Ellagic Acid as a Potential Inhibitor against the Nonstructural Protein NS3 Helicase of Zika Virus: A Molecular Modelling Study

Kullappan, Malathi; Benedict, Balakrishnan Anna; Rajajagadeesan, Anusha; Baskaran, Padmasini; Periadurai, Nanthini Devi; Ambrose, Jenifer Mallavarpu; Gandhamaneni, Sri Harshini; Nakkella, Aruna Kumari; Agarwal, Alok; Veeraraghavan, Vishnu Priya; Surapaneni, Krishna Mohan

doi:https://doi.org/10.1155/2022/2044577

BioMed Research International

On this page

Abstract Introduction Materials and Methods Results Discussion Conclusion Data Availability Consent Conflicts of Interest Authors’ Contributions Acknowledgments Supplementary Materials References Copyright Related Articles

Research Article | Open Access

Volume 2022 | Article ID 2044577 | https://doi.org/10.1155/2022/2044577

Ellagic Acid as a Potential Inhibitor against the Nonstructural Protein NS3 Helicase of Zika Virus: A Molecular Modelling Study

Malathi Kullappan,¹Balakrishnan Anna Benedict,²Anusha Rajajagadeesan,³Padmasini Baskaran,⁴Nanthini Devi Periadurai,⁵Jenifer Mallavarpu Ambrose,¹Sri Harshini Gandhamaneni,⁶Aruna Kumari Nakkella,⁷Alok Agarwal,⁸Vishnu Priya Veeraraghavan,⁹and Krishna Mohan Surapaneni¹⁰

Academic Editor: Mohd Saeed

Received14 Apr 2022

Revised27 May 2022

Accepted01 Aug 2022

Published21 Aug 2022

Abstract

Zika virus is a member of the Flaviviridae family and genus Flavivirus, which has a phylogenetic relationship with spondweni virus. It spreads to humans through a mosquito bite. To identify potential inhibitors for the Zika virus with biosafety, we selected natural antiviral compounds isolated from plant sources and screened against NS3 helicase of the Zika virus. The enzymatic activity of the NS3 helicase is associated with the C-terminal region and is concerned with RNA synthesis and genome replication. It serves as a crucial target for the Zika virus. We carried out molecular docking for the target NS3 helicase against the selected 25 phytochemicals using AutoDock Vina software. Among the 25 plant compounds, we identified NS3 helicase-ellagic acid (-9.9 kcal/mol), NS3 helicase-hypericin (-9.8 kcal/mol), and NS3 helicase-pentagalloylglucose (-9.5 kcal/mol) as the best binding affinity compounds based on their binding energies. To understand the stability of these complexes, molecular dynamic simulations were executed and the trajectory analysis exposed that the NS3 helicase-ellagic acid complex possesses greater stability than the other two complexes such as NS3 helicase-hypericin and NS3 helicase-pentagalloylglucose. The ADMET property prediction of these compounds resulted in nontoxicity and noncarcinogenicity.

1. Introduction

Zika virus (ZIKV) is an arbovirus related to the family of Flaviviridae [1]. It is a positive sense RNA, which is a constituent of the genus Flavivirus. This genus includes other viruses such as tick-borne encephalitis, yellow fever, tick-borne encephalitis, dengue, West Nile, Langat, St. Louis encephalitis, Modoc, Powassan, Japanese encephalitis, and Rio Bravo viruses [2]. The genome length of ZIKV is 10,794 kb, where the RNA has two noncoding regions such as 39 and 59 with a long open reading frame. ZIKV differs from other flaviviruses through the glycosylation spots located on the surface of the virus [3]. The ZIKV genome has three structural proteins and seven nonstructural proteins, which mediate the replication of the genome [4]. Here, the viral transmission follows the sylvatic cycle (nonhuman primate’s cycle), where the virus is disseminated by the bite of mosquitoes, namely, Aedes albopictus, Aedes aegypti, Aedes africanus, and Aedes hensilli, infected with ZIKV [5]. It can affect cells like dermal fibroblasts, skin keratinocytes, and dendritic cells. In vitro studies reveal that these cells demonstrated an increased level of infection after 24 to 48 hours of viral entry [6]. After biting an infected person, the mosquitoes consume the blood with ZIKV that initiate replication in the midgut epithelial cells and transfer to the salivary gland. After the incubation time of 10 days, the saliva is infected, which starts spreading by biting the other human. Through clathrin-mediated endocytosis, the viruses enter the host cell. When ZIKV enters the host skin, initially it affects the dermal fibroblasts that serve as a receptor for the viral attachment. Then, the membrane of the virus fuses with the endosomal membrane and would be delivered into the mammalian cells [7]. After completion of the replication, the virus gets multiplied and affects the lymph nodes and nervous system [8]. It also spreads through sexual intercourse and from the mother to the fetus through the placenta at the time of pregnancy [9]. These situations cause microcephaly and congenital Zika syndrome in the fetus. Neurological problems like Guillain-Barre syndrome (in adults—autoimmune disease), myelitis, and neuropathy were also reported in ZIKV infection [10]. The symptoms include mild rash, fever, joint and muscle pain, conjunctivitis, malaise, or headache, including symptomatic cases. It has been reported that terrorists are responsible for the vigorous transmission of ZIKV worldwide in around 60 nations [11]. In adults, the infection is associated with other defects such as thrombocytopenia, meningitis, multiorgan failure, and encephalitis. Mortality has been reported with sickle cell disease in children and cancer in adults. The persistence of ZIKV in the eye leads to the replication of the virus and causes maculopathy, uveitis, and conjunctivitis, which may tend to blindness in many cases [12]. Temperature is the greatest driver of vector-borne disease transmission; the mosquito transmission decreases with temperature. When the disease transmission occurs at 29°C, it was controlled under a cool temperature of <22°C. Also, the effectiveness of genome replication was reduced at cool temperatures [13]. The extensive outbreaks of ZIKV and its neurological defects have created a greater public health concern worldwide. There is an urgent need to identify the most promising inhibitor against ZIKV. Most of the existing antivirals are reported to have side effects. Hence, our current study exploited plant-derived phytochemicals against ZIKV, which has less toxicity as well as less chance to develop resistance. Usually, these plant-derived compounds were exploited to treat many diseases. Also, these plant compounds can inhibit replication. Alkaloids such as cepharanthine, fangchinoline, and tetrandrine isolated from Stephania tetrandra S. Moore have the ability to conceal the viral protein expression and thereby inhibit the human coronavirus OC43 replication in MRC-5 cells [14]. Synthesized or plant-derived quinolone alkaloids inhibit the activity of type II topoisomerase which leads to the inhibition of DNA replication [15]. Hence, our study has targeted the nonstructural ZIKV NS3 helicase, a crucial enzyme that participates in RNA unwinding of the replication process. The helicases of the Flavivirus are also informed to contribute to other vital roles like the splicing of pre-mRNA, ribosome biogenesis, export and degradation of RNA, maturation of RNA, and translation process. This NS3 helicase is a member of the superfamily SF2, which has a close identity with the Murray Valley encephalitis virus, dengue virus 2 (DENV2), and DENV4. The N terminal region of the protein possesses the activity of protease, and the C terminal region possesses the activity of helicase. Though the helicases of the Flavivirus family have conserved active site regions, ZIKV NS3 helicase has different RNA binding modes and motor domain movements [16]. All these crucial roles of ZIKV NS3 helicase inspired us to screen plant-derived phytochemicals. Various studies have been conducted to screen inhibitors against ZIKV; they used crucial targets such as NS3 helicase [17], NS2B/NS3 protease [18] [19, 20], and glycoprotein (gp) E [21]. The current study involves the use of docking protocols to analyze the binding mode of 25 plant compounds with the ZIKV NS3 helicase, molecular dynamic studies were further focused to analyze the protein-ligand complex stability, and finally, the oral bioavailability of the best binding energy (BE) ligands was also predicted. The overall study reveals the favorable BE, stability, and bioavailability of the plant compound ellagic acid. Our study will be helpful to the scientist who is involved in the inhibitor design for the ZIKV.

