key: cord-0953818-wg979q2d authors: Singh, Amit; Mishra, Abha title: Leucoefdin a potential inhibitor against SARS CoV-2 M(pro) date: 2020-06-17 journal: J Biomol Struct Dyn DOI: 10.1080/07391102.2020.1777903 sha: c9fb2e0bde084fc7f577a26b22a6f78cc47321b0 doc_id: 953818 cord_uid: wg979q2d Leucoefdin an important constituent of various fruits such as banana, raspberry, etc. was explored to target M(Pro) protease of SARS Co-V 2. Ligand was found to bind at active site of M(Pro) with large negative binding energies in molecular docking and simulation study. The docking results showed that Leucoefdin interacted with the M(Pro) by forming hydrogen bonds, at Leu 141, His163, His 164, and Glu 166. Other non-bonded interactions were seen at Met49, Pro52, Tyr54, Phe140, Leu141, Cys145 and Met165. Results of Leucoefdin was in coherence with the recently reported M(Pro) protease-inhibitor complex. It even displayed better binding energies (kcal/mol) in HTVS (-6.28), SP (-7.28), XP (-9.29) and MMGBSA (-44.71) as compared to the reference ligand [HTVS (-4.87), SP (-6.79), XP (-5.75) and MMGBSA (-47.76)]. Leucoefdin-M(Pro) complex on molecular dynamic simulation showed initial fluctuations in RMSD plot for a certain period and attained equilibrium which remained stable during entire simulation for 150 ns. RMSF of protein showed less secondary structure fluctuations and a greater number of H-bond formation with Leucoefdin during 150 ns simulation. Post simulation MMGBSA analysis showed binding energy of -45.98 Kcal/mol. These findings indicated the potential of Leucoefdin as lead compound in R&D for drug discovery and development against SARS CoV-2. Communicated by Ramaswamy H. Sarma Antiviral drugs target various stages of viral replication, like cell entry (fusion protein inhibitors), uncoating (ion channel blocker, capsid stabilizers), transcription of viral genome, translation of viral proteins, post-translational modification and release inhibitors. Till date, Food and Drug Administration (FDA) has not approved any specific drug against newly emerged SARS-CoV-2 virus . There are no clinical trial data supporting any prophylactic therapy. Many antiviral drugs are under clinical trial for the efficacy and safety against COVID-19. Several research and active clinical trials are underway especially on Protease inhibitors like Lopinavir, Ritonavir etc., RNA polymerase inhibitors such as Remdesivir, Favipiravir and Ribavirin etc, and antihelminths etc. Some other drugs like Baricitinib, Imatinib, Dasatinib, Cyclosporine are also undertaken for in vitro activity against SARS-CoV-2. Also, nonspecific anti-inflammatory and immunomodulators are being evaluated for its efficacy against COVID-19 (Agostini et al., 2018; Elfiky, 2020; Sarma et al., 2020; Wang et al., 2020; Wu et al., 2020) . The replicase gene of COVID-19 virus encodes two overlapping polyproteins, ppa and ppab which are vital for viral replication and transcription. A main protease (M Pro ) leads proteolytic processing of polyproteins and generate functional polypeptides. Therefore, M Pro can be an attractive target for antiviral drug design (Fearon et al., 2020; Jin et al., 2020) . The present study intended to find inhibitor of M Pro through molecular docking and simulation. The phytochemicals are rich source of many bioactive constituents which are essential for normal physiological functions, (USDA, 2017) . Leucoefdin, a member of catechol, and a polyphenol is colourless and linked to leucoanthocyanidins found in the Arachis hypogaea (Earth nut in seeds), Musa sp. (Banana, in the fruit), Phyllanthus emblica (Indian gooseberry), Vicia faba (bell-bean, in the seed), Nelumbo nucifera (lotusleaf) etc. (Drewes & Roux,1966; Ganguly et al., 1958) have been selected for the docking and simulation analysis. M Pro (PDB ID:5R84) (Fearon et al., 2020) and ligand Leucoefdin (Pubchem ID: 3081374) were obtained from RCSB protein data bank and Pubchem respectively. Structure based virtual screening and molecular docking was conducted using the Schr€ odinger package to find the best docked pose (LigPrep, 2015) . The docking was done with default parameters for grid and pose generation using HTVS, SP and XP mode. After XP study, the best interacting ligand was selected based on GlideScore (Friesner et al., , 2006 Halgren et al., 2004; Jacobson et al., 2004; Sastry et al., 2013) . SiteMap module of Schr€ odinger suite was used to forecast active sites of target protein structure. Potential sites of the target protein and site-score values were generated. Protein preparation wizard of Schr€ odinger suite was used to fix the drawbacks of protein. Further, protonation done by Epik module of Schr€ odinger suite. Water molecules more than 3 Å away from the ligand were removed before subjecting to protein minimization and hydrogen bond optimization. OPLSe force field employed in minimization studies. The structures of ligand were prepared with the LigPrep module, then energy minimisation was done in gas phase using Macromodel with the OPLSe force field. Glide-receptor grid generation module was used for docking studies (Halgren, 2007) . Co-crystallized ligand was first selected and around binding site M Pro ligand complex grid box was generated for molecular docking (Fearon et al., 2020) . Leucoefdin and known inhibitors from published sources Z31792168 was then docked with protein. The docking was done for high throughput virtual screening (HTVS), standard precision (SP) and extra precision (XP) to evaluate the binding energies with default parameters (Glide Schrodinger suite, 2015). First the co-crystallized ligand was extracted from the M Pro protein and re docked. Docking results of Leucoefdin and reference ligand was compared. The ligand with best XP score was further used for MD simulation studies. Molecular dynamics (MD) of Leucoefdin-M Pro complex was performed with Desmond module of Schrodinger Suite. The prepared complexes were solvated using orthorhombic solvent box and a solvent buffer (extended 10 Å outside the protein in all directions) and OPLS3 force field was applied. System was further neutralized by adding the desired number of counter ions with 0.15 M salt concentration. The six-stage system relaxation protocol of Schrodinger suite, was used before starting production runs (Harder et al., 2016; Macromodel, Schrodinger suite, 2015) . 2000 steps of steepest descent minimization in initial two stages were conducted with and without a restraint of 50 kcal/mol/Å2 on the solute atoms, then other four short MD simulations done with 12ps in NVT ensemble at 10 K with solute heavy atoms restrained with a force constant of 50 kcal/mol/Å2. Further 12ps MD simulation in NPT ensemble at10 K with the same restraint followed by 12ps MD simulation in NPT ensemble at 300 K with the same restraint and finally 24ps MD simulation in NPT ensemble at 300 K without any restraints. Unrestrained production simulations in the NPT ensemble was conducted for 150 ns at 300 K temperature and 1.01325 bar pressure. Relaxation time of 1ps was used with Nos e-Hoover chain thermostat and the isotropic Martyna-Tobias-Klein barostat was used with a relaxation time of 2ps. To evaluate short range interactions 9 Å cut off was used, while long range coulombic interactions were estimated with smooth particle mesh Ewald method (PME). r-RESPA integrator was used to study non-bonded interaction where every step was updated for short range forces and after every three steps the long range forces were updated. MD simulations for 150 ns was performed for M Pro -Leucoefdin complex using OPLS2005 (Optimized Kanhesia for Liquid Simulations) force fields. Trajectories were documented after every 4.8 ps, where 1.2 ps interval used for energy recording for 150 ns simulation. The trajectories were saved at 50.0 ps intervals for analysis (Macromodel, Schrodinger suite, 2015) . Prime-MM/GBSA (Molecular Mechanics/Generalized Born Surface Area) was used for the binding free energy calculation (Mobley & Dill, 2009; Prime, version 4.5, 2016; Kevin et al., 2006; Shivakumar et al., 2010; Wang, et al., 2015) . Almost all continents are having cases of SARS CoV-2, therefore, fast discovery of compounds against this problem is needed. Viral replication and transcription require a main protease M Pro for proteolytic processing of polyproteins to generate functional polypeptides of SARS CoV-2 replicase gene. M Pro being an attractive target for antiviral drug design as, it is not closely related to homologues in humans (Jin et al., 2020; Wu et al., 2020) . We have used database of 250 compounds (compiled from published literature on phytochemicals with antimicrobial properties) against M Pro protease and Leucoefdin was chosen for the present study based on its high binding energy and interaction property. Site Map analysis was made based on different parameters such as size, volume, enclosure, contact, amino acid exposure, hydrophobicity, hydrophilicity and donor/acceptor ratio for identification of potential sites. Results showed residues 41, 49, 52, 53, 54, 140, 141, 142, 144, 145, 163, 164, 165, 166, 187, 188, 189, 190 with a site score of 0. 89. The 306 amino acids length with a resolution of 1.83 Å of M Pro was taken for the study to find inhibitor from some selected phytochemicals. The M Pro (PDB ID:5R84) structure of COVID-19 virus was used to produce receptor grid for docking (Fearon et al., 2020) . Docking results showed ligand's interaction to the target with large negative binding energies (Table 1) Phe140, Leu141, Cys145 and Met165 (Liu et al., 2020) . M Pro Protein 6W63-complex with inhibitor X77 showed hydrogen bond interaction with residues Gly143, and Glu166 and similar non-bonded interactions at Tyr54, Phe140, Met165, Met49 and Leu141. Cys145 was making pi interactions with X77 (Mesecar, 2020). Leucoefdin-M Pro complex results was found to be comparable with other reported inhibitor. Molecular dynamics (MD) simulations represents dynamic fluctuations within the complex and interactions of ligand, lipid and water molecules. To study the dynamic interaction of viral protein with the docked Leucoefdin, molecular dynamics (MD) simulations of the M Pro complexed with Leucoefdin was performed using Desmond. It was used to evaluate the structural constancy and binding site adaptations to the docked ligand. MD simulations was run for the protein and ligand for 150 ns. The stability and fluctuations of the protein and ligand alone and in complex during the simulation were investigated and the resulting trajectory for complex was made with the backbone root mean square deviations (RMSDs) which is the average displacement of atoms from particular frame to a reference frame. Protein ligand binding overall was stable during entire simulation. Figure 2 displayed some fluctuation in Protein RMSD during initial 10 ns and remained stable for next 140 ns i.e. fluctuation remains within 1 A 0 . In case ligand initial fluctuation was observed till 70 ns thereafter it remained stable for next 80 ns that might be due to H-bond formation at Glu 166 (84% of the time) and other Hydrogen bonds formation as described in Figure 6a . It was evident from the Lig-fit-Prot of Figure 2 that the value was lower than the RMSD of the protein during entire simulation of 150 ns, thus it can be said that the ligand has not diffused away from the initial binding site and the interaction was stable (Figure 2) . Root mean square fluctuations (RMSF) was studied to recognize the critical residues involved in the major interactions with a ligand. The essential residues forming hydrogen bond with ligands being displayed by green bar during 150 ns simulation (Figure 3 ). The two major fluctuations in local domain of protein was seen in loop regions, first fluctuation was in between 1 and 10 residues and second after 300 residues. This might be due to the reason that N and C terminal fluctuate more than any other part of the protein (Figure 3 ). Figure 5 showed different types of contacts such as H bond, hydrophobic, ionic and water bridges between protein and ligand during simulation. Top panel indicated the sum of specific interaction protein makes with the ligand over the course of the trajectory. Specific residues interaction with the ligand in each trajectory frame was shown in bottom panel. Sometimes more than one specific contacts with ligands were formed with some residues which was represented by dark orange shade ( Figure 5 ). Furthermore, after 40 ns of simulation dark orange horizontal bar has been seen at Glu 166 which remained in hydrogen bond interaction with leucoefdin for entire 150 ns simulation. Figure 6a clearly indicated that Leucoefdin was interacting with protein throughout 150 ns simulation. This was further substantiated by the Figure 6b where the actual percentage of interaction was established between protein and ligand during simulation. Glu 166 was seen to interact 84% times with three hydrogen bonds formation. Ala 191 form 2 H bonds, other H bonds were also seen with Thr 190, Gln 192. Pro 168, Val 186 and Arg 188 were coordinated with water bridges. These interactions suggested that Leucoefdin-M Pro complex was stable throughout simulation. Complex was further subjected to Post simulation MM-GBSA analysis and binding energy found to be -45.98Kcal/mol. COVID-19 was first time reported in China, December 2019 subsequently more than 170 countries are affected with this disease. WHO declared COVID-19 as a pandemic. Lots of efforts being made by various researchers to develop drug against SARS CoV-2 virus. There are various targets like viral M Pro protease, RNA dependent RNA polymerase, host protease like TMPRSS2 and ACE2. We have performed computational study to find out inhibitor of M Pro protease. Study revealed that Leucoefdin has potential to inhibit M Pro protein, which is essential protease responsible for formation of functional viral polyprotein. Ligand protein complex showed good negative binding energies in HTVS, SP, XP and MMGBSA studies. MD simulations, study showed stable protein-ligand interaction. Thus, the present study gives an insight on a mechanism through which it interacts with viral protease. The finding also strengthens our understanding on these phytochemicals and their mechanism of action at molecular level. No potential conflict of interest was reported by the author(s). Coronavirus susceptibility to the antiviral remdesivir (GS-5734) is mediated by the viral polymerase and the proofreading Exoribonuclease A new flavan-3,4-diol from Acacia auriculiformis by paper ionophoresis SARS-CoV-2 RNA dependent RNA polymerase (RdRp) targeting: An in silico perspective PanDDA analysis group deposition-Crystal structure of COVID-19 main Glide: A new approach for rapid, accurate docking and scoring. 1. Method and assessment of docking Protein ligand interaction study during 150 ns simulation(green-Hbond; pink-ionic; blue-water bridges; violet-Hydrophobic); (b) Schematic ligand interaction with protein red charged, green hydrophobic, blue polar,ÁÀÁ salt bridge during simulation. accuracy Extra precision glide: docking and scoring incorporating a model of hydrophobic enclosure for protein-ligand complexes A study of leucoanthocyanidins of plants-I : Isomers of leucodelphinidin from karada bark and eucalyptus gum New method for fast and accurate binding-site identification and analysis Glide: A new approach for rapid, accurate docking and scoring. 2. Enrichment factors in database screening OPLS3: A force field providing broad coverage of drug-like small molecules and proteins A hierarchical approach to all-atom protein loop prediction The crystal structure of COVID-19 main protease in complex with an inhibitor N3. Protein DataBank Center for Structural Genomics of Infectious Diseases (CSGID) Structure of COVID-19 main protease bound to potent broad-spectrum non-covalent inhibitor X77 Binding of small-molecule ligands to proteins Schr€ odinger, LLC In-silico homology assisted identification of inhibitor of RNA binding against 2019-nCoV N-protein (N terminal domain Protein and ligand preparation: Parameters, protocols, and influence on virtual screening enrichments Prediction of absolute solvation free energies using molecular dynamics free energy perturbation and the OPLS force field Accurate and reliable prediction of relative ligand binding potency in prospective drug discovery by way of a modern free-energy calculation protocol and force field Remdesivir and chloroquine effectively inhibit the recently emerged novel coronavirus (2019-nCoV) in vitro A new coronavirus associated with human respiratory disease in China