key: cord-0852956-54xtztx4
authors: Rudrapal, Mithun; Celik, Ismail; Khan, Johra; Ansari, Mohammed Azam; Alarousy, Randa Mohammed Ibrahim Ismail; Yadav, Rohitash; Sharma, Tripti; Tallei, Trina Ekawati; Pasala, Praveen Kumar; Sahoo, Ranjan Kumar; Khairnar, Shubham J.; Bendale, Atul R.; Zothantluanga, James H.; Chetia, Dipak; Walode, Sanjay G.
title: Identification of bioactive molecules from Triphala (Ayurvedic herbal formulation) as potential inhibitors of SARS-CoV-2 main protease (Mpro) through computational investigations
date: 2022-01-10
journal: J King Saud Univ Sci
DOI: 10.1016/j.jksus.2022.101826
sha: e6507ef74a4c5a5853230a06357d7269b108885a
doc_id: 852956
cord_uid: 54xtztx4

Severe acute respiratory syndrome coronavirus disease (SARS-CoV-2) induced coronavirus disease 2019 (COVID-19) pandemic is the present worldwide health emergency. The global scientific community faces a significant challenge in developing targeted therapies to combat the SARS-CoV-2 infection. Computational approaches have been critical for identifying potential SARS-CoV-2 inhibitors in the face of limited resources and in this time of crisis. Main protease (M(pro)) is an intriguing drug target because it processes the polyproteins required for SARS-CoV-2 replication. The application of Ayurvedic knowledge from traditional Indian systems of medicine may be a promising strategy to develop potential inhibitor for different target proteins of SARS-CoV-2. With this endeavor, we docked bioactive molecules from Triphala, an Ayurvedic formulation, against M(pro) followed by molecular dynamics (MD) simulation (100 ns) to investigate their inhibitory potential against SARS-CoV-2. The top four best docked molecules (terflavin A, chebulagic acid, chebulinic acid, and corilagin) were selected for MD simulation study and the results obtained were compared to native ligand X77. From docking and MD simulation studies, the selected molecules showed promising binding affinity with the formation of stable complexes at the active binding pocket of M(pro) and exhibited negative binding energy during MM-PBSA calculations, indication their strong binding affinity with the target protein. The identified bioactive molecules were further analyzed for drug-likeness by Lipinski’s filter, ADMET and toxicity studies. Computational (in silico) investigations identified terflavin A, chebulagic acid, chebulinic acid, and corilagin from Triphala formulation as promising inhibitors of SARS-CoV-2 M(pro), suggesting experimental (in vitro/ in vivo) studies to further explore their inhibitory mechanisms.

