key: cord-0693846-bickmuv5
authors: Chikhale, Rupesh V.; Sinha, Saurabh K.; Patil, Rajesh B.; Prasad, Satyendra K.; Shakya, Anshul; Gurav, Nilambari; Prasad, Rupali; Dhaswadikar, Suhas R.; Wanjari, Manish; Gurav, Shailendra S.
title: In-silico investigation of phytochemicals from Asparagus racemosus as plausible antiviral agent in COVID-19
date: 2020-06-24
journal: J Biomol Struct Dyn
DOI: 10.1080/07391102.2020.1784289
sha: 413df35689085a42cb27fad1bb7faf3f80957487
doc_id: 693846
cord_uid: bickmuv5

COVID-19 has ravaged the world and is the greatest of pandemics in human history, in the absence of treatment or vaccine the mortality and morbidity rates are very high. The present investigation was undertaken to screen and identify the potent leads from the Indian Ayurvedic herb, Asparagus racemosus HIGHLIGHTS: Asparagus racemosus have antiviral potential. Phytochemicals of Shatavari showed promising in-silico docking and MD results. Asparaoside-C and Asparoside-F has good binding with target proteins. Asparagus racemosus holds promise as SARS-COV-2 (S) and (N) proteins inhibitor . Communicated by Ramaswamy H. Sarma

During the month of December 2019, a global pandemic emerged from the Wuhan city of China, called as the severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) and has become a major threat to the whole world due to the range of population it has infected throughout the world (Zhou et al., 2020) . The pandemic COVID 19 has affected around 212 countries with almost half million death throughout the globe. The corona viruses are highly enveloped single stranded RNA with a largest genome size (26-32-kilobases) that plays critical role in initial RNA synthesis of the infectious cycle, template for replication and transcription and also act as a substrate for packaging into the progeny virus. In all types of corona viruses, two-thirds of the genome encodes a replicase polyproteins, pp1ab, that consists of two overlapping open reading frames (ORFs) i.e. ORF1a and ORF1b and are processed by viral proteases to undergo cleavage into 16 different nonstructural proteins (NSPs), involved in transcription and replication (Boopathi et al., 2020; Cotten et al., 2013; Gupta et al., 2020) . Among these the NSP15 Endoribonuclease is a nidoviral RNA uridylate-specific endoribonuclease (NendoU) which comprises of C-terminal catalytic domain, specific for uridine, acts on single-as well as double stranded RNA (Elfiky, 2020; Mittal et al., 2020) . Studies have reported that the NendoU activity of NSP15 is mainly attributed to the observed protein interference with innate immune response and is therefore, a very essential segment for understanding the biological progression of coronavirus (Bhardwaj et al., 2008; Deng & Baker, 2018) . Another NSP that helps 2019-nCoV to get entry into the host cells involves a densely glycosylated, homotrimeric class I fusion spike (S) protein (Sarma et al., 2020) . This S protein present in a metastable prefusion conformation, undergoes structural rearrangement for the viral membrane to get fused with the membrane of host cell. The process is accelerated as a result of binding of S1 subunit to the host-cell receptor and transition of the S2 subunit to a highly stable postfusion conformation (Hasan et al., 2020) . Further, recent studies have reported that, spike (S) glycoprotein possess a human angiotensin-converting enzyme 2 (ACE2) binding site and has 10 to 20-fold higher binding affinity towards 2019-nCoV S as compared to SARS-CoVs (Chikhale et al., 2020; Choudhury, 2020; Han & Kr al, 2020; Ou et al., 2020; Pillaiyar et al., 2020; Sinha et al., 2020a) .

Even though, scientist throughout the world are trying their level best to identify or discover a new lead or repurposing of drugs which could be effective against SARS-CoV-2, it is a long way to go before we have a vaccine or drug molecule to treat COVID-19.The U.S. Food and Drug Administration has now issued an emergency use authorization for Remdesivir for the treatment of suspected or laboratory-confirmed adults and children hospitalized with severe symptoms of SARS-CoV-2 disease (FDA news release, 2020). The effectiveness of the treatment based on traditional medicinal plants has been reported during the 2003 SARS (Wen et al. 2011; Hoever et al,. 2005; Chen et al., 2008; Li et al., 2005; Lin et al., 2008; Ryu et al., 2010) . Therefore, studies on medicinal plants based on their history and traditional uses have already started among the scientific community (Al-Khafaji et al., 2020; Elmezayen et al., 2020; Enmozhi et al., 2020; Joshi et al., 2020; Sinha et al., 2020b; Wahedi et al., 2020; Zhang et al., 2020) .

