An In-Silico Investigation of Phytochemicals as Antiviral Agents Against Dengue Fever.

A virtual screening analysis of our library of phytochemical structures with dengue virus protein targets has been carried out using a molecular docking approach. A total of 2194 plant-derived secondary metabolites have been docked. This molecule set comprised of 290 alkaloids (68 indole alkaloids, 153 isoquinoline alkaloids, 5 quinoline alkaloids, 13 piperidine alkaloids, 14 steroidal alkaloids, and 37 miscellaneous alkaloids), 678 terpenoids (47 monoterpenoids, 169 sesquiterpenoids, 265 diterpenoids, 81 steroids, and 96 triterpenoids), 20 aurones, 81 chalcones, 349 flavonoids, 120 isoflavonoids, 74 lignans, 58 stilbenoids, 169 miscellaneous polyphenolic compounds, 100 coumarins, 28 xanthones, 67 quinones, and 160 miscellaneous phytochemicals. Dengue virus protein targets examined included dengue virus protease (NS2B-NS3pro), helicase (NS3 helicase), methyltransferase (MTase), RNA-dependent RNA polymerase (RdRp), and the dengue virus envelope protein. Polyphenolic compounds, flavonoids, chalcones, and other phenolics were the most numerous of the strongly docking ligands for dengue virus protein targets.

INTRODUCTION

Dengue fever is an acute viral infection caused by all four serotypes (1, 2, 3, or 4) of the dengue virus (DENV), and is one of the seventeen “Neglected Tropical Diseases” (NTDs) defined by the World Health Organization [1]. The highest incidence of dengue fever is in Southeast Asia, India, and the Neotropics; where the principal vector, the Aedes aegypti mosquito, is found [2]. There are two well defined manifestations of dengue virus infection in humans, dengue fever and severe dengue (dengue hemorrhagic fever / dengue shock syndrome, DHF/DSS) [3]. Dengue viral infections account for around 100 million cases of dengue fever and more than 500,000 cases of DHF/DSS each year [4]. India accounts for around 34% of the total global dengue infections while tropical America contributes around 14% [1, 5]. Delhi, India has had several epidemic outbreaks of dengue fever (1970, 1982, 1988, 1996, 2006, 2010, and 2013) [6]. In the 2013 outbreak, 71% of suspected patients tested positive for dengue virus, mainly DENV-2, but co-infection by different serotypes was also detected [7]. Rio de Janeiro, Brazil, has suffered similar epidemic outbreaks of dengue fever (1986, 1990, 2001-02, 2007-08, 2011, and 2012) [8, 9]. Between 2010 and 2012, dengue virus infection was confirmed in 47.5% of the cases; DENV-2 the prevalent serotype in 2010, but DENV-1 re-emerged as the major serotype in 2011 [9]. In addition, between 2010 and 2012, 7.3% of the cases were severe (DHF/DSS), due mainly to secondary dengue infection. In 2007, 53% of DHF/DSS cases occurred in children under 15 years of age [10]. Furthermore, a correlation between increased dengue infection risk and socioeconomic status has been demonstrated; poorer inhabitants with poor access to health services showed a significantly higher seroprevelence compared to wealthy inhabitants [10, 11]. A serological survey carried out in 2007-2008 in three neighborhoods of Rio de Janeiro showed more than half the population was seropositive (60.3% in Higienópolis, a densely populated urban neighborhood, 56.1% in Tubiacanga, a suburban residential neighborhood, and 77.4% in Palmares, a recently settled slum) [8].

Currently, there are no treatment options for dengue fever; there is no known cure or vaccine. There are, however, several dengue virus proteins that have been identified as potential drug targets [12]. Natural products have been and continue to be excellent sources of medicinal agents, as themselves or as templates for synthetic agents [13]. The people of developing nations often have limited access to health care and inadequate resources for pharmaceuticals, and traditional herbal medicines offer the potential for treatments of diseases in these regions. In this work, we have screened our virtual library of phytochemicals against several protein targets of dengue virus [14], included dengue virus protease (NS2B-NS3pro), helicase (NS3 helicase), methyltransferase (MTase), RNA-dependent RNA polymerase (RdRp), and the dengue virus envelope protein (Ep), using a molecular docking approach.