2. Materials and Methods

2.1. Preparation of the Target

The ZIKV NS3 helicase structure was obtained through the database Protein Data Bank (PDB: 5JRZ) [22, 23]. It is a monomer with a resolution of 1.62 Å. The X-ray diffraction technique was used to determine the crystal structure (Figure 1). It has a single chain with an amino acid length of 449. Nonresidue atoms coordinated with the crystal structure were removed. The missing residues and the atoms were modeled through the MODELLER 9.25 server [24]. Figure 2 depicts the study’s whole workflow. The figure was drawn with the help of the BioRender server [25].

2.2. Preparation of Active Sites

The amino acid information of the ATP binding site region in the ZIKV NS3 helicase (PDB: 5JRZ) structure was recovered from the literature [22].

2.3. Ligand Preparation

A total of twenty-five plant compounds with antiviral activity were collected from various studies [26–49], and their 3D structures were collected from the PubChem database [50]: apigenin (ID: 5280443), baicalein (ID: 5281605), berberine (ID: 2353), betulin (ID: 72326), chebulagic acid (ID: 442674), curcumin (ID: 969516), ellagic acid (ID: 5281855), epigallocatechin gallate (ID: 65064), fisetin (ID: 5281614), geraniin (ID: 3001497), glycyrrhizic acid (ID: 14982), hypericin (ID: 3663), hyperoside (ID: 5281643), kaempferol (ID: 5280863), lupeol (ID: 259846), mimusopic acid (ID: 6712545), mulberroside C (ID: 190453), myricetin (ID: 5281672), neoandrographolide (ID: 9848024), pentagalloylglucose (ID: 65238), piperine (ID: 638024), quercetin (ID: 5280343), rosmarinic acid (ID: 5281792), rutin (ID: 5280805), and torvoside (ID: 11018078). Also, N-(3-acetylphenyl) morpholine-4-carboxamide (ID: 671267), 2-(2-acetamidophenyl) acetic acid (ID: 14622178), and 2-(4-acetamidophenyl) acetamide (ID: 25862032) were obtained from the PubChem database.

2.4. Drug-Likeness Property Prediction for Phytochemicals

The drug candidates usually act upon the cellular targets; the molecule binds to the target and changes the cellular machinery and produces therapeutic action. A molecule before producing the pharmacodynamics effect on the human body should travel from the entry point to the active site, which is referred to as the pharmacokinetic properties. It was defined by the terms like absorption, distribution, metabolism, and excretion. To predict such pharmacokinetic properties of the drugs, Lipinski and coworkers have established a set of rules known as “Lipinski’s rule of five.” This rule states that the H-bond donors should be below 5, H-bond acceptors should not exceed 10, mlogP value should be below 5, and molecular weight should be below 500 Daltons [51]. These rules filter the compounds and make the drug candidates orally available. To analyze the oral bioavailability of the plant-derived compounds, the drug-likeness property of the compounds was predicted through the Molinspiration server [52].

2.5. Molecular Docking Studies

The interaction of small molecules with the target protein can be understood by molecular docking studies. For our study, we have utilized the AutoDock Vina, an open-source tool [53]. The structures of target NS3 helicase (PDB: 5JRZ) and 25 plant compounds were submitted as the input file. The protonation state of the amino acid side chains was predicted through the online server PropKa [54]. The stereochemistry of all the compounds was analyzed properly. Through the SWISS-PdbViewer, all the ligands were energy minimized [55]. SWISS-PdbViewer employs the steepest descent energy minimization through the GROMOS96 force field [56]. The grid size was fixed as with 0.375 Å spacing. Also, grid center XYZ coordinates were set as 39.629, 27.127, and 51.699, respectively. AutoDock Vina uses iterated local search global optimizer (a genetic algorithm with local gradient optimization), a stochastic global optimization algorithm, and hybrid scoring functions (a combination of empirical and knowledge-based scoring functions) [57]. AutoDock Vina produces the ten best binding modes for each run. The negative binding energy designates the prediction of ligand binding to a target protein. The higher negative value of the binding energy indicates the favorable interaction between the protein and ligand.

2.6. Molecular Dynamic Simulations

The simulation studies were executed to explore the stability of the selected bound complexes. From the molecular docking results, the three best BE complexes, NS3 helicase-ellagic acid, NS3 helicase-hypericin, and NS3 helicase-pentagalloylglucose were subjected to molecular dynamic simulations for 100 ns through GROMACS package in the Ubuntu environment [58] using Intel Xeon W-1270, 8 core, and 16 threaded processors. The complexes were located in the cubic box, and using transferable intermolecular Potential 3P (TIP3P) water, the system was solvated. The topology of the protein was produced through the Chemistry at Harvard Macromolecular Mechanics 36 (CHARMM36) force field, and for ligand, it was created by CHARMM General Force Field (CGenFF) 4.4 version [59]. The whole system was neutralized by adding the Na⁺ or Cl^- ions. Energy minimization of the complete system was carried out using the steepest descent algorithm. Next, two sets of equilibration were performed, one is NVT (canonical ensemble) run for 0.1 ns and another one is NPT (isothermal–isobaric ensemble) run for 0.1 ns. The bond length of the protein was restrained using the LINCS algorithm [60], and the Particle Mesh Ewald method was used to measure the electrostatic interactions. Finally, the production run of 100 ns was planned with 1.0 bar pressure and 300 K temperature. The trajectory analysis was made to produce different plots like Root Mean Square Deviation (RMSD), Root Mean Square Fluctuation (RMSF), Radius of Gyration (Rg), H-bond interactions, and Molecular Mechanics Poisson-Boltzmann Surface Area (MMPBSA) calculation of binding free energy which were accomplished to understand the stability of the complexes.

2.7. In Silico Absorption, Distribution, Metabolism, Excretion, and Toxicity (ADMET) Prediction

In the current study, we have used the admetSAR online tool for the ADME predictions [61]. An ADMET structure-activity relationship database has ADMET-related properties collected from the literature. It has 210000 data points for more than 96000 unique compounds. In absorption, BBB permeability, HIA, Caco-2 absorption, P-gp (drug transporter) inhibitor and substrate, and renal OCTs were analyzed. In metabolism, the cytochrome (CYP) P450 substrate 2C9, 2D6, 3A4 and CYP P450 inhibitors 1A2, 2C9, 2C19, 2D6, and 3A4 were analyzed. Also, toxicity and carcinogenicity of the small molecules were predicted by the admetSAR server.

3. Results

