International Journal of Clinical Biochemistry and Research

Print ISSN: 2394-6369

Online ISSN: 2394-6377

CODEN : IJCBK6

International Journal of Clinical Biochemistry and Research (IJCBR) open access, peer-reviewed quarterly journal publishing since 2014 and is published under auspices of the Innovative Education and Scientific Research Foundation (IESRF), aim to uplift researchers, scholars, academicians, and professionals in all academic and scientific disciplines. IESRF is dedicated to the transfer of technology and research by publishing scientific journals, research content, providing professional’s membership, and conducting conferences, seminars, and award more...

  • Article highlights
  • Article tables
  • Article images

Article statistics

Viewed: 154

PDF Downloaded: 41


Get Permission Velan: In silico investigation of therapeutic potential of phytosteroids against SARS-CoV2 infection via inhibition of viral multiplication and attenuation of host cytokine storm


Introduction

In December 2019, in the Chinese city of Wuhan, there have been reports of an outburst of clinical cases of patients with pneumonia-like symptoms as a result of an infection caused by a novel virus belonging to the coronavirus family.1 Molecular characterization as per World Health Organization (WHO) standards, categorized it as Severe Acute Respiratory Syndrome Coronavirus 2 (SARS-CoV2), based on which it was further named as novel coronavirus-2019 (2019-nCoV), and hence the disease was named as coronavirus disease-2019 (COVID-19).2 It has been recognized as a highly contagious and pathogenic virus that may cause a host of symptoms ranging from mild cold to fatal illnesses triggered by acute respiratory distress syndrome (ARDS). Cytokine storm, an important causal factor of ARDS, results in excessive and uncontrolled release of pro-inflammatory cytokines. Recent reports suggest that elevated cytokine profile that includes interleukin-1 (IL-1), interleukin-2 (IL-2), interleukin-10 (IL-10) and tumour necrosis factor-alpha (TNFα) are abnormally elevated in a COVID-19 patient’s blood serum.3 Cytokine storm results in abnormally high cell inflammation that leads to acute lung injury, thus resulting in ARDS, in addition to extrapulmonary systemic hyperinflammation syndrome. Hence, SARS-CoV2 infection causes cytokine storm that in turn triggers ARDS and multiple organ failure, resulting in mortality.4 In order to minimize the severity of SARS-CoV2 infection, attenuation of the cytokine storm may be an essential step. Among the recently reported methods, low-dose corticosteroids treatment inhibits secretion of pro-inflammatory markers including IL-1, IL-6, granulocyte colony stimulating factor (GCSF), TNF-α and cyclooxygenase (COX2) synthesis, thus minimizing severe lungs damage associated with COVID-19.5 However, inflammation inhibition alone cannot overcome COVID-19. Identification of therapeutic molecule that has dual properties of anti-inflammation and anti-viral multiplication may lead to an ideal drug treatment against COVID-19.

Figure 1

Chemical structure of the following compounds (A) Brassinolide; (B) 28-Homocastasterone; (C) 24-epibrassinolide (24-EB); and (D) Hydroxychloroquine (HCQ)

https://s3-us-west-2.amazonaws.com/typeset-prod-media-server/53370576-dd65-4f6b-b36c-201fe8af681dimage1.png

Brassinolide (BS), 28-Homocastasterone (28-HC) and 24-epibrassinolide (24-EB) are plant hormones belonging to brassinosteroid family (Figure 1). They are biosynthesized by CYP72B1 enzyme and regulates a wide range of physiological and metabolic processes including seed germination, senescence and stress response.6 Previous studies suggest that phytosteroids have antiviral potential against herpes simplex virus, Junin virus, arena viruses, and vesicular stomatitis viruses.7, 8 Mechanism of action of these phytosteroids mainly focuses on inhibition of viral multiplication.9 Conventional pharmacological methods of drug discovery are time-consuming, labour-intensive and expensive. An alternative concept of reverse pharmacology could be a major breakthrough in the field of drug discovery. The basis of reverse pharmacology includes in-silico analogue designing and ligand-receptor interaction.10

Figure 2

PDB structures of the following SARS-CoV proteins: (A) Spike protein; (B) Nucleocapsid protein; (C) Main protease

https://s3-us-west-2.amazonaws.com/typeset-prod-media-server/53370576-dd65-4f6b-b36c-201fe8af681dimage2.png

In the present in-silico study, we intend to investigate the anti-viral potential of specific phytosteroids, namely, BS, 28-HC and 24-EB, against SARS-CoV2 protein and human pro-inflammatory cytokines as an ideal drug candidate. These potential phytosteroids are hypothesized to have the ability to ct as inhibitor of viral proteins including Spike-protein (S-protein) subunit 1 (S1), nucleocapsid protein (N-protein), and main protease (Mpro) (Figure 2), in addition to the host pro-inflammatory proteins COX-2, prostaglandin synthase (PS), IL-1 and TNF-α. The docking studies was performed by using hydroxychloroquine (HCQ) as standard drug against the phytosteroids.