Dengue Virus Protease

DENV protease (NS2B-NS3pro) is a trypsin-like serine protease that cleaves the dengue polyprotein into individual proteins necessary for viral replication. DENV protease has been identified as a primary target for the development of dengue antiviral drugs [15]. There are two potential targets for DENV protease inhibition: (1) the active site of the protease [16], and (2) block association of NS3pro with its protein cofactor NS2B [17].

Dengue Virus Helicase

DENV NS3 helicase is an essential enzyme for viral replication and has been identified as a potential drug target. The enzyme shows both ATP-hydrolysis activity as well as RNA duplex unwinding activity. These two enzymatic activities are likely coupled with ATPase activity providing the chemical energy for RNA unwinding. Thus, there are two potential targets for DENV NS3 helicase, the RNA binding site [18] and the ATP binding site [12].

Dengue Virus Methyltransferase

DENV MTase is located on the N-terminus of NS5 and catalyzes the N-7 and 2ʹ-O methylations of the viral RNA cap [19]. These modifications are necessary for formation of the mature RNA cap structure. The MTase uses S-adenosylmethionine (SAM) as the methylating cofactor. There are binding sites for GTP and the guanine moiety of the RNA cap (the GTP pocket) and for SAM. DENV MTase is a valid antiviral target, and the GTP binding site has been identified as a valid site for development of antivirals [20]. Although S-adenosylmethionine associates tightly with DENV MTase and displacement of SAM may be difficult, the SAM binding site may, however, be capable of binding other ligands [12, 21].

Dengue Virus RNA-Dependent RNA Polymerase

DENV RdRp is essential for viral replication and there are no similar enzymes found in the host. This enzyme, then, is considered to be a prime target for antiviral development [12]. Nucleoside-based prodrug antivirals have been investigated. These agents are phosphorylated in vivo, and then target the active site of the RdRp enzyme [20]. Non-nucleoside small molecules can inhibit viral RdRp by binding to allosteric sites of the protein [19].

Dengue Virus Envelope Protein

DENV E protein mediates the binding of the virus to the host cell surface receptors [19]. Thus, interference of this interaction would be a potential avenue for antiviral development. A small hydrophobic channel has been identified in the DENV E protein and this site has been identified as a target for small-molecule inhibitors [12].

The purpose of this work was to examine whether higher plants may serve as a source of antiviral agents for treatment of dengue fever. A molecular docking approach was used to examine a selection of phytochemicals from our virtual library (2194 different compound structures) with the protein targets of dengue virus.

COMPUTATIONAL METHODS

