key: cord-0910458-euuii4r0
authors: Ghosh, Rajesh; Chakraborty, Ayon; Biswas, Ashis; Chowdhuri, Snehasis
title: Identification of polyphenols from Broussonetia papyrifera as SARS CoV-2 main protease inhibitors using in silico docking and molecular dynamics simulation approaches
date: 2020-08-07
journal: Journal of biomolecular structure & dynamics
DOI: 10.1080/07391102.2020.1802347
sha: 6f6ed901624e3af675024b88de1c0af3f7bd9218
doc_id: 910458
cord_uid: euuii4r0

The current COVID-19 pandemic is caused by SARS CoV-2. To date, ∼463,000 people died worldwide due to this disease. Several attempts have been taken in search of effective drugs to control the spread of SARS CoV-2 infection. The main protease (Mpro) from SARS CoV-2 plays a vital role in viral replication and thus serves as an important drug target. This Mpro shares a high degree of sequence similarity (>96%) with the same protease from SARS CoV-1 and MERS. It was already reported that Broussonetia papyrifera polyphenols efficiently inhibit the catalytic activity of SARS CoV-1 and MERS Mpro. But whether these polyphenols exhibit any inhibitory effect on SARS CoV-2 Mpro is far from clear. To understand this fact, here we have adopted computational approaches. Polyphenols having proper drug-likeness properties and two repurposed drugs (lopinavir and darunavir; having binding affinity −7.3 to −7.4 kcal/mol) were docked against SARS CoV-2 Mpro to study their binding properties. Only six polyphenols (broussochalcone A, papyriflavonol A, 3'-(3-methylbut-2-enyl)-3',4',7-trihydroxyflavane, broussoflavan A, kazinol F and kazinol J) had interaction with both the catalytic residues (His41 and Cys145) of Mpro and exhibited good binding affinity (−7.6 to −8.2 kcal/mol). Molecular dynamic simulations (100 ns) revealed that all Mpro-polyphenol complexes are more stable, conformationally less fluctuated; slightly less compact and marginally expanded than Mpro-darunavir/lopinavir complex. Even the number of intermolecular H-bond and MM-GBSA analysis suggested that these six polyphenols are more potent Mpro inhibitors than the two repurposed drugs (lopinavir and darunavir) and may serve as promising anti-COVID-19 drugs. Communicated by Ramaswamy H. Sarma

COVID-19 accounted for $8,760,000 infected cases worldwide while $463,000 people died between January to mid-June 2020. This highly contagious febrile respiratory illness was declared as a pandemic disease on January 30 2020, by the World Health Organization (WHO) (Cucinotta & Vanelli, 2020) . China was the epicenter of this disease, but it rapidly spread throughout the globe . The United States remains the most affected country with $2,300,000 infected cases and out of which $122,000 people died due to COVID-19. Fever, cough, sore throat, runny nose and difficulty in breathing remain the main symptoms but it has been reported to be asymptotic for some individuals which in turn, accelerates the spread of this disease (N. Chen et al., 2020; Ren et al., 2020; Zhu et al., 2020) . The unavailability of suitable drugs or therapies for effective treatment until now has transformed this disease into a dangerous and life-threatening.

A novel coronavirus, severe acute respiratory syndrome corona virus-2 (SARS CoV-2) has been identified as the etiological agent of the disease which belongs to the genus b-coronavirus (Zheng, 2020) . The whole-genome sequence of this RNA virus revealed that it is highly similar to that of SARS CoV-1 with a 79.6% sequence identity . However, the sequence similarities vary significantly for different viral proteins . For example, the sequence of spike proteins (S-protein) is quite divergent throughout different coronavirus species (Li, 2016) . This may be a consequence of rapid mutations and recombination across the species. Besides this, the binding propensities of these spike proteins towards the host receptors vary across the species (Lan et al., 2020) . For instance, both SARS CoV-1 and SARS CoV-2 use the same host receptor (ACE2) and show affinity to the same binding site but their binding affinities to ACE2 vary due to slight interface sequence variations (Lan et al., 2020) . On the other hand, the sequence of some proteins such as the main protease (Mpro) is highly conserved throughout coronavirus species (Mirza & Froeyen, 2020) . The Mpro from SARS CoV-2 shares more than 96% sequence similarity with the same protease from SARS CoV-1 and MERS (Supplementary Figure 1) . This makes Mpro an ideal target for broad-spectrum anti-CoV therapy. Mpro [also known as 3CLpro (chymotrypsin-like protease)] is a cysteine protease, which is an analog to the main picornavirus 3C protease (Rota et al., 2003) . Mpro plays an important role in the replication process of single-stranded RNA from SARS CoV-2. It helps in the proteolytic cleavage at 11 sites involving the Leu-Gln#(Ser, Ala, Gly) sequence of the viral polyprotein and resulting in the release of a total number of 16 nonstructural proteins (nsps) (Fan et al., 2004; Rota et al., 2003) . Each of the protomers of the homodimeric SARS CoV-2 Mpro protein consists of three domains (Supplementary Figure 1 ). Domain I (amino acid residues 8-101) and domain II (amino acid residues 102-184) form a chymotrypsin-like architecture and these two domains are connected to the domain III (amino acid residues 201-303) via a long loop (Jin et al., 2020) . Among them, domain I and II are essentially b-barrels while, domain III mainly consists of a-helices (Jin et al., 2020) . The catalytic site/active site/substrate binding site comprising of cysteine (Cys145) and histidine (His41) amino acid moiety is located at the cleft of domain I and domain II (Jin et al., 2020) . Cysteine145 serves as a common nucleophile and plays a vital role in the proteolytic functioning of Mpro (Anand et al., 2003; Chou et al., 2003; Hsu et al., 2005) . Deprotonation of Cys-thiol followed by nucleophilic attack of resulting anionic sulfur on the substrate carbonyl carbon is the first step in the proteolytic process of Mpro (Hsu et al., 2005) . As a result, a peptide product having an amine terminus is released whereas the deprotonated form of histidine is restored. Subsequently, the thioester product is hydrolyzed to produce a carboxylic acid and the free enzyme is generated in the final step (Hsu et al., 2005) . Therefore, proteolytic processing as the functional significance of Mpro in viral propagation makes this protease an important drug target. Besides this, lack of any human homolog of Mpro makes it an ideal target for the development of drugs against COVID-19 infection (Kim et al., 2016) . In recent times many computational studies have been carried out to repurpose various FDA approved antiviral drugs for COVID-19 treatment by targeting SARS CoV-2 Mpro (Arun et al., 2020; Elmezayen et al., 2020; Kandeel & Al-Nazawi, 2020; Mahanta et al., 2020; Muralidharan et al., 2020) . Apart from these drugs, various medicinal phytochemicals have been proposed as potent SARS CoV-2 Mpro inhibitors (Aanouz et al., 2020; Bhardwaj, Singh, Sharma, Rajendran, et al., 2020; Das et al., 2020; Enmozhi et al., 2020; Gyebi et al., 2020; Islam et al., 2020; Joshi et al., 2020; Umesh et al., 2020) . Even many plant-derived natural polyphenols (Oolonghomobisflavan-A from the tea plant Camellia sinensis L., epigallocatechin gallate, epicatechingallate and gallocatechin-3-gallate from green tea Camellia sinensis) having antiviral properties are effective SARS CoV-2 Mpro inhibitors (Bhardwaj, Singh, Sharma, Rajendran, et al., 2020; Ghosh et al., 2020) . Polyphenols having glucosidase inhibitory properties are reported to show anti-viral activity against various coronaviruses such as SARS CoV-1 and MERS (Dan et al., 2019; Williams & Goddard-Borger, 2020; Zhao et al., 2015) . A group of polyphenols from Broussonetia papyrifera (B. papyrifera) are known to possess a-glucosidase inhibition activity (H. W. Ryu et al., 2010; Hyung Won Ryu et al., 2012) . In 2017, Lee and coworkers have found that ten of them [broussochalcone B or bavachalcone (C1), broussochalcone A (C2), 4-hydroxyisolonchocarpin (C3), papyriflavonol A (C4), 3 0 -(3-methylbut-2-enyl)-3 0 ,4 0 ,7-trihydroxyflavane (C5), kazinol A (C6), kazinol B (C7), broussoflavan A (C8), kazinol F (C9) and kazinol J (C10)] (structures mentioned in Figure 1 ) show inhibitory effect against Mpro from SARS CoV-1 and MERS (Park et al., 2017) . In the same study, investigators have also experimentally revealed the mode of enzymatic inhibitionand calculated the IC50 values of these polyphenols against Mpro (Park et al., 2017) . But whether these polyphenols exhibit any antiviral activity against SARS CoV-2 by inhibiting the enzymatic/ proteolytic activity of Mpro is far from clear. Therefore, in this study, we have examined the inhibitory potency of these ten polyphenols from B. papyrifera against SARS CoV-2 Mpro with the aid of in-silico docking studies, molecular dynamics simulations and MM-GBSA analysis. This study has revealed that six polyphenols [broussochalcone A, papyriflavonol A, 3 0 -(3-methylbut-2-enyl)-3 0 ,4 0 ,7-trihydroxyflavane, broussoflavan A, kazinol F and kazinol J] exhibit stronger binding affinity with Mpro and are possibly effective Mpro inhibitors.

The structures of ten B. papyrifera polyphenols [broussochalcone B or bavachalcone (C1), broussochalcone A (C2), 4-hydroxyisolonchocarpin (C3), papyriflavonol A (C4), 3 0 -(3-methylbut-2enyl)-3 0 ,4 0 ,7-trihydroxyflavane (C5), kazinol A (C6), kazinol B (C7), broussoflavan A (C8), kazinol F (C9) and kazinol J (C10)] were downloaded from PubChem database server in SDF format (https://pubchem.ncbi.nlm.nih.gov). Using the PyMol (DeLano, 2002) , all the SDF files were converted to PDB format and each of the polyphenols structures was optimized with B3LYP/6-31G Ã basis set by using Gaussian09 software (Frisch & Clemente, 2009 ). Each of the optimized polyphenols structure was then inserted to AutoDock Tools and standard processes were used to obtain the pdbqt files.

The crystal structure of the SARS CoV-2 Mpro was taken from the RCSB Protein Data Bank (http://www.rcsb.org) (PDB ID: 6LU7) (Jin et al., 2020) . After correcting the improper bonds, missing hydrogens, side-chain anomalies etc (if any), pdbqt file for Mpro was made through AutoDock Tools (Morris et al., 2009; Morris et al., 2008) .

The docking of Mpro with two anti-HIV drugs and ten polyphenols from B. papyrifera was performed with the aid of AutoDock Vina (Morris et al., 2009; Morris et al., 2008) . The binding affinities of polyphenols-Mpro were determined and analyzed using it. As per the position of active site region, the grid box was assigned with a 10.0 Å radius throughout the initial inhibitor, where the ligands can easily be fitted and which covers the entire active site pocket. The same grid box size and other parameters were used for docking studies of two anti-HIV drugs along with all the ten polyphenols to obtain different docked conformations. The bestsuited conformations with the lowest root mean square deviation (RMSD) values along with the highest Vina score were selected for Mpro and B. papyrifera polyphenols complexes. The output from AutoDock Vina was rendered with PyMOL and DS visualizer softwares (Biovia, 2017; DeLano, 2002) .

The molecular dynamics (MD) simulations were performed using the GROningen MAchine for Chemical Simulations GROMACS 2019 (Abraham et al., 2015) . The GROMOS9653a6 force field and SPC water model were used for all the MD simulations (Oostenbrink et al., 2004) . The ligand topologies were obtained from the PRODRG server (Schuttelkopf & van Aalten, 2004) . All bond lengths of protein and anti-HIV drugs/polyphenols were constrained using the LINCS algorithm, while water molecules were restrained by SETTLE algorithm (Hess et al., 1997; Miyamoto & Kollman, 1992) . The system was accommodated in a cubic box with a total number of 30226, 30199, 30198, 30204, 30203, 30202, 30200, 30198, 30196 water molecules containing the unligated Mpro, Mpro-C2, Mpro-C4, Mpro-C5, Mpro-C8, Mpro-C9, Mpro-C10, Mpro-darunavir and Mpro-lopinavir complexes, respectively. Each system was energy-minimized using the steepest descent algorithm and equilibrated to achieve the appropriate volume. The leapfrog algorithm with time step 2 fs was used and at every 5 steps, the neighbor list was updated. The Particle Mesh Ewald method is used to treat the Longrange electrostatics with cut off 1.2 nm and with a Fourier grid spacing of 1.2 nm (Essmann et al., 1995) . Periodic boundary conditions were applied in all three directions. Equilibration of the systems was carried out in two main stages. First, the system was allowed to heat gradually to 300 K in NVT ensemble using the v-rescale algorithm for 10 ns. Then NPT ensemble was employed for 10 ns by positional restraining of the complexes (unligated Mpro, Mprodarunavir, Mpro-lopinavir and Mpro-polyphenol complexes) slowly allowing the solvent molecules to relax around it. We have used the Parrinello-Rahman method and Berendsen barostat to maintain the pressure and temperature, respectively (Berendsen et al., 1984; Parrinello & Rahman, 1981 ). For each system, the average temperature and pressure values remained close to the desired values. The equilibrated systems were then subjected to unrestrained production MD simulations of 100 ns each, maintaining target pressure (1 bar) and temperature (300 K). The root mean square deviation (RMSD), the total number of hydrogen bonds, root mean square fluctuation (RMSF), the radius of gyration (Rg), solvent accessible surface area (SASA) for each system were calculated from the MD trajectories Kamaraj & Purohit, 2016; Rajendran, 2016; Rajendran et al., 2018) .

To evaluate the theoretical free energies of binding of ligands to the receptor, generally, two methods are commonly used (a) the molecular mechanics generalized Born surface area (MM-GBSA) and (b) molecular mechanics Poisson-Boltzmann surface area (MM-PBSA). These two methods are equally efficient in predicting the correct binding affinities (Jianzhong Chen, 2016; J. Chen et al., 2015; Hou et al., 2011; Venugopal et al., 2020) . Here we have used the MM-GBSA method to calculate the relative binding free energies of anti-HIV drugs and B. papyrifera polyphenols to Mpro. The free energy of binding can be calculated as DG bind ¼DH À TDS.

DH ¼ DE elec þDE vdW þDG polar þDG non-polar , where E elec and E vdW are the electrostatic and van der Waal's contributions, and G ploar and G non-polar are the polar and non-polar solvation terms, respectively. The polar contribution of the free energy is estimated by a generalized Born model with an external dielectric constant of 80 and an internal dielectric constant of 1, while the non-polar energy contribution is calculated from the solvent accessible surface area (SASA). As similar types of ligands bind to the receptor, the entropic contribution is neglected here. Therefore, our calculated values referred to as relative binding free energies (DG bind ).

MM-GBSA is a popular method to calculate binding energy, which uses energy properties of free ligand, free receptor and receptor-ligand complex for binding affinity calculation. The prime module of the Schrodinger suite (Schr€ odinger Release 2020-1: Prime, Schr€ odinger, LLC, New York, NY, 2020) was used for all MM-GBSA calculations.

Swiss ADME and pkCSM-pharmacokinetics online softwares were used for the prediction of different pharmacokinetic properties of ten polyphenols from B. papyrifera (Daina et al., 2017; Pires et al., 2015) . Levels of toxicity along with the drug-likeness properties of these ten polyphenols such as absorption, distribution, metabolism and excretion parameters were mainly scrutinized.

The resolved crystal structure of Mpro in complex with i) a Michael inhibitor N3 and ii) inhibitor alpha-ketoamide provide useful information about the structural integrity of Mpro which we have elaborately discussed in the introduction section (Jin et al., 2020; Zhang et al., 2020) . These two crystal structures (PDB ID 6LU7 and 6Y2E) also laid foundation towards the structure-based drug design against Mpro. Islam et al., 2020; Joshi et al., 2020; Umesh et al., 2020) . Many anti-HIV drugs (darunavir, lopinavir, atazanavir, etc) also have a good binding affinity towards the active site of Mpro (Beck et al., 2020) . In recent past, many investigators have chosen darunavir and lopinavir (structures mentioned in Supplemental Figure 2 ) as standard substrates and have compared the binding affinity and/or binding modes between various small molecules with that of "Mpro-darunavir/lopinavir complex" (Bhardwaj, Singh, Sharma, Rajendran, et al., 2020; Gyebi et al., 2020; Mahanta et al., 2020) . Thus, we have also decided to take these two anti-HIV drugs as standard Mpro inhibitors for this study.

