key: cord-0724381-xoeokhrd
authors: Kumar, Vipul; Dhanjal, Jaspreet Kaur; Bhargava, Priyanshu; Kaul, Ashish; Wang, Jia; Zhang, Huayue; Kaul, Sunil C.; Wadhwa, Renu; Sundar, Durai
title: Withanone and Withaferin-A are predicted to interact with transmembrane protease serine 2 (TMPRSS2) and block entry of SARS-CoV-2 into cells
date: 2020-06-16
journal: J Biomol Struct Dyn
DOI: 10.1080/07391102.2020.1775704
sha: 6c515a645a55bd4ee1ee9471317289cd2aab959b
doc_id: 724381
cord_uid: xoeokhrd

Coronavirus disease 2019 (COVID-19) initiated in December 2019 in Wuhan, China and became pandemic causing high fatality and disrupted normal life calling world almost to a halt. Causative agent is a novel coronavirus called Severe Acute Respiratory Syndrome Coronavirus 2 (SARS-CoV-2/2019-nCoV). While new line of drug/vaccine development has been initiated world-wide, in the current scenario of high infected numbers, severity of the disease and high morbidity, repurposing of the existing drugs is heavily explored. Here, we used a homology-based structural model of transmembrane protease serine 2 (TMPRSS2), a cell surface receptor, required for entry of virus to the target host cell. Using the strengths of molecular docking and molecular dynamics simulations, we examined the binding potential of Withaferin-A (Wi-A), Withanone (Wi-N) and caffeic acid phenethyl ester to TPMRSS2 in comparison to its known inhibitor, Camostat mesylate. We found that both Wi-A and Wi-N could bind and stably interact at the catalytic site of TMPRSS2. Wi-N showed stronger interactions with TMPRSS2 catalytic residues than Wi-A and was also able to induce changes in its allosteric site. Furthermore, we investigated the effect of Wi-N on TMPRSS2 expression in MCF7 cells and found remarkable downregulation of TMPRSS2 mRNA in treated cells predicting dual action of Wi-N to block SARS-CoV-2 entry into the host cells. Since the natural compounds are easily available/affordable, they may even offer a timely therapeutic/preventive value for the management of SARS-CoV-2 pandemic. We also report that Wi-A/Wi-N content varies in different parts of Ashwagandha and warrants careful attention for their use. Communicated by Ramaswamy H. Sarma

Coronaviruses (family Coronaviridae; order Nidovirales) are large, enveloped, positive-stranded RNA viruses. Their genome, largest among all RNA viruses, is packed inside a helical capsid constituted of nucleocapsid protein (N), surrounded by an envelope comprised of three structural proteins; (i) the membrane protein-M, (ii) the envelope protein-E and (iii) the spike protein-S. Whereas M and E proteins are involved in virus assembly, the S protein (forms large protrusions from the virus surface, giving it an appearance of crown) mediates virus entry into target host cells and serves as a critical determinant of viral host range (Boopathi et al., 2020; Li, 2016) . Some coronaviruses also encode an envelope-associated hemagglutinin-esterase protein (HA). Among the four Coronaviruses genera (a, b, c and d), a and b coronaviruses infect mammals, c infect avian species and d infect both mammalian and avian species, and pose serious health threats (Li, 2016) . Severe acute respiratory syndrome coronavirus (SARS-CoV) and Middle East respiratory syndrome coronavirus (MERS-CoV) when transmitted from animals to humans, caused severe respiratory diseases (SARS and MERS in 2002 and 2012 with fatality rates of $10% and 36%, respectively). Porcine epidemic diarrhea coronavirus caused $100% fatality rate in piglets in 2013 in USA. A novel coronavirus, SARS-coronavirus 2 (SARS-CoV-2), closely related to SARS-CoV, has recently been causing infectious respiratory disease, coronavirus disease 19 . Originated and rapidly spreading (by human to human transmission) from Wuhan city of China to all over the world . The viral protein is synthesized as a precursor protein.

