key: cord-0747487-z2fdced4
authors: Kapoor, Neha; Ghorai, Soma Mondal; Kushwaha, Prem Kumar; Shukla, Richa; Aggarwal, Charu; Bandichhor, Rakeshwar
title: Plausible mechanisms explaining the role of cucurbitacins as potential therapeutic drugs against coronavirus 2019
date: 2020-11-25
journal: Inform Med Unlocked
DOI: 10.1016/j.imu.2020.100484
sha: 15f364f49ba6af3463da2f4584040a1c1bf944c0
doc_id: 747487
cord_uid: z2fdced4

In the year 2019, the potent zoonotic virus, severe acute respiratory syndrome coronavirus 2 (SARS-CoV2) began to rage globally, which resulted in the World Health Organization (WHO) declaring it as a pandemic on March 11th, 2020. Although extensive research is currently ongoing worldwide to understand the molecular mechanism and disease pathogenicity of SARS-CoV-2, there are still many nuances to elucidate. Therefore, developing an appropriate vaccine or therapeutic drug to combat coronavirus 2019 (COVID-19) is exceedingly challenging. Such scenarios require multifaceted approaches to identify suitable contenders for drugs against COVID-19. In this context, investigating natural compounds found in food, spices, and beverages can lead to the discovery of lead molecules that could be repurposed to treat COVID-19. Sixteen cucurbitacin analogues were investigated for activity against the SARS-CoV-2 main protease protein (Mpro), angiotensin-converting enzyme 2 (ACE2) binding receptor, nonstructural protein 12 (NSP12) RNA-dependent RNA polymerase (RdRp), NSP13 helicase, and Janus kinase 2 (JAK2)/signal transducer and activator of transcription 3 (STAT3) pathway using several relevant tools and simulated screening methods. All key proteins were found to bind efficiently only with cucurbitacin G 2-glucoside and cucurbitacin H with the lowest global energy. Further, the absorption, distribution, metabolism, and excretion (ADME) of all the cucurbitacins were analysed to explore their drug profiles. Cucurbitacin G 2-glucoside and H showed the best hits and all the analogues showed no adverse properties that would diminish their drug-likeness abilities. The encouraging results of the current study may lay the foundation for future research and development of effective measures and preventive medications against SARS-CoV-2.

Over the past two decades, humans have grappled with the spread of zoonotic pathogenic strains of many viruses such as the severe acute respiratory syndrome coronavirus (SARS-CoV) in 2002, followed by the Middle East respiratory syndrome-related (MERS)-CoV in 2012 with ~10% and ~36% mortality rates, respectively [1] . Although the fatality rate was lower, the advent of SARS-CoV-2 in 2020 heralded the spread of a highly infectious and pathogenic strain of betacoronavirus causing acute respiratory distress syndrome (ARDS), often leading to lung dysfunction, arrhythmia, and death. To date, there are no marketable antiviral agents or vaccines available against SARS-CoV-2. Developing therapies for the SARS-CoV-2-induced coronavirus disease 2019 (COVID-19) requires the complete elucidation of the pathology and molecular mechanisms underlying the viral mode of infection.

Therefore, it is crucial to understand and identify potential targets for drug design against both SARS-CoV-2 and human proteins, wherever they interact. Searching for biologically promising compounds from natural sources could be an alternate strategy for developing remedies against SARS-CoV-2. Although a few candidate drugs are currently under consideration as effective antivirals for SARS-CoV-2, we believe that the evaluation of small molecule drug candidates from natural compounds could potentially provide durable, broadspectrum host-directed therapies. Earlier studies have shown that cucurbitacins have remarkable inhibitory activity against the herpes simplex virus (HSV)-1 [2] and, thus, may have the potential to be developed as antiviral drugs against COVID- 19 .

Traditional belief as well as scientific evidences has proven that Cucurbita pepo, the oldest known cultivated species of the Curcurbiteae family, has immense medicinal and nutritional benefits [3] . Other species of the Curcurbiteae family, which are indigenous to the western hemisphere, include cucumbers, honeydew melons, cantaloupe, watermelons, and zucchini. To date, six cucurbitacins (B, D, E, E glucoside, Iso B, and Iso D) have been isolated from these species. Cucurbitacins are tetracyclic triterpenoids, and both the fruit and seeds of these plants are used as traditional medicines. Pulp with high carotenoid content and low fat content is known to have both beneficial physiological and immunomodulatory functions [4, 5] .

Pumpkin seeds have nutraceutical properties attributable to the rich content of various elements such as antioxidant vitamins, tocopherols, and carotenoids [6, 7] . They are also known to contain proteins, polyunsaturated fats, and phytosterols [8, 9] . Recently, cucurbitacins B, E, and D have been shown to exhibit potent antiviral activity against the bovine viral diarrhoea virus (BVDV), hepatitis C virus (HCV) [10] , and HSV [11] . In the current unprecedented times of the COVID-19 pandemic, there is an urgent need to explore the effects of natural compounds such as cucurbitacins, and carefully bridging the gap in our understanding of all the interacting proteins of both SARS-CoV-2 and humans would facilitate this process.

