key: cord-0879342-e510h0ob
authors: Solo, Peter; Maria doss, Arokia doss
title: Potential inhibitors of SARS-CoV-2 (COVID 19) spike protein of the and delta plus variant: In silico studies of medici delta nal plants of North-East India
date: 2021-10-22
journal: Curr Res Pharmacol Drug Discov
DOI: 10.1016/j.crphar.2021.100065
sha: 016e12a60e007e6e7de1c009a127f5307d176600
doc_id: 879342
cord_uid: e510h0ob

Phytochemicals of 38 Medicinal plants of North-East India, with anti-viral, anti-oxidant or anti-bacterial properties were screened for properties of drug likeness. 231 phytochemicals were screened with LIPINSKI rule of five to obtain 131 candidates, which were further screened with SWISS-ADME, to obtain 50 phytochemicals. These phytochemicals were docked with the spike protein of the Delta variant (B.1.617.2) and Delta-Plus (AY.1) variant of SARS-CoV-2 using Autodock Vina and MOE 09. The target proteins were constructed by homology modeling using Swiss-Model. Hydroxychloroquine, taken as a standard in docking analysis, exhibited a binding energy of −6.5 kcal/mol and −6.1 kcal/mol with respect to the Delta variant and Delta-Plus variant respectively. Among the 50 docked results most flavones showed very good docking scores. 3,5,8-Trimethoxy-6,7,4,5-bis(methylenedioxy)flavone, a Poly-Methoxyflavone, produced a highest docking score of −8.7 kcal/mol with respect to both the spike protein targets. Poly-Methoxyflavones and Poly-Ethoxyflavones exhibited good binding affinity for the target spike protein of SARS-CoV-2, and can be potential anti-viral drug candidates against the existing Delta variant of the SARS-CoV-2.

The current outbreak of severe acute respiratory syndrome coronavirus 2 (SARS CoV-2) has now undergone multiple significant mutations since its detection in 2019 in Wuhan, China. The spread of the Delta variant, which originated in India (Callaway, 2021) , has caused concern all over the world, with cases detected in over 96 countries (WHO, 2021) . The Delta variant has been denoted by the World Health Organisation as a "variant of concern". At present, the variant possess great threats to many countries like, the United States, Africa, Brazil Australia and Europe. India is still fighting a resurgence of the delta variant which appeared in the early part of 2021.

The Delta variant (B.1.617.2) has reported to be 60% more transmissible than the already highly infectious Alpha variant (B.1.1.7) (Callaway, 2021) , and is believed to spread faster than any other variants (Planas et al., 2021) . The recent studies on the variant has ignited fresh attention into how SARS-CoV-2 is able to adapt and mutate with the existing environment (Salvatore et al., 2021) . Another variant which is very similar to the Delta variant is the Delta plus variant (AY.1) which was first detected in Europe and was declared as a "variant of concern" by the U.K. governmental agency Health England. The delta plus variant is a sub lineage of the delta variant, with a notable difference of possessing K417N mutation in the spike protein. Most significant mutations in these variants have been occurring in the RBD region of the spike protein (Shu and McCauley, 2017: Khateeb et al., 2021) and these mutations corresponds to the increased transmissibility (Zhang et al., 2020; Erik et al., 2021) , increased immune evasiveness of the virus (Weisblum et al., 2020; Verma et al., 2021) and more flexibility to the spike protein to interact with the host receptors (Teruel N et al., 2021) .

With decrease in vaccine efficacy due to mutations (Noh et al., 2021) and the absence of strong anti-viral drug candidate against SARS-CoV-2, the world is still battling to overcome the Pandemic. In this situation, one can look into nature for a cure and a solution. Most traditionally used medicinal plants have phyto-constituents that are anti-viral, antiinflammatory, anti-oxidant and anti-microbial. These traditional plants can be investigated for potential anti-viral drug against SARS-CoV-2. A number of Insilco studies have been done with phytochemicals in pursuit of developing anti-viral drugs for SARS-CoV-2 (Pandey et al., 2020; Chojnacka et al., 2020; Basu et al., 2020; Kothandan et al., 2021; Prakrity et al., 2021) .

