In the present study, based on the importance of inhibition in reducing or stopping the activity of COVID19- ligand with the ability to inhibit the virus in the drug database has been used. For this purpose, and considering the importance of identifying ligand-receptor interactions, the effectiveness of each of these ligands on the related proteins has been evaluated during the molecular docking simulation process.
Abl kinase inhibitors
Abl kinases are non-receptor tyrosine kinases that are involved in various cellular processes, especially as mediators of viral infection and/or may be involved in T-cell signaling [14]. Imatinib is an Abl kinase inhibitor used to treat Philadelphia chromosome-positive chronic myelogenous leukemia (CML) and acute lymphocytic leukemia (ALL). This small-molecule inhibitor works by inhibiting Bcr-Abl tyrosine-kinase [15]. It has been shown, from previous studies, that an Abelson (Abl) kinase inhibitor, Imatinib, substantially decreases SARS-CoV and MERS-CoV viral titers. In the same study, they also investigated the bronchitis virus (IBV) to study the function of Abl kinase activity during coronavirus infection and found that Imatinib and two different Abl kinase inhibitors, GNF2 and GNF5, reduced IBV titers by blocking virus infection [1]. Previously, Imatinib was shown to block the entry of SARS-CoV or MERS-CoV S protein [1]. SARS-CoV-2 is highly homologous to SARS-CoV, so studying the effects of Abl kinase inhibitors on IBV, SARS-CoV and MERS-CoV may be useful in identifying the host cell pathways required for COVID-19 infection. It may also provide insights into potential strategies for the treatment of COVID-19 especially by using Imatinib. Virus–cell and cell–cell fusion induced by the coronavirus S protein has a very similar mechanism. Abl kinase activity plays a role in cytoskeletal rearrangement, regulating endothelial barrier, and junctional dynamics, and hence Abl kinase inhibitors might also be capable of interacting by interfering with the actin dynamics needed for virus – cell and cell – cell fusion in SARA-CoV-2. Furthermore, studies have indicated that Abl kinase regulates inflammatory signaling, NFκB signaling, and oxidant-induced epithelial cells injury caused by infection and ARDS which can be followed by COVID-19 [1, 16]. In epithelial cell injury H2O2 release leads to C-Abl activation and nuclear translocation. C-Abl inhibition by Imatinib increases the expression of antioxidant proteins such as catalase and glutathione peroxidase, which have been reduced due to oxidative stress. Therefore, treatment with Imatinib during ARDS may prevent the death of lung cells [16]. Coleman et al., investigated the Abl kinase inhibitors, including Imatinib, in SARS-CoV and MERS-CoV in vitro. In the early stages of infection, after internalization and endosomal trafficking, the anti-CoV activity of Imatinib is effected by inhibiting virion fusion in the endosomal membrane. Imatinib inhibits a step in virion replication before the genomic production of RNA. They also investigated the role of Abl2 in the replication of SARS-CoV and MERS-CoV. To knock down the Abl2 protein levels, siRNA was used. They demonstrated that Abl2 expression is essential for a productive viral replication and can be blocked by Imatinib [17]. Imatinib and its methane sulfonate derivative can be used in treating viral liver diseases, in particular viral hepatitis by inhibiting replication, transmission, or both, of hepatitis viruses or of other RNA viruses including respiratory syncytial virus, herpes virus, influenza virus, poxvirus, para influenza virus, rhinovirus, yellow fever virus, West Nile virus, and encephalitis virus in order to maintain or decrease RNA viral load. Although it is not meant to be limited to any particular mechanism of action or bound by definition, it is suspected that Imatinib's antiviral properties may be partly due to its ability to inhibit viral replication and transmission. Cellular signal transduction pathways are known to play an important role in viral infection, and cellular phosphorylation events during viral infection so they are required to effectively replicate and proliferate the virus. Several cellular signaling pathways, related to viral replication, have been investigated. Tyrosine kinase inhibitors can be used, for example, epidermal receptor growth factor (EGFR) inhibitors such as monoclonal antibodies and small-molecule inhibitors. EGFR inhibitors including monoclonal antibodies, such as IMC-C225 (Cetuximab), Trastuzumab (Herceptin), and others (ABX-EGF, EMD 72000), and tyrosine-kinase inhibitors, such as OSI-774 (Erlotinib, Tarceva), ZD1839 (Gefitinib, Iressa), and others (GW2016, CI-1033), can be used in combination therapy with Imatinib. Monoclonal antibodies can block extracellular ligand binding, but at the intracellular portion of the receptor, the small-molecule inhibitors may exert their effects to prevent tyrosine kinase phosphorylation and the activation of signal transduction pathways [18]. Consequently, combination therapy may be useful in COVID-19 treatment as well as for the other RNA viruses. According to the Hubei Anti-Cancer Association Chronic Myeloid Leukemia Standardized Management Collaboration Group research on 299 CML patients who responded optimally to anti-CML therapy using Imatinib and other tyrosine kinase inhibitors, 0.3% of patients were infected by SARS-CoV-2 and among those who failed to respond to CML treatment, 2% of patients were diagnosed with COVID-19. Therefore, patients who failed to receive an appropriate response to anti-CML therapy medications were more likely to get infected by SARS-CoV-2. Although more detailed clinical data and studies on the prevalence of COVID-19 in patients with CML is required, this idea may be consistent with several possibilities [19]. For instance, it has been shown that the total number of natural killer (NK) cells and Regulatory T cells was decreased markedly in patients with CML as well as COVID-19 infection [20, 21] while tyrosine kinase inhibitors are able to regulate the immune system by increasing the number of natural killer cells (NK) and Regulatory T cells [19]. In order to assess the potential effectiveness of Imatinib on COVID-19, by using molecular docking, we investigated the affinity and efficiency of Imatinib and possible intermediary proteins: Spike protein and RNA-dependent RNA polymerases (RdRps).According to the NCBI database and Basic Local Alignment Search Tool (BLAST), we found that SARS-CoV-2 spike protein (PDB: 6XR8) and SARS-CoV spike protein (PDB: 5 × 58) Query Cover is 76.16% with a GMQE of 0.80, which makes it promising to consider the proteins similar. The docking results are shown in Table 2; low energy indicates the optimum protein-ligand complexes. Accordingly, docking scores for SARS-CoV-2 and SARS CoV were − 9.6 and − 10.2 kcal/mol respectively which are low enough to show the appropriate protein-ligand complexes. Docking interactions of Imatinib based on docking studies are depicted in Fig. 1. Furthermore, the type of interaction with the number of active site amino acids is also considered effective. In terms of interactions, the presence of hydrogen bond interactions can be very important as they have critical contributions to the binding structures and binding free energies, although the van der Waals and Pi-interactions contributed to the stabilization of the binding structures. If these interactions take place in the active position of proteins, it will be much better and more desirable. Interaction Types and Amino Acids involved in the inhibition of proteins are shown in Table 6.
Table 2
The docking score of Imatinib to SARS-CoV-2 and SARS-CoV spike protein, SARS-CoV-2 RNA- dependent RNA polymerases (RdRps) (RMSD: 0.00)
Drug name
|
Drug bank ID
|
2D structure
|
Drug use
|
Binding affinity (total energy): SARS-CoV-2 spike protein
|
Binding affinity (total energy): SARS-CoV spike protein
|
Binding affinity (total energy): SARS-CoV-2 RNA-dependent RNA polymerases (RdRps)
|
Imatinib
|
DB00619
|
|
Anti-cancer
|
-9.6
|
-10.2
|
-8.1
|
Binding interactions of Imatinib and SARS-CoV-2 spike protein shows that Imatinib interacts by forming Hydrogen bond with residue SER B: 50 and SER B: 967. A pi-sigma interaction is also visible between the drug and amino acid THR B: 302 and Pi-Alkyl interactions with LEU C: 754 and CYS C: 760. The large number of Pi-sigma interactions which involves charge transfer and helps in intercalating the drug in the binding site of the receptor. Pi-alkyl bond also improves the hydrophobic interactions of the ligand in the binding pocket of the receptor. With regards to the van der Waals interactions, it should be mentioned that there are thirteen amino acids contributing with the ligand, including: LYS B: 304, HIS B: 49, ASN C: 764, CYS C: 738, ASP C: 737, VAL C: 736, SER C: 735, THR B: 315, GLN B: 314, SER B: 316, ARG B: 319, THR B: 274 and GLY C: 757. Also, Imatinib inhibits SARS-CoV-2 with an IC50 of 130 nM. But, although imatinib binds to the receptor-binding domain (RBD) of SARS-CoV-2 spike protein, it does not inhibit the spike RBD: ACE2 interaction, suggesting a Bcr-Abl kinase-mediated fusion inhibition mechanism is responsible for the inhibitory action [22]. Another small molecule tyrosine kinase inhibitor (TKI) is Abivertinib used in lung cancer treatment, targeting both mutant forms of the epidermal growth factor receptor (EGFR) and Bruton's tyrosine kinase (BTK). Abivertinib binds to the BTK receptor which results in receptor phosphorylation prevention. It also plays a powerful immunomodulatory role in vitro by inhibiting pro-inflammatory cytokine production that are associated with cytokine release syndrome (CRS) or cytokine storm and progression of COVID-19 infection such as IL-1beta, IL-6 and TNF-alpha in patients with acute respiratory distress syndrome (ARDS). It is worth noting that FDA clears Abivertinib for Phase 2 safety and efficacy study in hospitalized patients with moderateto severe COVID-19 (http://www.aceatherapeutics.com). This strongly indicates that the Abl kinase signaling pathway is a promising area to study for the development of antiviral therapies.
Janus Kinase Inhibitor (jak Inhibitor)
Janus kinase inhibitors are being used in cancer treatment and inflammatory diseases by inhibiting the activity of the Janus kinase family of enzymes (JAK1, JAK2, JAK3, TYK2), by interfering with the JAK-STAT signaling pathway [23]. Ruxolitinib is one of those inhibitors with an anti-inflammatory effect related to the inhibition of the release of cytokines. This drug is used for the treatment of myelofibrosis and polycythemia vera (PCV) [24]. Also, Ruxolitinib has been approved in Covid-19 patients with respiratory failure with no invasive assisted ventilation required [25]. Improvement in chest computed tomography and faster recovery from lymphopenia were seen in patients as well [26]. Another JAK inhibitor is Baricitinib used for rheumatoid arthritis. A randomized phase 2 trial for this drug has been licensed to the usual treatment of pneumonia in COVID-19 patients [25]. Transcription by interferon-activated JAK – STAT signaling pathway (mainly mediated by JAK1 and JAK2) contributes to the upregulation of several interferon-controlled genes which destroy viruses in infected cells rapidly. Many viruses have formed strategies to combat interferon effects by blocking their signaling pathways, and viral-encoded factors that antagonize the JAK – STAT pathway are important virulence determinants. So, baricitinib blocking of the JAK – STAT signal results in interferon-mediated antiviral response inhibition that has an impact on SARS-CoV-2 infection progression [27]. According to another cohort study, baricitinib in combination with remdesivir and hydroxychloroquine showed clinical improvement in patients [28]. Moreover, docking results between SARS-CoV-2 main protease and our two JAK inhibitors (Ruxolitinib and baricitinib) were − 6.2 and − 6.1 kcal/mol respectively (Table 3).
Table 3
The docking score of candidate inhibitors to SARS-CoV-2 main protease and SARS-CoV-2 spike protein.
Drug name
|
Drug bank ID
|
2D structure
|
Protein target
|
Drug use
|
Binding affinity (total energy)
|
Ruxolitinib
|
|
DB08877
|
SARS-CoV-2 main protease
|
treatment of intermediate or high-risk myelofibrosis
|
-6.2
|
baricitinib
|
|
DB11817
|
SARS-CoV-2 main protease
|
treatment of rheumatoid arthritis (RA)
|
-6.1
|
Sitagliptin
|
|
DB01261
|
SARS-CoV-2 spike protein
|
Treatment of diabetes mellitus type 2
|
-6.0
|
Dipeptidyl Peptidase 4 Inhibitors (dpp4i)
Besides the main viral entrance port, angiotensin converting enzyme 2 (ACE2), dipeptidyl peptidase 4 (DPP4) can be investigated as well. DPP4 is a type II trans membrane glycoprotein with its major role in glucose and insulin metabolism which is expressed in many tissues, such as the immune cells. Also, it plays a significant role in immune regulation by activating Tcells, modulating NF-jB pathway, and the expression of CD86. Dipeptidyl peptidase 4 inhibitors mainly Sitagliptin can be used to treat diabetes mellitus type 2. Moreover, it was identified as a functional receptor for the MERS-CoV spike protein and although SARS-CoV-2 spike protein does not necessarily need DPP4. In spite of Sitagliptin and SARS-CoV-2 spike protein docking result with the score of -6.0 kcal/mol and the possibility that it does not alter ACE2, the potential anti-inflammatory involvement of DPP4 inhibitors raises concerns about DPP4 modulation that might decrease the cytokine-mediated acute respiratory complications of COVID-19 infection [25, 29] (Table 3) .
Rna-dependent Rna Polymerase Inhibitors
RNA-dependent RNA polymerase (RdRP) is an enzyme that catalyzes the replication of RNA from an RNA template which is encoded in the genomes of all RNA viruses [30] including SARS-CoV-2. There are some drugs that are considered to be nucleotide analog inhibitors of RdRps. Remdesivir (RDV) is one of those investigational drugs that have a wide variety of antiviral activities against RNA viruses including coronaviruses [31]. Remdesivir suppresses viral replication and it was initially tested in clinical trials to prevent the 2014 Ebola outbreak. Later investigations indicated Remdesivir's ability to inhibit replication of coronavirus, including SARS-CoV-2 as well [32]. In another cohort study, medical progress was observed in 68% of patients taken to the hospital with severe COVID-19 treated with Remdesivir [33]. In Japan also Favipiravir, an antiviral drug that targets the influenza viral RNA-dependent RNA polymerase, has been used against SARS-Cov-2 [34]. Ribavirin and Galidesivir [3] are other recommended drugs that we studied using auto dock vina. Accordingly, Remdesivir (Fig. 1) has the best binding capability with the score of -9.0 kcal/mol. From our docking studies, with the better binding energy compounds, the identified active residues were Lys47, Ser784, Ser709, Tyr129, His133 and Thr141. The major interaction between Remdesivir and RdRP is characterized by hydrogen bonding between the oxygen with TYR A: 129. A Pi-Cation interaction of aromatic ring and LYS A: 780 and an Alkyl interaction with ALA A: 706 have been observed. The docked result of the displayed in Fig. 1 indicates the drugs has eight hydrogen bond interaction with six amino acids reported on Table 6. Other important interactions such as alkyl, Pi-Cation interactions were also reported Table 6. Remdesivir high affinity has also been correlated with the existence of Van der waal forces formed on the amide substituents backbone with the respective amino acids GLY A: 774,ASP A:135, ALA A: 46. THR A: 710, LYS A: 714, GLN A: 773, ASN A: 705 AND SER A: 784, which established a strong cohesive environment, thus stabilizing the formed complex. It exhibits effective in vitro activity against SARS-CoV-2 with an EC50 at 48 h of 0.77 µM in Vero E6 cells [5]. On the other hand, Favipiravir had the least affinity to RdRps despite the recorded efficiency based on clinical trials (Table 4).
Table 4
The docking score of candidate inhibitors to COVID-19 RNA-dependent RNA polymerases (RdRps)
Drug name
|
Drug bank ID
|
2D structure
|
Drug use
|
Binding affinity (total energy)
|
Remdesivir
|
DB14761
|
|
Antiviral
|
-9.0
|
Galidesivir
|
DB11676
|
|
Antiviral
|
-6.9
|
Ribavirin
|
DB00811
|
|
Antiviral
|
-6.3
|
Favipiravir
|
DB12466
|
|
Antiviral
|
-5.4
|
Papain-like Protease Inhibitors
Papain-like Protease (PLpro) is characterized in different coronaviruses, including SARS and MERS [35]. Also, the genome of SARS-CoV-2 encodes for different proteins including PLpro[36]. The SARS-CoV PLpro and SARS-CoV-2 PLpro protein sequences are similar, so protease inhibitors that have shown efficacy against SARS-CoV might be similarly effective against SARS-CoV-2. Papain-like protease (PLpro) has a crucial role in the viral life-cycle [37, 38]. Targeting PLpro with antiviral drugs may result in viral replication blockage and the deregulation of signaling cascades in infected cells inhibition [39]. Consequently, anti-HCV drugs (Telaprevir, Grazoprevir, and Boceprevir) that bind to the SARS-CoV-2 PLpro active site (contained residues Asp164, Val165, Arg166, Glu167, Met 208, Ala246, Pro247, Pro248, Tyr 264, Gly266, Asn267, Tyr 268, Gln269, Cys217, Gly271, Tyr273, Thr301 and Asp302), may therefore oppose viral replication [35, 40]. Similarly, our study showed that Telaprevir (Fig. 1), Grazoprevir, and Boceprevir (HCV protease inhibitor) may be effective in binding to SARS-CoV-2 papain-like protease (PLpro) active sites to prevent viral replication (Table 5)[41]. Regarding to the lowest binding energy, the best ligand was Telaprevir with the score of -9.9 kcal/mol. The results of docking analysis (Table 6) showed that Telaprevir forms Hydrogen bonds with the 6W9C amino acids LYS B:92, HIS B:89 and THR A:74. Also there can be seen that the ligand interacts with LYS B: 92 via Pi-Cation and TYR A: 171 and ILE B: 44 via Pi-Alkyl interactions.
Table 5
The docking score of candidate HCV protease inhibitors to SARS-CoV-2 Papain-like protease (PLpro)
Drug name
|
Drug bank ID
|
2D structure
|
Drug use
|
Binding affinity (total energy)
|
Telaprevir
|
DB05521
|
|
Anti-Hepatitis C
|
-9.9
|
Grazoprevir
|
DB11575
|
|
Anti-Hepatitis C
|
-8.7
|
Boceprevir
|
DB08873
|
|
Anti-Hepatitis C
|
-7.8
|
Table 6
Interaction types and Amino acids involved in the Inhibition of PDB: 6XR8, PDB: 6M71
Ligand
|
Protein
|
Conventional Hydrogen Bond
|
Carbon hydrogen bond
|
Pi-sigma and amide interaction
|
Alkyl Interaction
|
Pi- Cation interaction
|
Imatinib
|
PDB: 6XR8
|
SER B:50
SER B: 967
|
THR C: 761
|
THR B: 302
|
CYS C: 760
LEU C: 754
|
-
|
Remdesivir
|
PDB: 6M71
|
ASN A: 781
HIS A: 133
SER A: 709
TYR A: 129
LYS A: 47
ASP A: 711
|
SER A: 709
|
-
|
ALA A: 706
|
LYS A: 780
|
Telaprevir
|
PDB: 6W9C
|
LYS B: 92
|
HIS B: 89
THR A: 74
|
-
|
ILE B: 44
TYR A: 171
|
LYS B: 92
|
And PDB: 6W9C with the top three drug candidate. |