key: cord-335610-3v8140b6
authors: Prasanth, D. S. N. B. K.; Murahari, Manikanta; Chandramohan, Vivek; Panda, Siva Prasad; Atmakuri, Lakshmana Rao; Guntupalli, Chakravarthi
title: In silico identification of potential inhibitors from Cinnamon against main protease and spike glycoprotein of SARS CoV-2
date: 2020-06-22
journal: J Biomol Struct Dyn
DOI: 10.1080/07391102.2020.1779129
sha: 
doc_id: 335610
cord_uid: 3v8140b6

Cinnamon has been utilized to remedy a lot of afflictions of humans. Literary works illustrate that it possesses numerous biological activities. Our research study is intended to recognize the phyto-derived antiviral substances from Cinnamon against COVID-19 main protease enzyme and to understand the in silico molecular basis of its activity. In the present study, 48 isolates compounds from Cinnamon retrieved from the PubMed database, are subjected to docking analysis. Docking study was performed using Autodock vina and PyRx software. Afterwards, admetSAR, as well as DruLiTo servers, were used to investigate drug-likeness prophecy. Our study shows that the nine phytochemicals of Cinnamon are very likely against the main protease enzyme of COVID-19. Further MD simulations could identify Tenufolin (TEN) and Pavetannin C1 (PAV) as hit compounds. Utilizing contemporary strategies, these phyto-compounds from a natural origin might establish a reliable medication or support lead identification. Identified hit compounds can be further taken for in vitro and in vivo studies to examine their effectiveness versus COVID-19.

respectively (Fehr et al., 2017) . SARS originated in Guangdong province, China in 2002 and its eventual global expanded (de Wit et al., 2016; WHO, 2003) . Chinese horseshoe bats act as SARS-CoV reservoir hosts (Lau et al., 2005; Li et al., 2005) . Intermediate hosts such as civet cats and raccoon puppies, often marketed as food outlets in Chinese wet markets, enabled human transmission (Guan et al., 2003) . Currently, no new antivirals or licensed vaccinations are available to counter SARS, although traditional preventive steps, including travel bans although patient isolation, eventually stopped the SARS pandemic in 2002 and 2003. New respiratory infectious disease appeared in Wuhan, Hubei Province, China Wang et al., 2020; Zhu et al., 2020 ). An initial outbreak epidemic was connected to seafood center, possibly involved in animal exposure. Eventually, human-to-human infection emerged (Chan et al., 2020) and the now-called coronavirus disease 19 circulated exponentially in China. A unique coronavirus, SARS-coronavirus 2 (SARS-CoV-2), closely associated to SARS-CoV, was found in clinics and is suspected to be the emerging lung disease's etiologic agent (Chan et al., 2020) . On Feb 12, 2020, China recorded 44,730 lab-confirmed disorders, including 8,204 serious cases and 1,114 deaths (WHO, 2020) . Illnesses are also identified in 24 non-China nations, correlated with foreign travel. Whether sequencing correlations among SARS-CoV-2 and SARS-CoV translates into related biological properties, namely outbreak possibility, is currently unclear (Munster et al., 2020) .

Virus entry in cells is regulated by spike (S) glycoprotein; spike 1 (S1) surface unit helps the virus to be connected to cellular receptors. Cellular proteases cleave the S protein at the S1/S2 and S20 locations to allow viral particle entry. The viral capside is then fused with the cellular membrane, a subunit S2-guided process (Hoffmann et al., 2020) . It was established that angiotensin-converting enzyme 2 (ACE2) (15) mediates SARS-CoV entry and that serine protease TMPRSS2 is responsible for S protein cleavage (Glowacka et al., 2011; Li et al., 2003) (Hasan et al., 2020) . Analysis of receptor binding motif (RBM) sequences within the receptor binding domain reveals that it is accountable for binding to ACE2 and that residues have been preserved by SARS-CoV and SARS-CoV-2, indicating that binding to ACE2 could be identical while the same residues are absent in other coronaviruses (Ge et al., 2013; Menachery et al., 2019) . Some human ACE2 antibodies avoided SARS-CoV and SARS-CoV-2 infections (Boopathi et al., 2020) .