A protein-ligand docking and molecular dynamic simulations were conducted to locate the potential NS3 helicase inhibitors.

3.1. Active Site Prediction

According to Jain et al. (2016), the ZIKV NS3 helicase (PDB: 5JRZ) structure has three domains (domains 1, 2, and 3) with more or less similar sizes and two binding sites such as ATP binding site located between domains 1 and 2 and the RNA binding groove located between 1 and 3 (Figure 3). For the current study, we have retrieved the ATP binding site amino acids (Ala416, Glu231, Arg459, Glu286, Gln455, Arg202, Arg462, Pro196, and Lys200) to carry out the docking studies [22].

3.2. Drug-Likeness Property Prediction of Plant Compounds

All the 25 plant-derived antiviral compounds collected from the literature (Figure S1) were subjected to drug-likeness property prediction; out of which, a few compounds such as chebulagic acid, geraniin, glycyrrhizic acid, pentagalloylglucose, rutin, and torvoside showed 3 violations, where their molecular weight, hydrogen bond donor (HBD), and hydrogen bond acceptor (HBA) were found higher than the limited values. Lipinski’s rule of the five compounds violated only two rules. Similarly, the compounds epigallocatechin gallate and hyperoside had 2 violations. Here, we have observed several HBA and HBD counts. Similarly, betulin, lupeol, and myricetin had one violation, which was due to excess octanol partition coefficient values of botulin and lupeol, respectively; myricetin HBD count was found to be higher than the normal value (Table 1).

3.3. Molecular Docking Studies

Docking of the 25 plant compounds was carried out through the AutoDock Vina. The binding energies resulted in the range of -7.1 kcal/mol to -9.9 kcal/mol for the plant compounds against nonstructural protein NS3 helicase (PDB: 5JRZ) (Table 2). Here, favorable binding energy was observed for ellagic acid. The compounds ellagic acid, hypericin, pentagalloylglucose, epigallocatechin gallate, rutin, and glycyrrhizic acid fell within the binding energy ranges from -9.1 kcal/mol to -9.9 kcal/mol. Compounds such as mulberroside, myricetin, quercetin, baicalein, fisetin, torvoside, hyperoside, kaempferol, lupeol, neoandrographolide, and rosmarinic acid showed their binding energies within -8.0 kcal/mol to -8.9 kcal/mol. The remaining compounds like apigenin, chebulagic acid, mimusopic acid, berberine, betulin, curcumin, geraniin, and piperine fall within the binding energies between -7.1 kcal/mol and -7.9 kcal/mol.

The binding free energy was contributed by the interactions of electrostatic and nonelectrostatic forces. Here, the contribution of both the interactions was quantified by Coulomb’s potential and Lennard Jones potential, respectively. The contribution of each force is based on the charge and shape complementarity of the protein-ligand complex.

From Table 3, the most favorable interacting complexes were selected. Here, NS3 helicase-ellagic acid complex was ranked first as the interaction energy was observed to be -9.9 kcal/mol. It exhibited H-bond contacts with nine amino acids of Arg459 (bond length, 2.1 Å), Gln455 (2.4 Å), Glu231 (1.8 Å & 1.9 Å), Arg202 (2.1 Å), Thr201 (2.5 Å), Gly199 (2.8 Å), Ala198 (2.8 Å), Leu194 (2.0 Å), and His195 (2.1 Å) and hydrophobic interactions found with the amino acids of Arg776, Thr790, Leu858, Cys775, Leu777, Lys745, Asp855, Thr854, Met766, Leu844, Val726, Met793, Cys797, Ala743, Leu792, Leu718, Met1002, and Gly796 (Figure 4). These intermolecular H-bond and hydrophobic interactions are stabilizing the ligands in the binding sites of the NS3 helicase protein.

Hypericin was ranked second with -9.8 kcal/mol of binding energy and 5 amino acid interactions of Glu286 (2.7 and 3.4), Gly415 (2.6), Glu231, Thr201 (2.4), and Lys200 (2.2 and 2.3) and hydrophobic interactions of Gly415, Glu231, Ala416, Arg462, Arg202, Gly197, Asn417, and Gly199 (Figure 5).

Similarly, pentagalloylglucose was positioned third in the list with the binding energy of 9.5 kcal/mol and stabilized by 4 H-bonds with amino acids such as Asp410 (1.9 and 2.6), Met414 (2.0), Asp540 (1.9, 2.6, and 2.1), and Arg226 (2.0 and 2.1). It also displayed ten hydrophobic interactions with Asp291, Val227, Glu413, Arg388, Thr225, Ile411, Phe391, Cys262, Ala264, and Phe289 (Figure 6).

H-bonds are more important in defining the inhibitor’s activity against the target protein and in ensuring its stability with the protein. Here, the active site (ATP binding site) amino acids of NS3 helicases, such as Arg202, Arg459, Glu231, and Gln455, are involved in the interaction with ellagic acid; likewise, Glu286, Glu231, and lys200 amino acids were connected with the hypericin. None of the active site amino acids are found in the H-bond interactions between the NS3 helicase and pentagalloylglucose. Among the active site residues, Glu231 was found closely interacted with ellagic acid through 2 H-bonds with distances of 1.8 Å and 1.9 Å.

3.4. Molecular Dynamic Simulation Studies

The structural stability of the best binding energy complexes, NS3 helicase-ellagic acid, NS3 helicase-hypericin, and NS3 helicase-pentagalloylglucose complexes was analyzed by conducting molecular dynamic simulations for 100 ns, and the RMSD, RMSF, Rg, H-bond interactions, and MMPBSA calculation of binding free energy were analyzed.

3.4.1. RMSD

The RMSD is a measure of the difference between the initial structural confirmations of the protein backbone to its final position. The stability of the protein related to the structural conformation can be estimated by the deviations produced after the simulations. The RMSD plot for the Cα backbone atoms of the complexes NS3 helicase-ellagic acid, NS3 helicase-hypericin, and NS3 helicase-pentagalloylglucose was depicted in Figure 7(a). As shown in the figure, NS3 helicase-ellagic acid (blue) and NS3 helicase-hypericin (red) displayed slight variations in their RMSD values. At the beginning of the simulations, NS3 helicase-ellagic acid complex had the RMSD of 0.10 nm, which reached 0.12 nm until 10 ns. Similarly, the NS3 helicase-hypericin complex showed an RMSD value of 0.13 nm in the beginning and reached 0.15 nm till 10 ns; after which, both the complexes started showing deviation in their RMSD values. After 30 ns of simulation, both the complexes attained stability at 0.20 nm and 0.23 nm. Likewise, the NS3 helicase-pentagalloylglucose complex (black) showed fluctuation in the beginning and reached 0.25 nm after 10 ns. Further, as illustrated in the figure, we observed a fall in the RMSD value after 20 ns and attain stability at 0.27 nm. The unbound NS3 helicase has shown a stable RMSD of 0.42 nm after 30 ns (green). The overall comparison of the RMSDs of the three different complexes and unbound protein indicates the lower stability of the ellagic acid. Figure 7(b) depicts the RMSDs of the ligands ellagic acid (indigo), hypericin (cyan), and pentagalloylglucose (orange) after 100 ns. The ligands ellagic acid, hypericin, and pentagalloylglucose found their stability state at 0.02 nm, 0.06 nm, and 0.3 nm, respectively. The comparison of the RMSDs of the unbound NS3 helicase and free ligands with bound protein-ligand complex unfolds the impact of ligands binding with the NS3 helicase.

