key: cord-0789003-u2h5c916
authors: Swargiary, Ananta; Mahmud, Shafi; Saleh, Md. Abu
title: Screening of phytochemicals as potent inhibitor of 3-chymotrypsin and papain-like proteases of SARS-CoV2: an in silico approach to combat COVID-19
date: 2020-10-22
journal: Journal of biomolecular structure & dynamics
DOI: 10.1080/07391102.2020.1835729
sha: fe8cce48245c02a77ff0f7e9498ea4d96e19dd6d
doc_id: 789003
cord_uid: u2h5c916

COVID-19 and its causative organism SARS-CoV2 that emerged from Wuhan city, China have paralyzed the world. With no clinically approved drugs, the global health system is struggling to find an effective treatment measure. At this crucial juncture, screening of plant-derived compounds may be an effective strategy to combat COVID-19. The present study investigated the binding affinity of phytocompounds with 3-Chymotrypsin-like (3CLpro) and Papain-like proteases (PLpro) of SARS-CoV2 using in-silico techniques. A total of 32 anti-protease phytocompounds were investigated for the binding affinity to the proteins. Docking was performed in Autodock Vina. Pharmacophore descriptors of best ligands were studied using LigandScout. Molecular dynamics (MD) simulation of apo-protein and ligand-bound complexes was carried out in YASARA software. The druglikeness properties of phytocompounds were studied using ADMETlab. Out of 32 phytochemicals, amentoflavone and gallocatechin gallate showed the best binding affinity to 3CLpro (–9.4 kcal/mol) and PLpro (–8.8 kcal/mol). Phytochemicals such as savinin, theaflavin-3,3-digallate, and kazinol-A also showed strong affinity. MD simulation revealed ligand-induced conformational changes in the protein with decreased surface area and higher stability. The RMSD/F of proteins and ligands showed stability of the protein suggesting the effective binding of the ligand in both the proteins. Both amentoflavone and gallocatechin gallate possess promising druglikeness property. The present study thus suggests that Amentoflavone and Gallocatechin gallate may be potential inhibitors of 3CLpro and PLpro proteins and effective drug candidates for SARS-CoV2. However, the findings of in silico study need to be supported by in vivo studies to establish the exact mode of action. Communicated by Ramaswamy H. Sarma

COVID-19 pandemic is a rapidly spreading dreadful disease caused by a novel coronavirus SARS-CoV2 with more than 24.02 million positive cases across the world till date (2:44 p.m. CEST, 27 August 2020, WHO) . According to the latest WHO update, the global death rate of COVID-19 now stands at 3.42% (https://covid19.who.int/). SARS-CoV2 is a positivesense single-stranded RNA virus belonging to the family Coronaviridae that originated from the Wuhan city of China in 2019 (Tang et al., 2020) . It is considered to be the seventh human coronavirus (HCoV) under the same family after 229E, NL63, OC43, HKU1, SARS-CoV, and MERS-CoV (Benvenuto et al., 2020) . Like SARS and MERS coronaviruses, SARS-CoV2 is a highly pathogenic and contagious virus that causes severe respiratory complications and even death (Cui et al., 2019) . Sequence analysis reveals that all the human coronaviruses have a common origin (from bats?) except HKU1 and OC43 which are believed to be originated from rodents (Forni et al., 2017) . The attack and pathogenicity of SARS-CoV2 progress with four transmission stages-asymptomatic, moderate, extreme, and clinical which starts at the lower respiratory tract followed by invasion of pulmonary epithelial cells and hijacking the entire host cell machinery. The most common symptoms of COVID-19 include cough, fever, malaise, gastrointestinal symptoms, loss of smell, sore throat, heart failure, and acute kidney injury (Gandhi et al., 2020) .

Currently, there is no clinically approved drug/vaccine for the treatment of COVID-19 and SARS-CoV2. However, several clinical trials are in progress to find out possible therapies for COVID-19. Recently, hydroxychloroquine, remdesivir, and ivermectin have been positioned as a possible treatment for COVID-19 (Caly et al., 2020; Wang et al., 2020) . Several researchers are looking at the repurposing of existing antiviral drugs to test the efficacy on SARS-CoV2 (Elfiky, 2020; Kandeel & Al-Nazawi, 2020) . Meanwhile, it is also crucial to understand the pathogenic mechanism of SARS-CoV2 to develop effective drugs or vaccines against the virus. Based on the virus pathogenesis and clinical manifestations, several studies have put forwarded potential therapeutic drug targets of SARS-CoV2 (Omolo et al., 2020) . Among the several drug targets, coronavirus 3CLpro and PLpro are considered to be crucial targets of COVID-19 management (Zhou et al., 2020) . 3CLpro and PLpro are encoded by ORF-1 of virus genome which plays a major role in the processing of inactive polyproteins (pp1a and pp1ab) into active non-structural proteins (nsps). These nsps finally helps in the replication, generation, as well as infection of the virus (Zhou et al., 2019) . Inhibition of these proteases can be crucial for further replication and infection of SARS-CoV2 . Recent studies revealed that 3CLpro and PLpro share 96% and 83% sequence identity (at protein level) between SARS-CoV and -CoV2 (Morse et al., 2020) . However, no such distinct structural differences were found in the active-site pockets of the enzymes which therefore suggest that 3CLpro and PLpro could be a potential drug target in SARS-CoV2. Therefore, drugs or compounds that are used as protease inhibitors in SARS-CoV may also be investigated for its effectiveness against SARS-CoV2.

Plants have always been a rich source of bioactive compounds having tremendous medicinal values. Unlike synthetic drugs that are costly, and possess undesirable side-effects, phytocompounds are less expensive, readily available, and contain less side-effects. Several medicinal plants and isolated compounds have been established as anti-virus and protease inhibitors by many researchers Park et al., 2013) . reported the SARS-CoV protease inhibitor property of apigeni and luteolin isolated from Tripterygium regelii. Similarly, major bioactive compound rhamnazin 3-Orutinoside isolated from Sarcocornia fruticosa strongly inhibited the protease activity of Hepatitis-C virus (Hawas et al., 2019) . Several phytochemicals of known and unknown biological functions have been investigated by many authors in a computer-based approach to find out lead compounds that can effectively inhibit therapeutic targets in SARS-CoV2 including 3CLpro and PLpro (Chojnacka et al., 2020; Vardhan & Sahoo, 2020) . The present study investigated the binding affinity of anti-protease phytochemicals with 3-chymotrypsin-and papain-like proteases of SARS-CoV2 using molecular docking and molecular dynamics simulation approach. Furthermore, we also investigated the in silico druglikeness and ADMET properties of the phytochemicals.