Prior to conduct molecular docking studies, the polyphenols of B. papyrifera (structures shown in Figure 1) were screened based on their drug-likeness characteristics. Pharmacokinetics analysis revealed that the molecular weight (MW) of these seven polyphenols were less than 500 (ranging from 322 to 438) which suggested that all the polyphenols may easily transported, diffused and absorbed inside the body. The number of hydrogen bond donors (H-Do) was less than 5 and the number of hydrogen bond acceptors (H-Ac) was in range from 4 to 7 for these polyphenols, which are in accordance with Lipinski's rules. Furthermore, the topological polar surface area (TPSA) of all the polyphenols was found in the range of 55.76 to 131.36 Å 2 indicating good bioavailability of these polyphenols. The calculated intestinal absorption (IA) varies between 84 to 94% which signified good cell membrane permeability and oral bioavailability. All the polyphenols were negative towards AMES test and none of these polyphenols had shown hepatoxicity. However, pharmacokinetics analysis indicated a negative tolerance dose of three polyphenols (C1, C3 and C7) for humans (Table 1) . Such negative dose tolerance makes these three polyphenols susceptible to human use. As the rest of the polyphenols (C2, C4, C5, C6, C8, C9 and C10) harbor good and favorable pharmacokinetic properties (Table 1) , we selected them for testing their inhibition potency against Mpro.