(a)

(b)

3.4.2. RMSF

Fluctuations of each amino acid in the NS3 helicase (PDB: 5JRZ), when bound to compounds such as ellagic acid, hypericin, and pentagalloylglucose, are illustrated in Figure 8(a). The NS3 helicase and ellagic acid exhibited their amino acid interactions at Arg459, Gln455, Glu231, Arg202, Thr201, Gly199, Ala198, Leu194, and His195. Likewise, NS3 helicase-hypericin interacted with the amino acids Glu286, Gly415, Glu231, Thr201, and Lys200. Also, NS3 helicase-pentagalloylglucose interacted with Asp410, Met414, Asp540, and Arg226, respectively. The unbound NS3 helicase has shown the greater fluctuation (green) when compared to the other three complexes and the least fluctuation was observed for NS3 helicase-ellagic acid complex (blue). These fluctuation ranges revealed the higher stability of the NS3 helicase-ellagic acid complex than that of the other two complexes NS3 helicase-hypericin and NS3 helicase-pentagalloylglucose. Figure 8(b) depicts the RMSF values for three ligands, ellagic acid (indigo), hypericin (cyan), and pentagalloylglucose (orange). The highest fluctuation was observed for pentagalloylglucose, followed by hypericin and ellagic acid. Hence, the binding of phytochemicals with NS3 helicase has a greater impact on the rigidity of the protein-ligand complex.

(a)

(b)

3.4.3. Rg

The calculation of Rg revealed the compactness of the protein after binding to the ligand. Figure 9 shows the Rg values of the complexes NS3 helicase-ellagic acid (blue), NS3 helicase-hypericin (red), and NS3 helicase-pentagalloylglucose (black) and unbound NS3 helicase (green). From the beginning of the simulation to 30 ns, all the three complexes attained stability at the Rg values of 2.29 nm. Later, a decrease in Rg value was observed for all the three complexes after 30 ns. Both NS3 helicase-ellagic acid and NS3 helicase-hypericin were observed to gain stability at 2.24 nm. The lack of variations in the complex revealed the rigidity of the NS3 helicase protein after binding with the ligands ellagic acid and hypericin. As illustrated in the figure, the complex NS3 helicase-pentagalloylglucose gained stability at 2.27 nm and unbound NS3 helicase gained stability at 0.29 nm. The variation in the Rg values between these complexes and free protein revealed the rigidity of the complex after binding to the ligands.

3.4.4. H-Bond Interactions of Protein-Ligand Complexes

The H-bond formation between the protein and ligands generates specificity and directionality, which is the basic aspect of molecular recognition. The HBD and HBA share their energy in the binding site region. Also, for further stability of the complex, it must be complemented with the van der Waals interactions because of the shape complementarity. After the completion of 100 ns of MD simulations, the H-bond interaction between the complexes NS3 helicase-ellagic acid, NS3 helicase-hypericin, and NS3 helicase-pentagalloylglucose was analyzed. As depicted in Figure 10, the NS3 helicase-ellagic acid complex (blue) formed 7 strong H-bonds throughout the simulation. The NS3 helicase-hypericin complex (red) with 5 strong H-bond formations and NS3 helicase-pentagalloylglucose (black) with 3 H-bond formations were observed throughout the simulations. These H-bond formations stabilized the protein-ligand complexes.

3.4.5. MMPBSA Calculation of Binding Free Energy

The binding free energy for the phytochemicals ellagic acid, hypericin, and pentagalloylglucose with the target NS3 helicase (PDB: 5JRZ) was calculated by the MMPBSA technique. The total binding free energy is the combination of van der Waals energy, electrostatic energy, polar solvation energy, and solvent-accessible surface area (SASA). The van der Waals energy, electrostatic energy, and SASA contribute negative energy values, and polar solvation energy contributes positive energy values to the total BE. The comparison of three binding energies revealed the highest total BE of ellagic acid () than that of the other two complexes, which indicated the strong interaction of ellagic acid with the ZIKV NS3 helicase (PDB: 5JRZ). Here, the van der Waals energy contributed more than the electrostatic energy (Table 4).

3.5. In Silico ADMET Prediction

Drug designing and development is highly a complex, time-consuming, and costlier process which has increased attrition rates. The development of in silico ADMET predictions has reduced the failure rate in recent years. ADMET prediction is an essential component in the analysis of the efficacy and toxicity of the drug. It also determines the administration dose, route, and frequency. The physicochemical property of the compound candidate may affect the ADMET properties. Hence, the assessment of ADMET property is vital for the successful development of drugs. The initial evaluations of pharmacokinetic and toxic properties are essential to avoid compounds with ADME problems, which helps the scientist to prioritize particular compounds to synthesize and evaluate. Prediction of ADMET properties of these plant-derived compounds may save the time and cost of drug discovery. The blood-brain barrier (BBB) controls the flow of chemicals into the brain. The role of BBB is to separate the flow of blood from the central nervous system. This prevents harmful substances to enter the brain tissues. BBB controls the movement of external compounds to preserve the central nervous system at a steady state. Here, ellagic acid and pentagalloylglucose may cross the BBB (BBB_positive) and hypericin may not cross the BBB (BBB_negative). The phytochemicals ellagic acid and hypericin may get absorbed by the intestine (HIA_positive) and pentagalloylglucose may not be absorbed by the intestine (HIA_negative). The Caco-2 permeability was analyzed to identify the intestinal absorption of the drugs; the results showed the Caco-2 permeability of hypericin (Caco2_positive) and nonpermeability of ellagic acid and pentagalloylglucose (Table 5). The P-gp enzyme acts as a drug transporter, which controls the uptake and efflux of a wide variety of drugs. It also serves as an efflux pump. P-gp facilitates drug-drug interactions. The substrates of P-gp can also act as an inhibitor of the enzyme; inhibition of the enzyme results in greater drug bioavailability. Here, these three plant compounds, namely, ellagic acid, hypericin, and pentagalloylglucose, function as the substrate for P-gp. Similarly, cytochrome P450 enzymes such as 2D6, 2C9, and 3A4, hypericin, and the ellagic acid act as a nonsubstrate, and pentagalloylglucose serves as a substrate. The AMES toxicity and carcinogenic test revealed the nontoxicity and noncarcinogenic effects of ellagic acid, hypericin, and pentagalloylglucose.

4. Discussion