Asparagus racemosus (Willd.) is a well-known medicinal plant, grown in the tropical and subtropical regions of India, where its therapeutic importance is well documented in the Indian and British Pharmacopoeia along with various traditional system of medicines such as Ayurveda, Siddha and Unani (Gurav & Gurav, 2014; Singh & Geetanjali, 2016) . Traditionally, this plant is popularly known as 'Shatavari', and described in Ayurveda as a potential rasayana which prevents ageing, provides immunity, increases longevity, improves mental function and also helps in treatment of diseases related to female reproductive system, inflammation, dysentery, biliousness, tumor and diseases of the blood and eyes (Alok et al., 2013; Singh & Geetanjali, 2016) . The pharmacological studies carried out on the plant has revealed the potential adaptogenic, immunostimulant, anti-inflammatory, anti-microbial, anti-oxidant, anti-ulcer, galactogogue, phytoestrogenic, neuroprotective, aphrodisiac, anti-dyspepsia, anti-tussive and anti-cancerous activities (Alok et al., 2013; Boonsom et al., 2012; Singh & Geetanjali, 2016; Upadhyay et al., 2020) . The major phytoconstituents reported in the plant includes steroidal saponins such as Shatavarin (I-XI mainly I, IV and X), Aasparacosides (A-F), Racemosol, Asparagamine, Racemoside A-B, Asparagoside F, Quercetin, Rutin, Kaempferol and many more (Alok et al., 2013; Lalert et al., 2018; Singh & Geetanjali, 2016) .

Shatavari has been reported for its anti-bacterial, anti-candidal, anti-viral activities (Mandal et al., 2000; Raza et al., 2015; Shah et al., 2014; Uma et al., 2009) . Recent studies have reported that, the drug under investigation, if possess a potent immunomodulatory, anti-inflammatory and antioxidant potential could prove an effective agent in treatment of COVID 19 (Jayawardena et al., 2020; Patwardhan et al., 2020; Tillu et al., 2020) . Further, all these suggested properties are well equipped by the phytoconstituents from Asparagus racemosus.

Thus, keeping the above view into consideration, the present investigation was undertaken to determine the efficacy of 32 different bioactive molecules from the plant Asparagus racemosus against two different protein targets i.e. NSP15 Endoribonuclease and spike receptor-binding domain from SARS-CoV-2 with the help of molecular docking and molecular dynamics studies. 

To determine the potency and understand the insights into the possible mechanism of 32 different herbal-based ligands, molecular docking simulation study was carried out on two proteins i.e. spike receptor-binding domain and NSP15 Endoribonuclease from SARS-CoV-2.The spike receptor-binding domain from the virus gets entry into the host cell by combining with ACE-2 sreceptor, while NSP15 Endoribonuclease is said to be responsible for replication of virus genetic material. Docking simulation study was also performed on the remdesivir, a reference drug, on the same protein targets and the results were compared with the test ligands (Muralidharan et al., 2020) . Remdesivir was chosen as a reference drug as it is a well-known RNA dependent RNA polymerase, a nsp protein, inhibitor. Further, remdesivir was chosen as reference drug because it is being used as a agent under clinical investigation for the treatment of ncov-19 (Sinha et al., 2020b) . On the basis of the docking score (Table 1 ) and binding interaction (Tables 2  and 3) , only the best five tested ligands along with the reference drug with respect to different proteins have been discussed in our study. The docking scores of top scorer A. racemosus bio-actives are given in Table 1 whereas the docking scores of other phytochemicals are given in supplementary information Table S1 . The structures of all A. racemosus phytochemicals investigated in present investigation are given in supplementary information Figure S1 . Figure 1A ). The result outcomes from the docking of remdesivir with SARS-CoV-2 NSP15 Endoribonuclease are presented in (Table  1 and Figure 1B) . Remdesivir was connected with the binding site mainly through H-bonding and showed docking score -5.94 Kcal/mol. The hydroxyl group of tetrahydrofuran ring exhibited H-bonding with Val 292 and the oxygen of phosphoramide exhibited bi-furcated H-bonding with (protonated His) Hip 235, Thr 341 ( Figure 3A ).