All these seven polyphenols and two anti-HIV drugs (darunavir and lopinavir) were subjected to molecular docking studies to assess the polyphenol(s) exhibiting the higher or comparable binding energy to that of "Mpro-darunavir/lopinavir interaction". Darunavir interacted with Mpro via two hydrogen bonds [Gly143 (2.3 Å) and Glu166 (2.4 Å)], one Pi-sulfur bond (Met165) and multiple alkyl/Pi-alkyl bonds (Leu27, His41, Met49, Cys145 and His163) (Figure 2A and Table 2 ). It also formed many van der Waals interactions with different amino acid residues of Mpro ( Figure 2A and Table 2 ). Lopinavir formed only one hydrogen bond with Cys145 and several other non-covalent bonds with various important amino acid residues (such as Thr26, His41, Met49, Phe140, Glu166, Leu167, etc) within the active site of Mpro ( Figure 2B and Table 2 ). The binding affinity of darunavir and lopinavir towards Mpro was À7.4 and À7.3 kcal/mol, respectively (Table 2) . Interestingly, all the polyphenols (C2, C4, C5, C6, C8, C9 and C10) exhibited higher binding affinity (À7.6 to À8.2 kcal/mol) towards Mpro than that of darunavir and lopinavir ( Table 2 ). The highest binding affinity was observed for C5 (À8.2 kcal/mol) and the lowest one is for C6 (À7.6 kcal/mol).