The North-east India is part of the 'biodiversity hotspot' in the world; it shares both the Himalaya and the Indo-Burma region. The region constitutes a rich reservoir of medicinal plants (Roy et al., 2015) and has a rich variety of flora and fauna. The present research tries to explore the Phytochemicals of some of these medicinal plants of this regions in search of a potential anti-viral drug against SARS-CoV-2 using in silico studies.

Thirty-Eight indigenous medicinal plants, listed in supplementary table 1, of North-East India were selected for the study as they possessed anti-oxidant, anti-viral and anti-microbial properties based on researches carried out in literature (Chase & Singh, 2013; Bora et al., 2007; Panda et al., 2018) . 231 phytochemicals, listed in supplementary table 1, of theses medicinal plants were identified from literature and their sdf-structures were downloaded from Pubchem database, www.pubchem.ncbi.nlm.nih.gov (Kim et al., 2016) .

Two levels of screening for drug-likeness were performed on the phytochemicals that were selected. All the 231 compounds were screened with LIPINSKI rule of five (RO5) (Lipinski, 2004; Jayaram et al., 2012) and the results are listed in supplementary table 2. The selected compounds obtained after screening with Lipinski's rule were further screened with SWISS-ADME server, which analyses many important aspects of the drug, like; lipophilicity, water solubility, pharmacokinetics, druglikeness and the chemical properties (Diana et al., 2017) . The SDF structures of the compounds, downloaded from PUBCHEM, were fed into the SWISS-ADME server and were converted to SMILES format. 50 phytochemical compounds, listed in supplementary table 3, were finally selected as suitable candidates for docking with the spike protein of SARS-CoV-2 and their structures were converted into pdb format using PyMOL (Schrödinger and Delano, 2020).