Many strategic attempts have been made to interfere with the different steps in the coronavirus replication cycle. Previous outbreaks such as those of SARS-CoV and MERS have trained us to target key enzymes involved in viral replication such as proteases, polymerases, and helicases, as attractive therapeutic approaches [12] . A few antivirals with protease-and polymerase-related mechanisms of action are under investigation to counteract SARS-CoV-2 infection. One example is remdesivir, which is a nucleoside analogue and inhibitor of the RNAdependent RNA polymerase (RdRp) of coronavirus [13] .

Enzyme kinetics have shown that remdesivir triphosphate (RDV-TP) competes with natural ATP with four times more selectivity for position (i), thus, inhibiting RNA synthesis at position (i+5). This process delays chain termination and is a plausible mechanism for drug targeting in the SARS-CoV RdRp complex [14] [15] [16] . Thus, in the present study, we compared the cucurbitacin molecule with the RDV to assess its suitability as a target against the SARS-CoV RdRp complex. Cucurbitacin was also examined as a potential lead molecule to target the Mpro of SARS-CoV-2, which is the main enzyme mediating both viral replication and transcription [17] . Using combinatorial techniques of both structure-based virtual and high-throughput screening, innumerable pharmacologically active compounds, natural compound drug candidates, and repurposed drugs have been shown to inhibit Mpro proteases [18] .

The various envelope proteins such as S, M, and E on coronaviruses exhibit extensive variability, but the NSPs are quite conserved, including the protein helicase (NSP13), which was first studied in SARS-CoV and is an NTP-dependent 5' to 3' direction enzyme that coverts duplex oligonucleotides (RNA or DNA) by catalysing the unwinding into single strands [19, 20] . Thus, by targeting ATP-binding or direct NTPase activity, nucleic acids binding to the helicase, and blocking activities such as helicase translocation, inhibitors of NSP13 activity offer a potential therapeutic option against coronaviruses including SAR-CoV-2. Some natural compounds such as benzotriazole, imidazole, imidazodiazepine, phenothiazine, quinoline, anthracycline, triphenylmethane, tropolone, pyrrole, acridone, small peptide, and bananin derivatives are under consideration as potent inhibitors for NSP13 activity. Cucurbitacins have also been used as inhibitors of helicases.

Studies have proven that angiotensin-converting enzyme 2 (ACE2), a metallopeptidase, competently binds the S1 domain of the SARS-CoV "S" protein and mediates entry of their target cells. It has also been shown that the related enzyme ACE1 does not bind to the S1 domain and that the ACE2 receptor is the only functional receptor for SARS-CoV [21, 22] . The entry of SARS-CoV-2 S protein is known to be blocked by the anti-ACE-2 antibody [23] . Recently, plasma therapy with infused convalescent antibodies has been used to treat severe COVID-19 patients [24] . We conducted in silico analysis to determine if cucurbitacin disrupts the interaction between the virus and ACE2 receptors and, thus, might be a potential effective therapy for COVID-19.

Finally, we also attempted to study the signalling mechanism mediating the release of proinflammatory cytokines in SARS-CoV-2 infection. The release of interleukin (IL)-6, IL-1β, and IL-12 is known to cause cytokine storm, inducing multiple organ failure in patients with acute conditions. These cytokines are released from various innate immune cells (monocytes, neutrophils, and NK cells), which in turn, activate T-lymphocytes via the JAK/STAT pathway J o u r n a l P r e -p r o o f [25] . Thus, antagonists of the JAK/STAT pathway may be correlated in reducing the cytokine storm and, thus, saving lives. In the present study, cucurbitacins were also explored as inhibitors of relevant signalling pathways.