Protein-ligand docking studies were carried out based on the crystal structures (PDB 3U1J [22], 4M9K, 4M9M, and 4M9T [23]) and the NMR structures (PDB 2M9P and 2M9Q [24]) of dengue virus NS3 protease; crystal structures of dengue virus NS3 helicase (PDB 2BHR, 2BMF [25], 2JLQ, 2JLR, 2JLS, 2JLU, 2JLV, 2JLX [26], and 2WHX [27]); crystal structures of dengue virus NS5 RNA methyltransferase (PDB 3P8Z [28], 4CTJ, 4CTK, and 4R8S [21]); crystal structures of dengue virus NS5 RNA-dependent RNA polymerase (2J7U, 2J7W [29], 3VWS, and 4HHJ [30]); and crystal structures of dengue virus envelope protein (PDB 1OAN and 1OKE [31]). In order to examine selectivity for DEVN MTase, human RNA methyltransferase proteins (PDB 3BGV [32], 3EPP [33], and 4N49 [34]) were also used for docking. Similarly, in order to examine selectivity for DENV protease over mammalian serine proteases, we have carried out docking of the phytochemical ligands with crystal structures of rat trypsin (PDB 1BRB [35]), bovine trypsin (PDB 2ZDK [36] and 4MTB [37]), and porcine elastase (PDB 7EST [38]). Prior to docking, all solvent molecules and the co-crystallized ligands were removed from the structures. Molecular docking calculations for all compounds with each of the proteins were undertaken using Molegro Virtual Docker (version 6.0, Molegro ApS, Aarhus, Denmark) [39], with a sphere (15 Å radius) large enough to accommodate the cavity centered on the binding sites of each protein structure in order to allow each ligand to search. If a co-crystallized inhibitor or substrate was present in the structure, then that site was chosen as the binding site. If no co-crystallized ligand was present, then suitably sized (> 50 Å3) cavities were used as potential binding sites. Standard protonation states of the proteins based on neutral pH were used in the docking studies. Each protein was used as a rigid model structure; no relaxation of the protein was performed. Assignments of charges on each protein were based on standard templates as part of the Molegro Virtual Docker program; no other charges were necessary to be set. Overall, 2194 plant-derived secondary metabolites have been docked. This molecule set was comprised of 290 alkaloids (68 indole alkaloids, 153 isoquinoline alkaloids, 5 quinoline alkaloids, 13 piperidine alkaloids, 14 steroidal alkaloids, and 37 miscellaneous alkaloids), 678 terpenoids (47 monoterpenoids, 169 sesquiterpenoids, 265 diterpenoids, 81 steroids, and 96 triterpenoids), 20 aurones, 81 chalcones, 349 flavonoids, 120 isoflavonoids, 74 lignans, 58 stilbenoids, 169 miscellaneous polyphenolic compounds, 100 coumarins, 28 xanthones, 67 quinones, and 160 miscellaneous phytochemicals. While this list does not include all phytochemicals, it does represent a practical selection of phytochemical classes and structural types. Each ligand structure was built using Spartan ’14 for Windows (version 1.1.8, Wavefunction Inc., Irvine, California). For each ligand, a conformational search and geometry optimization

was carried out using the MMFF force field [40]. Flexible ligand models were used in the docking and subsequent optimization scheme. Variable orientations of each of the ligands were searched and ranked based on their re-rank score. For each docking simulation the maximum number of iterations for the docking algorithm was set to 1500, with a maximum population size of 50, and 100 runs per ligand. The RMSD threshold for multiple poses was set to 1.00 Å. The generated poses from each ligand were sorted by the calculated re-rank score.

RESULTS and DISCUSSION

In this analysis, we have rejected those phytochemical ligands that violate Lipinski’s rule of five [41]. That is, ligands with MW > 500 g/mol, hydrogen-bond-donating atoms > 5, hydrogen-bond-accepting atoms > 10, or ClogP > 5, even though they may have strong docking energies, were not considered.

The phytochemical ligands in this investigation were docked to the DENV NS2B-NS3 protease as well as with three different mammalian trypsin-like proteases (bovine trypsin, rat trypsin, and porcine elastase). Those ligands with the best docking for the DENV NS2B-NS3 protease (Edock ≤ -130 kJ/mol, comparable to the inhibitor (4S,5S)-4-amino-6,6,6-trifluoro-5-hydroxyhexyl-1-guanidinium) that were also selective for the virus protease over mammalian trypsin (at least 15 kJ/mol stronger docking) are listed in Table . The binding site of DENV NS2B-NS3 protease involves the catalytic triad of Ser196, His112, and Asp136. Additional amino acids at the active site are Glu46, Ile97, and Glu47. The binding site of this serine protease, as evidenced by the NMR structure with bound dipeptide Lys-Arg (PDB 2M9Q) [24], is bounded by Glu46 and Glu47, which both form electrostatic ionic interactions with the -NH3+ moiety of Lys-Arg. Other amino acids at the binding site include Val113, Ile97, Gln88, and Ile91.