ZIKV is one of the increasingly reported viral pathogens that sharply spreads beyond geographical borders. Since 2015, there are an increased number of microcephaly and Guillain-Barre syndrome incidences that have been observed worldwide. Specific vaccines or proper treatment option was not available to treat ZIKV till today. Few medications have been prescribed to release pain and fever related to the viral infection. Also, engineering methods were applied to identify potential peptides that penetrate the brain and cure infections. Apart from this, Ayurveda treatment is the most appropriate method to treat infections, which comprises only natural constituents and has no side effects. Hence, in our study, the molecular modeling protocols opted to identify natural phytochemicals as the potential inhibitors for ZIKV. Here, an active site-specific docking procedure was implemented and the active site residues were retrieved from the latest literature [22]. In the NS3 helicase (PDB: 5JRZ) structure, both RNA and ATP binding regions possess polar and hydrophobic characteristics and were suitable for high throughput screening; we utilized ATP binding sites for our study. A similar analysis conducted by Badshah et al. in 2019, which also utilized ATP binding site to carry out molecular docking for NS3 helicase against 1,4-benzothiazine derivatives to find out potential inhibitors for ZIKV [17]. According to “Lipinski’s rule of five,” the oral bioavailability of the natural constituents was assessed. This rule permits compounds with 2 violations. According to Benet et al.’s study in 2016, the rule of 5 is not applicable for natural compounds and natural compound derivatives [62]. To evaluate the nature of the compounds, the drug-likeness property prediction was carried out. Many computational studies report the use of screening techniques to identify potential inhibitors for various disease conditions. Kumari and Subbarao (2020) used FDA-approved drugs, natural products, and phytochemicals to screen against glutamine synthetase of Mycobacterium tuberculosis, an important enzyme of cell wall synthesis and nitrogen metabolism [63]. Sen et al. (2020) applied a structure-based virtual screening procedure to identify inhibitors for Severe Acute Respiratory Syndrome Coronaviruses 2 (SARS-CoV-2) [64]. They used phytochemicals retrieved from ginger, garlic, onion, peppermint, fenugreek, and chili to screen against the targets spike and main protease. Similarly, Gahlawat et al. (2020) screened inhibitors for the target main protease of SARS CoV-2; they used main protease inhibitors, natural products, and FDA-approved drugs [65]. The active site helps ligands to create enough contact sites to produce excellent interaction with the target protein by maintaining appropriate and desirable catalytic microenvironments [66–69]. In the ATP binding site of the NS3 helicase, the amino acids Glu231, Arg459, Gln455, Arg202, Glu286, Glu231, and Lys200 play a key role in the interaction with the ligands ellagic acid, hypericin, and pentagalloylglucose. From the docking results’ best binding energy compounds, ellagic acid (-9.5 kcal/mol), hypericin (-9.2 kcal/mol), and pentagalloylglucose (9.0 kcal/mol) were finally selected as the potential compounds. These docking results mimic the experimental study conducted by Acquadro et al. in 2020, which revealed the antiviral activity of leaf ethanolic extract of Punica granatum and ellagic acid against ZIKV [70]. The activity screening was carried out against the MR766 and HPF2013 strains of ZIKV. The study concluded that the compounds and the extracts were found to be active against the ZIKV; in particular, ellagic acid has shown EC₅₀ values of 30.86 μM and 46.23 μM for the strains MR766 and HPF2013, respectively. Another study revealed the antihuman rhinovirus (HRV) activity of ellagic acid isolated from the leaves of Lagerstroemia speciosa. The different strains of HRVs, namely, HRV-2, HRV-3, and HRV-4, have a 50% inhibitory concentration range of 38 μg/mL, 31 μg/mL, and 29 μg/mL, respectively, which was higher than that of the ribavirin drug. It also suggested that 50 μg/mL of ellagic acid inhibited viral replication by aiming at cellular components [31]. A randomized trial on ellagic acid and high-risk human papillomavirus- (HPV-) related low squamous intraepithelial lesion in women revealed a 74% HPV clearance in a group supplemented with 16 mg ellagic acid along with 100 mg of Annona muricata than the placebo group [71]. To find an inhibitor for Ebola virus cell entry, a high throughput assay was conducted with the extracts of 128 traditional medicine, which revealed that the crude extract and the compounds gallic acid and ellagic acid isolated from Rhodiola rosea were found to be effective against Ebola virus [72]. Ellagic acid showed favorable activity against many viruses. Hypericin, plant quinine, inactivates the human immunodeficiency virus type 1 (HIV-1), Sindbis virus, and murine cytomegalovirus [73]. Similarly, the pentagalloylglucose was also found to be effective against the influenza A virus, as it exhibited virus-induced hemagglutination in chicken red blood cells. Without affecting the nuclear transport of nucleoprotein or protein synthesis, pentagalloylglucose reduced the nucleoprotein accumulation in the plasma membrane during the replication cycle [44]. Jin et al.’s study in 2016 showcased the antiherpes simplex virus (HSV) activity of pentagalloylglucose, which slows down the process of nuclear transport in HSV-1 by blocking the upregulation of dynein [74]. Also, it affected the process of nucleocapsid egress in HSV-1. From the literature, it is evident that the plant compounds ellagic acid, hypericin, and pentagalloylglucose possess antiviral activity.

To understand the stability of the selected compounds, molecular dynamic studies were carried out for the complexes NS3 helicase-ellagic acid, NS3 helicase-hypericin, and NS3 helicase-pentagalloylglucose complexes, which revealed the higher stability of the complex; NS3 helicase-ellagic acid than the NS3 helicase-hypericin confirmed the nontoxicity and noncarcinogenicity of the compounds.

Medicinal plants serve as the most promising basis of natural antiviral substances against ZIKV. Bioactive substances isolated from medicinal plants have established their value against many life-threatening diseases like malaria, cancer, diabetes, heart disease, and Alzheimer’s diseases [75]. In the current study, a few antiviral plant compounds were screened against the NS3 helicase of ZIKV through molecular docking and molecular dynamic techniques and found ellagic acid as the more potential compound against ZIKV. Further, preclinical and clinical studies must be carried out to confirm the antiviral activity of ellagic acid against ZIKV.

5. Conclusion

The mosquito-transmitted ZIKV disease has now emerged as a public health threat worldwide due to its ability to cause neural infections. To identify potential inhibitors for ZIKV, we have screened medicinal plant-derived antiviral compounds against the nonstructural protein NS3 helicase of ZIKV. We conducted molecular docking and molecular dynamic studies and found ellagic acid as the most favorable binding energy compound against ZIKV NS3 helicase. Further clinical studies must be conducted on ellagic acid to ensure its antiviral activity against ZIKV.

Data Availability

The data that support the findings of this study are available from the corresponding author upon reasonable request.

Consent is not applicable.

Conflicts of Interest

The authors declare that there is no conflict of interest.

Authors’ Contributions

M.K., B.A.B., and S.K.M. were responsible for conceptualization; M.K., B.A.B., A.R., P.B., N.D.P., J.M.A., S.H.G., A.K.N., A.A., V.P.V., and S.K.M. were responsible for methodology; M.K. and S.K.M. were responsible for software; M.K., A.R., P.B., N.D.P., and S.K.M. were responsible for validation; M.K., B.A.B., A.R., P.B., N.D.P., J.M.A., S.H.G., A.K.N., A.A., V.P.V., and S.K.M. were responsible for data curation; M.K., B.A.B., A.R., P.B., N.D.P., J.M.A., S.H.G., A.K.N., A.A., V.P.V., and S.K.M. were responsible for writing—original draft preparation; M.K. and S.K.M. were responsible for writing—review and editing. All authors have read and agreed to the published version of the manuscript.

Acknowledgments

MK, JMA, and SKM would like to thank the management of PMCHRI for providing the necessary facilities to carry out this research work.

Supplementary Materials

