Molecular Docking, Dynamics Simulations, ADMET, and DFT Calculations: Combined <i>In Silico</i> Approach to Screen Natural Inhibitors of 3CL and PL Proteases of SARS-CoV-2

Mohanasundaram, Sugumar; Karthikeyan, Porkodi; Sampath, Venkatesan; Anbazhagan, M.; Prabhu, Sundramurthy Venkatesa; Khaled, Jamal M.; Thiruvengadam, Muthu

doi:https://doi.org/10.1155/2024/6647757

Cellular Microbiology

On this page

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

Research Article | Open Access

Volume 2024 | Article ID 6647757 | https://doi.org/10.1155/2024/6647757

Molecular Docking, Dynamics Simulations, ADMET, and DFT Calculations: Combined In Silico Approach to Screen Natural Inhibitors of 3CL and PL Proteases of SARS-CoV-2

Sugumar Mohanasundaram,¹Porkodi Karthikeyan,²Venkatesan Sampath,²M. Anbazhagan,³Sundramurthy Venkatesa Prabhu,^4,5Jamal M. Khaled,⁶and Muthu Thiruvengadam⁷

Academic Editor: Shiek Ahmed

Received19 Jun 2023

Revised15 Nov 2023

Accepted21 Dec 2023

Published16 Jan 2024

Abstract

Considering natural compounds for the antiviral effect is another opportunity for exploring novel drug candidates for severe acute respiratory syndrome coronavirus 2. The selected natural compounds were interacted using a molecular docking approach. The 3D structures of the main protease and papain-like protease were used for the virtual screening to detect the potent inhibitor against SARS-CoV-2. The top-scored compounds were further analyzed for absorption, digestion, metabolism, excretion, and toxicity properties and density functional theory analysis. Our results indicated that glycyrrhizin exhibited better docking scores of -9.5 kcal/mol with main protease and -9.7 kcal/mol with papain-like protease. Next to glycyrrhizin, rutin showed a better docking score of -9.1 kcal/mol and -9.2 kcal/mol with 3-chymotrypsin-like and papain-like proteases. Violaxanthin and naringin occupied the subsequent position in the docking score table with 3CL and PL proteases, respectively. In addition, the crucial properties like drug likeliness and pharmacokinetics of the compounds were determined. There is no significant toxicity identified.

1. Introduction

SARS-CoV-2, the first pandemic of this century, has been detected in 216 countries, with more than 756 million confirmed cases and 68,44,267 deaths as of February 20, 2023. The alteration happening in the genome of this virus resulting in different variants is the biggest challenge among the scientific community to predict its characteristic features for drug development. The UK found the first mutant variety (called alpha) in September 2020. In May 2020, South Africa found beta, and in November 2020, Brazil found gamma. By October 2020, India had reported the delta variant, and by December 2020, Peru had registered the lambda variant. In January of 2021, Colombia found the Mu species.