As a class, the strongest docking ligands were the chalcones (average docking energy = -121.2 kJ/mol) and the stilbenoids (average docking energy = -119.5 kJ/mol). The worst docking compounds were the terpenoids, monoterpenoids (average Edock = -63.8 kJ/mol), sesquiterpenoids (average Edock = -89.6 kJ/mol), diterpenoids (average Edock = -89.7 kJ/mol), and triterpenoids (average Edock = -91.8 kJ/mol). The lowest-energy docking poses of the chalcone ligands balsacone A, balsacone B, and balsacone C, place these compounds into both the binding site of Lys-Arg as well as the catalytic site occupied by the co-crystallized inhibitor M9P (see Fig. ).

Previous docking studies with DENV NS2B-NS3 protease by de Sousa and co-workers had focused on flavonoids that had shown protease inhibition activity [42]. The most active compound in that study was agathisflavone, but in our docking analysis, this compound docks to DENV protease with higher energy (Edock = -126.1 kJ/mol) than it does with rat trypsin (Edock = -138. 6 kJ/mol). Additionally, this biflavone has one rule-of-five violation (MW = 522.46 g/mol). Other flavonoid ligands examined by de Sousa et al.,

myricetin, quercetin, and kaempferol, were relatively weak docking ligands (Edock = -105.5, -102.7, and -91.5 kJ/mol, respectively) compared to those shown in Table .

ul Qamar and co-workers had carried out a docking analysis with DENV NS2B-NS3 protease using a library of 940 phytochemicals [43]. These workers found Garcinia phytochemicals to be their best hits, including gossypol, mangostenone C, garcidepsidone A, and dimethyl-calabaxanthone. We have docked these compounds in this present study, but they do not dock as well as the ligands listed in Table (Edock = -112.6, -116.8, -111.7 kJ/mol for gossypol, mangostenone C, and demethylcalabaxanthone, respectively). Garcidepsidone A had a docking energy of -127.4 kJ/mol, but this compound violates Lipinski’s rule of five (ClogP = 5.74).

A docking study with compounds from Murraya koenigii and DENV NS2B-NS3 protease revealed the bisindole alkaloid bismurrayafoline E to be a promising binding ligand [44]. This compound, however, violates Lipinski’s rule of 5 (MW = 724.98 g/mol), and was not considered in our study. In another docking study, panduratin A and 4-hydroxypanduratin A, two phenolic compounds from Boesenbergia rotunda with dengue-2 NS3 protease inhibitory activity [45], were found to dock well with the active site of DENV NS2B-NS3 protease [46]. In this present docking study, we find only modest docking energies for panduratin A and 4-hydroxypanduratin A (-97.1 and -90.9 kJ/mol, respectively) compared to those listed in Table 1, suggesting that much more potent dengue virus NS2B-NS3 protease inhibitors are available.

Table 1

Docking energies (kJ/mol) of phytochemical ligands with DENV NS2B-NS3 protease and with mammalian trypsin.

a M9P = the inhibitor (4S,5S)-4-amino-6,6,6-trifluoro-5-hydroxyhexyl-1-guanidinium [24].

Glycyrrhiza extracts have shown antiviral activity against human immunodeficiency virus [47, 48], Newcastle disease virus [49], human respiratory syncytial virus [50] and herpes simplex virus 1 [51]. Prenylated stilbenoids (3,3ʹ,4,5ʹ-tetrahydroxy-5-prenylbibenzyl, 3,3ʹ,5ʹ-trihydroxy-4-methoxy-5-prenylbibenzyl, 3-acetoxy-4ʹ,5-dihydroxy-3ʹ-prenyldihydrostilbene, licobenzofuran) and isoflavonoids (glycyrrhisoflavone, 4ʹ-O-methylglycyrrhisoflavone) chalcones (kanzonol Y) from G. glabra have demonstrated outstanding docking properties to DENV protease (Table 1). Prenylated phenolics from G. glabra have inhibited HIV giant cell formation [47]. In addition, kanzonol Y docked well with DENV RNA-dependent RNA polymerase and DENV envelope protein (see below). Glabraisoflavanone, also from G. glabra, docked well with DENV methyltransferase (see below).