Proteases are one of the important and best studied proteins in corona viruses (Khan, Jha, et al., 2020) . It is very essential for the processing poly proteins that are translated from viral RNA (Hilgenfeld, 2014) . Inhibition of this viral protease enzyme wills results into the blocking of viral replication. In this scenario, the development of inhibitors to the SARS-CoV-2 main protease gains importance in the drug discovery process against COVID-19. Recently Zhang et al solved the crystal structure of inhibitor bound SARS-CoV-2 main protease (Khan, Zia, et al., 2020; Zhang et al., 2020) .

The development of new therapeutics is an expensive and time consuming process. Normally it will take years to get the newly developed drug for the treatment (Elfiky, 2020; Enayatkhani et al., 2020; Muralidharan et al., 2020; Pant et al., 2020; Wahedi et al., 2020) . Drug repurposing is an efficient strategy in medicinal chemistry to bring faster and effective solutions to the unmet medical needs (Gupta et al., 2020) (Elfiky & Azzam, 2020; Sinha et al., 2020) . Repurposing of drugs derived from natural origin is considered as crucial therapeutic approach for the treatment of COVID-19 (Enmozhi et al., 2020) , considering the fast pace of its spread around the world (Elmezayen et al., 2020) .

In the present research, we have chosen Cinnamon, which is a potent source of antiviral agents (Aanouz et al., 2020) . Cinnamon is a traditional Indian medicine which is being used from hundreds of year to relieve various lung-related disorders includes pneumonia, infectious disease, as well as malignant pleural effusion (Lai et al., 2018; Townsend et al., 2013) . Recently, several studies also provided scientific data to support and unveil its antiparasitic, antihypertensive, antidiabetic, anti-hyperlipedimic, anti-oxidant, anti-inflammatory, analgesic, antimicrobial, antiviral, antitumor, anti-hypertension, anti-hyperlipemic, gastro-protective and immunomodulation activities.

Based on the above properties of Cinnamon, this research aimed to show a variety of active compounds across all Cinnamon varities and decide whether and how they interact with proteins i.e. main protease (Joshi et al., 2020) and spike protein, that are essential in the management of SARS-CoV-2.

Within this research, a dataset of active phytochemicals were acquired from Indian Medicinal Plants, Phytochemistry, and also Therapeutics data source (Chaudhuri et al., 2018; Mohanraj et al., 2018) .

The X-ray crystal structures of Main protease and spike receptor domain complexed with ACE2 (PDB ID: 6LU7, 6LZG) were downloaded from the RCSB PDB (Protein Data Bank) database (Islam et al., 2020; Sarma et al., 2020) . The Graphical User Interface program ''Auto-Dock Tools'' was used to prepare, run, and analyze the docking simulations. Kollman united atom charges, solvation parameters and polar hydrogen's were added to the receptor for the preparation of protein in docking simulation. Since ligands are not peptides, Gasteiger charge was assigned and then non-polar hydrogens were merged. AutoDock requires pre-calculated grid maps, one for each atom type, present in the ligand being docked as it stores the potential energy arising. This grid must surround the region of interest (active site) in the macromolecule. (Morris et al., 2009) (Figure 1 ).

