Next Article in Journal
A Comparison of Germany and the United Kingdom Indicates That More SARS-CoV-2 Circulation and Less Restrictions in the Warm Season Might Reduce Overall COVID-19 Burden
Next Article in Special Issue
Dithymoquinone Analogues as Potential Candidate(s) for Neurological Manifestation Associated with COVID-19: A Therapeutic Strategy for Neuro-COVID
Previous Article in Journal
Mineral Element Insiders and Outliers Play Crucial Roles in Biological Evolution
Previous Article in Special Issue
2-(3-Bromophenyl)-8-fluoroquinazoline-4-carboxylic Acid as a Novel and Selective Aurora A Kinase Inhibitory Lead with Apoptosis Properties: Design, Synthesis, In Vitro and In Silico Biological Evaluation
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
Life
Volume 12
Issue 7
10.3390/life12070952
life-logo
Submit to this Journal Review for this Journal Propose a Special Issue
► Article Menu

Article Menu

share Share announcement Help format_quote Cite question_answer Discuss in SciProfiles thumb_up
...
Endorse textsms
...
Comment
first_page
settings
Order Article Reprints
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Discovery of Small Molecules from Echinacea angustifolia Targeting RNA-Dependent RNA Polymerase of Japanese Encephalitis Virus

by
Pardeep Yadav
1,2,
Sherif A. El-Kafrawy
3,4,
Mai M. El-Day
3,4,
Wejdan T. Alghafari
3,5,
Arwa A. Faizo
3,4,
Saurabh Kumar Jha
1,6,7,
Vivek Dhar Dwivedi
2,8,* and
Esam I. Azhar
3,4,*
1
Department of Biotechnology, School of Engineering & Technology, Sharda University, Greater Noida 201310, India
2
Center for Bioinformatics, Computational and Systems Biology, Pathfinder Research and Training Foundation, Greater Noida 201308, India
3
Special Infectious Agents Unit-BSL-3, King Fahd Medical Research Center, King Abdulaziz University, Jeddah 21362, Saudi Arabia
4
Department of Medical Laboratory Sciences, Faculty of Applied Medical Sciences, King Abdulaziz University, Jeddah 21362, Saudi Arabia
5
Clinical Nutrition Department, Faculty of Applied Medical Sciences, King Abdulaziz University, Jeddah 21362, Saudi Arabia
6
Department of Biotechnology Engineering and Food Technology, Chandigarh University, Mohali 140413, India
7
Department of Biotechnology, School of Applied & Life Sciences (SALS), Uttaranchal University, Dehradun 248007, India
8
Institute of Advanced Materials, IAAM, 59053 Ulrika, Sweden
*
Authors to whom correspondence should be addressed.
Life 2022, 12(7), 952; https://doi.org/10.3390/life12070952
Submission received: 20 May 2022 / Revised: 13 June 2022 / Accepted: 20 June 2022 / Published: 24 June 2022
(This article belongs to the Special Issue New Trends in Pharmaceutical Science)
Browse Figures
Versions Notes

Abstract

:
The Japanese encephalitis virus (JEV), a mosquito-borne flavivirus that causes viral encephalitis leading to neural damage, is a major threat in most Asian countries. The RNA-dependent RNA polymerase (RdRp) present in the viral genome is the key component for genome replication, making it an attractive target for antiviral drug development. In this study, the natural products from Echinacea angustifolia were retrieved for structure-based virtual screening against JEV–RdRp. The top six compounds (Echinacoside, Echinacin, Rutin, Cynaroside, Quercetagetin 7-glucoside, and Kaempferol-3-glucoside) were obtained based on the highest negative docking score, ADMET (absorption, distribution, metabolism, excretion, and toxicity), and molecular interaction. The computational analysis of these selected compounds against the co-crystallized ligands, i.e., ATP and GTP, were performed. Further, 100 ns molecular dynamic simulation and post-free binding energy calculation of all the selected compounds complexed with JEV–RdRP were performed to check the stability of the complexes. The obtained results showed considerable stability and intermolecular interaction with native ligand-binding site residues of JEV–RdRp. Hence, selected natural compounds are admissible inhibitors of JEV–RdRp protein and can be considered for future antiviral drug development studies.

1. Introduction

Japanese encephalitis (JE) is the most commonly diagnosed epidemic encephalitis—a vector-borne acute central nervous system infection, caused by the Japanese encephalitis virus (JEV) that belongs to the genus Flavivirus of the family Flaviviridae and comprises five genotypes (GI–GV) [1,2,3]. JEV is transmitted to humans by the bite of infected mosquitoes (Culex spp.) and causes severe neurological manifestations [4,5]. JEV is an endemic in eastern and southern Asia countries and Australia, with fatality rates of 20–30% and life-long neurological impairments and sequelae among one-half of the survivors [6,7,8,9,10]. Nearly 68,000 clinical cases of JEV infection and approximately 20,000 deaths are reported annually in many Asian countries, and its transmission is endemic in 24 countries of Southeast Asia and Western Pacific regions [11,12,13]. According to World Health Organization (WHO), more than 3 billion people are at risk of JEV infection globally [11,12,13].
JEV is characterized as a positive-sense single-stranded RNA (+ssRNA) genome encapsulated by an icosahedral-shaped viral capsid. The RNA genome of JEV (≈11 kb) is modified at the 5′ end cap with an m7GpppAm structure and mimics cellular mRNAs, except for lacking the polyadenylated tail [14]. The viral genome encodes a single 370 kDa polyprotein, which is subsequently cleaved by a combination of viral and host proteases to yield three structural (capsid, membrane, and envelope) and seven nonstructural (NS: NS1, NS2A, NS2B, NS3, NS4A, NS4B, and NS5) proteins [14,15]. Among the seven NS proteins, NS5 is the largest (~105 kDa) and conserved viral protein—housing a C-terminal RNA-dependent RNA polymerase (RdRp) and an N-terminal guanylyltransferase (GTase) and methyltransferase (MTase) activity—and is associated with other viral proteins, host factors, and viral RNA to form the replication complex (RCs) [16,17]. The JEV NS5–RdRp domain (JEV–RdRp) replicates the genomic RNA into the uncapped minus-strand RNA, which is subsequently used as a template to produce a large excess of the positive-sense viral RNA genome [18]. The JEV NS5 MTase (JEVMTase) domain then performs both guanine N-7 and ribose 2′-OH methylations for the capping of the newly synthesized positive-sense RNAs [19,20]. Moreover, the JEV–RdRp domain is involved in the synthesis of both negative- and positive-strand RNAs during viral genome replication and translation [21]. As the JEV–RdRp domain is required for RNA synthesis in viral replication and lacks a similar protein in host cells, it has been demonstrated as a promising target for antiviral drug development against the JEV infection [14,22].
To tackle the JEV infection, both inactivated and live-attenuated vaccines have been progressively designed and utilized to prevent the JEV infection [23,24]. Unfortunately, regardless of a licensed vaccine to prevent JEV infection, infections strike annually due to a denial of complete coverage and access [25,26,27,28,29]. Several other efforts were also made in laboratories globally to identify an effective and safe therapeutic agent demonstrating an anti-JEV activity, including N-nonyl-deoxynojirimycin [30], dehydroepiandrosterone [31], N-methylisatin-beta-thiosemicarbazone derivative (SCH 16) [32], indirubin [33], manidipine [34], chlorpromazine [35], etanercept [36], minocycline [37], ouabain, and digoxin [38], but none of them is yet approved as antiviral against JEV infection treatment. The lack of efficient and approved antiviral agents in the prevention and treatment of the JEV infection, therefore, demands a potential antiviral to reduce JEV-infection-associated morbidity and mortality rates and continue to increase in global distribution.
In drug discovery pipelines, medicinal plants have become the popular source of natural compounds with high structural diversity against viral infections [39,40,41,42,43]. Several medicinal plants or natural compounds have shown promising effects against JEV [33,44,45,46,47]. Echinacea angustifolia (E. angustifolia) is one of the most well-known medicinal plants in its genus, which has been studied for medicinal properties, including antibacterial, anti-inflammatory, and antiviral activity properties [48,49,50,51,52,53]. However, to the best of our knowledge, the phytoconstituents of E. angustifolia are not yet demonstrated for the antiviral activity against JEV. Therefore, the present study was designed to explore the potential of natural compounds as therapeutic against the JEV infection via inhibition of the JEV–RdRp domain using molecular simulations, including structure-based virtual screening, molecular docking, and molecular dynamics (MD) simulations, as well as free energy calculations.