SARS-CoV-2 gets entry into the host cell through spike receptor-binding domain, which gets fused with ACE-2 receptor and its binding interaction depends on the binding of S1 subunit, having RBD. The RBD is consists of 5 twisted b sheets i.e. b 1 , b 2 , b 3 , b 4 and b 7 which are antiparallel to each other, where the extended inclusion is observed between b 4 and b 7 containing some a loops which designated as receptor binding motif (RBM). This RBM are said to contain almost all the residues that plays critical role in connectingn-COVID19 with ACE-2. Recently, it has been confirmed that, from all the residues between Arg 319 to Phe 541 only the 17 residues Lys 417, Gly (Figure 2A ). The bonding energy of Asparoside-C was found to be higher than all the reference drugs. All the oxane rings of was found to be properly buried into the outer concave surface of RBM and stabilized by 7H-binding interaction.

In accordance with our previous study (Sinha et al., 2020a) , in which we explained the importance of oxane rings of saikosaponins during binding. Figure 2C ). After structural elucidation of Asparoside-C and Asparoside-D we found that number of oxane ring was same in both the ligand, but aglycone part was differ. A hydroxyl group of aglycone part in Asparoside-C is naturally derivetized into methoxy in Asparoside-D, so we can say that Asparoside-D is the methoxy derivative of Asparoside-C and may be responsible for the lower binding affinity of Asparoside-D. Ser 494 residue was found to be common to form H-bond with both ligand and also found in the close proximity with other ligands in binding pocket. Figure 3 ).

The NSP15, bears a catalytic C-terminal domain, which is constitutes of a set of trimers i.e. hexamers and responsible for cutting the double-stranded (ds) RNA substrates with uridine specificity via the Mn þ 2-dependent endoribonuclease activity. Additionally, each monomeric unit is comprised of $345 amino acids which are folded into three domains viz. Racemoside-A (-5.993 Kcal/mol) was selected as a potent lead. All the selected ligand was found to be tightly fit into the binding pocket and interaction with the catalytic triad. These ligands also showed higher docking score as compare to the reference drugs. Hydroxyl group of oxane ring of Asparoside-C stabilized through H-bonding with Glu 234, Gly 230, Val 292, Hip 235, Asp 240 and hydroxymethyl group showed H-bonding with Val 292 ( Figure 4A ). Asparoside-F was found to interact with the binding pocket through trifurcated H-bonding with Glu 234, Gly 230, Ala 232, in which two are donor and one acceptor in nature, bi-furcated Hbonding with Hip 235, Asp 240 and conventional H-bonding with Val 339, Glu 234 ( Figure 4B ).

Rutin was found to be the least ranked ligand and get stabilized into the binding pocket through the chromenone and oxane ring. 

with SARS-CoV-2 main protease, papain like protease and RNA dependent RNA polymerase SARS-CoV2 main protease (M pro ), papain like protease (PL pro ) are other viral proteins which facilitates viral assembly by cleaving polyproteins (Mengist et al., 2020; Cardenas-Conejo et al., 2020) . RNA dependent RNA polymerase (RdRp) is another nonstructural protein which is vital for viral life cycle (Elfiky, 2020) . The docking study against these protein with the most promising bioactive ligand (Asparoside-C) against spike RBD and NSP 15 proteins was also performed. The crystal structure of SARS-Cov2 main protease (PDB ID: 4MDS) with 1.60 Å resolution was used in docking studies. All the previously mentioned protein refinement steps were performed. The aglycone part of Asparoside-C was properly positioned in to the binding site of main protease. The hydroxyl groups of sugar moieties form hydrogen bonds with the residues His41, Thr 190, Ala191, Glu 166, Asn 142, Asn 119, Thr 26 and Thr 24 ( Figure 5A ). Compared to the reference drug Remdesivir, which forms the hydrogen bond with Glu 166, Asparoside-C binds more favorably at the binding site. The docking score for Asparoside-C and Remdesivir are -7.431 and -6.099 respectively, which also suggests more binding affinity of Asaparoside-C. In case of SARS-CoV2 papain like protease (PL pro ), the experimentally solved crystal structure (PDB ID: 6WX4) with resolution 1.66 Å was used in docking study. The Asparoside-C binds at the binding site where the hydroxyl groups of sugar moieties form hydrogen bond with positively charged Lys 157, Arg 166, polar Thr 301, negatively charged Glu 161, Asp164, Glu 167 and hydrophobic Pro 248 ( Figure 5B ). On the other hand, the phenyl ring in Remdesivir forms p-cation interaction with positively charged Arg 166, while the pyrrole ring forms p-p stacking interaction with Thr 264 and only one hydrogen bond is formed with Thr 301. The docking scores for Asparoside-C and Remdesivir were -5.448 and -4.055 respectively.

RdRp crystal structure with resolution 2.90 Å (PDB ID: 6M71) was used to investigate possible mode of binding of Asparoside-C and compared with the binding mode of RdRp inhibitor, Remdesivir. The aglycone part of Asparoside-C occupies binding site in such a way that the major glycosidic part of this molecule accommodates in a binding cavity surrounded by hydrophobic Tyr 455, Tyr 458 and Val 166 residues. The hydroxyl groups of sugar moieties form hydrogen bond with Lys 798, Glu 167, Tyr 455, Arg 457, Asn 691, Asp 623, Asp 452, while the oxygen atom from glycosidic linkage forms hydrogen bond with positively charged Lys 621 residue ( Figure 5C ). Remdesivir, on the other hand forms hydrogen bond interactions with Asp 164, Lys 621 and Cys 622 residues. The charged Lys 621 also forms the p-cation interaction with pyrrole ring in Remdesivir. These docking scores for Asparoside-C and Remdesivir are -6.652 and -4.661 respectively. The docking study on these proteins suggest that Asparoside-C has favorable binding interactions and reasobaly good docking scores compared to reference drug Remdesivir.

Molecular dynamics simulations (MDS) provides with greater insight into binding of ligand molecules with the target protein under investigation. A workflow integrated with molecular docking, molecular dynamics simulations and free energy calculations has been employed by us in understanding the properties of various small molecule inhibitors (Chikhale et al., 2018; Kerzare et al., 2016; Pant et al., 2014) . The MDS studies of protein-ligand complex over a sufficiently long period allow exploring a wide space of conformational optimization and binding affinity. In this study a 100 ns MDS was performed to evaluate the conformational stability and binding affinity of the ligands bound to the viral spike receptor-binding domain and NSP15 endoribonuclease proteins. Six phytochemicals, Asparoside-C, Asparoside-D, Shatavarin-I in complex with SARS-CoV-2 spike protein receptor-binding domain and Asparoside-C, Asparoside-D, Asparoside-F in complex with NSP15 endoribonuclease protein were selected for the MDS based on their docking performance. Along with these phytochemicals, standard drug Remdesivir was studied for its interactions with spike protein receptor-binding domain and NSP15 endoribonuclease protein. MDS was performed with stages of three minimization steps, one heating step and two equilibration steps followed by 100 ns of production stage. MDS of these complexes provides for analyzing binding modes, formation of hydrogen bond, pi-pi interactions and van der Waals interactions. The MDS trajectories are analyzed study the following parameters; RMSD, RMSF, Ligand-RMSD and the binding free energies on the complexes with MMGBSA calculations.

SARS-CoV-2 spike receptor-binding domain binds with the host ACE-2 receptor leading to the viral entry in the host cell. Targeting this protein is a clear and direct approach to inhibit the COVID-19 infection. Three phytochemicals Asparoside-C, Asparoside-D, Shatavarin-I bound to spike protein were obtained from the docking studies. Post-MDS RMSD analysis showed Asparoside-C-spike protein complex achieves the conformational stability 12.5 ns with a RMSD of 2.5 Å to the initial structure and stays in this conformation for rest of the period of time (Figure 6, Black) . The Asparoside-D-spike protein complex shows optimum conformational stability around 2.5 Å between 15 to 60 ns of the MDS after which the RMSD increases to about 4 Å (Figure 6 , Red). The Shatavarin-I-spike protein complex is less stable conformationally with varying RMSD between 2-4 Å ( Figure 6 , Green), however, this is not a very large fluctuation. In case of Remdesivir, it had a stable RMSD convergence for first 25 ns around 0.5 Å. However, the RMSD increased gradually throughout the simulation and fluctuated to 6 Å towards the end, this could be attributed to the breaking of some interactions between the ligand and receptor active site residues.

The RMSF for individual amino acid (aa) reflects the trend observed in RMSD of complexes (Figure 7) . The RMSF per residue of Shatavarin-I-spike protein complex is higher compared to other tow complex. The RMSF of residue number between aa25 to aa50 is higher followed by residues between aa140 to aa160. To understand the conformational stability of the ligand bound in the active site ligand RMSD was studied for individual ligands; Asparoside-C, Asparoside-D, Shatavarin-I and Remdesivir.

The ligand RMSD was calculated to understand the ligand conformations during the period of MDS. The Asparoside-C (Figure 8 , Black) equilibrates for first 15 ns below 1 Å and then fluctuates by about 1.5 Å to stabilize for next 20 ns and then again fluctuates between 35 to 60 ns for about 1.5 Å. It equilibrates and stabilizes to 2 Å from the mean position after 60 ns. The Asparoside-D RMSD rises to around 2 Å and equilibrates for about 60 ns and then again increases to 4 Å after 60 ns (Figure 8, Red) . The RMSD of Shatavarin-I increases consistently over the MSD time form 0.5 Å to 4.5 and then equilibrates around 4.5 Å towards the end of MDS.

These ligands were studied further for their interactions with the individual residues port MDS. Asparoside-C bound to spike protein has two hydrogen bonds between hydroxyl groups of the Oxane rings and Gln493 and Ser494 ( Figure  9a ) before the start of production stage of the MDS. The Asparoside-C (O7-Ser494-O, Bond length 2.79 Å) hydrogen bond breaks during the MDS but Asparoside-C (O7-Gln493-NE2, Bond length 6.63 Å) is formed and stabilizes the complex (Figure 9b) .

The receptor surface analysis of the Asparoside-C bound to spike protein shows that the ligands are well accommodated in the active site during the MDS (Figure 9c) , the superposition of ligand conformation at the beginning and end of MDS shed light on the stability of the complex (Figure 9d) , this is also supported by the equilibrated ligand RMSD after 15 ns of the MDS. Three other phytochemicals Asparoside-C, Asparoside-D and Asparoside-F bound to SARS-CoV-2 NSP15 Endoribonuclease were selected from the docked complexes for MDS.

The binding mode of Remdesivir with the RBD domain of S protein is elaborated in Figure 10 . It was observed from the protein and ligand RMSD that there is gradual rise in RMSDs, on visual inspection of the MDS trajectory and analyzing the bond breaking and formation it was found that there are several conformational changes in Remdesivir during the MDS. These conformational changes contribute to the fluctuations in the RMSD. The hydrogen bond between Remdesivir and RBD domain (N12-Arg71) reduces in bond length from 3.26 to 2.90 Å, which explains the drop in RMSD in the last 5 ns (Figures 6 and 8) .

RMSD analysis of these complexes of phytochemicals and one of Remdesivir was performed the results are presented in Figure 11 . Asparoside-C and D-NSP15 endoribonuclease complexes (Figure 11 , black and Red respectively) shows similar pattern of RMSD with the initial equilibration to 2 Å for 35 ns. After 40 ns both the complexes shows increase in RMSD by another 1 Å and stays equilibrated in the same conformation till the end of MDS. NSP15 bound to Remdesivir also show similar trend and converges between 30 and 100 ns to RMSD between 2.5 to 3.0 Å (Figure 11 , Blue) Contrary to this the AsparosideF-NSP15 endoribonuclease complex remains around 2 Å throughout the MDS ( Figure  11 , Green).

To further understand these three complexes and the effect of ligand binding on the active site of the receptor RMSF of these complexes was studied ( Figure 12 ). The RMSF for Asparoside-C, Asparoside-D, Asparoside-F and Remdesivir bound D-NSP15 endoribonuclease protein amino acid residues had significantly low variations. The residues in the pocket forming the active site for binding of ligands had a low RMSF between 1-2 Å, this could be attributed to the low b-factor in that region owing to its structural confirmation.

The ligand RMSD for Asparoside-C, Asparoside-D, Asparoside-F and Remdesivir bound D-NSP15 endoribonuclease protein was studied to understand the conformational changes occurring during the MDS, the plot for ligand RMSD is presented in Figure 13 .

The phytochemical Asparoside-C had several conformational changes during the MDS which is evident from its evolving RMSD which fluctuates between 0 to 3 Å (Figure 13 , black). Asparoside-D equilibrated and the RMSD remains low for first 50 ns between 0.5 to 2 Å, after 50 ns the RMSD gradually rises to 2.5 and converges for good till the end of MDS (Figure 13 , Red). Asparoside-F, is an interesting case, the ligand RMSD changes from initial 1 Å to a narrow range of 1.75 to 2.0 Å, this reflects the limited access to conformational change due to either strong ligand binding due to hydrogen bonds or van Der Waals interactions or due to structural limitations. However, after 50 ns the molecule regains its original RMSD and fluctuates between 0.5 to 1.5 equilibrating till the end of MDS (Figure 13 , Green). More interesting changes were observed with Remdesivir, initially the there were fluctuation between 0-5 Å for first 10 ns and then the RMSD suddenly peaked to 20 Å and later stabilized for rest of the simulation period (Figure 13 , Blue). This high fluctuation in RMSD suggest major conformational and binding mode of Remdesivir.

To shed more light into this compound, interaction between receptor and ligand was explored by visualizing three-dimensional structure, given in Figure 14 . The initial conformation at the start of the MDS production phase was analyzed, the oxane ring of Asparoside-F formed a hydrogen bond with the His338 residue of the endonuclease protein active site (Asparoside-F-O14-His338-NE2, Bond length 2.91 Å) (Figure 14a ). After the completion of the MDS, this interaction is lost and the oxane rings form new bonds between Asparoside-F and Endoribonuclease (Asparoside-F-O12-Glu340-OE2, Bond length 2.65 Å) and (Asparoside-F-O20-Glu229-OE1, Bond length 2.65 Å) (Figure 14b ). The newly formed bonds and shorter bond Remdesivir-SARS-CoV-2 NSP15 Endoribonuclease shows high Ligand-RMSD but no major fluctuation in protein RMSD and RMSF. This was studied in detail and the trajectory was analyzed visually for the conformational evolution. It was found that the ligand breaks away from the complex around 12 ns breaking the interactions with Lys291, Val293 and Ser295 (Figure 15a and b) . During rest of the MDS, it does not form any further interactions with the receptor. Its positional and conformation changes (Figure 15c and d) contributed to the higher RMSD as observed in Figure 12 .

MM-GBSA analysis was performed on all of the six protein-ligand complexes to evaluate the affinity of ligands to the target protein receptors. The MM-GBSA based binding free energy (DG bind ) calculations were performed on the 100 ns long MDS trajectories. The binding energies measured by this method are more efficient than the GlideScore values for selection of protein-ligand complexes. The major energy components, such as H-bond interaction energy (DGbind_Hbond), coulomb or electrostatics interaction energy (DGbind_Coul), covalent interaction energy (DGbind_Cov), lipophilic interaction energy (DGbind_Lipo), electrostatic solvation free energy (DGbind_Solv) and van der Waals interaction energy (DGbind_vdW) altogether contribute to the calculation of MM-GBSA based relative binding affinity. The binding energies and the contributing factors calculated for the MDS trajectories are mentioned in the Table 4 . Out of the six complexes studied three complexes belonged to the SARS CoV-2 spike receptor-binding domain and other three bound to SARS-CoV-2 NSP15 Endoribonuclease. The Asparoside-D bound to SARS CoV-2 spike receptor-binding domain showed high binding free energies with DG bind of -66.49.41 Kcal/mol with high contribution form the DG gas (-115.98 kcal/mol) and DG GB (54.61 Kcal/mol). Asparoside-C had the second highest binding free energy (DG bind of -62.61 Kcal/mol). The DG Surf (-5.84 Kcal/mol) for Asparoside-C indicates its better surface accessibility compared to other molecules. Remdesivir showed least binding energy (DG bind of -18.45 Kcal/mol) in the studied molecules bound to the S protein. In the second group of compounds binding the SARS-CoV-2 NSP15 endoribonucelase, Asparoside-D shows a very good binding free energy (DG bind ¼ -72.46 Kcal/mol), second highest was Asparoside-F with binding free energy (DG bind ¼ -55.19 Kcal/mol). The high binding affinity is also justified by the fact that during MDS the molecule optimizes its conformation for better fitting with the receptor active site which is clearly evident with Asparoside-F bound to the SARS-CoV-2 NSP15 Endoribonuclease. The most significant finding was the comparison of Remdesivir binding to the SARS-CoV-2 NSP15 Endoribonuclease, it has a very poor binding free energy (DG bind ¼ -0.12 Kcal/mol). These results corroborate with results obtained from ligand RMSD analysis of Remdesivir. It suggests that binding of Remdesivir to SARS-CoV-2 NSP15 Endoribonuclease is weaker compared the phytochemicals under investigation.

