Abstract
COVID-19 is one of the members of the coronavirus family that can easily assail humans. As of now, 10 million people are infected and above two million people have died from COVID-19 globally. Over the past year, several researchers have made essential advances in discovering potential drugs. Up to now, no efficient drugs are available on the market. The present study aims to identify the potent phytocompounds from different medicinal plants (Zingiber officinale, Cuminum cyminum, Piper nigrum, Curcuma longa, and Allium sativum). In total, 227 phytocompounds were identified and screened against the proteins S-ACE2 and Mpro through structure-based virtual screening approaches. Based on the binding affinity score, 30 active phytocompounds were selected. Amongst, the binding affinity for beta-sitosterol and beta-elemene against S-ACE2 showed −12.0 and −10.9 kcal/mol, respectively. Meanwhile, the binding affinity for beta-sitosterol and beta-chlorogenin against Mpro was found to be −9.7 and −8.4 kcal/mol, respectively. Further, the selected compounds proceeded with molecular dynamics simulation, prime MM-GBSA analysis, and ADME/T property checks to understand the stability, interaction, conformational changes, binding free energy, and pharmaceutical relevant parameters. Moreover, the hotspot residues such as Lys31 and Lys353 for S-ACE2 and catalytic dyad His41 and Cys145 for Mpro were actively involved in the inhibition of viral entry. From the in silico analyses, we anticipate that this work could be valuable to ongoing novel drug discovery with potential treatment for COVID-19.
Introduction
The newly emerged pandemic disease of COVID-19 was first identified in Wuhan market, China in December 2019 (; ; ). At the same time, a group of people were confirmed as new cases and subsequently died from COVID-19 (; ). The precautions provided by the Chinese government such as intense quarantine, wearing masks, and social distance helped them control the rising rate of new cases (). Even with these precautions, COVID-19 spread throughout the world (). Unluckily, as of February 2021 over 10 million people were said to be infected and over two million people had died. Despite the wide availability of antiviral drugs, an appropriate medicine has not yet been found. Thus, the mortality and morbidity () rates continuously increase every day. Several researchers are currently trying to develop novel therapeutic drugs against COVID-19. Due to the shortcomings of existing drugs, they are not potent against COVID-19. Therefore, an alternative safe and effective antiviral therapy is urgently needed. In order to overcome any shortcomings, plant-based phytocompounds play an important role in treating COVID-19 (). Indian medicinal plants are rich in the source of phytocompounds that exhibit several biological properties including antiviral and immunomodulatory activities (; ). reported that the phytocompounds may exhibit antiviral properties against COVID-19. Moreover, the phytocompounds from various medicinal plants may help to boost the immune response and fight pathogens to combat viral disease (; ; ).
A computational method is a reliable approach to identify active phytocompounds from various databases. This strategy is currently employed in drug discovery research (). Recently, we have reported on the promising phytocompounds identified from different plant sources through the computational approaches for the treatment of COVID-19 (). In the drug discovery process with the help of a computational approach, the selection of structural proteins is one of the major tasks. Two proteins namely S-ACE2 and Mpro have been taken for the present study based on the higher resolution. SARS-CoV-2 contains a spike (S) protein that helps viral entry into the host cell and releases virus particles by attaching to a cell surface receptor called human angiotensin-converting enzyme 2 (ACE2) which facilitates viral transcription and replication (; ; ). It is one of the most common therapeutic drug targets because of its prominent role in host cell binding (). The spike protein of SARS-CoV-2 bound with ACE2 has a higher binding affinity than that of SARS-CoV and ACE2 (). The crystal structure of the S-ACE2 protein contains 603 amino acid residues with a resolution of 2.5Å. Moreover, this structure may be used to investigate the mechanism and binding strength of the viral complexes (). The main protease (Mpro) is the most attractive drug target to combat COVID-19 that can be correlated with viral replication and transcription of the virus life cycle (; ). The Mpro is responsible for the viral polyprotein proteolytic process and viral genome replication and transcription. Hence it is considered as an attractive drug target against COVID-19 (; ). The Mpro has a 1.3Å resolution crystal structure and a sequence similarity of 96% when compared to SARS-CoV. Several computational studies have been performed to identify a potent drug target against COVID-19, they suggested that the Mpro could be a significant drug target based on its binding mechanism and stability ().
In the present study, we selected five medicinal plants: Zingiber officinale (ginger), Cuminum cyminum (cumin), Piper nigrum (black pepper), Curcuma longa (turmeric), and Allium sativum (garlic) based on their different biological applications which include antioxidant, anti-inflammatory, anti-stress, and anti-cancer properties, relief from cardiovascular diseases and respiratory congestion, and anti-diabetic, immune-modulator, anti-microbial, anti-platelet, and cyclo-oxygenase-I (COX-1) inhibitory properties (). From the plants, a total of 227 phytocompounds (33 for Zingiber officinale; 45 for Cuminum cyminum; 42 for Piper nigrum; 54 for Curcuma longa, and 54 for Allium sativum) were screened against various SARS-CoV-2 protein targets (S-ACE2 and Mpro protein). Among them, the top binding affinity phytocompounds such as beta-sitosterol, beta-elemene and beta-chlorogenin were identified. The selected compounds were taken for further assessment using molecular dynamics simulation, ADME/T, and drug likeness properties.
Materials and Methods
The X-ray crystal structures of S-ACE2 with a resolution of 2.45 Å (PDB ID-6M0J) and Mpro with a resolution of 1.31 Å (PDB ID- 5R82) were retrieved from the Protein Data Bank (). The obtained crystal structures were prepared individually with the following steps: water molecules and bound ions were removed, missing hydrogen atoms were added, and Kollman charges were allocated. Then, the energy of the protein structure was minimized using GROMOS 43b1 force field with the help of Swiss-Pdb Viewer. Finally, the structures were saved in the.pdbqt format for molecular docking.
The five selected Indian medicinal plants (Zingiber officinale, Cuminum cyminum, Piper nigrum, Curcuma longa, and Allium sativum) contain 227 phytocompounds (Supplementary Table 1), which were retrieved from the Pubchem database1 and downloaded in the SDF file format ().
Structure-Based Virtual Screening
Structure-based virtual screening (SBVS) plays a vital role in identifying essential phytocompounds from different Indian natural medicinal plants for the development of drug-like compounds. The Vina wizard of the PyRx 8.0 virtual screening tool was employed to screen the phytocompounds against the S-ACE2 and Mpro target proteins (). All the phytocompounds were energy-minimized using the steepest descent gradient algorithm with an MMFF94 force field using the OpenBabel program. Subsequently, the prepared structures were converted into the.pdbqt file format. The 3D grid maps were generated to cover the active site of the target proteins.
Among the 227 phytocompounds, 30 phytocompounds were selected based on their top binding affinity against S-ACE2 and Mpro proteins. The selected 30 phytocompounds are highlighted in Supplementary Table 1. From the above-mentioned 30 compounds, the top three phytocompounds (beta-elemene, beta-sitosterol, and beta-chlorogenin) from each plant which exhibited higher binding affinities against S-ACE2 and Mpro were selected for further assessment. The interaction profiles of protein-ligand complexes were viewed using the Discovery Studio Visualizer 2019 (). The top three compounds were further validated with AutoDock 4.2 and AutoDock Vina tools to compare the binding affinity.
Prime MM-GBSA
Prime MM-GBSA (Molecular Mechanics, the Generalized Born model, and Solvent Accessibility) analysis for S-ACE2 with beta-elemene and beta-sitosterol; Mpro with beta-sitosterol and beta-chlorogenin complexes was performed to calculate the binding free energies using Schrodinger, LLC, NY (; ). The following equation was used to enumerate the binding free energy.
The Gcomplex indicates complex energy, Gprotein indicates the receptor energy, and Gligand indicates the unbound ligand energy.
ADME/T (Absorption, Distribution, Metabolism, Elimination, and Toxicity) Prediction
After the SBVS of the 227 phytocompounds against target proteins, based on the binding energy values, the phytocompounds, beta-elemene, beta-sitosterol, and beta-chlorogenin, were considered for ADME/T properties using the pkCSM server2 to predict their pharmacokinetic properties ().
Molecular Dynamics Simulation and Post-MM/GBSA Analyses
Furthermore, to evaluate the structural stability of the S-ACE2 with beta-elemene and beta-sitosterol; Mpro with beta-sitosterol and beta-chlorogenin complexes, Molecular dynamics simulation (MDS) was performed using the GROMACS package with a GROMOS96 43a1 force field (). The topology and parameters for the ligands were obtained from the PRODRG server3. The protein-ligand complexes were equilibrated under cubic periodic boundary conditions with the dimensions of 10 nm × 10 nm × 10 nm and solvated with an explicit SPC (simple point charge) water model. The proper counter ions Na+Cl– were added to neutralize the system. The energy minimization steps were performed using the steepest descent gradient algorithm to remove the weak van der Waals contacts. The NPT and NVT ensembles were employed for 50,000 steps in 100 ps, respectively. Then a constant temperature was maintained at 300 K with a coupling time of 0.1 ps. The constant pressure was also maintained with 1 bar for 100 ps. Further, the MDS was carried out for 100 ns for all the protein-ligand complexes. The trajectories were analyzed for docked complexes to calculate the root mean square deviation (RMSD), root mean square fluctuation (RMSF), and hydrogen bond using scripts included in the GROMACS package. Each trajectory was graphically visualized using the XMGRACE 2D plotting tool. In addition, the post MM/GBSA analysis was performed to calculate the binding energy using the prime MM/GBSA module in Schrodinger.
Results and Discussion
Structure-Based Virtual Screening
Firstly, the co-crystallized compound and known binders such as Lopinavir and Rilapladib were re-docked against S-ACE2 and Mpro. The co-crystallized compound of S-ACE2 showed a binding affinity of −6.9 kcal/mol, and a binding affinity of −7.0 kcal/mol for Mpro. The known binder result reveals that Lopinavir and Rilapladib against S-ACE2 showed −5.4 and −6.0 kcal/mol, respectively. On the other hand, Lopinavir and Rilapladib against Mpro showed −7.5 and −6.7 kcal/mol, respectively.
Based on the binding affinity of the above reported compounds, our study aimed to investigate the potential therapeutic drug candidates to combat COVID-19. The identified 227 phytocompounds from five different natural medicinal plants were screened against target proteins such as S-ACE2 and Mpro. Table 1 represents the binding affinity of the top phytocompounds against S-ACE2 and Mpro.
TABLE 1
| S-ACE2 (PDB ID—6M0J) | |||||||||
| Zingiber officinale—ginger | Cuminum cyminum—cumin | Piper nigrum—black pepper | Curcuma longa—turmeric | Allium sativum—garlic | |||||
| Compound name | Binding affinity (Kcal/mol) | Compound name | Binding affinity (Kcal/mol) | Compound name | Binding affinity (Kcal/mol) | Compound name | Binding affinity (Kcal/mol) | Compound name | Binding affinity (Kcal/mol) |
| Alpha-selinene | −10.9 | Beta-elemene | −10.9 | Beta-selinene | −7.8 | Beta-sitosterol | −12.0 | Saponin | −10.6 |
| Beta-sitosterol | −12.0 | Alpha-selinene | −10.9 | Beta-elemene | −10.9 | Cyclocurcumin | −9.1 | Beta-tocopherol | −10.5 |
| Beta-elemene | −10.9 | Beta-sitosterol | −12.0 | Beta-sitosterol | −12.0 | Demethoxycurcumin | −8.9 | Beta-chlorogenin | −8.4 |
| Mpro (PDB ID—5R82) | |||||||||
| Beta-sitosterol | −9.7 | Beta-sitosterol | −9.7 | Beta-sitosterol | −9.7 | Beta-sitosterol | −9.7 | Beta-chlorogenin | −8.4 |
| Kaempferol | −7.5 | Apigetrin | −8.0 | Alpha-tocopherol | −8.0 | Alpha-tocopherol | −8.0 | Apigenin | −7.5 |
| Alpha-selinene | −7.0 | Alpha-selinene | −7.2 | Kaempferol | −7.5 | Riboflavin | −7.2 | Kaempferol | −7.5 |
Binding affinity for the top phytoconstituents against S-ACE2 and Mpro proteins.
Based on their higher binding affinity, beta-sitosterol (−12.0 kcal/mol) and beta-elemene (−10.9 kcal/mol) against S-ACE2 and beta-sitosterol (−9.7 kcal/mol) and beta-chlorogenin (−8.4 kcal/mol) against Mpro were taken for further studies. Furthermore, the AutoDock and AutoDock Vina tools provided a close binding affinity score against S-ACE2 and Mpro when compared with the PyRx tool (Table 2).
TABLE 2
| Compound name | AutoDock (Kcal/mol) | AutoDock Vina (Kcal/mol) |
| S-ACE2 (6M0J) | ||
| Beta-sitosterol | −12.3 | −11.09 |
| Beta-elemene | −9.8 | −9.3 |
| Mpro (5R82) | ||
| Beta-sitosterol | −10.01 | −8.9 |
| Beta-chlorogenin | −8.8 | −7.1 |
Docking score for the top phytoconstituents against S-ACE2 and Mpro proteins.
The Binding Affinities of the Phytocompounds Into the S-ACE2 Active Site
The binding affinity for S-ACE2 with the beta-elemene and beta-sitosterol complexes were analyzed and shown in Figures 1A,B. Figure 1A displays the hydrogen bond, hydrophobic, and van der Waals interactions. Amino acid residues such as Glu37 and Arg403 actively participated in hydrogen bond interactions. Residues like Lys26, His34, Val93, Pro389, and Phe456 were actively involved with hydrophobic interactions. In van der Waals interactions, the Glu23, Thr27, Asp30, Asn33, Gln96, Ala387, Gln388, Phe390, Arg393, and Tyr505 residues actively interacted with beta-sitosterol. Polar amino acid residues such as Thr27, Gln96, and Tyr505 contributed to the interactions. Figure 1B reveals the hydrophobic and van der Waals interactions. Hydrophobic interactions involved the His34 residue. The van der Waals interaction was influenced by the Glu35, Glu37m and Asp38 residues along with hotspot residues Lys31 and Lys353. A recent study reports that the beta-elemene phytocompound can interact with the hotspot residues of Lys31 and Lys353 which play a crucial role for the inhibition of viral entry (). The interacting residues of these complexes are shown in Table 3.
FIGURE 1
TABLE 3
| Phytocompound name | H-bond | Hydrophobic interaction | Van der Waals interaction |
| Beta-sitosterol | Glu37 and Arg403 | Lys26, His34, Val93, Pro389, and Phe456 | Glu23, Thr27, Asp30, Asn33, Gln96, Ala387, Gln388, Phe390, Arg393, and Tyr505 |
| Beta-elemene | - | His34 | Glu35, Glu37, Asp38, Lys31, and Lys353 |
Interaction residues of beta-sitosterol and beta-elemene with S-ACE2 complexes.
The bolded Lys31 and Lys353 are the important hotspot residues for S-ACE2, and His41 and Cys145 are the essential catalytic dyads for Mpro. They are actively involved in the inhibition of viral entry.
Furthermore, the beta-sitosterol and beta-elemene were docked with S-ACE2 using the AutoDock (Figures 1C,D) and AutoDock Vina (Figures 1E,F) tools. The results showed that both the beta-sitosterol and beta-elemene complexes had similar binding affinity and interactions.
Generally, the mediating region of the S-protein of SARS-CoV-2 is highly responsible for contact with the ACE2 receptor through the receptor binding domain (RBD) which is present in host cells. SARS-CoV-2 enters into human cells through the ACE2 receptor, where targeting the receptor provides more insights on controlling SARS-CoV-2 infections. Thus, structure-based drug design (SBDD) approaches have been employed to identify the potential phytocompounds which could prevent the binding with the SARS-CoV-2 ACE2 receptor (). According to the structure-based virtual screening (SBVS) results, this study provides selected phytocompounds may be potential inhibitors against SARS-CoV-2.
The Binding Affinities of the Phytocompounds Into Catalytic Dyad of Mpro
Figures 2A,B shows the binding affinity of the Mpro complexed with beta-sitosterol and beta-chlorogenin. As shown in Figure 2A, beta-sitosterol had one hydrogen bond interaction with Thr24 and two hydrophobic interactions with residues Cys145 and Phe168. Also, the study observed that the van der Waals interactions involved the following residues: Thr25, Thr26, Leu27, His41, Ser46, Met49, Asn142, Gly143, Ser144, His163, Met165, Glu166, Leu167, and Gln189. As shown in Figure 2B, beta-chlorogenin formed two hydrogen bond interactions (Thr25 and Cys145 residues) and three hydrophobic interactions (Cys145, His163, and His172 residues). Van der Waals interactions with the residues Thr24, Thr25, Ser46, Phe140, Gly143, Ser144, His164, and Met165 were noticed. Polar amino acids such as Thr25 and Cys145 were influenced. The interaction residues of these complexes are shown in Table 4.
FIGURE 2
TABLE 4
| Phytocompound name | H-bond | Hydrophobic interaction | Van der Waals interaction |
| Beta-sitosterol | Thr24 and Thr26 | Cys145 and Pro168 | Thr25, Leu27, His41, Ser46, Met49, Asn142, Gly143, Ser144, His163, Met165, Glu166, Leu167, and Gln189 |
| Beta-chlorogenin | Thr25 and Cys 145 | Cys145, His163, and His172 | Thr24, Thr25, Ser46, Phe140, Gly143, Ser144, His164, and Met165 |
Interaction residues of beta-sitosterol and beta-chlorogenin with Mpro complexes.
The bolded Lys31 and Lys353 are the important hotspot residues for S-ACE2, and His41 and Cys145 are the essential catalytic dyads for Mpro. They are actively involved in the inhibition of viral entry.
In earlier reports, the catalytic dyad Cys145 and His41 residues of Mpro actively participated in the binding regions (; ). In our study, the aforementioned selected phytocompounds actively interacted with the binding pocket of the catalytic dyad (Cys145 and His41) of the substrate binding domain of the Mpro. The substantial binding affinity and formation of strong interactions between the above residues indicated the stability of the docked complexes.
The Mpro with beta-sitosterol and beta-chlorogenin complexes were further analyzed using AutoDock (Figures 2C,D) and AutoDock Vina (Figures 2E,F) tools. These showed similar binding affinity and interactions. Therefore, the molecular docking analysis suggested that beta-sitosterol and beta-chlorogenin play a vital role in replication and pathogenesis. Hence, it may inhibits COVID-19 Mpro enzyme activity which is taken further in molecular dynamics simulations.
Binding Free Energy
To evaluate the binding free energy of the complexes, the MM/GBSA approach was employed (Table 5). The results revealed that S-ACE2 with the beta-elemene and beta-sitosterol complexes showed binding free energies of −28.87 and −26.17 kcal/mol, respectively. The binding free energy of Mpro with the beta-sitosterol and beta-chlorogenin complexes were found to be −33.21 and −10.59 kcal/mol, respectively. The overall binding free energies indicate that the lead molecules such as beta-elemene, beta-sitosterol, and beta-chlorogenin may strongly bind in the binding region of taken targets to inhibit enzymatic activity as well as the innate immune response of the human body to prevent COVID-19.
TABLE 5
| S.No | Compound name | ΔGbind | ΔGcovalent | ΔGcoulomb | ΔGHBond | ΔGLipo | ΔGSolvGB | ΔGvdW |
| S-ACE2 (6M0J) | ||||||||
| 1 | Beta-sitosterol | −28.87 | 2.95 | −1.10 | −1.18 | −11.78 | 25.09 | −33.04 |
| 2 | Beta-elemene | −26.17 | 4.88 | −1.54 | −0.81 | −12.01 | 28.59 | −33.00 |
| Mpro (5R82) | ||||||||
| 3 | Beta-sitosterol | −33.21 | 5.57 | −5.28 | −1.01 | −15.17 | 21.19 | −38.50 |
| 4 | Beta-chlorogenin | −40.59 | 1.93 | −9.35 | −0.47 | −15.47 | 17.90 | −35.20 |
Binding free energy of top phytocompounds against S-ACE2 and Mpro.
ADME/T Properties
The efficiency of the therapeutic phytocompounds mainly depends on their bioactivity and pharmacokinetics properties. The ADME/T properties of beta-elemene, beta-sitosterol, and beta-chlorogenin were computed using the PkCSM web server, and are shown in Table 6. The selected phytocompounds showed an acceptable range including < 500 molecular weight, < 10 hydrogen bond acceptors, < 5 hydrogen bond donors, and logP values of < 5. Furthermore, the pharmacokinetic parameters such as human oral absorption, partition co-efficient (QPlogPo/w), and water solubility (QPlogS) values did not violate Lipinski’s RO5. Hence, the observed values for the selected phytocompounds can be considered as potent inhibitors for COVID-19.
TABLE 6
| Properties | Beta-elemene | Beta-sitosterol | Beta-chlorogenin |
| Molecular weight | 204.357 | 414.718 | 432.645 |
| TPSA | 94.774 | 187.039 | 188.008 |
| LogP | 4.7472 | 8.0248 | 4.7646 |
| H-bond acceptor | 0 | 1 | 4 |
| H-bond donor | 0 | 1 | 2 |
| Absorption | |||
| Water solubility (log mol/L) | −6.43 | −6.773 | −5.213 |
| Caco2 permeability (log Papp in10-6 cm/s) | 1.41 | 1.201 | 1.263 |
| Intestinal absorption (human) (% absorbed) | 94.359 | 94.464 | 96.823 |
| Skin permeability (log Kp) | −1.279 | −2.783 | −3.999 |
| P-Glycoprotein substrate | No | No | Yes |
| P-Glycoprotein I inhibitor | No | Yes | Yes |
| P-Glycoprotein II inhibitor | No | Yes | Yes |
| Distribution | |||
| VDss (human, log L/kg) | 0.601 | 0.193 | 0.192 |
| Fraction unbound (human) (Fu) | 0.157 | 0 | 0.037 |
| BBB permeability (logBB) | 0.809 | 0.781 | 0.004 |
| CNS permeability (log PS) | −1.714 | −1.705 | −1.592 |
| Metabolism | |||
| CYP2D6 substrate | No | No | No |
| CYP3A4 substrate | No | Yes | Yes |
| CYP1A2 inhibitor | No | No | No |
| CYP2C19 inhibitor | No | No | No |
| CYP2C9 inhibitor | No | No | No |
| CYP2D6 inhibitor | No | No | No |
| CYP3A4 inhibitor | No | No | No |
| Excretion | |||
| Total clearance (log ml/min/kg) | 0.251 | 0.628 | 0.346 |
| Renal OCT2 substrate | No | No | Yes |
| Toxicity | |||
| AMES toxicity | No | No | No |
| hERG I inhibitor | No | No | No |
| hERG II inhibitor | No | Yes | No |
| Hepatotoxicity | No | No | No |
| Skin sensitization | Yes | No | No |
ADME/T properties of the top phytocompounds using the PkCSM web server.
In addition, the phytocompounds displayed negative AMES toxicity and negative carcinogenicity which indicated that they are non-mutagenic and non-carcinogenic. From the results, it was noticed that the selected phytocompounds had an acceptable range except beta-sitosterol. It was because the LogP value was 8, hence beta-sitosterol had violated the drug-likeliness properties. Several reports suggested that the beta-sitosterol phytocompound possessed anti-HIV and anti-HBV (hepatitis B virus) activities (; ). In addition to this, beta-sitosterol has shown activity against SARS-CoV-2. An earlier study suggested that with the availability of beta-sitosterol in the drug industry, this molecule could be considered as a potential inhibitor of SARS-CoV-2 (; ). In line with this, beta-sitosterol has been considered to have potential antiviral activity against COVID-19.
Molecular Dynamics Simulations
A molecular dynamics simulation was performed to evaluate the stability and conformational changes of the docked complex. Figures 3A,B) represents the MDS of S-ACE2 with the beta-sitosterol and beta-elemene complexes; Mpro with the beta-sitosterol and beta-chlorogenin complexes, respectively. To examine the degree of conformational changes in the protein (backbone) and docked complexes, RMSD was calculated for 100 ns simulation trajectories for each complex. The RMSD of the backbone of S-ACE2 showed stability throughout the simulation (indicated in black) by maintaining a peak of 0.3 nm. S-ACE2 with the beta-sitosterol (indicated in red) and beta-elemene (indicated in green) complexes also occurred in the acceptable range and no major fluctuations were observed when compared to the backbone (Figure 3A). The backbone atoms of Mpro (indicated in black) maintained a peak with equilibration ranges between 0.25 and 0.4 nm and showed better stability until the end of the simulation. On the other hand, the Mpro with beta-sistosterol showed a similar acceptable range from 0.3 to 0.45 nm (indicated in red). The beta-chlorogenin equilibrated at 5 ns, fluctuated < 0.1 nm (on an average), and showed constant peaks throughout the simulation (indicated in green). Hence, all the three converged at 40 ns and were constant at the end of the simulation (Figure 3B). RMSD of the complexes were calculated as a function of time period to evaluate the conformational stability of the phytocompounds (). The trajectory analysis revealed that the protein-ligand complexes were stable and active in dynamic movement after attaining the equilibrium state. From the obtained results, it was observed that the taken complexes were stable throughout the simulation.
FIGURE 3
The RMSF graph provides information about the residual fluctuation of the protein-ligand complex. Figure 4 illustrates the RMSF of (A) S-ACE2 with beta-sitosterol and beta-elemene complexes; (B) Mpro with beta-sitosterol and beta-chlorogenin complexes for the period of 100 ns MD simulations. The average RMSF was maintained at 0.33 nm for S-ACE2 and 0.28 nm for Mpro complexes. For S-ACE2 with beta-sitosterol and beta-elemene complexes, the residues fluctuated between 100 and 120, 300, and 320. The following residues such as Asp615 (0.5 nm), Asp136 (0.6 nm), Asn338 (0.4 nm), and Val339 (0.4 nm) had a higher fluctuation with beta-sitosterol. Subsequently, the Asp615 (0.7 nm), Tyr613 (0.5 nm), Asn338 (0.4 nm), and Ala614 (0.4 nm) residues fluctuated in beta-elemene with the S-ACE2 complex. In the Mpro with beta-sitosterol and beta-chlorogenin complexes, fluctuations were observed between the amino acids 200 and 225. The Mpro with beta-sitosterol fluctuated for the following residues: Thr304 (0.6 nm), Phe223 (0.4 nm), Gln244 (0.4 nm), and Val303 (0.4 nm). The fluctuations for the Mpro with beta-chlorogenin complex were observed for residues Thr304 and Val303 (0.5 nm). The fluctuation of amino acids may be involved in the loop region which was monitored in the molecular dynamics simulations. From the outcome of the RMSF analysis, it was observed that the complexes were stable and the fluctuated residues did not have any impact on the protein-ligand complex.
FIGURE 4
Hydrogen bond interactions are the most important to determine the bonding strength between the protein-ligand complex. Figures 5A,B) indicates the number of hydrogen bonds formed during the simulation. The hydrogen bonds for S-ACE-2 with beta-sitosterol and beta-elemene complexes; Mpro with beta-sitosterol and beta-chlorogenin complexes increased during the simulations which indicated stable binding. Overall the simulation analysis indicated that no conformational changes were observed during the simulation. Therefore, the identified phytocompounds could be promising candidates for inhibiting S-ACE2 and Mpro.
FIGURE 5
Post MM/GBSA Analyses
The post MM/GBSA of the binding free energy calculation was analyzed with the generation of four-frame time intervals (25, 50, 75, and 100 ns) throughout the MDS as shown in Table 7. The results showed that all the complexes had better binding free energy at 100 ns when compared with the pre-MM/GBSA binding free energy.
TABLE 7
| S.No | Compound name | ΔGbind (25 ns) | ΔGbind (50 ns) | ΔGbind (75 ns) | ΔGbind (100 ns) |
| S-ACE2 (6M0J) | |||||
| 1 | Beta-sitosterol | −16. 09 | −19.09 | −23.96 | −31.48 |
| 2 | Beta-elemene | −12.24 | −19.98 | −21.90 | −28.09 |
| Mpro (5R82) | |||||
| 3 | Beta-sitosterol | −9.24 | −16.55 | −29.09. | −36.01 |
| 4 | Beta-chlorogenin | −13.00 | −20.24 | −29.99 | −45.77 |
Post MM/GBSA of top phytocompounds against S-ACE2 and Mpro.
Conclusion
Traditionally, drugs from medicinal plant sources have been widely used to control disease. At this critical stage of the COVID-19 infection, an effective drug is urgently needed. Plants contain a rich source of phytocompounds which may be an effective approach to combat COVID-19. The present study revealed the better binding affinities of beta-sitosterol, beta-elemene, and beta-chlorogenin in the active site of the S-ACE2 and Mpro proteins. In addition, MDS revealed that the amino acid fluctuation may be involved in the loop region that exhibits stability throughout the simulation, and hydrogen bonds for all the complexes increased during the simulations. Furthermore, the post-MM/GBSA showed that all the complexes had better binding free energy at 100 ns when compared with pre-MM/GBSA binding free energy. Moreover, the ADME/T properties confirmed that the identified phytocompounds could be considered as promising drug-like compounds. The identified phytocompounds (beta-sitosterol, beta-elemene, and beta-chlorogenin) are publicly available which will facilitate the rapid development of suitable effective therapeutic candidates to control viral replication as well as the innate immune response to combat COVID-19. Therefore, the current study will be highly useful for further research to design specific drugs against COVID-19. In future, the identified phytocompounds will be extended through experimental studies to confirm their status as novel compounds against COVID-19x.
Statements
Data availability statement
The original contributions presented in the study are included in the article/Supplementary Material, further inquiries can be directed to the corresponding author/s.
Author contributions
MS: investigation, methodology, and writing—original draft. BP: conceptualization, methodology, writing—review and editing, and supervision. NM: analyzing the result and discusses the manuscript. BR: discussing the in silico work and writing the manuscript. VR: editing the manuscript. PP: methodology design and editing. GS: developed the theoretical frame work. YW: review and editing, and designed the figure. BM: aided in interpreting the results and worked on the manuscript. BC: review and editing, and supervision. MB: review and editing the manuscript. SR: contributed to the interpretation of the results. JJ: contributed to the final version of the manuscript. GH: writing—review and editing, and supervision. SE: analyze the result and writing the manuscript. All authors contributed to the article and approved the submitted version.
Acknowledgments
The authors thankfully acknowledge the UGC-Innovative [No. F. 14-13/2013 (Inno/ASIST)], DST-FIST [SR/FST/LSI- 667/2016(C)], DST PURSE [SR/PURSE Phase 2/38 (G)], and MHRD-RUSA 2.0 [F.24/51/2014-U, Policy (TNMulti-Gen), Department of Education Government of India] for financial support and infrastructure facilities.
Conflict of interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Supplementary material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fmolb.2021.637122/full#supplementary-material
References
1
AdejoroI. S.BabatundeD. D.TolufasheG. T. (2020). Molecular docking and dynamic simulations of some medicinal plants compounds againstSARS-CoV-2: an in silico study.Taibah Univ. Sci.141563–1570. 10.1080/16583655.2020.1848049
2
AlexpandiR.De MesquitaJ. F.PandianS. K.RaviA. V. (2020). Quinolines-based SARS-CoV-2 3CLpro and RdRp Inhibitors and Spike-RBD-ACE2 Inhibitor for Drug-Repurposing Against COVID-19: an in silico analysis.Front. Microbiol.11:1796.
3
Arul MuruganN.ChitraJ. P.JeyakanthanJ. (2020). Computational investigation on Andrographis paniculata phytochemicals to evaluate their potency against sars-cov-2 in comparison to known antiviral compounds in drug trials.J. Biomol. Struct. Dyn.161–12. 10.1080/07391102.2020.1777901
4
BabichO.SukhikhS.ProsekovA.AsyakinaL.IvanovaS. (2020). Medicinal Plants to Strengthen Immunity during a Pandemic.Pharmaceuticals13:313. 10.3390/ph13100313
5
BaiC.WarshelA. (2020). Critical differences between the binding features of the spike proteins of SARS-CoV-2 and SARS-CoV.J. Phys. Chem. B1245907–5912. 10.1021/acs.jpcb.0c04317
6
BaigA. M.KhaleeqA.SyedaH. (2020). Docking prediction of amantadine in the receptor binding domain of spike protein of SARS-CoV-2.ACS Pharmacol. Transl. Sci31430–1433. 10.1021/acsptsci.0c00172
7
BattistiV.WiederO.GaronA.SeidelT.UrbanE.LangerT. (2020). A computational approach to identify potential novel inhibitors against the coronavirus SARS-CoV-2.Mol. Inform.39:2000090. 10.1002/minf.202000090
8
BenarbaB.PandiellaA. (2020). Medicinal plants as sources of active molecules against COVID-19.Front. Pharmacol.11:1189.
9
Ben-ShabatS.YarmolinskyL.PoratD.DahanA. (2020). Antiviral effect of phytochemicals from medicinal plants: applications and drug delivery strategies.Drug Deliv. Transl. Res.10354–367. 10.1007/s13346-019-00691-6
10
BharadwajS.LeeK. E.DwivediV. D.KangS. G. (2020). Computational insights into tetracyclines as inhibitors against SARS-CoV-2 Mpro via combinatorial molecular simulation calculations.Life Sciences.257118080. 10.1016/j.lfs.2020.118080
11
CavasottoC.Di FilippoJ. (2020). In silico drug repurposing for COVID-19: targeting SARS-CoV-2 proteins through docking and consensus ranking.Mol. Inform.40:2000115. 10.1002/minf.202000115
12
CavasottoC. N.LamasS. M.MagginiJ. (2021). Functional and druggability analysis of the SARS-CoV-2 proteome.Eur. J. Pharmacol.890:173705. 10.1016/j.ejphar.2020.173705
13
Chandra BabuT. M.RajeshS. S.BhaskarB. V.DeviS.RammohanA.SivaramanT.et al (2017). Molecular docking, molecular dynamics simulation, biological evaluation and 2D QSAR analysis of flavonoids from Syzygium alternifolium as potent anti-Helicobacter pylori agents.RSC Adv.718277–18292. 10.1039/c6ra27872h
14
ChengV. C.WongS. C.ChuangV. W.SoS. Y.ChenJ. H.SridharS.et al (2020). The role of community-wide wearing of face mask for control of coronavirus disease 2019 (COVID-19) epidemic due to SARS-CoV-2.J. Infect.81107–114. 10.1016/j.jinf.2020.04.024
15
ChikhaleR. V.GuravS. S.PatilR. B.SinhaS. K.PrasadS. K.ShakyaA.et al (2020). Sars-cov-2 host entry and replication inhibitors from Indian ginseng: an in-silico approach.J. Biomol. Struct. Dyn.20201–12. 10.1080/07391102.2020.1778539
16
ChirumboloS. (2012). Plant phytochemicals as new potential drugs for immune disorders and cancer therapy: really a promising path?J. Sci. Food Agric.921573–1577. 10.1002/jsfa.5670
17
ChoudharyS.MalikY. S.TomarS. (2020). Identification of SARS-CoV-2 cell entry inhibitors by drug repurposing using in silico structure-based virtual screening approach.Front. Immunol.11:1664. 10.3389/fimmu.2020.01664
18
ChowdhuryP. (2019). In silico investigation of phytoconstituents from Indian medicinal herb ‘Tinospora cordifolia (giloy)’ against SARS-CoV-2 (COVID-19) by molecular dynamics approach.J. Biomol. Struct. Dyn.71–18. 10.1080/07391102.2020.1803968
19
FerrazW. R.GomesR. A.NovaesA. L. S.Goulart TrossiniG. H. (2020). Ligand and structure-based virtual screening applied to the SARS-CoV-2 main protease: an in silico repurposing study.Future Med. Chem.201815–1828. 10.4155/fmc-2020-0165
20
GaoJ.ZhangL.LiuX.LiF.MaR.ZhuZ.et al (2020). Repurposing low–molecular-weight drugs against the main protease of severe acute respiratory syndrome coronavirus 2.J. Phys. Chem. Lett117267–7272. 10.1021/acs.jpclett.0c01894
21
GowrishankarS.MuthumanickamS.KamaladeviA.KarthikaC.JothiR.BoomiP.et al (2021). Promising phytochemicals of traditional Indian herbal steam inhalation therapy to combat COVID-19 – an in silico study.Food Chem. Toxicol.148:111966. 10.1016/j.fct.2020.111966
22
GuoR.ZhangY.LiX.SongX.LiD.ZhaoY. (2016). Discovery of ERBB3 inhibitors for non-small cell lung cancer (NSCLC) via virtual screening.J. Mol. Model22:135. 10.1007/s00894-016-3007-z
23
HalderT.DhorajiwalaT. M.SamantL. R. (2019). Multiple docking analysis and In silico absorption, distribution, metabolism, excretion, and toxicity screening of anti-leprosy phytochemicals and dapsone against dihydropteroate synthase of Mycobacterium leprae.Int. J. Mycobacteriol.8229–236. 10.4103/ijmy.ijmy_123_19
24
HuynhT.WangH.LuanB. (2020). In silico exploration of the molecular mechanism of clinically, oriented drugs for possibly inhibiting SARS-CoV-2′s main protease.J. Phys. Chem. Lett.114413–4420. 10.1021/acs.jpclett.0c00994
25
KhannaK.KohliS. K.KaurR.BhardwajA.BhardwajV.OhriP.et al (2020). Herbal immune-boosters: substantial warriors of pandemic Covid-19 battle.Phytomedicine85:153361. 10.1016/j.phymed.2020.153361
26
Linh NguyenH.Quoc ThaiN.Toan TruongD.Suan LiM. (2020). Remdesivir strongly binds to both RNA-Dependent RNA polymerase and main protease of SARS-CoV-2: evidence from molecular simulations.J. Phys. Chem. B12411337–11348. 10.1021/acs.jpcb.0c07312
27
Meyer-AlmesF.-J. (2020). Repurposing approved drugs as potential Inhibitors of 3CL-protease of SARS-CoV-2: virtual screening and structure based drug design.Comput. Biol. Chem.88:107351. 10.1016/j.compbiolchem.2020.107351
28
MukhtarM.ArshadM.AhmadM.PomerantzR. J.WigdahlB.ParveenZ. (2008). Antiviral potentials of medicinal plants.Virus Res.131111–120. 10.1016/j.virusres.2007.09.008
29
MuthumanickamS.IndhumathiT.BoomiP.BalajeeR.JeyakanthanJ.AnandK.et al (2020). In silico approach of naringin as potent phosphatase and tensin homolog (PTEN) protein agonist against prostate cancer.J. Biomol. Struct. Dyn.91–10. 10.1080/07391102.2020.1830855
30
ParvezM. K.RehmanM. T.AlamP.Al-DosariM. S.AlqasoumiS. IAlajmiM. F. (2019). Plant-derived antiviral drugs as novel hepatitis B virus inhibitors: cell culture and molecular docking study.Saudi Pharm. J.27:389. 10.1016/j.jsps.2018.12.008
31
QiF.QianS.ZhangS.ZhangZ. (2020). Single cell RNA sequencing of 13 human tissues identify cell types and receptors of human coronaviruses.Biochem. Biophys. Res. Commun.526135–140. 10.1016/j.bbrc.2020.03.044
32
RajagopalK.ByranG.JupudiS.VadivelanR. (2020). Activity of phytochemical constituents of black pepper, ginger, and garlic against coronavirus (COVID-19): an in silico approach.Int. J. Health Allied Sci.943–50.
33
RamachandranB.SrinivasadesikanV.ChouT.-M.JeyakanthanJ.LeeS.-L. (2020). Atomistic simulation on flavonoids derivatives as potential inhibitors of bacterial gyrase of Staphylococcus aureus.J. Biomol. Struct. Dyn.1–14. 10.1080/07391102.2020.1856184
34
RowaiyeA. B.OnuhO. A.SundayR. M.AbdulmalikZ. D.BurD.EmeterN. W.et al (2020). Structure-based virtual screening and molecular dynamic simulation studies of the natural inhibitors of SARS-CoV-2 main protease.J. Chem. Res.520–31.
35
SchoemanD.FieldingB. C. (2019). Coronavirus envelope protein: current knowledge.Virol. J.16:69.
36
SelvarajC.SinghS. K. (2014). Validation of potential inhibitors for SrtA against Bacillus anthracis by combined approach of ligand-based and molecular dynamics simulation.J. Biomol. Struct. Dyn.321333–1349. 10.1080/07391102.2013.818577
37
ShyrZ. A.GorshkovK.ChenC. Z.ZhengW. (2020). Drug discovery strategies for SARS-CoV-2.J. Pharmacol. Exp. Ther.375127–138. 10.1124/jpet.120.000123
38
TanakaT.IkedaT.KakuM.ZhuX. H.OkawaM.YokomizoK.et al (2004). New lignan glycoside and phenylethanoid glycosides from Strobilanthes cusia BREMEK.Chem. Pharm. Bull.521242–1245. 10.1248/cpb.52.1242
39
TariqA.MateenR. M.AfzalM. S.SaleemM. (2020). Paromomycin: a potential dual targeted drug effectively inhibits both spike (S1) and main protease of COVID-19.Int. J. Infect. Dis.98166–175. 10.1016/j.ijid.2020.06.063
40
TouretF.GillesM.BarralK.NougairèdeA.van HeldenJ.DecrolyE.et al (2020). In vitro screening of a FDA approved chemical library reveals potential inhibitors of SARS-CoV-2 replication.Sci. Rep.10: 70143.
41
WallsA. C.ParkY.-J.TortoriciM. A.WallA.McGuireA. T.VeeslerD. (2020). Structure, function, and antigenicity of the SARS-CoV-2 spike glycoprotein.Cell181281–292.e6.
42
WangJ. (2020). Fast identification of possible drug treatment of coronavirus disease-19 (COVID-19) through computational drug repurposing study.J. Chem. Inf. Model603277–3286. 10.1021/acs.jcim.0c00179
43
WhisenantJ.BurgessK. (2019). Blocking Coronavirus 19 Infection via the SARS-CoV-2 spike protein: initial steps.ACS Med. Chem. Lett.20201076–1078. 10.1021/acsmedchemlett.0c00233
Summary
Keywords
COVID-19, main protease, molecular dynamics simulation, natural medicinal plants, S-ACE2, structure-based virtual screening
Citation
Sankar M, Ramachandran B, Pandi B, Mutharasappan N, Ramasamy V, Prabu PG, Shanmugaraj G, Wang Y, Muniyandai B, Rathinasamy S, Chandrasekaran B, Bayan MF, Jeyaraman J, Halliah GP and Ebenezer SK (2021) In silico Screening of Natural Phytocompounds Towards Identification of Potential Lead Compounds to Treat COVID-19. Front. Mol. Biosci. 8:637122. doi: 10.3389/fmolb.2021.637122
Received
02 December 2020
Accepted
16 April 2021
Published
05 July 2021
Volume
8 - 2021
Edited by
Matteo Becatti, University of Firenze, Italy
Reviewed by
Bruce David Uhal, Michigan State University, United States; Claudio Norberto Cavasotto, Universidad Austral, Argentina; Kemal Yelekci, Kadir Has University, Turkey; Mahmoud Kandeel, King Faisal University, Saudi Arabia
Updates
Copyright
© 2021 Sankar, Ramachandran, Pandi, Mutharasappan, Ramasamy, Prabu, Shanmugaraj, Wang, Muniyandai, Rathinasamy, Chandrasekaran, Bayan, Jeyaraman, Halliah and Ebenezer.
This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Boomi Pandi, pboomi1983@gmail.comBalakumar Chandrasekaran, dhillbalu@gmail.com
†These authors have contributed equally to this work
This article was submitted to Molecular Diagnostics and Therapeutics, a section of the journal Frontiers in Molecular Biosciences
Disclaimer
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.