2. Methodology

2.1. Receptor and Ligand Collection

The three-dimensional (3D) crystal structure of the JEV–RdRp complexed with guanosine-5′-triphosphate (GTP) (PDB ID: 4HDG) [22] as a receptor was downloaded from the PDB database (https://www.rcsb.org/) [54]. Additionally, a total of 53 natural compounds, reported in Echinacea angustifolia, were collected from the literature, and their 3D atomic coordinates were retrieved from the PubChem database (https://pubchem.ncbi.nlm.nih.gov/) [55] to form the ligand library for further computational analysis.

2.2. Virtual Screening, Redocking Analysis, and ADME

For the structure-based virtual screening, the JEV–RdRp domain was processed and minimized using the Dock Prep tool in UCSF Chimera [56], with default parameters. Here, the receptor molecule was prepared by removing the native ligand, water molecules, heteroatoms, and non-polar hydrogen atoms, along with the addition of Gasteiger charges and polar hydrogen atoms, followed by structure minimization, using 100 and 50 steps of steepest descent and conjugate gradient methods, respectively. Furthermore, for the ligand library preparation, i.e., the collected natural compounds of Echinacea angustifolia and structure-based virtual screening, PyRx0.8 (Virtual Screening Tool) [57] was used, under default parameters. Herein, each ligand was initially added for Gasteiger partial charge, and energy was minimized under MMFF94 force field using default steepest descent and conjugate gradient with 1000 steps, followed by the addition of polar hydrogens, as reported earlier [58]. The prepared ligands were docked in the binding pocket of JEV–RdRp covering essential residues (Lys463, Lys471, Arg474, Asp541, Ser604, Asp668, Arg734, Arg742, Ser799, Try800, and Ser801) in the docking grid 39.3513 Å × 29.2100 Å × 28.2521 Å and center at −37.5873 Å, −5.1747 Å, −35.1844 Å along the X, Y, and Z axes.
Subsequently, the top 10 ranked ligands based on binding scores were selected for redocking analysis by comparison to the ATP and GTP as reference molecules using the Chimera-AutoDock Vina Plugin setup [56,59]. Herein, the processed protein structure was docked with the selected compounds in the same docking pocket as considered for virtual screening (center coordinates: −37.5873 Å, −5.1747 Å, −35.1844 Å; grid size: 39.3513 Å × 29.2100 Å × 28.2521 Å) under default parameters. Then, at least 10 docked poses were generated for each ligand, and the docked conformations with the highest negative docking scores and least root-mean-square deviation (RMSD) values (by default 0 in AutoDock Vina) were extracted for binding pose and intermolecular interaction analysis, under default parameters, in Maestro tool of Schrödinger suite 2018-4 [60]. All 3D and 2D figures were generated in Maestro tool of Schrödinger suite2018-4 [60].
SwissADME online server (http://www.swissadme.ch) was used to understand the properties related to the pharmacokinetics and drug likeliness of the selected compounds such as absorption, distribution, metabolism, and excretion (ADME) [61]. The bioavailability is another important parameter that makes any drug a promising therapeutic, and therefore, several parameters were observed during the ADME analysis such as blood–brain barrier (BBB), permeability, and cytochrome inhibition activities.

2.3. Molecular Dynamics Simulation

The molecular dynamics (MD) simulations for the selected best poses of the top six protein–ligand complexes were performed using with academic Maestro-Desmond v5.6 suite [62,63]. Initially, the docked complex was placed in the center of the orthorhombic box (10 Å  ×  10 Å  ×  10 Å) and amended with an explicit TIP4Pwater model using the system builder module. Thereafter, the complete system was neutralized using the counter sodium and chloride ions, while being aced at a distance of 20 Å around the docked ligand in the binding pocket of protein, and then minimized under default parameters using the minimization tool. Next, the whole system was simulated using an NPT ensemble, maintained by the Nose–Hoover thermostat and the Martyna–Tobias–Klein barostat method [64], with the temperature set at 300 K and pressure set at 1.013 bar under default parameters. For the simulations, the cutoff radius in Coulomb interactions was set at 9.0 Å, and the long-range electrostatic interactions were computed via the particle mesh Ewald method [65]. Each complex was simulated for 100 ns under optimized potentials for liquid simulations (OPLS)-2005 force-field parameters, and a total of 5000 frames were saved for analysis. The generated JEV–RdRp trajectories with selected ligands were analyzed for statistical parameters, including root-mean-square deviation (RMSD) and root-mean-square fluctuation (RMSF) analysis, and intermolecular interactions (protein–ligand contact mapping) as a function of 100 ns using the Simulation Interaction Diagram (SID) tool implemented in free academic Desmond module with Maestro-Schrödinger suite 2018–4 interface [62,63].

2.4. Post-Simulation Analysis

2.4.1. Essential Dynamics

The essential dynamics, defined in terms of principal component analysis (PCA), to study the dynamic motion of the protein was assessed in the presence of docked ligands from the respective MD simulation trajectories using the Bio3d package [66]. For this purpose, a total of 5000 frames of the Cα atoms from the respective MD simulation trajectories were superimposed on the initial pose to reduce the root-mean-square variations between similar residues by utilizing the fit.xyz function. These superimposed poses were then analyzed to generate respective plots of PCA components from each simulation trajectory using the pca.xyz function of the plot (pc) in the Bio3d package [66].

2.4.2. Prime MM/GBSA Binding Free Energy Calculations

To calculate the binding free energy and ligand strain energy for the docked poses, the endpoint MMGBSA calculation was applied to the complete 100 ns MD simulation trajectories using the thermal_mmgbsa.py python script in the Prime MMGBSA module of the Schrödinger suite (Schrödinger Release 2020-4: Prime, Schrödinger, LLC, New York, NY, USA, 2020), where each snapshot was treated for the removal of solvent and ions, as well as split into individual protein and ligand conformation for the free energy calculation. The net free binding energy (ΔG) was calculated using the following equation:
ΔGBind = ΔGcomplex(minimized) − (ΔGreceptor(minimized) + ΔGligand(minimized))
where ΔGBind denotes the binding free energy; ∆Gcomplex indicates the binding free energy of the complex; ∆GReceptor and ∆GLigand exhibit the energy for receptor and ligand, respectively. Finally, the computed binding free energy, along with energy dissociation components, was provided as average with the standard deviation for each simulation trajectory.

3. Results and Discussion

3.1. Structure-Based Virtual Screening and ADME

The natural products reported in the E. angustifolia were collected from the reported literature, and the respective 3D structures were collected retrieved from the PubChem database. The collected compounds were used as ligand libraries for the structure-based virtual screening against the nucleotide GTP-binding pocket in the crystal structure of JEV–RdRp (PDB ID: 4HDG). The top conformations of the screened ligands exhibited binding affinity ranging from −11.1 to −1.9 kcal/mol, targeting the binding pocket of the selected JEV–RdRp protein (Table S1). Thus, the top docked poses of the first six compounds—namely, Echinacoside, Echinacin, Rutin, Cynaroside, Quercetagetin 7-glucoside, and Kaempferol-3-glucoside—with significant binding scores (>−9 kcal/mol) were considered for the ADME/Tox and redocking analysis to find drug-likeness properties and the most ideal docking conformation, respectively, in the selected binding pocket of the viral protein.
Drug likeness and pharmacological characters are considered to be important factors for understanding the medicinal application of the compounds. Therefore, the ADME analysis of the top six compounds, i.e., Echinacoside, Echinacin, Rutin, Cynaroside, Quercetagetin 7-glucoside, and Kaempferol-3-glucoside, were predicted using SwissADME, resulting in the properties related to the pharmacokinetics and toxicity (Table S2). All of the selected compounds were found to be non-inhibitors of Cytochrome P450 2D6 (CYP2D6), as the inhibition of CYP2D6 leads to drug–drug interactions. Additionally, each of these natural compounds was found to be impermeable to the Blood–brain barrier (BBB) (Table S2). Based on these findings, as well as other properties such as pharmacokinetics and drug-likeness, the selected compounds have considerable medicinal properties.

3.2. Redocking and Molecular Contact Analysis

From the virtual screening, six potential natural compounds—Echinacoside, Echinacin, Rutin, Cynaroside, Quercetagetin 7-glucoside, and Kaempferol-3-glucoside—were selected for redocking in the selected binding pocket of JEV–RdRp against ATP and GTP as reference compounds using AutoDock Vina. Afterward, the docked poses with the highest negative docking energy values corresponding to zero RMSD values for each natural compound were considered for further computational analysis (Figure 1). In this study, among the selected natural compounds, Echinacoside docked with JEV–RdRp showed maximum docking energy (−11.1 kcal/mol), while JEV–RdRp–Kaempferol-3-glucoside docked complex was noted for lowest docking scores (10 kcal/mol) by comparison to reference compounds ATP (8.6 kcal/mol) and GTP (−9.0 kcal/mol). Additionally, a substantial number of intermolecular interactions, including hydrogen bond formation, π–π stacking, π–cation, hydrophobic, polar, negative, positive, glycine, and salt bridge interactions were also noted in the docked complexes by comparison to the reference docked complexes (Table 1, Figure S1). Therefore, the calculated docking scores and intermolecular contact profiling between the docked natural compounds and JEV–RdRp indicates the stability of the respective docked complexes.

3.3. Molecular Dynamics Simulation Analysis

A molecular dynamics simulation is used to predict the stability of protein–ligand complexes with respect to simulation intervals. In this study, the molecular dynamics simulation was analyzed in terms of the last pose from the simulation trajectory, protein and protein-fit-ligand root-mean-square deviation (RMSD), protein and protein-fit-ligand root-mean-square fluctuation (RMSF), and protein–ligand contact mapping as a function of 100 ns interval. Figure 2 shows the last poses of the MD trajectory were extracted and compared with the initially docked poses to monitor the relative occupation of the docked ligands in the protein structure. Notably, all of the docked natural compounds showed relatively substantial residence in the selective pocket of the JEV–RdRp as a function of 100 ns MD simulation interval, except for an acceptable deviation in the ligand conformations, which was noted against reference compounds (Figure S2). These results suggested the docked natural compounds as substantial inhibitors by comparison to the reference compounds, i.e., ATP and GTP.
Furthermore, the calculated RMSD for the protein in all of the docked complexes with natural compounds indicated acceptable deviations (<2 Å) with respect to 100 ns MD simulation by comparison to reference compounds (Figure 3). These observations were also supported by the calculated RMSF values (Figure S3) for the protein structures docked with respective compounds, except for occasional residual higher fluctuations (<5 Å) that were also noted in the terminal regions and residues forming direct contact or adjacent residues to the interacting residues with the docked ligands. Moreover, the calculated protein-fit-ligand RMSD values showed considerable deviations (<4.5 Å) (Figure 3) throughout the 100 ns simulation interval against reference compounds, i.e., ATP (<4.0 Å) (Figure 3g) and GTP (<4.0 Å) (Figure 3h). Notably, among the selected natural compounds, Rutin, Cynaroside, and Kaempferol-3-glucoside showed the most substantial stability and equilibrium in protein-fit-ligand values as a function of 100 ns against reference compounds, i.e., ATP and GTP. Moreover, the calculated RMSF values also showed variation within 2 Å (Figure S4) during the 100 ns, supporting the observed stability of the docked ligands with the JEV–RdRp during MD simulation.
Moreover, the protein–ligand contacts were also extracted from the 100 ns MD simulation trajectories for all of the docked complexes. Specifically, each type of intermolecular interaction, hydrogen bonding, hydrophobic interactions, water bridge formation, and ionic interaction was extracted and plotted as a fraction of the total molecular contact formed during the 100 ns interval for all the simulated complexes.
In the case of the Echinacoside–JEV–RdRp docked complex, it exhibited hydrogen bond formation for more than 100% of simulation time in Ser666 and Asp669 residues, while Trp800 and Ile802 residues showed 50% and 35% of total simulation time for the hydrophobic interaction with the docked ligand; Gly472, Cys714, and Arg474 residues participated in water bridges formation, in which Arg474 residue showed more than 100% of water bridge formation during the simulation time in addition to a significant H-bond. Additionally, Arg734 and Arg742 residues formed an ionic bond for more than 20% simulation time (Figure 4a). Likewise, the Echinacin–JEV–RdRp complex exhibited hydrogen bond interaction in Asp541 and Lys463 residue for 90% and 70% of total simulation time; Trp800 and Arg474 residues were also noted to form hydrophobic interaction for 100% of simulation time. Gly605 and Arg460 residues showed more than 80% of the total interaction fraction in water bridge formation during the simulation period (Figure 4b). Additionally, protein–ligand contact analysis of Rutin–JEV–RdRp showed a high contribution of Asp668, Asp669, and Ser604 in hydrogen bond formation (100%) during the total simulation interval. Tyr610 residue participated in hydrophobic interaction for more than 50% of the simulation time. A water bridge was formed by Gly412, Trp477, and Arg484 for more than 50% of simulation time and by Gly412 for 100% of simulation time (Figure 4c). Furthermore, the Cynaroside–JEV–RdRp complex exhibited a hydrogen bond with Glu510 and Asp668 residues for 100% of the simulation time. Tyr610 also showed the hydrophobic interaction with the docked ligand at a 100% of simulation interval. Asp541 and Asp669 residues were involved in water bridge formation for more than 50%, along with the hydrogen bond (Figure 4d). In comparison, in the Quercetagetin 7-glucoside–JEV–RdRp complex, Glu510 residue showed more than 100% of the total interaction fraction in hydrogen bond formation during the simulation interval, and Gly605 and Asp668 residues participated in hydrogen bond formation for more than 50% of simulation time; Asp668 also participated in water bridge formation (50%) during the simulation interval. Ile802 residue was noted for hydrophobic interaction with the docked ligand for more than 50% of the simulation time (Figure 4e). In the protein contact mapping analysis of the Kaempferol-3-glucoside–JEV–RdRp complex, Asp541 exhibited more than 100% of hydrogen bond formation during the 100 ns simulation time, whereas Tyr610, Trp800, and Ile802 residues showed 50% of hydrophobic interaction during the total simulation period. Notably, Asp668 residue showed more than 50% of the water bridge interaction along with the hydrogen bond (Figure 4f). In contrast, the protein–ligand mapping of the JEV–RdRp with its native ligand ATP exhibited more hydrogen bond formation for 100% of simulation time with Arg460, Lys463, Lys471, Arg474, and Arg734 residues, along with water bridge formation (50% of the interaction fraction). Ile802 and Lys471 residues participated in hydrophobic (50%) and ionic interaction (30%) during the 100 ns simulation period (Figure 4g). When the protein–ligand mapping of the JEV–RdRp with its native ligand GTP was analyzed, Arg460, Arg474, and Arg742 participated in hydrogen bond formation for more than 100%, in addition to the water bridge interaction of the total simulation time. Additionally, Lys463 and Arg734 exhibited hydrophobic and ionic bonds for more than 90% and more than 50% of the total interaction fraction (Figure 4h). The solvent accessible surface area (SASA) of the selected natural compound and the reference molecules were also studied to track the protein’s flexibility, stability, and folding in the presence and absence of ligands (Figure S5) [43].
Collectively, the analysis of the MD simulations supported the selected natural compounds as inhibitors of JEV–RdRp against the reference compounds, i.e., ATP and GTP. Hence, based on the statistical analysis and intermolecular interaction profiling of the docked natural compounds in the selected pocket of JEV–RdRp, the potential ligands can be placed in the order of Rutin, Cynaroside, Kaempferol-3-glucoside, Echinacoside, Quercetagetin 7-glucoside, and Echinacin, exhibiting the most stability with the JEV–RdRp protein.

3.4. Post-Simulation Analysis

3.4.1. Principal Component Analysis

Studying the essential dynamics based on principal component analysis provides a method to collect the domain dynamics and displacement of atoms in the protein, required for biological function. Figure 5 shows the percentage of variance (%) for the calculated mean square position variations in the covariance matrix as a function of extracted 20 eigenmodes from the respective MD simulation trajectories of the docked complexes. Notably, all of the JEV–RdRp structures docked with natural compounds relatively showed similar drops in eigen fraction values by comparison to the protein structures docked with reference compounds GTP but not ATP (Figure 5 and Figure S6). However, no significant changes in the eigen fraction were observed from 6 to 20 eigenvalues. These results indicated the significant conformational changes in the viral protein structure docked with natural compounds to attain the most stable complex formation against reference compounds.
Furthermore, the first three principal components were extracted and plotted to analyze the cluster motion in the protein structure as a function of 100 ns MD simulation. Notably, all of the complexes of JEV–RdRp docked with natural compounds showed higher cluster distributions of −40 and 40 along the direction of PC1, −30 and 30 along the direction of PC2, and −30 and 30 along PC3 by comparison to reference protein docked with ATP and GTP, which showed reduced and overlapped clusters along −30 and 20 and −40 and 20 coordinates for all three PCs (Figure 5 and Figure S6). These observations suggested that the docked natural compounds have the potential to induce conformational fluctuations in JEV–RdRp to distribute its biological function for the replication of the JE virus.

3.4.2. Binding Free Energy Calculation

All of the generated MD simulation trajectories created as a function of 100 ns interval were treated for free binding energy calculation using MM/GBSA method against reference complexes. Figure 6 shows high and low binding energy values for all the docked complexes with respect to time. Notably, Echinacoside, Echinacin, and Rutin were found to exhibit the highest binding energy values (>−100 kcal/mol) for some poses, with net binding free energy values of −80 ± 8.04, 81.67 ± 8.31, and 80.33 ± 5.54 kcal/mol, respectively. In comparison, the other three natural compounds showed mean binding free energy between 66 to 57 kcal/mol. Interestingly, the calculated binding free energy values were relatively higher than those of the reference complexes, i.e., ATP (69.62 ± 8.62 kcal/mol) and GTP (47.98 ± 11.41 kcal/mol) (Figure 6). Additionally, net energy dissociation components—namely, ΔGBind Coulomb, ΔGBind Covalent, ΔGBind Hbond, ΔGBind Lipo, ΔGBind Packing, ΔGBindSolvGB, and ΔGBind vdW—were determined using the Prime MM/GBSA method and were also studied on the complete 100 ns MD trajectories of each complex. Interestingly, substantial contributions of ΔGBind Coulomb, ΔGBind Lipo, and ΔGBind vdW interactions were noted for favorable energy, while ΔGBind Covalent and ΔGBindSolvGB exhibited unfavorable energy to the net binding free energy for all of the protein–ligand complexes during the 100 ns MD simulation (Table 2). Altogether, the natural compounds were noted for higher binding free energy, except Cynaroside, in the selective pocket of JEV–RdRp against reference compounds as a function of 100 ns interval.

4. Conclusions

Japanese encephalitis virus–RNA-dependent RNA polymerase (JEV–RdRp) is a potential target for antiviral drugs, as they are responsible for viral genome replication. The purpose of this study was to identify potent antiviral compounds from E. angustifolia by inhibiting JEV–RdRp using in silico approach. Based on the substantial docking energy (>−10 kcal/Mol) and pharmacokinetics analysis, six compounds—Echinacoside, Echinacin, Rutin, Cynaroside, Quercetagetin 7-glucoside, and Kaempferol-3-glucoside—were selected. Based on the hydrogen, hydrophobic, and other molecular interaction analyses, the stability of the respective compounds with JEV–RdRp’s binding site was studied. The final confirmation of the stability of the protein–ligand complex was made with the analysis of the MD simulation and post-simulation analysis. In conclusion, the results revealed that the selected compounds from E. angustifolia are acceptable inhibitors of the JEV–RdRp protein and can be utilized for further studies for developing a potential antiviral drug against Japanese encephalitis.

Supplementary Materials

The following supporting information can be downloaded at: https://www.mdpi.com/article/10.3390/life12070952/s1.

Author Contributions

Conceptualization: P.Y., V.D.D., E.I.A. and S.K.J.; methodology: P.Y., V.D.D., M.M.E.-D., S.A.E.-K., W.T.A. and A.A.F.; data curation and formal analysis: P.Y., M.M.E.-D., S.A.E.-K., W.T.A. and A.A.F.; investigation: P.Y., S.K.J., V.D.D. and E.I.A.; validation and visualization: P.Y., M.M.E.-D., S.A.E.-K., W.T.A. and A.A.F.; roles/writing—original draft: P.Y. and V.D.D.; writing—review and editing: V.D.D., E.I.A. and S.K.J.; Supervision: V.D.D., E.I.A. and S.K.J. All authors have read and agreed to the published version of the manuscript.

Funding

This research is supported by the charitable donation from the Late Sheikh Ibraheem Ahmed Azhar to cover the article processing charge (APC) as a contribution to the scientific research community.

Informed Consent Statement

Not applicable.

Data Availability Statement

The datasets used and/or analyzed during the current study are available from the corresponding author on reasonable request.

Acknowledgments

The authors acknowledge the generous charitable donation from the late Sheikh Ibraheem Ahmed Azhar as a contribution to the scientific research community. Additionally, the authors are highly thankful to Amaresh Kumar Sahoo, the Indian Institute of Information Technology, Prayagraj, India, for providing his kind support in binding free energy calculation in the Schrodinger 2020.4 suite.

Conflicts of Interest

The authors declare that they have no known competing financial interests or personal relationships that could have appeared to influence the work reported in this paper.

References

  1. Uchil, P.D.; Satchidanandam, V. Phylogenetic analysis of Japanese encephalitis virus: Envelope gene based analysis reveals a fifth genotype, geographic clustering, and multiple introductions of the virus into the Indian subcontinent. Am. J. Trop. Med. Hyg. 2001, 65, 242–251. [Google Scholar] [CrossRef] [PubMed]
  2. Li, M.H.; Fu, S.H.; Chen, W.X.; Wang, H.Y.; Guo, Y.H.; Liu, Q.Y.; Li, Y.X.; Luo, H.M.; Da, W.; Ji, D.Z.D.; et al. Genotype V Japanese Encephalitis Virus Is Emerging. PLoS Negl. Trop. Dis. 2011, 5, e1231. [Google Scholar] [CrossRef] [PubMed]
  3. Solomon, T.; Ni, H.; Beasley, D.W.; Ekkelenkamp, M.; Cardosa, M.J.; Barrett, A.D. Origin and evolution of Japanese encephalitis virus in southeast Asia. J. Virol. 2003, 77, 3091–3098. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  4. Turtle, L.; Solomon, T. Japanese encephalitis—The prospects for new treatments. Nat. Rev. Neurol. 2018, 14, 298–313. [Google Scholar] [CrossRef]
  5. Chen, K.M.; Tsai, H.C.; Sy, C.L.; Lee, S.S.J.; Liu, Y.C.; Wann, S.R.; Wang, Y.H.; Mai, M.H.; Chen, J.K.; Wu, K.S.; et al. Clinical manifestations of Japanese encephalitis in southern Taiwan. J. Microbiol. Immunol. 2009, 42, 296–302. [Google Scholar]
  6. Paul, W.S.; Moore, P.S.; Karabatsos, N.; Flood, S.P.; Yamada, S.; Jackson, T.; Tsai, T.F. Outbreak of Japanese encephalitis on the island of Saipan, 1990. J. Infect. Dis. 1993, 167, 1053–1058. [Google Scholar] [CrossRef]
  7. Wakai, S. Scourge of Japanese encephalitis in south-western Nepal. Lancet 1998, 351, 759. [Google Scholar] [CrossRef]
  8. Unni, S.K.; Ruzek, D.; Chhatbar, C.; Mishra, R.; Johri, M.K.; Singh, S.K. Japanese encephalitis virus: From genome to infectome. Microbes Infect. 2011, 13, 312–321. [Google Scholar] [CrossRef]
  9. Solomon, T.; Dung, N.M.; Kneen, R.; Gainsborough, M.; Vaughn, D.W.; Khanh, V.T. Japanese encephalitis. J. Neurol. Neurosurg. Psychiatry 2000, 68, 405. [Google Scholar] [CrossRef]
  10. Xu, C.; Zhang, W.; Pan, Y.; Wang, G.; Yin, Q.; Fu, S.; Li, F.; He, Y.; Xu, S.; Wang, Z.; et al. A Bibliometric Analysis of Global Research on Japanese Encephalitis From 1934 to 2020. Front. Cell Infect. Microbiol. 2022, 12, 833701. [Google Scholar] [CrossRef]
  11. WHO. Japanese Encephalitis. WHO Fact. Sheets. 2019. Available online: https://www.who.int/en/news-room/fact-sheets/detail/japanese-encephalitis (accessed on 30 September 2021).
  12. Campbell, G.L.; Hills, S.L.; Fischer, M.; Jacobson, J.A.; Hoke, C.H.; Hombach, J.M.; Marfin, A.A.; Solomon, T.; Tsai, T.F.; Tsu, V.D.; et al. Estimated global incidence of Japanese encephalitis: A systematic review. Bull. World Health Organ. 2011, 89, 766–774, 774A–774E. [Google Scholar] [CrossRef] [PubMed]
  13. Yun, S.I.; Lee, Y.M. Japanese encephalitis: The virus and vaccines. Hum. Vaccines Immunother. 2014, 10, 263–279. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  14. Kovanich, D.; Saisawang, C.; Sittipaisankul, P.; Ramphan, S.; Kalpongnukul, N.; Somparn, P.; Pisitkun, T.; Smith, D.R. Analysis of the Zika and Japanese Encephalitis Virus NS5 Interactomes. J. Proteome. Res. 2019, 18, 3203–3218. [Google Scholar] [CrossRef] [PubMed]
  15. Kumar, S.; Verma, A.; Yadav, P.; Dubey, S.K.; Azhar, E.I.; Maitra, S.S.; Dwivedi, V.D. Molecular pathogenesis of Japanese encephalitis and possible therapeutic strategies. Arch. Virol. 2022, 1–24. [Google Scholar] [CrossRef] [PubMed]
  16. Kapoor, M.; Zhang, L.W.; Ramachandra, M.; Kusukawa, J.; Ebner, K.E.; Padmanabhan, R. Association between Ns3 and Ns5 Proteins of Dengue Virus Type-2 in the Putative Rna Replicase Is Linked to Differential Phosphorylation of Ns5. J. Biol. Chem. 1995, 270, 19100–19106. [Google Scholar] [CrossRef] [Green Version]
  17. Uchil, P.D.; Kumar, A.V.; Satchidanandam, V. Nuclear localization of flavivirus RNA synthesis in infected cells. J. Virol. 2006, 80, 5451–5464. [Google Scholar] [CrossRef] [Green Version]
  18. Selisko, B.; Wang, C.; Harris, E.; Canard, B. Regulation of Flavivirus RNA synthesis and replication. Curr. Opin. Virol. 2014, 9, 74–83. [Google Scholar] [CrossRef] [Green Version]
  19. Hodel, A.E.; Gershon, P.D.; Shi, X.N.; Quiocho, F.A. The 1.85 angstrom structure of vaccinia protein VP39: A bifunctional enzyme that participates in the modification of both mRNA ends. Cell 1996, 85, 247–256. [Google Scholar] [CrossRef] [Green Version]
  20. Egloff, M.P.; Benarroch, D.; Selisko, B.; Romette, J.L.; Canard, B. An RNA cap (nucleoside-2′-O-)-methyltransferase in the flavivirus RNA polymerase NS5: Crystal structure and functional characterization. EMBO J. 2002, 21, 2757–2768. [Google Scholar] [CrossRef] [Green Version]
  21. El Sahili, A.; Soh, T.S.; Schiltz, J.; Gharbi-Ayachi, A.; Seh, C.C.; Shi, P.Y.; Lim, S.P.; Lescar, J. NS5 from Dengue Virus Serotype 2 Can Adopt a Conformation Analogous to That of Its Zika Virus and Japanese Encephalitis Virus Homologues. J. Virol. 2019, 94, e01294-19. [Google Scholar] [CrossRef]
  22. Surana, P.; Satchidanandam, V.; Nair, D.T. RNA-dependent RNA polymerase of Japanese encephalitis virus binds the initiator nucleotide GTP to form a mechanistically important pre-initiation state. Nucleic Acids Res. 2014, 42, 2758–2773. [Google Scholar] [CrossRef]
  23. Hoke, C.H.; Nisalak, A.; Sangawhipa, N.; Jatanasen, S.; Laorakapongse, T.; Innis, B.L.; Kotchasenee, S.; Gingrich, J.B.; Latendresse, J.; Fukai, K.; et al. Protection against Japanese encephalitis by inactivated vaccines. N. Engl. J. Med. 1988, 319, 608–614. [Google Scholar] [CrossRef] [PubMed]
  24. Xin, Y.Y.; Ming, Z.G.; Peng, G.Y.; Jian, A.; Min, L.H. Safety of a live-attenuated Japanese encephalitis virus vaccine (SA14-14-2) for children. Am. J. Trop. Med. Hyg. 1988, 39, 214–217. [Google Scholar] [CrossRef] [PubMed]
  25. Fan, Y.C.; Chen, Y.Y.; Chen, J.M.; Huang, C.; Huang, M.; Chiou, S.S. Effectiveness of Live-Attenuated Genotype III Japanese Encephalitis Viral Vaccine against Circulating Genotype I Viruses in Swine. Viruses 2022, 14, 114. [Google Scholar] [CrossRef] [PubMed]
  26. Ngwe Tun, M.M.; Kyaw, A.K.; Nwe, K.M.; Inoue, S.; Thant, K.Z.; Morita, K. Effectiveness of the SA 14-14-2 Live-Attenuated Japanese Encephalitis Vaccine in Myanmar. Vaccines 2021, 9, 568. [Google Scholar] [CrossRef] [PubMed]
  27. Tandale, B.V.; Khan, S.A.; Kushwaha, K.P.; Rahman, H.; Gore, M.M.; Japanese Encephalitis Vaccination Efficacy Case Control Study, G. Effectiveness of Japanese encephalitis SA 14-14-2 live attenuated vaccine among Indian children: Retrospective 1:4 matched case-control study. J. Infect. Public Health 2018, 11, 713–719. [Google Scholar] [CrossRef]
  28. Hennessy, S.; Liu, Z.; Tsai, T.F.; Strom, B.L.; Wan, C.M.; Liu, H.L.; Wu, T.X.; Yu, H.J.; Liu, Q.M.; Karabatsos, N.; et al. Effectiveness of live-attenuated Japanese encephalitis vaccine (SA14-14-2): A case-control study. Lancet 1996, 347, 1583–1586. [Google Scholar] [CrossRef]
  29. Gao, N.; Chen, W.; Zheng, Q.; Fan, D.Y.; Zhang, J.L.; Chen, H.; Gao, G.F.; Zhou, D.S.; An, J. Co-expression of Japanese encephalitis virus prM-E-NS1 antigen with granulocyte-macrophage colony-stimulating factor enhances humoral and anti-virus immunity after DNA vaccination. Immunol. Lett. 2010, 129, 23–31. [Google Scholar] [CrossRef]
  30. Wu, S.F.; Lee, C.J.; Liao, C.L.; Dwek, R.A.; Zitzmann, N.; Lin, Y.L. Antiviral effects of an iminosugar derivative on flavivirus infections. J. Virol. 2002, 76, 3596–3604. [Google Scholar] [CrossRef] [Green Version]
  31. Chang, C.C.; Ou, Y.C.; Raung, S.L.; Chen, C.J. Antiviral effect of dehydroepiandrosterone on Japanese encephalitis virus infection. J. Gen. Virol. 2005, 86, 2513–2523. [Google Scholar] [CrossRef] [Green Version]
  32. Sebastian, L.; Desai, A.; Shampur, M.N.; Perumal, Y.; Sriram, D.; Vasanthapuram, R. N-methylisatin-beta-thiosemicarbazone derivative (SCH 16) is an inhibitor of Japanese encephalitis virus infection in vitro and in vivo. Virol. J. 2008, 5, 64. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  33. Chang, S.J.; Chang, Y.C.; Lu, K.Z.; Tsou, Y.Y.; Lin, C.W. Antiviral Activity of Isatis indigotica Extract and Its Derived Indirubin against Japanese Encephalitis Virus. Evid.-Based Complement. Altern. Med 2012, 2012, 925830. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  34. Wang, S.; Liu, Y.; Guo, J.; Wang, P.; Zhang, L.; Xiao, G.; Wang, W. Screening of FDA-Approved Drugs for Inhibitors of Japanese Encephalitis Virus Infection. J. Virol. 2017, 91, e01055-17. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  35. Nawa, M.; Takasaki, T.; Yamada, K.I.; Kurane, I.; Akatsuka, T. Interference in Japanese encephalitis virus infection of Vero cells by a cationic amphiphilic drug, chlorpromazine. J. Gen. Virol. 2003, 84, 1737–1741. [Google Scholar] [CrossRef]
  36. Ye, J.; Jiang, R.; Cui, M.; Zhu, B.; Sun, L.; Wang, Y.; Zohaib, A.; Dong, Q.; Ruan, X.; Song, Y.; et al. Etanercept reduces neuroinflammation and lethality in mous.se model of Japanese encephalitis. J. Infect. Dis. 2014, 210, 875–889. [Google Scholar] [CrossRef]
  37. Mishra, M.K.; Basu, A. Minocycline neuroprotects, reduces microglial activation, inhibits caspase 3 induction, and viral replication following Japanese encephalitis. J. Neurochem. 2008, 105, 1582–1595. [Google Scholar] [CrossRef]
  38. Dwivedi, V.D.; Singh, A.; El-Kafraway, S.A.; Alandijany, T.A.; Faizo, A.A.; Bajrai, L.H.; Kamal, M.A.; Azhar, E.I. Mechanistic insights into the Japanese encephalitis virus RNA dependent RNA polymerase protein inhibition by bioflavonoids from Azadirachta indica. Sci. Rep. 2021, 11, 18125. [Google Scholar] [CrossRef]
  39. Dwivedi, V.D.; Bharadwaj, S.; Afroz, S.; Khan, N.; Ansari, M.A.; Yadava, U.; Tripathi, R.C.; Tripathi, I.P.; Mishra, S.K.; Kang, S.G. Anti-dengue infectivity evaluation of bioflavonoid from Azadirachta indica by dengue virus serine protease inhibition. J. Biomol. Struct. Dyn. 2021, 39, 1417–1430. [Google Scholar] [CrossRef]
  40. Dwivedi, V.D.; Tripathi, I.P.; Tripathi, R.C.; Singh, R.; Yadava, U.; Mishra, S.K. In silico docking of quercetin-3-O-β-D-glucoside from Azadirachta indica with NS2B-NS3 protease in Dengue virus. Online J. Bioinform. 2018, 19, 175–180. [Google Scholar]
  41. Ali, S.I.; Sheikh, W.M.; Rather, M.A.; Venkatesalu, V.; Muzamil Bashir, S.; Nabi, S.U. Medicinal plants: Treasure for antiviral drug discovery. Phytother. Res. PTR 2021, 35, 3447–3483. [Google Scholar] [CrossRef]
  42. Ganjhu, R.K.; Mudgal, P.P.; Maity, H.; Dowarha, D.; Devadiga, S.; Nag, S.; Arunkumar, G. Herbal plants and plant preparations as remedial approach for viral diseases. Virusdisease 2015, 26, 225–236. [Google Scholar] [CrossRef] [PubMed]
  43. Nedyalkova, M.; Vasighi, M.; Sappati, S.; Kumar, A.; Madurga, S.; Simeonov, V. Inhibition Ability of Natural Compounds on Receptor-Binding Domain of SARS-CoV2: An In Silico Approach. Pharmaceuticals 2021, 14, 1328. [Google Scholar] [CrossRef] [PubMed]
  44. Bhosale, S.; Kumar, A. Screening of phytoconstituents of Andrographis paniculata against various targets of Japanese encephalitis virus: An in-silico and in-vitro target-based approach. Curr. Res. Pharmacol. Drug Discov. 2021, 2, 100043. [Google Scholar] [CrossRef] [PubMed]
  45. Kajimura, K.; Takagi, Y.; Ueba, N.; Yamasaki, K.; Sakagami, Y.; Yokoyama, H.; Yoneda, K. Protective effect of astragali radix by oral administration against Japanese encephalitis virus infection in mice. Biol. Pharm. Bull. 1996, 19, 1166–1169. [Google Scholar] [CrossRef] [Green Version]
  46. Swarup, V.; Ghosh, J.; Mishra, M.K.; Basu, A. Novel strategy for treatment of Japanese encephalitis using arctigenin, a plant lignan. J. Antimicrob. Chemoth. 2008, 61, 679–688. [Google Scholar] [CrossRef] [Green Version]
  47. Johari, J.; Kianmehr, A.; Mustafa, M.R.; Abubakar, S.; Zandi, K. Antiviral activity of baicalein and quercetin against the Japanese encephalitis virus. Int. J. Mol. Sci. 2012, 13, 16785–16795. [Google Scholar] [CrossRef]
  48. Barnes, J.; Anderson, L.A.; Gibbons, S.; Phillipso.on, J.D. Echinacea species (Echinacea angustifolia (DC.) Hell., Echinacea pallida (Nutt.) Nutt., Echinacea purpurea (L.) Moench): A review of their chemistry, pharmacology and clinical properties. J. Pharm. Pharmacol. 2005, 57, 929–954. [Google Scholar] [CrossRef] [Green Version]
  49. Signer, J.; Jonsdottir, H.R.; Albrich, W.C.; Strasser, M.; Zust, R.; Ryter, S.; Ackermann-Gaumann, R.; Lenz, N.; Siegrist, D.; Suter, A.; et al. In vitro virucidal activity of Echinaforce(R), an Echinacea purpurea preparation, against coronaviruses, including common cold coronavirus 229E and SARS-CoV-2. Virol. J. 2020, 17, 136. [Google Scholar] [CrossRef]
  50. Aucoin, M.; Cooley, K.; Saunders, P.R.; Cardozo, V.; Remy, D.; Cramer, H.; Abad, C.N.; Hannan, N. The effect of quercetin on the prevention or treatment of COVID-19 and other respiratory tract infections in humans: A rapid review. Adv. Integr. Med. 2020, 7, 247–251. [Google Scholar] [CrossRef]
  51. Hudson, J.; Vimalanathan, S. Echinacea—A Source of Potent Antivirals for Respiratory Virus Infections. Pharmaceuticals 2011, 4, 1019–1031. [Google Scholar] [CrossRef] [Green Version]
  52. Bharadwaj, S.; El-Kafrawy, S.A.; Alandijany, T.A.; Bajrai, L.H.; Shah, A.A.; Dubey, A.; Sahoo, A.K.; Yadava, U.; Kamal, M.A.; Azhar, E.I.; et al. Structure-Based Identification of Natural Products as SARS-CoV-2 M(pro) Antagonist from Echinacea angustifolia Using Computational Approaches. Viruses 2021, 13, 305. [Google Scholar] [CrossRef] [PubMed]
  53. Hudson, J.B. Applications of the phytomedicine Echinacea purpurea (Purple Coneflower) in infectious diseases. J. Biomed. Biotechnol. 2012, 2012, 769896. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  54. Berman, H.M.; Westbrook, J.; Feng, Z.; Gilliland, G.; Bhat, T.N.; Weissig, H.; Shindyalov, I.N.; Bourne, P.E. The Protein Data Bank. Nucleic Acids Res. 2000, 28, 235–242. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  55. Kim, S.; Chen, J.; Cheng, T.; Gindulyte, A.; He, J.; He, S.; Li, Q.; Shoemaker, B.A.; Thiessen, P.A.; Yu, B.; et al. PubChem in 2021: New data content and improved web interfaces. Nucleic Acids Res. 2021, 49, D1388–D1395. [Google Scholar] [CrossRef] [PubMed]
  56. Pettersen, E.F.; Goddard, T.D.; Huang, C.C.; Couch, G.S.; Greenblatt, D.M.; Meng, E.C.; Ferrin, T.E. UCSF Chimera—A visualization system for exploratory research and analysis. J. Comput. Chem. 2004, 25, 1605–1612. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  57. Dallakyan, S.; Olson, A.J. Small-molecule library screening by docking with PyRx. Methods Mol. Biol. 2015, 1263, 243–250. [Google Scholar] [CrossRef]
  58. Bharadwaj, S.; Lee, K.E.; Dwivedi, V.D.; Yadava, U.; Nees, M.; Kang, S.G. Density functional theory and molecular dynamics simulation support Ganoderma lucidum triterpenoids as broad range antagonist of matrix metalloproteinases. J. Mol. Liq. 2020, 311, 113322. [Google Scholar] [CrossRef]
  59. Trott, O.; Olson, A.J. AutoDock Vina: Improving the speed and accuracy of docking with a new scoring function, efficient optimization, and multithreading. J. Comput. Chem. 2010, 31, 455–461. [Google Scholar] [CrossRef] [Green Version]
  60. Schrödinger Release 2020-4; Maestro, Schrödinger, LLC: New York, NY, USA, 2020.
  61. Daina, A.; Michielin, O.; Zoete, V. SwissADME: A free web tool to evaluate pharmacokinetics, drug-likeness and medicinal chemistry friendliness of small molecules. Sci. Rep. 2017, 7, 42717. [Google Scholar] [CrossRef] [Green Version]
  62. Bowers, K.J.; Chow, D.E.; Xu, H.; Dror, R.O.; Eastwood, M.P.; Gregersen, B.A.; Klepeis, J.L.; Kolossvary, I.; Moraes, M.A.; Sacerdoti, F.D. Scalable algorithms for molecular dynamics simulations on commodity clusters. In Proceedings of the SC’06: Proceedings of the 2006 ACM/IEEE Conference on Supercomputing, Tampa, FL, USA, 11–17 November 2006; p. 43. [Google Scholar]
  63. Schrödinger Release 2018-4: Desmond Molecular Dynamics System; D. E. Shaw Research: New York, NY, USA; Maestro-Desmond Interoperability Tools, Schrödinger: New York, NY, USA, 2018.
  64. Martyna, G.J.; Klein, M.L.; Tuckerman, M. Nose-Hoover Chains—The Canonical Ensemble Via Continuous Dynamics. J. Chem. Phys. 1992, 97, 2635–2643. [Google Scholar] [CrossRef]
  65. Toukmaji, A.Y.; Board, J.A. Ewald summation techniques in perspective: A survey. Comput. Phys. Commun. 1996, 95, 73–92. [Google Scholar] [CrossRef]
  66. Grant, B.J.; Rodrigues, A.P.; ElSawy, K.M.; McCammon, J.A.; Caves, L.S. Bio3d: An R package for the comparative analysis of protein structures. Bioinformatics 2006, 22, 2695–2696. [Google Scholar] [CrossRef] [PubMed] [Green Version]
Life 12 00952 g001 550
Figure 1. Molecular redocked poses of the selected natural compounds: (a) Echinacoside, (b) Echinacin, (c) Rutin, (d) Cynaroside, (e) Quercetagetin 7-glucoside, and (f) Kaempferol-3-glucoside, with comparison to the reference compounds (g) ATP and (h) GTP in the selected binding pocket of JEV–RdRp showing intermolecular interactions with residues extracted at 4 Å around the docked ligand.
Figure 1. Molecular redocked poses of the selected natural compounds: (a) Echinacoside, (b) Echinacin, (c) Rutin, (d) Cynaroside, (e) Quercetagetin 7-glucoside, and (f) Kaempferol-3-glucoside, with comparison to the reference compounds (g) ATP and (h) GTP in the selected binding pocket of JEV–RdRp showing intermolecular interactions with residues extracted at 4 Å around the docked ligand.
Life 12 00952 g001
Life 12 00952 g002 550
Figure 2. Representation of docked natural compounds: (a) Echinacoside, (b) Echinacin, (c) Rutin, (d) Cynaroside, (e) Quercetagetin 7-glucoside, and (f) Kaempferol-3-glucoside, in the selective binding pocket of JEV–RdRp before and after the 100 ns MD simulation.
Figure 2. Representation of docked natural compounds: (a) Echinacoside, (b) Echinacin, (c) Rutin, (d) Cynaroside, (e) Quercetagetin 7-glucoside, and (f) Kaempferol-3-glucoside, in the selective binding pocket of JEV–RdRp before and after the 100 ns MD simulation.
Life 12 00952 g002
Life 12 00952 g003 550
Figure 3. RMSD plots for the docked complexes of JEV–RdRP with selected natural compounds: (a) Echinacoside, (b) Echinacin, (c) Rutin, (d) Cynaroside, (e) Quercetagetin 7-glucoside, and (f) Kaempferol-3-glucoside, with comparison to the reference compounds (g) ATP and (h) GTP.
Figure 3. RMSD plots for the docked complexes of JEV–RdRP with selected natural compounds: (a) Echinacoside, (b) Echinacin, (c) Rutin, (d) Cynaroside, (e) Quercetagetin 7-glucoside, and (f) Kaempferol-3-glucoside, with comparison to the reference compounds (g) ATP and (h) GTP.
Life 12 00952 g003
Life 12 00952 g004 550
Figure 4. Protein–ligand interaction plots for the selected natural compounds: (a) Echinacoside, (b) Echinacin, (c) Rutin, (d) Cynaroside, (e) Quercetagetin 7-glucoside, and (f) Kaempferol-3-glucoside by comparison to the reference compounds (g) ATP and (h) GTP, docked with the JEV–RdRp; the plots were extracted from the 100 ns MD simulation interval.
Figure 4. Protein–ligand interaction plots for the selected natural compounds: (a) Echinacoside, (b) Echinacin, (c) Rutin, (d) Cynaroside, (e) Quercetagetin 7-glucoside, and (f) Kaempferol-3-glucoside by comparison to the reference compounds (g) ATP and (h) GTP, docked with the JEV–RdRp; the plots were extracted from the 100 ns MD simulation interval.
Life 12 00952 g004
Life 12 00952 g005 550
Figure 5. Principal component analysis for the generated molecular dynamics trajectories of JEV–RdRp docked with selected natural compounds: (a) Echinacoside, (b) Echinacin, (c) Rutin, (d) Cynaroside, (e) Quercetagetin 7-glucoside, and (f) Kaempferol-3-glucoside. The percentage of total mean square displacement of residual positional variations recorded in each dimension is categorized by equivalent eigenvalues or PCs.
Figure 5. Principal component analysis for the generated molecular dynamics trajectories of JEV–RdRp docked with selected natural compounds: (a) Echinacoside, (b) Echinacin, (c) Rutin, (d) Cynaroside, (e) Quercetagetin 7-glucoside, and (f) Kaempferol-3-glucoside. The percentage of total mean square displacement of residual positional variations recorded in each dimension is categorized by equivalent eigenvalues or PCs.
Life 12 00952 g005
Life 12 00952 g006 550
Figure 6. Binding free energy calculated on the 5000 poses generated during 100 ns MD simulation of the selected natural compounds: (a) Echinacoside, (b) Echinacin, (c) Rutin, (d) Cynaroside, (e) Quercetagetin 7-glucoside, and (f) Kaempferol-3-glucoside, (g) Adenosine Triphosphate, (h) Guanosine-5′-Triphosphate, docked with JEV–RdRp protein.
Figure 6. Binding free energy calculated on the 5000 poses generated during 100 ns MD simulation of the selected natural compounds: (a) Echinacoside, (b) Echinacin, (c) Rutin, (d) Cynaroside, (e) Quercetagetin 7-glucoside, and (f) Kaempferol-3-glucoside, (g) Adenosine Triphosphate, (h) Guanosine-5′-Triphosphate, docked with JEV–RdRp protein.
Life 12 00952 g006
Table
Table 1. Redocking scores and intermolecular interaction analysis for the potent natural compounds against JEV–RdRp extracted at 4 Å around the docked ligand.
Table 1. Redocking scores and intermolecular interaction analysis for the potent natural compounds against JEV–RdRp extracted at 4 Å around the docked ligand.
S. No.CompoundDocking Score (kcal/mol)H-Bondπ–π/π–Cation StackingHydrophobicPolarNegativePositiveGlycineSalt Bridge
a.Echinacoside−11.1Arg460(2), Glu461, Ser604(2), Ser666, Asp669Arg474Ile476, Val509, Tyr610, Phe713, Cys714, Trp800, Ile802Ser473, Ser604, Thr609, Asn613, Ser666, Ser801, His803Glu461,
Asp668, Asp669		Arg460, Lys462, Lys463, Lys471,
Arg474, Arg742		Gly412,
Gly472, Gly605,
Gly667		-
b.Echinacin−10.4Leu411, Tyr610, Cys714, Ser715, Arg734, Arg742,
Ser801		-Leu411, Ala413, Ile476, Val607, Tyr610, Cys714, Trp800, Ile802Ser604, Gln606, Thr609, Ser715, Asn716, Ser799 Ser801, His803Glu738Arg460, Arg474,
Arg734, Arg742,		Gly412-
c.Rutin−10.4Lyn463, Ala475, Ser604, Gln606, Ser666, Asp668Tyr610Ala413, Val414, Ala475, Ile476, Trp477, Val607, Tyr610, Cys714, Trp800, Ile802Ser604, Gln606, Thr609, Ser666, Ser715, Ser801, His803Asp668,
Asp669		Arg460, Lys471, Arg474, Arg742Gly412,
Gly605,
Gly667		-
d.Cynaroside−10.0Ser666, Asp668, Asp669, His803-Leu411, Ala413, Val607, Tyr610, Ala611, Cys714, Ile802Asn495, Ser604, Gln606, Thr609, Asn613, Ser666, Ser801, His803Glu510,
Asp668,
Asp669		Lys404, Lys471, Hip498Gly412,
Gly605,
Gly667		-
e.Quercetagetin 7-glucoside−10.0Ser666, Asp668,
Ser801		-Leu411, Ala413, Val607, Tyr610, Ala611, Cys714, Ile802Asn495, Ser604, Gln606, Thr609, Asn613, Ser666, Ser801, His803Glu510,
Asp668,
Asp669		Lys404, Lys471, Arg474Gly412,
Gly605,
Gly667		-
f.Kaempferol-3-glucoside−10.0Asp541, Ser604, Ser801-Leu411, Ala413, Trp540, Val607, Tyr610, Trp800, Ile802Asn495, Ser604, Gln605, Thr609, Asn613, Ser666, Ser801Asp541,
Asp668		Arg474Gly412,
Gly667,		-
GAdenosine triphosphate−8.6Ser604, Gln606, Ser715, Arg742, Trp800, Ile802--Leu411, Ala413, Val607, Tyr610, Trp800, Ile802Ser604, Gln606, Thr609, Ser715, Ser799, Ser801-Arg460, Lys463, Lys471, Arg474,
Arg734, Arg742		Gly412, Gly605Lys463, Arg474,
Arg734, Arg742
hGuanosine-5′-triphosphate−9.0Asp541, Asn613, Asp668, Ser801, Ile802Tyr610Trp540, Tyr610, Trp800, Ile802Ser604, Asn613, Thr609, Ser666, Ser801, His803Asp541, Asp668, Asp669Arg460, Lys463, Arg474, Arg742Gly667Arg460, Lys463,
Arg474, Arg742
Table
Table 2. Calculated net binding free energy for the selected docked poses of JEV–RdRp natural compounds snapshots from the last 10 ns interval of 100 ns MD simulation.
Table 2. Calculated net binding free energy for the selected docked poses of JEV–RdRp natural compounds snapshots from the last 10 ns interval of 100 ns MD simulation.
Energy (kcal/mol)JEV–RdRP Poses with Natural Compounds from Echinacea angustifolia
EchinacosideEchinacinRutinCynarosideQuercetagetin 7-GlucosideKaempferol-3-GlucosideAdenosine TriphosphateGuanosine-5′-Triphosphate
ΔGBind−80.15 ± 8.04−81.67 ± 8.31−80.33 ± 5.54−57.26 ± 4.13−66.24 ± 5.6−62.13 ± 5.97−69.62 ± 8.62−47.98 ± 11.41
ΔGBindCoulomb−57.04 ± 8.1−32.93 ± 8.82−41.87 ± 5.69−28.94 ± 6.36−30.43 ± 6.17−28.68 ± 5.86−242.97 ± 54.21−242.62 ± 52.45
ΔGBindCovalent4.62 ± 2.384.17 ± 2.115.19 ± 2.516.06 ± 2.14.29 ± 1.873.94 ± 1.677.83 ± 4.455.13 ± 2.45
ΔGBindHbond−7.11 ± 1−3.95 ± 0.55−6.11 ± 0.69−3.68 ± 0.79−4.56 ± 0.92−3.46 ± 0.63−15.77 ± 2.08−10.65 ± 2.24
ΔGBindLipo−18.26 ± 2.1−16.1 ± 1.22−17.43 ± 1.37−13.59 ± 1.49−16.4 ± 1.04−16.22 ± 1.34−5.09 ± 0.58−3.56 ± 0.61
ΔGBindPacking−2.04 ± 0.63−7.16 ± 1.27−0.71 ± 0.17−3.35 ± 0.62−2.07 ± 0.35−2.38 ± 0.76−0.15 ± 0.28−1.48 ± 0.61
ΔGBindSolvGB66.69 ± 6.1141.52 ± 3.646.21 ± 3.4442.08 ± 3.1842.84 ± 4.7239.37 ± 3.82224.6 ± 49.37236.56 ± 49.6
ΔGBindvdW−67 ± 4.05−67.22 ± 3.38−65.54 ± 3.71−55.81 ± 3.28−59.9 ± 3.17−54.7 ± 3.41−38.07 ± 4.69−31.35 ± 4.63
Ligand Strain Energy11.92 ± 3.737.91 ± 2.3710.64 ± 3.466.0 ±2.135.2 ± 2.195.57 ± 2.517.55 ± 4.063.99 ± 2.95
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Yadav, P.; El-Kafrawy, S.A.; El-Day, M.M.; Alghafari, W.T.; Faizo, A.A.; Jha, S.K.; Dwivedi, V.D.; Azhar, E.I. Discovery of Small Molecules from Echinacea angustifolia Targeting RNA-Dependent RNA Polymerase of Japanese Encephalitis Virus. Life 2022, 12, 952. https://doi.org/10.3390/life12070952

AMA Style

Yadav P, El-Kafrawy SA, El-Day MM, Alghafari WT, Faizo AA, Jha SK, Dwivedi VD, Azhar EI. Discovery of Small Molecules from Echinacea angustifolia Targeting RNA-Dependent RNA Polymerase of Japanese Encephalitis Virus. Life. 2022; 12(7):952. https://doi.org/10.3390/life12070952

Chicago/Turabian Style

Yadav, Pardeep, Sherif A. El-Kafrawy, Mai M. El-Day, Wejdan T. Alghafari, Arwa A. Faizo, Saurabh Kumar Jha, Vivek Dhar Dwivedi, and Esam I. Azhar. 2022. "Discovery of Small Molecules from Echinacea angustifolia Targeting RNA-Dependent RNA Polymerase of Japanese Encephalitis Virus" Life 12, no. 7: 952. https://doi.org/10.3390/life12070952

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop