key: cord-0825147-xs44uinq
authors: Joshi, Chinmayi; Chaudhari, Armi; Joshi, Chaitanya; Joshi, Madhvi; Bagatharia, Snehal
title: Repurposing of the herbal formulations: molecular docking and molecular dynamics simulation studies to validate the efficacy of phytocompounds against SARS-CoV-2 proteins
date: 2021-05-14
journal: Journal of biomolecular structure & dynamics
DOI: 10.1080/07391102.2021.1922095
sha: 0b774722fdfb7b40f5ea05dbf84b91eee3989c7a
doc_id: 825147
cord_uid: xs44uinq

Herbal formulations mentioned in traditional medicinal texts were investigated for in silico effect against SARS-COV-2 proteins involved in various functions of a virus such as attachment, entry, replication, transcription, etc. To repurpose and validate polyherbal formulations, molecular docking was performed to study the interactions of more than 150 compounds from various formulations against the SARS-CoV-2 proteins. Molecular dynamics (MD) simulation was performed to evaluate the interaction of top scored ligands with the various receptor proteins. The docking results showed that Liquiritic acid, Liquorice acid, Terchebulin, Glabrolide, Casuarinin, Corilagin, Chebulagic acid, Neochebulinic acid, Daturataturin A, and Taraxerol were effective against SARS-COV-2 proteins with higher binding affinities with different proteins. Results of MD simulations validated the stability of ligands from potent formulations with various receptors of SARS-CoV-2. Binding free energy analysis suggested the favourable interactions of phytocompounds with the recpetors. Besides, in silico comparison of the various formulations determined that Pathyadi kwath, Sanjeevani vati, Yashtimadhu, Tribhuvan Keeratiras, and Septillin were more effective than Samshamni vati, AYUSH-64, and Trikatu. Polyherbal formulations having anti-COVID-19 potential can be used for the treatment with adequate monitoring. New formulations may also be developed for systematic trials based on ranking from these studies. Communicated by Ramaswamy H. Sarma.

Coronavirus disease (COVID-19) is a pandemic disease caused by the novel Coronavirus and has been declared as a Global Public Health Emergency by World Health Organization (WHO). As of January 17, 2021, the disease has spread worldwide with 93,217,287 cumulative cases and 2,014,957 deaths, and in India, 10,557,985 cumulative confirmed cases and 1,52,274 cumulative deaths are reported (World Health Organization, 2020) . This pandemic situation becomes more challenging to handle in low and middle-income countries because a large proportion of individuals may be at increased risk of infection and face difficulties in receiving quality health services (Dong et al., 2020; Gupta et al., 2020) . The governments have prompted various strategies for the prevention and management of disease, but it has become one of the forefronts of health challenges due to the lack of the approved targeted therapy . Modern molecules such as Remdesivir, Lopinavir/Ritonavir, Favipiravir/Umiferovir, and anti-malarial drugs such as chloroquine, hydroxychloroquine have been reported to have anti-COVID-19 potential, however, reports about their ineffectiveness against COVID-19 is also published (Geleris et al., 2020; Marzolini et al., 2020; Wang et al., 2020) . Researchers are actively involved in searching for novel drug molecules alternatives to existing drugs, and pluralistic knowledge systems available worldwide may offer a bright ray of hope (Rastogi et al., 2020) . Many of the traditional systems such as traditional Chinese medicine, Kampo, traditional Korean medicine has been practiced in some areas of the world for the treatment of SARS-CoV diseases (Lee et al., 2020; Komuro, 2017; Yang et al., 2020) .