Molecular docking studies further depicted that all these polyphenols efficiently interacted with different amino acid residues of domain I and II of Mpro (Figure 3-4 and Table 2 ). C2 formed hydrogen bonds with five amino acid residues (Thr26, Gly143, Ser144, Cys145 and Glu166) of Mpro ( Figure  3A , Table 2 ). Many other amino acid residues including His41 were also involved in different non-covalent interactions (van der Waals, Pi-sulfur, Pi-sigma and alkyl/Pi-alkyl) with C2 ( Figure 3A ). When C4 was docked into the active site of Mpro, three hydrogen bondinteractions [Leu141 (2.3 Å), Cys145 (2.3 Å) and Arg188 (2.1 Å)] and thirteen other noncovalent interactions (van der Waals and alkyl/Pi-alkyl) were evidenced ( Figure 3B and Figure 3C and Table 2 ). In case of Mpro-C6 complex, we observed two hydrogen bond interactions [Ser144 (2.4 Å, 2.5 Å)], twelve van der Waals interactions (Thr25, Phe140, Leu141, Asn142, Gly143, His163, Met165, Glu166, Arg188, Gln189, Thr190 and Gln192) and two other types of non-covalent interactions (alkyl and Pi-alkyl) ( Figure 3D and Table 2 ). When the other three polyphenols (C8, C9 and C10) were docked individually to Mpro, these complexes were stabilized by three number of hydrogen bond interactions and many non-covalent (C-H bond, van der Waals, Pi-alkyl etc) interactions ( Figure 4A -C and Table 2 ). Gly143 of Mpro formed two and Glu166 formed one hydrogen bond with C8 ( Figure 4A , Table 2 ).