Its cleavage by host cell proteases including transmembrane protease serine 2 (TMPRSS2), endosomal cathepsin L protease or human airway trypsin-like protease (HAT) is an essential step for virus fusion to the host cell membrane, infectivity and propagation (B€ ottcher-Friebertsh€ auser et al., 2010; Yamaya et al., 2015) . Based on these information, protease inhibitors are considered as drugs for treatment/prevention of viral infections. Indeed, overexpression of TMPRSS2 and HAT was shown to promote the growth of different subtypes of human and avian influenza viruses (B€ ottcher et al., 2009; Shen et al., 2017) . On the other hand, serine protease inhibitors abrogated the virus proliferation (B€ ottcher et al., 2009; B€ ottcher-Friebertsh€ auser et al., 2011) . Peptide-conjugated phosphorodiamidate morpholino oligomers (PPMO) that sterically blocked TMPRSS2 transcription caused significant suppression of viral titers (B€ ottcher-Friebertsh€ auser et al., 2011) . Benzylsulfonyl-d-arginine-proline-4-amidinobenzylamide, a potent inhibitor of HAT and TMPRSS2 was also shown to cause strong suppression of virus propagation (B€ ottcher-Friebertsh€ auser et al., 2012 (B€ ottcher-Friebertsh€ auser et al., , 2013 . Camostat mesylate (an inhibitor of serine protease, clinically used to treat chronic pancreatitis) has also been suggested as a candidate antiviral drug to prevent or suppress TMPRSS2-dependent infection by SARS-CoV (Kawase et al., 2012) . Cells constitutively expressing TMPRSS2 showed higher susceptibility and larger syncytia when infected with MERS-CoV. Camostat mesylate although completely inhibits the syncytium formation, could only partially block virus entry into cells. The simultaneous treatment with inhibitors of cathepsin L and TMPRSS2 was shown to completely block virus entry into Vero-TMPRSS2 cells (Shirato et al., 2013) . Cysteine protease inhibitor K11777 (used for the treatment of tropical parasitic disease, American trypanosomiasis) and its closely related vinylsulfones have been shown to target cathepsin-mediated cell entry and act as broad-spectrum antivirals. It was shown to inhibit SARS-CoV and Ebola virus entry into the host cells. Activity of K11777, its derivative and Camostat mesylate were compared in pathogenic animal model of SARS-CoV infection. It was demonstrated that the viral pathogenesis of SARS-CoV is dependent on serine rather than cysteine proteases and hence Camostat mesylate was recommended as a potential therapeutic intervention for respiratory coronavirus infections (Yamaya et al., 2015; Zhou et al., 2015) . TMPRSS2 was recently proven to be crucial for hemagglutinin cleavage of some human influenza viruses. Since the catalytic sites of the diverse serine proteases linked to influenza, parainfluenza and coronavirus activation are structurally similar, active site inhibitors of these airway proteases were predicted to possess broad therapeutic applicability against multiple respiratory viruses (Limburg et al., 2019) . Alternatively, superior selectivity could be achieved with allosteric inhibitors of TMPRSS2 or another critical protease (Laporte & Naesens, 2017) . Recently, SARS-CoV-2 was shown to use the ACE2 receptor for entry and the serine protease TMPRSS2 for S protein priming. Camostat mesylate blocked the entry of SARS-CoV-2 virus into the cells (Hoffmann et al., 2020) .

Genomic characterization of the SARS-CoV-2, its variance, evolution, transmission dynamics (human to human transmission by droplets or direct contact) have also recently been reported (Benvenuto et al., 2020; Ceraolo & Giorgi, 2020; Lu et al., 2020; Paraskevis et al., 2020) . It is an enveloped single-stranded RNA b-coronavirus similar to the SARS and MERS viruses and closely related to two bat-derived SARS-like coronaviruses (CoVZC45 and bat-SL-CoVZXC21; Lai et al., 2020) . Lacking any specific treatment modalities for the new virus, the current treatment options are tailored from earlier SARS-CoV and MERS-CoV experiences and repurposing of existing drugs. These include broad-spectrum antibiotics, corticosteroids, interferons, ribavirin, lopinavir-ritonavir or mycophenolate mofetil and have not been subjected to clinical trials as yet. Shirato et al. (2017) examined the mechanism of cell entry used by a pseudo typed virus bearing the HCoV-229E spike (S) protein in the presence or absence of protease inhibitors. They showed that compared with a laboratory strain isolated in 1966, the recent clinical isolates of HCoV-229E were less likely to use endosomal protease cathepsin L due to two viral mutations (amino acid substitutions, R642M and N714K) that altered their sensitivity to a cathepsin L inhibitor. They suggested that the cell surface TMPRSS2 is a better target than endosomal cathepsin for therapeutic agents for SARS and MERS (Shirato et al., 2017) . Explorations on repurposing of drugs continues worldwide to find immediate therapeutic strategies for deadly SARS-CoV-2. Computational approaches for quick screening of available compounds is gaining attention among the scientific community for the same. In a recent study, Indinavir and Remdesivir were investigated for their docking potential to SARS-CoV-2 protease. Though no binding was observed within the catalytic site of the protease, the active form (CHEMBL2016761) of Remdesivir could dock in the overlapping region of the NTP binding motif urging its validation through clinical trials for COVID-19 infection (Chang et al., 2020) . Remdesivir has also been tested for its potential to bind to RNA dependent RNA polymerase of SARS-CoV-2 (Nejadi Babadaei, Hasan, Vahdani, et al., 2020) . The U.S. Food and Drug Administration (FDA) has now issued an emergency use authorization for Remdesivir for the treatment of suspected or laboratory-confirmed adults and children hospitalized with severe symptoms of SARS-CoV-2 disease (FDA news release, 2020; Hendaus, 2020). In line with this, many research groups across the globe are exploring various drugs for potential repurposing to find an immediate and effective therapeutic candidate for COVID-19 pandemic, by mainly targeting SARS-CoV-2 main protease (Aanouz et al., 2020; Al-Khafaji et al., 2020; Das et al., 2020; Enmozhi et al., 2020; Gyebi et al., 2020; Islam et al., 2020; Joshi et al., 2020; Khan, Jha, et al., 2020; Khan, Zia, et al., 2020; Kumar, Kumari, et al., 2020; Lobo-Galo et al., 2020; Muralidharan et al., 2020; Pant et al., 2020; Umesh et al., 2020) , E and N protein (Gupta et al., 2020; Sarma et al., 2020) , RNA-dependent RNA polymerase (Elfiky, 2020) , NSP15 and prefusion 2019-nCoV spike glycoprotein (Sinha et al., 2020) and cell surface receptors (ACE-2 and TMPRSS2) of host cells (Abdelli et al., 2020; Elmezayen et al., 2020; Thuy et al., 2020; Wahedi et al., 2020) .

Indian Ayurvedic herb, Ashwagandha (Withania sominifera) and honey bee propolis have been heavily used in traditional home medicine systems. Both of these are trusted to boost the immune function, possess a variety of prophylactic and therapeutic activities (Latheef et al., 2017; Munagala et al., 2011) . Several studies have provided evidence that honeybee propolis inhibits variety of viruses including herpes simplex virus, sindbis virus, parainfuenze-3 virus, human cytomegalovirus, dengue virus type-2, influenza virus A1 and rhinovirus (Amoros et al., 1994; Kwon et al., 2019; Munagala et al., 2011; Serkedjieva et al., 1992) . Withaferin-A (Wi-A), a bioactive withanolide from Ashwagandha, was shown to possess inhibitory activity for HPV and influenza viruses (Cai et al., 2015) . In light of these literature, we examined the potential of three natural compounds, i.e. Wi-A and Withanone (Wi-N) derived from Ashwagandha, and caffeic acid phenethyl ester (CAPE) from honey bee propolis for their ability to dock and inhibit TMPRSS2, while comparing with Camostat mesylate by molecular docking and molecular dynamics (MD) simulations analyses. We found that both Wi-A and Wi-N could bind and stably interact at the catalytic site of TMPRSS2. Wi-N showed stronger interactions with TMPRSS2 catalytic residues than Wi-A and was also able to induce changes in the allosteric site of TMPRSS2 similar to the ones reported for Camostat mesylate. We also examined if Wi-N could affect the endogenous expression of TMPRSS2 and found its transcriptional downregulation in Wi-N treated cells. The data predicts dual action of Wi-N to block SARS-CoV-2.

The 3D structure of TMPRSS2 is not yet available in the Protein Data Bank, however, a structure modeled using a homologous protein, serine protease hepsin (5ce1.1), is available in the Swiss model repository (O15393). The identity between the sequence of TMPRSS2 and hepsin is 33.82%, and the Q mean of the modeled structure was -1.62. Also, the catalytic domain in both the proteins was found to be well conserved with identical catalytic residues His296, Asp345 and Ser441 (Afar et al., 2001) . This modeled structure was further prepared for docking using the protein preparation wizard of the maestro from Schrodinger suite (Schr€ odinger, 2020). The preparation steps mainly included the addition of missing disulfide and hydrogen bonds, filling of missing amino acids side chains and optimization of hydrogen bonds. Then OPLS3e forcefield was used for restrained minimization until the average root mean square deviation (RMSD) of the non-hydrogen atoms converged to 0.30 Å (Harder et al., 2016) . The structure of Camostat mesylate (CID 5284360), Wi-A (CID 265237), Wi-N (CID 21679027) and CAPE (CID 5281787) was downloaded from the PubChem database. The structure of all the ligands shown in Figure 1 was prepared using the LigPrep module of the Schrodinger suite (Schr€ odinger, 2020).

For the ligand preparation, OPLS3e force field was used for energy minimization and a pH of 7.0 ± 2.0 was set to generate all the possible ionization states. Following this, the ligands were desalted, their tautomers were generated, while retaining the chiralities, possible stereoisomers were also generated (Sastry et al., 2013) .

The grid was generated at the catalytic site of TMPRSS2, mainly covering His296, Asp345 and Ser441 residues for performing molecular docking studies. Extra precision flexible docking protocol was used for all the dockings (Friesner et al., 2006 ).

Desmond package with OPLS3e force field from Schrodinger was used to study the stability of the protein-ligand system in the presence of explicit water molecules (Schr€ odinger, 2020). Each protein-ligand complex was solvated with the TIP4P water model in an orthorhombic periodic boundary box. To prevent interaction of the protein complex with its own periodic image, the distance between the complex and the box wall was kept 10 Å. The system was then neutralized by the addition of ions. Energy of the prepared systems was minimized for 5000 steps using the steepest descent method or until a gradient threshold of 25 kcal/mol/Å was achieved. It was followed by L-BFGS (Low-memory Broyden-Fletcher-Goldfarb Shanno quasi-Newtonian minimizer) until a convergence threshold of 1 kcal/mol/Å was met (Harder et al., 2016 ). The equilibrated system was then subjected to 100 ns simulation in NPT ensemble with 300 K temperature, constant pressure of 1 atm and a time step of 2 fs.

The MD trajectories were analyzed post MD simulations using the Desmond Simulation event analysis tool. RMSD of conformations acquired by all protein-ligand complexes during the 100-ns simulation run with respect to their docked pose was analyzed and compared with the behavior of Camostat mesylate-TMPRSS2 complex as a function of time to investigate the stability of ligands within the active site pocket of TMPRSS2. Root mean square fluctuation (RMSF) of TMPRSS2 residues was also analyzed to investigate the flexibility of the protein region interacting with the ligands. The interaction pattern of all three ligands-Wi-A, Wi-N and CAPE within the catalytic domain of TMPRSS2 was studied and compared with the interactions of Camostat mesylate with TMPRSS2. Various residues residing in the catalytic domain of the protein were investigated for their time of contact with the different ligands to calculate their occupancy percentage for interaction during the simulation time. The radius of gyration and RMSD of ligands along the simulation trajectory was calculated and analyzed.

A total of 100 frames representing the protein-ligand conformations spanned in between 40 and 100 ns of simulation run were extracted to calculate the MM/GBSA (molecular mechanics energies combined with the 906generalized Born and surface area continuum solvation) free binding energy using the Prime module of Schrodinger suite (Schr€ odinger, 2020).

The equation used for the calculation is:

DG complex , G receptor and G ligand represent the free energies of the complex, receptor and ligand, respectively.

2.6. Cell culture and viability assays MCF7 (human breast carcinoma known to express TMPRSS2) cells were bought from The Japanese Collection of Research Bioresources Cell Bank (Tokyo, Japan) and maintained in Dulbecco's Modified Eagle's Medium (DMEM) (Invitrogen, Carlsbad, California, USA) supplemented with 10% fetal bovine serum and 1% penicillin/streptomycin in a humidified incubator (37 C and 5% CO 2 ). Five thousand cells/well were plated in 96-well plate and allowed to settle overnight, followed by treatment with varying doses of Wi-N. Control cells were treated with dimethyl sulfoxide (DMSO). After 48 h, 10 lL of MTT f3-(4,5-dimethylthiazol-2-yl)-2,5-diphenyltetrazolium bromideg (Sigma-Aldrich, M2003-1G) in DMEM was added to each well and incubated for 4 h. The media and MTT from the wells were aspirated out and replaced with 100% DMSO, followed by measurement of optical density at 570 nm using Tecan infinite M200V R Pro using a microplate reader (Infinite 200 PRO; Tecan Group Ltd., Mannedorf, Switzerland). For Quantitative Crystal Voilet (QCV) staining, after 48 h treatment, cells were fixed, stained with CV and incubated with de-staining solution (40% methanol:10% glacial acetic acid:50% water) for 5-10 mins. The optical density of de-staining solution was measured by Tecan infinite M200V R Pro using a microplate reader at 570 nm. The data represents mean ± SD from three experiments. Statistical significance was calculated by Unpaired t test (GraphPad Prism, GraphPad Software, San Diego, CA; ÃÃ represents p value < 0.05 for significant change). Phase contrast images of control and treated cells were captured on bright field microscopy (10X) Nikon ECLIPSE TE300 (Nikon Instruments Inc. Tokyo, Japan).

MCF7 cells were treated with non-toxic (determined by independent experiments) dose (20-40 lM) of Wi-N for 48 h. Control and Wi-N treated cells were harvested by trypsin-EDTA (Wako, Japan). Total RNA was extracted using TRIzol TM (Thermo Fisher Sci., Japan). cDNA synthesis was performed using QuantitTect Reverse Transcription Kit (Qiagen, Germany) followed by real-time RT-PCR on the ECO TM real time system (Illumina, San Diego, CA, USA) using SYBRV R Green Select Master Mix (Applied Biosystem, Japan). The expression was normalized against 18S as an internal control. Single PCR product amplification was confirmed by meltcurve analysis 2 -DDCT . Primer sequences used in this study were -TMPRSS2 [forward] 5 0 -GAGGACGAGAATCGGTGTGT-3 0 and [reverse] 5 0 -TCCAGTCGTCTTGGCACA-3 0 (Foulds et al., 2010) and internal control 18S

[forward] 5 0 -CAGGGTTCGATTCCGTAGAG-3 0 and [reverse] 5 0 -CCTCCAGTGGATCCTCGTTA-3 0 . The data represents mean ± SD from three independent experiments. Statistical significance was calculated by unpaired t test (GraphPad Prism, GraphPad Software, San Diego, CA; p values ÃÃÃ < 0.001 represent highly significant change).

Dry powder of roots, stems and leaves from organically raised Ashwagandha were extracted in b-cyclodextrin-water and the amount of Wi-N and Wi-A was evaluated in each by HPLC, as described earlier . The analyses were performed using Develosil (150 mm L Â 4.6 mm i.d. 5 lm) C30-UG column (Nomura Chemical Co. Ltd., Japan) on Shimadzu HPLC CBM-20A/20Alite (Shimadzu Corporation, Tokyo, Japan) equipped with UV detector (SPD-20A/20AV) and Lab solutions Software. The separation was carried out using a linear gradient elution program with gradient 

3.1. TMPRSS2, an attractive therapeutic target for COVID-19

TMPRSS2, a single-pass type 2 membrane protein, has been reported as the activator of S protein of various strains of SARS-CoV (Belouzard et al., 2012; Bertram et al., 2011) . A very recent report has confirmed that SARS-CoV-2 entry to human cells also depends on the S protein priming by TMPRSS2 (Hoffmann et al., 2020) . It is also reported that SARS-CoV-2 can use the endosomal cysteine proteases CatB/L for S priming in TMPRSS2cells, however, TMPRSS2 but not CatB/L is essential for the virus entry and its spread (Iwata-Yoshikawa et al., 2019). The main catalytic residues of TMPRSS2 involved in the proteolytic activity are His296, Asp345 and Ser441 (Paoloni-Giacobino et al., 1997; Wilson et al., 2005) . Further, TMPRSS2 has been shown to be dispensable in mice knockout model that showed no significant effect on development, survival or normal organ structure or functions, suggesting its functional redundancy with other transmembrane serine proteases or its non-vital function in case of stressed or disease scenario (Kim et al., 2006) . In these premises, TMPRSS2 has emerged as an attractive therapeutic target against COVID-19. Camostat mesylate (a serine protease inhibitor), well documented for the treatment of cancer, pancreatitis and liver fibrosis, has recently been shown to inhibit TMPRSS2 and entry of SARS-CoV-2 virus into the cells (Hoffmann et al., 2020) . In earlier studies, it was shown to significantly reduce H1N1 influenza virus titers in the culture supernatant (Yamaya et al., 2015) and protect mice from H1N1 influenza virus infection (Bahgat et al., 2011; Lee et al., 1996) . It was also shown to inhibit TMPRSS2 and yield 10-fold reduction in SARS-CoV virus titre (Shulla et al., 2011) and protected mice from the infection thus increasing the survival rate to 60% (Zhou et al., 2015) . Its action against MERS-CoV has also been reported. Concentration of 10 lM could reduce the entry of virus in Vero-TMPRSS2 cells by 15folds. Two hundred and seventy-fold decrease in the amount of viral RNA in the culture supernatant of Calu-3 cells was also observed three days post MERS-CoV infection with 100 lM dose of Camostat (Shirato et al., 2017) . Camostat mesylate, similar to the other non-specific serine protease inhibitors, targets a highly conserved catalytic domain comprising of His-Asp-Ser triad. The protease catalytic domain of TMPRSS2 is constituted of amino acid residues His296, Asp345 and Ser441, corresponding to His57, Asp102 and Ser195 of chymotrypsinogen, another well-known serine protease (Lin et al., 1999; Park, 2011) . Although the exact molecular mechanism of action of camostat remains elusive, to be effective against different serine proteases, it must target the active site harboring the conserved catalytic triad residues. Based on the anti-viral activity of Camostat mesylate, we decided to use it as a reference molecule. We first modelled its interactions with the catalytic residues of TMPRSS2 and then compared them with our small molecules; Wi-A, Wi-N and CAPE. The crystal structure of TMPRSS2 has not been reported yet. Therefore, we used a homology-based structural model from Swiss model repository (O15393). This structure was modeled using the structural coordinates of human hepsin as a template. Both TMPRSS2 and hepsin belong to the same family of the proteases i.e. type-II membrane-associated serine peptidase (Afar et al., 2001) . Both Hepsin and TMPRSS2 are reported to share the same arrangement of the catalytic triad amino acids -His296, Asp345 and Ser441 (Afar et al., 2001) . Based on this assessment of model confidence, we further processed this structure for carrying out molecular docking and MD simulations through restrained minimization. It allowed hydrogen atoms to be freely minimized. It also allowed sufficient movement of heavy atom to relax strained bonds, angles and clashes.

In these premises, we first investigated the docking score for Camostat mesylate, a known inhibitor, against the substrate- binding pocket of TMPRSS2. The docking score for the best pose of TMPRSS2/Camostat mesylate complex was found to be -5.90 Kcal/mol. Gly464 was mainly contributing to this by making a hydrogen bond with the ligand. Wi-A, Wi-N and CAPE were docked with TMPRSS2 using the same grid. All the ligands were binding at the catalytic site of the protein. For TMPRSS2/CAPE complex, the best pose had the binding energy of -6.20 Kcal/mol, and Gly464 and Ser436 were involved in the hydrogen bond interactions. Complex of TMPRSS2 with Wi-A had the docking score of -5.60 Kcal/mol. Here, Glu299 and Lys342 were mainly involved in the polar interactions. And in case of Wi-N, the docking score was -4.30 Kcal/mol and Gly462 of TMPRSS2 was involved in the polar interactions. The detailed list of residues involved in polar and non-polar interactions with all the four ligands -Camostat mesylate, Wi-A, Wi-N and CAPE is shown in Table 1 . As molecular docking alone is not sufficient to give the real estimate of binding pose and binding affinity, the best-docked pose for each complex was further subjected to MD simulations and MM/GBSA free binding energy calculations for better understanding of the molecular basis of interactions.

To investigate the dynamic behavior of each protein-ligand complex and account for the stability of the interactions we performed 100 ns long MD simulations in water solvated state. Visual analysis of the trajectory showed that TMPRSS2 could not retain CAPE in its substrate binding pocket for long. It lost its interaction and moved out within 37 ns of the simulation run. Wi-A and Wi-N, on the other hand, were found to be stably interacting at the docked site throughout the 100 ns of simulations and was very similar to Camostat mesylate. The overall conformation of TMPRSS2-Camostat mesylate, TMPRSS2-Wi-A and TMPRSS2-Wi-N were found stable and did not show much fluctuations (Figure 2(A) ). The RMSF plot of the complexes was found to be aligned with TMPRSS2(Apo) residues, and no significant fluctuations were observed in the catalytic triad of TMPRSS2 (Figure 2(B) ). To get further insights of the conformational changes, the simulated complexes were superimposed on TMPRSS2(Apo). The RMSD of TMPRSS2-Camostat mesylate with respect to TMPRSS2(Apo) was estimated to be 1.76 Å, for TMPRSS2-Wi-A it was 2.0 Å and for TMPRSS2-Wi-N it was 2.4 Å. In further scrutiny, it was found that the loop region from Gly245 to Gly265 had an outward conformation in the natural compound complexes, whereas it had an inward orientation in Camostat mesylate complex, when compared to TMPRSS2(Apo) structure. Further, it was found that Wi-N and Camostat mesylate were inducing the same conformation changes in TMPRSS2. Both these molecules forced the region from Arg316 to Tyr322 to acquire a helical secondary structure, while in case of Wi-A and TMPRSS2(Apo) this region was found to be a loop ( Figure  2(C) ). Furthermore, as all the ligands were binding within the same pocket of TMPRSS2, interacting with almost the same residues, suggested a similar mechanism of their action (Figure 2(D) ). The stable binding of the ligands and the compactness of Wi-A and Wi-N were also seen through RMSD and radius of gyration calculations. All the calculations with the average values and standard deviation throughout the 100 ns of the simulations are shown in Table 2. 3.4. Wi-N binding to Ser441 of TMPRSS2 is likely to disrupt substrate binding at the catalytic site Hydrogen bonds predominantly contribute to the specificity of molecular recognition and are therefore exploited the most for the design of effective drug molecules. In view of this, we further investigated the hydrogen bond interactions made by the natural compounds during the simulation run. Hydrogen bonds were calculated taking the criteria that the distance between hydrogen bond donor and acceptor should be less than 2.5 Å, the donor angle should be less than or equal to 120 degrees between donor-hydrogen acceptor atoms and the acceptor angle should be less than or equal to 90 degrees between the hydrogen acceptorbonded atoms. It was found that Wi-N was making significant polar contact (for more than 50% of the simulation time) with Ser441, a residue directly involved in catalytic activity of TMPRSS2, whereas Wi-A was making water-mediated contact with Ser 441 (Figure 3(A,B) ). In both Wi-A and Wi-N, this interaction with the catalytic residues was maintained by the lactone ring. Though Camostat mesylate had numerous polar and non-polar contacts within the binding pocket of TMPRSS2, it was not making polar interaction with any of the catalytic triad residues (His296, Asp435, Ser441; Figure 3 (C)). Overall analyses showed that at any fraction of Val275, Gln276, Val278, His279, Val280, Glu 299, Lys300, Pro301, Leu302, Tyr337, Asp345, Glu389, Lys390, Gly391, Lys392, Thr393, Asp435, Ser436, Cys437, Asp440, Asp458, Thr459, Ser460, Trp461, Ser463, Arg470, Pro471, Val473, Ser436, Cys437, Ser441, Thr459, Ser460, Trp461, Gly462, Gly464, Tyr474 Ser206, His296, Glu299, Lys342, Glu388, Glu389, Lys390, Gly391, Asp435, Ser436, Cys437, Gln438, Gly439, Asp440, Ser441, Thr459, Ser460, Trp461, Gly464, Lys467, Arg470, Pro471, Val473, Gly439, Asp440, Ser441, Gly462, Ser463, Gly464 His296, Glu299, Lys342, Asp435, Ile381, Ala386, Glu389, Lys390, Val434, Ser436, Cys437, Gln438, Gly439, Asp440, Ser441, Thr459, Ser460, Trp461, Gly462, Ser463, Gly464, Cys465, Ala466, Lys467, Tyr474 time, all the ligands had a similar number of polar and nonpolar interactions with TMPRSS2 while in terms of consistent interactions with catalytic residues of TMPRSS2, natural molecules were better. A detailed analysis of the interacting residues during the course of MD simulations has been shown in Table 3 . The most significant interactions between the protein-ligand complexes are illustrated in Figure 4 (A-C). This analysis taking Camostat mesylate as a reference molecule thus suggests that Wi-N and Wi-A have a good affinity towards the substrate-binding pocket of TMPRSS2 and could probably be a natural and readily available drug for the inhibition of S protein priming of SARS-Cov2.

3.5. Wi-N shows dual inhibitory activity towards TMPRSS2-inhibition of its catalytic activity and transcriptional downregulation MM/GBSA is a method commonly used or estimating ligandbinding affinities in protein systems. We used 100 snapshots from the conformational space of each protein-ligand complex in between 40 to 100 ns for calculating the MMGBSA free binding energy. The average MMGBSA free binding energy of Camostat mesylate, Wi-N and Wi-A with TMPRSS2 was estimated to be -65.20 ± 1.64 Kcal/mol, À46.80 ± 5.13 Kcal/mol and -37.80 ± 5.80 Kcal/mol, respectively. Though the free binding energy of Camostat mesylate was more in comparison to Wi-N and Wi-A, the free binding energy of Wi-N was in range of -50 Kcal/mol clearly indicating equitable affinity of this natural molecule towards TMPRSS2. Camostat mesylate, an established inhibitor of TMPRSS2, has been shown to cause 30%-40% reduction in SARS infection (virus titer; Kawase et al., 2012) . It was recently demonstrated that SARS-CoV-2 follows the mechanism similar to SARS for invasion into host cell and requires TMPRSS2 for its protein priming. When Vero cells expressing TMPRSS2 were treated with Camostat, the infection by SARS-CoV-2 was reduced to 40%-50% (Hoffmann et al., 2020) . However, the mechanism of inhibition of TMPRSS2 by Camostat mesylate, per se, has not been clearly elucidated. Azouz et al. (2020) showed decrease in TMPRSS2 protein in Camostat mesylate-treated cells. On the other hand, Yamaya et al. (2015) reported that the expression levels of TMPRSS2 mRNA remained unaffected by Camostat mesylate.

Downregulation of endogenous TMPRSS2 expression clearly holds promise as an approach towards preventing transmission of coronavirus. Earlier, peptide-conjugated PPMO, single-stranded-DNA-like antisense agents that readily enter cells and sterically blocking cRNA, targeted to TMPRSS2 caused significant suppression of viral titers (B€ ottcher-Friebertsh€ auser et al., 2011) . Most recently, Shen et al. (2020) reported an inhibition of influenza A virus propagation by benzoselenoxanthenes by downregulation of TMPRSS2 expression. We next considered and examined if Wi-N, in addition to the predicted inhibition of TMPRSS2, have such capability. We used human breast carcinoma (MCF7) cells, known to express TMPRSS2. Non-toxic doses of Wi-N were first determined by MTT and QCV assays. Expression of TMPRSS2 was then determined in control and Wi-N treated cells by real time RT-PCR using specific primers for TMPRSS2 and 18S (as an internal control). As shown in Figure 5 , 40%-50% reduction in TMPRSS2 transcript was observed in cells treated with non-toxic doses of Wi-N. These data suggested that Wi-N may have dual impact on TMPRSS2, i.e. inhibition of its catalytic activity as predicted by molecular modelling data and transcriptional downregulation as shown in Figure 5 . Both of these may collectively block the cleavage of viral protein and entry of the virus to host cells, effectively.

In light of the above findings, we considered to examine the best source of Wi-N. The alcoholic extract of Ashwagandha leaves and Wi-N have earlier been shown to possess anticancer potential, and safe for normal cells (Widodo et al., 2008) . The normal cells were shown to be protected from oxidative, DNA damage, chemical and amnesia stresses in cell culture and mice models of study, respectively (Gao et al., 2014; Konar et al., 2011; Wadhwa et al., 2016; Widodo et al., 2009) . It was earlier reported that Ashwagandha leaves possess higher amount of Wi-A and Wi-N, as compared to the roots . Although there was a high ratio of Wi-N/Wi-A in roots, the total content was extremely (10 times) low as compared in leaves . Methods of water-based extraction (using cyclodextrins) for these bioactive components, and hydroponic cultivation yielding leaves with enriched level of Wi-N have also been established . The use of vermicompost as an organic fertilizer in Ashwagandha cultivation was seen to yield not only the higher leaf and root mass, but also higher level of Wi-N and Wi-A (Kaur et al., 2018) . In view of the current findings that Wi-A and Wi-N may be particularly useful for the management of COVID-19 disease, we investigated their amount in the synchronous Ashwagandha sources and their extracts. Root, stem and leaves from the same Ashwagandha harvest were carefully separated, dried, powdered, extracted and accessed for Wi-A and Wi-N contents (Foulds et al., 2010) . In multiple replicates, we found that b-CD assisted water extract from three parts of the plant Figure 6 . Content of Wi-A and Wi-N in root, stem and leaves of Ashwagandha as detected by high pressure liquid chromatography. Dry powders of each of the sample obtained from the same plant/harvest of Ashwagandha were extracted with b-CD water . Content of Wi-A and Wi-N was determined by HPLC with respect to the standards shown on the top. X-axis and Y-axis show minutes and absorbance units, respectively. contained both Wi-A and Wi-N. Interestingly, whereas the amount of Wi-A was highest in leaves, Wi-N was enriched in stem ( Figure 6 ). Recently, Wi-N has also been predicted to possess the ability to interact with and inhibit the main protease of SARS-CoV-2 (Kumar, Dhanjal, et al., 2020) , and may hence act as two-way sword, namely (i) inactivating the virus and (ii) blocking its entry to the target cells. Interestingly, an independent bioinformatics screening of thousands of phytochemicals selected Wi-N as a potential inhibition of ACE2, another cell surface receptor protein required for entry of virus to the target cells (Balkrishna et al., 2020) . Taken together, these findings propose that (i) Wi-N could be a strong candidate for the management of COVID-19 and (ii) Ashwagandha stem and its Wi-N rich extract may provide preventive and therapeutic benefits in the current pandemic.

The entry of SARS-CoV-2 into the host depends on the interaction with ACE-2 receptor of the host cell surface. After the successful attachment with ACE2 receptor, the S protein of SARS-CoV-2 needs to be cleaved and activated for its fusion into the host cell membrane, facilitated by the TMPRSS-2 enzyme of the host (Hoffmann et al., 2020) . In the present study, the inhibitory potential of three natural compounds, Wi-A, Wi-N and CAPE for TMPRSS2 in comparison to its known inhibitor Camostat mesylate was investigated by molecular modeling tools. We demonstrate that both Wi-A and Wi-N could bind and stably interact at the catalytic site of TMPRSS2. The interactions of Wi-N with catalytic residues were stronger and significant than Wi-A. Wi-N also able to induce changes in the allosteric site of TMPRSS2 similar to the ones reported for Camostat mesylate. Furthermore, cells treated with Wi-N showed remarkable downregulation of TMPRSS2 suggesting dual potential of Wi-N to block TMPRSS2 and hence the entry of SARS-CoV-2 to host cells. While these preliminary data necessitate further in vitro and in vivo experimental evidence, we demonstrate that Wi-N/Wi-A (i) vary in their content in different parts of Ashwagandha plant and hence their careful recruitment may be helpful in the management of current outbreak of fatal SARS-CoV-2 and (ii) may provide a reliable direction to new drug development.

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

In silico study the inhibition of Angiotensin converting enzyme 2 receptor of COVID-19 by Ammoides verticillata components harvested from western Algeria

Catalytic cleavage of the androgen-regulated TMPRSS2 protease results in its secretion by prostate and prostate cancer epithelia

Using integrated computational approaches to identify safe and rapid treatment for SARS -CoV-2

Comparison of the anti-herpes simplex virus activities of propolis and 3-methyl-but-2-enyl caffeate

Alpha 1 antitrypsin is an inhibitor of the SARS-CoV2-priming protease TMPRSS2. bioRxiv

Inhibition of lung serine proteases in mice: A potentially new approach to control influenza infection

Withanone from Withania somnifera may inhibit novel coronavirus (COVID-19) entry by disrupting interactions between viral S-protein receptor binding domain and host ACE2 receptor

Mechanisms of coronavirus cell entry mediated by the viral spike protein

The 2019-new coronavirus epidemic: Evidence for virus evolution

Cleavage and activation of the severe acute respiratory syndrome coronavirus spike protein by human airway trypsin-like protease

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

MDCK cells that express proteases TMPRSS2 and HAT provide a cell system to propagate influenza viruses in the absence of trypsin and to study cleavage of HA and its inhibition

Cleavage of influenza virus hemagglutinin by airway proteases TMPRSS2 and HAT differs in subcellular localization and susceptibility to protease inhibitors

Activation of influenza viruses by proteases from host cells and bacteria in the human airway epithelium. Pathogens and Disease

Hemagglutinin activating host cell proteases provide promising drug targets for the treatment of influenza A and B virus infections. Vaccine

Inhibition of influenza virus infection in human airway cell cultures by an antisense peptide-conjugated morpholino oligomer targeting the hemagglutinin-activating protease TMPRSS2

Promising anti-influenza properties of active constituent of Withania somnifera ayurvedic herb in targeting neuraminidase of H1N1 influenza: Computational study

Genomic variance of the 2019-nCoV coronavirus

Potential therapeutic agents for COVID-19 based on the analysis of protease and RNA polymerase potential therapeutic agents for COVID-19 based on the analysis of protease and RNA polymerase docking

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

SARS-CoV-2 RNA dependent RNA polymerase (RdRp) targeting: An in silico perspective

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

Coronavirus (COVID-19) update: FDA Issues emergency use authorization for potential COVID-19 treatment

Research resource: Expression profiling reveals unexpected targets and functions of the human steroid receptor RNA activator (SRA) gene

Extra precision glide: Docking and scoring incorporating a model of hydrophobic enclosure for protein-ligand complexes

Withanone-rich combination of Ashwagandha withanolides restricts metastasis and angiogenesis through hnRNP-K

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

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

OPLS3: A force field providing broad coverage of drug-like small molecules and proteins

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

Remdesivir in the treatment of Coronavirus Disease 2019 (COVID-19): A simplified summary

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

TMPRSS2 contributes to virus spread and immunopathology in the airways of murine models after coronavirus infection

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

Novel methods to generate active ingredients-enriched Ashwagandha leaves and extracts

Organic cultivation of Ashwagandha with improved biomass and high content of active Withanolides: Use of

Simultaneous treatment of human bronchial epithelial cells with serine and cysteine protease inhibitors prevents severe acute respiratory syndrome coronavirus entry

Targeting SARS-CoV-2: A systematic drug repurposing approach to identify promising inhibitors against 3C-like proteinase and 2'-O-ribose methyltransferase

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

Phenotypic analysis of mice lacking the Tmprss2-encoded protease

Protective role of Ashwagandha leaf extract and its component Withanone on scopolamine-induced changes in the brain and brain-derived cells

Understanding the binding affinity of noscapines with protease of SARS-CoV-2 for COVID-19 using MD simulations at different temperatures

Withanone and caffeic acid phenethyl ester are predicted to interact with main protease (Mpro) of SARS-CoV-2 and inhibit its activity

Antiviral effects and possible mechanisms of action of constituents from Brazilian propolis and related compounds

Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) and coronavirus disease-2019 (COVID-19): The epidemic and the challenges

Airway proteases: An emerging drug target for influenza and other respiratory virus infections. Current Opinion in Virology

Immunomodulatory and prophylactic efficacy of herbal extracts against experimentally induced chicken infectious anaemia in chicks: Assessing the viral load and cell mediated immunity

Evaluation of antiinfluenza effects of camostat in mice infected with non-adapted human influenza viruses

Structure, function, and evolution of coronavirus spike proteins

Transmission dynamics and evolutionary history of 2019-nCoV

Potential of large "first generation" human-to-human transmission of 2019-nCoV

TMPRSS2 is the major activating protease of influenza A virus in primary human airway cells and influenza B virus in human type II pneumocytes

Prostate-localized and androgen-regulated expression of the membrane-bound serine protease TMPRSS2

FDA-approved thiol-reacting drugs that potentially bind into the SARS-CoV-2 main protease, essential for viral replication

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

Withaferin A induces p53-dependent apoptosis by repression of HPV oncogenes and upregulation of tumor suppressor proteins in human cervical cancer cells

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

The expression level of angiotensin-converting enzyme 2 determine the severity of COVID-19: Lung and heart tissue as targets

Development of Remdesivir repositioning as a nucleotide analog against COVID-19 RNA dependent RNA polymerase

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

Cloning of the TMPRSS2 gene, which encodes a novel serine protease with transmembrane, LDLRA, and SRCR domains and maps to 21q22.3

Full-genome evolutionary analysis of the novel corona virus (2019-nCoV) rejects the hypothesis of emergence as a result of a recent recombination event

TMPRSS2 (transmembrane protease, serine 2)

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

Schr€ odinger Release 2020-1: Maestro 019-3 SR, Glide, LigPrep, Protein Preparation Wizard, Prime, Desmond Molecular Dynamics System

Anti-influenza virus effect of some propolis constituents and their analogues (esters of substituted cinnamic acids)

TMPRSS2: A potential target for treatment of influenza virus and coronavirus infections

Inhibition of Influenza A virus propagation by benzoselenoxanthenes stabilizing TMPRSS2 Gene G-quadruplex and hence downregulating TMPRSS2 expression

Clinical isolates of human coronavirus 229E bypass the endosome for cell entry

Middle East respiratory syndrome coronavirus infection mediated by the transmembrane serine protease TMPRSS2

A transmembrane serine protease is linked to the severe acute respiratory syndrome coronavirus receptor and activates virus entry

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

Investigation into SARS-CoV-2 Resistance of compounds in garlic essential oil

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

Nootropic potential of Ashwagandha leaves: Beyond traditional root extracts

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

Clinical characteristics and therapeutic procedure for four cases with 2019 novel coronavirus pneumonia receiving combined Chinese and Western medicine treatment

WHO Coronavirus Disease (COVID-19) Dashboard

Deceleration of senescence in normal human fibroblasts by Withanone extracted from ashwagandha leaves

Selective killing of cancer cells by leaf extract of Ashwagandha: Components, activity and pathway analyses

The membrane-anchored serine protease, TMPRSS2, activates PAR-2 in prostate cancer cells

The serine protease inhibitor camostat inhibits influenza virus replication and cytokine production in primary cultures of human tracheal epithelial cells

Protease inhibitors targeting coronavirus and filovirus entry

A novel coronavirus from patients with pneumonia in China

The computations were performed at the Bioinformatics Centre supported by the Department of Biotechnology (DBT) Govt. of India at IIT Delhi.

All data generated or analyzed during this study are included in this published article.

No potential conflict of interest was reported by the authors.

This study was supported by the funds granted by AIST (Japan) and DBT (Government of India).