The crystal 3D structure of the following active compounds of Cinnamon were retrieved from PubChem database (O'Boyle et al., 2011) . Drug-likeliness properties of ligands were analyzed for the selected active compounds using DruLiTo software (Pangastuti et al., 2016) .

Autodock 4.0 methodology was validated with the respective co-crystallized ligands of target proteins to ensure the virtual screening process. Autodock 4.0 represents valid RMSD score and accurate binding with target receptor. In this context, Mpro (PDB ID: 6LU7) was tested with its co-crystalized inhibitor N3 (n-[(5-methylisoxazol-3-yl)carbonyl]alanyl-l-valyln$1$-((1r,2z)-4-(benzyloxy)-4-oxo-1-f[(3r)-2-oxopyrrolidin-3yl]$9#methylgbut-2-enyl)-l-leucinamide) (Agostino et al., 2009; Cosconati et al., 2010) .

Accurate prediction of active sites is an important tool in bioinformatics. In this study, for spike protein (PDB ID: 6LZG), the active site was predicted by using Biovia Drug discovery studio visualizer 2020 (Design, 2014) .

Molecular testing of all the compound libraries was performed using PyRx software by autodock wizard as the engine for docking (Dallakyan & Olson, 2015) . During the docking period, the ligands were considered to be flexible, and the protein was supposed to be rigid. The configuration file for the grid parameters was generated using Grid box for 6LU7 (x ¼ À17.59, y ¼ 15.81, z ¼ 63.53) and 6LZG (x ¼ À26.73, y ¼ 8.32, z ¼ À14.07) in PyRx respectively (Hall & Ji, 2020) . The application was also used to know/predict the amino acids in the active site of the protein that interact with the ligands. The results less than 1.0 Å in positional root-meansquare deviation (RMSD) were considered ideal and clustered together for finding the favourable binding. The highest binding energy (most negative) was recognised as the ligand with maximum binding affinity. Visual examination of the docking site was performed using Biovia Drug discovery studio 2019, and the results were validated using Autodock Vina (Seeliger & de Groot, 2010 ).

ADMET of the ligands is pharmacokinetic properties calculation that is required to be examined to establish their function inside the body. The ADMET inheritance of the ligands was studied, making use of admetSAR (Cheng et al., 2012; Yang et al., 2019) .

Prophecy of Cinnamon for antiviral activity was created with the assistance of software, PASS. PASS is a computer system based program utilized for the prognosis of various sorts of physiological actions for multiple compounds consisting of phytoconstituents. The estimated activity of a substance is predicted as probable activity (Pa) and probable inactivity (Pi). The substances revealing Pa higher than Pi are actually the only components thought about as feasible for a specific medical activity ( 

The crystal structure of Main protease (6LU7) and Spike receptor-binding domain complexed with its receptor ACE2 (6LZG) with selected top ligands identified from docking analysis such as Tenufolin (TEN) and Pavetannin C1 (PAV) were subjected to molecular dynamics using gromacs GPU enabled package. Ligand topology was selected from Prodrug server. The pdb2gmx, a module of GROMACS was used to add hydrogens to the heavy atoms. Prepared systems were first vacuum minimized for 1500 steps using the steepest descent algorithm. Then the structures were solvated in a cubic periodic box with a water simple point charge (SPCE) water model. The complex systems were subsequently maintained with an appropriate salt concentration of 0.15 M by adding suitable numbers of Na þ and Cl À counter ions. The system preparation was referred on the basis of previously published paper (Gangadharappa et al., 2019) . Each resultant structure from the NPT equilibration phase was subjected for final production run in NPT ensemble for 50 ns simulation time. Trajectory analysis was performed by using the gromacs simulation package of RMSD and RMSF (Mohankumar et al., 2020) . Molecular Mechanics Poisson-Boltzmann surface area (MM-PBSA) approach was employed to understand the binding free energy (DG binding) of an inhibitor with protein over simulation time. A GROMACS utility g_mmpbsa was employed to estimate the binding free energy (Kumari et al., 2014) . To obtain an accurate result, we computed DG for the last 20 ns in with dt 1000 frames.

For the main protease, target structure (PDB-6LU7) with reported inhibitor as co-crystallized ligand was retrieved from PDB for docking purpose. To further validate the docking methodology, RMSD value was calculated. RMSD score of the co-crystalized (internal) ligand and extracted internal ligand of the docked target protein-ligand complex structure served as control docking model as represented in Figure 2 . The docking outcome showed that Autodock 4.0 determined the optimal orientation of the co-crystallized ligand. RMSD value 0.271 indicates the methodology is accurate to predict the binding affinity for unknown ligands. For spike protein, target structure with reported inhibitor is not available in the PDB collection or literature. We have utilized the active site prediction option of Discovery Studio Visualizer (DSV) 2020. After loading the protein structure (PDB-6LZG) to DSV, it reads the protein and highlights the probable active site i.e. residue information in yellow color. For the following spike protein, three residues i.e. HIS374; HIS378; GLU402 were recognized as crucial amino acids at the active pocket. Considering all the three residues, grid file was generated with following dimensions (x ¼ À26.73, y ¼ 8.32, z ¼ À14.07) using PyRx program and progressed for molecular docking. Further Molecular Dynamics integrated with MM-PBSA calculations has validated the hit molecules identified from docking.

In order to identify a prospective candidate for managing COVID-19, molecular docking was executed over 48 phytoconstituents acquired from different species of Cinnamon on the binding pocket of enzyme COVID-19 (PDB ID: 6LU7; 6LZG). From the literature it was clearly indicated that virus enters into the host cell with the assistance of ACE2 receptor. S-protein binds directly to the Angiotensin Converting Enzyme 2 (ACE2) receptor of the human host cell surface -thus enabling virus entry and replication. We have taken the PDB ID: 6LZG where spike protein is complexed with ACE2 to investigate the binding affinity of cinnamon derivatives with the complex (Sinha et al., 2020) .

All these 48 compounds were docked against the target enzyme COVID-19 and ranked based on their dock score. Compounds possessing dock score of À7.0 or even less are thought about a better representative for restraint of the COVID-19. A comprehensive evaluation can be done by referring to table in supplementary file. This table exemplifies the list of active molecules acquired after docking studies. Those active molecules possess dock score value of À7.0 or lower were picked. Total of 7 compounds was selected based on the binding interactions with 6LU7 and 6LZG (Figures 3-10) . Out of the seven compounds, Tenufolin exhibited the bestdocked score (À8.8 Kcal/mol) with SARS-CoV2 Main Proteases (6LU7) and Pavetannin C1 exhibited the best-docked score (À11.1 Kcal/mol) with SARS-CoV2 spike protein (6LUZ).

The rigid docking results were envisioned utilizing Discovery studio for evaluation of communications. The best binding postures of protein-ligand communications were envisioned and charted in Tables 1 and 2. Main Protease, involved in severe acute respiratory syndrome (SARS), displayed best docking score of À8.8 kcal/mol with Tenuifolin, among the phytochemicals. Rest of the compounds 6-Glucopyranosylprocyanidin B1, Cinnamtannin-B1, Kaempferol 3-alpha-L-arabinofuranoside-7-rhamnoside, Pavetannin-C1, Proanthocyanidin-A2 and Procyanidin-B7 also interacted well with main Protease showing moderate binding energies of À7.6 kcal/mol, À8.4 kcal/mol, À8.1 kcal/mol, À7.3 kcal/mol, À8 kcal/mol and À8.2 kcal/mol respectively (Table 1 , Figures  1-4) . Meticulous in silico analysis revealed that all the phytoconstituents interacted decently at the active site of main protease, reflecting probable inhibitory tendencies against the COVID-19.

Tenufolin has the best correlation with main protease (6LU7) protein complexes. The main protease with tenufolin Spike protein (6LZG), associated with SARS was found to exhibit the best possible interaction with Pavetannin C1 (À11.1 kcal/mol) among the phytochemicals ( Table 2) . Rest of the phytoconstituents 6-Glucopyranosylprocyanidin B1, Cinnamtannin-B1, Kaempferol 3-alpha-L-arabinofuranoside-7rhamnoside, Proanthocyanidin-A2, Procyanidin-B7, Tenuifolin displayed moderate binding with spike protein at the active site with binding energies of À9.9 kcal/mol, À10.2 kcal/mol, À8.7 kcal/mol, À9.4 kcal/mol, À9.6 kcal/mol and À8.7 kcal/ mol respectively. Meticulous interface profiling showed promising phytochemicals inhibitory activity against COVID-19's Spike Protein (Figures 6-9) .

Pavetannin C1 has the best correlation with main protease (6LZG) protein complexes. 

The physicochemical properties of the chosen seven active compounds were studied on DruLiTo software. As all the compounds are of natural source, none of them obeyed Lipinski's rule (Table 1) . Kaempferol 3-alpha-L-arabinofuranoside-7-rhamnoside and Tenufolin shows higher TPSA (63.22, 52.6) and AMR (315.59, 172.9) (Table 1 and Figure 11 ). TPSA, as well as AMR, are fundamental physicochemical properties mostly entailed in drug absorption, transport and penetration mechanism (Ertl et al., 2000) .

The ADMET properties of the ligands were assessed, making use of admetSAR. ADMET properties for the substances in the research study were evaluated, making use of admetSAR. All the substances revealed excellent human intestinal absorption (HIA), blood-brain barrier (B.B.B.) infiltration. None of the compound was found carcinogenic. All the Table 4 .

The biological activity spectra of previously identified phytoconstituents were obtained by online PASS version. These predictions were interpreted and used in a flexible manner and given in Table 5 .

RMSD analysis of the protein gives insights and variations in its structural confirmation during the simulations, which - Figure 11 . The three dimensional structures of in silico active ligands. A: 6-Glucopyranosylprocyanidin B1; B: Cinnamtannin B1; C: Kaempferol 3-alpha-L-arabinofuranoside-7-rhamnoside-3; D: Pavetannin_C_1; E: Proanthocyanidin-A2; F: Procyanidin_B7; G: Tenufolin. includes the stability of protein and confirms whether simulation has equilibrated. The 6LU7 complex with TEN and 6LZG complex with PAV RMSD plot of backbone atoms was shown in Figures 12 and 13 . RMSD was plotted for 6LU7-TEN and 6LGZ-PAV complex structures that converged during the 50 ns MD simulation. Both complex structures were found stable after 10 ns. The average values of 6LU7-TEN was $0.45 nm and 6LGZ-PAV was $0.24 nm of 50 ns simulation. This indicates the protein complex structures were stable during the md simulation. RMSF analysis measures the fluctuations of each residue during simulation. Ligand binding poses energy and interaction is directly dependent on residual fluctuation (RMSF) values. The 6LU7 complex with TEN and 6LZG complex with PAV RMSF plot for each residue was shown in Figure 13 . It was observed that residues in loop region are more fluctuated during the simulation. This indicates that the protein was not fluctuated in the 50 ns simulation periods.

RBD domain superimpose analysis was carried out with open and close confirmation of the protein. Initial and last confirmation of the dynamics trajectory was analyzed using the Pymol software. RMSD ¼ 2.192 (572 to 572 atoms) indicates that there are some structural changes between the protein structure. Superimposed structure was depicted in Figure 14 .

We have reported the total binding energy (deltaG value) of complex contributed from various energy terms after MD simulations. The binding free energy value of 6LU7_TEN was found to be À123.949 þ/À 16.613 kJ/mol which indicates the inhibitor has higher affinity with the 6LU7 complex and binding free energy value of 6LZG_PAV was found to be À158.870 þ/À 30.378 kJ/mol which indicates that the inhibitor has higher affinity with 6LZG complex. The 6LU7_TEN and 6LZG_PAV systems were submitted for residue level MM-PBSA binding energy calculation. Hotspots in interaction with the inhibitor (TEN) complex with 6LU7 are as follows:

LYS137 À1.392 kcal/mol, ASP197 À1.088 kcal/mol, TYR239 À2.051 kcal/mol and GLU288 À1.103 kcal/mol. The total of hotspot interaction binding energy was À5.634 kcal/mol. Hotspots in interaction with the inhibitor (PAV) complexed with 6LZG are as follows:

SER44 À 1.061 kcal/mol, HIS345 À0.247 kcal/mol, ALA348 À3.580 kcal/mol, ASP350 À9.864 kcal/mol, ARG514 À0.212 kcal/mol and TYR515 À6.672 kcal/mol. The total of hotspot interaction binding energy was À21.636 kcal/mol. MD trajectories analysis and MM-PBSA results revealed that TEN complex with 6LU7 showed moderate binding with key residues and PAV complex with 6LZG formed strong binding with key residues towards COVID-19 targets (Table 6 ).

Coronaviruses have a long history of infecting humans and animals and causing respiratory, digestive, liver and central nervous system diseases . A novel emerged SARS-CoV-2 is presenting major threats to human health nowadays (Zhu et al., 2020) . The primary focus has been on clinical management which includes the prevention of infection, control measures and supportive care. Currently, no specific clinical therapeutics are available for the treatment of SARS-CoV-2-mediated infections (Zhou et al., 2020) . Thus, Table 5 . Results of PASS calculations for antiviral activity of isolated phyto-constituents from ginger.

Main predicted activity by PASS online the need of the hour is to identify and characterize novel drug candidates to overcome the health loses caused by SARS-CoV-2. With this new breakthrough of Mpro structure in COVID-19, it has offered an astounding possibility to recognize the prospective drug candidates for the effective therapy of coronavirus. In this context, natural products have gained importance as potent anti-viral agents during recent years (Lin et al., 2014; Martinez et al., 2015) . Considering the immediate need of therapeutics against COVID-19 and services of natural products in drug discovery, we have screened phytoconstituents from Cinnamon as novel hit molecules against Mpro, of SARS-CoV-2 for the identification of Mpro inhibitors to provide natural scaffolds for drug development.

Generally speaking, after virtual screening, recognizing lead compounds focuses on three main criteria: binding intensity, molecule linkage associations and therapeutic characterization. Lead compounds are more definitely molecules with very small binding energies, intense hydrogen and hydrophobic bond interactions, and relatively good ADMET properties. Therefore, top seven compounds were picked as suitable hit molecules for further study based on their low binding energies and relatively good poses inside the active site pocket after virtual screening. Out of 48 candidates, seven compounds displayed a higher binding affinity and least binding energy with the main protease enzyme and spike protein. Tenuifolin has exhibited highest dock score with least binding energy of À8.8 Kcal/mol and found to make six hydrogen bonds with six amino acids, i.e. ARG A:105 (6.34), GLN A:107 (4.08), GLN A:110 (3.43, 5.02), THR A:111 (3.39), ASP A:295 (3.97) and one hydrophobic interaction with PHE A:294 (4.35). Except with, ARG A:105, remaining, the bond length of hydrogen bonds is <5 A , which indicates the bonding is stronger and formed stable complexe. Whereas, Gingerenone A has a binding energy of À6.5 kcal/mol, associates with three hydrogen bonds with GLN A: 107 (4.04), GLN A:110 (3.79), THR A:111 (4.55) and hydrophobic interactions with VAL A:104 (4.98), ILE A:106 (5.48), PHE A:294 (4.73). These two compounds have the least binding affinity in comparison with other ligands due to the formation of more hydrogen bonds with the proteins. All the ligands are involved in hydrophobic interactions; mainly, two amino acids were involved, i.e. PHE A:294, VAL A:104 and ILE A:106.

Physicochemical characteristics were analyzed in light of the traditional drug-like laws that are implemented in the early stages of drug development, namely Lipinski's Five law (RO5). Lipinski's rule of Five (RO5) indicates that a medication molecule appears to show good oral bioavailability, smooth membrane permeability, and strong gastrointestinal absorption in human intestine while its log p 5; MW 500 dalton; HBAs 10 and HBDs 5. In our study, all the selected candidates with strong dock score have violated the Lipinski's rule particularly in two main properties i.e. molecular weight >500 and hydrogen bonding acceptor >10. A significant class of drug molecules omitted from the initial 'ruleof-five' study are natural products. The high efficiency of natural product-based drug development can be due to the fact that evolutionary selection, which determines structural prerequisites for protein binding, has biologically pre-validated their chemical structures. It also has the benefit over conventional semi-synthetic methods to refine a structurally diverse natural product lead without the concomitant rise in molecular size. Therefore, synthetic biology may be a great lead optimization method to preserving or growing ligand effectiveness.

In the In silico measurement of intestinal absorption, both phytochemicals were theoretically strongly soluble in the gastrointestinal tract (Table 5 ). Despite effective distribution of medications in the human body, significant toxic effects such as carcinogenicity (Benigni & Bossa, 2011; Liebler & Guengerich, 2005) became important concerns. Accordingly, in the early phases of drug growth, identification and evaluation of novel product candidates is proposed to be important. None of the compounds have demonstrated in silico toxicity in this analysis.

The biological activity spectra of all selected compounds were determined by using an online version of PASS software. Results obtained (Table 6 ) were interpreted and used in a flexible manner. Thus, we anticipate that the consumption of Cinnamon has the potential to boost immunity to fight against COVID-19 infections.

In this research study, we have actually made use of Bioinformatics resources, PyRx and also Autodock-Vina and GROMACS to identify the potent molecules from Cinnamon against COVID-19 Main Proteases and spike protein, which participate in a vital part in Corona virus propagation. Final in silico validation by molecular dynamics integrated with free energy calculations and per residue analysis have demonstrated stable complexes with main protease and spike protein. Overall results proposed that tenufolin and Pavetannin C1 are the potent compounds which has good binding efficacy with main protease and spike proteins. Comparatively, Pavetannin C1 has strong affinity towards spike protein and which need to verified experimentally. Natural substances are deemed key targets of drug discovery and lead greatly to innovative product design and production initiatives. As the compounds are bulky in nature with high molecular weight, so further optimization is required. Further, its efficiency can be proved by its in-vitro, in-vivo and clinical studies. 

Moroccan medicinal plants as inhibitors of COVID-19: Computational investigations

Molecular docking of carbohydrate ligands to antibodies: Structural validation against crystal structures

Mechanisms of chemical carcinogenicity and mutagenicity: A review with implications for predictive toxicology

Novel 2019 Coronavirus structure, mechanism of action, antiviral drug promises and rule out against its treatment

A familial cluster of pneumonia associated with the 2019 novel coronavirus indicating person-to-person transmission: A study of a family cluster

Innovation and trends in the development and approval of antiviral medicines: 1987-2017 and beyond

admetSAR: A comprehensive source and free tool for assessment of chemical ADMET properties

Coronaviruses as the cause of respiratory infections

Virtual screening with AutoDock: Theory and practice

Small-molecule library screening by docking with PyRx

SARS and MERS: Recent insights into emerging coronaviruses

Pharmacophore and ligand-based design with Biovia Discovery StudioV R

Natural products may interfere with SARS-CoV-2 attachment to the host cell

Novel guanosine derivatives against MERS CoV polymerase: An in silico perspective

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

Reverse vaccinology approach to design a novel multi-epitope vaccine candidate against COVID-19: An in silico study

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

Fast calculation of molecular polar surface area as a sum of fragment-based contributions and its application to the prediction of drug transport properties

Middle East respiratory syndrome: Emergence of a pathogenic human coronavirus

Structural insights of metallo-beta-lactamase revealed an effective way of inhibition of enzyme by natural inhibitors

Isolation and characterization of a bat SARS-like coronavirus that uses the ACE2 receptor

Evidence that TMPRSS2 activates the severe acute respiratory syndrome coronavirus spike protein for membrane fusion and reduces viral control by the humoral immune response

PASS-assisted exploration of new therapeutic potential of natural products

Isolation and characterization of viruses related to the SARS coronavirus from animals in southern China

In-silico approaches to detect inhibitors of the human severe acute respiratory syndrome coronavirus envelope protein ion channel

A search for medications to treat COVID-19 via in silico molecular docking models of the SARS-CoV-2 spike glycoprotein and 3CL protease

A review on the cleavage priming of the spike protein on coronavirus by angiotensin-converting enzyme-2 and furin

From SARS to MERS: Crystallographic studies on coronaviral proteases enable antiviral drug design

SARS-CoV-2 cell entry depends on ACE2 and TMPRSS2 and is blocked by a clinically proven protease inhibitor

Clinical features of patients infected with 2019 novel coronavirus in Wuhan, China. The Lancet

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

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

Targeting SARS-Cov-2: A systematic drug repurposing approach to identify promising inhibitors against 3C-like Proteinase and 2'-O-RiboseMethyltransferase

Identification of chymotrypsin-like protease inhibitors of SARS-CoV-2 via integrated computational approach

PASS assisted prediction and pharmacological evaluation of novel nicotinic analogs for nootropic activity in mice

g_mmpbsa-a GROMACS tool for high-throughput MM-PBSA calculations

Clinical practice guidelines for diagnosis and management of cough-Chinese Thoracic Society (CTS) Asthma Consortium

Severe acute respiratory syndrome coronavirus-like virus in Chinese horseshoe bats

Structure of SARS coronavirus spike receptor-binding domain complexed with receptor

Angiotensin-converting enzyme 2 is a functional receptor for the SARS coronavirus

Elucidating mechanisms of drug-induced toxicity

Antiviral natural products and herbal medicines

Antiviral drug discovery: Broad-spectrum drugs from nature

Trypsin treatment unlocks barrier for zoonotic bat coronavirus infection

PASSassisted exploration of antidepressant activity of 1, 3, 4-trisubstitutedb-lactam derivatives

Design and molecular dynamic Investigations of 7,8-dihydroxyflavone derivatives as potential neuroprotective agents against alpha-synuclein

IMPPAT: A curated database of Indian medicinal plants

AutoDock4 and AutoDockTools4: Automated docking with selective receptor flexibility

A novel coronavirus emerging in China-key questions for impact assessment

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

Open Babel: An open chemical toolbox

Natural bioactive compound from Moringa oleiferaagainst cancer based on in silico screening

Peptide-like and small-molecule inhibitors against Covid-19

In-silico homology assisted identification of inhibitor of RNA binding against 2019-nCoV N-protein (N terminal domain)

Ligand docking and binding site analysis with PyMOL and Autodock/Vina

An in-silico evaluation of different Saikosaponins for their potency against SARS-CoV-2 using NSP15 and fusion spike glycoprotein as targets

From SARS coronavirus to novel animal and human coronaviruses

Effects of ginger and its constituents on airway smooth muscle relaxation and calcium regulation

Stilbene-based natural compounds as promising drug candidates against COVID-19

A novel coronavirus outbreak of global health concern

Summary of probable SARS cases with onset of illness from 1

Coronavirus disease 2019 (COVID-19): Situation report

admetSAR 2.0: Web-service for prediction and optimization of chemical ADMET properties

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

Network-based drug repurposing for novel coronavirus 2019-nCoV/ SARS-CoV-2

A novel coronavirus from patients with pneumonia in China

Authors are thankful to K L College of Pharmacy, Koneru Lakshmaiah Education Foundation and M.S. Ramaiah University of Applied Sciences for providing necessary support to conduct the research and Siddaganga Institute of Technology, Tumakuru, India, for their encouragement and support in MD simulations.

All the authors declare no conflict of interest.

Vivek thank Karnataka Biotechnology and Information Technology Services (KBITS), Bangalore for providing the grant to establish the computational facility under Biotech Policy-II, Biotechnology Finishing School and Biotechnology Skill Enhancement Programme (BiSEP) at the Department of Biotechnology, Siddaganga Institute of Technology, Tumakuru -572103, Karnataka.