key: cord-1023523-eg8dih0o
authors: Tripathi, Manish Kumar; Singh, Pushpendra; Sharma, Sujata; Singh, Tej P.; Ethayathulla, A. S.; Kaur, Punit
title: Identification of bioactive molecule from Withania somnifera (Ashwagandha) as SARS-CoV-2 main protease inhibitor
date: 2020-07-08
journal: Journal of biomolecular structure & dynamics
DOI: 10.1080/07391102.2020.1790425
sha: 38f795b33ce33a4cdc53ea8522aeb30a948c1f70
doc_id: 1023523
cord_uid: eg8dih0o

SARS-CoV-2 is the causative agent of COVID-19 and has been declared as pandemic disease by World Health Organization. Lack of targeted therapeutics and vaccines for COVID-2019 have triggered the scientific community to develop new vaccines or drugs against this novel virus. Many synthetic compounds and antimalarial drugs are undergoing clinical trials. The traditional medical practitioners widely use Indian medicinal plant Withania somnifera (Ashwagandha) natural constituents, called withanolides for curing various diseases. The main protease (M(pro)) of SARS-CoV-2 plays a vital role in disease propagation by processing the polyproteins which are required for its replication. Hence, it denotes a significant target for drug discovery. In the present study, we evaluate the potential of 40 natural chemical constituents of Ashwagandha to explore a possible inhibitor against main protease of SARS-CoV-2 by adopting the computational approach. The docking study revealed that four constituents of Ashwagandha; Withanoside II (-11.30 Kcal/mol), Withanoside IV (-11.02 Kcal/mol), Withanoside V (-8.96 Kcal/mol) and Sitoindoside IX (-8.37 Kcal/mol) exhibited the highest docking energy among the selected natural constituents. Further, MD simulation study of 100 ns predicts Withanoside V possess strong binding affinity and hydrogen-bonding interactions with the protein active site and indicates its stability in the active site. The binding free energy score also correlates with the highest score of −87.01 ± 5.01 Kcal/mol as compared to other selected compounds. In conclusion, our study suggests that Withanoside V in Ashwagandha may be serve as a potential inhibitor against M(pro) of SARS-CoV-2 to combat COVID-19 and may have an antiviral effect on nCoV. Communicated by Ramaswamy H. Sarma

The novel strain of coronavirus was first identified at the end of December 2019 and emerged during an outbreak in Wuhan, China. Coronavirus disease (COVID-19) is a transmissible infectious disease initiated by severe acute respiratory syndrome coronavirus-2 (SARS-CoV-2) and has become a severe threat to human health . On 30 January 2020, World Health Organization WHO declared this outbreak as a Public Health Emergency of international concern and later on 11 February 2020 named it as Coronavirus . It subsequently declared COVID-19 as a global pandemic on 11 March 2020 (Green, 2020) . Till 28 May 2020, WHO reported worldwide 5,593,631 confirmed cases and 353,334 deaths due to SARS-CoV-2 (WHO, 2020).

The term coronavirus was coined from the Latin term corona that means "crown" like shape (Rabi et al., 2020) . The suffix 2019 novel coronavirus or "2019-nCoV" coronavirus was given by WHO on 12 January 2020. SARS-CoV-2 transmission is via human-to-human, mainly through direct contact and/or by contact with the infected surface with touching of the face. Additionally, it also spreads through respiratory droplets released by coughs or sneezes . Common symptoms of coronavirus (SARS-CoV-2) include respiratory symptoms, cough, fever, shortness of breath, and more severe symptoms of this infection can cause pneumonia, severe acute respiratory syndrome, and even death (Qiu et al., 2020) . SARS-CoV-2 coronavirus commonly belongs to the beta-CoVs category and has 60-140 nm diameter with a round or elliptical shape (Fehr & Perlman, 2015) .

SARS-CoV-2 is a positive-stranded RNA virus wherein part of its genome encodes, 16 non-structural proteins (Nsp1-16), viral RNA dependent RNA polymerase, RNA synthesis materials, and two large non-structural polyproteins. The remaining genome encodes four structural proteins, spike (S), envelope (E), membrane (M), nucleocapsid (N), and the other accessory proteins Lu et al., 2020) . A maximum number of non-structural proteins are known to play a significant role in the SARS-CoV-2 replication.