The balsacones A-C are antibacterial dihydrochalcones from Populus balsamifera [52]. We are unaware of any antiviral properties of these compounds, but they do exhibit excellent docking properties to DENV NS2B-NS3 protease (Table ) and DENV helicase (see below, Table ). The prenylated flavonoids glabranine and 7-O-methylglabranine, isolated from Tephrosia spp., had shown in-vitro antiviral activity against dengue virus serotype 2 [53]. Although these two flavonoids were not among the strongest docking, they did show docking selectivity with DENV protease (Edock = -113.4 and -118.8 kJ/mol, respectively). Two diprenylated flavanones, euchrestaflavanone A from Macaranga pleiostemona [54] and flemiflavanone D from Flemingia stricta [55] showed notably strong docking with DENV NS2B-NS3 protease (Table ). These compounds had shown antibacterial activity, but to our knowledge have not been screened for anti-dengue activity.

In addition to docking energies, it is useful to consider proximity of electrophilic ligands (e.g., epoxides or enones) to nucleophilic residues in the active site. A search of the docked ligands with DENV protease has revealed three sesquiterpenoid α-methylene lactones that dock in close proximity to the –OH group of Ser196. Ineupatorolide A, bigelovin, and 2-de-ethoxy-2-methoxyphantomolin (Fig. ) dock with the methylene carbon within 3.6 Å of the serine –OH. These ligands can presumably undergo conjugate addition with the –OH group and form a covalent intermediate thereby inhibiting the protease.

The strongest docking (most exothermic) ligands that do not violate Lipinski’s rule of five are listed in Table along with their docking energies to the RNA binding site and the ATP binding site. These ligands were selected based on their docking energies compared to the docking energies of the co-crystallized ligands in the ATP binding site, ANP (phosphoaminophosphonic acid adenylate ester, -156.9 kJ/mol) and ADP (adenosine 5’-diphosphate, -142.1 kJ/mol). Note that the phytochemical ligands showed a notable docking preference for the ATP site of DENV helicase over the RNA binding site. The ATP site of DENV helicase is dominated by protonated basic amino acids (Arg463, Lys199, Lys201, Arg460) that form ionic interactions and hydrogen bonds with the phosphate groups of ATP and the co-crystallized ligands ADP or ANP. Additional hydrogen-bonding amino acids include Gly196, Gly198, and Thr200. Lys201 and Arg418 form hydrogen-bonds with the adenine moiety.

As was the case with DENV protease, the strongest docking class of phytochemicals with DENV helicase ATP site were the chalcones (average Edock = -126.1 kJ/mol) and the stilbenoids (average Edock = -117.6 kJ/mol) while the worst were the terpenoids. The strongest docking ligands were curcumin II (Edock = -164.0 kJ/mol), caffeoyl-p-coumaroyltartaric acid (Edock = -156.2 kJ/mol), balsacone C (Edock = -153.6 kJ/mol), curcumin I (Edock = -153.4 kJ/mol), and angusticornin B (Edock = -152.5 kJ/mol). These ligands dock into the phosphate site of the ATP binding site and not in the adenine site of co-crystallized ADP or ANP. The common structural features of the best docking ligands are two phenolic groups with a flexible linker between them. One of the phenolic groups occupies the phosphate site while the other occupies a site flanked on either side by Arg463 and Arg418 (Figure 3). Several cinnamic acid derivatives from black cohosh (Cimicifuga racemosa) [56] have shown exceptional docking to the ATP site of DENV helicase, including cimicifugic acids B and G, cimiciphenol, and cimiracemate D (Table ).

Fig. (3)