However, C9 formedhydrogen bonds with Leu141, Gly143 and Met165 of Mpro and the amino acid residues of Mpro engaged in the hydrogen bonding with C10 were Ser144, His163 and Thr190, respectively ( Figure 4B -C, Table 2 ). C10 interacted with His41 and Cys145 of Mpro via Pi-sigma and Pi-alkyl interactions, respectively ( Figure 4B ). These two amino acid residues of Mpro formed Pi-Pi (His41) and Pi-alkyl (Cys145) bonds with C9 ( Figure 4A ). It was also evidenced that the other four polyphenols (C2, C4, C5, and C8) were involved in interactions with His41 and Cys145 via hydrogen bonding or other non-covalent forces (Figure 3 and 4) . Only polyphenol C6 had no interaction with these two key catalytic residue(s) of Mpro protease ( Figure 4C ). Thus, it can be concluded that all these polyphenols (except C6) may possibly inhibit the proteolytic activity of Mpro and may be effective candidates for the treatment of COVID-19 disease. We selected six Mpro-polyphenol complexes (Mpro-C2, Mpro-C4, Mpro-C5, Mpro-C8, Mpro-C9 and Mpro-C10) for performing the subsequent studies. 

In order to get an idea about the structural stability, conformational fluctuations, compactness and folding behavior of Mpro alone/Mpro (unligated) and Mpro complexed with two anti-HIV drugs as well as six polyphenols, we performed MD simulations for 100 ns using GROMOS9653a6 force field. The analysis of RMSD usually provides important information about the stability of the protein-ligand complex. Thus, we first estimated the RMSD of backbone alpha carbon atoms of all these systems ( Figure 5 ). The RMSD of Mpro (unligated) maintained a constant value ($0.21-0.22 nm) from 2 ns to 17 ns. Thereafter the RMSD value gradually increased till 25 ns and reached $0.35 nm. Then, the value .was slightly decreased and persisted at $0.31 nm from 65 ns till the end of the MD run. The RMSD values for both Mpro-darunavir and Mpro-lopinavir complexes were found to remain almost constant ($0.36-0.37 nm) from 10 ns to 100 ns with some marginal fluctuations ( Figure 5 ). The RMSD profiles of Mpropolyphenol complexes revealed that most of the complexes were stabilized quite quickly ( Figure 5) The magnitude of RMSD corresponding to three Mpro-polyphenol complexes (Mpro-C2, Mpro-C9 and Mpro-C10) attained an equilibrium value after 40 ns ($0.24 nm for Mpro-C2 and $0.26-0.27 nm for other two systems) and remained almost the same throughout the 100 ns simulation time span. On the contrary, the saturation of the RMSD curve for two other Mpro-polyphenol complexes (Mpro-C4 and Mpro-C8) was observed after 30 and 20 ns, respectively ( Figure 5 ). For Mpro-C5 complex, the RMSD value from 2 ns to 34 ns oscillated between $0.19 and 0.24 nm. Then, within next 8 ns, the value progressively increased to $0.35 nm and remained almost the same till 50 ns, Again the value started increasing and reached to $0.46 nm after 58 ns and remained almost the same till the end of MD run with some fluctuations ( Figure 5 ). The average RMSD values for Mpro (unligated), Mpro-darunavir complex and Mpro-lopinavir complex were found to be $0.31 nm, $0.36 nm and $0.37 nm, respectively, which are in agreement with previously reported values (Table 3 ) (Bhardwaj, Singh, Sharma, Rajendran, et al., 2020) . While the same for all the Mpro-polyphenol complexes were ranging from $0.23 nm to $0.34 nm with the lowest RMSD for the Mpro-C8 system and the highest one for Mpro-C5 system (Table 3) . Thus, it can be concluded that all Mpro-polyphenol complexes are stable. However, the stability of Mpro-C5 is least among all of them. These findings additionally indicated that the stability of these six Mpro-polyphenol complexes is comparable or relatively more than that of the Mpro-darunavir/Mpro-lopinavir complex. The conformational stability of these six Mpro-polyphenol complexes was further analyzed by estimating the total number of intermolecular hydrogen bonds formed during the entire 100 ns simulation time span (Table 3 ). The average number of intermolecular hydrogen bonds in the Mpro (unligated) system was 547. In Mpro-darunavir and Mpro-lopinavir complex, the existence of more number of intermolecular hydrogen bonds was evidenced (550 and 551, respectively). We also noticed a higher number of intermolecular hydrogen bonds in all selected Mpro-polyphenol complexes (552-564) ( Table 3) . Among these six complexes, the highest number of intermolecular hydrogen bonds (564) was observed when C8 was complexed with Mpro. These results support the RMSD data obtained from MD simulations. Based on these findings we can convincingly say that the complex originated from the binding of each selected polyphenol to Mpro is quite stable.