The diverse pharmacological activities of Asparagus racemosus are well documented. However, in current outbreak of SARS-CoV-2 its anti-viral activities are immensely important. All together due to its anti-viral, immunomodulatory, antiinflammatory and antioxidant potentials, this plant and its phytoconstituents could become vital therapeutic alternative in COVID 19. On this line, in the present work, the potential of major phytochemicals from Asparagus racemosus was explored through in silico methods. The docking studies pointed out the phytoconstituents viz. Asparoside-C, Shatavarin-I, Asparoside-D, Shatavarin-X and Racemoside-A having the lowest binding free energies with SARS-CoV-2 spike receptor-binding domain and Asparoside-C, Asparoside-F, Rutin, Asparoside-D, Racemoside-A with NSP15 endoribonuclease. In order to gain deeper insights of binding affinities and binding modes of top phytoconstituents, 100 ns molecular dynamics simulations and MM-GBSA binding free energy calculations were carried out. It was found that the Asparoside-C-spike protein complex achieves the conformational stability at 12.5 ns with comparably low RMSD of 2.5 Å. In case of NSP15 endoribonuclease, Asparoside-F was found better stabilized with RMSD of 2 Å throughout the MDS. Compared to docking studies, MM-GBSA binding free energy calculations and estimates of various DG parameters are more accurate. These calculations suggest that the Asparoside-D and Asparoside-C bound to SAR CoV-2 spike receptor-binding domain with DG bind of -66.49 and -62.61 Kcal/mol respectively are most favorable while Asparoside-D and Asparoside-F bound to NSP15 endoribonucelase with DG bind -72.46 and -55.19 Kcal/mol respectively are most favorable in terms of binding interactions and À0.12 (6.89) Ã All energies are in Kcal/mol with standard deviation in parenthesis. DE VDW ¼ van der Waals contribution from MM; DE ELE ¼ electrostatic energy as calculated by the MM force field; DG GB ¼ the electrostatic contribution to the solvation free energy calculated by GB; DG Surf ¼ solvent-accessible surface area; DG Sol ¼ solvation free energy; DG gas ¼ gas phase interaction energy; DG bind ¼Binding free energy. binding free energies. In this study we also found that the standard drug Remdesivir, which is being used as treatment for COVID-19 has less binding affinity to the S protein compared to other phytochemicals. In case of NSP15 endoribonuclease, it shows very low binding affinity compared to the phytochemicals. The binding affinity of Remdesivir is more towards the SARS CoV-2 S protein compared to NSP15 endoribonuclease which suggest its site of binding to be RBD domain of S protein. The future scope of our findings lies in the experimental verification of these results towards the potential of Asparoside-D, Asparoside-C and Asparoside-F in inhibiting key SAR CoV-2 proteins. These results could be exploited further to design future ligands from natural sources.

No potential conflict of interest was reported by the author(s). 

Using integrated computational approaches to identify safe and rapid treatment for SARS -CoV-2

Plant profile, phytochemistry and pharmacology of Asparagus racemosus (Shatavari): A review

Structural and functional analyses of the severe acute respiratory syndrome coronavirus endoribonuclease Nsp15

Molecular analysis of the genus Asparagus based on matK sequences and its application to identify A. racemosus, a medicinally phytoestrogenic species

Novel 2019 coronavirus structure, mechanism of action, antiviral drug promises and rule out against its treatment

An exclusive 42 amino acid signature in pp1ab protein provides insights into the evolutive history of the 2019 novel human-pathogenic coronavirus (SARS-CoV-2)

Toona sinensis Roem tender leaf extract inhibits SARS coronavirus replication

SARS-CoV-2 host entry and replication inhibitors from Indian ginseng: In silico approach