The main protease (M pro ) also known as 3-chymotrypsinlike protease (3CL pro ) is a key protein required for the proteolytic maturation of the virus (Morse et al., 2020; Stobart et al., 2012) . Thus, targeting main protease can provide an effective treatment against SARS-CoV-2 by inhibition of the viral polypeptide cleavage . It is one of the best-characterised drug targets among structural and non-structural proteins of coronavirus and currently Xray structures are reported as both liganded and un-liganded SARS-CoV-2 (Jin et al., 2020; Zhang et al., 2020) . Since this protein is absent in humans, it forms an excellent target for drug discovery. Moreover, there is presently no specific medicine or vaccine available for the treatment of COVID-19 infection caused by SARS-CoV-2.

Initial studies have indicated lopinavir, ritonavir, nelfinavir, antimalarial drug chloroquine, hydroxychloroquine and many other antiviral drugs, which have been screened as potential inhibitors against SARS-CoV-2 (Jin et al., 2020; Pant et al., 2020; Wang et al., 2020) . Bioactive natural compounds have diverse bioavailability, with relatively less non-toxic property. In the bioactive natural compounds, various phytoconstituents like phenols, steroids, flavonoid molecules are present, that have been reported and screened for their possible therapeutic effects against various viral diseases including COVID-19 (Bhardwaj et al., 2020; Das et al., 2020; Thuy et al., 2020) . Natural compounds such as crocin, digitoxigenin, b-eudesmol, withanone, caffeic acid phenethyl ester and other various bioactive molecules have also been reported to interact with SARS-CoV-2 main protease (Aanouz et al., 2020; Bhardwaj et al., 2020; Gyebi et al., 2020; Islam et al., 2020; .

Thus, in the present pandemic situation, there is an urgent need to investigate bioactive compounds which have the potency to combat this viral infection and provide immunity and strength to our body to combat SARS-CoV-2. Indian medicinal plants have been widely used for various diseases to develop drugs as well as afford strength to our immune system. There are countless benefits of herbal medicine such as naturally available, alleviates hormones and metabolism, anti-inflammatory, and the boosting up of the immune system. Perceiving the significance of immunity boosting actions during the COVID-19 situation, it is very vital to consume supplements in the form of immune nutrients such as vitamin C, B-complex, zinc and copper that will support our body to fight against the SARS-CoV-2. India's Ministry of AYUSH has also stated, "Ayurveda's immunity-boosting measures for self-care during COVID 19 crisis" (Ministry of Ayush, 2020). Withania somnifera (Ashwagandha) is one of the most valued Indian ayurvedic medicinal plants that are revered as a biologically active immunomodulator and its ability to balance, energise, rejuvenate, and revitalise (Vetvicka & Vetvickova, 2011) . Ashwagandha comprises a rich source of various phytoconstituents such as, Withaferin A, steroidal lactones of the withanolide series, steroidal alkaloids of the withanoside series and some other compounds (Glotter et al., 1973 (Glotter et al., , 1966 Hirayama et al., 1982; Matsuda et al., 2001; Tong et al., 2011) . Active constituents of Ashwagandha are known to have a promising anti-influenza activity, against chikungunya, inhibit infectious Bursal Disease virus and are beneficial for the treatment of genital disease caused by Herpes Simplex Virus among African tribes, etc. (Cai et al., 2015; Ganguly et al., 2018; Jain et al., 2018) .

In the current study, the aim is to explore the medicinal potential of Withania somnifera against the main protease of SARS-CoV-2 by adopting computational approaches. The data generated is very encouraging and suggests that Withania somnifera has the capability of being effective in the treatment of SARS-CoV-2.

The 40 natural chemical constituents of Withania somnifera (Ashwagandha) were obtained from reported literature and used as ligand dataset for this study ( Table 1 ). The 2D structures of these identified compounds are shown in supplementary table (Table S1 ). Further, these identified natural compounds were retrieved from PubChem database (https:// pubchem.ncbi.nlm.nih.gov/) as .sdf format and imported into the maestro module. Ligprep module was used for the preparation and optimization of these ligands (Schr€ odinger, 2018) and prepared ligands were used for further studies.

The glide module of Schr€ odinger Maestro was used to perform the molecular docking study to analyse the binding mode of identified phyto-compounds of Withania somnifera (Ashwagandha) into the active site of main protease SARS-CoV-2. The crystal structure of SARS-CoV-2 main protease was obtained from the protein data bank (PDB ID:6LU7) (Jin et al., 2020) . The crystal structure was further pre-processed, refined, optimized and minimized by the protein preparation wizard module. The missing side chain and loops in the crystal structure were included using the Prime module. Water molecules in the active site more than 5 Å from the hetero atoms including ligands were deleted. Further, the crystal structure was optimized by PROPKA method and minimized by using the OPLS3 force field (Jorgensen et al., 1996) . The receptor grid generation module was used for the generation of a grid around the centroid of co-crystallized ligand (N3). Re-docking of the co-crystallized ligand in the active site of main protease of SARS-CoV-2 was performed to validate the docking protocol. The docked complex was superimposed and aligned with the original crystal structure by employing the superimposing tool to calculate the RMSD value. The generated grid was used to perform the molecular docking with the prepared ligands of Ashwagandha obtained from the ligprep module. After docking, the binding affinity of the docked complex was predicted on the basis of the glide score and results analyzed using glide XP visualizer.

The classical molecular dynamics simulation was carried out for the docked complexes using the Desmond module of Schrodinger suite to affirm the stability and binding interaction of these docked complexes (Bowers et al., 2006) . The docked complexes were solvated using TI3P solvent model (specifies a three-site rigid water molecule with charges) in an orthorhombic box of 15 Â 15 Â 15 Å 3 size. For neutralizing the system, appropriate cation (Na þ ) or anion (Cl -) along with a salt concentration of 0.15 mol/L were added. The prepared system was energy minimized using steepest descent conjugate gradient method, where maximum interaction was set to 2000, and the convergence threshold was set to 1.0 kcal/mol/Å. Prior to a production run, the NPT simulation were performed at 300 K temperature set by Nose-Hoover thermostat and 1.01325 bar pressure to relax the complexes (Kaczor et al., 2015) . Finally, a production run of 100 ns was performed and trajectory was recorded at every 20 ps time. After the completion of simulation obtained trajectories were analysed by using various parameters of MD simulation, root mean square deviation (RMSD), protein root mean square fluctuation (RMSF) and protein ligand (PL) contacts with the help of simulation interaction analysis module of Desmond.

Molecular mechanics generalized Born surface area (MM-GBSA) calculations are widely performed in computational drug discovery process to estimate the binding free energy prediction of protein-ligand complexes. A total of 100 frames of protein ligand conformations spanned between 60 to 100 ns simulation run were extracted for binding energy calculation. The prime module of Schrodinger was used to perform MM-GBSA calculations (Hou et al., 2011) . The binding energy calculation (DG bind) occurs as per the below equation

where G complex, G receptor and G ligand represent the free energies of the complex, receptor and ligand respectively.

The reported literature was utilized to identify 40 phyto-compounds of Ashwagandha, and a database was created of these compounds by downloading their structures from PubChem database. All the identified compounds were docked in the active site of main protease of SARS-CoV-2. Molecular docking yielded four compounds which were selected on the basis of their docking score that was higher than the co-crystallized ligand N3 (-8.12 Kcal/mole). The docking score of the co-crystallized ligand N3 subsequently served as a control in this study. The docking scores of all Withasomnine 442877 (Schr€ oter et al., 1966) the remaining compounds are provided in Supplementary material Table S1 .

The main protease of SARS-COV-2 is a cysteine protease and is a homodimeric protein where each monomer is designated as a promoter. Each promoter comprises three domains each, domain I, domain II and domain III containing residues 8-101, residues 102-184, and residues 201-303 respectively. Domain II is connected to domain III through along loop from by residues 185 to 200 .

The residues 1 to 7 comprise the N-finger wherein the N-finger of one promoter (Promoter A) sits in the pocket formed in second promoter (Promoter B). The active site dyad amino acid residues, Cys-145 and His-41, are located between the gap formed between domain I and domain II. The binding site is surrounded by hydrophilic as well as hydrophobic amino acid residues with two negatively charged residues, Glu-166 and Asp-187, and one positively charged residue, Arg-188 (Jin et al., 2020; Yang et al., 2003) . The hydrophilic and hydrophobic amino acid surrounding the binding site comprises Thr-24, Leu-27, His-41, Met-49, Asn-142, His-163, His-164, Phe-140, Met-145, Gln-189, Pro-168 and Cys-145. The main protease protein is typically conserved amongst the various CoVs, and the absence of the protein in the humans makes it an attractive target for drug design. The initial docking protocol was validated by cross-docking the inhibitor N3 into the binding site of the protein, main protease of SARS-CoV-2 (PDB ID:6LU7). The docking protocol revealed a root mean square deviation (RMSD) of 1.0 Å between the positions of the ligand N3 observed in the crystal structure and the docked ligand. Subsequently, molecular docking study was performed within the region of main protease protein where the inhibitor bound with the validated grid. The obtained results after the docking study revealed that from all the selected compounds used in the present study, four compounds, Withanoside II (PubChem CID-101168811), Withanoside IV (PubChem CID-71312551) Withanoside V (PubChem CID-10700345) and Sitoindoside IX (PubChem CID-189586) possessed a higher docked score with a range of À8.37 to À11.30 kcal/mole as compared to the standard compound N3 whose docked score was found to À8.12 kcal/mole. The structures of these identified compounds are given in Table 2 .

The ligand Withanoside II interacts with the catalytic dyad amino acid residues, His-41 and Cys-145 with a docking score of À11.30 Kcal/mole (Table 3) . Withanoside II forms hydrogen bond interaction with His-41 and hydrophobic interaction with Cys-145 amino acid residues. Further, interaction analysis revealed that the compound formed hydrogen bond with Thr-26, Ser-46, Asn-142 and Gln-189 amino acid residues. Withanoside II also interacts through hydrophobic interactions with residues Met-49, Cys-145, Met-165, Pro-168 and charged interaction with His-164, Glu-166 amino acid residues ( Figure 1A,B) .

Withanoside IV docked into the protein active site and it formed interactions with the catalytic dyad residues His-41 and Cys-145 with a docking score of À11.02 Kcal/mole (Table 3) . It is also involved in hydrogen bonded interaction with the amino acid residues Thr-26, Leu-141, Gln-189, Thr-190 and Gln-192 as well as a number of hydrophobic interactions with Met-49, Pro-168, Leu-167, Met-165, Leu-27. Withanoside IV was further stabilized through polar interaction with amino acid residues His-41, His-163, His-164, Asn-142 and Glu-192, respectively (Figure 2A,B) .

Like Withanoside II and IV, Withanoside V also interacted with the active site dyad amino acid residues His-41 and Cys-145 present between gap of the two domains, domain I and domain II. The docking score of this docked complex was found to be À8.96 Kcal/mole (Table 3 ). The docked protein-ligand complex was stabilised through hydrogen bonding interactions with residues Thr-24, Thr-25, Thr-26, Gly-143, His-163, Glu-166. A number of hydrophobic interaction were observed between the ligand and protein comprising residue Met-49, Cys-44, Leu-141, Phe-140, Met-165, Leu-167 and Pro-168, respectively ( Figure 3A ,B).

The docking energy of Sitoindoside IX-protease docked complex was À8.37 Kcal/mole (Table 3) . Sitoindoside IX disclosed polar interaction with His-41 and hydrogen bonding interaction with Cys-145 amino acid of the catalytic dyad residues. Sitoindoside IX showed hydrogen bond interaction with Thr-24, Thr-26, Gly-143 and Cys-145 amino acid residues. The ligand formed hydrophobic interactions with amino acid residues in the substrate binding region ( Figure 4A,B) . From docking studies, we infer that the residues of the catalytic dyad Cys-145 and His-41 together with other amino acid residues have the highest binding capability for the selected compounds. This is comparable to the binding observed of the co-crystallized ligand N3 with the protein main protease protein ( Figure 5A,B) .

The molecular docking studies pinpoint that these are the most significant residues required for ligand binding situated in the substrate-binding region of the main protease enzyme. The compounds bind with the catalytic dyad residues and occupy the substrate binding region. A superposition of all compounds indicates that they adopt a similar mode of binding except for compound Withanoside IV ( Figure 6A,B) .

Further, we performed MD simulation study of 100 ns time scale of all these selected docked complex (Withanoside II, Withanoside IV, Withanoside V, Sitoindoside IX) and co-crystallized ligand N3 to check their stability in the active site of main protease SARS-COV-2 protein. The results obtained after the simulation were analysed from the generated equilibrium trajectories by using the various protocol of MD simulation, i.e. root mean square deviation (RMSD), root mean square fluctuations (RMSF) and interaction analysis of hydrogen bond and their stability during the simulation period.

Root mean square deviation measures the stability of the simulation system and calculates its conformational perturbations which occur in the protein backbone during the simulation time (Sargsyan et al., 2017) . After the simulation, obtained RMSD results were plotted against the simulation time for the analysis of results (Figure 7) .

From the RMSD graph, we predict that the protein backbone of main protease was stable during the entire simulation time with a mean value of 2 Å. These lower fluctuations indicate the stable behaviour of protein backbone during the entire simulation. Comparing the docked complex RMSD trajectories, we inferred that Withanoside II, Withanoside IV and Sitoindoside IX manifested fluctuations initially up to 30 ns time scale and after 30 ns the compound Withanoside II did not reveal any significant RMSD fluctuations as compared to the docked complex trajectories of the other compounds. The Withanoside V docked complex trajectory was comparatively more stable with a mean value of 3 Å as compared to the trajectories of other docked complexes. Moreover, the trajectory of the Withanoside V docked complex did not indicate any fluctuating behaviour; this indicates its stable complex formation and least conformational changes during simulation. Thus, from the RMSD trajectory analysis of all docked complexes, we observed that Withanoside V trajectory is steady throughout the simulation and changes in backbone RMSD are within the acceptable range.

Further, to understand the residue wise fluctuations in the docked complexes, root mean square fluctuations (RMSF) analysis was performed. The RMSF behaviour of all the docked complexes was generated from the stable trajectory behaviour and plotted against the amino acid residue in x-axis vs RMSF(Å) in the y-axis. The RMSF plot for Withanoside V, yielded the least fluctuations at the catalytic active site dyad, His-41 and Cys-145, as well as the substrate-binding pocket residues Glu-166, Asp-187 and Arg-188 with flexibility range of 0.6 to 3.0 Å (Figure 8) .

The RMSF pattern of Withanoside V is found to be comparable to the co-crystalized ligand N3. The other selected compounds (Withanoside II, Withanoside IV and Sitoindoside IX) demonstrated higher fluctuations behaviour (Supplementary material Figure S1 ).

The protein-ligand interaction analysis was also assessed for all the docked complexes formed during the simulation time to check the stability of docked complexes. Thus, the hydrogen-bonding pattern of all the identified hits was evaluated and plotted against the simulation time scale along with the percentage interaction of active site residues. The interaction pattern of Withanoside II obtained from MD simulation, divulged hydrogen bonding interaction with His-41 amino acid residue apart from Cys-145 in the catalytic dyad region. The ligand-protein contact analysis revealed that catalytic dyad residue Cys-145 formed a hydrogen bond through a solvated molecule with the main protease protein residues. The histogram representation of interaction fraction and timeline of percentage index is represented in Supplementary material Figure S2 . Further, analysis of result at the substrate binding region revealed, interactions with residues Asn-142, Glu-166 and Gln-189. The RMSD trajectory analysis of Withanoside II was also not found to be stable. The interaction pattern of the selected compound, Withanoside IV obtained from MD simulation indicated hydrogen bonding interactions with the catalytic dyad amino acid residues His-41 and Cys-145 wherein the interaction with Cys-145 residue was maintained during the entire simulation time. The RMSD pattern of this docked complex showed greater fluctuations as compared to all the other docked complexes. Withanoside IV like Withanoside II possessed interactions with the specific substrate-binding region amino acid residues Asn-142, Glu-166 and Gln-189 (Supplementary material Figure S3 ).

The MD simulation interaction analysis of Withanoside V with main protease protein was also evaluated and plotted as a percentage index of total contacts during the simulation time. The hydrogen bond interactions with catalytic dyad residues His-41 and Cys-145 were seen to be retained during the entire simulation with more than 53% and 40% simulation time. The polar interaction with His-41 and hydrophobic interaction with Cys-145 amino acid residues formed in the docked complex were also maintained throughout the simulation. Moreover, a stable RMSD behaviour of its trajectory was observed during simulation indicating the stability of this protein-ligand complex in the active site. The histogram, in combination with the 2D simulation interaction diagram between Withanoside V with main protease protein is illustrated in Figure 9 . The interaction analysis revealed that in addition to the stable interactions present at the catalytic dyad interaction, substrate binding site residues Thr-24, were also observed to be stable during the entire simulation time.

The MD simulation of Sitoindoside IX-protein complex revealed that the interaction with His-41 in the catalytic dyad region was not adequately retained during the simulation. The histogram in combination with the 2D simulation interaction diagram indicated that the interactions present at the substrate binding region with residues Asn-142 and Glu-166 were not stable during the entire simulation scale (Supplementary material Figure S4 ).

The interaction pattern for the co-crystallized ligand N3, obtained from MD simulation showcased hydrogen bonding interactions with catalytic dyad amino acid residues His-41 and Cys-145. However, its interaction with Cys-145 residue was not observed during the entire simulation time. The histogram representation of interaction fraction illustrated interactions with the residues in the specific substrate-binding region, Asn-142, Glu-166 and Gln-189 (Supplementary material Figure S5 ).

The binding free energy calculation was carried out for all the systems (Withanoside II, Withanoside IV, Withanoside V, Sitoindoside IX) along with the co-crystallized ligand N3. The prime MM/GBSA (molecular mechanics energies combined with the generalized Born and surface area continuum solvation) free energy calculation is widely accepted and commonly used to estimate the ligand-binding affinities in the protein systems. To evaluate the binding free energy, 100 conformations that spanned between the 60 to 100 ns time interval of the MD production run were extracted.

The MM-GBSA DG bind score obtained after the calculations is provided in Table 4 . The calculated binding free energy of these compounds Withanoside II (-62.50 ± 5.25 Kcal/mole), Withanoside IV (-81.29 ± 4.78 Kcal/mole), Withanoside V (-87.01 ± 5.01 Kcal/mole), Sitoindoside IX (-49.90 ± 4.15 Kcal/mole) and for the co-crystallized ligand (-86.60 ± 6.16 Kcal/mole) were then determined. The negative values of binding energy indicate that all the hit compounds interact favourably with the receptor. These results clearly suggest that Withanoside V possessed the maximum negative binding energy (-87.01 ± 5.01 Kcal/mol) as compared to the other selected compounds. The binding energy of Withanoside II, and Sitoindoside IX was not found to be comparable with the co-crystallized ligand N3. Thus, Withanoside V with the least binding energy and better binding affinity against the target could consequently prove to be a potential inhibitor of main protease of SARS-CoV-2. The comparison of the interaction patterns of all the selected hit compounds pinpoints Withanoside V, as the compound of choice as it revealed stable interactions with both the residues of catalytic dyad (His-41 and Cys-145) in main protease of SARS-CoV-2. Though Withanoside IV exhibited stable interactions with the catalytic dyad residues, however, its RMSD behaviour was not within the acceptable range as it presented higher fluctuation and its binding energy was also found to be À81.29 ± 4.78 Kcal/mole. In contrast, Withanoside V displayed stable interactions with the substrate binding site amino acid residues, Thr-24, Asn-142, His-163, Glu-166, Gln-189 throughout the MD simulation. The co-crystallized ligand N3 did not retain the interaction with Cys-145 throughout the simulation but its interaction with substrate binding regions was present throughout. The residue Glu-166 in main protease is known to be involved in the formation of its functional dimeric form . The ability of Withanoside V to bind to the catalytic dyad along with the substrate binding region which houses the Glu-166 residues prevents the dimerization of main protease as the residue Glu-166 of one monomer/promoter is unable to interact with the N-finger residues of other monomer. The MM-GBSA binding energy of Withanoside V was also found to be most negative, À87.01 ± 5.01 Kcal/mole, compared to other selected compounds. Thus, based on these observations, we conclude that Withanoside V present in Withania somnifera (Ashwagandha) can be a potent and efficacious biomolecule which interacts with SARS-CoV-2 and can be further explored as a drug candidate.

The molecular docking and dynamics simulation approach was applied for the natural chemical constituents of Withania somnifera (Ashwagandha) to identify their potential against the main protease protein of SARS-CoV-2, a pandemic COVID-19 that is presently wreaking havoc across the world. Our study suggested that out of the forty constituents of Ashwagandha, Withanoside V was the most potent natural inhibitor. Moreover, the RMSD trajectories of Withanoside V exhibited a stable behaviour with comparatively lesser fluctuations as compared to other ligands. The ligand displayed stable interactions with the catalytic dyad residues (His-41 and Cys-145) of main protease as well the substrate binding region which promotes the dimerization of the protein.

Further, the MM-GBSA score of Withanoside V presented the maximum negative binding energy (-87.01 ± 5.01 Kcal/mol). M pro promotes activation for the polyactive site of main protease, which is important for the replication of the SARS-CoV-2. Hence, our study concludes that Withanoside V which is present in the plant Withania somnifera is a potent inhibitor of main protease of SARS-CoV-2. The status of Withania somnifera as an immunity enhancer, anti-diabetic agent and blood pressure and hormonal regulation is well known. The present study suggests that additionally, the Ayurveda herb Ashwagandha could prove to be an alternative to available treatments of COVID-19. Moreover, the backbone structure of Withanoside V could also be further exploited to develop more potent inhibitors of SARS-CoV-2 M pro . However, further in vitro and in vivo studies are necessary for validation of this compound to combat COVID-19.

Moroccan medicinal plants as inhibitors against SARS-CoV-2 main protease: Computational investigations

New withanolides and other constituents from the fruit of Withania somnifera

The withanolides of Withania somnifera chemotypes I and II

Comparative study of withanolide production and the related transcriptional responses of biosynthetic genes in fungi elicited cell suspension culture of Withania somnifera in shake flask and bioreactor

8 -Withania somnifera: Ethnobotany, pharmacology, and therapeutic functions

Coronavirus main proteinase (3CLPro) structure: Basis for design of anti-SARS drugs

Bioactive metabolites from the leaves of Withania adpressa

Withanolide C, A chlorinated withanolide from Withania somnifera

Identification of bioactive molecules from Tea plant as SARS-CoV-2 main protease inhibitors

Scalable algorithms for molecular dynamics simulations on commodity

Promising anti-influenza properties of active constituent of Withania somnifera ayurvedic herb in targeting neuraminidase of H1N1 influenza: Computational study

Fortschritte der Chemie organischer Naturstoffe Progress in the Chemistry of Organic Natural Products

Analysis of withanolides in root and leaf of Withania somnifera by HPLC with photodiode array and evaporative light scattering detection

Withanone, an active constituent from Withania somnifera, affords protection against NMDA-induced excitotoxicity in neuron-like cells

An investigation into the identification of potential inhibitors of SARS-CoV-2 main protease using molecular docking study

Somniferine, a novel dimeric opium alkaloid

Coronaviruses: An overview of their replication and pathogenesis

Nitric oxide induced by Indian ginseng root extract inhibits Infectious Bursal Disease virus in chicken embryo fibroblasts in vitro

Constituents of Withania somnifera Dun-XIII: The withanolides of chemotype III

Constituents of Withania somnifera Dun. Part VIII. A new steroidal lactone, 27-deoxy-14-hydroxywithaferin A

Did the hesitancy in declaring COVID-19 a pandemic reflect a need to redefine the term? The Lancet

Potential inhibitors of coronavirus 3-chymotrypsin-like protease (3CLpro): An in silico screening of alkaloids and terpenoids from African medicinal plants

Withanolide E sensitizes renal carcinoma cells to TRAIL-induced apoptosis by increasing cFLIP degradation

Stereoselective synthesis of withafein A and 27-deoxywithaferin A1

Assessing the performance of the MM/PBSA and MM/GBSA methods. 1. The accuracy of binding free energy calculations based on molecular dynamics simulations

A molecular modeling approach to identify effective antiviral phytochemicals against the main protease of SARS-CoV-2

In vivo evaluation of Withania somnifera-based Indian traditional formulation (Amukkara Choornam), against Chikungunya virus-induced morbidity and arthralgia

Growth inhibition of human tumor cell lines by withanolides from Withania somnifera leaves

Structure of Mpro from SARS-CoV-2 and discovery of its inhibitors

Development and testing of the OPLS all-atom force field on conformational energetics and properties of organic liquids

Comparative molecular field analysis and molecular dynamics studies of a/b hydrolase domain containing 6 (ABHD6) inhibitors

Withanolides Q and R, two new 23-hydroxy-steroidal lactones

Restoration of stress-induced altered T cell function and corresponding cytokines patterns by Withanolide A

Withanoside IV and its active metabolite, sominone, attenuate Abeta(25-35)-induced neurodegeneration

Identification of phytochemical inhibitors against main protease of COVID-19 using molecular modeling approaches

Withanone and caffeic acid phenethyl ester are predicted to interact with main protease (Mpro) of SARS-CoV-2 and inhibit its activity

Barrier protective effects of withaferin A in HMGB1-induced inflammatory responses in both cellular and animal models

Research and development on therapeutic agents and vaccines for COVID-19 and related human coronavirus diseases

Withaferin A-related steroids from Withania aristata exhibit potent antiproliferative activity by inducing apoptosis in human tumor cells

Genomic characterisation and epidemiology of 2019 novel coronavirus: Implications for virus origins and receptor binding

Structures of withanosides

Ayurveda's immunity boosting measures for self care during COVID 19 crisis

Bak compensated for Bax in p53-null cells to release cytochrome c for the initiation of mitochondrial signaling during Withanolide D-induced apoptosis

Learning from the past: Possible urgent prevention and treatment options for severe acute respiratory infections caused by 2019-nCoV

Screening of promising chemotherapeutic candidates from plants against human adult T-cell leukemia/lymphoma (II): Apoptosis of antiproliferactive principle (24,25-dihydrowithanolide D) against ATL cell lines and structure-activity relationships with withanolides isolated from solanaceous plants

Peptide-like and small-molecule inhibitors against Covid-19

Transmission routes of 2019-nCoV and controls in dental practice

Clinical and epidemiological features of 36 children with coronavirus disease 2019 (COVID-19) in Zhejiang, China: An observational cohort study

SARS-CoV-2 and coronavirus disease 2019: What we know so far

Molecular docking and dynamic simulation studies evidenced plausible immunotherapeutic anticancer property by Withaferin A targeting indoleamine 2,3-dioxygenase

An in silico guided identification of nAChR agonists from Withania somnifera

How molecular size impacts RMSD applications in molecular dynamics simulations

Withasomnine. A pyrazole alkaloid from Withania somnifera Dun

Secondary metabolites of medicinal plants: Ethnopharmacological properties, biological activity and production strategies

Silencing of sterol glycosyltransferases modulates the withanolide biosynthesis and leads to compromised basal immunity of Withania somnifera

Impaired motor learning attributed to altered AMPA receptor function in the cerebellum of rats with temporal lobe epilepsy: Ameliorating effects of Withania somnifera and withanolide A

Temperature-sensitive mutants and revertants in the coronavirus nonstructural protein 5 protease (3CLpro) define residues involved in long-distance communication and regulation of protease activity

Ashwagandhanolide, a bioactive dimeric thiowithanolide isolated from the roots of Withania somnifera

Investigation into SARS-CoV-2 resistance of compounds in garlic essential oil

Chlorinated Withanolides from Withania somnifera

Withania somnifera induces cytotoxic and cytostatic effects on human T leukemia cells

Immune enhancing effects of WB365, a novel combination of Ashwagandha (Withania somnifera) and Maitake (Grifola frondosa) extracts

Withanone as an inhibitor of survivin: A potential drug candidate for cancer therapy

Remdesivir and chloroquine effectively inhibit the recently emerged novel coronavirus (2019-nCoV) in vitro

Natural withanolides in the treatment of chronic diseases

Coronavirus disease (COVID-19) Situation Report -129

Analysis of therapeutic targets for SARS-CoV-2 and discovery of potential drugs by computational methods

The crystal structures of severe acute respiratory syndrome virus main protease and its complex with an inhibitor

Crystal structure of SARS-CoV-2 main protease provides a basis for design of improved a-ketoamide inhibitors

Withanolide Derivatives from the Roots of Withania somnifera and Their Neurite Outgrowth Activities

A pneumonia outbreak associated with a new coronavirus of probable bat origin

The authors have no conflict of interest regarding the publication of this article.

http://orcid.org/0000-0002-7358-3716