Next, we looked into the flexibility of different regions of Mpro by calculating the RMSF of alpha carbon atoms for all systems (Figure 6 ). It was quite evident from the RMSF profiles that all systems experience higher conformational fluctuations in domain III. Moreover, the fluctuations for the amino acid residues pertaining to domain III of the Mpro-C5 complex were highest among all the studied systems. In the case of Mpro (unligated) system, we additionally observed higher fluctuations (up to $0.6 nm) in a certain portion of domain I (residues 45-60). In fact, most of the amino acid residues within the domain I and II of this system had RMSF fluctuation below 0.3 nm. The average RMSF value for Mpro (unligated) system was 0.194 nm ( Table 3 ). The Mpro-darunavir and Mpro-lopinavir system experienced more or less similar conformational fluctuations to that of Mpro (unligated) system ( Figure 6 ). In fact, the fluctuations for the residues 45-60, were reduced upon the binding of lopinavir to Mpro (up to 0.35 nm). For both Mpro-darunavir and Mpro-lopinavir complexes, the average RMSF value was 0.195 nm (Table 3) . Furthermore, it was evidenced after analyzing the RMSF profiles of four Mpro-polyphenol systems (Mpro-C2, Mpro-C4, Mpro-C8 and Mpro-C10) that these complexes exhibited lower fluctuations (especially in domain I and II) as compared to the Mpro (unligated) system. The average RMSF values of these Mpro-polyphenol complexes were ranging between $0.141 nm to $0.157 nm (Table 2) . Even in Mpro-C9 complexes, the fluctuations of most of the amino acid residues (residues 134-144 in Mpro-C9 complex) residing at the domain I and II were reduced ( Figure 6 ). The average RMSF value of this complex is $0.159 nm (Table 3 ). The RMSF plot of the Mpro-C5 complex reflected that very few amino acid residues within domain I and II have an RMSF value of more than 0.25 nm ( Figure 6) . Interestingly, the RMSF values of several stretches within these two domains of this Mpro-polyphenol system (residues 21-24, 85-109 and 130-136) were more compared to that of Mpro (unligated) system, The average RMSF value of Mpro-C5 complex was $0.207 nm (Table 3) . Most importantly, the fluctuations of many key residues of the binding region of Mpro were ceased down after binding to these six polyphenols. These findings suggested that the conformational fluctuations of these six Mpro-polyphenol complexes are comparable or relatively less than that of the Mpro-darunavir/Mpro-lopinavir complex.