Materials and Methods

Compound preparation

In-silico structures of four molecules, namely, BS (CID:115196), 28-HC (CID:11038340), 24-EB (CID:443055) and HCQ (CID:3652) were downloaded from PubChem database (http://pubchem.ncbi.nlm.nih.gov/) as .SDF file format. The files were converted to .mol2 format using Open Babel software version 2.4.1. and optimized by means of ligand preparation script in AutoDock ver. 4.0. program. The inhibitors were prepared for docking as detecting root, torsion tree were set and saved in .pdbqt file format.

Protein structure and preparation

3D structures of S-protein S1 subunit (PDB ID: 6acc), N-protein (PDB ID: 4kxj) and Mpro (PDB ID: 6lu7), COX-2 (PDB ID: 4m1), PS (PDB ID: 1z9h), IL-1(PDB ID: 2nvh) and TNF-α (PDB ID: 2az5) were retrieved from RCSB protein databank (http://www.rcsb.org). AutoDock Tools (http://autodock.scripps.edu/resources/adt) facility supported the protein set up for docking. The script procedure employs the following ADT routines including removal of water molecules, addition of polar hydrogen atoms, assignment of Kollman charges and conversion of the protein files in. pdb format for molecular docking.

Grid box generation

3D structures of selected proteins and chemical structures were complexed together to forrm a grid. Therefore, the centroid of the small molecule in complex structures was chosen to generate grid points X = 60, Y = 60 and Z = 60 axis set for docking analysis. The grid file was generated by means of “grid generation panel” in AutoDock software ver. 4.0.11

Molecular docking simulation

Protein-ligand docking simulations were performed using AutoDock4 tool. For each ligand (chemical structure), 100 docking runs with default parameters were performed by treating protein as rigid and the ligand as flexible. The results were visualized using PyMol (The PyMOL Molecular Graphics System, Version 1.5.0.4 Schrodinger, LLC.), wherein all the conformations for each of the ligand was found to be within the cavity of protein indicating that the docking run was free from errors. The conformational clusters with lowest binding energy were considered for further analysis.10

Evaluation of the total binding energy

The AutoDock ver. 4.0 algorithm was applied to evaluate the total binding energy of phytosteroids against various proteins. Various docked conformations were obtained and the ones with lowest binding energy towards ligand binding cavity of protein were selected as the possible binding conformation, and considered for further interaction analysis. The final evaluations of the interactions between the target compounds and amino acid residues of the ligand binding cavity of protein were analyzed using online Protein Ligand Interaction Profiler tool (https://plip.biotec.tu-dresden.de/plip-web/plip/index).

Results

The SARS-CoV2 S1 subunit (S1) binding interaction with dietary phytosteroids were studied by in-silico docking simulations. Based on docking energies obtained, S1 exhibited highest binding affinity towards BS, followed by 28-HC, and 24-EB, thereby exhibiting -7.24, -8.27 and -8.2 Kcal/mole respectively, as compared to HCQ that exhibited binding affinity of -6.36 Kcal/mole. In S1-BS interaction, BS interacted with LYS929 residues via hydrogen bond and PRO710, ASP932, VAL 933, GLU999, ILE1000, LEU1006 residues via hydrophobic interaction (Table 1 and Figure 3). Similarly, in S1-28-HC interaction, GLY981 and GLN984 amino acid residues form hydrogen bonds with 28-HC, while THR943, LEU944, GLN947, GLN984, THR988 and LEU1006 residues form hydrophobic interactions (Table 1). Likewise, the binding studies performed between 24-EB and S1, indicates 24-EB interacting with amino acid residues GLY981, GLN984, SER985 via hydrogen bond, and with amino acid residues THR943, LEU944, GLN947, GLN984, THR988, TYR989 via hydrophobic interaction (Table 1 and Figure 3). HCQ, on the other hand, interacts with S1 amino acid residues GLN992 and ARG996 via hydrogen bond, while hydrophobic interactions are observed with ALA940, THR943, THR988, TYR989 amino acid residues with lowest binding affinity (-5.22 Kcal/mole) as compared to the three dietary phytosteroids (Table 1 and Figure 3).

Table 1

Molecular docking analysis ofSpike (S) protein S1 subunit interactions with various ligands

Spike (S) protein S1 subunit (PDB:6acc)

G-Score (Kcal/mole)

Hydrogen bounds

Hydrophobic interaction

BS (CID: 115196)

-7.24

LYS929.

PRO710, ASP932, VAL 933, GLU999, ILE1000, LEU1006.

28-HB (CID: 11038340)

-8.27

GLY981, GLN984.

THR943, LEU944, GLN947, GLN984, THR988.

24-EBS (CID: 443055)

-8.26

GLY981, GLN984, SER985.

THR943, LEU944, GLN947, GLN984, THR988, TYR989.

Hydroxychloroquine (CID: 3652)

-6.36

GLN992, ARG996.

ALA940, THR943, THR988, TYR989.

Table 2

Molecular docking analysis ofnucleocapsid protein interactions with various ligands

Nnucleocapsid protein (PDB: 4kxj)

G-Score (Kcal/mole)

Hydrogen bounds

Hydrophobic interaction

π-Stacking

Salt Bridges

BS (CID: 115196)

-11.24

VAL59, SER64, ARG164.

PHE57, PRO61, HIS104, TYR124.

TYR63, TYR124

-

28-HB (CID: 11038340)

-10.80

SER64, TYR102.

PHE57, PRO61, TYR102, HIS104, TYR124.

TYR63, TYR124

-

24-EBS (CID: 443055)

-11.23

TYR62, SER64, TYR102

PHE57, VAL59, PRO61, TYR63, TYR102, HIS104, TYR124

HIS104, TYR124

-

Hydroxychloroquine (CID: 3652)

-7.16

ASN58, VAL59.

PHE57, PRO61, TYR102, TYR124.

-

-

Table 3

Molecular docking analysis ofmain protease interactions with various ligands

Main protease protein (PDB:6lu7)

G-Score (Kcal/mole)

Hydrogen bounds

Hydrophobic interaction

π-Stacking

Halogen Bonds

BS (CID: 115196)

-10.36

GLY 143, SER144, CYS145, HIS163.

HIS41, MET49, MET165, ASP187.

HIS41.

-

28-HB (CID: 11038340)

-9.25

LEU141, HIS163, HIS164, GLU166.

HIS41, MET49, MET165, GLU166, GLN189.

-

-

24-EBS (CID: 443055)

-10.16

PHE140, GLY143, SER144, CYS145, HIS163.

HIS41, ASP187.

HIS41.

-

Hydroxychloroquine (CID: 3652)

-8.82

HIS41, GLU166, GLN189.

GLY143, SER144, CYS145, HIS164.

HIS41.

-

Table 4

Molecular docking analysis ofinterleukin-1 cytokine interactions with various ligands

Interleukin-1 protein (PDB:2nvh)

G-Score (Kcal/mole)

Hydrogen bounds

Hydrophobic interaction

π-Stacking

Salt Bridges

Brasslinolide (CID: 115196)

-9.02

THR79, LEU80, LEU134.

THR79, LEU80, PHE133.

-

LYS77

28-homocastasterone (CID: 11038340)

-13.83

TYR24, LEU26, VAL132.

GLU25, LYS77, LEU80, PHE133, LEU134.

-

-

24-epibrasslinolide (CID: 443055)

-9.97

PRO78, LEU134, GLY136, ASP142.

LYS77, LEU80, TRP120, PHE133.

LYS77.

-

Hydroxychloroquine (CID: 3652)

-4.99

LEU134, GLY136.

LYS77, TRP120, PHE133.

PHE133.

-

Table 5

Molecular docking analysis oftumor necrosis factor-α cytokine interactions with various ligands

Tumor necrosis factor-alpha protein (PDB:2nvh)

G-Score (Kcal/mole)

Hydrogen bounds

Hydrophobic interaction

π-Stacking

Salt Bridges

Brassinolide (CID: 115196)

-9.96

VAL17, ARG32, ALA33, ASN34, VAL91.

ALA18, ALA33, ASN34, SER147, GLY148, GLN149, VAL150.

-

-

28-homocastasterone (CID: 11038340)

-6.07

THR77, THR79, ILE97, GLU146.

LYS65, THR77, TYR115, ASP143, ALA145, GLU146, GLN149.

-

LYS65

24-epibrassinolide (CID: 443055)

-5.51

LYS65, THR77, TYR115, ASP143, ALA145, GLU146, GLN149.

THR77, ILE97, GLU146.

-

LYS65

Hydroxychloroquine (CID: 3652)

-7.51

ALA18, ASN34, SER147, GLY148, VAL150

VAL17, ASN34.

-

-

Table 6

Molecular docking analysis ofcyclooxygenase-2 interactions with various ligands

Cyclooxygenase-2 protein (PDB:4m1)

G-Score (Kcal/mole)

Hydrogen bounds

Hydrophobic interaction

π-Stacking

Salt Bridges

Brasslinolide (CID: 115196)

-12.25

ALA199, PHE200, ALA202, PHE210, TYR385, TRP387, LEU390, LEU391

28-homocastasterone (CID: 11038340)

-13.77

ALA199, PHE200, ALA202, TRP387, LEU390, LEU391.

ALA199, THR206, ASN382.

HIS207, HIS386, HIS388.

HIS207, HIS386.

24-epibrasslinolide (CID: 443055)

-13.22

ALA199, ASN382.

ALA199, PHE200, ALA202, PHE210, TYR385, TRP387, LEU390, LEU391.

HIS207, HIS386.

HIS207.

Hydroxychloroquine (CID: 3652)

-10.44

PHE310, TYR385, TRP387.

TYR148, THR206, THR212, ASN382, TRP387.

-

-

Table 7

Molecular docking analysis ofprostaglandin synthase interactions with various ligands

Prostaglandin synthase protein (PDB:1z9h)

G-Score (Kcal/mole)

Hydrogen bounds

Hydrophobic interaction

π-Stacking

Salt Bridges

Brasslinolide (CID: 115196)

-8.01

ASP164, SER165, GLU204, TYR209.

PHE112, TYR145, LYS 147.

LYS208

LYS208

28-homocastasterone (CID: 11038340)

-8.95

GLN161, ASN163, ASP164, THR187.

LEU162, VAL167, ILE186, TYR189, PHE205.

-

-

24-epibrasslinolide (CID: 443055)

-9.48

VAL148, ASP164, SER165, LYS208, TYR209.

PHE112, LYS147, ASP239.

PHE112, LYS147.

-

Hydroxychloroquine (CID: 3652)

-6.52

LYS115, VAL148, ASP164, SER166.

PHE112, LYS147, ASP239.

PHE112, LYS147

-

In the current in-silico study, N-terminal domain (NTD) were docked with the dietary phytosterols and HCQ. The docking scores of BS, 28-HC, and 24-EB against NTD of N-protein were -11.24 Kcal/mole, -10.80 Kcal/mole and -11.23 Kcal/mole, respectively compared to HCQ, that scored -7.16 Kcal/mole of binding energy with NTD of SARS-CoV2 (Table 2 and Figure 4). The BS forms hydrogen bond with VAL59, SER64 and ARG164 amino acid residues, while it forms hydrophobic interaction with PHE57, PRO61, HIS104, and TYR124 amino acid residues of NTD of N-protein. Likewise, 28-HC forms hydrogen bond with SER64 and TYR102 residues, while it forms hydrophobic interaction with PHE57, PRO61, TYR102, HIS104, and TYR124 amino acid residues of NTD of N-protein (Table 2). Both BS and 28-HC ligands form π-stacking interaction with TYR63, TYR124 amino acids residues. Correspondingly, 24-EB forms hydrogen bond with TYR62, SER64, and TYR102 amino acid residues and hydrophobic interaction with PHE57, VAL59, PRO61, TYR63, TYR102, HIS104, and TYR124 residues of NTD of N-protein (Table 2 and Figure 4). Hence, π-stacking interaction of N-protein’s NTD with TYR63 and TYR124 amino acids residues is common to all the three phytosterols. Docking studies between HCQ ligand and NTD of N-protein revealed hydrogen bonding between HCQ and amino acid residues ASN58 and VAL59, while PHE57, PRO61, TYR102, TYR124 are involved in hydrophobic interaction with HCQ.

Figure 3

The molecular docking of SARS-CoV-2 Spike (S) protein S1 subunit interaction with ligands. (A, C, E, G) The putative binding site of BS and 24-BC 24-Epi and HOC on SARS-CoV2 protein 3-D ribbon structure of SARS-CoV2 protein. (B, D, F, H). The interacting amino acid residues of SARS-CoV-2 with BS, 24-BC, 24-Epi and HOC.

https://s3-us-west-2.amazonaws.com/typeset-prod-media-server/53370576-dd65-4f6b-b36c-201fe8af681dimage3.png
Figure 4

The molecular docking of SARS-CoV2 Nucleocapsid protein interaction with ligands. (A, C, E, G) The putative binding site of BS and 24-BC 24-Epi and HOC on SARS-CoV2 protein 3-D ribbon structure of SARS-CoV2 protein. (B, D, F, H). The interacting amino acid residues of SARS-CoV-2 with BS, 24-BC, 24-Epi and HOC.

https://s3-us-west-2.amazonaws.com/typeset-prod-media-server/53370576-dd65-4f6b-b36c-201fe8af681dimage4.png
Figure 5

The molecular docking of SARS-CoV2 Main protease protein interaction with ligands. (A, C, E, G) The putative binding site of BS and 24-BC 24-Epi and HOC on SARS-CoV2 protein 3-D ribbon structure of SARS-CoV2 protein. (B, D, F, H). The interacting amino acid residues of SARS-CoV-2 with BS, 24-BC, 24-Epi and HOC.

https://s3-us-west-2.amazonaws.com/typeset-prod-media-server/53370576-dd65-4f6b-b36c-201fe8af681dimage5.png

Docking scores of BS, 28-HC, and 24-EB against Mpro were -10.36, -9.25, and -10.16 Kcal/mole respectively (Table 3 and Figure 5), whereas HCQ scored -7.90 Kcal/mole of binding energy against Mpro. Among the three phytosterols, BS had the highest docking score and consequently highest binding affinity for Mpro, thereby forming hydrogen bond with GLY 143, SER144, CYS145 and HIS163 residues, while HIS41, MET49, MET165 and ASP187 residues were involved in hydrophobic interaction. It also formed π-Stacking bond formed with HIS41 amino acid of Mpro (Table 3 and Figure 5). In comparison to BS, HCQ formed hydrogen bonds with PHE140, GLY143, SER144, CYS145 and HIS163 residues, while HIS41, MET49, MET165 and ASP187 of Mpro formed hydrophobic interaction.

Docking scores of BS, 28-HC and 24-EB against pro-inflammatory protein, IL1, are -9.02 Kcal/mole, -13.83 Kcal/mole, and -9.97 Kcal/mole, respectively (Table 4). Whereas, the same series of phytosterols has a docking score of -9.96 Kcal/mole, -6.07 Kcal/mole, and -5.51 Kcal/mole respectively against TNF-α (Table 5). Docking scores of phytosterols against pro-inflammatory COX-2 are -12.25 Kcal/mole, -13.77 Kcal/mole, and -13.22 Kcal/mole (Table 6), whereas the same series of phytosterols when docked against PS are -8.01 Kcal/mole, -8.95 Kcal/mole, and -9.48 Kcal/mole (Table 7). Similarly, HCQ’s docking score against IL1 and TNF-α are -4.99 Kcal/mole and -7.51 Kcal/mole respectively (Table 4, Table 5). Consequently, HCQ’s docking score against pro-inflammatory enzymes, COX-2 and PS are -10.44 Kcal/mole and -6.52 Kcal/mole (Table 6, Table 7) respectively.

The common molecular interaction between BS, 28-HC, and 24-EB against IL1 and TNF-α amino acid residues are (LYS77, LEU80, PHE133, LEU134) and (LYS65, THR77, TYR115, ASP143, ALA145, GLU146, GLN149) respectively whereas HCQ does not show significant similarity of molecular interaction with two proteins. Similarly, common molecular interaction between BS, 28-HC, 24-EB against the two targets, COX2 and PS, include amino acid residues (ALA199, PHE200, ALA202, PHE210, TYR385, TRP387, LEU390, LEU391) and (PHE112, VAL148, ASP164, SER165, LYS208, TYR209) respectively.

Discussion

Phytosteroids are consumed as part of functional food and herb-based folk medicine. 12 Assimilation of phytosteroids into vital organs via blood circulation induces metabolic changes in mammalian cells that includes antihyperglycemic and alterations of pro-inflammatory cytokines including IL-1, TNF- α, and COX-2. 13, 14 Three-dimensional structure of IL-1 consists of open-ended barrel shape composed beta pleated strands with two IL-1 receptor binding sites. 15 3D-structure of TNF-α consists of three monomers to form a compact bell-shaped trimer. TNF-α monomer forms an elongated, antiparallel beta-pleated sheet with a jelly-roll topology. 16 COX2 converts arachidonate into prostaglandins and thromboxane, thus, contributing to acute and chronic inflammatory modulatory effects. 17

Structurally, SARS-CoV2 is a nanometre-sized entity ranging between 65-125 nm in diameter. It has non-segmented positive direction RNA of 26-32kbs in length with crown-shaped spikes on the outer surface. It contains four structural proteins, namely, S-protein, E-protein, M-protein, and N-protein. 18 S-protein is a type I transmembrane glycoprotein of approx. 150 kDa molecular weight and is composed of two subunits, namely, S1 and S2, protruding from the viral surface. The S1 subunit (Figure 2 a) binds to angiotensin 2 converting enzyme (ACE2) host cellular receptor through its receptor-binding domain (RBD), followed by conformational changes in the S2 subunit, which allows the fusion peptide to insert into the host cell membrane. 19 Similarly, N-protein (Figure 2 b), a major structural protein that plays various roles in the viral replication and forms a ribonucleoprotein complex with the viral RNA through its NTD. The SARS-CoV2 main protease (SARS-CoV2 Mpro) consist of three domains that processes polypeptide with the help of a catalytic dyad and the active site of the protease is located between the domains I and II (Figure 2 c). 20 The spikes are responsible for the attachment of the viruses on the host surface and their subsequent entry into the host cells. Similarly, the N-protein and main protease (Mpro) participates in packaging and multiplication of the virus. Disruption of virus entry into the host cells and virus replication machinery may be potential therapeutic targets. Development of S-protein, N-protein and Mpro-targeted drugs against coronaviruses may be probable therapeutic option. 21 So far, no specific drug against COVID 19 has been approved by the Food and Drug Administration Agency (USFDA). 22 Initiation of SARS-CoV2 infection requires successful penetration of the virus in host cell membrane barrier followed by modulation host cellular metabolic stimuluses that may recognise the virus via a number of physicochemical and cell-signalling methods. Host cellular modulation aids in virus multiplication by influencing the homeostatic regulatory processes that controls cell and tissue functions. It ensures equilibrium dynamics of factors such as rate limiting steps, circulatory partitioning, product inhibition, protein synthesis and gene expression that affects respiratory, hepatic, gastrointestinal and neurologic homeostatic processes.13

Natural compounds from plant sources are among the potential candidates that are currently being investigated for their influence on viral cellular metabolic processes, thus acquiring early translational knowledge of such compounds as potential therapeutic interventions against COVID-19. The current study was undertaken to identify probable phytochemical candidates that could bind to the spike S1 subunit of S-protein, N-protein and Mpro, all of which are understood to be potential drug-targets against SARS-CoV2. Among the various phytochemicals tested against these three protein targets, we identified three dietary phytosteroid molecules that could bind with the inhibitor binding domain of the target proteins, as compared to an anti-malarial drug, HCQ, which is currently being repurposed for managing COVID-19, due to lack of any standard drug. 23 HCQ as the choice of positive control standard drug may be debatable. However, our choice of HCQ was based on recent evidences of COVID-19 response to HCQ treatment. One of the recent reports published in The Lancet Respiratory Medicine Journal suggest that in absence of any significant alternative, clinical trials are being conducted on HCQ with some significant results, especially in critically ill COVID-19 patients with history of pneumonia. 24 Another report of experiments conducted on COVID-19 infected non-human primates suggest that the cytokine storm was significantly attenuated by the HCQ treatment. 25 The molecular mechanism of action of HCQ on inhibition of cytokine storm in non-COVID-19 patients, however, is well established, especially in rheumatoid arthritis and Systemic lupus erythematosus. 26, 27

Comparatively low binding affinity of HCQ with S1 may be due to sterically bulkier carbon chain which could disrupt the interaction of HCQ with the binding site, thus contributing to the significantly lower binding affinity. 28 Both phytosteroids, 28-HC and 24-EB, forms hydrogen bond with GLY981 and GLN984 of S1 subunit, while they both form of hydrophobic interactions with amino acid residues THR943, LEU944, GLN947, GLN984, THR988 of the S1 subunit protein. HCQ forms lesser number of hydrophobic interaction and hydrogen bonds with the amino acid residues of S1 subunit, thus, further suggesting the reason behind low binding affinity of HCQ in comparison to phytosteroids, 28-HC and 24-EB.

N-protein binds to the viral RNA genome, forming a helical ribonucleoprotein (RNP) or nucleocapsid structure complex. RNP maintains RNA conformation that is necessary for replication and transcription of viral genome, in addition to cell-cycle progression via actin reorganization. N-protein consists of three domains, namely, N-terminal RNA-binding domain (NTD), C-terminal dimerization domain (CTD) and a central Ser/Arg (SR)-rich linker region. 21 Previous studies suggests that the N- and C-terminal domains of the N proteins are responsible for RNA binding and oligomerization, respectively. 28

Mpro, a 306 amino acids enzyme, induces the spike protein that further helps in viral penetration. 29 The current study was undertaken to identify phytosterols as a potential drug candidate against COVID-19 that could significantly bind and inhibit Mpro enzymes. Except for 28-HC, all other ligands formed π-Stacking interaction with HIS-41 amino acid residue, indicating that the amino acid HIS-41 may play an important role in inhibitor-protein binding affinity, thus, displaying lower binding energy.

Figure 6

Schematic illustration of the proposed mechanism of action of phytosteroids against COVID-19 prophylaxis. Phytosteroids may have multiple levels of mechanism of action. They may prevent viral entry into the host by preventing binding of viral subunit-S1 with ACE2 receptor. mechanism in inhibition of virus’ host entry and multiplication. If the virus is able to bypass this inhibition, these phytosteroids may prevent viral multiplication by inhibiting protein synthesis and vial components assembling. In case, the second level of defence also fails, the virus may trigger upregulation and activity of pro-inflammatory cytokines resulting in cytokine storm. Phytosteroids may attenuate this cytokine storm as the last resort to overcome manifestations of SARS-CoV2 infection.

https://s3-us-west-2.amazonaws.com/typeset-prod-media-server/53370576-dd65-4f6b-b36c-201fe8af681dimage6.png

IL-1 IL-6, GCSF and TNF-α are some of the pro-inflammatory cytokines produced by activated macrophages. Low-dose corticosteroids treatment lowers expression levels of these pro-inflammatory markers, in addition to quenching of COX-2 synthesis. These biomarkers are otherwise abnormally elevated as part of the natural immune response and may have a role in severe lung damage associated with COVID-19. 30 Therefore, to overcome these acute and chronic inflammatory modulation effects, in-silico molecular docking studies of potential anti-inflammatory phytosteroids were undertaken to identify potential alternative therapeutic molecules against COVID-19. The docking scores of BS, 28-HC, 24-EB and HCQ against two pro-inflammatory markers and two pro-inflammatory enzyme targets have been compared. BS has the highest docking score against all four pro-inflammatory targets suggesting that BS could have superior anti-inflammatory activity as compared to other two phytosteroids and HCQ.

Inhibition of host cell entry and viral replication machinery may be ideal therapeutic targets, but a parallel defence mechanism via controlling of cytokine storm gives therapeutic edge at two levels against SARS-CoV2. In case, a virus escapes firstly line of defence and is able to replicate inside host cell, a second level of defence may be triggered which attenuates cytokine storm from getting triggered. Hence, we identified seven such target proteins, S-, N-, and Mpro proteins, IL-1, TNF-α, COX2 and PS that may be potentially targeted by these three phytosteroids. Hence, based on the above hypothesis and its corresponding in silico results, we proposed a plausible molecular mechanism of action of these phytosteroids wherein viral entry points and multiplication may be hindered. If at all the virus succeeds in multiplication and trigger cytokine storm, these steroids may counter inhibit the expressions and activity of the corresponding cytokines (Figure 6). In summary, SARS-CoV2-related proteins are involved in host cell entry and viral replication, while the cytokine proteins are involved in triggering cytokine storm. All the seven target proteins were docked against the three phytosterols, BS, 28-HC, and 24-EB in comparison to the current repurposed drug, HCQ, wherein the docking scores and in turn the binding affinity of BS is highest, thus, suggesting that BS may be an effective alternative to HCQ against COVID-19.

Conclusions

The present in-silico study suggests that BS is an effective alternative to HCQ since it may potentially inhibit COVID-19 at multiple levels including host entry, viral multiplication and attenuation of cytokine storm. The outcome of our in-silico data may be the basis for in vivo and in vitro studies against COVID-19, simultaneously, establishing the concept of reverse pharmacology.

Declarations

The authors declare that there are no conflicts of interest, whatsoever, among themselves.

Source of Funding

None.

Conflict of Interest

None.

Acknowledgments

The authors acknowledge the Department of Biochemistry and Molecular Biology, Pondicherry University, Puducherry for provided facilities to carry out in this work.

References

1 

C Wang P W Horby F G Hayden G F Gao A Novel Coronavirus Outbreak of Global Health ConcernLancet2020395102233018530194

3 

P Conti G Ronconi A Caraffa C Gallenga R Ross I Frydas Induction of Pro-Inflammatory Cytokines (IL-1 and IL-6) and Lung Inflammation by Coronavirus-19 (COVI-19 or SARS-CoV-2): Anti-Inflammatory StrategiesJ Biol Regul Homeost Agents202034232731

4 

W Zhang Y Zhao F Zhang Q Wang T Li Z Liu The Use of Anti-Inflammatory Drugs in the Treatment of People with Severe Coronavirus Disease 2019 (COVID-19): The Perspectives of Clinical Immunologists from ChinaClin Immunol2020214108393

5 

Y Wang W Jiang Q He C Wang B Wang P Zhou Low-Dose and Short-Term Application of Corticosteroid Treatment in Patients with Severe COVID-19 Pneumonia: Single-Center Experience from Wuhan, China.medRxiv20209374292302

6 

EM Turk S Fujioka H Seto Y Shimada S Takatsuto S Yoshida Inactivates Brassinosteroid Hormones: An Intersection between Photomorphogenesis and Plant Steroid Signal TransductionPlant Physiol20031334164353

7 

C Eugster G Rivara A Biglino R Cavallo P Gioannini B Forno Phytosterol Compounds Having Antiviral Efficacy. Panminerva Med19973911220

8 

MB Wachsman JA Ramírez LB Talarico LR Galagovsky CE Coto Antiviral Activity of Natural and Synthetic Brassinosteroids.Curr Med Chem: Anti-Infective Agents20043216379https://doi.org/10.2174/1568012043354026

9 

V Castilla M Larzábal NA Sgalippa MB Wachsman CE Coto Antiviral Mode of Action of a Synthetic Brassinosteroid against Junin Virus ReplicationAntiviral Res20056828895

10 

XY Meng HX Zhang M Mezei M Cui Molecular Docking: A Powerful Approach for Structure-Based Drug DiscoveryCurr Comput Aided-Drug Des20117214657

11 

DV Der Spoel E Lindahl B Hess G Groenhof AE Mark HJ Berendsen GROMACS: Fast, Flexible, and FreeJ Comput Chem2005261617018

12 

VN Zhabinskii NB Khripach VA Khripach Steroid Plant Hormones: Effects Outside Plant KingdomSteroids2015978797

13 

V Athithan K Srikumar 28-Homocastasterone down Regulates Blood Glucose, Cholesterol, Triglycerides, SREBP1c and Activates LxR Expression in Normal & Diabetic Male RatChem Biol Interact2017277820

14 

J Leiro D García JA Arranz R Delgado ML Sanmartín F Orallo An Anacardiaceae Preparation Reduces the Expression of Inflammation-Related Genes in Murine MacrophagesInt Immunopharmacol2004489911003

15 

SZ Ben-Sasson J Hu-Li J Quiel S Cauchetaux M Ratner I Shapira IL-1 Acts Directly on CD4 T Cells to Enhance Their Antigen-Driven Expansion and DifferentiationProc Natl Acad Sci200910617711924

16 

MJ Eck SR Sprang The Structure of Tumor Necrosis Factor-Alpha at 2.6 A Resolution. Implications for Receptor BindingJ Biol Chem19892642917595605

17 

E Ricciotti GA Fitzgerald Prostaglandins and InflammationArterioscler Thromb Vasc Biol20113159861000

18 

AR Fehr S Perlman An Overview of Their Replication and PathogenesisCoronaviruses: Methods and Protocols20151282123

19 

DS Hui EI Azhar TA Madani F Ntoumi R Kock O Dar The Continuing 2019-NCoV Epidemic Threat of Novel Coronaviruses to Global Health - The Latest 2019 Novel Coronavirus Outbreak in WuhanChina Int J Infect Dis2020912646

20 

D Forni R Cagliani M Clerici M Sironi Molecular Evolution of Human Coronavirus GenomesTrends Microbiol20172513548

21 

L Zhang D Lin X Sun U Curth C Drosten L Sauerhering Crystal Structure of SARS-CoV-2 Main Protease Provides a Basis for Design of Improved α-Ketoamide InhibitorsScience2020368648940912

22 

A Zumla JFW Chan EI Azhar DSC Hui KY Yuen Coronaviruses-Drug Discovery and Therapeutic OptionsNat Rev Drug Dis201615532747

23 

MS Khuroo Chloroquine and Hydroxychloroquine in Coronavirus Disease 2019 (COVID-19). Facts, Fiction and the Hype: A Critical AppraisalInt J Ant Agents2020563106101

24 

FS Taccone J Gorham JL Vincent Hydroxychloroquine in the Management of Critically Ill Patients with COVID-19: The Need for an Evidence Base. Lancet Respir Med202091301729

25 

P Maisonnasse J Guedj V Contreras S Behillil C Solas R Marlin Hydroxychloroquine Use against SARS-CoV-2 Infection in Non-Human PrimatesNature202058558487

26 

C Ponticelli G Moroni Hydroxychloroquine in Systemic Lupus Erythematosus (SLE)Drug Safety Eval20171614119

27 

TW Nystad AM Fenstad LI Havelin O Furnes AK Skredderstuen BT Fevang SAT0061 Reduction in Joint Replacement Surgery in Patients with Reumatoid ArthritisAnn Rheum Dis201435126113

28 

X Xue H Yu H Yang F Xue Z Wu W Shen Structures of Two Coronavirus Main Proteases: Implications for Substrate Binding and Antiviral Drug DesignJ Virol2008825251527

29 

JK Mille GR Whittaker Host Cell Entry of Middle East Respiratory Syndrome Coronavirus after Two-Step, Furin-Mediated Activation of the Spike ProteinProc Natl Acad Sci USA201411142152149

30 

E Auriel K Regev AD Korczyn Nonsteroidal Anti-Inflammatory Drugs Exposure and the Central Nervous SystemHandbook of Clinical Neurology201457784



jats-html.xsl


This is an Open Access (OA) journal, and articles are distributed under the terms of the Creative Commons Attribution-NonCommercial-ShareAlike 4.0 License, which allows others to remix, tweak, and build upon the work non-commercially, as long as appropriate credit is given and the new creations are licensed under the identical terms.

Article type

Original Article


Article page

162-170


Authors Details

Athithan Velan


Article History

Received : 18-09-2024

Accepted : 03-10-2024


Article Metrics


View Article As

 


Downlaod Files