The omicron variant has been reported in more than 50 countries (coronavirus disease (COVID-19) pandemic, https://www.who.int/emergencies/diseases/novel-coronavirus-2019) [1]. Practical diagnostic approaches have been developed, including RT-PCR, CT imaging, and interleukin-based serological screening [2].

Since the abrupt prerequisite for an effective therapeutic option against SARS-CoV-2 is in the relentless phase of the explosion, evaluating the prevailing antiviral agents for repurposing is obligatory. Many countries conduct in vivo and human trials with existing drugs for possible repurposing. Hence, virtual screening of antiviral medications against SARS-CoV-2 will institute a platform for the lead identification process.

Over 250 genome sequences have been deposited into GenBank (GenBank, taxid: 2697049). This β-coronavirus shares an 89% resemblance to the bat coronavirus [3, 4]. A 29 kb size, single-stranded RNA encoding 9,860 amino acids, is present. Since December 2019, researchers have identified several proteins as targets. Main protease, RNA polymerase, and spike protein have attracted outstanding attention owing to their significant role in interactions with the host (human) and the subsequent replication.

Viral proteases are endopeptidases that catalyze the cleavage of viral polyproteins [5]. Coronaviridae genomes carry two polyproteins, which get hewed and changed into nonstructural proteins (NSP) by the actions of the 3CL protease and the PL protease [6]. The NSP has a role in viral replication in host cells [7]. As a result, targeting those proteases to mimic the transformation of NSP into biologically functional proteins is worthwhile. The SARS papain-like protease cleaves ISG15 and polyUb [8, 9], whereas 3CL protease is indispensable to the viral life cycle [10].

Food and Drug Administration (FDA) of the USA has not approved any antiviral medicine to treat SARS-CoV-2 except remdesivir under emergency use authorization (EUA). The focus of research has shifted to evaluating the available antiviral drugs. An efficient, naturally occurring drug must be identified and capable of targeting ACE-2 receptors, spike protein, RdRp, or SARS-CoV-2 proteases (https://www.guidetopharmacology.org/coronavirus.jsp).

Literature studies indicate that saquinavir, remdesivir, and darunavir and two natural drugs, flavone and coumarin derivatives, were employed to inhibit 3CL protease [11]. Herbal treatment was demonstrated as effective in suppressing infectious diseases during the 2003 SARS outbreak [12]. Significant similarities can be found between SARS-CoV-2 and MERS-CoV [13, 14]. Thus, screening for antiviral compounds in herbals previously used to treat viral pneumonia will expedite the identification of a potential SARS-CoV-2 antiviral drug. In India, Andrographis paniculata was widely used to control the dengue outbreak in 2006, following the Ministry of Ayush’s recommendations. Traditionally used by the Tamil people of South India, “Nilavembu Kudineer” is an aqueous preparation of Andrographis paniculata. It was administered to dengue-infected patients for treatment and to healthy individuals as a prophylactic immune booster. Its antidengue properties have been extensively explored and demonstrated [15]. In India, papaya leaves were used to treat dengue [16]. An investigation was conducted with 49 naturally occurring metabolites derived from 23 distinct medicinal plant species. Among the total of 49 compounds that were subjected to screening, only 7 exhibited drug-like properties. The study employed molecular docking, dynamics simulation, ADME/T analysis, and DFT analysis to identify naringenin’s most promising ligand [17]. Another similar study confirmed that artificially synthesized acyl phloroglucinols showed notable inhibitory potential against SARS-CoV-2 PL protease [17].

Other plants considered include Azadirachta indica (Neem), Eclipta prostrata, Carica papaya, and Zingiber officinale, all extensively examined for their antiviral activity [18]. The preceding data concluded that natural antiviral chemicals are critical in mimicking infectious diseases such as SARS, dengue, and chikungunya. As a result, antiviral compounds are chosen for this investigation from the sources listed above.

The objectives of this study were the effective interaction of selected natural compounds with the chosen protein targets by in silico docking analysis followed by dynamics simulation analysis and conceptual density functional theory (DFT) computation. Finally, the absorption, digestion, metabolism, excretion, and toxicity (ADMET) properties were investigated to determine their drug-like properties.

2. Materials and Methods

2.1. Retrieval of Protein Receptor

The 3D structures of 3CL protease (PDB ID: 7C8T) and PL protease (PDB ID: 6WX4) were downloaded from the Protein Data Bank based on the resolution of crystal structure (Figure 1).

(a)

(b)

2.2. Retrieval of Ligands

3D structures of antiviral compounds (Table 1) such as 6-gingerol, catechin, chlorogenic acid, coumaric acid, caricaxanthin, dasyscyphin C, glycyrrhizin, hyperoside, kaempferol, naringin, nimbaflavone, quercetin, rutin, violaxanthin, zeaxanthin, and zingiberene were downloaded from PubChem.

2.3. Ligand and Protein Preparation

The ligand molecules are prepared by energy minimization using ArgusLab software [19]. The energy minimization has been done by the restricted Hartree-Fock (closed-shell). Until attaining the lowest potential energy, all the ligand molecules were energy minimized. The energy-minimized ligands were used for docking studies.

The 3D structures of chosen protein molecules were directly downloaded from PDB (PDB ID: 7C8T and PDB ID: 6WX4) and used for docking. The water molecules and ligand (inhibitor) were removed from the structure, and the essential polar hydrogen was added to the structure using Discovery Studio Visualizer (version 20.1.0.19295) developed by Dassault Systems BIOVIA Corporation. The prepared protein structure was further finalized using an online server of the Centre for Molecular and Biomolecular Informatics, Radboud University, Netherlands (https://swift.cmbi.umcn.nl/servers/html/index.html), and then saved for docking in the form of PDB format [20].

2.4. Molecular Docking

PyRx (version 0.8) was used to dock the protein receptors with the selected ligands, and the final output files were obtained as PDBQT files. Individual poses are prepared by a vina split [21]. The amino acids involved in different interactions were observed.

2.5. Conceptual DFT

DFT helps analyze a structure by molecular orbital energies [22]. It may help clarify the structure-activity relationship of compounds. DFT calculation will help link biological activities to structure [23]. The Hohenberg and Kohn theorems underpin it [24]. This work used a subset of DFT called conceptual DFT to investigate chemistry using electron density ideas [25]. The electron density of a molecule yields ten separate molecular characteristics, including the electrophilicity index, chemical potential, molecular dipole moment, and electronegativity.

2.6. Prediction of ADME Properties

All the ligands used in this study were not approved as drugs by the FDA. As a result, it is critical to understand the drug likelihood of the best ligands targeting both 3CL and PL proteases. The ligands’ ADME properties were foreseen by SwissADME (https://www.swissadme.ch). The ADME properties were observed [26] for the ligands that showed the best docking score with 3C and papain-like proteases of SARS-CoV-2.

2.7. Prediction of Toxicity

The selected best ligands with significant docking scores were subjected to toxicity prediction using the pkCSM (http://biosig.unimelb.edu.au/pkcsm/). The canonical SMILES of the target ligands, taken from the PubChem website, were utilized on the pkCSM server under toxicity mode for prediction [27].

2.8. Molecular Dynamics Simulation

The glycyrrhizin, rutin, and violaxanthin ligands were simulated using the GROMACS 5.1.2 software and the GROMOS96 54a7 force field. The study mentioned above employed traditional algorithmic techniques and methodologies. The PDBQT files of the docked complexes’ main protease and drug molecules were converted to PDB files using Discovery Studio 2020. These PDB data were subsequently refined using Swiss PDB Viewer. GROMACS generated the core protease topology file, while PRODRG2 was responsible for creating the topologies of the drug candidates. Topology and data for protein-ligand complexes were generated. Four sodium ions were introduced into a dodecahedron container with a volume of 1 nanometer to preserve electrostatic charge equilibrium, which was subsequently filled with water. The energy was reduced after 1,000 iterations using the steepest descent algorithm. The drug candidates and main protease were later constrained, allowing the solvent to diffuse without restriction. The conditions required for achieving the constant number of particles, volume, and temperature (NVT) and continuous number of particles, pressure, and temperature (NPT) equilibration were as follows. The output analogue coordinates consisted of a total of 500 steps. The computational technique employed in this study was the Particle Mesh Ewald approach, with a cut-off distance of 1.2 nm for the electrostatic interactions. Following the completion of NVT (canonical ensemble) and NPT (isothermal-isobaric ensemble) simulations, production molecular dynamics (MD) runs were conducted at a temperature of 300 K and a pressure of 1 bar, with a time step of 2 femtoseconds. A 20,000-picosecond simulation was performed [20, 28].

3. Results

3.1. Analysis of Docking with 3CL Protease

To understand the molecular interaction, all compounds were docked with the 3CL protease (PDB ID: 7C8T). Our results indicated that glycyrrhizin had a better docking value of -9.5 kcal/mol and coumaric acid had a lower docking value of -5.3 kcal/mol (Figure 2). Table 2 lists the amino acids of the 3CL protease that form chemical bonds with each ligand. The docking image (Figure 3) of glycyrrhizin with 3C protease revealed seven hydrogen bonds with CYS145, ARG131, LYS137, THR199, LEU287, ASP289, and GLU288 and one electrostatic interaction with LYS137. Compared to glycyrrhizin, rutin (Figure 4) and violaxanthin (Figure 5) demonstrated significant docking scores with 3C protease. Rutin formed six hydrogen bonds and three additional contacts, while violaxanthin formed three hydrophobic connections (Table 2).

3.2. Docking Analysis with PL Protease

The PL protease (PDB ID: 6WX4) was docked individually with each identified ligand. Glycyrrhizin was discovered with the highest docking value of -9.7 kcal/mol (Figure 6). Zingiberene had the lowest docking value of -4.8 kcal/mol. The amino acid details of each ligand’s bonding with the protein receptor are listed in Table 3. According to the Discovery Studio Visualizer (Figure 7), glycyrrhizin forms three hydrogen bonds with PL protease via LEU211, TYR305, and GLY209. Compared to glycyrrhizin, rutin (Figure 8) and naringin (Figure 9) had significant docking scores. Rutin formed four hydrogen bonds and seven other interactions with the PL protease, whereas naringin formed four hydrogen bonds and four other interactions.

3.3. Prediction of ADME Properties

The docking investigation revealed that glycyrrhizin and rutin have a better binding affinity for 3CL and PL proteases, respectively. It is worth noting that glycyrrhizin has a high potential for interaction with the 3CL and PL proteases. Compared to glycyrrhizin, rutin has a high potential for interaction with both proteases. ADME properties of the top four ligands are determined by their docking score for both proteases. As an outcome, glycyrrhizin, rutin, violaxanthin, and naringin were chosen for ADME property prediction using the SwissADME server.

The SwissADME results (Table 4) indicate that none of the ligands can inhibit the cytochrome P450 system, which is a significant finding. All ligands have a negligible level of intestinal absorption rather than none. None of the ligands are capable of passing through the blood-brain barrier. Rutin and naringin are found to be water-soluble, whereas glycyrrhizin and violaxanthin are found to be insoluble and very poorly soluble in water. Compared to the other ligands, glycyrrhizin and violaxanthin do not contain heavy aromatic atoms. To be regarded as a drug-like molecule in ADMET prediction, any molecule must adhere to Lipinski’s rule of 5. Here, the final chosen molecules violate some of Lipinski’s rules, which must be evaluated before evaluating any ligand as a prospective drug-likeness compound based on its ADME features.

3.4. Conceptual DFT

The phytochemicals were optimized using the B3LYP function [29, 30] with a 6-31G (2p, d) basis in Gaussian 16 [31] to determine their molecular properties [32]. The highest occupied molecular orbital (HOMO) energy () and the least unoccupied molecular orbital (LUMO) energy () were calculated, which are critical descriptors that refer to a molecule’s ability to donate as well as accept electrons, correspondingly. We developed and examined maps reflecting the electron density in various locations of the molecules at HOMO and LUMO. The and values for the phytochemicals are listed in Table 5.

Table 6 shows the electron density maps of selected phytochemicals’ molecular orbitals. We built and evaluated density maps of molecular orbitals. The energy gap () was determined like . The energy gap is proportionate to the molecules’ reactivity [33]. As a result, with an value of -1.22 eV, the rutin protease had the smallest energy gap. Glycyrrhizin, on the other hand, showed the largest, measuring 4.97 eV.

A molecule’s molecular dipole moment can be related to its chemical reactivity [34]. Rutin protease and glycyrrhizin exhibited dipole moments greater than 7.0 Debye. Using predicted and scores, electronegativity and other properties of the ligands were determined. The electron density of a compound is significant because it influences the electronegativity. If the electronegativity is low, the inhibitory activity will be higher and effective [35]. Glycyrrhizin has the lowest electronegativity value of -3.76 eV, while the rutin protease has a value of 4.44 eV. Correlations are found between the values of descriptors for phytochemicals derived using conceptual DFT and the docking analysis. From this, glycyrrhizin was the highest scoring molecule with huge electronegativity than rutin, which is consistent with the docking analysis.

3.5. Prediction of Toxicity

The toxicity of the final selected ligands, glycyrrhizin, rutin, naringin, and violaxanthin, was predicted and analyzed [36]. The results (Table 7) demonstrate unequivocally that none of the four ligands contain AMES toxicity and is neither mutagenic nor carcinogenic. The maximum recommended tolerated dose in humans for all ligands is less than 0.477, which is considered very low. None of the ligands were found to be hepatotoxic or skin-sensitizing. Glycyrrhizin and violaxanthin showed no inhibitory potential on hERG genes encoding for potassium channels. The LD₅₀ and minnow toxicity values are significantly high for all the ligands, which shows that only at very high doses, ligands may tend to cause toxicity (http://biosig.unimelb.edu.au/pkcsm/static/help/pkcsm_theory.pdf).

3.6. Molecular Dynamics Simulation

Following the molecular docking and ADMET prediction analysis, the best three phytochemicals were subjected to research the stability of ligands with 3CL and PL proteases using a molecular dynamics simulation study. Root-mean-square deviation (RMSD), root-mean-square fluctuation (RMSF), and protein-ligand contacts were calculated. Glycyrrhizin, rutin, and violaxanthin molecules show good binding affinity and dug-likeness properties in this docking and ADMET analysis. RMSD, RMSF, Rg, and solvent-accessible surface are shown in Figures 10–12. The average RMSD value of glycyrrhizin, rutin, and violaxanthin molecules with 3CL and PL proteases shows 0.19 nm, 0.18 nm, and 0.18 nm and 0.2 nm, 0.21 nm, and 0.2 nm, respectively. The stability of the 3CL with glycyrrhizin, rutin, and violaxanthin in radius of gyration (total and around axes) shows at 2.35 Rg (nm) up to 20000 ps.

The results obtained in this investigation are consistent and comparable with the findings reported in the earlier study conducted by Jang et al. [37]. The results demonstrate the stability of the relationship between the natural metabolites and the protein structure. The results obtained for the root-mean-square deviation (RMSD) and root-mean-square fluctuation (RMSF) in our study demonstrate comparable or superior performance compared to the findings reported in earlier studies. The molecular dynamics simulation conducted for 2000 picoseconds showed encouraging interactions between the protein and the selected metabolites. In the investigation mentioned above, the compounds that were chosen had protein-based root-mean-square deviation (RMSD) values of 0.5, 0.4, and 0.3 nm, respectively (refer to Figure 3(d)). Additionally, their root-mean-square deviation (RMSF) values displayed comparable patterns. Moreover, the in vitro assessment demonstrated that this interaction has a higher degree of stability, indicating inhibitory potential. The findings closely align with our results.

4. Discussions

Virtual screening of the existing antiviral drugs and the naturally available compounds with antiviral properties is essential for the treatment. The widely used drug, remdesivir, is still in trial and is authorized to be used only for emergency management. A recent report [38, 39] stated that combining azithromycin and hydroxychloroquine significantly aided in emergency cases. The inhibition of proteases was observed due to the increase in pH. However, the adverse effects of hydroxychloroquine and azithromycin must be considered, and a side-effect-free alternative is urgently needed [39].

Indian spices are of great importance because of their traditional use as a source of medicine. Many bioactive metabolites with therapeutic applications have been identified [40]. Plants are well characterized and proven for their noteworthy basis of high-value bioactive metabolites, including alkaloids, phenols, flavonoids, coumarins, steroids, and lignans. Plants have widened the window of opportunity to discover natural compounds with therapeutic potential [41–43]. According to a study [44], the plant-derived chemicals carnosol and rosmanol may be employed as therapeutic candidates against SARS-CoV-2. Another study [45] also asserts that three natural substances, ursolic acid, carvacrol, and oleanolic acid, are possible inhibitors.

In this study, four compounds, glycyrrhizin, rutin, violaxanthin, and naringin, were finally chosen from a pool of sixteen compounds based on their docking score.

Glycyrrhizin, the compound with the highest docking score, has been reported to be a moderately potent antiviral agent against HIV and herpes simplex virus [46]. It was recently found as a possible inhibitor of HIV’s replicative mechanism [47]. Glycyrrhizin inhibits many viruses, including herpes, corona-alpha, flaviviruses, HIV, Epstein–Barr virus, influenza A virus, vaccinia, and polio type I viruses [48–53]. Our findings correlated with earlier research, demonstrating a potent 3CL and PL protease inhibitors. The drug is also effective against SARS in vitro studies and may show promise for treating COVID-19 infections. More promising, some derivatives of glycyrrhizin are shown to have a much higher antiviral activity than glycyrrhizin itself [54]. The in silico results of this paper are similar to a recent study [55] that has clearly explained the suitability of glycyrrhizin against the 3CL protease of SARS-CoV-2 through in vitro assays, and EC₅₀ value was found to be 0.44 mg/mL.

Additionally, conceptual DFT analyses demonstrate that all molecular characteristics indicate that this compound is a potent inhibitor. The combination of glycyrrhizin and boswellic acid demonstrated efficacy in reducing death rates, accelerating the healing process, and enhancing prognosis or reducing clinical status scores on a 7-point scale. The laboratory findings indicate notable disparities in the serum levels of C-reactive protein (CRP) and the percentage of lymphocytes between the placebo group and the intervention group, providing evidence to support the efficacy of GR+BA in promoting improvement. The combination above has qualities of safety and affordability, rendering it a potential candidate for addressing mild to moderate infections caused by SARS-CoV-2 or its variations associated with COVID-19 [56]. Glycyrrhizin had superior inhibitory effectiveness against SARS-CoV-2 ( mg/L) during the processes of viral adsorption and penetration, in comparison to other antiviral agents such as ribavirin, 6-aziridine, pyrazofurin, and mycophenolic acid, while also demonstrating a higher selectivity index [57].

Rutin has the second-highest docking score. It is a glycosidic derivative of a flavonoid available in natural resources. According to a study [58], the compound’s 50% effective concentration (EC₅₀) against the cucumber mosaic virus was 394.78 g/mL. A prior work [59] demonstrated that sodium rutin sulfate could be a sulfated rutin analogue that can inhibit HIV-1. Additionally, rutin has been shown [45] to significantly inhibit the avian influenza strain H5N1 in the Madin-Darby canine kidney.

Additionally, there is evidence [60] that rutin significantly inhibits the viral replication cycle during canine distemper virus infection during adsorption and diffusion. The current results of the molecular docking study demonstrate the compound’s antiviral activity and suggest that it may be a significant inhibitor of both 3CL and PL proteases. It has been studied [61] that rutin can inhibit 3CL protease. The findings of the conceptual DFT calculation also indicate that it can be a probable inhibitor next to glycyrrhizin.

Violaxanthin is a naturally occurring orange-coloured xanthophyll pigment found in plants. In silico analysis demonstrated a moderate antidengue potential for this compound [62]. There is scant evidence for violaxanthin’s antiviral properties; however, there is evidence for its anti-inflammatory, antiproliferative, and antioxidant properties. Thus, more research and in vitro testing are required to demonstrate its antiviral efficacy.

The fourth ligand with an exceptionally high docking score, naringin, a flavanone 7-O-glycoside, is found in citrus fruits, particularly grapefruits. Naringin was found to have antiadsorption properties against dengue virus—type 2 in a study [63], with an IC₅₀ of 168.2 g/mL and a corresponding selectivity index of 1.3. Another study [64] reported that naringin interacted successfully with the chikungunya virus’s nsp2, inhibiting viral replication at the initiation stage. Naringin also suppresses human rotavirus (Wa-1 strain) [65, 66]. The docking research results indicate that it can interact strongly with 3CL and PL proteases. Using conceptual DFT calculations, violaxanthin and naringin are placed third and fourth in potency.

Molecular dynamics simulation validated the conversational stability of glycyrrhizin, rutin, and violaxanthin molecules with the main protease. The RMSD of these ligands showed minor fluctuation with each protein conformation at 20 ns. Also, the binding nature of glycyrrhizin, rutin, and violaxanthin ligands with 3CL did not produce any conformational stability of the protein. Furthermore, the RMSF analysis shows that the best ligands were firmly bound to the active site of the 3CL protein and did not change the loop sections. The gyration radius was also examined to illustrate the protein’s structural alterations. Figures 10–12 show the time development of the Rg of 3CL subjected to applied electric fields of various strengths and directions. The value of Rg rises with the power of the applied electric field V/nm. When 1BBL is exposed to electric fields, the importance of Rg becomes increasingly exceptional V/nm as the intensity of the electric field increases due to the protein’s extension along the electric field. The protein-ligand association was stable during MD modelling, with slight backbone changes.

According to the toxicity prediction, glycyrrhizin and violaxanthin exhibited no toxicity with a higher LD₅₀ value. Additionally, rutin and naringin inhibited the hERG II gene. The chosen ligands are prospective antiviral medicines due to their effective interaction with proteases. As per the results of ADME property prediction, glycyrrhizin is the best choice for an antiviral molecule with specific breaches of Lipinski’s rule. Additionally, no substantial toxicity was detected using toxicological prediction. Interactions with SARS-CoV-2 proteases in silico demonstrate their ability to work against various viruses. While glycyrrhizin, rutin, and naringin all exhibit promising antiviral effects, their drug-likeness violations should not be overlooked when considering them as drug lead compounds. As a result, glycyrrhizin is the best antiviral drug and may be a probable alternative treatment option.

5. Conclusion

Repurposing existing medications and utilizing natural antiviral compounds can be the best scientific option until the development of a new drug. The primary target of most viral infections, including HIV, dengue, SARS, influenza, and n-COVID-19, is a protease. The in silico screening of sixteen natural antiviral medicines for their capacity to inhibit the SARS-CoV-2 protease suggests that glycyrrhizin and rutin may be deemed more powerful antiviral agents than others. Additionally, their ADME and toxicology features indicate their potential as anti-SARS-CoV-2 agents. Also, molecular dynamics simulation analysis reveals the strong interaction of top ligands with the protein targets. As a result, additional research and trials are necessary to establish them as antiviral agents.

Data Availability

The data used to support the findings of this study are included in the article.

All authors agree to the publishing of the paper.

Conflicts of Interest

The authors declare no conflict of interest.

Authors’ Contributions

Sugumar Mohanasundaram was responsible for the methodology and investigation and wrote the original draft. Porkodi Karthikeyan was responsible for the software and supervision and wrote, reviewed, and edited the manuscript. Venkatesan Sampath was responsible for the conceptualization and supervision. Anbazhagan. M was responsible for the methodology and wrote, reviewed, edited, and validated the manuscript. Sundramurthy Venkatesa Prabhu was responsible for the conceptualization and supervision and wrote, reviewed, and edited the manuscript. Jamal M. Khaled was responsible for the funding acquisition and validation. Muthu Thiruvengadam was responsible for the conceptualization and supervision.

Acknowledgments

The authors extend their appreciation to the Deputyship for Research and Innovation, “Ministry of Education” in Saudi Arabia for funding this research (IFKSUOR3-457-2).

References

“Coronavirus disease (COVID-19) pandemic,” 2019, http://www.who.int/emergencies/diseases/novel-coronavirus-2019.
View at: Google Scholar
A. Sharma, F. I. Ahmad, and S. K. Lal, “COVID-19: a review on the novel coronavirus disease evolution, transmission, detection, control and prevention,” Viruses, vol. 13, no. 2, p. 202, 2021.
View at: Publisher Site | Google Scholar
J. F. W. Chan, K. H. Kok, Z. Zhu et al., “Genomic characterization of the 2019 novel human-pathogenic coronavirus isolated from a patient with atypical pneumonia after visiting Wuhan,” Emerging Microbes & Infections, vol. 9, no. 1, pp. 221–236, 2020.
View at: Publisher Site | Google Scholar
R. Lu, X. Zhao, J. Li et al., “Genomic characterisation and epidemiology of 2019 novel coronavirus: implications for virus origins and receptor binding,” Lancet, vol. 395, no. 10224, pp. 565–574, 2020.
View at: Publisher Site | Google Scholar
J. Konvalinka, H. G. Kräusslich, and B. Müller, “Retroviral proteases and their roles in virion maturation,” Virology, vol. 479-480, pp. 403–417, 2015.
View at: Publisher Site | Google Scholar
R. Hilgenfeld, “From SARS to MERS: crystallographic studies on coronaviral proteases enable antiviral drug design,” The FEBS Journal, vol. 281, no. 18, pp. 4085–4096, 2014.
View at: Publisher Site | Google Scholar
V. Thiel, K. A. Ivanov, A. Putics et al., “Mechanisms and enzymes involved in SARS coronavirus genome expression,” The Journal of General Virology, vol. 84, no. 9, pp. 2305–2315, 2003.
View at: Publisher Site | Google Scholar
V. G. Bhoj and Z. J. Chen, “Ubiquitylation in innate and adaptive immunity,” Nature, vol. 458, no. 7237, pp. 430–437, 2009.
View at: Publisher Site | Google Scholar
M. K. Isaacson and H. L. Ploegh, “Ubiquitination, ubiquitin-like modifiers, and deubiquitination in viral infection,” Cell Host & Microbe, vol. 5, no. 6, pp. 559–570, 2009.
View at: Publisher Site | Google Scholar
P. Mukherjee, F. Shah, P. Desai, and M. Avery, “Inhibitors of SARS-3CLpro: virtual screening, biological evaluation, and molecular dynamics simulation studies,” Journal of Chemical Information and Modeling, vol. 51, no. 6, pp. 1376–1392, 2011.
View at: Publisher Site | Google Scholar
S. A. Khan, K. Zia, S. Ashraf, R. Uddin, and Z. Ul-Haq, “Identification of chymotrypsin-like protease inhibitors of SARS-Cov-2 via integrated computational approach,” Journal of Biomolecular Structure and Dynamics, vol. 39, no. 7, pp. 2607–2616, 2020.
View at: Publisher Site | Google Scholar
Z. Chen and T. Nakamura, “Statistical evidence for the usefulness of Chinese medicine in the treatment of SARS,” Phytotherapy Research, vol. 18, no. 7, pp. 592–594, 2004.
View at: Publisher Site | Google Scholar
H. Y. Tian, “2019-nCoV: new challenges from coronavirus,” Chinese with Abstract in English, vol. 54, p. E001, 2020.
View at: Publisher Site | Google Scholar
P. Zhou, X. L. Yang, X. G. Wang et al., “A pneumonia outbreak associated with a new coronavirus of probable bat origin,” Nature, vol. 579, no. 7798, pp. 270–273, 2020.
View at: Publisher Site | Google Scholar
K. Krishnasamy, S. Ramalingam, S. Karupannan et al., “Anti-dengue activity of Andrographis paniculata extracts and quantification of dengue viral inhibition by SYBR green reverse transcription polymerase chain reaction,” An International Quarterly Journal of Research in Ayurveda, vol. 39, no. 2, pp. 87–91, 2018.
View at: Publisher Site | Google Scholar
P. Senthilvel, P. Lavanya, K. M. Kumar et al., “Flavonoid from Carica papaya inhibits NS2B-NS3 protease and prevents dengue 2 viral assembly,” Bioinformation, vol. 9, no. 18, pp. 889–895, 2013.
View at: Publisher Site | Google Scholar
W. E. Prasetyo, H. Purnomo, M. Sadrini, F. R. Wibowo, M. Firdaus, and T. Kusumaningsih, “Identification of potential bioactive natural compounds from Indonesian medicinal plants against 3-chymotrypsin-like protease (3CLpro) of SARS-CoV-2: molecular docking, ADME/T, molecular dynamic simulations, and DFT analysis,” Journal of Biomolecular Structure and Dynamics, vol. 41, no. 10, pp. 4467–4484, 2023.
View at: Publisher Site | Google Scholar
B. N. Dhawan, “Antiviral activity of Indian plants,” Proceedings of the National Academy of Sciences, India Section B, vol. 82, no. 1, pp. 209–224, 2012.
View at: Publisher Site | Google Scholar
M. A. Thompson, “Molecular docking using ArgusLab, an efficient shape-based search algorithm and AScore scoring function,” in Proceedings of the ACS Meeting, Philadelphia, Pa, USA, 2004, 172: CINF 42.
View at: Google Scholar
S. Sivakumar, S. Mohanasundaram, N. Rangarajan, V. Sampath, and M. V. Dass Prakash, “In silico prediction of interactions and molecular dynamics simulation analysis of Mpro of severe acute respiratory syndrome caused by novel coronavirus 2 with the FDA-approved nonprotein antiviral drugs,” Journal of Applied Pharmaceutical Science, vol. 12, no. 5, pp. 104–119, 2022.
View at: Publisher Site | Google Scholar
R. T. Kroemer, “Structure-based drug design: docking and scoring,” Current Protein & Peptide Science, vol. 8, no. 4, pp. 312–328, 2007.
View at: Publisher Site | Google Scholar
P. Geerlings, E. Chamorro, P. K. Chattaraj et al., “Conceptual density functional theory: status, prospects, issues,” Theoretical Chemistry Accounts, vol. 139, no. 2, p. 36, 2020.
View at: Publisher Site | Google Scholar
S. K. Nagarajan, S. Babu, H. Sohn, P. Devaraju, and T. Madhavan, “Toward a better understanding of the interaction between somatostatin receptor 2 and its ligands: a structural characterization study using molecular dynamics and conceptual density functional theory,” Journal of Biomolecular Structure and Dynamics, vol. 37, no. 12, pp. 3081–3102, 2019.
View at: Publisher Site | Google Scholar
P. Hohenberg and W. Kohn, “Inhomogeneous electron gas,” Physics Review, vol. 136, no. 3B, pp. B864–B871, 1964.
View at: Publisher Site | Google Scholar
L. R. Domingo, M. Ríos-Gutiérrez, and P. Pérez, “Applications of the conceptual density functional theory indices to organic chemistry reactivity,” Molecules, vol. 21, no. 6, p. 748, 2016.
View at: Publisher Site | Google Scholar
A. Daina, O. Michielin, and V. Zoete, “SwissADME: a free web tool to evaluate pharmacokinetics, drug-likeness and medicinal chemistry friendliness of small molecules,” Scientific Reports, vol. 7, no. 1, article 42717, 2017.
View at: Publisher Site | Google Scholar
D. E. Pires, T. L. Blundell, and D. B. Ascher, “PkCSM: predicting small-molecule pharmacokinetic and toxicity properties using graph-based signatures,” Journal of Medicinal Chemistry, vol. 58, no. 9, pp. 4066–4072, 2015.
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,” Software X, vol. 1-2, pp. 19–25, 2015.
View at: Publisher Site | Google Scholar
A. Becke, “Density-functional exchange-energy approximation with correct asymptotic behavior,” Physical Review A, vol. 38, no. 6, pp. 3098–3100, 1988.
View at: Publisher Site | Google Scholar
C. Lee, W. Yang, and R. Parr, “Development of the Colle-Salvetti correlation-energy formula into a functional of the electron density,” Physical Review B, vol. 37, no. 2, pp. 785–789, 1988.
View at: Publisher Site | Google Scholar
M. Frisch, G. Trucks, H. Schlegel et al., Gaussian 16 (version revision B.01) [Linux], Gaussian, Inc., Wallingford CT, 2016.
K. Fukui, “The role of frontier orbitals in chemical reactions (Nobel lecture),” Angewandte Chemie International Edition in English, vol. 21, no. 11, pp. 801–809, 1982.
View at: Publisher Site | Google Scholar
R. Bostan, S. Varvara, L. Găină, and L. M. Mureşan, “Evaluation of some phenothiazine derivatives as corrosion inhibitors for bronze in weakly acidic solution,” Corrosion Science, vol. 63, pp. 275–286, 2012.
View at: Publisher Site | Google Scholar
B. Mert, M. Erman Mert, G. Kardas, and B. Yazıcı, “Experimental and theoretical investigation of 3-amino-1,2,4-triazole-5-thiol as a corrosion inhibitor for carbon steel in HCl medium,” Corrosion Science, vol. 53, no. 12, pp. 4265–4272, 2011.
View at: Publisher Site | Google Scholar
C. Zhan, J. Nichols, and D. Dixon, “Ionization potential, electron affinity, electronegativity, hardness, and electron excitation energy: molecular properties from density functional theory orbital energies,” The Journal of Physical Chemistry, vol. 107, no. 20, pp. 4184–4195, 2003.
View at: Publisher Site | Google Scholar
M. Wang, R. Cao, L. Zhang et al., “Remdesivir and chloroquine effectively inhibit the recently emerged novel coronavirus (2019-nCoV) in vitro,” Cell Research, vol. 30, no. 3, pp. 269–271, 2020.
View at: Publisher Site | Google Scholar
J. H. Kim, Y.-I. Park, M. Hur et al., “Inhibition by components of Glycyrrhiza uralensis of 3CLpro and HCoV-OC43 proliferation,” Journal of Enzyme Inhibition and Medicinal Chemistry, vol. 38, no. 1, p. 1, 2023.
View at: Publisher Site | Google Scholar
Y. K. Bokam, C. Guntupalli, S. N. K. R. Gudhanti et al., “Importance of pharmacists as a front line warrior in improving medication compliance in Covid-19 patients,” Indian Journal of Pharmaceutical Sciences, vol. 83, no. 2, pp. 398–401, 2021.
View at: Publisher Site | Google Scholar
S. K. Enmozhi, K. Raja, I. Sebastine, and J. Joseph, “Andrographolide as a potential inhibitor of SARS-CoV-2 main protease: an in silico approach,” Journal of Biomolecular Structure and Dynamics, vol. 39, no. 9, pp. 1–7, 2020.
View at: Publisher Site | Google Scholar
K. R. A. Sachan, S. Kumar, K. Kumari, and D. Singh, “Medicinal uses of our spices used in our traditional culture,” Journal of Medicinal Plant Studies, vol. 6, no. 3, pp. 116–122, 2018.
View at: Google Scholar
B. Jee, S. Kumar, R. Yadav, Y. Singh, A. Kumar, and M. Sharma, “Ursolic acid and carvacrol may be potential inhibitors of dormancy protein small heat shock protein 16.3 of Mycobacterium tuberculosis,” Journal of Biomolecular Structure and Dynamics, vol. 36, no. 13, pp. 3434–3443, 2018.
View at: Publisher Site | Google Scholar
A. Kumar, R. Kumar, M. Sharma, U. Kumar, M. N. P. Gajula, and K. P. Singh, “Uttarakhand medicinal plants database (UMPDB): a platform for exploring genomic, chemical, and traditional knowledge,” Data (Data), vol. 3, no. 1, p. 7, 2018.
View at: Publisher Site | Google Scholar
S. S. Panchangam, M. Vahedi, M. J. Megha et al., “Saffron‘omics’: the challenges of integrating omic technologies,” Avicenna Journal of Phytomedicine, vol. 6, no. 6, pp. 604–620, 2016.
View at: Google Scholar
D. K. Umesh, C. Selvaraj, S. K. Singh, and V. K. Dubey, “Identification of new anti-nCoV drug chemical compounds from Indian spices exploiting SARS-CoV-2 main protease as target,” Journal of Biomolecular Structure and Dynamics, vol. 39, no. 9, pp. 1–9, 2020.
View at: Publisher Site | Google Scholar
A. Kumar, G. Choudhir, S. K. Shukla et al., “Identification of phytochemical inhibitors against main protease of COVID-19 using molecular modeling approaches,” Journal of Biomolecular Structure and Dynamics, vol. 39, no. 10, pp. 3760–3770, 2021.
View at: Publisher Site | Google Scholar
G. Williams, C. Leclerq, and A. G. A. C. Knaap, WHO Food Additives Series 54: Safety and Evaluation of Certain Food Additives, WHO, Geneva, 2006.
M. Ito, A. Sato, K. Hirabayashi et al., “Mechanism of inhibitory effect of glycyrrhizin on replication of human immunodeficiency virus (HIV),” Antiviral Research, vol. 10, no. 6, pp. 289–298, 1988.
View at: Publisher Site | Google Scholar
J. Cinatl, B. Morgenstern, G. Bauer, P. Chandra, H. Rabenau, and H. W. Doerr, “Glycyrrhizin, an active component of liquorice roots, and replication of SARS-associated coronavirus,” Lancet, vol. 361, no. 9374, pp. 2045-2046, 2003.
View at: Publisher Site | Google Scholar
J. M. Crance, N. Scaramozzino, A. Jouan, and D. Garin, “Interferon, ribavirin, 6-azauridine and glycyrrhizin: antiviral compounds active against pathogenic flaviviruses,” Antiviral Research, vol. 58, no. 1, pp. 73–79, 2003.
View at: Publisher Site | Google Scholar
S. Harada, “The broad antiviral agent glycyrrhizin directly modulates the fluidity of plasma membrane and HIV-1 envelope,” The Biochemical Journal, vol. 392, no. 1, pp. 191–199, 2005.
View at: Publisher Site | Google Scholar
T. Hattori, S. Ikematsu, A. Koito et al., “Preliminary evidence for inhibitory effect of glycyrrhizin on HIV replication in patients with AIDS,” Antiviral Research, vol. 11, no. 5-6, pp. 255–261, 1989.
View at: Publisher Site | Google Scholar
G. Lampis, D. Deidda, M. Pinza, and R. Pompei, “Enhancement of anti-herpetic activity of glycyrrhizic acid by physiological proteins,” Antiviral Chemistry and Chemotherapy, vol. 12, no. 2, pp. 125–131, 2001.
View at: Publisher Site | Google Scholar
J. C. Lin, “Mechanism of action of glycyrrhizic acid in inhibition of Epstein-Barr virus replication in vitro,” Antiviral Research, vol. 59, no. 1, pp. 41–47, 2003.
View at: Publisher Site | Google Scholar
H. Al-Kamel and O. Grundmann, “Glycyrrhizin as a potential treatment for the novel coronavirus (COVID-19),” Mini-Reviews in Medicinal Chemistry, vol. 21, no. 16, pp. 2204–2208, 2021.
View at: Publisher Site | Google Scholar
L. van de Sand, M. Bormann, M. Alt et al., “Glycyrrhizin effectively inhibits SARS-CoV-2 replication by inhibiting the viral main protease,” Viruses, vol. 13, no. 4, p. 609, 2021.
View at: Publisher Site | Google Scholar
A. A. Gomaa, Y. A. Abdel-Wadood, and M. A. Gomaa, “Glycyrrhizin and boswellic acids, the golden nutraceuticals: multitargeting for treatment of mild-moderate COVID-19 and prevention of post-COVID cognitive impairment,” Inflammopharmacology, vol. 30, no. 6, pp. 1977–1992, 2022.
View at: Publisher Site | Google Scholar
S. Banerjee, S. K. Baidya, N. Adhikari, B. Ghosh, and T. Jha, “Glycyrrhizin as a promising kryptonite against SARS-CoV-2: clinical, experimental, and theoretical evidences,” Journal of Molecular Structure, vol. 1275, article 134642, 2023.
View at: Publisher Site | Google Scholar
H. Yu, Y. Ding, D. Xie et al., “Design, synthesis, and antiviral activity of novel rutin derivatives containing 1, 4-pentadien-3-one moiety,” European Journal of Medicinal Chemistry, vol. 92, pp. 732–737, 2015.
View at: Publisher Site | Google Scholar
J. Tao, Q. Hu, J. Yang et al., “In vitro anti-HIV and -HSV activity and safety of sodium rutin sulfate as a microbicide candidate,” Antiviral Research, vol. 75, no. 3, pp. 227–233, 2007.
View at: Publisher Site | Google Scholar
A. K. Ibrahim, A. I. Youssef, A. S. Arafa, and S. A. Ahmed, “Anti-H5N1 virus flavonoids from Capparis sinaica Veill,” Natural Product Research, vol. 27, no. 22, pp. 2149–2153, 2013.
View at: Publisher Site | Google Scholar
B. Rizzuti, F. Grande, F. Conforti et al., “Rutin is a low micromolar inhibitor of SARS-CoV-2 main protease 3CLpro: implications for drug design of quercetin analogs,” Biomedicine, vol. 9, no. 4, p. 375, 2021.
View at: Publisher Site | Google Scholar
O. V. Carvalho, C. V. Botelho, C. G. Ferreira et al., “In vitro inhibition of canine distemper virus by flavonoids and phenolic acids: implications of structural differences for antiviral design,” Research in Veterinary Science, vol. 95, no. 2, pp. 717–724, 2013.
View at: Publisher Site | Google Scholar
N. Radhakrishnan, K. W. Lam, and M. E. Norhaizan, “Molecular docking analysis of Carica papaya Linn constituents as antiviral agent,” International Food Research Journal, vol. 24, no. 4, pp. 1819–1825, 2017.
View at: Google Scholar
K. Zandi, B. Teoh, 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, no. 1, p. 560, 2011.
View at: Publisher Site | Google Scholar
P. K. Tripathi, J. Singh, N. Gaurav, D. K. Garg, and A. K. Patel, “In-silico and biophysical investigation of biomolecular interaction between naringin and nsP2 of the chikungunya virus,” International Journal of Biological Macromolecules, vol. 160, pp. 1061–1066, 2020.
View at: Publisher Site | Google Scholar
L. A. Savi, T. Caon, A. P. de Oliveira et al., “Evaluation of antirotavirus activity of flavonoids,” Fitoterapia, vol. 81, no. 8, pp. 1142–1146, 2010.
View at: Publisher Site | Google Scholar

Copyright

Copyright © 2024 Sugumar Mohanasundaram 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

753

Downloads

257

Citations