Afterward, we estimated the Rg to assess the compactness of all the complexes (Figure 7 and Table 3 ). The average Rg value for Mpro (unligated) and the other two complexes (Mpro-darunavir and Mpro-lopinavir) was almost identical ($2.20 nm). We observed a slightly higher average Rg value for six Mpro-polyphenol systems ($2.21-2.23 nm) (Table 3) . Thus, it can be suggested that all these Mpro-polyphenol complexes were slightly less compact in comparison to free Mpro/ Mpro (unligated) and other two Mpro-HIV drug complexes. Thereafter, SASA was employed to assess the extent of expansion of protein volume in each system (Figure 8 and Table 3 ). The average SASA values of Mpro-darunavir complex (151.154 nm 2 ) and Mpro-lopinavir complex ($151.283 nm 2 ) were found to be in the similar range with Mpro (unligated) ($151.448 nm 2 ) ( Table 3) . Higher SASA values were observed in all the Mpro-polyphenol complexes ($152.775-$156.399 nm 2 ). This increased SASA values indicated marginal expansion of the Mpro upon interaction with these six polyphenols.

Finally, we estimated the binding energy of "Mpro-polyphenols interaction" as well as "Mpro-HIV drugs interaction" using the MM-GBSA method. We utilized the docking conformation having the highest AutoDock Vina energy values. Table 4 , illustrated very high MM-GBSA binding free energies of our Mpropolyphenol complexes, which are comparable with the docking results. These higher MM-GBSA free energy results signify greater stability of these Mpro-polyphenol complexes. The MM-GBSA free energy values of Mpro-darunavir and Mpro-lopinavir complexes were found to be À35.65 kcal/mol and À40.39 kcal/ mol, respectively (Table 4) . On the contrary, all the six Mpropolyphenol complexes exhibited higher MM-GBSA binding energy than that of Mpro-darunavir complex and Mpro-lopinavir complex. The MM-GBSA free energies of six Mpro-polyphenol complexes were ranging from À41.32 kcal/mol to À56.23 kcal/mol (Table 4) .

Among all of them, the Mpro-C9 complex exhibited the highest binding free energies, while the Mpro-C8 system showed the lowest binding free energies. It is quite evident from these MM-GBSA values that all six polyphenols efficiently interacted with Mpro with higher binding free energy than that of "Mpro-darunavir/lopinavir interaction." The higher MM-GBSA values (DG bind ) in the case of these six Mpro-polyphenol complexes were mostly contributed by the coulombic interactions, SASA and hydrophobic interactions. Among these six Mpro-polyphenol complexes, Mpro-C2 showed the highest, whereas Mpro-C10 showed the lowest coulombic interactions. The maximum free energy contribution from SASA was experienced by Mpro-C9 and Mpro-C10 complexes, whereas the rest of the Mpro-polyphenol complexes showed similar SASA values. The Mpro-C9 showed the highest hydrophobic interactions, whereas the lowest hydrophobic interactions was observed in the case of Mpro-C2 complex. Other Mpro-polyphenol complexes showed similar hydrophobic interactions.

This study is aimed to test the inhibition potency of B. papyrifera polyphenols against SARS CoV-2 Mpro using a computational approach. The polyphenols which possess favorable drug-likeness characteristics [broussochalcone A (C2), papyriflavonol A (C4), 3 0 -(3-methylbut-2-enyl)-3 0 ,4 0 ,7-trihydroxyflavane (C5), kazinol A (C6), broussoflavan A (C8), kazinol F (C9) and kazinol J (C10)] including two anti-HIV drugs (darunavir and lopinavir) were subjected to molecular docking studies. All these polyphenols had higher AutoDock Vina energy values than the darunavir and lopinavir. Among them, six polyphenols (C2, C4, C5, C8, C9 and C10) had interaction with both the key catalytic residues (His41 and Cys145) of Mpro. The RMSD and RMSF profiles corresponding to these six Mpro-polyphenol complexes clearly suggested that they (complexes) are highly stable and experience less conformational fluctuations. The Rg and SASA analysis revealed that all Mpro-polyphenol complexes are slightly less compact and expand marginally. The existence of a higher number of intermolecular hydrogen bonds in the complexes with B. papyrifera polyphenols (C2, C4, C5, C8, C9 and C10) than in Mpro-darunavir/lopinavir complex suggesting greater stability of these polyphenols in the binding pockets of Mpro. These findings were further corroborated by MM-GBSA analysis. This analysis revealed that all Mpro-polyphenol complexes were more stable than Mpro-darunavir and Mprolopinavir complex. Therefore, it can be concluded that broussochalcone A, papyriflavonol A, 3 0 -(3-methylbut-2-enyl)-3 0 ,4 0 ,7-trihydroxyflavane, broussoflavan A, kazinol F and kazinol J were more effective Mpro inhibitors than earlier recommended anti-HIV drugs (darunavir and lopinavir). However, their inhibitory effectiveness and usage as anti-COVID-19 drugs should be thoroughly examined using various experimental studies.

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

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

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