Figure S1: Plant-derived antiviral compounds: (a) apigenin, (b) baicalein, (c) berberine, (d) betulin, (e) chebulagic acid, (f) curcumin, (g) ellagic acid, (h) epigallocatechin gallate, (i) fisetin, (j) geraniin, (k) glycyrrhizic acid, (l) hypericin, (m) hyperoside, (n) kaempferol, (o) lupeol, (p) mimusopic acid, (q) mulberroside C, (r) myricetin, (s) neoandrographolide, (t) pentagalloylglucose, (u) piperine, (v) quercetin, (w) rosmarinic acid, (x) rutin, and (y) torvoside. (Supplementary Materials)

References

D. Musso and D. J. Gubler, “Zika virus,” Clinical Microbiology Reviews, vol. 29, no. 3, pp. 487–524, 2016.
View at: Publisher Site | Google Scholar
R. Zhao, M. Wang, J. Cao et al., “Flavivirus: from structure to therapeutics development,” Life (Basel), vol. 11, no. 7, p. 615, 2021.
View at: Publisher Site | Google Scholar
S. S. Kazmi, W. Ali, N. Bibi, and F. Nouroz, “A review on Zika virus outbreak, epidemiology, transmission and infection dynamics,” Journal of Biological Research (Thessaloniki), vol. 27, no. 1, p. 5, 2020.
View at: Publisher Site | Google Scholar
D. Sirohi and R. J. Kuhn, “Zika virus structure, maturation, and receptors,” The Journal of Infectious Diseases, vol. 216, suppl_10, pp. S935–S944, 2017.
View at: Publisher Site | Google Scholar
K. Gorshkov, S. A. Shiryaev, S. Fertel et al., “Zika virus: origins, pathological action, and treatment strategies,” Frontiers in Microbiology, vol. 9, p. 3252, 2019.
View at: Publisher Site | Google Scholar
F. Noorbakhsh, K. Abdolmohammadi, Y. Fatahi et al., “Zika virus infection, basic and clinical aspects: a review article,” Iranian Journal of Public Health, vol. 48, no. 1, pp. 20–31, 2019.
View at: Google Scholar
A. Agrelli, R. R. de Moura, S. Crovella, and L. A. C. Brandão, “ZIKA virus entry mechanisms in human cells,” Infection, Genetics and Evolution, vol. 69, pp. 22–29, 2019.
View at: Publisher Site | Google Scholar
R. Hamel, O. Dejarnac, S. Wichit et al., “Biology of Zika virus infection in human skin cells,” Journal of Virology, vol. 89, no. 17, pp. 8880–8896, 2015.
View at: Publisher Site | Google Scholar
P. Pielnaa, M. Al-Saadawe, A. Saro et al., “Zika virus-spread, epidemiology, genome, transmission cycle, clinical manifestation, associated challenges, vaccine and antiviral drug development,” Virology, vol. 543, pp. 34–42, 2020.
View at: Publisher Site | Google Scholar
World Health Organization, https://www.who.int/news-room/fact-sheets/detail/zika-virus (10.01.2022).
P. J. Imperato, “The convergence of a virus, mosquitoes, and human travel in globalizing the Zika epidemic,” Journal of Community Health, vol. 41, no. 3, pp. 674–679, 2016.
View at: Publisher Site | Google Scholar
T. C. Pierson and M. S. Diamond, “The emergence of Zika virus and its new clinical syndromes,” Nature, vol. 560, no. 7720, pp. 573–581, 2018.
View at: Publisher Site | Google Scholar
B. Tesla, J. S. Powers, Y. Barnes, S. Lakhani, M. D. Acciani, and M. A. Brindley, “Temperate conditions limit Zika virus genome replication,” Journal of Virology, vol. 96, no. 10, article e0016522, 2022.
View at: Publisher Site | Google Scholar
D. E. Kim, J. S. Min, M. S. Jang et al., “Natural bis-benzylisoquinoline alkaloids-tetrandrine, fangchinoline, and cepharanthine, inhibit human coronavirus OC43 infection of MRC-5 human lung cells,” Biomolecules, vol. 9, no. 11, p. 696, 2019.
View at: Publisher Site | Google Scholar
S. Heeb, M. P. Fletcher, S. R. Chhabra, S. P. Diggle, P. Williams, and M. Cámara, “Quinolones: from antibiotics to autoinducers,” FEMS Microbiology Reviews, vol. 35, no. 2, pp. 247–274, 2011.
View at: Publisher Site | Google Scholar
D. Kumar, N. Sharma, M. Aarthy, S. K. Singh, and R. Giri, “Mechanistic insights into Zika virus NS3 helicase inhibition by epigallocatechin-3-gallate,” ACS Omega, vol. 5, no. 19, pp. 11217–11226, 2020.
View at: Publisher Site | Google Scholar
S. L. Badshah, N. Ahmad, A. Ur Rehman et al., “Molecular docking and simulation of Zika virus NS3 helicase,” BMC chemistry, vol. 13, no. 1, p. 67, 2019.
View at: Publisher Site | Google Scholar
H. Choudhry, F. A. Alzahrani, M. A. Hassan et al., “Zika virus targeting by screening inhibitors against NS2B/NS3 protease,” BioMed Research International, vol. 2019, 3947211 pages, 2019.
View at: Publisher Site | Google Scholar
L. R. Bowen, D. J. Li, D. T. Nola et al., “Identification of potential Zika virus NS2B-NS3 protease inhibitors via docking, molecular dynamics and consensus scoring-based virtual screening,” Journal of Molecular Modeling, vol. 25, no. 7, p. 194, 2019.
View at: Publisher Site | Google Scholar
H. Lee, J. Ren, S. Nocadello et al., “Identification of novel small molecule inhibitors against NS2B/NS3 serine protease from Zika virus,” Antiviral Research, vol. 139, pp. 49–58, 2017.
View at: Publisher Site | Google Scholar
S. M. Telehany, M. S. Humby, T. D. McGee Jr., S. P. Riley, A. Jacobs, and R. C. Rizzo, “Identification of Zika virus inhibitors using homology modeling and similarity-based screening to target glycoprotein E,” Biochemistry, vol. 59, no. 39, pp. 3709–3724, 2020.
View at: Publisher Site | Google Scholar
R. Jain, J. Coloma, A. García-Sastre, and A. K. Aggarwal, “Structure of the NS3 helicase from Zika virus,” Nature Structural & Molecular Biology, vol. 23, no. 8, pp. 752–754, 2016.
View at: Publisher Site | Google Scholar
S. K. Burley, H. M. Berman, G. J. Kleywegt, J. L. Markley, H. Nakamura, and S. Velankar, “Protein Data Bank (PDB): the single global macromolecular structure archive,” Methods in Molecular Biology, vol. 1607, pp. 627–641, 2017.
View at: Publisher Site | Google Scholar
B. Webb and A. Sali, “Protein structure modeling with MODELLER,” Methods in Molecular Biology, vol. 1137, pp. 1–15, 2014.
View at: Publisher Site | Google Scholar
Adapted from “Zika Virus”, by http://BioRender.com (2022). Retrieved from https://app.biorender.com/biorender-templates.
X. Lv, M. Qiu, D. Chen, N. Zheng, Y. Jin, and Z. Wu, “Apigenin inhibits enterovirus 71 replication through suppressing viral IRES activity and modulating cellular JNK pathway,” Antiviral Research, vol. 109, pp. 30–41, 2014.
View at: Publisher Site | Google Scholar
T. C. Hung, A. Jassey, C. H. Liu et al., “Berberine inhibits hepatitis C virus entry by targeting the viral E2 glycoprotein,” Phytomedicine, vol. 53, pp. 62–69, 2019.
View at: Publisher Site | Google Scholar
N. I. Pavlova, O. V. Savinova, S. N. Nikolaeva, E. I. Boreko, and O. B. Flekhter, “Antiviral activity of betulin, betulinic and betulonic acids against some enveloped and non-enveloped viruses,” Fitoterapia, vol. 74, no. 5, pp. 489–492, 2003.
View at: Publisher Site | Google Scholar
P. Li, R. Du, Y. Wang et al., “Identification of chebulinic acid and chebulagic acid as novel influenza viral neuraminidase inhibitors,” Frontiers in Microbiology, vol. 11, p. 182, 2020.
View at: Publisher Site | Google Scholar
R. K. Thimmulappa, M. N. Kumar, C. Shivamallu et al., “Antiviral and immunomodulatory activity of curcumin: a case for prophylactic therapy for COVID-19,” Heliyon, vol. 7, no. 2, article e06350, 2021.
View at: Publisher Site | Google Scholar
S. W. Park, M. J. Kwon, J. Y. Yoo, H. J. Choi, and Y. J. Ahn, “Antiviral activity and possible mode of action of ellagic acid identified in Lagerstroemia speciosa leaves toward human rhinoviruses,” BMC Complementary and Alternative Medicine, vol. 14, pp. 1–8, 2014.
View at: Google Scholar
Y. Zhong, C. M. Ma, and F. Shahidi, “Antioxidant and antiviral activities of lipophilic epigallocatechin gallate (EGCG) derivatives,” Journal of Functional Foods, vol. 4, no. 1, pp. 87–93, 2012.
View at: Publisher Site | Google Scholar
S. Y. Lyu, J. Y. Rhim, and W. B. Park, “Antiherpetic activities of flavonoids against herpes simplex virus type 1 (HSV-1) and type 2 (HSV-2) in vitro,” Archives of Pharmacal Research, vol. 28, no. 11, pp. 1293–1301, 2005.
View at: Publisher Site | Google Scholar
Y. Yang, L. Zhang, X. Fan, C. Qin, and J. Liu, “Antiviral effect of geraniin on human enterovirus 71 in vitro and in vivo,” Bioorganic & Medicinal Chemistry Letters, vol. 22, no. 6, pp. 2209–2211, 2012.
View at: Publisher Site | Google Scholar
Z. G. Sun, T. T. Zhao, N. Lu, Y. A. Yang, and H. L. Zhu, “Research progress of glycyrrhizic acid on antiviral activity,” Mini Reviews in Medicinal Chemistry, vol. 19, no. 10, pp. 826–832, 2019.
View at: Publisher Site | Google Scholar
P. Miskovsky, “Hypericin--a new antiviral and antitumor photosensitizer: mechanism of action and interaction with biological macromolecules,” Current Drug Targets, vol. 3, no. 1, pp. 55–84, 2002.
View at: Publisher Site | Google Scholar
P. Halder, U. Pal, P. Paladhi et al., “Evaluation of potency of the selected bioactive molecules from Indian medicinal plants with M^Pro of SARS-CoV-2 through in silico analysis,” J. Ayurveda Integr. Med., vol. 13, no. 2, article 100449, 2022.
View at: Publisher Site | Google Scholar
M. Behbahani, S. Sayedipour, A. Pourazar, and M. Shanehsazzadeh, “In vitro anti-HIV-1 activities of kaempferol and kaempferol-7-O-glucoside isolated from Securigera securidaca,” Res. Pharm. Sci., vol. 9, no. 6, pp. 463–469, 2014.
View at: Google Scholar
F. M. Tolo, G. W. Rukunga, F. W. Muli et al., “The antiviral activity of compounds isolated from Kenyan Carissa edulis (Forssk.) Vahl,” Journal of Medicinal Plant Research: Planta Medica, vol. 4, pp. 1517–1522, 2010.
View at: Google Scholar
N. P. Sahu, K. Koike, Z. Jia, and T. Nikaido, “Triterpenoid saponins from _Mimusops elengi_,” Phytochemistry, vol. 44, no. 6, pp. 1145–1149, 1997.
View at: Publisher Site | Google Scholar
Y. Cao, E. Lei, L. Li et al., “Antiviral activity of Mulberroside C against enterovirus A71 _in vitro_ and _in vivo_,” European Journal of Pharmacology, vol. 906, article 174204, 2021.
View at: Publisher Site | Google Scholar
S. Pasetto, V. Pardi, and R. M. Murata, “Anti-HIV-1 activity of flavonoid myricetin on HIV-1 infection in a dual-chamber in vitro model,” PLoS One, vol. 9, no. 12, article e115323, 2014.
View at: Publisher Site | Google Scholar
A. K. Jadhav and S. M. Karuppayil, “Andrographis paniculata (Burm. F) Wall ex Nees: antiviral properties,” Phytotherapy Research, vol. 35, no. 10, pp. 5365–5373, 2021.
View at: Publisher Site | Google Scholar
G. Liu, S. Xiong, Y. F. Xiang et al., “Antiviral activity and possible mechanisms of action of pentagalloylglucose (PGG) against influenza A virus,” Archives of Virology, vol. 156, no. 8, pp. 1359–1369, 2011.
View at: Publisher Site | Google Scholar
A. Nag and R. R. Chowdhury, “Piperine, an alkaloid of black pepper seeds can effectively inhibit the antiviral enzymes of dengue and Ebola viruses, an in silico molecular docking study,” Virus, vol. 31, no. 3, pp. 308–315, 2020.
View at: Publisher Site | Google Scholar
K. Zandi, B. T. Teoh, S. S. Sam, P. F. Wong, M. R. Mustafa, and S. AbuBakar, “Antiviral activity of four types of bioflavonoid against dengue virus type-2,” Virology Journal, vol. 8, pp. 1–1, 2011.
View at: Google Scholar
Y. Tsukamoto, S. Ikeda, K. Uwai et al., “Rosmarinic acid is a novel inhibitor for hepatitis B virus replication targeting viral epsilon RNA-polymerase interaction,” PLoS One, vol. 13, no. 5, article e0197664, 2018.
View at: Publisher Site | Google Scholar
D. D. Orhan, B. Ozçelik, S. Ozgen, and F. Ergun, “Antibacterial, antifungal, and antiviral activities of some flavonoids,” Microbiological Research, vol. 165, no. 6, pp. 496–504, 2010.
View at: Publisher Site | Google Scholar
A. K. Ibrahim, A. I. Youssef, A. S. Arafa, and S. A. Ahmed, “Anti-H5N1 virus flavonoids fromCapparis sinaicaVeill,” Natural Product Research, vol. 27, no. 22, pp. 2149–2153, 2013.
View at: Publisher Site | Google Scholar
S. Kim, P. A. Thiessen, E. E. Bolton et al., “PubChem substance and compound databases,” Nucleic Acids Research, vol. 44, no. D1, pp. D1202–D1213, 2016.
View at: Publisher Site | Google Scholar
C. A. Lipinski, “Lead- and drug-like compounds: the rule-of-five revolution,” Drug Discovery Today: Technologies, vol. 1, no. 4, pp. 337–341, 2004.
View at: Publisher Site | Google Scholar
T. B. Hadda, V. Rastija, F. AlMalki et al., “Petra/Osiris/Molinspiration and molecular docking analyses of 3-hydroxy-indolin-2-one derivatives as potential antiviral agents,” Current Computer-Aided Drug Design, vol. 17, no. 1, pp. 123–133, 2021.
View at: Publisher Site | Google Scholar
O. Trott and A. J. Olson, “AutoDock Vina: improving the speed and accuracy of docking with a new scoring function, efficient optimization, and multithreading,” Journal of Computational Chemistry, vol. 31, no. 2, pp. 455–461, 2010.
View at: Publisher Site | Google Scholar
M. Rostkowski, M. H. Olssonm, C. R. Sondergaard, and J. H. Jensen, “Graphical analysis of pH-dependent properties of proteins predicted using PROPKA,” BMC Structural Biology, vol. 11, no. 1, p. 6, 2011.
View at: Publisher Site | Google Scholar
N. Guex and M. C. Peitsch, “SWISS-MODEL and the Swiss-Pdb Viewer: an environment for comparative protein modeling,” Electrophoresis, vol. 18, no. 15, pp. 2714–2723, 1997.
View at: Publisher Site | Google Scholar
T. Schwede, J. Kopp, N. Guex, and M. C. Peitsch, “SWISS-MODEL: an automated protein homology-modeling server,” Nucleic Acids Research, vol. 31, no. 13, pp. 3381–3385, 2003.
View at: Publisher Site | Google Scholar
T. F. Vieira and S. F. Sousa, “Comparing AutoDock and Vina in ligand/decoy discrimination for virtual screening,” Applied Sciences, vol. 9, no. 21, p. 4538, 2019.
View at: Publisher Site | Google Scholar
M. J. Abraham, T. Murtola, R. Schulz et al., “GROMACS: high performance molecular simulations through multi-level parallelism from laptops to supercomputers,” SoftwareX, vol. 1-2, pp. 19–25, 2015.
View at: Publisher Site | Google Scholar
K. Vanommeslaeghe, E. Hatcher, C. Acharya et al., “CHARMM general force field: a force field for drug-like molecules compatible with the CHARMM all-atom additive biological force fields,” Journal of Computational Chemistry, vol. 31, no. 4, pp. 671–690, 2010.
View at: Publisher Site | Google Scholar
B. Hess, “P-LINCS: a parallel linear constraint solver for molecular simulation,” Journal of Chemical Theory and Computation, vol. 4, no. 1, pp. 116–122, 2008.
View at: Publisher Site | Google Scholar
H. Yang, C. Lou, L. Sun et al., “admetSAR 2.0: web-service for prediction and optimization of chemical ADMET properties,” Bioinformatics, vol. 35, no. 6, pp. 1067–1069, 2019.
View at: Publisher Site | Google Scholar
L. Z. Benet, C. M. Hosey, O. Ursu, and T. I. Oprea, “BDDCS, the rule of 5 and drugability,” Advanced Drug Delivery Reviews, vol. 101, pp. 89–98, 2016.
View at: Publisher Site | Google Scholar
M. Kumari and N. Subbarao, “Virtual screening to identify novel potential inhibitors for glutamine synthetase ofMycobacterium tuberculosis,” Journal of Biomolecular Structure & Dynamics, vol. 38, no. 17, pp. 5062–5080, 2020.
View at: Publisher Site | Google Scholar
D. Sen, P. Debnath, B. Debnath, S. Bhaumik, and S. Debnath, “Identification of potential inhibitors of SARS-CoV-2 main protease and spike receptor from 10 important spices through structure-based virtual screening and molecular dynamic study,” Journal of Biomolecular Structure & Dynamics, vol. 40, pp. 1–22, 2020.
View at: Google Scholar
A. Gahlawat, N. Kumar, R. Kumar et al., “Structure-based virtual screening to discover potential lead molecules for the SARS-CoV-2 main protease,” Journal of Chemical Information and Modeling, vol. 60, no. 12, pp. 5781–5793, 2020.
View at: Publisher Site | Google Scholar
M. Tasleem, A. Alrehaily, T. M. Almeleebia et al., “Investigation of antidepressant properties of yohimbine by employing structure-based computational assessments,” Current Issues in Molecular Biology, vol. 43, no. 3, pp. 1805–1827, 2021.
View at: Publisher Site | Google Scholar
M. Y. Alshahrani, K. M. Alshahrani, M. Tasleem et al., “Computational screening of natural compounds for identification of potential anti-cancer agents targeting MCM7 protein,” Molecules, vol. 26, no. 19, p. 5878, 2021.
View at: Publisher Site | Google Scholar
M. Tasleem, R. Ishrat, A. Islam, F. Ahmad, and I. Hassan, “Structural characterization, homology modeling and docking studies of ARG674 mutation in MyH8 gene associated with trismus-pseudocamptodactyly syndrome,” Letters in Drug Design & Discovery, vol. 11, no. 10, pp. 1177–1187, 2014.
View at: Publisher Site | Google Scholar
M. Akhter, M. Tasleem, M. M. Alam, and S. Ali, “_In silico_ approach for bioremediation of arsenic by structure prediction and docking studies of arsenite oxidase from _Pseudomonas stutzeri TS44_,” International Biodeterioration & Biodegradation, vol. 122, pp. 82–91, 2017.
View at: Publisher Site | Google Scholar
S. Acquadro, A. Civra, C. Cagliero et al., “Punica granatum leaf ethanolic extract and ellagic acid as inhibitors of Zika virus infection,” Planta Medica, vol. 86, no. 18, pp. 1363–1374, 2020.
View at: Publisher Site | Google Scholar
M. Le Donne, M. Lentini, A. Alibrandi et al., “Antiviral activity of ellagic acid and Annona muricata in cervical HPV related pre-neoplastic lesions: a randomized trial,” Journal of Functional Foods, vol. 35, pp. 549–554, 2017.
View at: Publisher Site | Google Scholar
Q. Cui, R. Du, M. Anantpadma et al., “Identification of ellagic acid from plant Rhodiola rosea L. as an anti-Ebola virus entry inhibitor,” Viruses, vol. 10, no. 4, p. 152, 2018.
View at: Publisher Site | Google Scholar
J. B. Hudson, I. Lopez-Bazzocchi, and G. H. Towers, “Antiviral activities of hypericin,” Antiviral Research, vol. 15, no. 2, pp. 101–112, 1991.
View at: Publisher Site | Google Scholar
F. Jin, K. Ma, M. Chen et al., “Pentagalloylglucose blocks the nuclear transport and the process of nucleocapsid egress to inhibit HSV-1 infection,” Japanese Journal of Infectious Diseases, vol. 69, no. 2, pp. 135–142, 2016.
View at: Publisher Site | Google Scholar
S. R. Ahmed, A. Banik, S. M. Anni, and M. M. H. Chowdhury, “Inhibitory potential of plant-derived metabolites against Zika virus: a computational-aided approach,” Phytomedicine Plus, vol. 1, no. 4, article 100129, 2021.
View at: Publisher Site | Google Scholar

Copyright

Copyright © 2022 Malathi Kullappan et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

PDF Download Citation

Download other formats

Order printed copies

Views

365

Downloads

408

Citations