In India, The Ministry of Ayurveda, Yoga, and Naturopathy, Unani, Siddha and Homoeopathy (AYUSH), a federal government organization that is actively involved in research and development on Indian traditional medicine (http://ayush. gov.in/) presented a plausible plan of action for Ayurvedic intervention which is not limited to prophylaxis alone but also includes the therapeutic and integrative model of care.

The crystal structures of target proteins were downloaded in pdb (gz) format from the RCSB protein data bank (PDB) (https://www.rcsb.org/). A list of the selected protein targets is given in the supplementary material (Table S1 ). Protein structures were prepared by removing water molecules and heteroatoms followed by the conversion in .pdb format using BIOVIA Discovery Studio Visualizer v20.1.0.19295 (Dassault Syst emes, San Diego, CA, USA). The protein structures were prepared for molecular docking by adding polar hydrogen atoms and Kollman charges on them. This procedure was followed by the conversion of the pdb structure of protein into pdbqt and grid preparation using AutoDock Tools version 1.5.7 (ADT; Scripps Research Institute, La Jolla, San Diego, USA).

Chemical structures of more than 150 compounds from various polyherbal formulations and 14 FDA-approved drugs were retrieved from the PubChem database (https://pubchem.ncbi.nlm.nih.gov/) in the Spatial Data File (.sdf) format). Chemical structures of ligands in the .sdf format were converted to the .pdb format using BIOVIA Discovery Studio Visualizer. Ligand structures were prepared by adding nonpolar hydrogens, Gasteiger changes, and rotatable bonds and converted in .pdbqt format using AutoDock Tools. A list of formulations is given in supplementary material (Table S2 ).

Target proteins and ligands were docked using AutoDock Vina (The Scripps Research Institute, La Jolla, San Diego, USA) (Trott and Olson, 2010) . A three-dimensional grid box to define binding pocket was set using MGL tools. Binding pockets were optimized for each protein. Docking of each ligand was done to each receptor with grid coordinates and grid boxes of certain sizes for each receptor using AutoDock Vina installed in the ToolShed Galaxy server (The Galaxy Team, 2014) was used for molecular docking. Ten binding poses were generated for each ligand and protein. Results were obtained and ranked by the binding affinity (Shukla et al., 2019) , which was predicted as negative Gibbs free energy (DG) scores (kcal/mol). The docked ligand poses and protein were saved as a single file using PyMOL (DeLano Scientific LLC, Palo Alto, CA, USA and the interaction of proteins and ligands were analyzed using the BIOVIA Discovery Studio Visualizer. Protein-ligand docking poses were analyzed by studying the hydrogen-bond interactions, hydrophobic interactions, and bonding distances. The interaction of each ligand with the target proteins was studied and favorable confirmations were selected.

The protein-ligand complex structure of SARS-CoV-2 proteins and top 4 ranked ligands were subjected to MD simulation using Desmond Molecular Dynamics System, version 6.2, D. E. Shaw Research, New York (Bowers et al., 2006) . MD simulation of Remdesivir with various receptors was also performed. A 100 ns simulation was carried out for each complex. Protein-ligand complexes were prepared for simulation by adding hydrogens, filling missing side chains or whole residues. Refinement of prepared structure file was done. In this step, protonation and tautomeric states of variable residues were changed manually, which was followed by the generation of a solvated system for simulation. TIP3P water model was used for the solvation of complexes followed by neutralization by adding ions. After neutralization, energy minimization was performed. OPLS_2005 force field parameters were used for simulations. To perform simulations, the temperature was 300 K, the pressure was 1.01325 bar, and cut off radius was 5 Å. Root mean square deviation (RMSD), root mean square fluctuation (RMSF), the number of hydrogen bonds, and the radius of gyration (Rg) were calculated.

Protein-ligand confirmations and major global motions upon ligand binding were studied by principle component analysis (PCA). Bio3D, an R package was used to perform PCA. After removal of the rotational and translational movements, the positional covariance matrix of atomic coordinates and its eigenvectors were computed by overlaiding coordinates onto a reference structure followed by the diagonalization on the estimated symmetric matrix by creating the diagonal matrix of eigenvalues. Columns were the eigenvectors corresponding to the direction of motion relative to the initial coordinates in the diagonal matrix (Grant et al., 2006) .

The molecular mechanics energies combined with the generalized Born and surface area continuum solvation (MM/GBSA) method was used to calculate the free energy of the binding of ligands to proteins. The free energy of solvation and MM potential energy were calculated using Prime 4.0 (Schr€ odinger, LLC, New York, NY). The last 10 ns of the MD trajectories were taken for MM/GBSA (Jacobson et al., 2002 (Jacobson et al., , 2004 .

Docking analysis generated binding affinities of compounds with SARS-CoV-2 proteins. A higher negative value of binding affinity was considered as the most favorable interaction. The compounds were showing varied docking scores with different target proteins, and the average binding affinity of compounds was calculated to get a holistic picture of the docking results (Related file, Docking results of compounds). Docking of drugs such as lopinavir/ritonavir, remdesivir, favipiravir, chloroquine and hydroxychloroquine, favipiravir, darunavir, ivermectin was also done. Among the drugs, darunavir showed the highest average docking score followed by azithromycin. Docking results of existing drugs are shown in Table 1 .

On the basis of the docking results, the top 10 compounds were Liquiritic acid (-13.6), Terchenulin (-13), Liquorice acid (-12.7), Glabrolide (-11.6), Casuarinin (-11.6), Corilagin (-11), Chebulagic acid (-10.8), Neochebulinic acid (-10.9), Daturataturin A (-10.7) and Taraxerol (-10.6). Besides -7.0 -6.2 -6.9 Lopinavir/Ritonavir -6.5 -5.2 -4.1 -9.4 -5.4 -6.1 -3.8 -8.1 -5.8 -6.7 -7.7 -8.4 -7.2 -7.3 -6.9 -6.0 -6.6 Azithromycin -9.4 -5.4 -6.3 -12.2 -8.5 -10.2 -6.5 -10.4 -9.3 -11.9 -11.2 -10.6 -11.0 -10.2 -9.1 -8.9 -9.1 Ivermactin -8.8 -6.9 -8.4 -9.5 -10.5 -9.2 -5.9 -9.6 -7.7 -9.7 -10.3 -10.9 -10.6 -9.5 -7.6 -6.4 -8.3 Aspirin -5.6 -5.3 -5.5 -5.8 -6.1 -4.9 -3.8 -5.9 -5.6 -6.6 -5.3 -6.1 -6.7 -5.9 -5.1 -4.7 -5.5 Darunavir -10.5 -9.2 -9.0 12.1 -10.2 -10.7 -7.3 -12.6 -9.4 -12.7 -11.0 -11.5 -12.1 -11.2 -9.5 -8.7 -10.2 Galidesivir -6.3 -4.5 -4.8 -7.5 -5.6 -6.3 -4.3 -7.4 -6.2 -7.5 -7.0 -7.8 -7.5 -7.5 -5. the top 10 compounds, numerous compounds have also shown a good binding affinity with target proteins ranging from -10.5 to -7.5. A list of the top 25 compounds having binding affinity ranging from -13.6 to -9.9 is given in Table 2 . Given the interactions of compounds with the individual proteins, Liquiritic acid showed the highest binding affinities with all the proteins selected for the study. Docking scores of Liquiritic acid ranged from -16 to -8.9, where the highest affinity was observed with NSP15 Endoribonuclease (-16) of SARS-COV-2 followed by the nucleocapsid protein (-15.1). Liquiritic acid showed good interactions with NSP12, main protease, spike glycoprotein, RNA-dependent RNA polymerase, NSP3 also. After Liquiritic acid, Terchebulin also showed good scores with target proteins. It showed the highest score with RBD/ACE2-B0AT1 complex (-15.5), main protease (-15), followed by spike glycoprotein (-14.1).

The ranking of plants based on the binding affinity of compounds is given in Table 3 . Plant-wise sorting of docking results showed that the compounds from the plants i.e. Glycyrriza glabra (Yashtimadhu), Terminalia chebula (Harde), Tinospora cordifolia (Galo), Phyllanthus emblica (Amla), Semecarpus anacardium (Bhilamo), Andrographis paniculata (Kariyatu), Dhatura metel (Ganthalo Daturo), Alstonia scholaris (Saptaparna), Aconitum ferox (Vachnag), Azadirachta indica (Neem), Withania somnifera (Ashwgandha) were showing good binding affinities ranging from -13.6 to -9.6. Compounds from Tinospora cordifolia, Piper longum, Azadirachta indica, Phyllanthus emblica, Oroxylum indicum, Aegele marmelos and Ocimum sanctum were also showing good interactions with proteins in terms of binding affinity.

Formulation-wise segregation of the docking results showed that Yashtimadhu, Pathyadi kwath, Sanjeevani vati, Septillin, and Tribhvan Keerti ras have plants and phytochemicals in formulations showing higher binding scores in silico. A list of formulations containing the top 10 compounds is given in Table 4 .

Total eight ligand-protein complexes (including reference ligand Remdesivir) were subjected to MD simulation. A box of water with Na þ and Clions were used to prepare proteinligand complexes for the simulation. The stability and confirmation of the protein-ligand complexes were assessed by simulating for 100 ns and analyzed by plotting the RMSD, RMSF, a radius of gyration, and H-bonds graph for 100 ns.

RMSD plots of ligands and receptors showed that the ligands were having multiple binding orientations. Figure 1 shows the RMSD values of a protein (left Y-axis), and ligand (right Y-axis). Lig fit Prot and Lig fit Lig show the RMSD of a ligand when the protein-ligand complex is first aligned on the protein backbone of the reference and its reference confirmation, respectively. The RMSD plot of Liquiritic acid-6VSB complex showed the changes in RMSD values larger than 3 Å, which indicates the large conformational change in the protein structure during the simulation. However, the RMSD values of protein and ligand remained intact between 3.

Yashtimadhu -12.7 -11.3 -11.6 -14.7

-12.5 -12.5 -8.4 -14.6 -13.2 -13.5 -15.0 -14.3 -11.8 -13.2 -12.4 -12.7 -12.1

Yashtimadhu -11.6 -11.8 -13.0 -12.0 -12.9 -11.0 -8.8 -11.0 -10.0 -12.2 -11.4 -10.6 -13.1 -13.2 -11.8 -11.0 -12.4

Harde -11.6 -11.0 -11.8 -13.0 -12.0 -11.5 -7.1 -13.2 -12.0 -12.8 -13.0 -13.2 -13.6 -9.9 -9.6 -9.8 -12.3

6.

Amla -11.0 -10.0 -12.0 -11.9

-11.6 -9.2 -7.4 -12.0 -11.2 -13.6 -11.8 -12.0 -11.7

-10.4 -11.5 -9.9 -9.8

Neochebulinic acid Terminalia chebula -12.0 -11.5 -11.1 -10.9

-12.6 -9.5 -11.0 -11.3

10.

Gokarni -10.6 -12.1 -9.8 -10.2 -9.3 -9.7 -8.9 -9.8 -10.3 -10.4 -10.6 -11.4 -11.3 -11.2 -12.3 -10.9 -11.3

Hispaglabridin A

Yashtimadhu -10.5 -8.6

-9.4 -11.2 -9.8 -9.9 -6.5 -11.2 -10.2 -10.9 -11.3 -11.0 -11.6 -12.0 -12.3 -11.1 -11.7

Chebulinic acid -9.8 -6.4 -13.2 -11.6 -11.6 -11.2 -12.1 -11.2 -10.7 -9.9 -9.7 -8.6

17. 0-15 ns, 20-40 ns, 60-70 ns, and 85-90 ns of the simulation (Figure 2(A) ). The RMSD plots of Tinocordioside-6M3M, Daturataturin-6NUR, and Corilagin-6VWW complexes showed that the complexes equilibrated after 80, 60, and 50 ns respectively (Figure 2(B-D) ). Noticeably, in all the complexes except liquiritic acid-6VSB complex, stabilized RMSD values indicated that the system had equilibrated after 100 ns of simulation. RMSD plots of reference compound (Remdesivir) with the selected proteins (Figure 3) , indicated that the reference compound formed more stable complexes as compared to ligands. However, in the case of Remdesivir-6VSB complex, RMSD values were found to be increasing even after 100 ns simulation.

The RMSF of the protein-remdesivir complexes and the protein-ligand complexes were calculated and shown in Figure 4 . RMSF plots were generated for all the complexes to analyze the residual atomic fluctuations of the protein atoms in the presence of ligands. In RMSF plots, fluctuations were observed during the simulation which was indicated by peaks. It is commonly known that the helices and sheets show less flexibility, whereas loops and turns show higher flexibility (Shukla et al., 2019) . Here in the present study, the average RMSF value for Remdesivir-6VSB, 6NUR, 6M3M, and 6VWW complex was 6.57, 6.09, 20.79, and 5.26 Å, respectively, whereas the average RMSF value for liquiritic acid-6VSB, daturataturin-6NUR, tinocordioside-6M3M, and the colilagin-6VWW complex was 12.41, 2.63, 12.73, and 4.62 Å, respectively. The RMSF values indicate that Remdesivr showed lesser fluctuations than daturataturin and corilagin. RMSF plots of Daturataturin-6NUR and Corilagin-6VWW complexes were found to fluctuate more from their mean structure. However, Corilagin-6VWW complex was stable after the 50 ns simulation, whereas initially RMSF value of Liquiritic acid-6VSB complex was decreased after 35 ns and was found to be almost stable between 60 ns simulation.

Intermolecular H-bonds are also important factors to be considered to check the stability of protein-ligand complexes during the simulations. The number of H-bonds was also analyzed using simulation trajectories. H-bonds of ligands (reference and phytochemicals) are shown in Figure 5 . Among the four ligands (liquiritic acid, corilagin, tinocordioside, daturataturin), corilagin achieved the highest H-bonds (average $4). Reference ligand showed the lesser number of H-bonds as compared to phytochemicals. In total, the constant RMSD of protein-ligand complexes during the simulation period was achieved by the contribution of the number of H-bonds. From the results, it can be said that ligand-binding affinity alters with the change of the number of H-bonds.

In addition to H-bonds, hydrophobic and ionic interactions also play an important role in protein-ligand interactions. A schematic of detailed ligand atom interactions with the protein residues is given in Figure 6 . Liquiritic acid docked with 6VSB receptor showed the interactions to 3 residues THR302 (47%), LEU303 (48%), LYS304 (45%), whereas remdesivir showed H-bond with ILE312 (65%). In the case of 6M3M receptor, Tinocordioside showed interaction with VAL159 with 32% occupancy, whereas Remdesivir showed interaction with ASN76, SER79 and ILE75. In the case of 6NUR, Daturataturin showed a better interaction as compared to remdesivir. The radius of gyration (Rg) values of the protein and ligand complexes were also calculated to analyze the compactness of protein-ligand complexes (Figure 7 ). Figure 7(A) shows the Rg values of reference compounds. The average Rg values of remdesivir with 6VSB, 6M3M, 6NUR, and 6VWW were 0.52, 0.45, 0.48, and 0.46, respectively. The reference complexes showed less Rg values compared to phytocompound-protein complexes, suggested that it generates more compact complexes compared to the phytocompounds.

In the case of 6VSB and liquiritic acid, the Rg value of protein and ligand was found to be 0.90 and 0.88 nm respectively after 40 ns simulation. Protein 6VSB showed fluctuation from $1.43 to $0.91 nm, whereas ligand liquiritic acid showed fluctuation from $1.36 to $0.91 nm (Figure 7(B) ). In the case of 6M3M and tinocordioside, Rg values were found to be increased after 100 ns simulation from $1.0 to $1.20 nm (Figure 7 (C)), whereas, in the case of 6NUR and daturataturin, the Rg value of protein and ligand was found to be 1.87 and 1.86 nm respectively after 100 ns simulation. Both the protein 6NUR and ligand daturataturin showed fluctuation from $1.87 to $1.80 nm (Figure 7(D) ). Rg values of 6VWW and corilagin were found to fluctuate from $1.63 to $1.67 nm (Figure 7(E) ). However, the fluctuation pattern of the Rg values was the same for ligands and their complex proteins and overall Rg plot analysis showed that ligands Rg curves significantly followed the fluctuation pattern of the respective proteins, indicating the ability of the ligand to stabilize the proteins.

PCA method was used to study the concerted motions after ligand binding. PCA plot and Porcupine plot of reference complex and liquiritic acid-6VSB complex are given in Figure 8 . It shows the statistically meaningful confirmations in both the complexes. The principal motions and the vital motions required for conformational changes were identified in both the complexes. Two different groupings along the PC1 plane were observed, whereas the groupings observed along the PC2 and PC3 planes were not separated. These results indicate the presence of a non-periodic conformational change and periodic global motions. However, the porcupine plot showed the difference between the two complexes more accurately. The direction and magnitude of the motion were presented by the arrows on the protein-ligand complexes. Differences in the motions were observed in the Remdesivir-6VSB complex as compared to Liquiritic acid-6VSB complex. Prominent motions were observed in the Remdesivir-6VSB complex, whereas a large deviation pattern was observed in the Liquiritic acid-6VSB complex.

The principal motions and the vital motions required for conformational changes in other protein-ligand complexes are shown in PCA plots (Figure 9 ). The graph showed that Remdesivir-protein complexes showed a much stable cluster as compared to all complexes. PCA plots were also generated for the rest of the six complexes. Periods jumps were observed in almost all complexes. 6NUR_Daturataurin and 6NUR_Remdesivir were showing positive distinct correlated motions (blue), may be defined these complexes have stable inhibitors for respected protein.

All the protein-ligand complexes with reference ligand were subjected for calculation of binding free energy by using MM-GBSA method. All the complexes including phytocompound-protein complexes and Remdesivir-protein complexes showed a negative binding affinity, indicating favorable interactions. Additionally, other interactions such as electrostatic energy, Van der Waals interactions, nonpolar solvation energy, etc. were also calculated and presented in Table 5 . In all the cases, the reference compound showed a better binding affinity as compared to phytocompounds. Other interactions i.e. solvation energy, contributed to the total interaction energy adversely, whereas electrostatic(coloumb), van der Waals energy contributed to the binding energy favorably. However, the results obtained from binding free energy analysis of phytocompounds-protein complexes indicate that the phytocompounds from various formulations can bind the SARS-CoV2 proteins, and they can generate stable complexes with proteins. Figure 10 , showing electrostatic interactions which were calculated using the PRIME energy calculation module of Schrodinger, whose surface depicting electrostatic energy ranging from values -30 kcal/mole (blue) to þ30 kcal/mole(red). In 6VSB-Remdesivir overall electrostatic energy was contributing more compare to 6VSB_Liquiritic acid. Where ligand seems to contribute electrostatic potential positive in both cases. Minimized structure obtained through MMGBSA, as shown in Figure 10 . 6VSB_Liquiritc acid seems to have more hydrogen bond interactions with amino acids GLN957, ASN953, and THR827 

Despite extensive advancement in modern medicine, the effectiveness of drugs such as remdesivir, hydroxychloroquine, favipiravir was found to be decreased due to the rapid emergence of SARS-CoV-2 drug resistance. These drugs can target cellular functions such as viral replication, terminal glycosylation of ACE2, etc., and affect the viral life cycle. Moreover, Padhi et al. (2021) performed a mutational mapping and provided insight into the functional outcomes of mutations in the remdesivir-binding site in the nsp12 subunit of RdRp. Their study suggested that very few mutations in nsp12 of SARS-COV-2 can lead the resistance against remdesivir. The drug resistance takes place due to the evolutionary pressure on the viral population, and the pressure provides a favorable condition to the subpopulation that has a relatively better fitness with reduced fidelity and emergence of mutation. In the presence of evolutionary pressure, the mutated viral subpopulation will become predominant. In face of such difficulties in the management of diseases, currently, a therapeutic requirement of alternative treatment is high. In this line, natural products can be trusted due to their chemical diversity, biological diversity, and drug-like properties. The traditional medicinal system uses natural products, which have been practiced around the globe for several years. In India, Ayurveda is one of the traditional medicinal systems, in which single or multiple herbs are being used for the treatment. A single herb or polyherbal formulation may contain more than one phytoconstituent which may act synergistically to achieve extra therapeutic effectiveness by maintaining healthy conditions and improving immunity. Simultaneously these compounds may target multiple functions of pathogens at the various stages of the infection. In total, traditional medicine can use the "multi-drugs and multi-targets" mode for the treatment of complex diseases which may be more effective than the individual drug (Padhi et al., 2021; Parasuraman et al., 2014) .

As finding solutions for the management of COVID-19 disease is the highest priority in the current pandemic situation, molecular docking studies were planned. SARS-CoV-2 genome encodes for 14 open reading frames and 16 protein replicase-transcriptase that consists of multiple enzymes essential for attachment and entry of the virus into host cells, replication and pathogenicity, virus transcription, infection to human cells as well as transcription process (Cheng, 2007; Dinesh et al., 2020; Gao, 2020; Kirchdoerfer & Ward, 2019; Yan, 2020) . In this study, we aimed to evaluate the anti-COVID-19 potential of ayurvedic polyherbal formulations that contain several natural compounds and plant extracts that have been reported for their antiviral potential (Churiyah et al., 2015; Ganjhu, 2015; Saha & Ghosh, 2012; Tan et al., 2013) . In the present study, the selection of polyherbal formulations was done based on the advisory of the Ministry of AYUSH for the management of COVID-19 spread in India. Various interventions for the prevention of disease have been proposed by AYUSH which are being used as prophylactic and symptomatic management for disease management (Guidelines for Ayurveda practitioners for covid 19; WHO, 2020).

The objective was to study the interactions of ligands present in polyherbal formulations with the proteins involved in COVID-19 and to analyze the stability of protein-ligand complexes. Compounds from polyherbal/polycomponent formulations were docked into the active site of the SARS-CoV-2 target proteins. Amongst the top five compounds, liquorice acid has been used in the treatment of hepatic disease for more than 40 years in Japan (Li et al., 2014) . It has also been reported for its potential to inhibit SARS-coronavirus (SARS-CoV) replication (Hoever, 2005) . Another compound casuarinin was also found to be effective against herpes simplex type 2 (HSV-2) in vitro (Cheng et al., 2002) .

Plant-wise sorting indicated that compounds from G. glabra, T. chebula, T. cordifolia, P. emblica were found to be effective in silico. In Ayurveda, these plants are extensively used for the treatment of various disorders. Yashtimadhu has been prescribed for respiratory and digestive disorders, and to improve immunity. Glycyrriza glabra reported to have activity against herpes simplex, Varicella zoster, Japanese encephalitis, influenza virus, vesicular stomatitis virus, type A influenza virus (Adam, 1997; Pompei et al., 1979 Pompei et al., , 1980 . In addition to that, authors (Cinatl et al., 2003) demonstrated in vitro effect of compounds present in Glycyrrhiza glabra against two clinical isolates of coronavirus (FFM-1 and FFM-2) from patients with SARS admitted to the clinical center of Frankfurt University, Germany. Likewise, T. chebula has also been reported for various biological properties such as anticancer, anti-inflammatory, antioxidant, anti-protozoal, antimicrobial activity. Badmaev and Nowakowski (2000) reported the protective activity of T. chebula against influenza A virus. Various studies reported the antiviral activity of T. chebula on different viruses (Ajala et al., 2014; Kolla et al., 2017; Oyuntsetseg et al., 2014) . Chebulagic acid from Terminalia chebula is reported to have potential as broad-spectrum antivirals for controlling viral infections which engage host cell surface glycosaminoglycans that play a role in initial binding to the host cell (Lin, 2013) . Plants such as Tinospora cordifolia, Piper longum, Azadirachta indica, Phyllanthus emblica, Oroxylum indicum, Aegele marmelos, and Ocimum sanctum have been used in various polyherbal formulations of Ayurveda and also reported for their antiviral properties against numerous viruses (Ghoke et al., 2018; Jiang, 2013; Lalita, 2002; Tiwari et al., 2010; Vellingiri et al., 2020) .

The results of formulation-wise sorting demonstrate the effectiveness of ayurvedic formulations. Few of these formulations and the plants used in formulations have been also mentioned in the Indian System of Medicine such as Siddha and Unani. Overall, the study demonstrates the efficacy of polyherbal formulations using modern tools.

In the context of epidemics and pandemics, complementary and alternative medicine recommends preventive and prophylactic measures to enhance immunity. The traditional medicinal systems recognize the effect of ecological and environmental conditions on the state of health of humans and make a great effort on the prevention of diseases.

Few plants such as Glycyrriza glabra (Yashtimadhu), Aegele marmelos (Bilva), Terminalia chebula (Harde), Tinospora cordifolia (Guduchi), etc. also reported having anti-inflammatory activity (Bag, 2013; Nirmala & Selvaraj, 2011; Rajaram, 2018; Upadhyay et al., 2010) . In severe COVID-19 patients, the uncontrolled release of pro-inflammatory cytokines is very common and the progression of the disease results in the loss of immune regulation due to exacerbation of the inflammatory components (Garcia, 2020) . Zhang et al. (2020) shared their experience about the anti-inflammation treatment of patients with severe COVID-19 and described the importance of anti-inflammatory drugs for the treatment in patients suffering from symptoms including multiple organ failure and acute respiratory distress syndrome (ARDS).

In the current outbreak of COVID-19, one other major concern is the D-dimer production. Few studies revealed the increase of D-dimer levels in severe COVID-19 patients Zhou et al., 2020) . Wichmann et al. (2020) performed autopsies on 12 patients with COVID-19. They showed deep vein thrombosis and pulmonary embolism (PE) in seven and four patients respectively. In this context, Ayurveda has also described the management of D-dimer formation. Herbs of the formulations such as Yashtimadhu, Pathyadi Kwath, Sanjeevani Vati, Samshamni Vati were reported to have antithrombotic activity. Additionally, compounds such as Glycyrrhizin, Corilagin, Withanoferin A also reported for antithrombotic activity (Ku & Bae, 2014; Lugun, 2018; Mendes-Silva, 2003; Ri & Ho Ju, 2018; Saleem, 2019) .

Molecular dynamics simulation results indicate that RMSD values of proteins and ligands followed a similar pattern except for liquiritic acid-6VSB complex. Overall, protein-ligand complexes were stable and RMSD fluctuations were observed due to the small size of ligands and partial binding pocket occupancy. The number of H-bonds indicated the higher stability of protein-ligand binding conformations (Londhe et al., 2019) . Peele et al. (2020) carried out 20 ns MD simulations on the main protease and inhibitors complex. Their study validated the stability of lopinavir, amodiquine, and theaflavin digallate in the protein binding pocket as potent binders. Similar types of studies have been reported about the screening of natural compounds having anti-COVID potential from Traditional Chinese Medicine (Romeo et al., 2020; Selvaraj et al., 2020) .

Overall, a collective approach of docking and MD simulation is beneficial for repurposing ayurvedic formulations for the management of the disease. Docking results and their correlation with in vitro studies reported in the literature generate reliable evidence to use these ayurvedic formulations for the management of COVID-19. Ancient literature such as Siddha, Unani, Ayurveda has the potential to prevent and treat COVID-19 which can be further tested. According to guidelines published by AYUSH, these formulations can be used for the management of the disease. This study provides evidence of the interaction of the natural compounds with the proteins involved in the disease using computational tools.

According to the present study, the formulations tested in this study had significantly higher binding efficacy against their SARS-CoV-2 targets. In silico results reveal that these formulations may be effective inhibitors of SARS-CoV-2 through their binding to the spike glycoprotein, RNAdependent RNA polymerase, protease, which can be further studied in vitro. Interactions of compounds with the target proteins suggest that these compounds will be efficacious in preventing both viral attachment and replication, as well as these formulations, can be directly used for patients having mild to moderate symptoms of COVID-19. SARS-CoV-2 is reported to have a high predilection in the pharyngeal epithelial cells. As the deliverability of these extracts to the pharyngeal regions is feasible by suitable oral formulations, they will be useful in the medical management of SARS-CoV-2 infections.

In vitro antiviral activity of indigenous glycyrrhizin, licorice and glycyrrhizic acid (Sigma) on Japanese encephalitis virus

Hepatitis C virus inhibitory hydrolysable tannins from the fruits of Terminalia chebula

Protection of epithelial cells against influenza A virus by a plant derived biological response modifier Ledretan-96

Anti-inflammatory, anti-lipid peroxidative, antioxidant and membrane stabilizing activities of hydroalcoholic extract of Terminalia chebula fruits

Dissemination of scientific software with Galaxy ToolShed

Scalable algorithms for molecular dynamics simulations on commodity clusters < SE-END>

Proceedings of the 2006 ACM/IEEE Conference on Supercomputing, November)

Antiherpes simplex virus type 2 activity of casuarinin from the bark of Terminalia arjuna Linn

Severe acute respiratory syndrome coronavirus as an agent of emerging and reemerging infection

Antiviral and immunostimulant activities of Andrographis paniculata

Glycyrrhizin, an active component of liquorice roots, and replication of SARS-associated coronavirus

Structural basis of RNA recognition by the SARS-CoV-2 nucleocapsid phosphoprotein

An interactive web-based dashboard to track COVID-19 in real time

Herbal plants and plant preparations as remedial approach for viral diseases

Structure of the RNA-dependent RNA polymerase from COVID-19 virus

Immune Response, Inflammation, and the Clinical Spectrum of COVID-19

Observational study of hydroxychloroquine in hospitalized patients with Covid-19

Evaluation of antiviral activity of Ocimum sanctum and Acacia arabica leaves extracts against H9N2 virus using embryonated chicken egg model

Bio3d: An R package for the comparative analysis of protein structures

Clinical characteristics of coronavirus disease 2019 in China

Guidelines for Ayurveda practitioners for COVID 19

The need for COVID-19 research in low-and middle-income countries

Antiviral activity of glycyrrhizic acid derivatives against SARS À coronavirus

On the role of crystal packing forces in determining protein sidechain conformations

A hierarchical approach to all-atom protein loop prediction

Anti-HBV active constituents from Piper longum

Structure of the SARS-CoV nsp12 polymerase bound to nsp7 and nsp8 co-factors

Terminalia chebula Retz.-an important medicinal plant

Japanese Kampo Medicines for the Treatment of Common Diseases: Focus on Inflammation

Antiplatelet, anticoagulant, and profibrinolytic activities of withaferin A

In vitro antiviral activity of bael (Aegle marmelos Corr) upon

A consensus guideline of herbal medicine for coronavirus disease 2019

Glycyrrhizic acid in the treatment of liver diseases: Literature review

Broad-spectrum antiviral activity of chebulagic acid and punicalagin against viruses that use glycosaminoglycans for entry

Investigation of molecular details of Keap1-Nrf2 inhibitors using molecular dynamics and umbrella sampling techniques

Evaluation of antithrombotic activities of Solanum xanthocarpum and Tinospora cordifolia

Effect of systemic inflammatory response to SARS-CoV-2 on lopinavir and hydroxychloroquine plasma concentrations

Antithrombotic effect of Glycyrrhizin, a plantderived thrombin inhibitor

Anti-inflammatory and anti-bacterial activities of Glycyrrhiza glabra L

Evaluation of direct antiviral activity of the Deva-5 herb formulation and extracts of five Asian plants against influenza A virus H3N8

High-throughput rational design of the remdesivir binding site in the RdRp of SARS-CoV-2: Implications for potential resistance. iScience

Polyherbal formulation: Concept of ayurveda

Molecular docking and dynamic simulations for antiviral compounds against SARS-CoV-2: A computational study

Glycyrrhizic acid inhibits virus growth and inactivates virus particles

Antiviral activity of glycyrrhizic acid

Anti-inflammatory profile of Aegle marmelos (L) Correa (Bilva) with special reference to young roots grown in different parts of India

COVID-19 Pandemic: A pragmatic plan for Ayurveda Intervention

The antithrombotic activity of corilagin purified from Korean herb-Phyllanthus ussuriensis

Targeting the SARS-CoV-2 spike glycoprotein prefusion conformation: Virtual screening and molecular dynamics simulations applied to the identification of potential fusion inhibitors

Tinospora cordifolia: One plant, many roles

Antithrombotic Potential of Ajwa Dates and Piper Nigrum: In Vitro &

Structure-based virtual screening and molecular dynamics simulation of SARS-CoV-2 Guanine-N7 methyltransferase (nsp14) for identifying antiviral inhibitors against COVID-19

Structural and energetic understanding of novel natural inhibitors of Mycobacterium tuberculosis malate synthase

Evaluation of antiviral activities of four local Malaysian Phyllanthus species against herpes simplex viruses and possible antiviral target

Abnormal coagulation parameters are associated with poor prognosis in patients with novel coronavirus pneumonia

In vitro antiviral activity of neem (Azardirachta indica L.) bark extract against herpes simplex virus type-1 infection

AutoDock Vina: Improving the speed and accuracy of docking with a new scoring function, efficient optimization, and multithreading

Tinospora cordifolia (Willd.) Hook. f. and Thoms. (Guduchi)-validation of the Ayurvedic pharmacology through experimental and clinical studies

COVID-19: A promising cure for the global panic

Remdesivir in adults with severe COVID-19: A randomised, double-blind, placebo-controlled, multicentre trial

Autopsy findings and venous thromboembolism in patients with COVID-19: A prospective cohort study

Clinical management of severe acute respiratory infection (SARI) when the COVID-19 disease is suspected: Interim Guidance (Version 1.2), 13 March 2020 released by World Health Organization

Structural basis for the recognition of SARS-CoV-2 by fulllength human ACE2

Traditional Chinese medicine in the treatment of patients infected with 2019-new coronavirus (SARS-CoV-2): A review and perspective

The traditional medicine and modern medicine from natural products

The use of anti-inflammatory drugs in the treatment of people with severe coronavirus disease 2019 (COVID-19): The experience of clinical immunologists from China

Clinical course and risk factors for mortality of adult inpatients with COVID-19 in Wuhan, China: a retrospective cohort study

The authors are thankful to Zuber Saiyed for providing technical support for the installation of software and tools for molecular docking and molecular dynamics simulations.

The authors declare no competing interests.