Cucurbitacin was selected for screening for activity against SARS-CoV-2, and its threedimensional (3D) structure was retrieved from PubChem (https://pubchem.ncbi.nlm.nih.gov/) in the SDF format. The 3D structure of cucurbitacin was minimized with retained specified chirality using the default force field OPLS3 of ligprep/maestro and epik to generate the possible state at the default pH. The molecular enzymes of SARS-CoV-2 NSP12 (Protein Databank [PDB] Id 6NUR) [14] with bound cofactor NSP7 and NSP8, the main protease (PDB Id 6LU7), JAK2 (PDB Id 4GFM), ACE2 (PDB Id 6VW1), and NSP13 helicase (6ZSL) were targeted by the selected cucurbitacin to inhibit the viral infection of COVID-19.

These protein structures were retrieved from the PDB (http://www.rcsb.org/pdb) and prepared from protein preparation using the wizard function of Maestro 12.4 in Schrodinger 2020-2. The 3D structures of the proteins were pre-processed by choosing the default option and filling the missing side chains and loops with prime. Further, the structures were modified by removing hets/water within 5 Å, and finally, refined by assigning the H-bonds, removing water within 3 Å, and performing retrained minimization by choosing the OPLS3 force field.

The protein-binding site "sitemap" was constructed using the standard default parameter setting of the maestro suite. The sitemap also facilitated the characterization of hydrophobic, hydrophilic, hydrogen donor, and hydrogen acceptor residues in the binding site. The top-ranked potential binding sites were identified, and the best predicted binding site was chosen based on a Dscore value > 1.

The receptor grid was generated by using default parameter settings from maestro suite. Predicted sitemap binding sites were used for receptor grid generation, and further predicted receptor grids were used for ligand docking. Docking calculations were performed using the standard default parameter setting of the ligand-docking task of Maestro in which Cucurbitacin was docked into the predicted receptor grid with extra precision along with XP descriptor information. The ligand sampling was kept flexible while the proteins were considered as rigid structures and epik state penalties were applied. Finally, for the output file, a pose viewer file was chosen and post-docking minimization was performed. The best dock score was identified as a "hit" Pymol 2.4.0 was also used for visualisation and figure generation.

Drug performance and pharmacological efficacy are critically influence by four major parameters: absorption, distribution, metabolism, and excretion (ADME). Prior knowledge of the ADME and toxicological (Tox) parameters of drugs enables the control of their pharmacological activity and pharmacokinetics. Thus, pharmacology and effectiveness are mainly measured through the analysis of factors that influence the kinetics of drug doses and contact with the tissues in an organism. In this study, we used the qikprop function of maestro 12.4 to determine the ADME/Tox properties of all 16 of the cucurbitacin analogues. The server carefully predicts the toxicity endpoints by not only predicting the 2D similarity to compounds with known median lethal doses (LD 50 ), but also by drawing parallels on fragment and molecular similarity and fragment tendency.

Considering the already known therapeutic activities of cucurbitacins and using XP docking methods, we identified all known cucurbitacins with structures that could be used to inhibit multiple targets and steps in the life cycle of coronaviruses (Table 1 ). In the present study, we selected 16 of the best cucurbitacins with known structural configurations designated as A, B, C, D, E, F, G, H, I, J, K, L, O, P, Q, R, and S ( Fig. 1 ). Based on binding affinity and structurebased potential, we virtually screened using an effective procedure to identify ligand hits and conduct lead optimization for cucurbitacins as drug candidates. We docked 16 potential candidates and arranged a congeneric series for binding affinity prediction using methods such as macrocycle pose prediction by engaging polypeptides, predicting protein-ligand complex geometries, and free energy perturbation or MM-GBSA.

Consequently, we identified key molecules likely to bind well to the cucurbitacins. Thus, based on structure-based virtual screening and the induced fit docking protocols for the most effective ligand docking, we determined the most important functional cellular machinery for SARs-CoV2 and host interaction. These were the SARS-CoV-2 polymerase, NSP12 (known for direct supply of positive sense RNA, which increases the rate of protein synthesis); the active viral protease, Mpro; COVID-19 protein, NSP13 helicase (known for repairing cellular components) and the host's ACE2 receptor (on the epithelial and Vero E6 cells that open door for J o u r n a l P r e -p r o o f entry of coronavirus); and the host signalling machinery, JAK2 (activates T-lymphocytes via JAK/STAT pathway.

The sitemap showed that when all the selected cucurbitacins were bound to the kinase domain of JAK2, the binding mode and best pose (Emodel) of cucurbitacin G 2-glucoside were mostly similar to those of the other cucurbitacins ( Fig. 2a and b) . The only difference was observed in the H-bond energy of G, which was slightly higher than that of the remaining cucurbitacins. Further, G formed hydrogen bonds with the Asp 976, Glu 1015, ARG 980, Asp 939, Ser 936, and Leu 855 residues with an excellent binding affinity (-8.06 kcal/mol). These results predicted that this molecule could be used as an inhibitor of SARS-CoV2 infection. Similarly, among all the analogues, cucurbitacin G 2-glucoside with a Dscore of -9.43 kcal/mol ( Fig. 2 c and d) showed the best hit and appeared to be the most potent inhibitor of COV 2 Mpro protein. Thus, cucurbitacins are also active against the target Mpro of SARS-CoV-2.

The docking of all cucurbitacins showed very good interaction with the RBD of ACE2 and spike (S1) glycoproteins of CoV2 with a similar interaction profile. This was further corroborated by the sitemap, which determined that the RBD of ACE2 and spike glycoproteins had Dscore > 1. In particular, cucurbitacin G 2-glucoside showed the best binding affinity for ACE2 and formed hydrogen bonds with Tyr510, Tyr516, Glu402, Asn51, and His378 residues ( Fig. 2e and f) with hydrogen bond energy of -9.92 kcal/mol, which was much better than that the interaction of other reference molecules to date, including hydroxychloroquine.

The binding mode and best pose (Emodel) of cucurbitacin G 2-glucoside was slightly higher than that of the remaining cucurbitacins. However, for RdRp, cucurbitacin-K formed hydrogen bonds with Phe396, Thr394, Hip256, and Ala252 residues with binding affinity -9.57 kcal/mol ( Fig. 2g and h) , whereas remdesivir was reported to form a hydrogen bond with binding affinity of -5.91 kcal/mol (Table 1 ). Cucurbitacins are triterpenes that may potentially bind to the allosteric site of RdRp and probably provide key acuities to NSP12 polymerase catalysis and reliability that may act as a template for the design of novel antiviral therapeutics.

The sitemap was used to determine the possible nucleic acid binding site (Dscore = 1.07), which is sandwiched between RecA1 (241-443 amino acids [aa]) and RecA2 (444-596 aa) domain of the NSP13-helicase of COVID-19. In the nucleic acid binding site, Arg560 from motif IV was the only residue that might interact with the base and sugar moiety of the single strand DNA (ssDNA). Docking of the cucurbitacins led to the identification of the best possible hit for J o u r n a l P r e -p r o o f the NSP13 helicase of COVID-19 (Table 2 ). Our study clearly showed that one of the hydroxy groups of cucurbitacin H and cucurbitacin J interacted very strongly with Arg460 and revealed the best hit with a binding affinity of -10.069 kcal/mol (Fig. 3) .

The drug likeness of all 16 cucurbitacins was estimated using different ADME-related properties, namely physicochemical parameters, pharmacokinetics, lipophilicity, water solubility, and their medicinal chemistry. Analysis of the inhibitory effects of different cucurbitacins predicted that cucurbitacin G 2-glucoside and H are good candidates for antiviral drugs. The fast infection rate of pathogenic diseases has motivated drug discovery for these conditions, and recent technological advancements have increased the rate of drug discovery with prior knowledge of the drug candidate. We determined the anti-COVID-19 activity and drug-likeness of selected cucurbitacins by performing in-silico ADME/Tox calculations (Table 3) .

For cucurbitacin H, the solvent accessible surface area (SASA), predicted water-gas partition coefficient (QPlogPw), octanol/water partition coefficient (QPlogPo/w), aqueous solubility (QPlogS), half-maximal inhibitory concentration (IC 50 ) value for blockade of the human ether-a-go-go-related gene (HERG) K + channels (QPlogHERG), apparent Caco-2 cell permeability (QPPCaco), brain/blood partition coefficient (QPlogBB), skin permeability 

Reports by experts had already predicted in 2018 that a flu pandemic was only a matter of time and there could be numerous undiscovered viruses that could kill millions worldwide. In the past, we witnessed MERS and SARS and there is no shortage of articles regarding their mechanism of action. Extensive research are now been followed to understand the molecular mechanism and disease pathogenicity of SARS-CoV-2, but still developing an appropriate vaccine or therapeutic drug to combat coronavirus 2019 (COVID-19) is a challenging task. A multifaceted approach is required to rein this virus and we believed that the natural compounds found in food, spices, and beverages can lead to the discovery of lead molecules that could be repurposed to treat COVID-19. In our present study, we investigated sixteen cucurbitacin analogues for activity against the SARS-CoV-2 main protease protein (Mpro), angiotensin-converting enzyme 2 (ACE2) binding receptor, nonstructural protein 12 (NSP12) RNAdependent RNA polymerase (RdRp), NSP13 helicase, and Janus kinase 2 (JAK2)/signal transducer and activator of transcription 3 (STAT3) pathway using several relevant tools and simulated screening methods. We thus, could identify the host ACE2 receptor, CoV-2 RdRp polymerase, CoV-2 Mpro protease, CoV-2-helicase protein, and JAK2 signalling machinery as the best hits for specific cucurbitacins. All key proteins were found to bind efficiently only with cucurbitacin G 2-glucoside and cucurbitacin H with the lowest global energy. Also, the absorption, distribution, metabolism, and excretion (ADME) of all the cucurbitacins showed that Cucurbitacin G 2-glucoside and H showed the best hits and all the analogues showed no adverse properties that would diminish their drug-likeness abilities.

Plant-derived natural compounds have played an essential role in disease therapy. Cucurbitacins are triterpenes with various chemical structures found in the members of the Cucurbitaceae and numerous other plant families, such as Trichosanthes, Cucurbita, Cucumis, and Citrullus, which possess numerous medicinal properties [26] . Cucurbitacin B has been shown to inhibit HSV-1 and, thus, holds promise in the development of antiviral drugs [2] . The current COVID-19 pandemic has precipitated an urgent need for scientific strategies to repurpose already known lead compounds to develop possible therapies with activity against the SARS-CoV-2 virus.

Although cucurbitacins are known toxic compounds, they have been shown to possess various bioactivities and pharmacological efficiency against cancer, inflammation, diabetes, and atherosclerosis. There have been instances of severe poisoning in cattle that consumed bitter fruits of the Cucumis and Cucurbita species. The dose range of the toxicity of cucurbitacins based on the few in vivo toxicity reports has been reported to be between 2 and 12.5 mg/kg. Double bonds at C-23 and the acetyl group at C-25 have been found to modulate the toxicity of cucurbitacins, and their bioactive and toxic dose have been found to be close [3] . Among the various bioactivities shown by cucurbitacins, plant species such as Momordica are held in high esteem in many countries as traditional medicine for certain metabolic diseases such as diabetes [27] .

The antidiabetic property of cucurbitacins is mediated by activation of the AMP-activated protein kinase (AMPK) pathway (a major regulatory pathway for GLUT4 translocation) [28, 29] . These compounds also show anti-atherosclerosis by inhibiting the generation of lipid oxidation products such as malonaldehyde (MAD) and 4 hydoxynonenal (4 HNE) [3] . Cucurbitacin D improves blood circulation through inhibition of Na+/K+ ATPase [30] . The concentration of cucurbitacin C in the leaves of plants has also been reported to be an important parameter in the spider mite resistance of C. sativus, perhaps by acting as an antagonist of the spider mite ecdysteroid receptor [31] . However, the most promising therapeutic properties of cucurbitacins are their anticancer [32, 33] and anti-inflammatory [34, 35] activities.

The diverse pharmacological activity of cucurbitacins is due to the variety of their side chain derivatives [36, 37] , which can be used to group these compounds and their derivatives into 12 main categories [33] . There are 16 main molecules designated cucurbitacin A to T, and hundreds of derivatives of the parent molecule. Chemically, cucurbitacins are derived from the basic cucurbitane ring skeleton, a triterpene hydrocarbon (International Union of Pure and Applied Chemistry [IUPAC] name: 19 [10-9β]-abeo-5α-lanostane), which is modified by groups containing oxygen and double bonds to produce numerous cucurbitacins with distinctive characteristics [38] . Pure cucurbitacins and their glycosidic forms have been isolated from plant tissues using non-polar solvents to remove undesired components [39] . Most cucurbitacins are soluble in non-polar solvents, but are insoluble in ether, whereas they are only slightly soluble in water. The majority of cucurbitacins are usually crystalline or exist as needle-shaped particles at 68 to 77°F (20 to 25°C) except for cucurbitacin H, which is an amorphous solid. Most cucurbitacins are derived as 2-O-glucoside, except for cucurbitacin D, which lacks the acetyl group at 25-OH, but is the most ubiquitous cucurbitacin isolated. Cucurbitacin G 2-glucoside and H have the same structures but differ in the configuration of their hydroxyl groups at position C-24, which has not been established [33] . Cucurbitacin R was identified to be 23,24dihydrocucurbitacin D (DHCD), hence, has been moved to the cucurbitacin D group based on this description [40] . Similarly, Cucurbitacin J and K differ only in the configuration of their hydroxyl groups at position 24 and their chemical validation is still pending [41] .

The earliest studies during the onset of the pandemic were conducted to understand whether the new SARS-CoV2 binds to the host cells in a similar manner to that of the known SARS-CoV. Stefan Pohlmann and his team were among the first to discover that both viruses used ACE2 surface receptors to penetrate cells and the same protease to become infectious [42] . Binding of the CoV-2 homotrimeric spike glycoprotein (comprising a S1 subunit and S2 subunit in each spike monomer) to the host ACE2 receptor is the first critical step in the entry of SARS-CoV into target cells [43] [44] [45] . In addition to the lungs, ACE2 receptors are expressed in many extrapulmonary tissues, including the heart, kidney, endothelium, and intestine [46] [47] [48] [49] [50] . The receptor-binding domain (RBD) of the S1 subunit is responsible for the interaction between the CoV spike and ACE2 receptor [49] . Structural analysis of various spike residues has identified highly conserved or identical side chain residues in the SARS-CoV-2 RBD to those of the SARS-CoV RBD, suggesting convergent evolution and improved binding to the ACE2 receptor [51] . Notably, ACE2 is widely expressed on the luminal surface of intestinal epithelial cells, functioning as a co-receptor for nutrient uptake, particularly for amino acid resorption from food [52] . Therefore, the intestine could be a potent entry point for CoV-2, which gave us the idea of using oral formulations of natural products such as cucurbitacins as antagonist of the ACE2 receptor. Docking results showed that most of the selected cucurbitacins interacted well with the RBD of ACE2 and S1 glycoprotein of CoV2, thus, providing evidence supporting its promising role as a target site.

In the last decade, the world witnessed the emergence of numerous coronavirus infections, such as SARS and MERS, which showed that SARS-CoVs possess the largest RNA genomes and require fidelity for the RNA synthesis complex to replicate. The new emerging viruses of diverse CoV families require innovative research approaches to generate antiviral strategies that target conserved components of the viral life cycle. This includes strategies mainly J o u r n a l P r e -p r o o f focusing on the viral machinery responsible for the replication and transcription of the positivestrand viral RNA genome. NSP12 is a complex RNA synthesis machinery of CoV formed as a product of the ORF1a and ORF1ab viral polyproteins [53] . Sequence conservation and homology modelling studies have revealed mostly highly conserved surfaces on the NSP12 CoV family of viruses at the template entry, template-primer exit, NTP tunnels and the polymeraseactive site [14] . Interestingly, SARS-CoV NSP12 contains two metal-binding sites with two zinc atoms, which are distal to known active sites as well as protein-protein and protein-RNA interactions. Coronavirus NSPs usually have bound zinc atoms in NSP3, NSP10, NSP13, and NSP14 that point to the extensive utilisation of zinc ions for folding proteins of the viral replication complex [54] [55] [56] . To date, the nucleotide-based prodrug "remdesivir" is the only known inhibitor of COVID-19, which inhibits the NSP12 CoV-RdRp by creating a covalent bond between the 5′hydroxy of remdesivir and 3′hydroxy of the RNA ribose in the RNA-binding site [57] . Cucurbitacin K forms a high-affinity hydrogen bond with NSP12 polymerase, thus, targeting these NSPs may hold the key to blocking viral replication.

In continuation of experiments targeting viral NSPs, subsequent studies in SARS/MERS did not clearly specify the function of NSPs such as NSP13 or helicases. Based on in silico modelling and virtual screening, Mirza and Froeyen [58] predicted a similarity in the amino acid/protein profile of SARS-CoV-2 NSP13 to those of SARS-and MERS-NSP13. NSP13 was thought to probably interact with other NSPs such as NSP7, NSP8, and NSP12 to contribute to viral replication [59, 60] . Adedeji et al. [61] demonstrated that the helicase activity of NSP13 is enriched in the presence of RdRp NSP12, suggesting that these proteins interact in the functional replication complex and that the protein-protein interaction contributes to the efficiency of viral replication. SARS-CoV NSP13 catalyses the NTP-dependent 5′ to 3′ unwinding reaction and converts double-stranded oligonucleotides into single strands [19, 20] . It has been proposed that molecules with the potency to inhibit putative ATP binding sites associated with NTPase activity could be used for future therapeutics. NTPase active site residues of SARS-CoV-2, including Lys288, Ser289, Asp374, Glu375, Gln404, and Arg567 are conserved throughout the coronavirus family. Thus, many natural compounds that inhibit ATP hydrolysis could be considered potential inhibitors of SARS-CoV2 replication. The present study also showed that one of the hydroxy groups of cucurbitacins H and J interacted very strongly with Arg460 and, thus, these compounds have the potential to inhibit SARS-CoV2.

SARS-CoV-2 replication and transcription requires encoding of two overlapping polyproteins (pp1a and pp1ab), which undergo extensive autolytic cleavage and proteolytic processing at 11 conserved sites through the action of 33.8 kDa Mpro (also known as 3C-like protease) to release from the polyproteins. The 3CL protease of coronaviruses facilitates viral assembly by cleaving polyproteins and this functional importance of Mpro in the viral life cycle makes it an attractive target for antiviral drug design [18, 62] . The crystal structure of the Mpro of SARS-CoV-2 is highly similar to that of SARS-CoV Mpro, with only a 0.53 Å r.m.s difference between the two free enzymes [63] . In our study, cucurbitacin G showed potent inhibitory activity against the CoV2 Mpro protein.

Viral entry into the host normally activated inflammatory pathways, creating an anti-viral environment to stop the virus infection pathogenesis; however, an aberrant or uncontrolled host immune response may cause host cell and tissue damage [64] . The severity of COVID-19 infection is due to an exaggerated inflammatory response with the release of a large amount of pro-inflammatory cytokines. Most studies of SARS-CoV-2 have related to this cytokine storm to lung injury and multiple organ failure, which is activated via the JAK/STAT pathway [65] [66] [67] [68] [69] . Phosphorylation of JAK signalling molecules activates STAT proteins, which induces a cascade of cytokines [70, 71] . Cucurbitacin B was reported to be an anti-proliferative agent in pancreatic cancer by inhibiting the activation of JAK2, STAT3, and STAT5 components [72] . Suppression of the JAK/STAT pathway increased the expression of p21WAF1 independently of p53 activity, thus, inducing cell cycle arrest. This mechanism could be exploited to stop the proliferation and recruitment of T-lymphocytes, thereby inhibiting cytokine production. The mechanisms of the anticancer activity of cucurbitacins are possibly mediated by inhibition of the JAK/STAT3 signalling pathway [73, 74] . Interestingly, the same pathway is used by coronaviruses to activate T-lymphocytes at the later stage of infection, which leads to cytokine storm and subsequent mortality of the infected person. Moreover, cucurbitacins are also used as immunosuppressants, and they have been shown to inhibit the expression of the surface markers CD69 and CD25 required for the activation of T lymphocytes [75] . In the present study, we attempted to show that cucurbitacins can inhibit the JAK/STAT3 pathway and, thus, could inhibit SARS-CoV2 infection. Intriguingly, cucurbitacins have been recognised to possess anti-inflammatory activity and inhibit the expression of TNF and proinflammatory mediators such as nitric oxide (NO) synthase 2 and cyclooxygenase 2-mediated inhibition of NO generation [34, 35] . A similar strategy could be employed to curb the surging production of cytokines in SARS-CoV2 infection, where cucurbitacins could be a critical source of future anti-inflammatory drugs for severe viral infections.

The past two decades have seen a surge in the emergence of novel viral diseases such as the global incidence of viral infections caused by SARS-CoV and MERS, which belong to the coronavirus family. In addition, repeated outbreaks of other coronaviruses have affected the human population such as the 229E, HKU1, NL63, and OC43. Incidentally, the current SARS-CoV2 has an 80% nucleotide identity with the previous infective viruses, and there is a fairly good understanding of effective therapies against this family of viruses. The information acquired from investigations of these coronaviruses suggest that a promising therapeutic strategy would be to target the fairly conserved key enzymes involved in the replication cycle of the virus. This suggests that known therapies that target these enzymes, even in closely related coronavirus species, could be used in SARS-CoV-2 infection. Proteases and polymerases of SARS-CoV2 as targets for antiviral drugs, have gained attention in recent weeks in clinical trial of drugs such as indinavir, saquinavir, lopinavir/ritonavir, and remdesivir for COVID-19 with. In contrast, helicases have not been explored as much as the proteases and polymerases for this purpose. Moreover, the structural proteins of coronaviruses such as spike (S) proteins are known to show considerable variability among coronavirus species, and targeting the spike proteins has been an extremely laborious task. Key NSPs such as the helicase NSP13 have shown similar homology to helicases of other coronavirus species. Our present perspective is to present the potential application of helicase as one of the target enzymes that could be inhibited using natural compounds. Moreover, currently, there are very limited known studies of the use of natural compounds as repurposed drugs with antiviral activity administered as a single native molecule or in combination therapy. There are multiple steps in the life cycle of a virus and in the present study, cucurbitacin analogues were targeted against the virus at every possible process, such as attachment, entry, replication, morphogenesis, and exit from the host cell. The results suggest that the natural molecule cucurbitacin has potential as a therapeutic drug against COVID-19, especially considering the very encouraging outcome observed with cucurbitacin G 2-glucoside, cucurbitacin-H and cucurbitacin-J. However, all investigational agents for SARS-CoV-2 are under strict regulation by the World Health Organization (WHO) and this study urges the scientific community to continue conducting further in vivo trials for the experimental validation of our findings.

The authors declare that they have no conflict of interest 

Middle East respiratory syndrome and severe acute respiratory syndrome: current therapeutic options and potential targets for novel therapies

Cucurbitacin B interacts synergistically with antibiotics against Staphylococcus aureus clinical isolates and exhibits antiviral activity against HSV-1

Cucurbitacins-an insight into medicinal leads from nature. Pharmacognosy reviews

Antidiabetic potentials of Momordica charantia: multiple mechanisms behind the effects

A review on the medicinally important plants of the family Cucurbitaceae

Oil and tocopherol content and composition of pumpkin seed oil in 12 cultivars

Amino acid, mineral and fatty acid content of pumpkin seeds (Cucurbita spp) and Cyperus esculentus nuts in the Republic of Niger. Plant foods for human nutrition

Phytosterol composition of nuts and seeds commonly consumed in the United States

Fatty acid composition of seed and leaf oils of pumpkin, walnut, almond, maize, sunflower and melon. Chemistry of Natural Compounds

Anticancer and Antiviral Activities of Cucurbitacins Isolated From Cucumis Prophetarum var. Prophetarum Growing in the Southwestern Region of Saudi Arabia

Effects of Z-venusol, isolated from Gunnera perpensa L., on IL-6 and cAMP activity in human breast cancer cells in vitro

Lipid-lowering therapy and reninangiotensin-aldosterone system inhibitors in the era of the COVID-19 pandemic

The antiviral compound remdesivir potently inhibits RNA-dependent RNA polymerase from Middle East respiratory syndrome coronavirus

Structure of the SARS-CoV nsp12 polymerase bound to nsp7 and nsp8 co-factors

Structural and molecular basis of mismatch correction and ribavirin excision from coronavirus RNA

One severe acute respiratory syndrome coronavirus protein complex integrates processive RNA polymerase and exonuclease activities

Structures of two coronavirus main proteases: implications for substrate binding and antiviral drug design

Structure of M pro from SARS-CoV-2 and discovery of its inhibitors

Structure and mechanism of helicases and nucleic acid translocases

Delicate structural coordination of the Severe Acute Respiratory Syndrome coronavirus Nsp13 upon ATP hydrolysis

Angiotensin-I-converting enzyme and its relatives

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

Blocking of the high-affinity interactionsynapse between SARS-CoV-2 Spike and human ACE2 proteins likely requires multiple high-affinity antibodies: An Immune perspective

Treatment of COVID-19 with convalescent plasma: lessons from past coronavirus outbreaks. Clinical Microbiology and Infection

The COVID-19 cytokine storm; what we know so far

New cucurbitacins from Phormium tenax and Marah oreganus

Anti-diabetic properties and phytochemistry of Momordica charantia L.(Cucurbitaceae)

AMP-activated protein kinase and malonyl-CoA: targets for treating insulin resistance? Drug Discovery Today: Therapeutic Strategies

AICAR administration causes an apparent enhancement of muscle and liver insulin action in insulin-resistant high-fat-fed rats. diabetes

Active ingredients in Chinese medicines promoting blood circulation as Na+/K+-ATPase inhibitors

Role of cucurbitacin C in resistance to spider mite (Tetranychus urticae) in cucumber (Cucumis sativus L.)

Lead compounds from Cucurbitaceae for the treatment of cáncer. Phytochemicals-Isolation, Characterization and Role in Human Health

Cucurbitacins and cucurbitane glycosides: structures and biological activities

Search for a novel antioxidant, antiinflammatory/analgesic or anti-proliferative drug: Cucurbitacins hold the ace

Cucurbitacin E exhibits anti-inflammatory effect in RAW 264.7 cells via suppression of NF-κB nuclear translocation

Cucurbitacins with unusual side chains from Picrorhiza kurroa

Cucurbitacins are insect steroid hormone antagonists acting at the ecdysteroid receptor

Biogene Gifte-Biologie, Chemie

Method of isolating cucurbitacin

Constituents of Ecballium elaterium XXIII: Cucurbitacins and hexanorcucurbitacin

The Merc index-An encyclopedia of Chemicals and drugs

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

Structure, function, and antigenicity of the SARS-CoV-2 spike glycoprotein

Functional assessment of cell entry and receptor usage for SARS-CoV-2 and other lineage B betacoronaviruses

Back to the spring of Wuhan: facts and hope of COVID-19 outbreak

Tissue distribution of ACE2 protein, the functional receptor for SARS coronavirus. A first step in understanding SARS pathogenesis

Organ distribution of severe acute respiratory syndrome (SARS) associated coronavirus (SARS-CoV) in SARS patients: implications for pathogenesis and virus transmission pathways

Multiple organ infection and the pathogenesis of SARS

Essential role for collectrin in renal amino acid transport

Angiotensin-converting enzyme 2 is an essential regulator of heart function

Receptor recognition by the novel coronavirus from Wuhan: an analysis based on decade-long structural studies of SARS coronavirus

ACE2 links amino acid malnutrition to microbial ecology and intestinal inflammation

The coronavirus replicase

Crystal structure of Middle East respiratory syndrome coronavirus helicase. PLoS pathogens

Structural basis and functional analysis of the SARS coronavirus nsp14-nsp10 complex

Crystal structure of nonstructural protein 10 from the severe acute respiratory syndrome coronavirus reveals a novel fold with two zinc-binding motifs

Remdesivir is a direct-acting antiviral that inhibits RNA-dependent RNA polymerase from severe acute respiratory syndrome coronavirus 2 with high potency

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

Genome-wide analysis of protein-protein interactions and involvement of viral proteins in SARS-CoV replication

Analysis of intraviral protein-protein interactions of the SARS coronavirus ORFeome. PloS one

Mechanism of nucleic acid unwinding by SARS-CoV helicase

Designing of improved drugs for COVID-19: Crystal structure of SARS-CoV-2 main protease M pro

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

Immunity to viruses

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

Clinical predictors of mortality due to COVID-19 based on an analysis of data of 150 patients from Wuhan, China. Intensive care medicine

Clinical and immunologic features in severe and moderate forms of Coronavirus Disease

Diagnostic utility of clinical laboratory data determinations for patients with the severe COVID-19

Analysis of clinical features of 29 patients with 2019 novel coronavirus pneumonia. Zhonghua jie he he hu xi za zhi= Zhonghua jiehe he huxi zazhi= Chinese journal of tuberculosis and respiratory diseases

Cucurbitacin Q: a selective STAT3 activation inhibitor with potent antitumor activity

Blockade of constitutively activated Janus kinase/signal transducer and activator of transcription-3 pathway inhibits growth of human pancreatic cancer. Cancer letters

Cucurbitacin B induces apoptosis by inhibition of the JAK/STAT pathway and potentiates antiproliferative effects of gemcitabine on pancreatic cancer cells

Cucurbitacin B inhibits STAT3 and the Raf/MEK/ERK pathway in leukemia cell line K562. Cancer letters

Cucurbitacin E inhibits Huh7 hepatoma carcinoma cell proliferation and metastasis via suppressing MAPKs and JAK/STAT3 pathways

Cucurbitacin IIb exhibits anti-inflammatory activity through modulating multiple cellular behaviors of mouse lymphocytes

The authors would like to acknowledge the Principal, Hindu College, University of Delhi, and the Director of Birla Institute of Technology, Mesra, Ranchi for the technical support provided for this research work.

The authors would like to acknowledge the Schrödinger for the technical support provided for this research work.J o u r n a l P r e -p r o o f