Keywords
SARS-CoV-2, Main protease, ADMET analysis, Molecular docking, Molecular dynamics simulation, Essential oil.
-
Views
-
Downloads
How to cite this article
Hasib RA, Ali MC, Rahman MS et al. A computational biology approach for the identification of potential SARS-CoV-2 main protease inhibitors from natural essential oil compounds. [version 1; peer review: 1 approved]. F1000Research 2021, 10:1313 (https://doi.org/10.12688/f1000research.73999.1)
NOTE: it is important to ensure the information in square brackets after the title is included in all citations of this article.
Export Citation
Sciwheel
EndNote
Ref. Manager
Bibtex
ProCite
Sente
▬
✚
Research Article
A computational biology approach for the identification of potential SARS-CoV-2 main protease inhibitors from natural essential oil compounds.
[version 1; peer review: 1 approved]
PUBLISHED 23 Dec 2021
OPEN PEER REVIEW
REVIEWER STATUS
This article is included in the Bioinformatics gateway.
This article is included in the Emerging Diseases and Outbreaks gateway.
This article is included in the Coronavirus collection.
Abstract
Corresponding author:
Mohammad Abu Hena Mostofa Jamal
Competing interests:
No competing interests were disclosed.
Grant information:
This work was supported by the Ministry of Science and Technology, Government of the People’s Republic of Bangladesh under grant number [ES-361, 2020-21].
The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript.
Copyright:
© 2021 Hasib RA et al.
This is an open access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.
How to cite: Hasib RA, Ali MC, Rahman MS et al. A computational biology approach for the identification of potential SARS-CoV-2 main protease inhibitors from natural essential oil compounds. [version 1; peer review: 1 approved]. F1000Research 2021, 10:1313 (https://doi.org/10.12688/f1000research.73999.1)
First published: 23 Dec 2021, 10:1313 (https://doi.org/10.12688/f1000research.73999.1)
Latest published: 15 Jun 2022, 10:1313 (https://doi.org/10.12688/f1000research.73999.2)
The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript.
There is a newer version
of this article available.
of this article available.
Introduction
A new strain of coronavirus termed severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) is responsible for a variety of respiratory diseases that emerged in Wuhan, Hubei province, China towards the end of 2019, and since then has spread globally.1,2 SARS-CoV-2 is a positive sense, enveloped RNA virus that belongs to the β genus of the coronaviridae family.3 Live animals like bats sold at the local Huanan seafood wholesale market are thought to be a possible point of origin of SARS-CoV-2 as most of the fatalities are inhabited nearby it.4 Patients with SARS-CoV-2 infection suffer from a fever at the preliminary stage of the disease as a clinical symptom with some other symptoms including headache, dry cough, difficulty breathing, and pneumonia.5 Progressive respiratory failure due to alveolar damage with the initiation of the disease can even cause the death of the patient.6 Considering the potential health risks associated with this disease, a therapeutic strategy to reduce the transmission rate is urgently needed, but no specific drugs have yet been discovered.7 Several medications like ritonavir, lopinavir, oseltamivir, and ganciclovir have been assayed as yet for the inhibition of SARS-CoV-2. All such drugs are repurposed, and though they may have shown some effectiveness against SARS-CoV-2, but their adverse effects and poor success rate against emerging diseases cannot be neglected.8 Simultaneously, lopinavir-ritonavir based treatment in hospitalized adult patients with severe SARS-CoV-2 has not shown any significant effect.9 Combined doses of hydroxychloroquine and azithromycin exhibited more satisfactory results in the removal of viral load of SARS-CoV-2 than hydroxychloroquine alone, as claimed by Gautret et al., in their study10 but was refuted by another study mentioning the clinical failure of these two drugs in severe cases.11
A protease enzyme called main protease (Mpro) or chymotrypsin-like protease (3CLPro) and a papain-like protease act in the maneuvering of two polyproteins (1a and 1ab) into non-structural proteins (nsps), which are liable for viral replication.12,13 The SARS-CoV-2 main protease specifically cleaves polypeptide sequences after a residue of glutamine as like Leu-Gln↓ (Ser, Ala, Gly) (↓ indicates the cleavage site), which represents the main protease as a suitable drug target since no human host-cell proteases are identified with such substrate specificity.14–17 SARS-CoV-2 main protease maintains an amino acid sequence analogy of 96% in contrast to SARS-CoV-1 and varies merely at 12 out of 303 places in the amino acid sequence.18 Moreover, SARS-CoV-2 main protease comprises three major domains, namely chymotrypsin-like domain I (residues 10–99), picornavirus 3C protease-like domain II (residues 100–184); these two catalytic domains contain a six-stranded antiparallel β-barrel and the α-helical domain III (residues 201–303); necessary for dimerization. Domain II and domain III are linked together through an extended loop region between residues 185-200.19 A cleft between domains I and II holds the substrate-binding site, and the Cys145 and His41 catalytic dyad residues are also included here.20,21
As mentioned earlier, there is no appropriate medication for SARS-CoV-2; however, several previous studies indicated that phytochemicals and essential oils secreted from medicinal plants may have been used to prevent viral replication.22 Essential oils have been reported to exhibit antiviral effects along with their antibacterial and antifungal activity.23 Compounds derived from essential oils have been proven efficient against the herpes simplex virus (HSV-I and HSV-II) and SARS-CoV by numerous experiments.24,25 An in-silico investigation against SARS-CoV-2 reported that several antiviral agents or other alleviation of SARS-CoV-2 symptoms could be amplified by essential oils.26 So, through this research, we hope to provide evidence that our selected essential oil compounds could be effectively used as a treatment against the deadly SARS-CoV-2 infection. The goal of this study was to investigate the effects of essential oils and phytochemicals obtained from different plants against SARS-CoV-2 main protease using using absorption, distribution, metabolism, elimination, and toxicity (ADMET) profiling, molecular docking, and molecular dynamics simulation.
Methods
Essential oil library generation
In our study, we considered essential oil compounds as a possible anti-SARS-CoV-2 agent. In August 2020, we conducted a randomized study to identify relevant essential oil compounds for the current study. From previous research and observations, we formulated a repository of the bioactive compounds of essential oils. Through checking relevant studies in PubMed, PubChem, Google Scholar, the Web of Science, and Scopus databases, we built a catalog of biologically active essential oil molecules. As inclusion criteria, the compounds were evaluated for their antimicrobial and antiviral properties. Our research included the compounds that had previously shown these two characteristics. Compounds were excluded from the study if they did not have these two key features. The 3D conformation in spatial data file (SDF) format, InChI key, and canonical simplified molecular input line entry systems (SMILEs) of all the compounds were collected from the PubChem database once identified from the search databases.27,28 This research identified a total of 114 essential oil compounds, which were then investigated to see whether they had antiviral properties against SARS-CoV-2 (See extended data file for plant information of all the compounds including their binding affinity scores).29
ADMET analysis
Based on their ADMET properties, all the essential oil compounds were filtered. Qikprop determined the absorption, distribution, metabolism, elimination (ADME) properties of the compounds.30 Physicochemical properties, pharmacokinetics, lipophilicity, water-solubility, medicinal chemistry are some of the ADME properties estimated to analyze drug profiles of top drug candidates.31 In this study, the rule of five, rule of three, solvent accessible surface area (SASA), aqueous solubility (QPlogS), human oral absorption, Caco-2 cell permeability, and other pharmacological and druglikeness properties like predicted octanol/water partition coefficient (QPlogPo/w), predicted water/gas partition coefficient (QPlogPw), and conformation-independent predicted aqueous solubility (CIQPlogS) were analyzed using Qikprop and the outputs were compared to the threshold values. The compounds were filtered in a stepwise process. The ligands were initially preprocessed for ADME profiling using the LigPrep module of the Schrödinger suite 2020–3 (https://www.schrodinger.com/products/ligprep). The processed ligands were then directly inserted as an input file (.maegz) into Qikprop (https://www.schrodinger.com/products/qikprop). For the pharmacokinetic properties mentioned above, an output file (.csv) was generated automatically.
First, compounds that pursued the rule of five and the rule of three were included. Following that, these compounds that passed the rule of five and three tests were assessed based on their maximum oral absorption rate (cut-off value is equal to three). In addition, the toxicity of the compounds was analyzed via the online platform admetSAR (http://lmmd.ecust.edu.cn/admetsar2).32 admetSAR is a free tool that offers a user-friendly interface to search for ADMET properties. It is estimated that there are more than 210,000 ADMET specified datasets for more than 96,000 specific compounds in admetSAR, with 45 different kinds of ADMET-associated properties, which were collated from a wide range of previous studies.33 The toxicity assessment results were generated using the Canonical SMILES of the compounds obtained from PubChem as input data into admetSAR. Here, only negative outcomes in the hERG inhibition, AMES test, and carcinogenic test were taken into account. In summary, the following parameters were used to filter the compounds from the initial 114 identified: Rule of Five = 0 violations, Rule of Three = 0 violations, Human Oral Absorption = 3, hERG inhibition = non-inhibitor, AMES Toxicity = non-AMES Toxic, and Carcinogenic Property = non-carcinogen.
Preparation of receptor
The crystallized structure of the SARS-CoV-2 main protease (PDB ID 6LU7) was compiled from the reserves of the Research Collaboratory for Structural Bioinformatics Protein Data Bank (RCSB PDB) (https://www.rcsb.org).34–36 The RCSB PDB is the world's first open-access digital resource for experimentally determined 3D structures of biological macromolecules and their complexes.37 Using the search, visualisation, and analysis tools of RCSB PDB, users can easily explore chemical interactions that stabilize macromolecules and play a significant role in their interactions and functions.38 During receptor preparation, the water molecules and ligands had been eliminated from it, and the steepest descent and conjugate gradient techniques were utilized for energy minimization. The final receptor was configured using GROMACS 96 43B1 algorithm in Swiss-PdbViewer (version 4.0.4) (https://spdbv.vital-it.ch)39 and Chimera (Amber Force field, version 1.14) (https://www.cgl.ucsf.edu/chimera).40 Swiss-PdbViewer is a user-friendly platform that combines the functions for protein structure modeling, assessment, and manipulation. It has been improved with numerous features that allow users to identify and search for general structural motifs in a set of structures.41 In addition, UCSF Chimera is an interactive tool for visualizing molecular structures and is completely free for non-commercial purposes.42
Ligand preparation and molecular docking
The ligand preparation was carried out by utilizing the open-source virtual screening tool PyRx (https://pyrx.sourceforge.io). To minimize the energy of the ligands, the mmff94 (Merck molecular force field) force field,43 which is available in PyRx, was employed.44 For ligand optimization, the steepest descent technique was used, with a total number of 200 steps. After that, the ligands were transformed into a protein data bank (PDB) format by implementing the open-source chemical toolbox Open Babel (http://openbabel.org/wiki/Main_Page). For this, the SDF format of the ligands retrieved from the PubChem database was used as inputs in Open Babel. Then the autodock protein data bank (PDB), partial charge (Q), and atom type (T) [PDBQT] format of the ligands was created via PyRx for further study and analysis.45 The grid box generation and molecular docking study were performed by utilizing the autodock wizard of PyRx, version 0.8 as the free version (https://pyrx.sourceforge.io). PyRx provides an autodock wizard with a simple user interface, making it a useful tool for computer-aided drug design (CADD).46 During the molecular docking study, ligands were considered flexible, and the protein was assumed to be rigid.44,47 The grid parameter configuration file was initiated utilizing the auto grid engine in autodock (grid box size X, Y, Z; 31.93, 21.49, and 30.66, respectively). As potential inhibitors, ligands with the most negative docking scores were taken. Low-affinity ligand binding entails less intermolecular force between the ligand and its receptor than high-affinity ligand binding does. The lower the binding affinity of the ligand is, the better the binding of the ligand and the receptor will be. Finally, using Biovia discovery studio visualizer (v 4.5) (https://discover.3ds.com/discovery-studio-visualizer), molecular interactions between ligands and receptor were visualized.48
Molecular dynamics simulation
To attain detailed information about the predicted molecular interactions, the best-docked complexes were analyzed for 50 ns of MD simulations utilizing the GROMOS96 43al force field of the ‘GROMACS’ macromolecular simulation system (https://simlab.uams.edu/ProteinWithLigand/protein_with_ligand.html).49 In the beginning, topology parameters for the selected ligands were obtained from the ‘PRODRG’ webserver (http://davapc1.bioch.dundee.ac.uk/cgi-bin/prodrg). PRODRG delivers high-throughput protein-ligand crystallography with quick, automated, and consistent access to small-molecule topologies and parameters. PRODRG produced topologies are often of higher quality than topologies generated by other methods and are acceptable for MD simulation in the GROMACS environment.50
In this study, Mpro-ligand complexes were solvated in a cubic box (0.15 M) with the SPC (spc216) water model. Cubic boxes are generally utilized in molecular dynamics simulations for the sake of geometric simplicity.51 Despite considering polarization, the SPC water model is one the most reliable three-center waters models. It provides accurate water density and dielectric permittivity in molecular dynamics simulations.52
For each simulation method, the steepest gradient technique (5000 cycles) was implemented to minimize the energy. Under periodic boundary conditions with a predefined pressure of 1 bar, the equilibration period was accomplished for 100 ps at a constant temperature of 300 K. Finally, the simulation outputs like root mean square deviation (RMSD), root mean square fluctuations (RMSF), radius of gyration (Rg), solvent accessible surface area (SASA), and the number of hydrogen bonds were calculated employing the built-in features of "GROMACS". The lower RMSD, Rg, and SASA value represents the compactness of the system where the high RMSF value indicates the structural flexibility. Likewise, intermolecular hydrogen bonds formed during protein-ligand MD simulation are crucial for the conformational stability of the complex.
Results
ADMET profiling
The propensity of a drug to acquire pharmacologically active concentration at targeted therapies can be evaluated by a set of ADME parameters. Evaluation of ADME may be efficient in combination with toxicity appraisement in silico models.53 The ADMET properties and drug-likeness of each essential oil compound were determined in this study (See underlying data file, ADME.xlsx).54 We did not consider the ability of the compounds to cross the blood-brain barrier (BBB) because SARS-Cov-2 is not directly linked with our central nervous system. Out of 114 compounds tested, 41 were found to violate the rule of five, rule of three, or both. Just one of the remaining 73 compounds had an oral absorption rate of less than three. As a result, the remaining 72 compounds were screened based on the hERG inhibition, AMES, and carcinogenic tests (See underlying data file, Toxicity.xlsx).54 Finally, 53 essential oil compounds met the druglikeness and ADMET profiles requirements outlined in the methods. Following that, a molecular docking approach was then used to examine these 53 filtered compounds.
Molecular docking
Since there is no specific medication to date for the treatment of SARS-CoV-2, it is very urgent to identify feasible biological compounds which block SARS-CoV-2 main protease and act as potential anti-SARS-CoV-2 agents. In this investigation, we docked 53 essential oil molecules against SARS-CoV-2 main protease (PDB ID: 6lu7) that passed the ADMET criteria and evaluated the binding affinity of all these probable inhibitor compounds to the receptor-binding site (See underlying data file for all binding affinity scores).54 We chose the top three potential compounds based on their higher binding affinity scores. The PubChem Id, molecular formula, molecular weight, IUPAC name, 2D structures, and binding affinity of the selected hits including the positive control α-ketoamide are shown in Table 1.
Table 1. PubChem Id, molecular formula, molecular weight, IUPAC name, 2D structure, and binding affinity of selected compounds and the control.
Among the compounds, bisabololoxide B (obtained from Matricaria recutita) disclosed the highest binding affinity of −6.6 kcal/mol, and both eremanthin (obtained from Laurus nobilis) and leptospermone (obtained from Leptospermum scoparium) exhibited a binding affinity of −6.3 kcal/mol during docking analysis. Zhang et al.,53 have shown in a recent analysis that the positive inhibitor α-ketoamide interacts with the main protease residues His41, Gly143, Ser144, Cys145, His163, His164, Glu166, Pro168, and Gln189.55 The Cys-His catalytic dyad (Cys-145 and His-41) is the target region for the inhibitors to bind in. Our selected three hit compounds interact with either Cys-145 or His-41 or both by forming different types of bonds.
Bisabololoxide B is demonstrated in Figure 1(b), it interacted with A: HIS-41 and A: MET-49 by forming alkyl bonds. Eremanthin, demonstrated in Figure 1(c), formed conventional hydrogen bonds with A: SER-144 and A: GLY-143. At the same time, eremanthin interacted with A: CYS-145 by forming both conventional hydrogen and alkyl bonds, whereas it formed only alkyl bonds with A: HIS-163 and A: HIS-172. Leptospermone, demonstrated in Figure 1(d), interacted with A: CYS-145, A: SER-144, and A: GLY-143 through strong hydrogen bonds. It also interacted with A: HIS-41 and A: MET-49 by forming alkyl bonds and with A: MET-165 by carbon-hydrogen bond. Such findings indicated that our selected essential oil compounds could be effective in treating SARS-CoV-2.
Figure 1. Molecular interactions between selected essential oil compounds and SARS-Cov-2 main protease (Mpro).
(a) Molecular interactions of positive control α-ketoamide with Mpro, (b) interactions between Bisabololoxide B and Mpro, (c) between Eremanthin and Mpro, and (d) between Leptospermone and Mpro.
Figure 2. Molecular dynamics simulation analysis of (a) root mean square deviation (RMSD); (b) root mean square fluctuations (RMSF); (c) radius of gyration (Rg); (d) solvent accessible surface area (SASA); and (e) number of hydrogen bonds.
Molecular dynamics simulation
The molecular interactions can describe biological systems in a dynamic fashion, but usually, molecular docking defines only a single whiff of protein-ligand interaction.56 Hence, we simulated the dynamic nature of all the top three complexes to capture their different conformations attained in the solvated condition. Five parameters, RMSD, RMSF, Rg, SASA, and the number of hydrogen bonds were utilized in our study to analyze the dynamic behavior of the top three protein-ligand complexes. The average RMSD values for α-ketoamide-Mpro, bisabololoxide B-Mpro, eremanthin-Mpro and leptospermone-Mpro complex were 0.34 Å, 0.33 Å, 0.32 Å, and 0.23 Å, respectively where leptospermone showed the lowest RMSD value (shown in Table 2). The positive control α-ketoamide had a stable conformation from 1 to 20 ns, but from 21 ns to the rest of the simulation, there was significant structural distortion. From 1 ns to 20 ns and 40 ns to 50 ns, the bisabololoxide B complex was relatively stable; however, it fluctuated noticeably from 21 ns to 40 ns. Eremanthin remained almost stable throughout the simulation process; only a small amount of structural drifting was found at 8 ns to 10 ns and from 48 ns to 50 ns. The leptospermone complex remained stable from 15 ns to 20 ns and 35 ns to 45 ns but exhibited significant structural instability during the rest of the time. Furthermore, to evaluate the local residual transition, the RMSF benchmark was used. There is a proportional relationship between the RMSF value and the flexibility of a system. If the RMSF value is greater, the system's flexibility will be greater.57 In comparison to the positive control α-ketoamide-Mpro complex (0.21 Å), all three hits had a lower average RMSF value. Eremanthin- Mpro complex showed the lowest average RMSF value (0.15 Å) compared to the other systems. Besides this, both bisabololoxide B-Mpro and leptospermone-Mpro complex displayed an average RMSF value of 0.16 Å (shown in Table 2). Consequently, the lower Rg value evaluates the compactness of the protein structure.58 In this study, bisabololoxide B-Mpro and eremanthin- Mpro complex showed the lowest Rg value (2.13 Å). On the other hand, leptospermone-Mpro complex exhibited the same Rg value (2.16 Å) as the positive control α-ketoamide-Mpro complex (shown in Table 2). A minimal fluctuation rate is always expected during a simulation process. Here also, we calculated the SASA for all systems to determine the compactness of the system (shown in Table 2). Both eremanthin-Mpro (134.31 nm2) and bisabololoxide B-Mpro complex (134.44 nm2) showed lower average SASA value than the positive control α-ketoamide-Mpro complex (139.03 nm2). However, the average SASA value of leptospermone-Mpro complex was 139.75 nm2 which is almost similar to the positive control α-ketoamide-Mpro complex (shown in Table 2). Furthermore, we calculated the total average number of hydrogen bonds for Mpro-ligand complexes, as intermolecular hydrogen bonds between protein-ligand complex plays a vital role in its conformational stability. Here, the positive control α-ketoamide-Mpro complex formed the highest average number of hydrogen bonds (213.79) during the simulation process. Similarly, the eremanthin-Mpro complex produced almost the same number (213.45) of hydrogen bonds as the positive control α-ketoamide-Mpro complex, whereas bisabololoxide B-Mpro and leptospermone-Mpro complex formed an average of 211.03 and 207.13 hydrogen bonds respectively (shown in Table 2).
Discussion
In reality, SARS-CoV-2 is a big concern for the scientific community due to its high infectivity.59 On 30 December 2019 and 30 January 2020, the World Health Organization (WHO) raised warnings and announced this viral infection a ‘public health emergency of international concern’. Since the first report of a coronavirus-related pneumonia outbreak in December 2019, the virus SARS-CoV-2 that causes the infection has emerged into a pandemic, with >100 million individuals infected in over 210 countries and around two million people dying from COVID-19 as of today.60
Nevertheless, no sufficient breakthrough in action against SARS-CoV-2 has been developed yet,61 but researchers are profoundly focused on implementing vaccines and methodologies to eradicate the disease.62,63 However, clinical trials of many candidates are still going on, and some of them pose contentious issues.31,64–66 For instance, An adenosine analogue named Remdesivir which inhibits viral RNA polymerases was known to be an advantageous antiviral agent against a diverse species of RNA viruses along with Ebola, SARS-CoV, and MERS-CoV and exhibited remedial activity in non-clinical models of these coronaviruses but had several undesirable effects in case of SARS-CoV-2 and failed in phase three clinical trials.55,67–69 Failure of hydroxychloroquine in preventing illness of patients with high or moderate-risk exposure to SARS-CoV-2 was also reported previously.70 Therefore, we desperately need potential anti-SARS-CoV-2 drugs to keep people's lives safe.71 New plant-based therapies compensate for around 60% of total Western medicines, and now they are being extensively studied to develop a cure for combating SARS-CoV-2.72,73 This study assayed several plant essential oil compounds through computational approaches to find out effective SARS-CoV-2 Mpro inhibitors.
Since prehistoric periods, plants and their derivatives have been utilized to treat diseases in humans, and they were reported to be clinically safe for humans.74 Plant essential oils are considered to have a range of useful antimicrobial, antibacterial, antiviral, antiparasitic, and insecticidal properties, and now widely studied as a consequential source of drug development for their safety profiles and advantages over synthetic drugs.69,75,76 Some essential oils have shown to be effective against multiple RNA, and DNA viruses, including avian influenza A virus, herpes simplex virus, dengue virus, and poliovirus, as well as are being successful in the treatment of both infectious and chronic diseases.77–81
The protocol of new drug development is complicated, time-intensive, costly, and involves both preclinical and clinical stages. With the aid of computational biology and bioinformatics, manufacturing operations of new drugs are altered and optimized, and hence development costs are minimized.82 This approach focuses on ligands interacting with their target proteins to evaluate the clinical efficacy of specific compounds.83–85
The current study can help to explain the attributes of druglikeness and medicinal chemistry of the key metabolites against SARS-CoV-2. ADME properties estimated by the computational approach provide crucial perceptions into how the body treats or accepts a drug. Computational methods perform a vital role in predicting possible issues with ADME and toxicity and reducing the number of experiments that demand animal trials. Thus, supreme drug candidates were utilized to estimate their drug profiles with ADME analysis.86 In this experiment, all 114 compounds were firstly filtered based on their different ADMET properties. The rule of five and the rule of three were perfectly implemented in 72 compounds out of 114. There were no unfavorable outcomes discovered that could lead to a decrease in their ADME properties and druglikeness. The druglikeness property is assessed by Lipinski's rule of five, whereas Jorgensen's rule of three is more likely to suggest a compound's bioavailability.87,88
Furthermore, when developing new drugs for oral administration, solubility is one of the most significant parameters influencing absorption, and it improves many drug development operations.89,90 The 72 filtered compounds are all water-soluble and have a high rate of oral absorption. Consequently, toxicity tests showed that 53 of the 72 compounds were non-hERG inhibitors and found to be non-Ames toxic and non-carcinogenic in the Ames test and carcinogenic test, respectively.
The 53 screened compounds were then evaluated by molecular docking study. We considered the compounds as a set of ligands in our current analysis and selected the top three hits (Bisabololoxide B, Eremanthin, and Leptospermone) based on their higher binding affinity scores over the other compounds. Results revealed that bisabololoxide B showcased the maximum binding affinity score of -6.6 and the remaining two metabolites (eremanthin and leptospermone) exhibited a binding affinity score of -6.3 kcal/mol, respectively. Bisabololoxide B possesses anti-inflammatory properties, and its pharmacological significance was previously noted.91 Antioxidant properties of eremanthin were previously assayed in STZ induced diabetic rats.92 In an earlier study, eremanthin exhibited notable cytotoxic effects on HeLa cervical cancer cells and hence acted as a promising anti-cancer agent against human cervical cancer.93 Leptospermone is an active compound of the β-triketone portion of manuka oil that exhibited good antiviral effects against HSV-1 and HSV-2 prior to infection. It works by interrupting the structure of the virion envelope and mask viral particles requires for adsorption.94 Leptospermone also showed antimicrobial effects against two bacterial strains, namely C. difficile and C. perfringens in an earlier study.95
Moreover, molecular dynamics simulation was employed to define the exact physical movement and flexibility of ligand-receptor complexes to further evaluate the molecular docking performance.96 RMSD indicates the stability of ligand-protein complexes, where residual flexibility is characterized by RMSF.97 All of our hit compounds showed lower RMSD and RMSF values compared to the positive α-ketoamide control, which depicts their stable conformation with SARS-CoV-2 main protease. The lowest RMSD value was asserted by leptospermone-Mpro complex, which represents its structural compactness. Compared to other systems, eremanthin-Mpro complex displayed a lower RMSF value that denotes its conformational stability. In addition, we also computed the Rg value for each system to determine their compactness. Both bisabololoxide B-Mpro and eremanthin-Mpro complex showed a lower Rg value (2.13 Å), however, the leptospermone-Mpro complex depicted the same Rg value (2.16 Å) as the positive control α-ketoamide-Mpro complex. These outputs suggested that bisabololoxide B-Mpro and eremanthin-Mpro complex were more compact than the positive control α-ketoamide-Mpro complex. Likewise, to evaluate the nature of protein expansion, we estimated the SASA values for all complexes. Here, the eremanthin-Mpro (134.31 nm2) and bisabololoxide B-Mpro complex (134.44 nm2) disclosed a smaller SASA value compared to the positive control α-ketoamide-Mpro complex (139.03 nm2). Consequently, the leptospermone-Mpro complex exhibited a nearly identical but a bit higher SASA value (139.75 nm2) to that of the positive control α-ketoamide-Mpro complex, indicating that the Mpro expanded slightly when it interacted with leptospermone. Furthermore, we calculated the average mean value of total intermolecular hydrogen bonds generated during the simulation process to analyze the structural consistency of the complexes. The bisabololoxide B-Mpro and leptospermone-Mpro complex formed 211.03 and 207.13 hydrogen bonds during the molecular dynamics simulation, while the number of intermolecular hydrogen bonds for eremanthin-Mpro complex (213.45) is nearly equivalent to the positive control α-ketoamide-Mpro complex (213.79). All such outcomes obtained from molecular dynamics simulation indicate the conformational stability and compactness of our selected target compounds with the SARS-CoV-2 main protease.
Throughout the current investigation, we unveiled the molecular interactions of the leading drug candidates with SARS-CoV-2 main protease. The catalytic domain of SARS-CoV-2 main protease owned the ligand-binding sites for each ligand, and His41 and Cys145 residues serve as a site for substrate recognition among the other active residues as it comprises the catalytic dyad.98–100 All of the top three drug candidates were well adjusted into the pocket of SARS-CoV-2 main protease in which several hydrophobic amino acid residues such as Met-49, Gly-143, His-41, Ser-144, His-163, Cys-145, and Met-165 constitute a fairly hydrophobic atmosphere to maintain conformational stabilization of the compounds.100 Thus, based on the findings of this study, we can claim that each of our lead compounds has the potential to act as an inhibitor of SARS-CoV-2, and we hope that the study will pave the way for the development of novel pharmaceutical drugs.
Conclusion
The rapid spreading of novel coronavirus threatens human life terribly, which needs to be stopped for reducing the mortality rate. The main objective of this study was to identify novel inhibitors against the SARS-CoV-2 main protease. To discover novel Mpro inhibitors derived from natural essential oil compounds, ADMET analysis, molecular docking, and MD simulation were successfully conducted here. A total number of 114 essential oil compounds were screened by ADMET profiling and molecular docking study, and the efficacy of the top three compounds was validated by MD simulation. During the 50ns simulation, all the hits (bisabololoxide B, eremanthin, and leptospermone) appeared stable and possessed high affinities against SARS-CoV-2 main protease but further in vivo studies are needed for experimental validation. We conclude that the findings of this study could lead to the development of promising natural pharmaceutical agents for SARS-CoV-2 in the future.
Data availability
Underlying data
Figshare: Underlying data for ‘A computational biology approach for the identification of potential SARS-CoV-2 main protease inhibitors from natural essential oil compounds.’
https://doi.org/10.6084/m9.figshare.16879777.v1.54
This project contains the following underlying data:
• ADME.xlsx (spreadsheet of ADME properties for all the 114 compounds).
• Toxicity.xlsx (spreadsheet of toxicity for screened 72 compounds).
• Binding Affinity.xlsx (spreadsheet of binding affinity for screened 53 compounds).
Data are available under the terms of the Creative Commons Attribution 4.0 International license (CC-BY 4.0).
Extended data
Figshare: Underlying data for ‘A computational biology approach for the identification of potential SARS-CoV-2 main protease inhibitors from natural essential oil compounds.’
https://doi.org/10.6084/m9.figshare.16879777.v1.54
This project contains the following extended data:
• Extended Data 1.docx (document file for the plant source with their corresponding binding affinity of all the 114 compounds).
Data are available under the terms of the Creative Commons Attribution 4.0 International license (CC-BY 4.0).
Authors contribution
RAH and MCA conceived the plan of this research. RAH wrote the manuscript. RAH analyzed the data and made the figures. MSR, MMR, MAI, and FFA edited the manuscript. MCA did the molecular dynamics simulation, and MAHMJ supervised the whole project. All authors revised and approved the manuscript for final submission.
Download
Export To
metrics
| Views | Downloads | |
|---|---|---|
| F1000Research | - | - |
|
PubMed Central
Data from PMC are received and updated monthly.
|
- | - |
Citations
CITE
how to cite this article
Hasib RA, Ali MC, Rahman MS et al. A computational biology approach for the identification of potential SARS-CoV-2 main protease inhibitors from natural essential oil compounds. [version 1; peer review: 1 approved] F1000Research 2021, 10:1313 (https://doi.org/10.12688/f1000research.73999.1)
NOTE: it is important to ensure the information in square brackets after the title is included in all citations of this article.
track
receive updates on this article
Track an article to receive email alerts on any updates to this article.
Open Peer Review
Current Reviewer Status:
Key to Reviewer Statuses
VIEW
HIDE
ApprovedThe paper is scientifically sound in its current form and only minor, if any, improvements are suggested
Approved with reservations
A number of small changes, sometimes more significant revisions are required to address specific details and improve the papers academic merit.
Not approvedFundamental flaws in the paper seriously undermine the findings and conclusions
Version 1
VERSION 1
PUBLISHED 23 Dec 2021
Views
29
How to cite this report:
Nguyen TTB. Reviewer Report For: A computational biology approach for the identification of potential SARS-CoV-2 main protease inhibitors from natural essential oil compounds. [version 1; peer review: 1 approved]. F1000Research 2021, 10:1313 (https://doi.org/10.5256/f1000research.77704.r138870)
The direct URL for this report is:
https://f1000research.com/articles/10-1313/v1#referee-response-138870
https://f1000research.com/articles/10-1313/v1#referee-response-138870
NOTE: it is important to ensure the information in square brackets after the title is included in this citation.
Reviewer Report 31 May 2022
Approved
VIEWS 29
Article "A computational biology approach for the identification of potential SARS-CoV-2 main protease inhibitors from natural essential oil compounds".
The manuscript is clearly presented and well substantiated and I recommend it for indexing in F1000Research with the following minor ... Continue reading
The manuscript is clearly presented and well substantiated and I recommend it for indexing in F1000Research with the following minor ... Continue reading
CITE
HOW TO CITE THIS REPORT
Nguyen TTB. Reviewer Report For: A computational biology approach for the identification of potential SARS-CoV-2 main protease inhibitors from natural essential oil compounds. [version 1; peer review: 1 approved]. F1000Research 2021, 10:1313 (https://doi.org/10.5256/f1000research.77704.r138870)
The direct URL for this report is:
https://f1000research.com/articles/10-1313/v1#referee-response-138870
https://f1000research.com/articles/10-1313/v1#referee-response-138870
NOTE: it is important to ensure the information in square brackets after the title is included in all citations of this article.
-
Author Response 15 Jun 2022Mohammad Abu Hena Mostofa Jamal, Biotechnology and Genetic Engineering, Islamic University, Kushtia, 7003, Bangladesh15 Jun 2022Author ResponseThank you for your insightful feedback to improve our manuscript. Regarding your valuable feedback, we have addressed all the comments individually as follows:
Comment 1: What is the reason ... Continue reading
COMMENTS ON THIS REPORT
-
Author Response 15 Jun 2022Mohammad Abu Hena Mostofa Jamal, Biotechnology and Genetic Engineering, Islamic University, Kushtia, 7003, Bangladesh15 Jun 2022Author ResponseThank you for your insightful feedback to improve our manuscript. Regarding your valuable feedback, we have addressed all the comments individually as follows:
Comment 1: What is the reason ... Continue reading
Comments on this article Comments (0)