(A) Lowest-energy docking poses of angusticornin B (green), balsacone C (yellow), caffeoyl-p-coumaroyltartaric acid (aqua), and curcumin II (magenta) with DENV helicase (PDB 2JLX) [26]. (B) Curcumin II in the active site of DENV helicase (PDB 2JLX) showing key intermolecular contacts (hydrogen-bonds are shown as blue dashed lines). The ADP co-crystallized ligand is shown as a yellow wire figure.

DENV MTase has a relatively large binding cavity (around 425-515 Å3) that accommodates the viral RNA cap (the GTP site) and SAM (the SAM site). The strongest docking phytochemical ligands all docked at the SAM site. In addition, it is beneficial for antiviral agents to show selectivity for DENV MTase over human RNA MTase. Table 3 lists those compounds that docked with DENV MTase with docking energies comparable or more exothermic than SAM and also with docking energies significantly more exothermic with DENV MTase than with Hs RNA MTase. Not surprisingly, the strongest docking ligands were relatively large molecules (MW around 680-802 g/mol), and therefore violate Lipinski’s rule of five. Based on this molecular docking analysis, the ligands in

Table 3

MolDock docking energies (kJ/mol) of phytochemical ligands with dengue virus methyltransferase and human RNA methyltransferase.

a SAM and SFG are co-crystallized ligands.

Table are predicted to show inhibition of DENV MTase by competition with SAM and to selectively inhibit viral MTase over human MTase. Additionally, there are several phytochemical ligands that showed more exothermic docking energies than SAM for the SAM site of DENV MTase, more exothermic docking with DENV MTase than with human RNA MTase, and no Lipinski’s rule violations: Isoborreverine, diplacone, styracifolin B, neosilyhermin A, and 8β-[4-hydroxy-5-(5-hydroxytigloyloxy)tigloyl]-sant-amarin.

Like DENV MTase, DENV RdRp has a very large binding cavity (ca. 3000 Å3), so there are several sites available for docking. Eight phytochemical ligands from our docking analysis showed excellent docking energies to DENV RdRp (Table ). These ligands meet the rule-of-five criteria for drug likeness and have significantly lower (more exothermic) docking energies than the co-crystallized ligands guanosine 5ʹ-triphosphate (GTP) [29] or 5-{[(4-chlorophenyl)sulfonyl]amino}-2-methyl-1-benzofuran-3-carboxylic acid (VWS) [30]. Four of the strongly docking ligands docked preferentially at the GTP site of DENV RdRp, dimethylisoborreverine, drummondin D, flinderole B, and pungiolide A. The other four ligands docked preferentially at sites removed from either the GTP binding site or the VWS site. Flinderole A and 4,6-dibenzoyl-2-[phenylhydroxymethyl]-3(2H)-benzofuranone preferentially docked in a hydrophobic pocket surrounded by Trp302, Phe354, Val358, Val577, Val579, and Gly599. Neosilyhermin B and kanzonol Y docked preferentially in a hydrophobic pocked formed by Ala406, Ala407, Asn492, Glu507, Val603, Tyr606, and Ile797.

Curcuma longa extracts have shown activity against hepatitis B virus replication [57] and dengue-2 virus protease [58]. In addition, curcumin, one of the major components of C. longa, has shown inhibition of HIV-1 and HIV-2 proteases [59] and in-vitro infectivity of dengue virus type 2 [60]. García Ariza and co-workers have suggested that dengue virus RNA-dependent RNA polymerase inhibition by curcumin, based on molecular docking studies, may be responsible for the anti-dengue activity of curcumin [61]. In this present docking study, several curcuminoids (1,7-bis(4-hydroxy-3-methoxyphenyl)-1,4,6-heptatrien-3-one, curcu-min I, curcumin II, curcumin III, and dihydrocurcumin) showed notably strong docking energies to the ATP site of dengue virus helicase (Table ). These ligands also showed docking energies somewhat stronger than those of the co-crystallized ligands for DENV RdRp (e.g., curcumin I, Edock = -119.3 kJ/mol). Curcumin I also showed excellent docking properties to the DENV envelope protein hydrophobic pore (see below, Table ).