Drug repurposing against SARS-CoV-2 using E-pharmacophore based virtual screening, molecular docking and molecular dynamics with main protease as the target

Predicting commercially available antiviral drugs that may act on the novel coronavirus (SARS-CoV-2) through a drug-target interaction deep learning model

Molecular dynamics with coupling to an external bath

A new insight into protein-protein interactions and the effect of conformational alterations in PCNA

Structural based study to identify new potential inhibitors for dual specificity tyrosine-phosphorylation-regulated kinase

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

Discovery studio modeling environment

Probing the structure-function relationship of Mycobacterium leprae HSP18 under different UV radiations

Drug resistance mechanisms of three mutations V32I, I47V and V82I in HIV-1 protease toward inhibitors probed by molecular dynamics simulations and binding free energy predictions

A comparative insight into amprenavir resistance of mutations V32I, G48V, I50V, I54V, and I84V in HIV-1 protease based on thermodynamic integration and MM-PBSA methods

Epidemiological and clinical characteristics of 99 cases of 2019 novel coronavirus pneumonia in Wuhan, China: A descriptive study

Binding mechanism of coronavirus main proteinase with ligands and its implication to drug design against SARS

WHO declares COVID-19 a pandemic

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

Benzonate derivatives of acetophenone as potent a-glucosidase inhibitors: synthesis, structure-activity relationship and mechanism

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

The PyMOL Molecular Graphics System

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

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

A smooth particle mesh Ewald method

Biosynthesis, purification, and substrate specificity of severe acute respiratory syndrome coronavirus 3C-like proteinase

Gaussian 09, Revision A. 01

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

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

LINCS: A linear constraint solver for molecular simulations

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

Mechanism of the maturation process of SARS-CoV 3CL protease

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

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

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

Mutational analysis on membrane associated transporter protein (MATP) and their structural consequences in Oculocutaeous albinism Type 4 (OCA4)-A molecular dynamics approach

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

Reversal of the progression of fatal coronavirus infection in cats by a broad-spectrum coronavirus protease inhibitor

Structure of the SARS-CoV-2 spike receptor-binding domain bound to the ACE2 receptor

Structure, function, and evolution of coronavirus spike proteins

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

Potential anti-viral activity of approved repurposed drug against main protease of SARS-CoV-2: An in silico based approach

Structural elucidation of SARS-CoV-2 vital proteins: Computational methods reveal potential drug candidates against main protease, Nsp12 polymerase and Nsp13 helicase

Settle: An analytical version of the SHAKE and RATTLE algorithm for rigid water models

AutoDock4 and AutoDockTools4: Automated docking with selective receptor flexibility

Using AutoDock for ligandreceptor docking

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

Conformational perturbation, hydrophobic interactions and oligomeric association are responsible for the enhanced chaperone function of Mycobacterium leprae HSP18 under pre-thermal condition

A biomolecular force field based on the free enthalpy of hydration and solvation: The GROMOS force-field parameter sets 53A5 and 53A6

Evaluation of polyphenols from Broussonetia papyrifera as coronavirus protease inhibitors

Polymorphic transitions in single crystals: A new molecular dynamics method

pkCSM: predicting small-molecule pharmacokinetic and toxicity properties using graphbased signatures

Structural analysis of oncogenic mutation of isocitrate dehydrogenase 1

Pathological role of a point mutation (T315I) in BCR-ABL1 protein-A computational insight

Identification of a novel coronavirus causing severe pneumonia in human: A descriptive study

Characterization of a novel coronavirus associated with severe acute respiratory syndrome

Polyphenols from Broussonetia papyrifera displaying potent alpha-glucosidase inhibition

Inhibition of xanthine oxidase by phenolic phytochemicals from Broussonetia papyrifera

PRODRG: A tool for highthroughput crystallography of protein-ligand complexes

Identification of new anti-nCoV drug chemical compounds from Indian spices exploiting SARS-CoV-2 main protease as target

Effect of hydrophobic and hydrogen bonding interactions on the potency of ss-alanine analogs of G-protein coupled glucagon receptor inhibitors

alpha-glucosidase inhibitors as host-directed antiviral agents with potential for the treatment of COVID-19

COVID-19 transmission through asymptomatic carriers is a challenge to containment

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

Inhibition of endoplasmic reticulum-resident glucosidases impairs severe acute respiratory syndrome coronavirus and human coronavirus NL63 spike protein-mediated entry by altering the glycan processing of angiotensin I-converting enzyme 2

SARS-CoV-2: An emerging coronavirus that causes a global threat

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

A novel coronavirus from patients with pneumonia in China

RG acknowledges IIT Bhubaneswar for providing fellowship. The authors thank IIT Delhi HPC facility for computational resources.

The authors declare that they have no conflicts of interest with the contents of this article.

http://orcid.org/0000-0002-1155-7862