Coronavirus disease 2019 (COVID-19) is a highly infectious disease with a high rate of reproducibility. Severe acute respiratory syndrome coronavirus disease-2 (SARS-CoV-2) is the causative agent for COVID-19 . This deadly disease has infected nearly every country on the planet, with 266,946,588 confirmed cases and 5,282,046 fatalities, according to the latest reports of World Health Organization (https://covid19.who). With the majority of the world under lockdown and the alarming threat of thousands of deaths every day, rapid drug discovery for the identification of anti-COVID drugs is an urgent need of the hour. Conventional drug discovery processes could take years. To aid in the rapid discovery of new drugs, computational approaches have been effective in identifying potential drugs/medications for COVID-19 therapy. The use of computational algorithms such as molecular docking, molecular dynamics (MD) simulation and density functional theory (DFT) method has made the search for drugs of natural or synthetic origin that target the viral system, particularly SARS-CoV-2 main protease more feasible while causing fewer side effects with improved therapeutic performance in human populations. In this context, repurposing the traditional herbal remedies such as Ayurvedic formulations can play a significant role in the development of effective therapeutic entities. Moreover, the uses of phytochemicals over synthetic drugs are advantageous since they are safer and a higher dose is well tolerated.

Ayurvedic formulation Triphala, which literary means Tri (three) and Phala (fruit) in Sanskrit, is a 'tridoshic rasayan' that balances and rejuvenates the "tridosha", i.e., vata, pitta, and kapha, as well as promotes immunity and health. It is a polyherbal formulation made up of a proportional amount of dried fruit powder from three plants: Phyllanthus emblica (Indian gooseberry or Amalaki), Terminalia bellirica (Bibhitaki), and Terminalia chebula (Haritaki) (Peterson et al., 2017) . In the clinical practice of Ayurveda, the formulation (Triphala) is consumed as a 'health tonic'. It is also used to treat digestive problems such as gastric acidity, constipation and poor food assimilation. Additionally, several diseases such as fever, cough, jaundice, anemia, asthma, chronic ulcers, pyorrhea and leucorrhoea are also treated with Triphala. Additionally, obesity, inflammation, fatigue, cardiovascular disorders, liver dysfunction and ophthalmic problems are some other health ailments that Triphala is being recommended for. Moreover, it demonstrates cardiotonic properties by lowering serum cholesterol levels, reducing myocardial necrosis, enhancing blood circulation, and strengthening capillaries. Flavonoids, tannins, gallic acid, pectin, vitamin C, corilagin, and phyllaemblic compounds are present in in P. emblica fruits (Habib-ur-Rehman et al., 2007) . Terminalia bellirica extract contains β-sitosterol, tannins and phenolic acids such as ellagic acid, gallic acid, galloyl glucose, ethyl gallate, and chebulaginic acid (Kumar and Khurana, 2018) . The main bioactive active constituents present in T. chebula are chebulic, chebulinic acid, 2,4-chebulyl-D-glucopyranose, terflavin A, terchebin, ellagic acid, gallic acid, tannic acid, ethyl gallate, arjungenin, arjunglucoside I, punicalagin and luteolin. Recently, computational techniques are being extensively used to investigate the inhibitory potential of phytocompounds against SARS-CoV-2 (Singh et al., 2021; Sharma et al., 2021; Bhardwaj et al., 2021) . Numerous studies have been conducted to determine the efficacy of traditional medicine as anti-SARS-CoV-2 agents (Verma et al., 2020; Chakraborty et al., 2021; Dutta et al., 2021) . Furthermore, effective design of SARS-CoV-2 inhibitors can be accomplished through the selection of druggable targets. Main protease (M pro ) or 3chymotrypsin-like protease 3 (CLpro) cleaves non-structural proteins (Nsp) at eleven different sites and leads to the formation of Nsp4-Nsp16, which is essential for viral replication . The S protein and Mpro are important targets for designing SARS-CoV-2 inhibitors (Idroes et al., 2021) . No adverse toxic effect is likely to be observed by inhibiting Mpro as there are no known similar enzymes in humans with similar cleavage specificity (Khater and Nassar, 2021) . Papain-like protease (PL pro ) is another enzyme similar to M pro ; it generates Nsp1-Nsp3 by breaking at various sites, thereby generating essential proteins that suppress the immune system by speeding up the infection process in the host cell (Harcourt et al., 2004) . One of the important targets is RNA-dependent RNA polymerase (RdRp) as the viral genome depends on it for its replication and transcription (Borgio et al., 2020; Rudrapal et al., 2021) . In literature (Ahmad et al., 2021) , the possible potential of Indian medicinal plants and formulations against SARS-CoV-2/ COVID-19 has been mentioned based upon the recommendation (https://www.ayush.gov.in/docs/ayush-Protocol-covid-19.pdf; https://www.ayush.gov.in/docs/ayurved-guidlines.pdf) of AYUSH (Ayurveda, Yoga, Naturopathy, Unani, Siddha, Sowa-Rigpa and Homoeopathy).The usefulness of Triphala formulation and other Ayurvedic herbal remedies as therapeutic supplements in COVID-19 has also been reported in some recent literature (Ozah, 2020) . Khan et al. (2021) reported the inhibitory potential of Triphala against spike glycoprotein of SARS-CoV-2 by molecular docking. In view of these information, in the urge to design new SARS-CoV-2 inhibitors by blockage of its main protease (M pro ), we performed structure-based molecular docking studies of bioactive molecules from Triphala formulation. The selected candidates were subjected to MD simulations for confirmation of the SARS-CoV-2 inhibitory potential of bioactive molecules. Further, in silico screening of the selected compounds was conducted with the application of Lipinski's filters to determine their drug-likeness and ADMET properties.

The crystal structure of the Mpro of SARS-CoV-2 bound with non-covalent inhibitor X77 (N-(4tert-butylphenyl)-N-[(1R)-2-(cyclohexylamino)-2-oxo-1-(pyridin-3-yl) ethyl]-1H-imidazole-4carboxamide) at a resolution of 2.1 Å (PDB id: 6W63) was retrieved from RCSB Protein data bank (https://www.rcsb.org/). The maestro's protein preparation wizard (version 2021-2 Schrodinger suite) was used to prepare proteins by removing all water within 5Å distance from the ligand, filling amino acid side chains, adding missing hydrogen atoms, creating disulfide bonds, and optimizing hydrogen bonds. For restrained minimization, OPLS4 force field was used until the average RMSD (root-mean-square deviation) of the non-hydrogen atom converged to 0.30 Å, allowing heavy atoms to move sufficiently to release strained bonds, angles, and clashes (Roos et al., 2019; Kumar et al., 2021) .

The Triphala formulation's bioactive molecules were chosen from the literature (Prasad and Srivastava, 2020) , and their chemical structures were obtained from a popular online database, PubChem (https://pubchem.ncbi.nlm.nih.gov/). LigPrep was used to generate the ligands' threedimensional (3D) coordinates (Chen and Foloppe, 2010) . Ionization state (pH 7.4), tautomer, and chirality were predicted by Epik module. The OPLS4 force field was employed for energy minimization until an energetically stable confirmation was obtained (Roos et al., 2019) .

The selected set of phytoconstituents was subjected for molecular docking against Mpro with the Extra Precision (XP) protocol of the glide docking module of Maestro (Friesner et al., 2006) . The glide-receptor grid generation module was used with the default settings to identify the binding site, and the centroid of the bound ligand was taken into account for grid generation (Halgren et al., 2004) . The single best pose for each molecule was saved and the compounds were then classified based on their docking scores, which represent free binding energies.

Following that, the top four compounds with the best binding free energies, namely terflavin A, corilagin, chebulagic acid, and chebulinic acid, were selected to further ensure the flexibility and stability at the binding site of M pro . A 100 ns MD simulation was performed using GROnin-gen MAchine for Chemical Simulation (GROMACS) 2019.2 software (Abraham et al., 2015) . The topologies of the protein and ligand were determined using the 'pdb2gmx' script and the GlycoBioChem PRODRG2 server (http://davapc1.bioch.dundee.ac.uk/cgi-bin/prodrg), respectively. Energy minimization was accomplished using the GROMOS96 54a7 force field. The complex was solvated using the SCP water model, electrically neutralized by adding sodium ions and chloride ions using 'gmxgenion' script, equilibrated under an 0.3 ns isothermalisochoric (NVT) and 0.3 ns isothermal-isobaric (NPT) ensembles and subjected to MD simulations (100 ns) (Ghosh et al., 2021) . The pressure and temperature of the system were maintained at 1 atm and 300 K, respectively. The leap frog algorithm was used to integrate the motion equation in triplicate at a time step of 2 femtoseconds. Finally, the protein-ligand complexes were put to generate trajectories of RMSD, RMSF, Rg and H bond by GROMACS scripts. The trajectories were analyzed and visualized by VMD 1.9.3 (Humphrey et al., 1996) .

The binding free energies of protein-ligand complexes were quantified by the Molecular Mechanics Poisson Boltzmann Surface Area (MM-PBSA). This technique allows us to study the biomolecular interactions between the target and bioactive molecules. The g_mmpbsa script was utilized for MM-PBSA binding free energy calculations (Kumari et al., 2014) . Van der Waals, polar solvation, and electrostatic contributions to energy were calculated. Solvent accessible surface area (SASA) was used to calculate non-polar energy contribution.

The selected bioactive molecules were further checked for drug-likeness properties using Lipinski's rule of five (RO5) (Hussain et al., 2021; Rudrapal and Sowmya, 2019) . The ADMET properties including absorption, distribution, metabolism, excretion and toxicity of the selected were determined by the admetSAR server (Cheng et al., 2012; Ghost et al., 2021) . Oral in silico toxicity prediction was performed using web-based server ProTox-II (Banerjee et al., 2018) .

Triphala (Ayurvedic herbal formulation) is made up of dried fruits from three different herbal plants in equal parts (Jantrapirom et al., 2021; Tarasiuk et al., 2018) . Following a thorough review of the literature on the three plants T. chebula, T. bellirica, and P. emblica, 19 bioactive molecules were chosen for molecular docking against the M pro of SARS-CoV-2. The binding pocket around the inhibitor X77 of the target protein M pro was used as a reference for determining the grid box. Docking analysis generated many configurations, which were scored to determine the most favorable binding modes. The bioactive molecules docked against M pro of SARS-CoV-2 exhibited excellent binding, as evidenced by high docking scores that exceeded X77. The validation of the docking method was carried out by redocking of co-crystal ligand, X77 in the original ligand binding site of Mpro. The redocking was successful. The measured RMSD value of redocking was 0.4493 Å between the experimental docking pose of X77 and the natural binding mode of the native ligand. The superimposed conformation of the experimental ligand with the original co-crytal ligand, X77 is displayed in Fig. 2 . The docking results and the binding free energies of 19 compounds are given in Table 1. Terflavin A achieved the highest docking score, with a value of -14,521 kcal/mol. Corilagin, chebulagic acid, chebulinic acid and X77 had XP Gscore values of -12.845 kcal/mol, -11,359 kcal/mol, -11,217 kcal/mol, and -5,966 kcal/mol, respectively. Chemical structures of the four best-docked compounds and the native ligand X77 are depicted in Fig. 1 . Chebulagic acid has a lower XP Gscore value than chebulinic acid, but the binding energy value of chebulinic acid is better (-64.05 kcal/mol). Furthermore, despite having a higher XP Gscore, the dG bind value of the X77 ligand is -85,045 kcal/mol, lower than corilagin, chebulagic acid, and chebulinic acid.

The interactions with M pro of the top four docked compounds terflavin A, chebulagic acid, chebulinic acid, corilagin, including the native ligand X77 are visualized in Fig. 3. Terflavin A, chebulagic acid, chebulinic acid, corilagin are types of hydrolysable tannins, which are typically found on Terminalia species leaves (L.). Terflavin A and corilagin are ellagitannin, while chebulagic acid and chebulinic acid are benzopyran tannin. Tannins can inhibit SARS-CoV-2 through binding with catalytic dyad residues of M pro (Chang et al., 2019; Das et al., 2020) . Tannic acid even inhibits both M pro and TMPRSS2 protease, effectively suppressing SARS-CoV-2 (Wang et al., 2020) . The catalytic residues of 3CL pro (or M pro in the case of SARS-CoV-2) are Cys145 and His41, which are surrounded by other residues that confer substrate specificity Singh et al., 2020 (1.66 Å), and Glu166 (1.78 Å) then H-carbon bonds at Met165 (2.06 Å) and pi-alkyl and pisulphur interactions at residues Met49 (2.13 Å) and Cys145 (2.12 Å), respectively. As evidenced, H-bond interactions exist between corilagin and Cys145 (1.72 Å) and His41 (1.08 Å). As a native ligand, the X77 forms H-bonds at residues Asn142 (2.16 Å), Gly143 (1.89 Å), His163 (1.67 Å), and Glu166 (2.05 Å), then H-carbon bonds at Thr26 (2.10 Å) and Met165 (2.90) and pi-pi stacked and pi-sulphur interactions at residues His41 (2.17 Å), Met49 (2.04 Å), and Cys145 (1.14 Å), respectively. These findings support previous findings that the selected compounds bind to the catalytic residues of M pro (Silva et al., 2021) . The substrate-binding cleft contains a catalytic dyad composed of His41 and Cys145, which is located between domain I and II at the N-terminus, and a hydrogen-bonded water molecule bound to His41 acts as a third component of the catalytic triad (Tan et al., 2005) .

MD simulations and other related techniques are on the verge to act as routine techniques for the discovery of drugs (De Vivo et al., 2016) . In this study, the stability and dynamic behavior of the four ligand-protein (bound) complexes were verified by MD simulations for a 100 ns time. The RMSD, root mean square fluctuation (RMSF), radius of gyration (Rg) and solvent accessible surface area (SASA) graphs were plotted of the apoprotein, co-crystal ligand complex and four ligand-protein complexes. The RMSD analysis provides information insight into the structural deviation and conformational stability and convergence during the time of the simulation. Thereby, in an equilibrated, stable system, the RMSD towards the end of MD simulations should be within 1-3 Å for a globular protein. The time-dependent RMSD graph was plotted for the backbone of the apoprotein, co-crystal ligand-protein complex and ligand (terflavin A, corilagin, chebulagic acid, chebulinic acid) protein (M pro ) complex. Fig. 4 (A) shows the RMSD values of apoprotein and co-crystal ligand-protein complex fluctuate in the beginning from 0.2 nm and at the end at around ~0.3 nm, confirming the stability of the protein during 100 ns of simulation. At first, the RMSD values of the apoprotein and the protein-ligand complex are different, but after 50 nm, the RMSD values become stable and remain stable until 100 ns. Fig. 4 (B) shows the RMSD values of protein (M pro ) and ligand (terflavin A, corilagin, chebulagic acid, chebulinic acid) complex, indicated with colour code including M pro -TerA (green), M pro -Cor (pink), M pro -Cgic (cyan), and M pro -Cnic (indigo). The RMSD of M pro -TerA fluctuates at the beginning until 70 ns (0.2 nm to ~0.6 nm) and at 100 ns, it fluctuates at around ~0.3 nm. The RMSD values of both the Protein-Corilagin (M pro -Cor) and Protein-Chebulagic acid (M pro -Cgic) complex are stable from the beginning (0.2 nm), and stay stable until 100 nm (0.25 nm), confirming the stability of both protein-ligand complexes during 100 ns of simulation. Whereas, Protein-Chebulinic acid complex fluctuats (between 0.2 nm to 0.3 nm) up to 70 ns, but after 70 nm, it becomes stable and stays stable (0.2 nm) until 100 nm. The overall RMSD results reveal that the binding of the ligands (corilagin, chebulagic acid, chebulinic acid) in the active site of M pro is stable. The RMSF study gives data regarding the amino acid residues fluctuation from the reference structure. Fig. 4 (C) shows that the co-crystal ligand fluctuation is almost similar to apoprotein, whereas Fig. 4 (D) shows the fluctuation of the ligand's binding to active site residues Thr24, Thr25, His41, Cys44, Met49, Tyr54, Phe140, Asn142, Gly143, Cys145, His163, His164, Met165, Glu166, Leu167, Pro168, Asp187, Arg188, Gln189, and Thr190 in each of M pro TerA, M pro -Cor, M pro -Cgic, and M pro -Cnic indicates that there is nearly similar RMSF fluctuation compared to apoprotein and co-crystal ligand-protein residues. Thus, analysis of the RMSF suggested a stable protein-ligand complex since there are no significant changes at the active site's reference location after successful binding by the apoprotein. Further, the 'sasa' function was used to calculate the SASA to examine the accessibility of solvent molecules on the protein surface. The SASA value of apoprotein, co-crystal ligandprotein complex, and protein-ligand complexes M pro -TerA, M pro -Cor, M pro -Cgic, and Mpro-Cnic were stable during 100 ns MD. Thus, it indicates that ligand-protein complexes have comparable stability to apoprotein and co-crystal ligand-protein complexes [ Fig. 4 (E) ]. During MD simulations, a study on the radius of gyration (Rg) indicates the protein compactness. Rg of the C α atoms of the Mpro SARS-CoV-2 protein was computed to ascertain conformational changes in it. The Rg graph was plotted for protein using 100 ns trajectories as represented in Fig. 4 (F) . The Rg plot indicates that the Rg value apoprotein (mean: 2.20 nm) and protein-ligand complex (mean: 1.23 nm) are nearly identical. Finally, Rg results suggest that all secondary structure of protein complexes is closely packed during 100 ns MD simulation. Previous literature has highlighted the potential impact of medicinal plants and their phytoconstituents against SARS-CoV-2 . Studies have reported that thioflavonoids and epicatechin-3-O-gallate interacted with (binding sites) Mpro and PLpro, and were considered as SARS-CoV-2 inhibitors Ismail et al., 2021) . Thus, bioactive molecules bind to or interact with the active residues which are present at the binding cavity of Mpro and PLpro.

Python script MmPbSaStat.py provided in the g_mmpbsa was utilized to calculate average BFE and its standard deviation error. The energy released during the interaction of protein-ligand is represented in the form of BFE. It is a cumulative sum of Van der Waals, polar salvation, electrostatic, and SASA energy. Precisely, the MM-PBSA method estimates the binding free energy of a ligand-protein complex as the difference between the free energy of the complex and the free energies of the unbound components (receptor and ligand). More negative the BFE, the more effectively the ligand binds to the protein ( Table 2 ). The results indicate that all the forms of energy contributed fairly to the interactions between M pro and the selected molecules. M pro -TerA complex showed the least binding BFE (-233.782 kJ/mol) among the selected bioactive molecules followed by M pro -Cgic (-221.916 kJ/mol). The binding free energy (BFE) of the X77 (-160.625 kJ/mol) is less negative compared to the selected bioactive molecules. A comparison of the BFE of all the complexes reveals that, except for corilagin, all of the other selected bioactive molecules outperform X77 in inhibiting the protein M pro of SARS-CoV-2.

Lipinski's rule of five (Ro5) is a set of rules used to assess the drug-likeness of a chemical compound that has a pharmacological activity that could lead to it being an orally active drug for human use (Lipinski et al., 2001; Kalita et al., 2020) . The Ro5 specifies that the molecular weight must be less than 500 Da, the number of H-bond donors must be less than 5, the number of H-bond acceptor must be less than 10, the xlogP must be less than 5, and the molar refractivity must be between 40 and 130 (Fatimawali et al., 2021; Kashyap et al., 2016; Sefren Geiner et al., 2020) . Regarding this statement, Table 3 indicates that only compound X77 appears to meet the Ro5 criteria. However, strict adherence to Ro5 is not required for naturally occurring compounds or phytocomponents. Lipinski specifically states that the natural products and substrates of biological transporters are not covered by Ro5 (Lipinski et al., 2001; Rudrapal et al., 2017; Rudrapal et al., 2018) . It is, therefore important to note here that the desired pharmacological efficacy of Ayurvedic/ herbal formulations is possibility due to the multi-component nature (different components) of such formulations.

On its way to the target site, the drug candidate undergoes pharmacokinetic events such as absorption, distribution, metabolism, and excretion (ADME). ADME parameters, in conjunction with toxicity, are critical for the discovery and development of potential drug candidates (Guan et al., 2019) . Not only should a high-quality drug candidate be efficacious against the therapeutic target, but it should also exhibit appropriate ADMET properties at the dosage for therapeutic purposes. Based on the results of the admetSAR and ProTox analyses (Table 4 ), all of the compounds were from any signs of toxicity such as hepatotoxic, carcinogenic, or mutagenic. The acute oral toxicity analysis reveals that all compounds are classified as class III (500 mg/kg LD 50 5000 mg/kg), which means that they are not toxic, except for chebulagic acid, which is classified as category II (50 mg/kg LD 50 500 mg/kg), indicating a moderate level of toxicity (Herkenne et al., 2008; Li et al., 2014) . Typically, animal experiments are used to determine a compound's toxicity, carcinogenicity, and mutagenicity. However, this takes a long time and necessitates the sacrifice of animals. Thus, prediction of toxicity within in silico technique is a faster and less expensive alternative, as it is based on existing data on a compound's toxicity. As a result, it is used to develop models capable of predicting the toxicity, carcinogenicity, and mutagenicity of novel compounds.

Bioavailability is a pharmacological term that refers to the level and amount to which a bioactive compound can be absorbed and made available (Kim et al., 2014) . A compound's low oral bioavailability is a frequent issue during the drug design process. As a result, a drug may fail in clinical trials despite its high efficacy in vitro and in vivo experiments (Fukunishi et al., 2014) . A compound is predicted to have low oral bioavailability if its bioavailability score is lower than 0.5. If the score is greater than 0.5, the compound is expected to have a high oral bioavailability (Rath et al., 2021) . According to the findings, all compounds, except for X77, have a high oral bioavailability. Toxic doses are frequently expressed as LD 50 (mg/kg body weight). LD50 (median lethal dose) is a particular dose wherein 50% of the participants die after being exposed to a compound (Erhirhie et al., 2018) . A compound is considered highly toxic if the LD 50 value is higher than 300 mg/kg. If the LD50 is between 300 and 1,000 mg/kg, a compound is regarded as moderately toxic. If, on the other hand, the LD 50 value of a compound is between 1,000 and 5,000 mg/kg, it is considered to be less or slightly toxic. Thus, chebulagic acid, chebulinic acid, and X77 are classified as moderately toxic, whereas corilagin and terflavin A are classified as less toxic.

To develop drugs from ancient Indian knowledge against SARS-CoV-2, we made an attempt to analyze the binding modes of the phytoconstituents of Triphala formulation by docking against the active site of M pro . Terflavin A, corilagin, chebulagic acid and chebulinic acid demonstrated promising binding affinity with the target protein, M pro . Molecular dynamics simulation was carried out for these selected bioactive molecules to further elucidate their binding affinity towards the M pro and explored the obtained trajectories for RMSD, RMSF, Rg, and SASA. The trajectories were found to have increased flexibility and higher stability. Further, strong binding affinity between the bioactive molecules and M pro were confirmed by low binding free energies, dynamic behaviors of complexes as well as residue contribution energies obtained by MM-PBSA. This study identifies four bioactive molecules from Triphala formulation as promising SARS-CoV-2 M pro inhibitors; hence these bioactive molecules can be further examined using experimental studies for confirmation of their therapeutic efficacies.

This research received no external funding.

The authors declare no conflicting interests.

Data and material will be made available after proper request to the corresponding author.

MR developed the idea of the work. IC prepared the research methodology. IC, TET, MR, RY, TS, JHZ, DC and SGW were involved in software-based investigations, data analysis and interpretation of results. RY and TS prepared the original draft. MR, IC, PKP, SJK and DC reviewed the draft. JK, MAA, RMIIA, ARB and RKS were involved in funding acquisition. JHZ and MR completed the manuscript writing and editing. All authors read and approve the manuscript. 

GROMACS: High performance molecular simulations through multi-level parallelism from laptops to supercomputers

Indian medicinal plants and formulations and their potential against COVID-19-preclinical and clinical research

ProTox-II: a webserver for the prediction of toxicity of chemicals

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

State-of-the-art tools unveil potent drug targets amongst clinically approved drugs to inhibit helicase in SARS-CoV-2

Bromelain a Potential Bioactive Compound: A Comprehensive Overview from a Pharmacological Perspective

A Comprehensive Review of the Structure Elucidation of Tannins from Terminalia Linn. Evidence-Based Complement

Drug-like Bioactive Structures and Conformational Coverage with the LigPrep/ConfGen Suite: Comparison to Programs MOE and Catalyst

Plants of the Genus Terminalia: An Insight on Its Biological Potentials, Pre-Clinical and Clinical Studies

Role of Molecular Dynamics and Related Methods in Drug Discovery

Appraisals of the Bangladeshi Medicinal Plant Calotropis gigantea Used by Folk Medicine Practitioners in the Management of COVID-19: A

Advances in acute toxicity testing: strengths, weaknesses and regulatory acceptance

Data on the docking of phytoconstituents of betel plant and matcha green tea on SARS-CoV-2. Data Br

Extra Precision Glide: Docking and Scoring Incorporating a Model of Hydrophobic Enclosure for Protein−Ligand Complexes

Prediction of Synthetic Accessibility Based on Commercially Available Compound Databases

Design, Molecular Docking, Druglikeness and Molecular Dynamics Studies of 1,2,4-Trioxane Derivatives as Novel P. falciparum Falcipain-2 (FP-2) Inhibitors. Biotechnologia

Identification of potential plant-based inhibitor against viral proteases of SARS-CoV-2 through molecular docking, MM-PBSA binding energy calculations and molecular dynamics simulation

ADMET-score -a comprehensive scoring function for evaluation of chemical drug-likeness

Studies on the chemical constituents of Phyllanthus emblica

Glide: A New Approach for Rapid, Accurate Docking and Scoring. 2. Enrichment Factors in Database Screening

Identification of Severe Acute Respiratory Syndrome Coronavirus Replicase Products and Characterization of Papain-Like Protease Activity

In Vivo Methods for the Assessment of Topical Drug Bioavailability

VMD: Visual molecular dynamics

Bioactive Antidiabetic Flavonoids from the Stem Bark of Cordia dichotoma Forst.: Identification, Docking and ADMET Studies

Suhendrayatna, 2021. The study of Calotropis Gigantea leaf metabolites from Ie Brouk geothermal area Lamteuba-Aceh Besar using molecular docking

Targeting multiple conformations of SARS-CoV2 Papain-Like Protease for drug repositioning: An in-silico study

Pharmacological Benefits of Triphala: A Perspective for Allergic Rhinitis

Design, Synthesis, Antimalarial Activity and Docking Study of 7-Chloro-4-(2-(substituted benzylidene)hydrazineyl)quinolines

Synthesis, Antimalarial Activity Evaluation and Drug likeness Study of Some New Quinoline-Lawsone Hybrids

Discovery of Potential Inhibitors of the Receptor-binding Domain (RBD) of Pandemic Disease-causing SARS-CoV-2 Spike Glycoprotein from Triphala Through Molecular Docking

In silico molecular docking analysis for repurposing approved antiviral drugs against SARS-CoV-2 main protease

Critical Evaluation of Human Oral Bioavailability for Pharmaceutical Drugs by Using Various Cheminformatics Approaches

Phytochemistry and medicinal potential of the terminalia bellirica roxb. (bahera)

Cerebroprotective Effect of Aloe Emodin: In Silico and In Vivo Studies

g_mmpbsa -A GROMACS Tool for High-Throughput MM-PBSA Calculations

In silico Prediction of Chemical Acute Oral Toxicity Using Multi-Classification Methods

Experimental and computational approaches to estimate solubility and permeability in drug discovery and development settings 1PII of original article: S0169-409X(96)00423-1. The article was originally published in Advanced

Triphala: a useful therapeutic supplement during COVID-19 pandemic

Therapeutic Uses of Triphala in Ayurvedic Medicine

Oxidative Stress and Cancer: Chemopreventive and Therapeutic Role of Triphala

Design, synthesis and molecular modeling studies of novel mesalamine linked coumarin for treatment of inflammatory bowel disease

OPLS3e: Extending Force Field Coverage for Drug-Like Small Molecules

Novel series of 1,2,4-trioxane derivatives as antimalarial agents

In silico screening of phytopolyphenolics for the identification of bioactive compounds as novel protease inhibitors effective against SARS-CoV-2

Coronavirus disease-2019 (COVID-19): An updated review

Design, Synthesis, Drug-Likeness Studies and Bio-Evaluation of Some New Chalconeimines

Newer series of trioxane derivatives as potent antimalarial agents

The potential of leaf extract of Pangium edule Reinw as HIV-1 protease inhibitor: A computational biology approach

An in-silico evaluation of different bioactive molecules of tea for their inhibition potency against non structural protein-15 of SARS-CoV-2

Drug repurposing and computational modeling for discovery of inhibitors of the main protease (M pro ) of SARS-CoV-2

In-silico evaluation of bioactive compounds from tea as potential SARS-CoV-2 nonstructural protein 16 inhibitors

A comprehensive review on promising anti-viral therapeutic candidates identified against main protease from SARS-CoV-2 through various computational methods

pH-dependent Conformational Flexibility of the SARS-CoV Main Proteinase (Mpro) Dimer: Molecular Dynamics Simulations and Multiple X-ray Structure Analyses

Triphala: current applications and new perspectives on the treatment of functional gastrointestinal disorders

Anti-SARS-CoV Natural Products With the Potential to Inhibit SARS-CoV-2 (COVID-19)

Tannic acid suppresses SARS-CoV-2 as a dual inhibitor of the viral main protease and the cellular TMPRSS2 protease

Design of Wide-Spectrum Inhibitors Targeting Coronavirus Main Proteases

Chebulagic acid -221

All MD simulations presented in the study were carried out with resources provided by TÜBTAK (Turkish Scientific and Technological Research Council), ULAKBIM (Turkish Academic Network and Information Centre), and TRUBA (High Performance and Grid Computing Centre).