The 3D structure of the spike glycoprotein of the Delta (B.1.617.2) and Delta plus (AY.1) variants were not available at RCSB protein data bank, but were constructed with Swiss Model server (Waterhouse et al., 2018) . The target sequence of the variants with GenBank IDs QWE52264.1 (https://www.ncbi.nlm.nih.gov/protein/QWE52264.1) and QVO78843.1 (https://www.ncbi.nlm.nih.gov/protein/QVO78843.1) were downloaded in FASTA format from NCBI (https://www.ncbi.nlm.nih.gov/) website for the Delta and Delta-plus variants respectively. The selected Gen sequence for the delta variant (QWE52264.1) contains all the significant mutations like; L452R, T478K, P681R and D614G, Fig. 1(a) . The Gene sequence for the Delta plus (QVO78843.1) also has the significant mutations such; D614G, K417N, P681N, T478K and L452R, Fig. 1(b) . The Amino acid sequences of the surface glycoproteins were pasted in the Swiss Model server for construction of 3D model. BLAST and HHblits methods were used to search for template structure against the SWISS-MODEL template library. The same template with PDB ID 7krs.1.A was selected for both the delta (B.1.617.2) and delta plus (AY.1) variants for building the 3D models. The description of the template used for building the 3D model is given in table 1. The Structure obtained for the spike glycoprotein of the two variants has a QMEAN score of -1.61 and -1.57 for the Delta and Delta-plus variant respectively. In both the spike protein models chain A is slightly up compared to the other two chains which could be due to mutations (Akbulut, 2021) , Fig. 1 (a) and Fig. 1(b) . 

The modelled protein structures were prepared for docking using Autodock Vina (Trott and Olson, 2010). Ligands were removed, polar hydrogen bonds and Gasteiger charges were added and the structures were saved in pdbqt format. The 50 phytochemical ligands were also prepared using Autodock Vina and their structures were saved in pdbqt formats. The binding site for the target spike protein was chosen by performing a blind docking with Hydroxychloroquine as the ligand. The grid box contained both the S1 and S2 subunits of the spike protein. The most probable active site of the spike protein was identified in the receptor binding domain (319-541aa). For docking with MOE 09, the sdf-structure of the 50 phytochemical ligands were converted to pdb format using Pymol and the structures were optimized for docking using Amber12:EHT Forcefield. 3D-protonation and energy minimization using Amber12:EHT was used to prepare the pdb-structure of the target proteins that were constructed using homology modelling.

Docking studies was carried out using Autodock Vina. A grid box of 62X68X40 Å was assigned in the active site of the spike glycoprotein for both delta and delta-plus variant. The set of 50 ligands were screened with the receptor using Autodock Vina. Hydroxychloroquine was used as a standard in comparing the docking scores. The docking results, listed in Table 2 , were viewed and analysed with Discovery studio visualizer (BIOVIA, Dassault Systèmes, 2021). MOE docking tool was also used to dock the 50 ligands into the receptor binding domain of the spike proteins (319-541aa). The 5 finest docked postures were created applying a scoring J o u r n a l P r e -p r o o f 

The ligand, 3,5,3'-Trimethoxy-6,7:4',5'-bis(methylenedioxy)flavone, was taken to perform MD simulation with both the Target variants, B.1.617.2 and AY.1. VMD, NAMD and CHARMM-GUI was used to perform molecular dynamics simulation. CHARMM-GUI was used to obtain CGenFF topology file for the ligand (Jo, S. et al., 2008) . VMD was used to generate the topology file for the both the target proteins, which was then merged with the ligand and the complex was solvated (Humphrey et al., 1996) . NAMD was then used to perform simulation for the Protein-Ligand complex with CHARMM36m force field. (James et al., 2020) . All the Simulation analysis plots, i.e., RMSD, RMSF and H-bonds were prepared using QtGrace.

In this study 231 phytochemicals were initially selected from 38 medicinal plants, which were reported to exhibit anti-oxidant, anti-viral or anti-microbial properties. Screening of these phytochemicals for drug likeness was performed using LIPINSKI rule of five and SWISS-ADME. A list of 50 phytochemicals were finally selected for docking with the target spike glycoprotein of SAR-CoV-2 (Delta and Delta-plus variant). The target proteins, i.e., spike protein of the Delta and Delta-plus SARS-CoV-2 variant were constructed using homology modelling with SWISS MODEL. Hydroxychloroquine, with a binding affinity of -6.5 kcal/mol with the delta variant and -6.1 kcal/mol with the delta-plus variant, was used as a standard in the docking analysis. Most phytochemicals exhibited higher binding score than the standard. The top 15 phytochemicals in the binding score, listed in Table 2 , are all flavones. The binding affinity of the phytochemicals towards the two targets are quite similar with only a few exceptions.

The docking score using MOE 09 were almost in correspondence with the scores using Autodock Vina, accept that the binding energy computed using MOE are generally lesser for all the 50 ligands (supplementary table 4).

The docking analysis using Autodock Vina identified 3,5,3'-Trimethoxy-6,7:4',5'bis(methylenedioxy)flavone (Fig. 2) as having the highest binding affinity of -8.7 kcal/mol for both the target proteins, in comparison to the other 49 phytochemicals. This particular phytochemical is present in Nicotiana plumbaginifolia (Shajib et al., 2018; Nishtha et al., 2017) . The same ligand performed very well when docked with MOE, producing the highest score with the delta variant and the second highest score with the delta plus variant. It is quite interesting to note that, all the hits produced using MOE docking tools are all present in the top 16 hits produced using Autodock Vina, accept for the ligand Reticuline, which did not perform well in Autodock Vina.

The interactions of the ligand (3,5,3'-Trimethoxy-6,7,4',5'-bis(methylenedioxy)flavone) with the amino acids ( Fig. 3 & Fig. 4 ) at the active site are very similar for both the target variants. GLU340, VAL341, PHE347, SER349, ALA352, LYS356, ALA344, TYR351, TYR451 and ALA348 are the common Amino residues in both the docking sites. The interactions of 3,5,3'-Trimethoxy-6,7,4',5'-bis(methylenedioxy)flavone with the delta variant spike protein has an additional amino acid residue (SER399) at the binding site (Fig. 3 ). This could indicate that there is very less structural difference between the two target proteins at the active site. Rhamnocitrin -8.0 -6.1 -7.9 -6.4 12 1,5-Dihydroxy-3-methoxy-7methylanthracene-9,10dione -8.0 -5.6 -7.4 -6.6

13 Oxyberberine -8.4 -5.7 -8.1 -7.1 14

Rugosinone -8.1 -6.3 -7.8 -6.5 15 Jaceosidin -7.8 -5.7 -7.7 -6.9 Hydroxychloroquine -6.5 -6.1 -6.1 -6.2

Among the top 15 hits, Afzelechin, 5,7-dihydroxy-2-(4-hydroxyphenyl) chroman-4-one, Dihydrokaempferol, Kaempferol, Apigenin, Isorhamnetin, Chrysosplenetin, Resutin, Rhamnocitrin, and 3,5,3'-Trimethoxy-6,7:4',5'-bis(methylenedioxy)flavone have more structural similarities and have better binding energy. In all these ligands, A344, A346, A347, A348 and B165 are the most common amino acid residues in the binding site. The Amino acid residues A344, A347 and A348 are involved in Pi-Pi T-shaped and Pi-alkyl interactions with the two aromatic rings of the ligands. The amino acid residue A346 is involved in conventional hydrogen bond with the carbonyl oxygen between the aromatic rings and B165 forms hydrogen bond with the aromatic hydroxyl group (Fig. 3 & Fig. 4 ). The receptor binding domain (RBD) situated in the S1 subunit initiates viral entry into the host by binding to the cell receptor ACE2 in the aminopeptidase N (Huang, Y. et al., 2020) . All the common amino acid residues actively involved in the interaction with the top hit ligands are in the RBD, which suggest that the ligands can disrupt the viral binding to ACE2 receptor.

The RMSD plots of the MD simulation in Fig. 5 (A, B) suggest that the protein-ligand complexes are stable within acceptable deviation. The systems after 70 to 80 ns, becomes more stable in the course of the simulation. 3,5,8-Trimethoxy-6,7,4,5-bis(methylenedioxy)flavone with the Delta variant spike protein exhibited an average RMSD of 1.348Å, and the same ligand exhibited an average RMSD of 1.577Å with the Delta Plus Spike protein. Thus, these results indicate that the protein-ligand complex system remains stable amidst internal motions and fluctuations during the course of the simulation.

Both protein-ligand complexes showed similar RMSF plot, Fig. 5 (C, D), which also indicates their structural similarity. The stability of the complexes around the active RBD (319-541a) is noticeable in the plot, which could be attributed to the ligand recognition in the active site. Both the complexes showed an average RMSF of 1Å with significant fluctuations in the region1000-1100. These results show that the target proteins are stabilized by the binding of the ligand, 3,5,8-Trimethoxy-6,7,4,5-bis(methylenedioxy)flavone, at the active site. Based on the structural similarities of the tops hits, a ligand based pharmacophore model is generated by PharmaGist (Dina et al., 2008) which is shown in figure 6 . The tops 15 hits were submitted to PharmaGist in mol2 format. The five acceptor groups and the three aromatic pockets are all involved in binding with the active sites of the spike proteins. 

In conclusion, the study reports that many medicinal plants rich in Poly-Methoxyflavones and Poly-Ethoxyflavones, especially Nicotiana Plumbaginifolia, Bergenia Ciliata, Cardiocrinum Gigantum Wall, Cauletya Spicata Sm, Centella Asiatica Linn, Houttuynia Cordata Thunb, Mentha spicata, Paris Polyphylla Sm, Thalicrum foliolosum and Wedelia chinensis Osbeck, could be used to treat the outbreak of COVID-19. Further in vivo and in vitro studies on Poly-Methoxyflavones and Poly-Ethoxyflavones are needed to confirm the theory proposed by this paper. We assume that more Insilco studies on phytochemicals will reveal the potentials that are inbuilt in many traditional plants, to inhibit SARS-CoV-2.

Peter Solo: Methodology, Software, Investigation, Resources, Data curation, Visualization, Writing-Original draft preparation. M. Arockia doss: Conceptualization, Validation, Supervision, Writing-Reviewing and Editing. 

Mutations in the SARS CoV-2 spike protein may cause functional changes in the protein quaternary structure

Molecular docking study of potential phytochemicals and their effects on the complex of SARS-CoV2 spike protein and human ACE2

Medicinal plants used by the people of Northeast India for curing malaria

Delta coronavirus variant: scientists brace for impact

Ethnomedical plants used by Angami tribe of Nagaland

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

SwissADME: a free web tool to evaluate pharmacokinetics, druglikeness and medicinal chemistry friendliness of small molecules

PharmaGist: a webserver for ligand-based pharmacophore detection

Evaluating the Effects of SARS-CoV-2 Spike Mutation D614G on Transmissibility and Pathogenicity

Structural and functional properties of SARS-CoV-2 spike protein: potential antivirus drug development for COVID-19

VMD: visual molecular dynamics

Hydrogen Bond: ARG346, ASN354; Pi-Pi T-shaped

Alkyl: ARG346; Pi-Alkyl: ALA348, ALA344, LYS356 6 5,7-dihydroxy-2-(4-hydroxyphenyl) chroman-4-one -7.7 -5.4434 Hydrogen bond: ASN165, ARG346; Pi-Pi T-shaped: PHE347; Pi-Alkyl: ALA344

Hydrogen bond: ASN165, ARG346, VAL341; Pi-Pi T-shaped: PHE347; Pi-Alkyl: ALA344

Dihydrokaempferol -7.9 -5.6789 Hydrogen bond: ASN165, ARG346; Pi-Pi T-shaped: PHE347; Pi-Alkyl: ALA344

GLU 340, SER399, ARG346; Pi-Pi T-shaped: PHE347; Pi-Alkyl: ALA344

Stacked: PHE347; Pi-Alkyl: ARG346

Hydrogen bond: ASN165, ARG346; Pi-Pi T-shaped: PHE347

Stacked: PHE347; Pi-Alkyl: ARG346, ALA344, ALA348, LYS356 9 Apigenin -7.8 -5.8967 Hydrogen bond: ASN165, ARG346, SER399; Pi-Pi Tshaped

Stacked: PHE347; Pi-alkyl: ALA344

Hydrogen bond: ARG346, ASN354; Pi-Pi T-shaped: PHE347

Stacked: PHE347

PHE347; Alkyl: ALA344, LYS356; Pi-Alkyl: ALA344 -6.7 -5.2339 Hydrogen bond: ASN354; Pi-Sigma: PHE347

Hydrogen bond: ASP405; Alkyl: ARG408, PRO412, ALA411; Pi-Alkyl: TYR380 8-Dimethyl-4a,5-dihydroxy-5-isopropyl-4a

TYR380 13 3,8-Dimethyl-5-isopropyl-4a

Carbon Hydrogen bond: PHE347; Alkyl: ALA344, ARG346

Hydrogen bond: ARG346; Carbon Hydrogen bond: PHE347, SER399

Hydrogen bond: ASP405, GLN409, LYS417

Hydrogen bond: GLU340, VAL341, ARG346

Hydrogen bond: ASN165, GLU340, ARG346

Stacked: PHE347; Alkyl: ALA344, LYS356

Hydrogen bond: LYS378

Hydrogen bond: ASN165, VAL341, ARG346; Carbon Hydrogen bond: ALA352

Hydrogen bond: ASP405, GLN409, ASN417; Carbon Hydrogen bond: GLU406; Pi-Cation: ARG403; Pi-Alkyl: ALA352, ALA348, ALA344, LYS356; Pi-Alkyl: ALA348

Anion: ASP405; Alkyl: ALA411

Hydrogen bond: LYS378, GLU406, GLN409, LYS417; Carbon Hydrogen bond: GLY416; Pi-Cation: LYS378

Pi-Donor Hydrogen bond: LYS417; Pi-Alkyl: LYS417

Hydrogen bond: GLN409, ASN417; Pi-Cation: LYS378; Pi-Anion: GLU406

Hydrogen bond: ASP405, ASN417

Hydrogen bond: ASN417

Pi-Cation: ARG403; Pi-Anion: ASP405; Pi-Sigma: ALA411

Carbon Hydrogen bond: GLU340, VAL341

ALA348; Alkyl: ALA352, ALA344; Pi-alkyl: ALA344, PHE347

Hydrogen bond: ASN165, ARG346, SER399

T-shaped: PHE347; Amide-Pi Stacked: PHE347; Alkyl: ALA348

ALA348; Alkyl: ALA352

Carbon Hydrogen bond: GLU340, VAL341, PHE347, SER349; Alkyl: ALA352, LYS356, ALA344; Pi-LYS356, ALA344; Pi-Alkyl: TYR351, TYR451, ALA348

Carbond Hydrogen bond: LYS378, TYR380; Alkyl: PRO384, LEU455; Pi-Alkyl: PHE456

Hydrogen bond: LYS378, GLN409, ASN417; Pi-Anion: GLU406; Pi-Pi T-shaped: TYR380; Alkyl: PRO412

Alkyl: LYS417; Pi-Alkyl: TYR380

Hydrogen bond: SER469; Pi-Donor Hydrogen bond: GLN474; Alkyl: LYS458, PRO491; Pi-Alkyl: LYS458 26 Harman -6.8 -3.8876 Hydrogen bond: SER383, LYS417; Pi-Sigma: LYS417, LEU455

ALA344, ALA348 28 Harmol -6.6 -3.7002 Hydrogen bond: SER399; Pi-Pi T-shaped: PHE347; Pi-Alkyl: ALA344

Hydrogen bond: ASN165, SER399

T-shaped: PHE347; Pi-Alkyl: ALA344

Hydrogen bond: SER371, ARG454

Alkyl: PRO384; Pi-Alkyl: PRO412 Pi-Alkyl: PRO491, ARG457, LYS458, ALA372 30 n-cis-Feruloyltyramine

Hydrogen bond: GLN409, GLN414, ASP427; Carbon Hydrogen bond: TYR380, ASP405; Pi-Alkyl: ARG408

Hydrogen bond: GLU340, LYS356

Alkyl: LYS356; Pi-Alkyl: ALA344

Hydrogen bond: ASN165, VAL341; Alkyl: ALA348; Pi-Alkyl: ALA344

Hydrogen bond: GLU406, GLN409; Alkyl: PRO384; Pi-Alkyl: TYR421

Hydrogen bond: TYR369, GLU406, ARG408

Alkyl: ALA372; Pi-Alkyl: ALA372

Hydrogen bond: LYS378, SER383, LYS417; Pi-Pi Tshaped: TYR380; Alkyl: PRO384; Pi-Alkyl: LYS417

Hydrogen bond: ARG403, GLN409, GLN414, PHE429; Pi-Cation: ARG403; Pi-Anion: ASP405; Pi-Sigma: ALA411

Alkyl: ALA411; Pi-Alkyl: TYR380

Carbon Hydrogen bond: CYS379; Pi-Alkyl: LYS417

PHE347; Pi-Alkyl: ALA344

Hydrogen bond: GLN474; Pi-Donor Hydrogen bond: GLN474

Carbon Hydrogen bond: GLY416; Pi-Donor Hydrogen bond: LYS417; Pi-Sigma: LYS417

Carbon Hydrogen bond: ASP427; Pi-Cation:ARG403

Pi-Anion: ASP405; Pi-Sigma: ALA411; Pi-Pi T-shaped: TYR380; Alkyl: LYS417

Hydrogen bond: ARG403; Carbon Hydrogen bond: LYS378, ASP427

Pi-Anion: ASP405; Pi-Pi T-shaped: TYR380; Alkyl: PRO412, ALA411; Pi-Alkyl: ALA411 38 Noroxyhydrastinine -6.9 -5.2435 Hydrogen bond

Carbon Hydrogen bond: VAL341; Pi-Pi T-shaped: PHE347; Alkyl: ALA344, ALA348; Pi-Alkyl: PHE347, ALA344 -6.9 -5.7754 Hydrogen bond: GLU340; Carbon Hydrogen bond: VAL341; Pi-Pi Tshaped: PHE347

Carbon Hydrogen bond: LYS378, PRO384

Hydrogen bond: LYS417; Alkyl: LYS378, LEU455, PRO384; Pi-Alkyl: LYS417

Hydrogen bond: ARG403; Carbon hydrogen bond: ARG403; Pi-Cation: LYS378; Pi-Anion: GLU406; Pi-Pi T-shaped: TYR380

Alkyl: PRO412; Pi-Alkyl: TYR453

Carbon Hydrogen bond: GLU340, VAL341; Pi-Cation: LYS356; Pi-Anion: GLU340

Alkyl: ALA344, ALA348; Pi-Alkyl: LYS356, ALA344

Hydrogen bond: ARG403; Carbon Hydrogen bond: TYR380; Pi-Pi Tshaped: TYR380; Alkyl: ARG408

Hydrogen bond: GLN409, LYS417; Carbon Hydrogen bond: GLY416; Pi-Cation: LYS417; Pi-Pi T-shaped: TYR380

Hydrogen bond: ASN354; Carbon Hydrogen bond: ARG346, ASN450; Pi-Donor Hydrogen bond: SER349; Pi-Sigma: ALA348

PRO412 LYS356; Pi-Alkyl: ALA348

Carbon Hydrogen bond: GLU340, VAL341; Alkyl: ALA344, ALA348

Hydrogen bond: ARG403; Pi-Pi Tshaped: TYR380; Alkyl: ALA411

Hydrogen bond: ARG346, SER399; Pi-Pi Stacked: PHE347; Alkyl: ALA344, ALA348; Pi-Alkyl: ARG346

SER399; Alkyl: ALA344, ALA348; Pi-Alkyl: ALA344

Hydrogen bond: ASN165, ARG346, ALA348; Alkyl: ALA344; Pi-Alkyl: ARG346

Hydrogen bond: ASN165, ARG346, ALA348

Alkyl: ALA344; Pi-Alkyl: LYS356

Carbon Hydrogen bond: ASN370, ARG457; Pi-Alkyl: ARG457

Hydrogen bond: ASP405; Carbon Hydrogen bond: ASP427; Pi-Cation: LYS378

Hydrogen bond: LYS378, ARG403; Pi-Pi T-shaped: TYR380; Alkyl: PRO412, ALA411; Pi-Alkyl: PRO412, TYR380 -7.2 -5.2180 Hydrogen bond: GLU406; Pi-Pi Tshaped: TYR380

TYR380, PRO412, ALA411 47 Cleroindicin-E -5.0 -3.4657 Carbon Hydrogen bond: ARG457

Hydrogen bond: LYS378, SER383, LYS417; Carbon Hydrogen bond: LYS378; Pi-Pi T-shaped: TYR380

Hydrogen bond: ARG403, GLN414, ASN417; Pi-Cation: ARG403; Pi-Anion: ASP405; Pi-Sigma: ALA411; Pi-Pi T-shaped: TYR380

TYR421, LYS417

Carbon Hydrogen bond: GLY416; Pi-Sigma: LYS417

The authors are grateful to the authorities of St. Joseph's College (autonomous) for providing the space and the equipment for performing computational analysis.

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