A total of 32 phytocompounds that are reported to have protease inhibitor property are collected from different pieces of literature. The name of the chemicals and the publication details are presented in Table 1 . The sdf files of each chemical were downloaded from the Pubchem database (http://pubchem.ncbi.nlm.nih.gov/). The files were converted into pdb file using OpenBabel software (O'Boyle et al., 2011) . Finally, chemicals were converted into pdbqt file using Autodock tool (Trott & Olson, 2010 ).

Three-dimensional structures of the 3-Chymotrypsin-like protease (3CLpro) and Papain-like protease (PLpro) of SARS-CoV2 were downloaded from RCSB-PDB database (https:// www.rcsb.org/). The pdb IDs of the crystal structures of proteases were 6M2N for 3CLpro and 7JN2 for PLpro. Both the protein structures were processed by removing the water and other hetatms. Next, polar hydrogens and Kollman charges were added to the pdb files, and finally converted into pdbqt file format using AutoDock Tools (Trott & Olson, 2010 ).

After the compounds and enzymes were prepared docking was carried out using AutoDock Vina (Trott & Olson, 2010) . The amino acid residues interacting with the co-crystallized ligands of both the 6M2N and 7JN2 were taken as the active site residues and accordingly docking grid parameters were set. The active pocket amino acids residues in 6M2N were His41, Cys44, Met49, Leu141, Asn142, Gly143, Ser144, Cys145, His164, Met165, Glu166, Asp187, and Gln189 while the active pocket of 7JN2 contains amino acid residues such as Leu162, Gly163, Asp164, Glu167, Pro247, Pro248, Tyr264, Tyr268, Gln269, Tyr273, and Thr301. The grid parameters were set as x, y, z size-and centre-coordinates: 40, 50, 72 and -32.939, -65.379, and 39 .480 for 3CLpro and for PLpro the parameters were 48, 54, 42, and 51.130, 31.061, and 1.325 for x, y, z sizeand centre-coordinate, respectively. The docking algorithm was carried out by keeping the default exhaustiveness at 8. After docking, the best poses scoring the lowest binding energy (kcal/mol) for both 3CLpro and PLpro were selected and visualize in Discovery Studio.

To validate the docking study and evaluate the change in protein conformation, molecular dynamics simulation was implemented. The YASARA dynamics commercial package (Krieger et al., 2003) was used to simulate the complex where AMBER14 force field was applied (Dickson et al., 2014) . The water molecules were added and with 0.9% NaCl salt at 310 K temperature (Krieger et al., 2006) . The periodic boundary condition was maintained where the Particle Mesh Ewald method was applied to calculate long range electrostatic interaction. A cubic simulation cell was created to simulate the protein which was 20 Å bigger than the protein. The Berendsen thermostat was used to control the temperature of the simulation system. The protein complex and apo structure was initially cleaned and minimized by steepest gradient approaches (5000 cycles). The normal simulation time was maintained with a time step of 1.25 fs (Krieger & Vriend, 2015) . The simulation trajectory was saved for every 100ps to analyze the result. The simulation result was incorporated with the default script of YASARA. Finally, root mean square deviation (RMSD), root mean square fluctuation (RMSF), radius of gyration (Rg), solvent accessible surface area (SASA), and Hydrogen bonds (H-bonds) were analyzed (Islam et al., 2019; Khan et al., 2020; Mahmud et al., 2020a Mahmud et al., , 2020b . To have more correct result from molecular dynamics simulation, each complex were run thrice (n¼3) and average result was used for analysis. To calculate the binding free energy, the Molecular Mechanics Poisson-Boltzmann Surface Area method was applied. The edited YASARA script was employed for this calculation. The following equation had been used for the calculation of MM-PBSA (Gilson & Honig, 1988) .

Here, DG MM is the sum of van der Waals and electrostatic interaction, DG PB and DG SA is the polar and non-polar solvation energies, TDS is the entropic contribution (Razzaghi-Asl et al., 2018) . The 1000 trajectory files were considered for MM-PBSA calculation.

Pharmacophore features of the two best binding compounds were determined by LigandScout software which demonstrated Structural Activity Relationship (Muchtaridi et al., 2017) and also suggested compound's descriptors necessary for optimal molecular interactions with the biological targets.

The optimized chemical structures (pdb) of all the compounds were loaded into Ligandscout working space and key pharmacophore features were determined including Hbond donor, H-bond acceptor, hydrophobic, aromatic, halogen bond donor, positively and negatively ionizable groups.

Top ten best binding phytochemicals with the 3CLpro and PLpro were selected and druglikeness properties were studied using SwissADME (Daina et al., 2017) and PubChem database. The druglikeness properties of selected drugs were evaluated based on Lipinski's rule (Lipinski, 2004) . Similarly, the absorption, distribution, metabolism, excretion, and toxicity (ADMET) properties of drugs were predicted using ADMETlab (Dong et al., 2018 ).

Statistical calculation were carried out in Microsoft excel and OriginPro softwares. Values are expressed in mean ± standard deviation. 

The search for drugs and vaccines for COVID-19 pandemic is the most urgent and challenging task for global researchers at this moment. In this drive of drug discovery, we have investigated 32 phytochemicals from 17 medicinal plants that are reported to have anti-protease property. Table 1 showed the name of the plants, isolated compounds, and the target proteins of the phytocompounds reported from different plants. All the selected phytocompounds are reported to possess in vitro protease inhibition properties. The compounds were found to have inhibitory activity against 3CLpro and PLpro enzymes of many viruses such as SARS-CoV, HIV, and Hepatitis C virus.

Three-dimensional structures of 3CLpro and PLpro have been retrieved from PDB database. Structurally, 3CLpro and PLpro contains 306 and 312 amino acid residues, respectively. The total numbers of atoms present in the protein structures were 4682 and 4898, respectively. The secondary structure of 3CLpro has 17.26% alpha helix and 25.01% beta strands and the protein secondary structure elements (SSE) were found to be 42.27%. Similarly, PLpro has 21.9% and 28.23% alpha helix and beta strand, respectively forming the total SSE 50.13%. The active site of the 3CLpro protein is found to be located at the cleft region between domain I and II region. Similarly, the active site of PLpro is found to be located at the groove region between the Thumb domain and Finger domain of the protein.

All the 32 phytocompounds were analyzed for the binding affinity with two key proteases of SARS-CoV2. Figure 1 showed the binding energies (kcal/mol) of all the 32 compounds with 3CLpro and PLpro proteins of the SARS-CoV2. In silico docking studies found that the binding affinities of anti-protease phytocompounds range from -5.1 to -9.4 kcal/ mol in 3CLpro and -5.5 to -8.8 kcal/mol in PLpro proteins, respectively. The average binding energies were found to be -7.34 kcal/mol and -7.48 kcal/mol for 3CLpro and PLpro proteins, respectively. Half of the compounds were found to have -7.3 kcal/mol or higher binding affinity against 3CLpro while 18 compounds out of 32 were found to have stronger binding affinity than the average energy. Amentoflavone showed the strongest affinity to 3CLpro while cinnamic amide showed weakest binding energy. Similarly, gallocatechin gallate and sinigrin showed the strongest and weakest binding affinities to PLpro protein. The two-dimensional display of binding interactions, ligand-protein interacting residues, hydrophobicity as well as Ramachandran plots are presented in Figure 2 . The threedimensional surface view structure with ligand binding sphere view shows that ligand amentoflavone binding site is in between the domain-I and -II region of the 3CLpro protein ( Figure 2(a) ). Thirteen amino acid residues of 3CLpro such as His41, Cys44, Met49, Phe140, Leu141, Asn142, Cys145, Met165, Glu166, Val186, Asp187, Arg188, and Gln189 were found to make interactions with the ligand. Out of 13 active site amino acid residues, 10 residues were found to make interaction with the ligand. Five H-bonds were formed between ligand and His41, Met49, Glu166, Arg188, and Gln189 amino acids of 3CLpro (Figure 2(b) ). Figure 2 (c) showed the amino acid distribution in the Ramachandran plot of complete 3Clpro protein while Figure 2 (f) compares distribution of the amino acid residues that are interacting to the ligand only. Most of the amino acid residues (seven) that are interacting with the ligand amentoflavone were from coil-structure of 3CLpro while two amino acids each from helix, turn, and b-sheet structure that made the interaction with the ligand. Out of 13 ligand-interacting residues, six are strongly hydrophilic and remaining seven are hydrophobic in nature. Figure 3 showed the two-dimensional display of binding interactions, ligand-protein interactions, hydrophobicity, and Ramachandran plots of the ligand-interacting residues. The three-dimensional surface view structure with ligand binding sphere view showed that ligand gallocatechin gallate binding site lies between thumb domain and palm domain region of PLpro protein (Figure 3(a) ). Fourteen amino acid residues of PLpro Asp164, Arg166, Glu167, Met208, Pro247, Pro248, Tyr264, Gly266, Asn267, Tyr268, Gln269, Cys270, Tyr273, and Thr301 were found to make interactions with the ligand. Out of 11 active site amino acid residues, nine residues were found to make interactions with the ligand. Five amino acids formed conventional H-bond with the ligand while other residues formed van der Waals interaction (11 residues), pibonds, and other non-covalent interactions (Figure 3(b) ). Ramachandran plots displayed the amino acid distribution and conformation of amino acid residues in the plot ( Figure  3(c,f) ). Nine amino acids residues of PLpro are found to be distributed outside the allowed region of the plot. Out of 14 amino acid residues interacting with the ligand, 6 belongs to b-sheet structure, 4 to coil structure, and 2 each from turn and helix structure of the protein. Amino acid Gln269 formed unfavorable bond with the ligand and is place outside the allowed region in the Ramachandran plot (Figure 3(b,f) ). The H-bond donor or acceptor residues as well as hydrophobicity property of the ligand-interacting amino acid residues are presented in Figure 3 (d,e). The ligand surrounding amino acids are found to be entirely consist of hydrophilic residues except two amino acids Met209 and Cys270 that are hydrophobic in nature (Figure 4 ).

The RMSD of 3CLpro apo-protein (ligand-free) and ligandbound complexes, Rg fluctuations, SARA values, and hydrogen bond profile during the period of simulation (100 ns) are presented in Figure 5 . Figure 5(a) showed that the apo-protein and ligand bound 3CLpro and PLpro exhibit initial rise till 5 ns and thereafter showed stable nature. But after 40 ns all four complexes showed unstable due to the fluctuations in RMSD profile. Interestingly, the RMSD descriptors of the complexes never exceeded 2.5 Å value, which denotes the structural integrity of the proteins. The average RMSD descriptors of 3CLpro apo-protein and ligand-complex were found to be 1.13 ± 0.9 Å and 1.191 ± 0.2 Å, respectively. Similarly, the RMSD of apo and ligand-bound PLpro were found to be 1.11 ± 0.01 Å and 1.06 ± 0.02 Å, respectively. The MD simulation result from the 3CLpro and PLpro complex suggested that upon ligand binding no significant deviations or conformational changes were taken place in the protein structure. The radius of gyration of the protein complex denotes the degree of compactness and rigidness of the protein. Greater Rg value indicates higher flexibility and conformation of the protein whereas lower Rg value denotes more rigidity. Figure 5 (b) showed that the 3CLpro protein has higher stability compared to PLpro. It has been observed that both the proteases have almost similar stability in apoform as well as ligand bound-complex form. The average Rg values of 3CLpro were found to be 22.37 ± 0.90 Å and 22.36 ± 0.90 Å for apo and ligand bound complex. Similarly, Rg values were found to be 23.44 ± 0.30 Å and 23.25 ± 1.30 Å for apo-and ligand-bound PLpro protein, respectively. Figure  5 (c) showed that the solvent-accessible surface areas (SARA) of 3CLpro and PLpro both in apo form and ligand-bound complex form. In both the proteases, ligand-bound complexes were found to have smaller SASA value compared to apo-form of the proteins. The average SASA values decreased from 14345.04 ± 0.20 Å 2 in apo-form to 14168.25 ± 0.90 Å 2 in ligand-bound complex of 3CLpro. However, after 60 ns both the apo and ligand-bound 3CLpro complex showed almost similar surface area. Like Rg, PLpro was found to have slightly bigger solvent accessible surface area than 3CLpro. The SASA values were found to be 15033.04±.02 Å 2 and 14804.37 ± 0.9 Å 2 for apo and ligand- bound structures of PLpro, respectively. Moreover, to understand the hydrogen bond and their contributions in gaining stability of the four systems, quantitative assessment of hydrogen bond was done. From Figure 5 (d) it can be observed that, all four complexes did not deviate and similar number of hydrogen bond were found for apo-protein and ligand-bound complexes. The average numbers of hydrogen bonds were found to be 534.05 ± 0.1, 526.04 ± 0.9, 575.14 ± 0.30, and 573.77 ± 1.10 for 3CLpro apo, 3CLpro-Ligand, PLpro apo, and PLpro-ligand complexes, respectively. Upon binding of ligands to the proteins there were slight changes in the number of H-bonds. The fluctuations of amino acids of both the protease enzymes in apo-and ligand-bound forms, fluctuations of ligand-atoms, and binding energies during the period of simulations (100 ns) were presented in Figure 6 . The root mean square fluctuation (RMSF) is useful for characterizing local changes along the protein chain. The amino acid residues Ser1 (gamma-turn), Glu55 (helix-strand), Ile59 (helix-strand), Gln244 (helix-strand), and Gln306 (gamma-turn) from SARS-Cov-2 3CLpro exhibited more fluctuations. While in case of PLpro protein, higher RMSF peaks were observed in His47 (beta-turn), Asn48 (beta-turn), Val66 (helix-strand), Leu113 (helix-strand), Lys190 (beta-turn), Lys218 (helix-strand), Lys228 (beta-turn), Gln229 (helix-strand), Cys270 (beta-turn), and Lys315 (helix-strand) amino acid residues. Overall, less flexible nature of residues stabilized both the protein-ligand complexes. Figure 6 (c) showed the RMSDs both the ligands of 3CLpro and PLpro proteases. Ligand RMSD indicates how stable the ligand is with respect to the protein and its binding pocket. It was observed that there was a slight increase in fluctuation of both the ligands till 10 ns and reached almost stable state without much fluctuation. After 20 ns simulation time, RMSD profile reaches its stability and no instability was observed for 3CLpro complex. However, PLpro complex followed the similar trend at the initial phase but they stabilized after 40 ns simulation time. The comparative analysis between these two complexes suggested that 3CLpro complex has less ligand flexibility than the PLpro ligand. The binding free energy of the two docked complexes was calculated and the average binding energy was found to be -69.77 ± 1.10 kJ/mol and -142.56 ± 0.90 kJ/mol for 3CLpro and PLpro complexes, respectively. The 3CLpro complex had more favourable binding with the ligand and higher binding energy represents more tight binding. From Figure 6 (d) it was observed that PLpro complex had similar MM-PBSA profile till 40 ns and some aberration were observed from 40 to 60 ns simulation time. However, during rest of the simulation period the complex, binding energy remained in similar pattern. On the other hand, 3CLpro complex had also some deviation in binding energy in 50 to 60 ns period and except that all simulation trajectory confirms its stability in MM-PBSA descriptor.

To understand the change and binding stability in the ligand-docked complexes, we also studied the protein-ligand interaction profile of the two complexes. Table 2 and Figures  7 and 8 showed the binding properties of protein-ligand complexes and the amino acid residues involved at different time intervals. Figures 7(a-e) and 8(a-e) showed the surface view of the topological changes of the protein-ligand complexes at 0, 25, 50, 75, and 100 ns. The study showed that binding of ligand amentoflavone to the 3CLpro protein is stabilized by 5, 2, 2, 5, and 5 hydrogen bonds at 0, 25, 50, 75, and 100 ns during simulation. The amino acid residues involved in the interactions, nature of bonding, and bond length were shown in Table 2 and Figures 7(f-j) and 8(f-j). Amino acid residue Glu166 was found to be a crucial residue forming H-bond with the ligand throughout the simulation period. Most of the residues were found to be repeated at different time intervals forming non-covalent interactions. After 100 ns of simulation, the amino acid residues involved in H-bond were found to be Glu166, Phe140, His172, Leu141 and His41 in 3CLpro protein. However, one amide-pi-stacked bond was observed at Asn142 and one pi-alkyl bond at Met49 position for 3CLpro complex. Interestingly, non-covalent interaction at the active groove of 3CLpro (His41, Leu141, Asn142, Met49, and Glu166) was observed after 100 ns simulation which correlates with possible inhibition of 3CLpro by the ligand molecules. However, at 100 ns simulation time, PLpro complex had two hydrogen bond at Glu167, Asp164, one electrostatic bond at Arg166, one Pi-pi-stacked bond at Tyr268, two pi-alkyl bond at Pro247 and Pro248 residue with gallocatechin gallate ligand. This PLpro complex also had non-bonded interaction at the active sites (Asp164, Pro248, and Tyr268) after dynamics simulation which is required for the inhibition of this protein. Moreover, complex from 25, 50, and 75 simulation trajectories also followed the similar trend binding pattern. Figure 9 showed the superimposition between pre-MD and post-MD structures of both the enzymes which showed structural changes with RMSD values 1.799 Å and 1.78 Å for PLpro and 3CLpro protein complexes, respectively. The RMSD profile from superimposition also confirmed the inflexibility of the docked complexes.

ADMET profile of a compound is an important property that predicts the druglikeness of a compound. Table 3 showed the various druglikeness properties of phytocompounds. According to Lipinski's rule, a compound may be regarded to have druglikeness property if the molecular weight of a compound is 500 Da, LogP value 5, number of H-bond donor 5, and H-bond acceptor is 10. Table 2 showed that most of the compounds possess druglikeness property except few compounds such as rhamnazin-3-O-rutinoside and theaflavin-3,3-digallate which violated the Lipinski's rule. The top binding ligands, amentoflavone and gallocatechin gallate showed violation in two parameters. Amentoflavone violated the rule in MW and HBD, while gallocatechin gallate violated in HBD and HBA. Topological polar surface areas of most of the compounds were below the 100Å 2 . Smaller the TPSA value higher the permeability of the molecule through the cell membrane, and therefore better distribution. The heat map of the ADMET properties of the top ten best binding phytocompounds with 3CLpro and PLpro are presented in the Figure 10 . It is observed from the study that almost all the phytocompounds have considerable ADMET properties. The absorption and permeability of all the top ten ligands showed moderate to high vales except few compounds such as rhamnazin-3-O-rutinoside and theaflavin-3-3-digallate which are predicted to have to permeability and absorption though human intestine. Amentoflavone showed low BBB property suggesting easy permeability and distribution to the brain cells while gallocatechin gallate showed high BBB value. Cytochrome-P450 is an important enzyme family necessary for detoxification and metabolism of drugs and toxic substances. Inhibition of this enzyme complex allows the drug molecules to reach the target site. 3CLpro top binding compound amentoflavone showed strong inhibitory property to Cyto-P450 enzyme family along with other phytocompounds. Gallocatechin gallate, on the other hand showed moderate cyto-P450 inhibitory property. The 3CLpro and PLpro best binding ligands were predicted to have low toxicity and carcinogenic property. Similarly, other phytocompounds are also predicted to have moderate to low toxicity activity.

The pandemic and severity of COVID-19 and its causative organism SARS-CoV2 have tremendously affected the health and wealth of the global community. Millions of lives have been lost so far and tens of millions are still undergoing treatment. With no clinically approved drugs, the current situation requires speedy investigation in the line of drug design to develop an effective drug to combat COVID-19. The molecular docking and molecular dynamics simulation are the most common and frequently used techniques to assess the interaction between ligand-protein complexes at the atomic level. By using these techniques novel inhibitors against disease-causing biological targets can be designed (Gogoi et al., 2019; Swargiary et al., 2020) . In the present study, we have investigated the 3CLpro and PLpro inhibitory property of 32 phytocompounds reported to have anti- protease activity by several researchers. The study found that amentoflavone and gallocatechin gallate isolated from Torreya nucifera and Pichia pastoris possess strong binding affinity to 3CLpro and PLpro proteins, respectively. Amentoflavone is a biflavonoid constituent of a number of plants including Torreya nucifera, Chamaecyparis obtusa, Ginkgo biloba, and many others (Pan et al., 2005; . The plant compound is known to possess several biological activities such as neuroprotector (Shin et al., 2006) , anti-inflammatory activity (Banerjee et al., 2002) , anticancer , and many others. Amentoflavone have also been reported by several researchers as inhibitor of viral protease enzyme (Sun et al., 1995; Jayadevappa et al., 2020) . Furthermore, rhamnazin-3-O-rutinoside, theaflavin-3,3-digallate, and kazinol-A also showed good binding affinity to the target proteins. Several recent studies have reported many phytochemicals that are predicted to have considerable binding affinity to proteases and other target protein of SARS-CoV2 (Qamar et al., 2020; Pan et al., 2005) . We have found that the highest binding energy of 3CLpro and PLpro was -9.4 kcal/mol and -8.8 kcal/mol with the top ligands. Similarly, most of the research findings showed that the binding energies of protein and ligand averages about -7.0 to -8.0 kcal/mol (Shin et al., 2006) . In a similar study, in silico analysis by Strodel et al. (2020) revealed amentoflavone as one of the top SARS-CoV2 3Clpro (Mpro) binding ligand (-9.28 kcal/mol). The strong binding property of amentoflavone has been attributed to the presence of extended ring structures combined with functional group. It was also reported that one of the phenolic ring make H-bond with Ser46 residue of protein and 11 van der Waals interactions. However, in the present study we have found that all the five H-bonds (Glu166, Arg188, Gln189, Met49, and His41) were formed by the phenolic rings of the 3CLpro, and Ser46 did not form any contact with the protein. Similarly, Mishra et al. (2020) also revealed that amentoflavone made strongest interaction with the SARS-CoV2 main protease (-9.13 kcal/mol). In a 10 ns MD simulation they have reported that the peak RMSF was observed in amino acid residue-106, 181, 46 and 151 which is different from our findings. We observed that the amino acid residues Ser1, Glu55, Ile59, Gln244, and Gln306 from SARS-CoV2 3CLpro exhibited more fluctuations. Gallocatechin gallate is a derivative of catechin, a natural phenol having several medicinal values. In the present study gallocatechin gallate have been found to bind strongly with PLpro. Nguyen et al. (2012) also reported the Figure 10 . ADMET properties of the best-binding phytocompounds with 3CL protease and PL protease of SARS-CoV2. HIA, human intestinal absorption; BBB, blood-brain barrier; H-HT, human hepatotoxicity; AMES, Ames mutagenicity, SkinSen, skin sensitization; DILI, drug induced liver injury; blue colored chemicals; PL, protease docked compounds.

inhibition property of epigallocatechin gallate and gallocatechin gallate against 3CLpro protein of SARS-CoV. Recent in silico study by Ghosh et al. (2020) revealed that gallocatechin gallate exhibited strong binding affinity to 3CLpro (-9.0 kcal/ mol). In another study, a similar compound, epigallocatechin gallate showed strong binding interaction with PLpro (-13.91 kcal/mol). Amino acid residues His89, Lys92, Trp93, Trp106, Asp108, Cys155, Asn156, and Lys157 were found to involve in H-bonding with the ligand (Laskar & Choudhury, 2020) . In the present study, gallocatechin gallate is found to make five H-bonds with PLpro pre-MD simulation.

It is also observed that the binding affinity of compounds depends upon the pharmacophore descriptors of the molecule. MD simulation studies have revealed that two best ligand molecules based on binding affinity showed more stability than the ligand-free protein structure. The higher RMSD values of apo structures revealed that, the 3CLpro and PLpro of SARS-CoV-2 have higher rigidity after the binding of the ligand suggesting the stability and effective binding of the ligands. Also, radius of gyration and solvent accessible surface area suggest the similar trend of the docked complex where two hit compounds exhibit less mobile nature and structural compactness. The tight packaging of the docked complex was illustrated and also no expansion of protein volume was maintained for both docked structure. Insilico druglikeness and ADMET properties are important parameters for the screening of the possible drug candidate. By using in silico tools it is now possible to predict the druglikeness of a compound based on its structure, physical, and chemical properties (Banerjee et al., 2002) . Lipinski's rule of five predicts whether a compound is orally absorbed or not and it depends on four main characteristics of a compound, molecular weight, lipophilicity, H-bond donor, and acceptor. According to the rule, a molecule will not be orally active if it violates two or more of the four rules (Lipinski, 2004) . The present study found that the best binding ligands possess the considerable property of druglikeness. A good drug candidate must also possess other properties such as high absorption by gastrointestinal tract, easy distribution, metabolism, excretion or elimination, and less toxicity called ADMET profile . Similarly, the present study showed that amentoflavone and gallocatechin gallate possesses good ADMET properties suggesting the possibility of potential drug candidates against SARS-CoV2. The high absorbability by GI, high permeability through cell membrane, resistant to on route metabolism, and less cytotoxicity properties suggest therapeutic potentiality of the phytocompounds. Other phytocompounds such as iguesterin, psoralidin, rhamnazin 3-o-rutinoside, theaflavin-3,3-digallate, and kazinol-A also showed promising ADMET properties with less side effects.

Plants have been the source of healthcare needs since ancient times. In this crucial juncture of COVID-19 pandemic, a speedy approval of effective drug and medicine is indeed the need of the hour. Plants, being a rich source of phytocompounds may be an effective strategy to combat SARS-CoV2. The present study revealed the binding affinities of phytocompounds to the active sites of two crucial proteins 3-chymotrypsin-and papain-like proteases of SARS-CoV2. Phytocompounds amentoflavone and gallocatechin gallate not only bound strongly with the target proteins but also stabilized the three-dimensional conformations of the protein structures after binding. Furthermore, insilico druglikeness and ADMET profiling of the compounds showed promising therapeutic potential. The present study therefore may be taken as the basis for further studies and to finally establish the efficacy of the phytocompounds.

Inhibition of TNFa-induced cyclooxygenase-2 expression by amentoflavone through suppression of NF-jB activation in A549 cells

The 2019-new Coronavirus epidemic: Evidence for virus evolution

The FDA-approved drug ivermectin inhibits the replication of SARS-CoV-2 in vitro

Inhibition of SARS-CoV 3C-like protease activity by theaflavin-3,3'-digallate (TF3)

Phytochemicals containing biologically active polyphenols as an effective agent against Covid-19-inducing coronavirus

Origin and evolution of pathogenic coronaviruses

SwissADME: A free web tool to evaluate pharmacokinetics, drug-likeness and medicinal chemistry friendliness of small molecules

Lipid14: The amber lipid force field

The price of innovation: New estimates of drug development costs

ADMETlab: A platform for systematic ADMET evaluation based on a comprehensively collected ADMET database

Anti-HCV, nucleotide inhibitors, repurposing against COVID-19

Molecular evolution of human coronavirus genomes

Mild or moderate covid-19

Evaluation of green tea polyphenols as novel corona virus (SARS CoV-2) main protease (Mpro) inhibitors -an in silico docking and molecular dynamics simulation study

Calculation of the total electrostatic energy of a macromolecular system: Solvation energies, binding energies, and conformational analysis

Energetically favorable anti-electrostatic hydrogen bonded cationic clusters in Ni (II) 3,5-dimethylpyrazole complexes: Anticancer evaluation and theoretical studies

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

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

In vitro inhibition of Hepatitis C virus protease and antioxidant by flavonoid glycosides from the Saudi costal plant Sarcocornia fruticosa

Assessment of structurally and functionally high-risk nsSNPs impacts on human bone morphogenetic protein receptor type IA (BMPR1A) by computational approach

Investigation of plant flavonoids as potential dengue protease inhibitors

Virtual screening and repurposing of FDA approved drugs against COVID-19 main protease

Comparative molecular investigation of the potential inhibitors against SARS-CoV-2 main protease: A molecular docking study

Phenolic phytochemical displaying SARS-CoV papain-like protease inhibition from the seeds of Psoralea corylifolia

YASARA -Yet another scientific artificial reality application

Fast empirical pKa prediction by Ewald summation

New ways to boost molecular dynamics simulations

Search for therapeutics against COVID 19 targeting SARS-CoV-2 papain-like protease: An in silico study. ResearchSquare (Preprint)

Prenylisoflavonoids from Erythrina senegalensis as novel HIV-1 protease inhibitors

Fatty acid synthase inhibition by amentoflavone suppresses HER2/neu (erbB2) oncogene in SKBR3 human breast cancer cells

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

Lead-and drug-like compounds: The rule-of-five revolution

Scutellaria baicalensis extract and baicalein inhibit replication of SARS-CoV-2 and its 3C-like protease in-vitro

Exploring the potent inhibitors and binding modes of phospholipase A2 through in silico investigation

Molecular docking and dynamics study of natural compound for potential inhibition of main protease of SARS-CoV-2

Natural compounds as potential inhibitors of novel coronavirus (COVID-19) main protease: An in silico study

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

Molecular Docking and 3d-pharmacophore modeling to study the interactions of chalcone derivatives with estrogen receptor alpha

Flavonoid-mediated inhibition of SARS coronavirus 3C-like protease expressed in Pichia pastoris

Ellagic acid & gallic acid from Lagerstroemia speciosa L. inhibit HIV-1 infection through inhibition of HIV-1 protease & reverse transcriptase activity

Open Babel: An open chemical toolbox

Use of artemisinin and its derivatives for the treatment of malaria in children

Update on therapeutic approaches and emerging therapies for SARS-CoV-2 virus

Amentoflavone and its derivatives as novel natural inhibitors of human Cathepsin B

Targeting SARS-CoV-2 spike protein of COVID-19 with naturally occurring phytochemicals: An in silico study for drug development

Diarylheptanoids from Alnus japonica inhibit papain-like protease of severe acute respiratory syndrome coronavirus

Dieckol, a SARS-CoV 3CL(pro) inhibitor

Chalcones isolated from Angelica keiskei inhibit cysteine proteases of SARS-CoV

Evaluation of polyphenols from Broussonetia papyrifera as coronavirus protease inhibitors

Structural basis of SARS-CoV-2 3CLpro and anti-COVID-19 drug discovery from medicinal plants

Identification of COX-2 inhibitors via structure-based virtual screening and molecular dynamics simulation

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

SARS-CoV 3CLpro inhibitory effects of quinone-methide triterpenes from Tripterygium regelii

Polyphenol amentoflavone affords neuroprotection against neonatal hypoxic-ischemic brain damage via multiple mechanisms

Papain-like protease (plpro) inhibitory effects of cinnamic amides from Tribulus terrestris fruits

High throughput virtual screening to discover inhibitors of the main protease of the coronavirus SARS-CoV-2

Test of the biflavones from Selaginella moellendorffii on the in vitro inhibition of HIV-1 protease

Antioxidant and antiproliferative activity of selected medicinal plants of lower Assam, India: An in vitro and in silico study. Anti-Cancer Agent in Medical Chemistry

On the origin and continuing evolution of SARS-CoV-2

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

In silico ADMET and molecular docking study on searching potential inhibitors from limonoids and triterpenoids for COVID-19

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

Specific plant terpenoids and lignoids possess potent antiviral activities against severe acute respiratory syndrome coronavirus

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

Potential therapeutic targets and promising drugs for combating SARS-CoV

Identification of novel proteolytically inactive mutations in coronavirus 3C-like protease using a combined approach

Authors thank Head, Department of Zoology, Bodoland University, India for providing necessary facility for carrying out this work. Authors also acknowledge the Head, Department of Genetic Engineering and Biotechnology, University of Rajshahi, Bangladesh, for providing necessary facilities.

Authors declare no conflict of interest.

This research did not receive any specific grant from funding agencies in the public, commercial, or not-for-profit sectors.

http://orcid.org/0000-0001-9594-3666 Shafi Mahmud http://orcid.org/0000-0003-1604-8626