2021 Volume 18 Pages 226-240
More than one and half years have passed, as of August 2021, since the COVID-19 caused by the novel coronavirus named SARS-CoV-2 emerged in 2019. While the recent success of vaccine developments likely reduces the severe cases, there is still a strong requirement of safety and effective therapeutic drugs for overcoming the unprecedented situation. Here we review the recent progress and the status of the drug discovery against COVID-19 with emphasizing a structure-based perspective. Structural data regarding the SARS-CoV-2 proteome has been rapidly accumulated in the Protein Data Bank, and up to 68% of the total amino acid residues encoded in the genome were covered by the structural data. Despite a global effort of in silico and in vitro screenings for drug repurposing, there is only a limited number of drugs had been successfully authorized by drug regulation organizations. Although many approved drugs and natural compounds, which exhibited antiviral activity in vitro, were considered potential drugs against COVID-19, a further multidisciplinary investigation is required for understanding the mechanisms underlying the antiviral effects of the drugs.
There is still a strong requirement of effective therapeutics for overcoming the COVID-19 pandemic. Up to 68% of the total amino acid residues encoded in the SARS-CoV-2 genome have been currently covered by the structural data. The recent activities in drug discovery against COVID-19 are reviewed in this article with emphasizing a protein structure-based approaches.
The novel coronavirus SARS-CoV-2, which emerged in 2019 in China and rapidly spread worldwide, has been causing the COVID-19 pandemic [1]. As of July 2021, there were more than 189 million confirmed cases and 4 million deaths globally. While several vaccines have been developed and the vaccination has been begun worldwide, there remains a great threat to public health and an urgent need for safety and effective therapeutic development [2]. More than one and half years since the pandemic situation, global efforts for drug discovery for combat COVID-19 are continuing. In this review, we overview the recent progress in structural analyses of the SARS-CoV-2 proteins and structure-based drug repurposing against COVID-19.
SARS-Cov-2 is an enveloped (+) single-stranded RNA coronavirus, and its genome size is ~29.9 kb in length with a 7-methyl-G 5' cap, a 3' poly-A tail, and more than 14 open reading frames (ORFs) as shown in Figure 1A. The orf1ab encodes a polyprotein of non-structural proteins, which were proteolytically processed by self-encoded two cysteine proteases, namely Papain-like protease (PLpro) in Nsp3 and main protease (Mpro or 3CLpro), into distinct 16 non-structural proteins (Nsp1 to Nsp16) essential for viral replication. The genome encodes four structural proteins, namely, spike glycoprotein (S protein), membrane (M) protein, envelope (E) protein, and nucleocapsid phosphoprotein (N protein). The genome also has at least nine ORFs for the accessory proteins (3a, 3b, 6, 7a, 7b, 8, 9b, 9c, and 10).
The genome structure of SARS-CoV-2 and map of 3D assignment. A: The genomic structure and coding proteins in the SARS-CoV-2 RNA genome. The coding proteins are colored by categories as; non-structural polyprotein (pale-blue), structural proteins (orange), and accessory proteins (green). B: Coverage of each coding region of the genome by the 3D structures deposited in the PDB. The color indicates the different types of proteins as; SARS-CoV-2 proteins (red), SARS-CoV or homologous proteins (orange), Model structures predicted by AlphaFold or CASP_commons (blue). C: The coverage depth of each genomic region by the 3D structures.
The non-structural proteins compose the viral replication and transcription complex (RTC) that are responsible for the RNA-processing, modifying, and proofreading functions maintaining the integrity of the large viral genome [3]. In addition to the two proteases, several non-structural proteins, namely, Nsp12, Nsp13, Nsp14, Nsp15, and Nsp16, have enzymatic activities for the viral replication, thus, were considered as drug targets. Nsp12 is the centerpiece of the RTC machinery and acts as RNA-dependent RNA polymerase (RdRp), which is the target protein of the approved drug Remdesivir. Nsp13, which encodes an RNA helicase, utilizes the energy derived from the hydrolysis of nucleoside triphosphates (NTPs) to unwind double-stranded (ds) RNA. Nsp14 has two catalytic domains where the N-terminus is 3'-to-5' exoribonuclease for RNA-proofreading, and the C-terminus is N7-methyltransferase (N7-MTase) for RNA modification processes, respectively. Nsp15 is a uridine-specific endoribonuclease with a C-terminal catalytic domain belonging to the EndoU nuclease family, which is highly conserved in coronaviruses. Nsp16 acts as a 2'-O-methyltransferase catalyzing the final step of the mRNA maturation where it catalyzes methylation at the ribose 2'-O position of the first nucleotide of the RNA with S-adenosyl-L-methionine (SAM) as the methyl group donor. Nsp10 plays a pivotal role for Nsp14 and Nsp16 enzymatic activities upon binding to the proteins to ensure structural stability.
The interaction of S protein, through the receptor-binding domain (RBD), to the cellular entry receptor angiotensin-converting enzyme 2 (ACE2) is the mandatory step of SARS-CoV-2 infection. Therefore, finding molecules that effectively inhibit the interaction between the virus-host proteins, including antibodies, is thought to be most promising in drug development.
Functions of the accessory proteins are largely unknown, but some of them are considered to be involved in evasion from host immune response [3–5].
When we began to study the knowledge-based structural modeling for SARS-CoV-2 proteins in January 2020, there was no available 3D structure from the virus. But the structure data of SARS-CoV, which caused the SARS outbreak of 2003, were available in the Protein Data Bank (PDB) [6,7]. Hence, we utilized those structures as the template for homology modeling (Figure 1B). We obtained structural models for 17 proteins out of a total 26 proteins (the total number was based on the initial annotation in the reference genome of SARS-CoV-2 isolated in Wuhan city, China [8]. Note that additional 3 ORFs, namely 3b, 9b and 9c, were also considered to be expressed currently).
Since the structure of SARS-CoV-2 3CLpro, the first 3D structure of the virus, was determined by X-ray crystallography in February 2020, the 3D structure data of SARS-CoV-2 proteins have been rapidly accumulated in the PDB. According to our surveillance of the database in June 2021, more than 1300 entries had been currently released in the database, including redundant entries which cover the same protein regions, and 68% of the amino acid residues of SARS-CoV-2 were covered by at least one 3D structure of its own protein during past one and half years. Figure 1C shows the accumulation of PDB entries that cover each genome region. The largest number of structures was deposited for the RBD domain, followed by 3CLpro and Macrodomain of Nsp3, indicating the global focus on these particular domains for drug development [9,10].
Owing to the breakthrough in cryo-electron microscopy (cryo-EM) techniques, many large protein complex structures have been determined to atomic or semi-atomic resolutions. For the RTC and the trimer of the S protein of SARS-CoV-2, 3D structures were also revealed by cryo-EM method.
Computational approaches for protein structure predictions also contributed to the completion of the SARS-CoV-2 protein structures. In March 2020, the community of protein structure prediction (CASP; Critical Assessment of Techniques for Protein Structure Prediction) had launched the SARS-CoV-2 structure modeling initiative for a challenge to predict the 3D structures of the ten viral proteins and domains, namely, Nsp2, the C-terminal domain of Nsp3, Nsp4, Nsp6, ORF3a, M protein, ORF6, ORF7b, ORF8, and ORF10, for which no apparent template structure was available (https://predictioncenter.org/caspcommons/). In these challenges, about 40 participants, including groups from Japan, deposited their modeled structures of the target proteins to the CASP site by the designated deadline (all the challenges were finished by the end of May 2020). After the challenge, the structures of ORF3a and ORF8 have been experimentally determined. Particularly noteworthy was that the model structure of ORF3a deposited by AlphaFold, which utilized a deep-learning approach developed by DeepMind [11] and demonstrated startling accuracy at CASP14, showed very high similarity to that of the experimentally determined one. The other protein models by AlphaFold also assisted in the structure determinations of ORF8 and Nsp2 proteins by X-ray crystallography and cryo-EM methods, respectively, indicating that the protein structure prediction with machine learning opened the door toward the next generation of structural biology. They also provided the model structures of Nsp4, Nsp6, and M protein, of which 3D structures were not determined yet, at the website (https://deepmind.com/research/open-source/computational-predictions-of-protein-structures-associated-with-COVID-19).
When we built knowledge-based models of SARS-CoV-2 proteins with approved drugs, only a small set of structural data of SARS-CoV or other related homologous proteins bound to small compounds or ligands were available [12]. Here we re-examined the PDB for the SARS-CoV-2 proteins and found that some novel structures in complex with drugs/compounds were deposited as summarized in Table 1. It clearly demonstrated the consequence of the increasing interest in the mechanism of action of potential drugs against the SARS-CoV-2 proteins, as described in detail below.
Drug | KEGG | DrugBank | Target | PDB codes | Clinical Status |
---|---|---|---|---|---|
Remdesivir | D11472 | DB14761 | Nsp12 (RdRp) | 7c2k, 7bv2, 7b3b, 7b3c, 7l1f | EUA for COVID-19 |
Favipiravir | D09537 | DB12466 | Nsp12 (RdRp) | 7aap, 7ctt | Phase 3 for COVID-19 |
Suramina) | D00808 | DB04786 | Nsp12 (RdRp) | 7d4f | Phase 2 for acute kidney injury |
Bamlanivimab | D11936 | DB15718 | S | 7l3n | EUA for COVID-19 |
Etesevimab | D11944 | DB15897 | S | 7c01 | EUA for COVID-19 |
Casirivimab | D11938 | DB15941 | S | 6xdg | EUA for COVID-19 |
Imdevimab | D11939 | DB15940 | S | 6xdg | EUA for COVID-19 |
PF-00835231 | None | None | 3CLpro | 6xhm | Phase 1 |
GC-376a) | None | None | 3CLpro | 6wtj, 7c8u, 7cb7, 7cbt, 7jsu, 7k0f | n.a. |
Boceprevira) | D08876 | DB08873 | 3CLpro | 6wnp, 6xqu, 6zru, 7brp, 7c6s, 7com, 7k40 | Approved for HCV |
Telaprevir | D09012 | DB05521 | 3CLpro | 6xqs, 6zrt, 7c7p, 7k6d, 7k6e, 7lb7 | Approved for HCV |
Baicaleina) | C10023 | DB16101 | 3CLpro | 6m2n | Phase 2 for influenza infection |
Carmofura) | D01784 | DB09010 | 3CLpro | 7buy | n.a. |
Shikonina) | C17412 | None | 3CLpro | 7ca8 | n.a. |
Ebselena) | None | DB12610 | 3CLpro | 7bak, 7bal, 7bfb | Phase 2 for COVID-19 |
Myricetina) | C10107 | DB02375 | 3CLpro | 7b3e | n.a. |
Ascorbic acid | D00018 | DB00126 | 3CLpro | 7mpb | Phase 3 for COVID-19 |
Ebselena) | None | DB12610 | Nsp3 (PLpro) | 7m1y | Phase 2 for COVID-19 |
GRL0617 | None | DB08656 | Nsp3 (PLpro) | 7jrn, 7cjm, 7cmd, 7jir | n.a. |
GS-441524 | C22275 | DB15686 | Nsp3 (PLpro) | 7bf6 | Phase 1 for COVID-19 |
Tipiracil | D10467 | DB09343 | Nsp15 | 6wxc | Approved for colorectal cancer |
Sinefungin | D05846 | DB01910 | Nsp16 | 6wkq, 6yz1 | n.a. |
a) The inhibitory activity was measured and confirmed by high-throughput in vitro assays. EUA: Emergency used authorization for COVID-19 by FDA.
Main protease 3CLpro is one of the most-studied proteins for drug discovery since the past SARS outbreak in 2003. Therefore, many efforts have been made for structural analyses of the protein and its potential inhibitors. After the first 3D determination of 3CLpro of SARS-CoV-2, more than 300 structures were determined by X-ray crystallography. Because of the demand for understanding the inhibitory mechanism, many of the crystal structures were determined in complex with peptide mimetic compounds, natural compounds, or others.
The following inhibitor compounds were co-crystalized with 3CLpro: PF-00835231 (see Table 1 for PDB ID) was originally developed more than 15 years ago against 3CLpro of SARS-CoV [13]. GC-376 was initially developed for a broad-spectrum inhibitor of 3C proteases in Picornaviruses, Noroviruses, and Coronaviruses [14]. Boceprevir and Telaprevir were originally developed as Hepatitis C virus (HCV) protease inhibitors and approved by the Food and Drug Administration (FDA) in U.S.A. These peptide-mimetic drugs, all of which were covalent inhibitors, exhibited high inhibitory activity for the 3CLpro in SARS-CoV-2 in vitro [13,15] and were utilized as lead compounds for further drug development by increasing efficacy and bioavailability [16]. Natural compounds such as baicalein, myricetin, shikonin, and ascorbic acid also bound to 3CLpro in the crystals [17] as shown in Figure 2A. These compounds, except ascorbic acid, showed inhibitor potency by in vitro high-throughput screenings [18–20] but no clinical trial conduction was announced yet for them. Despite no direct evidence for the 3CLpro inhibitory activity, ascorbic acid (vitamin C) is evaluated in phase 3 clinical trial for COVID-19 because it has anti-inflammatory properties, which implies influence of SARS-CoV-2 infection [21].
Structures of SARS-CoV-2 proteases. (A) The crystal structures of 3CLpro complexed with compounds, baicalein (upper left, [PDB ID: 6m2n], myricetin (upper right, [PDB ID: 7b3e]), shikonin (lower left, [PDB ID: 7ca8]), ascorbic acid (lower right, [PDB ID: 7mpb]). The catalytic residues (H41 and C145) are represented as stick models. (B) Superposition of the crystal structures of PLpro domains of Nsp3 complexed with its inhibitor GRL0617 (color in beige, [PDB ID: 7jrn]) and complexed with ISG15 (color in blue, [PDB ID: 6xa9]). The catalytic residues are represented by stick models.
Another protease, PLpro in Nsp3, was also crucial for viral replication and evasion from host immunity. The three distinct substrates of PLpro, namely the viral polyprotein, degradative Lys48-polyubiquitin, and antiviral ISG15 signals by deubiquitination or deISGylation, make PLpro an excellent candidate for pharmacological intervention. One promising compound was GRL0617, originally discovered as a noncovalent inhibitor specifically for PLpro of SARS-CoV in 2008 [22]. The co-crystal structure of SARS-CoV-2 PLpro with GRL0617 revealed that the compound resided in a pocket in the palm region, apart from the catalytic site. When the structure was compared with another co-crystal structure of PLpro with human ISG15, GRL0617 was implied to make a steric clash to the C-terminus of ISG15, suggesting that the compound would block the C-terminus of ISG15 from binding to the catalytic site of PLpro (Figure 2B). An in vitro experiment demonstrated that GRL0617 inhibited the deubiquitination and deISG15ylation activity of PLpro, which indicated that GRL0617 was a noncovalent inhibitor of PLpro [23].
RdRp—Remdesivir/Favipiravir/SuraminRNA-dependent RNA polymerase (RdRp/Nsp12) is one of the major drug targets of SARS-CoV-2. Remdesivir (approved drug) and Favipiravir (evaluated in phase 3 clinical trial) target RdRp. Recent advances in the cryo-EM revealed the RdRp-drug interaction modes. Remdesivir is a prodrug converted to the active drug in the triphosphate form (remdesivir-triphosphate; RTP) within cells. The complex structure of RdRp-RNA with RTP revealed that when RTP was incorporated into a new RNA product, the incorporated Remdesivir moiety made a barrier to further RNA translocation by steric hindrance, consequently stalling the RNA polymerase reaction [24,25]. The structural details of the protein-drug interaction would contribute to developing more effective drugs compared to Remdesivir. Some of the new nucleotide-analog drugs, including oral drug Molnupiravir with higher bioavailability, entered the phase 2/3 clinical trial (ClinicalTrials.gov identifier NCT04575597).
Suramin, which has been used to treat African sleeping sickness caused by trypanosomes for about 100 years, is now evaluated in phase 2 clinical trial for acute kidney injury (ClinicalTrials.gov identifier NCT04496596). Recently suramin has been shown to inhibit infection of SARS-CoV-2 in cell culture by preventing cellular entry of the virus [26]. The cryo-EM structure of SARS-CoV-2 RdRp bound to suramin revealed two binding sites. One site directly blocked the binding of the RNA template strand, and the other site made clashes with the RNA primer strand near the RdRp catalytic site [27]. The same group also demonstrated that suramin and a suramin derivative were at least 20-fold more potent than Remdesivir by a biochemical assay.
Nsp15—TipiracilNsp15 is a uridine-specific endoribonuclease with the C-terminal catalytic domain belonging to the EndoU family that are highly conserved in coronaviruses. This enzyme acts in homo-hexamer. It is well known that uridylate-specific nucleolytic activity of Nsp15 on single-stranded and double-stranded (ds) RNA limits the formation of dsRNA intermediates. Thus Nsp15 compromises the ability of specific cytoplasmic viral RNA sensors to activate the IFN-I response of innate immunity to infection [28]. Although a biochemical screening identified NSC95397, which was validated as an inhibitor of the endoribonuclease, this compound did not inhibit SARS-CoV-2 growth in Vero E6 cells [29].
Recently, the structure of SARS-CoV-2 Nsp15 was determined by X-ray crystallography [30]. It formed hexamer and the overall structure resembled SARS-CoV Nsp15 because of the high sequence identity (95.7%) between the two proteins. The structure of SARS-CoV-2 Nsp15 was also determined with uracil analogs, including Tipiracil, which was a thymidine phosphorylase inhibitor approved for the treatment of refractory metastatic colorectal cancer [31]. Biochemical and whole-cell assays demonstrated that Tipiracil inhibited SARS-CoV-2 Nsp15 activity by interacting with the uridine binding pocket in the active site, suggesting that uracil and its derivatives might represent a plausible lead for nucleotide-like drug development [30].
Nsp16—SinefunginThe viral mRNA capping is essential for efficient viral protein production, protection of viral RNA from host degradation, and evasion from the host’s innate immune responses [32]. Therefore, the enzymes in the RNA capping are now considered one of the potential targets for drug development, while it has been less focused at the earlier stage of the pandemic. Nsp16 is an enzyme that acts in the final step of mRNA maturation as described in the previous section. Nsp16 is activated upon binding to Nsp10, which is also a cofactor of Nsp14, forming a heterodimer.
We previously built a homology model of the SARS-CoV-2 Nsp16-Nsp10 heterodimer using that of SARS-CoV as a template [33]. The template structure of SARS-CoV Nsp16 bound Sinefungin, a natural product from Streptomyces griseolus and experimentally used as antibiotics [34]. Utilizing this drug-protein complex structure, we built the model structure of SARS-CoV-2 Nsp16-Nsp10 heterodimer complexed with Sinefungin (Figure 3). Recently, the crystal structure of the Nsp16-Nsp10 heterodimer from SARS-CoV-2 was solved [35,36]. Similar to the case of SARS-CoV, the solved structure bound Sinefungin, and the complex structure is similar to our model [12]. However, because of the close structural similarity of Sinefungin to the natural metabolite SAM, it was thought to be a broad-spectrum inhibitor of methyltransferase, and inappropriate for therapeutics. Therefore attempts to design a compound that inhibits Nsp16 more specifically were reported [37].
Comparison of structures of Nsp10-Nsp16 complex. (A) The homology model of SARS-CoV-2 Nsp10 (beige) and Nsp16 (blue) complex bound to Sinefungin and m7GpppA, which were based on the structure of SARS-CoV as the template in our previous study (YP_009725311_v02_1.pdb from BSM-Arc [85] entry BSM-ID BSM00015). (B) The crystal structure of the Nsp10-Nsp16 complex from SARS-CoV-2 bound to Sinefungin [PDB ID: 6yz1].
The S protein is essential for the viral entry to host cells. In the viral entry process, the receptor-binding domain (RBD) interacts with the host cell membrane receptor, angiotensin-converting enzyme 2 (ACE2). The mutations in the RBD affect the virus-host interaction and some of the mutants are assigned as variant of concern (VOC), which would have higher infectivity and/or severer clinical consequence. The interest of many researchers on this target is how to interfere with the host-viral interaction, and the most direct solution is developing neutralizing antibodies. Hence, many structural data of S protein complexed with various antibodies were deposited to the PDB. The entries contain neutralizing antibodies approved for emergency use, such as Bamlanivimab [38], Etesevimab, Casirivimab, and Imdevimab (Figure 4). Because these antibodies do not entirely share the binding interfaces, the cocktail antibody drugs are expected to be effective for viral variants, which have accumulated several mutations.
Experimentally determined structures of the RBD of S proteins from SARS-CoV-2 complexed with the approved neutralizing antibodies. (A) The cryo-EM structure of RBD complexed with Bamlanivimab [PDB ID: 7l3n]. (B) The crystal structure of RBD complexed with Etesevimab [PDB ID: 7c01]. (C) The cryo-EM structure of RBD complexed with Casirivimab and Imdevimab [PDB ID: 6xdg]. The RBD and neutralizing antibodies are represented as surface model and ribbon model, respectively. The surfaces colored yellow in RBD depict the antibodies binding regions.
Once an atomic structure of the target protein becomes available to a high resolution, a computational molecular docking methodology enables us to virtually screen the candidate compounds that potentially inhibit or attenuate the protein functions from libraries of hundreds of thousands of compounds [39]. If an approved drug, the one with proven acceptable safety for humans, was found to be effective for COVID-19, it could make a shortcut to the early-stage of clinical trials. Therefore, many studies have employed virtual screening strategies to find potentially effective compounds from existing approved drugs for repurposing.
Many of the target viral proteins are essential enzymes for the viral replication as represented by proteases (3CLpro and PLpro) [40–47], RNA replicase (RdRp and Nsp13) [43,46,48–52], methyltransferases (Nsp14 and Nsp16) [53–56], and endoribonuclease (Nsp15) [57–59]. In the reported virtual screening studies listed in Table 2, DrugBank [60], ChEMBL [61], PubChem [62], and/or ZINC databases [63], which contained FDA-approved or clinical stage drugs, were used as compound libraries. Some studies also used databases containing natural compounds, medicinal plants, or phytochemicals such as KNApSAcK [64], MPD3 [65], IMPPAT [66], or NANPDB [67].
Target | Drug Libraries | Proposed potential drugs | Reference |
---|---|---|---|
3CLpro | Phytochemicals with antiviral activity from 11 plants retrieving from literature and PubChem | Quercetin 3-vicianoside, Absinthin, Delphinidin 3-O-glucoside, Chrysoeriol 8-C-glucoside, Piperolactam A | Joshi et al., 2020 [40] |
DrugBank | Proflavine, Chloroxine, Demexiptiline, Fluorouracil, Oteracil | Gao et al., 2020 [41] | |
FDA approved drugs from DrugBank, ZINC, Selleckchem, Enamine subsets, e-Drugs and BindingDB | Metyrapone, Rufinamide, Zonisamide, Lacosamide, Apatinibe | Federico et al., 2021 [42] | |
FDA approved drugs from ZINC | Tetracycline, Dihydroergotamine, Ergotamine, Dutasteride, Nelfinavir | Gul et al., 2020 [43] | |
FDA approved drugs from ChEMBL, DrugBank, DrugCentral, Selleckchem | ENMD-981693, Felypressin, Brilacidin, Ritonavir, Saquinavir | Cavasotto et al., 2021 [44] | |
FDA approved drugs from ChEMBL, DrugBank, ZINC, Selleckchem | Talampicillin, Lurasidone | Elmezayen et al., 2020 [45] | |
KNApSAcK | Caribine, Cryptopleurine, Justicidin D, Diphyllin, Quercetin | Gani et al., 2021 [46] | |
FDA approved drugs from DrugBank and ZINC | Ivermectin, Diosmin, Selinexor, Bromocriptine, Elbasvir, Dihydroergotamine | Yuce et al., 2021 [47] | |
PLpro | FDA approved drugs from ChEMBL, DrugBank, DrugCentral, Selleckchem | Anatibant, Pilaralisib, Tiracizine, Zabofloxacin, Picotamide | Cavasotto et al., 2021 [44] |
RBD of S protein |
FDA approved drugs from ChEMBL, DrugBank, DrugCentral, Selleckchem | Pralatrexate, Carumonam, Aclerasteride, Granotapide | Cavasotto et al., 2021 [44] |
KNApSAcK | (–)-beta-sitosterol, 10-Methoxycamptothecin, Inoxanthone, Emetine, Alpha-cedrene | Gani et al., 2021 [46] | |
Nsp12 (RdRp) |
FDA approved drugs from ZINC database | Eltrombopag, Tipranavir, Ergotamine, Conivaptan | Gul et al., 2020 [43] |
FDA approved drugs from e-Drugs3D | Quinupristin, Cetrorelix, Dactinomycin, Rifampin, Sirolimus | Pokhrel et al., 2020 [48] | |
ZINC15 database | Ivermectin, Rifabutin, Rifapentine, Fidaxomicin, 7-methyl-guanosine-5'-triphosphate-5'-guanosine | Parvez et al., 2020 [49] | |
KNApSAcK | Justicidin D, 10-Methoxycamptothecin, Inoxanthone, 3-O-Caffeoylquinic acid | Gani et al., 2021 [46] | |
Phytochemicals from Indian medicinal plants | Cordifolide A, Sitoindoside IX | Koulgi et al., 2021 [50] | |
Nsp13 | Natural compounds from MPD3 database | ZINC257223845 | Ahmad S. et al., 2021 [51] |
Natural compounds from IMPPAT database | (+)-Epiexcelsin, Euphorbetin, Isorhoeadine, Picrasidine M, Picrasidine N | Vivek-Ananth et al., 2021 [52] | |
Nsp14 (N7-MTase) |
Compounds from Traditional Chinese Medicine database | TCM20111, TCM31007, TCM3495, TCM5376, TCM57025 | Selvaraj et al., 2021 [53] |
FDA approved, worked-out-FDA or investigational-only drugs from ZINC database | Hypericin, Olysio, Sovaprevir, Celsentri, Saquinavir | Liu et al., 2021 [54] | |
Nsp14 (ExoN) |
FDA approved, worked-out-FDA or investigational-only drugs from ZINC database | Hypericin, Bromocriptine, Tanespimycin, Idarubicin, Emend | Liu et al., 2021 [54] |
Nsp15 | FDA approved drugs from ZINC | Citrate, Dihydroergotamine, Ergotamine, Glisoxepine, Idarubicin | Chandra et al., 2020 [57] |
FDA approved drugs from DrugBank | Elbasvir, Paritaprevir | Sixto-Lopez Y et al., 2021 [58] | |
Selleckchem Natural compound library | Oleuropein, Thymopentin | Vijayan R & Gourinath S 2021 [59] | |
Nsp16 | FDA approved drugs from the PyRx 0.8 virtual screening tool | Digitoxin, Dihydroergotamine, Irinotecan, Sinefungin, Teniposide | Sharma K et al., 2020 [55] |
Natural compounds from North African Natural Products database | Citrinamide A, 4,5-Di-p-trans-coumaroylquinic acid, Genkwanin-6-C-beta-glucopyranoside, Paraliane diterpene | Mohammad A et al., 2021 [56] |
All the target proteins except for PLpro were examined in two or more independent studies. Intriguingly, the drugs raised by each study for the same target showed limited overlap with each other. It might be due to the differences in employed compound libraries, docking tools, or different atomic coordinates of the target proteins. On the other hand, several compounds were commonly raised in the independent studies; for example, Dihydroergotamine, originally discovered as a serotonin receptor antagonist and was approved for aborting or preventing vascular headaches, was commonly found as a potential 3CLpro inhibitor in two independent studies [43,47]. Other examples were quercetin and its derivative, which were sort of flavonoids and previously reported as the candidates of 3CLpro inhibitor for SARS-CoV [68].
Some compounds were commonly raised as candidates for multiple target proteins of SARS-CoV-2. For instance, Elbasvir, which was originally discovered as an HCV NS5A inhibitor, was selected as the candidate for 3CLpro and Nsp15. Also, Ergotamine, a similar drug as Dihydroergotamine, was selected as a drug for 3CLpro, Nsp15, and RdRp in different studies, suggesting that these drugs have the potency to bind several target proteins.
In general, the drug candidates raised by virtual screening should also be verified experimentally via in vitro biochemical or cell-based assays. To our best knowledge, no drug of COVID-19 found only in the virtual-screening studies has successfully passed to clinical trials to date.
Several independent groups performed in vitro biochemical assays for screening compounds from chemical libraries that can inhibit each viral enzyme. To identify drug-repurposing candidates, most of the studies used chemical libraries containing FDA-approved or late-stage clinical trial drugs provided by non-profit organizations such as Drug Repurposing Hub [69] and ReFRAME library [70], or purchased from commercial suppliers. These in vitro assays targeted mainly seven enzymatic activities of six SARS-CoV-2 proteins, namely, 3CLpro (protease) [15,17,19], PLpro (protease) [71], RdRp (RNA replicase) [72], Nsp13 (helicase) [73], Nsp14 (exonuclease) [74], Nsp14 (N7-MTase) [75], and Nsp15 (endonuclease) [29,76] (Table 3).
Target | Drug Libraries | Proposed candidate | Reference |
---|---|---|---|
3CLpro | ~10,000 compounds of approved or clinical used drugs and natural products from commercial libraries e.g. Target Mol, Selleck, Shanghai Institute for Advanced Immunochemical Studies | Ebselen, Disulfiram, Tideglusib, Carmofur Shikonin | Jin et al., 2020 [17] |
18 protease inhibitors including HCV or HIV proteases | Boceprevir, GC-376 | Fu et al., 2020 [15] | |
Various flavonoids | Baicalin, Herbacetin, Pectolinarin | Jo S. et al., 2020 [20] | |
Dompe “Safe-In-Man” proprietary collection (containing 607 drug candidates), EU-OPENSCREEN, DrugRepurposingHub) | Thioguanosine, MG-132, Myricetin | Kuzikov et al., 2021 [19] | |
PL pro | 5,576 compounds (3,727 approved drugs and late-stage clinical drug candidates) | No candidate found | Klemm T et al., 2020 [71] |
RdRp | ~5000 Commercial libraries (Sigma, Selleck, Enzo, Tocris, Calbiochem, Symansis) | GSK-650394, C-646, BH3I-1, Suramin, Cefsulodin | Bertolin et al., 2021 [72] |
Nsp13 | ~5000 Commercial libraries (Sigma, Selleck, Enzo, Tocris, Calbiochem, Symansis) | FPA124, Suramin | Zeng et al., 2021 [73] |
Nsp14 (N7-MTase) |
~5000 Commercial libraries (Sigma, Selleck, Enzo, Tocris, Calbiochem, Symansis) | Trifluperidol, PF-03882845, Inauhzin, Lomeguatrib | Basu et al., 2021 [75] |
Nsp14 (ExoN) | ~5000 Commercial libraries (Sigma, Selleck, Enzo, Tocris, Calbiochem, Symansis) | Patulin, Aurintricarboxylic acid | Canal et al., 2021 [74] |
Nsp15 | ~5000 Commercial libraries (Sigma, Selleck, Enzo, Tocris, Calbiochem, Symansis) | NSC95397 | Canal et al., 2021 [29] |
ReFRAME drug repurposing library | Exebryl-1 | Choi R et al., 2021 [76] |
Although both virtual-screening and in vitro assay approaches utilized similar drug libraries in the cases of 3CLpro, most of the compounds raised by the virtual screening studies were not re-discovered by in vitro high-throughput inhibitory assays. Exceptionally, compounds with flavonoid backbone, e.g., baicalein, quercetin, myricetin, and baicalin, were commonly observed in both lists. This result might be partly due to the limited number of drugs/compounds listed as candidates in the computational studies.
Among the drugs discovered by those in vitro assays, we confirmed that Ebselen and Disulfiram, both exhibited inhibitory activities against 3CLpro, have been entered into phase 2 clinical trials for COVID-19 (ClinicalTrials.gov identifiers NCT04484025 and NCT04485130, respectively).
The cell-based assay is most widely used for the drug screening because of its high-applicability to throughput assays compared to biochemical approaches. Usually, the cytopathogenic effect (CPE), the changes in host cell morphology caused by the target infecting virus, is measured as infection level in this method. The antiviral activity can be evaluated by measuring a change of the degree of CPE with/without a compound. At the earlier stage of the COVID-19 pandemic, many approved drugs, including Chloroquine/Hydroxychloroquine [77], Ivermectin [78], Nelfinavil [79] and Cepharanthine [79,80], were identified as potent anti-COVID-19 drugs by the cell-based assays. Many of these drugs were entered into clinical trials including ongoing ones. The FDA granted an emergency use authorization (EUA) for chloroquine and hydroxychloroquine for the treatment of hospitalized patients with COVID-19 in March 2020. However, the emergency use was revoked in June based on the follow-up assessment of EUA drug and the clinical trial results, in which it was demonstrated that the drugs did not show any effectiveness in reducing either mortality or morbidity [81]. For Ivermectin, many clinical trials are ongoing in several countries including Japan, while WHO has recommended its limited use within clinical trials because of less than sufficient evidence for the drug efficacy against COVID-19 [82]. Nelfinavir, which is a viral protease inhibitor and approved for HIV infection, has been entered into a clinical trial in Japan [83].
With the cell-based assays, compounds can be evaluated under the presence of human and viral genome/proteome, which endorses promising feasibility for the drug candidates. On the other hand, due to the nature of this method, target proteins and inhibition mechanisms of the identified potential drugs in host cells tend to be uncertain. A recent report demonstrated that some antiviral drugs accumulated in an intracellular compartment, such as endosomes or lysosomes, can directly or indirectly inhibit lipid processing. Consequently, the accumulation caused toxic phospholipidosis, which is known to be associated with inhibition of coronavirus replication. It was suggested to be the inhibition mechanism of several candidates found in the cell-based assays. In this case, the potential drugs do not directly interfere with viral proteins, and likely cause cell-toxic side effects [84]. In fact, Hydroxychloroquine or similar drugs, which caused phospholipidosis in in vitro assays, eventually failed clinical trials. Also it should be noted that phospholipidosis induced by drugs could be overlooked in the period of clinical trials because it is typically an in vivo side effect that appears after chronic administration.
Although not all the drugs found in cell-based assays would have the phospholipidosis effect, mechanisms underlying the antiviral activities of the drugs remain largely elusive. In order to verify the underlying mechanisms of drugs, a multidisciplinary approach combining in silico, in vitro, and in vivo assays is promising.
Owing to the global efforts of experimental determination of SARS-CoV-2 protein structures, nearly 70% of the viral protein structurome has been covered within one and half years. Notably, the recent startling success of the machine learning-based structural prediction significantly contributed to the progress of the structurome. Consequently, an enormous number of in silico virtual screenings based on the viral protein structures or in vitro biochemical screening studies have been reported for repurposing approved drugs in a very short period. However, only a few of the detected drug candidates have successfully passed clinical trials and been approved at this point of time, highlighting a large gap lying between in silico/in vitro and in vivo/clinical studies for drug repurposing. Remdesivir is the sole approved COVID-19 therapeutic drug directory-targeted to the SARS-CoV-2 proteins. This fact might represent a potential limitation in repurposing the existing drugs, which were not particularly designed for targeting the SARS-CoV-2 proteins. Accordingly, the global efforts would be shifted to a specific design of new chemical entities against COVID-19. In fact, a couple of 3CLpro inhibitors, which have been newly designed by pharmaceutical companies including one in Japan, are being evaluated in phase 1 clinical trials. In such cases, the accumulation of the viral protein structures in complex with known drugs would contribute as primers to a new drug design not only in developing more effective therapeutic against COVID-19, but also preparing weapons against emerging new viruses in the near future.
This work was partly supported by the Platform Project for Supporting Drug Discovery and Life Science Research (Basis of Supporting Innovative Drug Discovery and Life Science Research (BINDS)) from AMED (JP20am0101111 and JP20am0101069), and Grants-in-Aid for scientific research from the Ministry of Education, Culture, Sports, Science and Technology-Japan (JP21H03547 and JP19K12211).
The authors declare no competing financial interests.
A.H., C.S., M.S., S.N., M.O., S.K., T.S. wrote the manuscript.