Design, synthesis and anticancer studies of novel aminobenzazolyl pyrimidines as tyrosine kinase inhibitors

Fragment tailoring strategy to design novel chemical entities as potential binders of novel corona virus main protease

Transmission and evolution of the Middle East respiratory syndrome coronavirus in Saudi Arabia: A descriptive genomic study

An "Old" protein with a new story: Coronavirus endoribonuclease is important for evading host antiviral defenses

SARS-CoV-2 RNA dependent RNA polymerase (RdRp) targeting: An in silico perspective

Drug repurposing for coronavirus (COVID-19): In silico screening of known drugs against coronavirus 3CL hydrolase and protease enzymes

Andrographolide as a potential inhibitor of SARS-CoV-2 main protease: An in silico approach

Coronavirus (COVID-19) update: FDA issues emergency use authorization for potential COVID-19 treatment

In-silico approaches to detect inhibitors of the human severe acute respiratory syndrome coronavirus envelope protein ion channel

Herbal drug microscopy

Computational Design of ACE2-Based Peptide Inhibitors of SARS-CoV-2

A review on the cleavage priming of the spike protein on coronavirus by angiotensin-converting enzyme-2 and furin

Antiviral activity of glycyrrhizic acid derivatives against SARS-coronavirus

Enhancing immunity in viral infections, with special emphasis on COVID-19: A review. Diabetes and Metabolic Syndrome

Discovery of potential multi-target-directed ligands by targeting host-specific SARS-CoV-2 structurally conserved main protease

Design, synthesis, pharmacological evaluation and molecular docking studies of substituted oxadiazolyl-2-oxoindolinylidene propane hydrazide derivatives

Neuroprotective effect of Asparagus racemosus root extract via the enhancement of brain-derived neurotrophic factor and estrogen receptor in ovariectomized rats

TanIdentification of natural compounds with antiviral activities against SARS-associated coronavirus

Anti-SARS coronavirus 3C-like protease effects of Isatis indigotica root and plant-derived phenolic compounds

Evaluation of antibacterial activity of Asparagus racemosus Willd. root

Designing of improved drugs for COVID-19: Crystal structure of SARS-CoV-2 main protease Mpro

Insights into structural dynamics of allosteric binding sites in HCV RNA-dependent RNA polymerase

Computational studies of drug repurposing and synergism of lopinavir, oseltamivir and ritonavir binding with SARS-CoV-2 protease against COVID-19

Characterization of spike glycoprotein of SARS-CoV-2 on virus entry and its immune cross-reactivity with SARS-CoV

LEDGF/p75 IN interaction inhibitors: In silico studies of an old target with novel approach

Ayurveda rasayana in prophylaxis of COVID-19

Recent discovery and development of inhibitors targeting coronaviruses

Antiviral and immune boosting activities of different medicinal plants against Newcastle disease virus in poultry

Biflavonoids from Torreya nucifera displaying SARS-CoV 3CL(pro) inhibition

In-silico homology assisted identification of inhibitor of RNA binding against 2019-nCoV N-protein (N terminal domain)

Antibacterial activity of chemical constituents isolated from asparagus racemosus

Asparagus racemosus: A review on its phytochemical and therapeutic potential

Identification of bioactive compounds from glycyrrhiza glabra as possible inhibitor of SARS-CoV-2 spike glycoprotein and non-structural protein-15: A pharmacoinformatics study

An in-silico evaluation of different Saikosaponins for their potency against SARS-CoV-2 using NSP15 and fusion spike glycoprotein as targets

Public Health Approach of Ayurveda and Yoga for COVID-19 Prophylaxis

Anticandidal activity of Asparagus racemosus

Asparagus racemosus bZIP transcription factor-regulated squalene epoxidase (ArSQE) promotes germination and abiotic stress tolerance in transgenic tobacco

Effect of ethanolic extract of Solanum virginianum Linn. on neuropathic pain using chronic constriction injury rat model and molecular docking studies

Stilbene-based natural compounds as promising drug candidates against COVID-19

Traditional Chinese medicine herbal extracts of Cibotium barometz, Gentiana scabra, Dioscorea batatas, Cassia tora, and Taxillus chinensis inhibit SARS CoV replication

In silico screening of Chinese herbal medicines with the potential to directly inhibit 2019 novel coronavirus

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