A number of alkylated flavanones, the chartaceones A-F, isolated from the bark of Cryptocarya chartacea, have shown potent activity against dengue virus NS5 RdRp [62]. The chartaceones C-F had shown the highest potency and these ligands also showed very exothermic docking energies (ranging from -125 to -139 kJ/mol). These compounds, however, all violate Lipinski’s rule (MW = 660.75 g/mol). The lower molecular weight chartaceones A-B, however, docked with DENV RdRp with docking energies (Edock = -101 to -117 kJ/mol) comparable to GTP (Edock = -107.7 kJ/mol). Chartaceone A2 showed selective docking to DENV protease (Table 1), chartaceone B3 docked strongly with the ATP site of DENV helicase (Table ), and chartaceone B4 docked well with DENV MTase (Table ).

Several phytochemical compounds have been found that dock into the small hydrophobic channel of DENV envelope protein with docking energies significantly more exothermic than the co-crystallized inhibitor, n-octyl-β-d-glucoside [12] (Table ). Additionally, these phytochemical ligands completely occupy the channel of the envelope protein (see Fig. ). Thus, for example, the flavonoid canniflavin A docks inside the cavity pore, which is largely hydrophobic. There are, however, key hydrogen-bonding contacts between the flavonoid ligand and Thr48 and Thr280 (Fig. ). An obvious structural feature of strongly docking phytochemical ligands to the hydrophobic pore is the ability to form an extended structure, particularly involving hydrophobic groups such as prenyl and geranyl.

A number of flavonoids from Silybum marianum have shown exceptional docking to the hydrophobic pore of DENV Ep (2,3-dehydroxilybin, cisilandrin, isosilybins A-D, silandrins A-B, and silybin (Table ). Silybin (silibinin, silymarin) has shown antiviral activity against hepatitis C virus [63] and chikungunya virus [64]. C-geranylated flavanones from Paulownia tomentosa fruits have shown antibacterial [65, 66], cytotoxic [67], and anti-inflammatory [68] activities as well as SARS-CoV protease inhibitory activity [69]. Diplacone showed excellent docking properties with DENV MTase (Table ) and both 3ʹ-O-methyldiplacone and 3ʹ-O-methyldiplacol docked well with the ATP site of DENV helicase (Table ). Diplacone, bonannione A, 3ʹ-O-methyldiplacone and 4ʹ-O-methyldiplacone docked well with the hydrophobic channel of DENV Ep (Table ). Likewise, the C-geranylated flavone cannflavin A, from Cannabis sativa [70], also docked well with DENV Ep (Table ).

SUMMARY AND CONCLUSION

A molecular docking analysis of 2194 plant-derived secondary metabolites with dengue virus protein targets has been carried out. The analysis has revealed 24 compounds that docked strongly to dengue virus NS2B-NS3 protease (Edock < -130 kJ/mol) but not as strongly with mammalian trypsin, 21 compounds that docked strongly to the ATP binding site of DENV NS3 helicase (Edock < -140 kJ/mol), 13 compounds that docked strongly to DENV methyltransferase (Edock < -130 kJ/mol) but not as strongly with human RNA methyltransferase, 8 phytochemicals that showed notable docking to DENV RNA-dependent RNA polymerase, and 32 compounds that showed excellent docking properties with the hydrophobic pore of DENV envelope protein. The results of this study mirror previous docking studies that showed polyphenolic phytochemicals to be excellent docking ligands to dengue protein targets [42, 43, 46, 61]. These results and the current study underscore the importance of natural products from higher plants in drug discovery and may provide potential avenues for development of chemotherapeutic agents for the treatment of dengue fever.

Table 2

Lowest-energy docking scores (re-rank scores) for phytochemical ligands with dengue virus helicase.

Table 4

Lowest MolDock docking energies for phytochemical ligands with dengue virus RNA-dependent RNA polymerase.

Table 5

MolDock docking energies (kJ/mol) for phytochemical ligands with dengue virus envelope protein hydrophobic channel.