Next Article in Journal
SARS-CoV-2/COVID-19: Viral Genomics, Epidemiology, Vaccines, and Therapeutic Interventions
Next Article in Special Issue
Broad-Spectrum Host-Based Antivirals Targeting the Interferon and Lipogenesis Pathways as Potential Treatment Options for the Pandemic Coronavirus Disease 2019 (COVID-19)
Previous Article in Journal
Glucose-Regulated Protein 78 Interacts with Zika Virus Envelope Protein and Contributes to a Productive Infection
Previous Article in Special Issue
Repurposing Antiviral Protease Inhibitors Using Extracellular Vesicles for Potential Therapy of COVID-19
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Potential Drugs Targeting Early Innate Immune Evasion of SARS-Coronavirus 2 via 2’-O-Methylation of Viral RNA

by
José Antonio Encinar
1,* and
Javier A. Menendez
2,3,*
1
Institute of Research, Development and Innovation in Biotechnology of Elche (IDiBE) and Molecular and Cell Biology Institute (IBMC), Miguel Hernández University (UMH), 03202 Alicante, Spain
2
Program Against Cancer Therapeutic Resistance (ProCURE), Metabolism and Cancer Group, Catalan Institute of Oncology, 17005 Girona, Spain
3
Girona Biomedical Research Institute, 17007 Girona, Spain
*
Authors to whom correspondence should be addressed.
Viruses 2020, 12(5), 525; https://doi.org/10.3390/v12050525
Submission received: 22 April 2020 / Revised: 5 May 2020 / Accepted: 8 May 2020 / Published: 10 May 2020
(This article belongs to the Special Issue Drug-Repositioning Opportunities for Antiviral Therapy)

Abstract

:
The severe acute respiratory syndrome coronavirus-2 (SARS-CoV-2) causing the COVID-19 respiratory disease pandemic utilizes unique 2′-O-methyltransferase (2′-O-MTase) capping machinery to camouflage its RNA from innate immune recognition. The nsp16 catalytic subunit of the 2′-O-MTase is unusual in its requirement for a stimulatory subunit (nsp10) to catalyze the ribose 2′-O-methylation of the viral RNA cap. Here we provide a computational basis for drug repositioning or de novo drug development based on three differential traits of the intermolecular interactions of the SARS-CoV-2-specific nsp16/nsp10 heterodimer, namely: (1) the S-adenosyl-l-methionine-binding pocket of nsp16, (2) the unique “activating surface” between nsp16 and nsp10, and (3) the RNA-binding groove of nsp16. We employed ≈9000 U.S. Food and Drug Administration (FDA)-approved investigational and experimental drugs from the DrugBank repository for docking virtual screening. After molecular dynamics calculations of the stability of the binding modes of high-scoring nsp16/nsp10–drug complexes, we considered their pharmacological overlapping with functional modules of the virus–host interactome that is relevant to the viral lifecycle, and to the clinical features of COVID-19. Some of the predicted drugs (e.g., tegobuvir, sonidegib, siramesine, antrafenine, bemcentinib, itacitinib, or phthalocyanine) might be suitable for repurposing to pharmacologically reactivate innate immune restriction and antagonism of SARS-CoV-2 RNAs lacking 2′-O-methylation.

Graphical Abstract

1. Introduction

As of 4 May 2020, the pandemic of coronavirus disease 2019 (COVID-19) respiratory disease caused by the pathogenic severe acute respiratory syndrome coronavirus-2 (SARS-CoV-2) has led to more than 3,500,000 confirmed cases and more than 250,000 deaths worldwide [1,2,3,4,5]. Laboratory-based studies using the nucleotide analog remdesivir—a pan-inhibitor of viral RNA-dependent RNA polymerases—and preliminary clinical reports with (hydroxy)chloroquine—an approved, anti-inflammatory drug used to treat malaria, lupus, and rheumatoid arthritis—suggest their potential benefit against SARS-CoV-2 infection and the possible amelioration of viral shedding [6,7,8,9,10,11,12,13]. Accordingly, clinical trials evaluating the survival or time to clinical improvement in severely ill adult patients hospitalized for COVID-19 after adding remdesivir or hydroxychloroquine to standard supportive care, and clinical trials exploring hydroxychloroquine for preventing secondary SARS-CoV-2 transmission following initial contact exposure, are either recruiting or underway. However, no antiviral drugs are yet available with proven efficacy for SARS-CoV-2 treatment or prophylactic strategies to successfully protect individuals at high risk for COVID-19 infection (e.g., close contacts, households, and healthcare workers).
The current development of novel therapeutics to counteract SARS-CoV-2 infection can be categorized into at least four different strategies, namely: (a) broad-spectrum anti-virals (e.g., remdesivir, ribavirin, cyclophilin, and interferon) [14,15]; (b) drugs targeting the proinflammatory hypercytokinemia (termed “cytokine storm”) driving the transition from first COVID-19 symptoms to acute respiratory distress syndrome (e.g., IL-6 antibody blockers, IL-1 receptor antagonists, and JAK inhibitors) [16,17,18,19,20]; (c) inhibitors of host cell proteases that participate in the priming of the viral Spike (S) glycoprotein [21,22,23,24]; and (d) therapeutics targeting the host–virus interface linking the viral S protein to the angiotensin-converting enzyme 2 (ACE2) receptor in host cells [25,26,27,28,29,30,31,32,33]. In the current pandemic, identifying new targets for already approved drugs (drug repurposing) might shorten the development time and reduce the cost compared with de novo discovery of new compounds targeting one or several of the repertoire of viral proteins (up to 29) [34,35]. The bulk of the drug repurposing efforts seem to be directed toward pharmacologically targeting 3CLpro/nsp5-dependent viral replication [36,37], RdRp/nsp12-driven viral RNA synthesis, and S protein-driven viral cellular entry [22].
SARS-CoV-2 RNAs are capped at the 5′ end to impede degradation by 5′ exoribonucleases, ensure efficient translation, and evade recognition by the host cell innate immune system [38,39,40,41,42]. Interestingly, the SARS-CoV-2 2′-O-methyltransferase (2′-O-MTase) nsp16 protein is an RNA cap-modifying enzyme that is devoid of enzymatic activity and is activated by nsp10, which interacts with nsp16 and selectively confers upon it 2′-O-MTase activity on N7-methyl guanine RNA caps [43,44,45,46,47]. Thus, the methylation process follows an ordered sequence whereby RNA cap guanine-N7-methyltransferase (N7-MTase, nsp14)-mediated N7-guanine methylation precedes nsp16/nsp10-catalyzed RNA 2′-O-methylation [45,46,47,48,49,50] (Figure 1A). Nsp10 binds nsp16 through a ≈930 Å activation surface in nsp10, a molecular event that promotes nsp16 binding to the capped RNA substrate and the methyl donor S-adenosyl-l-methionine (SAM), stabilizing the SAM-binding pocket and extending the capped RNA-binding groove [45,46,47]. The requirement of nsp10 for nsp16 to execute its 2′-O-MTase activity is a unique feature of SARS-CoV-2 that has not been found in any other virus or host cell. The recently described crystal structure of the nsp16/nsp10 heterodimer has revealed that the nsp16/nsp10 interface and the RNA substrate binding sites may represent better drug targets than the MTase active site for developing highly specific anti-SARS-CoV-2 drugs [45,46]. Crucially, the absence of 2′-O-MTase activity results in a significant attenuation of SARS-CoV infection, which is characterized by decreased viral replication and limited breathing difficulties in animal models [51]. Therefore, pharmacological exploitation of 2′-O-MTase activity might open new treatment and prevention avenues to restore viral RNA recognition and activate intrinsic cell immunity against SARS-CoV-2 (Figure 1A). However, despite the evident therapeutic implications, there have been almost no studies aimed at designing and developing highly specific anti-SARS-CoV-2 drugs impeding the functioning of the RNA cap 2′-O-MTase nsp16/nsp10 protein complex.
Here, we used a virtual screening approach of molecular docking of ≈9000 U.S. Food and Drug Administration (FDA)-approved investigational and experimental (discovery-phase) drugs to identify potential candidates that can be directed to the RNA cap 2′-O-MTase nsp16/nsp10 complex (Figure 1B). To provide a better prediction of the binding modes of the selected drugs to the nsp16/nsp10 protein complex and to explore the stability of the nsp16/nsp10–drug complexes, we performed 100 ns molecular dynamics (MD) simulations and monitored the convergence of the root mean square deviation (RMSD) plots of the nsp16/nsp10 protein complex backbone and the non-hydrogen atoms of the drugs. To further explore the stability of the binding modes of the predicted candidates to the nsp16/nsp10 protein complex, MD simulations were employed for free energy calculations using molecular mechanics Poisson–Boltzmann surface area (MM/PBSA) approximation approaches. Finally, we searched the literature to consider the pharmacological overlapping of the selected, high-scoring nsp16/nsp10-targeting drugs with key functional modules of the virus–host interactome that are relevant for the viral lifecycle, as well as with clinical and laboratory features of COVID-19.

2. Materials and Methods

2.1. Drug Structures

The 3D structures of all the tested compounds belonging to DrugBank Release v.5.1.5. were downloaded from the DrugBank website as SDS files [52] and converted to MOL2 format using Marvin Suite 6.0 tools. Briefly, the script executed from a console (Windows or Linux) was as follows: /ChemAxon/MarvinBeans/bin/molconvert mol 2 /path-SDF-file/3D-structures.sdf -o/path-mol2-files/DBnumber.mol2 -m -g -F. The MOL2 files were then converted to PDBQT files using Open Babel 2.4.1 suite by executing the following script from a console (Windows or Linux): /OpenBabel-2.4.1./obabel/path-mol2-files/DBnumber.mol2 -O /path-pdbqt-files/DBnumber.pdbqt –gen3d -xh -p 7.4.

2.2. Viral Proteins Structures

Crystallographic coordinates of the nsp16/nsp10 complex 6W4H with a suitable resolution for docking studies (1.8 Å) were downloaded from the Protein Data Bank (PDB) in the pdb format. The specific edition of protein structures was made using PyMol 2.0 software (PyMOL Molecular Graphics System, v2.3.3 Schrödinger, New York, NY, USA, LLC, at http://www.pymol.org/) without further optimization. The appropriate pH 7.4 protonation state of the nsp17/nsp10 side chains was created using the YASARA structure v19.12.14 software [53].

2.3. Molecular Docking Simulation

To perform the docking studies with AutoDockVina (v1.1.2, San Diego, CA, USA), the structures of nsp16 or the nsp16/nsp10 complex were transformed to the PDBQT format, including the atomic charges and atom-type definitions [54,55]. These preparations were performed using the AutoDock/Vina plugin with scripts from the AutoDock Tools package [56].
Molecular dockings simulations were carried out in three different regions of the nsp16 protein, namely: the SAM-binding site, the nsp16/nsp10 interface, and the RNA-binding groove. For each region, the grid dimensions were as follows: 23 × 23 × 23 (SAM-binding site), 18 × 25 × 15 (nsp16/nsp10 interface), and 30 × 23 × 30 (RNA groove). AutoDock/Vina was set up on a Linux cluster (see Acknowledgments section), which generates a file for each drug containing the molecular coordinates of up to 20 poses, as well the Gibbs free energy variation (ΔG, kcal/mol) for each pose [57].

2.4. Molecular Dynamics Simulations

YASARA dynamics v19.9.17 (Vienna, Austria) was employed to carry out all the MD simulations with AMBER14 as a force field. The simulation cell was allowed to include 20 Å surrounding the protein that was filled with water at a density of 0.997 g/mL. Initial energy minimization was carried out under relaxed constraints using steepest descent minimization. Simulations were performed in water under constant pressure and constant temperature (25 °C) conditions. To mimic a physiological environment, counter ions were added to neutralize the system (Na+ or Cl were added as a replacement for water to give a final NaCl concentration of 0.9% and the pH was maintained at 7.4). Hydrogen atoms were added to the protein structure at the appropriate ionizable groups according to both the calculated pKa and the simulation pH (i.e., a hydrogen atom was added if the computed pKa was higher than the pH). The pKa was computed for each residue according to the Ewald method [58]. All simulation steps were run using a pre-installed macro (md_run.mcr) within the YASARA suite. Data were collected every 100 ps. The molecular mechanics Poisson–Boltzmann surface area (MM/PBSA) calculations were used to determine the alchemical binding free energy of drug candidates against nsp16 using the YASARA macro md_analyzebindenergy.mcr, as previously described [59,60,61].

2.5. Evaluation of nsp16-Drug Interactions

The best poses of the predicted drug candidates with the best affinity scores were merged with the optimized protein structure in the pdb format using PyMol 2.0 software (New York, NY, USA). Contact residues of the merged structures were analyzed using the protein–ligand interaction profiler (PLIP) algorithm.

3. Results

We performed structure-based virtual screening involving molecular docking and MD strategies to select putative candidates among 8696 drugs in the DrugBank database. The target employed was 6W4H, the 1.80 Å resolution crystal structure of the SARS-CoV-2 nsp16/nsp10 complex that was recently released (2020-03-18) by the Center for Structural Genomics of Infectious Diseases (www.csgid.org). The database was docked into three different binding sites of the nsp16/nsp10 complex, namely: the SAM-binding site, the nsp16/nsp10 interface, and the RNA-binding groove (Figure 2). Once the conformational search algorithm probed the energy landscape of each drug in each binding site [54,55] (http://shaker.umh.es/Publications/PDF/nsp16/DrugBank-nsp16-docking.zip), we selected the top 20 high-scoring compounds (which therefore have potentially higher affinity) for each site in a first filtering step. The selected compounds exhibited binding energy values below −10 kcal/mol in the SAM-binding site, below −8 kcal/mol in the nsp16/nsp10 interface, and below −9 kcal/mol in the RNA binding groove (Table S1), thus representing KD (KD = expΔG/RT) values in the nanomolar to low-micromolar range. The compounds were then grouped based on interactions (i.e., predicted bioactivity) with the same (or different) binding sites (Figure 1B).
Based on the results of the molecular docking calculations, twelve drugs were predicted to occupy the SAM binding site with a high affinity (Table S2), twelve were predicted to occupy the nsp16/nsp10 interface with a high affinity (Table S3), fourteen were predicted to occupy the RNA-binding groove with a high affinity (Table S4), three were predicted to occupy both the nsp16/nsp10 interface and the SAM binding site with a high affinity (Table S5), one was predicted to dually target the SAM-binding site and the RNA-binding groove (Table S6), one was predicted to dually target the interface and the RNA-binding groove (Table S7), and four were predicted to occupy all the tested cavities of the nsp16/nsp10 protein complex with a high affinity (Table S8).
Because the flexibility of the target-binding site during the molecular recognition process is an essential but frequently overlooked aspect to be considered in molecular docking, we approached such flexibility issues using MD techniques. Accordingly, we performed MD simulations for 100 ns for each of the high-scoring nsp16/nsp10–drug complexes to confirm the kinetic stability of the poses and validate the binding poses obtained by docking. The nsp16/nsp10 protein backbone RMSD plots of the drugs’ heavy atoms, measured after superimposing nsp16/nsp16 on its reference structure during the MD simulation, were prepared in parallel. A second filter stage considered the stable behaviors of the drug–protein complexes when RMSD values did not exceed 10–15 Å and the starting and ending poses of the drug were similar throughout the MD simulation. Those drugs that did not meet such requirements were discarded as potential candidates (Table S9). To further verify the stability of the nsp16/nsp10–drug complexes, we evaluated the MM/PBSA parameters, which estimate the free energy of the binding of small ligands to biological macromolecules, using MD simulations of the receptor–ligand complex; the results showed good correlations with the experimentally obtained values despite excluding the conformational entropy or the number and free energy of water molecules in the binding site [62,63]. In this third filter, a well-recognized scoring function that commonly serves as a powerful tool in the correct ranking of drug candidates, which is the fact that more positive energies indicate the occurrence of stronger binders. We finally combined the second and third filtering criteria to generate a definitive list of candidates that were classified as either “strong” (>25 kcal/mol) or “weak” (<25 kcal/mol) (Table S10).
Twenty different compounds were cataloged as strong candidates, with three of them were predicted to stably target the SAM-binding site, fourteen of them were predicted to stably target the nsp16/nsp10 interface, and four of them were predicted to stably target the RNA binding groove (Figure 2).
To facilitate the understanding of the predicted binding modes of the candidates selected in Table S2, we have illustrated the detailed interactions of nsp16 (or the nsp16/nsp10 complex) with the SAM-binding-site-targeting drugs antrafenine, entrectinib, and lifitegrast (Figure 3), with the nsp16/nsp10 interface-targeting drugs TG-100801, LY-2456302, tegobuvir, sonidegib, siramesine, and lifirafenib (Figure 4), and with the RNA binding groove-targeting drugs DB02449, DB03067, and PF-04457845 (Figure 5). Such depictions include the MS simulations at 0 and 100 ns, the time evolution of RMSD relative to the initial structure of nsp16/nsp10 in the absence and presence of drug candidates, MM/PBSA binding energy analyses calculated from the trajectory of the 100 ns (or last 30 ns) MD simulation, and identification of amino acid residues participating in the drug-binding pocket.
As internal controls, RMSD-based evaluation of the protein’s thermodynamic stability during the MD simulation included not only SAM but also the nsp16-SAM and/or nsp16/nsp10–SAM complex; given that SAM is a methyl donor for the methylation reaction that operates as an nsp16 coenzyme, it is not expected to stabilize the protein structure. The calculated MM/PBSA values of the nsp16–SAM complex in the absence of nsp10 notably increased in its presence, thereby in silico confirming that the partnering of nsp10 stabilizes the interaction of SAM with the catalytic nsp16 [45,46]. Drug candidates such as antrafenine, entrectinib, and lifitegrast are predicted to compete with SAM for the SAM-binding pocket of nsp16 (Figure 3), thereby impeding its methyltransferase activity. Drug candidates such as TG-100801, LY-2456302, tegobuvir, sonidegib, siramesine, and lifirafenib are predicted to bind tightly (mostly via hydrophobic interactions) to the activating interface of nsp16 (Figure 4), thereby impeding the obligatory formation of the active nsp16/nsp10 complex. Drug candidates such as DB02449, DB03067, and PF-04457845 are predicted to stably occupy the RNA groove, mostly via hydrophobic interactions and hydrogen bonds (Figure 5), thereby preventing a proper positioning of the viral RNA that enables its capping.

4. Discussion

The SARS-CoV-2 causing the current pandemic of COVID-19 respiratory disease tends to have a longer incubation period (2–11 days) than the influenza virus (1–4 days). This long prodromal phase, which significantly increases the pre- and oligosymptomatic virus transmission capability of infected individuals, is largely due to its ability to efficiently evade immune detection and to dampen human immune responses. This adaptation of SARS-CoV-2 to escape recognition and activation of the host immune response involves the RNA cap modifying enzyme 2′-O-MTase, which provides the viral mRNA the ability to camouflage itself from the host cell [64]. Accordingly, human coronavirus mutants lacking 2′-O-MTase are exquisitely sensitive to a higher expression of induced type-I interferons [40,41,42]. However, despite the strong rationale for the development of anti-SARS-CoV-2 interventions, such as live attenuated vaccines (e.g., employing 2′-O-methylation RNA mutants) and small-molecule inhibitors, very few attempts have been made to reactivate the innate immune restriction and antagonism of SARS-CoV-2 RNAs lacking 2′-O-methylation via the pharmacological blockade of the coronaviral 2-O′-MTase.
Here, we provide a first-in-class computational basis for drug repositioning or de novo drug development based on three differential traits of the intermolecular interactions of the RNA cap 2′-O-MTase nsp16/nsp10 protein complex, namely: the (nsp10-stabilized) SAM-binding pocket of nsp16, the (nsp10-extended) RNA-binding groove of nsp16, and the unique nsp16/nsp10 interaction interface that is obligatorily required by nsp16 to execute its 2′-O-MTase activity [44,45,46,47]. We provide a careful description and a ranked list of compounds that can be tested experimentally. Some of the identified small molecules might be repurposed or be considered as potential lead molecules for further optimization and drug development against SARS-CoV-2. Of the initially selected 47 docking candidates, only 20 were considered further as strong candidates upon review of the MD/RMSD stability criteria (Figure 1B). We acknowledge that such conventional MD simulations of nanosecond time scales, which will require enhanced sampling methods to provide a complete picture of the key residues used for drug candidates for effective nsp16/nsp10 inhibition, should be viewed as a preliminary approach to identify moieties critical for the retention of key nsp16/nsp10–drug interactions. Nevertheless, they may be useful when analyzing the predicted candidates as lead molecules for optimization and drug development against the nsp16/nsp10 complex. We also acknowledge that some of the predicted drugs, especially those for the treatment of tumors, are expected to produce serious side-effects, and in some cases, our data inferred that high concentrations are required to inhibit nsp16/nsp10 effectively. Accordingly, these drugs may not be considered safe therapies for a population of infected/sick people. Furthermore, while computational approaches offer novel testable hypotheses for systematic drug repositioning, it can be argued that single virus proteins will carry a high risk of drug resistance given the rapid evolution of virus genomes.
Because viruses require host cellular factors for successful replication, the systematic identification of the so-called virus–host protein–protein interactome [34,65] would be expected to provide key insights into effective molecular targets for developing broadly acting effective treatments against SARS-CoV-2 and other upcoming coronavirus strains. A recent landmark systematic mapping of the interaction landscape between SARS-CoV-2 proteins and human proteins provided the interaction map of nsp10 [66]. Unfortunately, the authors were unable to generate data for the nsp16 protein. Despite this, and with the intention of providing a broader mechanistic context to allow investigators to experimentally or clinically pursue drugs directly, we decided to carefully re-evaluate the selected anti-nsp10/nsp16 candidates in terms of not only their pharmacological overlapping with key functional modules of the SARS-CoV-2-human interactome relevant for the viral lifecycle but also with clinical and laboratory features of COVID-19 (Figure 6). Based on such a comprehensive evaluation, we prioritized several candidates, as described below.

4.1. Virus Cell Entry and Replication

4.1.1. Tegobuvir

Tegobuvir (GS-9190), a novel inhibitor of the hepatitis C virus (HCV) RNA replication in vitro, has shown potential antiviral activity in patients chronically infected with genotype 1 HCV [67,68,69,70]. A recent docking-based repurposing approach predicted the ability of tegobuvir to target the SARS-CoV-2 core RNA-synthesis machinery by occupying the interface between the nsp12 polymerase and its essential co-factor nsp7 [71]. Our present in silico approach predicted the ability of tegobuvir to occupy the activating interface of the nsp16 2′-O-MTase and its essential co-factor nsp10. It is noteworthy that ribavirin, a prodrug that is metabolized to nucleoside analogs to block viral RNA synthesis and viral mRNA capping proposed as an anti-viral therapy against SARS-CoV-2 based on in vitro data [72,73], was previously demonstrated to exhibit limited efficacy against coronaviruses in vivo [74]. This is because the coronaviral nsp14 protein, a bifunctional enzyme carrying RNA cap guanine N7-MTase and 3′-5′ exoribonuclease activities, associates with the low-fidelity nsp12 viral RNA polymerase, thereby generating a robust RNA synthesis and proofreading pathway in which the antiviral ribavirin 5′-monophosphate is significantly incorporated but also readily excised from RNA. Tegobuvir might be worthy of further investigation as a non-nucleoside inhibitor of RNA synthesis and 2′-O-methylation of viral mRNA to attenuate SARS-CoV-2 replication.

4.1.2. Siramesine

The gold standard sigma-2 receptor agonist siramesine (formerly known as Lu-28-179) was among the high-scoring candidates predicted to stably occupy the activating surface occurring between nsp10 and nsp16. Siramesine is a direct lysosomal toxin that disrupts the lysosomal pH gradient, which is likely done by targeting vacuolar V-ATPase in the endosomal/lysosome membrane, thereby triggering lysosomal dysfunction/destabilization and the accumulation of autophagosomes to counteract the lysosomal cell death pathway [75,76,77,78,79]. Coronaviruses, such as SARS-CoV-2, interact differentially with the autophagy machinery to utilize components for virus replication while attenuating autophagy per se to promote viral production [80]. Compounds such as siramesine, which can reinstate cytoprotective autophagy, might reduce viral production. Moreover, siramesine-related drugs, such as the sigma-1 receptor agonist chloroquine [81,82] and the endo/lysosomal V-ATPase inhibitor bafilomycin A1, may also promote viral trapping with perinuclear vacuoles and interference with terminal glycosylation of the SARS-CoV-2 cellular receptor ACE2 [83,84]. Siramesine belongs to the group of lysosomotropic agents with potential therapeutic efficacy against COVID-19 by impacting cellular trafficking through the endosome/lysosome endocytic pathways and by promoting the accumulation of autophagosomes, which associate with a rapid and sustained inhibition of mTOR, a negative regulator of autophagy [76,77].
It should be noted that other SARS-CoV-2 non-structural proteins, such as the virus replication component nsp6, have been shown to interact with sigma receptors [66], which are also targeted (albeit with lower affinity) by therapeutic candidates for COVID-19 treatment and the prevention of secondary SARS-CoV-2 transmission, such as (hydroxy)chloroquine. Future work should evaluate the apparently central role of the sigma-receptor-regulated endoplasmic reticulum stress/unfolded protein response pathway [85,86,87], not only as a major aspect of the coronavirus-host interactome but also in the mechanistic context of sigma-receptor-targeting drugs with therapeutic potential against SARS-CoV-2.
Siramesine was originally designed to treat anxiety and depression [88,89]. Although clinical trials failed to show satisfactory efficacy in the treatment of psychiatric disorders, the drug was well-tolerated and non-toxic in humans. The lack of side-effects; the strong and immediate effects of siramesine on lysosome acidification, trafficking, and autophagy activation; and its predicted capacity to target the nsp16/nsp10 interface might open new avenues for the therapeutic consideration of this lysosomotropic agent and/or closely related sigma-2-receptor-targeting drugs [66,90,91,92].

4.2. Host Immune Response

4.2.1. Bemcentinib

Bemcentinib (BGB324/R-428), a first-in-class highly selective oral AXL tyrosine kinase inhibitor in phase II development, was also predicted to stably occupy the nsp16/nsp10 activating interface. Although mainly recognized for its ability to block epithelial–mesenchymal transition phenomena in various cancer types [93,94,95], bemcentinib is now known for its anti-inflammatory and anti-fibrotic roles in non-cancer diseases [96,97], which are two phenotypic aspects that might be relevant in the therapeutic management of COVID-19. Moreover, AXL is a family member of the TAM receptor tyrosine kinases that operates both as a pleiotropic inhibitor of the innate immune response and a docking site during the cellular entry of some viruses [98,99,100]. As AXL receptor signaling participates in the downmodulation of interferon-related host immune responses to promote viral infection, inhibiting AXL function might be a new therapeutic avenue toward reducing the transmission of some of the most troublesome emerging enveloped viruses, including dengue, West Nile, Ebola, lentivirus, and most recently, Zika [99,100]. If experimentally confirmed as a blocker of the nsp16/nsp10-driven viral innate immunoevasion via suppression of viral RNA 2′-O-methylation, bemcentinib might also stimulate an advantageous increase in the early expression of interferon-stimulated genes that are capable of promoting a greater attenuation of SARS-CoV-2 replication. Bemcentinib was recently shown to block lysosomal acidification and promote the accumulation of autophagosomes in cancer cells independently of its inhibitory effects on AXL [101,102], thereby suggesting an additional overlapping mechanism with the above-mentioned lysosomotropic agents.

4.2.2. Sonidegib

Sonidegib (Odomzo®), an FDA-approved oral, selective smoothened (SMO) inhibitor used to treat basal cell carcinomas via regulation of the hedgehog (Hh) signaling pathway [103,104,105] was among the best drug candidates predicted to impede the formation of the obligatory nsp16/nsp10 complex by stably occupying the nsp16 activating surface. The Hh signaling is targeted by numerous viral pathogens to control the local infected environment [106]. Due to the key roles of Hh signaling in tissue repair processes and immunity (e.g., immune cell proliferation and activation), the ability of viral pathogens (e.g., influenza, Epstein-Barr virus, Hepatitis B and C virus, human immunodeficiency virus) that cause damage to the host tissue can exacerbate such Hh-driven responses (e.g., activation of STAT3 and upregulation of IL-6 expression) to drive a detrimental outcome, including fibrosis and an imbalanced immune response [106]. Accordingly, the therapeutic use of molecules inhibiting Hh signaling, such as sonidegib, might be expanded for use as potent, broad-spectrum inhibitors of pathogenic infections by targeting highly-conserved host targets also shared by SARS-CoV-2/COVID-19.

4.2.3. Itacitinib

Itacitinib is a JAK1-selective inhibitor of the JAK/STAT pathway with >20-fold selectivity over JAK2 and a >100-fold selectivity over JAK3 and TYK2 [107,108]. Cell and animal models and phase II studies have shown potent, pre- and clinical efficacy of itacitinib against psoriasis and rheumatoid arthritis in a JAK2-independent manner [109,110,111,112,113]. As with other JAK inhibitors, such as baricitinib or tofacitinib, itacitinib is being studied for the prevention of the so-called cytokine release syndrome (CRS) or cytokine storm syndrome (CSS) [114], an acute systemic inflammatory syndrome characterized by fever and multiple organ dysfunction [115]. It is now known that a subgroup of patients with severe COVID-19 have had a CSS that occurred 7–10 days after the disease onset, coinciding with the peak of respiratory distress [14,15,116,117]; immune dysregulation, rather than the level of peak viremia, is responsible for the SARS-CoV-2 infection driving a too-little-and-too-late type-I-interferon response accompanied by aberrant proinflammatory cytokine secretion by alveolar macrophages and subsequent CD4+ and CD8+ T-cell dysfunction. In this life-threatening clinical scenario of hyperinflammation with dysregulated and profuse immune responses causing lung damage, the use of antivirals as single-agents may be insufficient to stop the cytokine storm, pulmonary destruction, and respiratory distress in patients. Accordingly, pharmacological suppression of hypercytokinemia using a selective cytokine blockade and JAK inhibition could improve mortality. Targeted immunomodulation with the dual JAK1/JAK2 inhibitor baricitinib is currently being evaluated in clinical trials aimed at reducing the SARS-CoV-2-driven CCS, ameliorating pulmonary inflammation, and improving mortality [118,119]. Because different JAK inhibitors modulate distinct cytokine pathways to varying degrees but fail to potently or continuously inhibit an individual cytokine pathway [111,113], further clinical studies are warranted to explore whether nsp16/nsp10 inhibition using a JAK inhibitor, such as itacitinib, might provide a means to first re-activate an innate immune response (via nsp16/nsp10 targeting) and later provide prophylactic prevention against the CSS.

4.3. Neurological Effects

4.3.1. MK-3207, PF-04457845, and Antrafenine

The small-molecule inhibitor of the CGRP receptor belonging to the “gepant” class of compounds for acute migraine therapy MK-3207 [120,121,122,123,124,125] was among the high-scoring candidates predicted to target the nsp16/nsp10 2′-O-MTase complex. A significant number of patients with COVID-19 show neurologic manifestations, including not only stroke and encephalopathy but also headache, seizures, nausea, olfactory disorders, and vomiting [126,127,128,129], which are also common symptoms of acute migraines. Recent reports of profound myocarditis and fatal arrhythmias have suggested the potential influence of SARS-CoV-2 on the cardiovascular system [130,131,132]. The neuroinvasive potential of SARS-CoV-2, which has been shown to invade the central nervous system and to infect the brainstem [133,134], may play an unexpected role in the respiratory distress and cardiovascular impact of COVID-19. Coronaviruses, such as SARS-CoV-2, can travel retrograde along neuronal synapses from the mechanoreceptors and chemoreceptors present in the lung and lower respiratory airways to the central cardiorespiratory center in the brain [126]. It is thus noteworthy that neuropeptides, such as CGRP and their G-protein coupled receptors (GPCRs), are present and play a prominent role not only in the lung but also in the brain area involved in ventilatory and cardiovascular regulation [135,136,137,138].
Although few studies have evaluated the integrative physiological actions of CGRP, it might be of interest to evaluate whether CGRP receptor-targeted gepants might provide not only symptomatic alleviation of COVID-19-related neurological symptoms but also beneficial central actions on ventilatory and cardiovascular disturbances accompanying COVID-19. Moreover, the recent realization that acute lung injury induces the release of CGRP, which in turn is implicated in COVID-19 symptoms, such as cough, fever, pain, and the release of IL-6, might illuminate a new scenario in which “gepants” might be used to target the SARS-CoV-2 nsp16/nsp10 complex while treating excessive lung inflammation caused by SARS-CoV-2. Because the neuronal system is responsible for pain transmission, inflammation, and immunomodulation throughout the entire pulmonary system, including the localized release of COVID-19-related cytokines [139], anti-inflammatory drugs previously employed for reducing inflammatory and neuroinflammatory pain, such as the fatty acid amide hydrolase inhibitor PF-04457845 [140,141,142] and the cyclooxygenase inhibitor antrafenine [143], might be repurposed to dually target the SARS-CoV-2 nsp16/nsp10 complex while reducing the pulmonary innervation system in ARDS associated with advanced COVID-19 disease.

4.3.2. TG-100801 and Lifitegrast

The best predicted candidate in terms of binding energy and stability targeting the unique nsp16/nsp10 interaction interface was TG-100801, the first topically applied, multitargeted VEGFR/Src kinase inhibitor to advance into the clinic. TG-100801 is administered noninvasively as an eye drop and has been designed to suppress VEGF-mediated leakage and additional kinase targets associated with inflammation, edema, and angiogenesis in several pathological eye conditions, including viral infections [144,145]. Lifitegrast, an FDA-approved drug for the treatment of dry eye syndrome [146,147], was also among the final list of strong nsp16/nsp10-targeting candidates. While the potential extra-respiratory transmissible routes of SARS-CoV-2 remain unclear, clinical entities, such as conjunctivitis, anterior uveitis, retinitis, and optic neuritis, have been documented in feline and murine coronavirus models. Indeed, SARS-CoV-2 might be transmitted via aerosol contact with the conjunctiva, causing conjunctival congestion as the initial symptom of COVID-19, and viral particles are present in ocular secretions [148,149,150]. Thus, because exposed mucous membranes and unprotected eyes appear to significantly increase the risk of SARS-CoV-2 transmission, the possibility of ocular SARS-CoV-2 infection cannot be ignored and the increase in awareness of eye protection is necessary. Indeed, dry eye, blurred vision, and foreign body sensation have been shown to rank as the top three COVID-19-related ocular symptoms and diseases [151]. The computational prediction of TG-100801 and lifitegrast as anti-SARS-CoV-2 drugs might have not only pharmacophoric implications for designing more specific compounds targeting the nsp16/nsp10 interface but could also provide a pharmacological approach of blocking eye-related transmission routes or treating SARS-CoV-2-driven ocular manifestations in patients.

4.4. Decontamination and Self-Disinfection

Phthalocyanine

Phthalocyanine is an archetypal member of so-called photodynamic antimicrobials, in which the combination of a sensitizing drug and visible light causes selective destruction of viruses, bacteria, and other pathogens when applied to consumer or medical products [152,153,154,155]. Although the length of time that SARS-CoV-2 can survive on inanimate surfaces [156,157] varies depending on numerous factors, the risk of transmission of COVID-19 in social contact environments and non-healthcare work settings can be largely minimized through good general hygiene, including handwashing and wearing a face mask. A different scenario occurs in the so-called “infectious clean” following the discharge or transfer of patients with COVID-19 in healthcare facilities, which requires both thorough cleaning and disinfection for environmental decontamination. Currently, such procedures involve physical cleaning with detergent followed by disinfection with a hospital-grade disinfectant with activity against some viruses, or a chlorine-based product, such as sodium hypochlorite. Photochemotherapy using metal-derivatives of phthalocyanines, which might become activated even under weak indoor light, has been investigated for the disinfection of lipid envelope viruses, including HIV, hepatitis B virus, herpes simplex virus 1, and influenza A virus [H1N1] [152,153,154,155]. Our current prediction of phthalocyanine as a putative pan-inhibitor of the SAM-binding site, the nsp16/nsp10 interface, and the RNA-binding groove of the SARS-CoV-2 2′-O-MTase, might provide an additional biological layer of interest to the usage of phthalocyanine derivatives as biocidal molecules that can be incorporated into self-disinfecting materials, such as fabrics, films, and coatings, which can play a vital role in preventing the transmission of SARS-CoV-2, especially in healthcare settings.
We acknowledge that other candidates predicted to target the nsp16/nps10 2′-O-MTase might provide additional clinico-molecular advantages against SARS-CoV-2/COVID-19. Lifirafenib, an orally available dual RAF kinase/EGFR inhibitor with an acceptable risk–benefit profile [158], and more particularly, entrectinib (Rozlytrek, RDX-1; Genentech Inc.: South San Francisco, CA, USA), a potent oral inhibitor of the tyrosine kinase TRK A/B/C, ROS1, and ALK with a good safety profile [159,160,161,162], were predicted to highly stably interact with the nsp16/nsp10 interface and the (nsp10-stabilized) SAM-binding pocket, respectively. At a moment in which the oncology community is carefully balancing treatment decisions for lung cancer patients in light of the COVID-19 pandemic (e.g., withholding immunotherapy from patients whose tumors do not have known driver mutations; interactions of COVID-19 treatments in clinical trials, such as remdesivir and hydroxychloroquine, with immunotherapeutics and tyrosine kinase inhibitors), if COVID-19-positive oncologic patients would be treated with tyrosine kinase inhibitors, such as entrectinib, associating clinical outcomes with treatment with these drugs may merit investigation.

5. Conclusions

The COVID-19-causing virus SARS-CoV-2 has evolved distinct 2′-O’-MTase machinery to autonomously modify the 5′ cap of its RNA. Such 2′-O-methylation is a crucial mechanism used to subvert innate host antiviral responses by escaping interferon-stimulated proteins that are capable of detecting non-2′-O-methylated viral RNA and inhibiting its translation [40]. Although the previous biochemical and structural delineation of distinctive activating traits of nsp16 by nsp10 immediately suggested a new way to design and develop highly SARS-CoV-2-specific antivirals [45,46], very few attempts have been made to reactivate the innate immune restriction and antagonism of SARS-CoV-2 RNAs by pharmacologically preventing their RNA cap 2′-O-methylation. Here, we aimed to systematically identify repurposable drugs that we predict can target the RNA cap 2′-O-MTase nsp16/nsp10 protein complex but also other molecular and clinical traits specific to SARS-CoV-2 and COVID-19. Although we acknowledge that our final computational prediction might appear too wide for pre-clinical and clinical validation, we expect that by providing a multi-layer network framework (Figure 6), our findings should help to select the narrowest list of candidates that can be rapidly tested before evaluating their in vivo efficiency and side-effects. Indeed, some of the predicted drugs (e.g., tegobuvir, sonidegib, siramesine, antrafenine, bemcentinib, itacitinib, and phthalocyanine) might be suitable for rapid repurposing to pharmacologically reactivate the innate immune restriction and antagonism of SARS-CoV-2 RNAs lacking 2′-O-methylation.

Supplementary Materials

The following are available online at https://www.mdpi.com/1999-4915/12/5/525/s1. Table S1: Binding energies and dissociation constants of top-scoring nsp16/nsp10-targeting candidates. Table S2: Drug candidates targeting the SAM-binding site of the nsp16/nsp10 protein complex. Table S3: Drug candidates targeting the nsp16/nsp10 interface of the nsp16/nsp10 protein complex. Table S4: Drug candidates targeting the RNA-binding groove of the nsp16/nsp10 protein complex. Table S5: Drug candidates targeting the SAM-binding site and nsp16/nsp10 interface of the nsp16/nsp10 protein complex. Table S6: Drug candidates targeting the SAM-binding site and RNA-binding groove of the nsp16/nsp10 protein complex. Table S7: Drug candidates targeting the nsp16/nsp10 interface and RNA-binding groove of the nsp16/nsp10 protein complex. Table S8: Drug candidates targeting the SAM-binding site, nsp16/nsp10 interface, and RNA-binding groove of the nsp16/nsp10 protein complex. Table S9: Filtering of nsp16/nsp10-targeting candidate drugs based on MD simulations for drug-protein complexes. Table S10: Selection of nsp16/nsp10-targeting candidate drugs based on MD simulations for drug-protein complexes.

Author Contributions

J.A.E. and J.A.M. conceived, initiated, and coordinated the project; J.A.E. performed the virtual screening; J.A.E. and J.A.M. analyzed the computational data; J.A.M. drafted the manuscript; J.A.E. and J.A.M. discussed the results. All authors have read and agreed to the published version of the manuscript.

Funding

Work in the Menendez laboratory is supported by the Spanish Ministry of Science and Innovation (grant SAF2016-80639-P, Plan Nacional de I+D+I, founded by the European Regional Development Fund, Spain) and by an unrestricted research grant from the Fundació Oncolliga Girona (Lliga catalana d’ajuda al malalt de càncer, Girona). The Spanish Ministry of Economy and Competitiveness (MINECO, Project RTI2018-096724-B-C21) and the Generalitat Valenciana (PROMETEO/2016/006) supports work in the Encinar laboratory.

Acknowledgments

We are grateful to the High-Performance Computing Center Stuttgart (HLRS) for letting us take advantage of the Hawk computer cluster (https://hlrs.de/systems/hpe-apollo-hawk/), which has provided ≈70% of the calculation time necessary to prepare this manuscript. We wish to especially thank Thomas Beisel for his kind assistance and cooperation. We are also grateful to the Cluster of Scientific Computing (http://ccc.umh.es/) of the Miguel Hernández University (UMH) for providing computing facilities. This manuscript was written during the second and third weeks of the lockdown in Spain due to the SARS-CoV-2 coronavirus pandemic in March–April 2020. This paper is dedicated to the memory of all those who have succumbed to COVID-19 during those hard days. The authors would like to thank Kenneth McCreath for editorial support.

Conflicts of Interest

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest. Ethics approval was not required for this study as per the local legislation. The authors have not filed for patent protection on the use of predicted drugs for targeting SARS-CoV-2 2′O-MTase to ensure all this information is freely available to accelerate the discovery of a treatment for COVID-19.

References

  1. Wu, F.; Zhao, S.; Yu, B.; Chen, Y.M.; Wang, W.; Song, Z.G.; Hu, Y.; Tao, Z.W.; Tian, J.H.; Pei, Y.Y.; et al. A new coronavirus associated with human respiratory disease in China. Nature 2020, 579, 265–269. [Google Scholar] [CrossRef] [Green Version]
  2. Li, Q.; Guan, X.; Wu, P.; Wang, X.; Zhou, L.; Tong, Y.; Ren, R.; Leung, K.S.M.; Lau, E.H.Y.; Wong, J.Y.; et al. Early Transmission Dynamics in Wuhan, China, of Novel Coronavirus-Infected Pneumonia. N. Engl. J. Med. 2020, 382, 1199–1207. [Google Scholar] [CrossRef] [PubMed]
  3. Wilder-Smith, A.; Chiew, C.J.; Lee, V.J. Can we contain the COVID-19 outbreak with the same measures as for SARS? Lancet Infect. Dis. 2020, 20, e102–e107. [Google Scholar] [CrossRef] [Green Version]
  4. Gates, B. Responding to Covid-19—A Once-in-a-Century Pandemic? N. Engl. J. Med. 2020. [Google Scholar] [CrossRef] [PubMed]
  5. Novel Coronavirus (2019-nCoV) Situation Reports. Available online: https://www.who.int/emergencies/diseases/novel-coronavirus-2019 (accessed on 4 May 2020).
  6. Sheahan, T.P.; Sims, A.C.; Leist, S.R.; Schäfer, A.; Won, J.; Brown, A.J.; Montgomery, S.A.; Hogg, A.; Babusis, D.; Clarke, M.O.; et al. Comparative therapeutic efficacy of remdesivir and combination lopinavir, ritonavir, and interferon beta against MERS-CoV. Nat. Commun. 2020, 11, 222. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  7. Wang, M.; Cao, R.; Zhang, L.; Yang, X.; Liu, J.; Xu, M.; Shi, Z.; Hu, Z.; Zhong, W.; Xiao, G. Remdesivir and chloroquine effectively inhibit the recently emerged novel coronavirus (2019-nCoV) in vitro. Cell Res. 2020, 30, 269–271. [Google Scholar] [CrossRef] [PubMed]
  8. Martinez, M.A. Compounds with therapeutic potential against novel respiratory 2019 coronavirus. Antimicrob. Agents Chemother. 2020, 64. [Google Scholar] [CrossRef] [Green Version]
  9. Liu, J.; Cao, R.; Xu, M.; Wang, X.; Zhang, H.; Hu, H.; Li, Y.; Hu, Z.; Zhong, W.; Wang, M. Hydroxychloroquine, a less toxic derivative of chloroquine, is effective in inhibiting SARS-CoV-2 infection in vitro. Cell Discov. 2020, 6, 16. [Google Scholar] [CrossRef] [Green Version]
  10. Gautret, P.; Lagier, J.C.; Parola, P.; Hoang, V.T.; Meddeb, L.; Mailhe, M.; Doudier, B.; Courjon, J.; Giordanengo, V.; Vieira, V.E.; et al. Hydroxychloroquine and azithromycin as a treatment of COVID-19: Results of an open-label non-randomized clinical trial. Int. J. Antimicrob. Agents. 2020, 105949. [Google Scholar] [CrossRef]
  11. Yao, X.; Ye, F.; Zhang, M.; Cui, C.; Huang, B.; Niu, P.; Liu, X.; Zhao, L.; Dong, E.; Song, C.; et al. In Vitro Antiviral Activity and Projection of Optimized Dosing Design of Hydroxychloroquine for the Treatment of Severe Acute Respiratory Syndrome Coronavirus 2 (SARS-CoV-2). Clin. Infect. Dis. 2020, ciaa237. [Google Scholar] [CrossRef] [Green Version]
  12. Touret, F.; de Lamballerie, X. Of chloroquine and COVID-19. Antivir. Res. 2020, 177, 104762. [Google Scholar] [CrossRef] [PubMed]
  13. Devaux, C.A.; Rolain, J.M.; Colson, P.; Raoult, D. New insights on the antiviral effects of chloroquine against coronavirus: What to expect for COVID-19? Int. J. Antimicrob. Agents. 2020, 105938. [Google Scholar] [CrossRef] [PubMed]
  14. Totura, A.L.; Bavari, S. Broad-spectrum coronavirus antiviral drug discovery. Expert Opin. Drug Discov. 2019, 14, 397–412. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  15. Sheahan, T.P.; Sims, A.C.; Graham, R.L.; Menachery, V.D.; Gralinski, L.E.; Case, J.B.; Leist, S.R.; Pyrc, K.; Feng, J.Y.; Trantcheva, I.; et al. Broad-spectrum antiviral GS-5734 inhibits both epidemic and zoonotic coronaviruses. Sci. Transl. Med. 2017, 9, eaal3653. [Google Scholar] [CrossRef] [Green Version]
  16. Zhou, F.; Yu, T.; Du, R.; Fan, G.; Liu, Y.; Liu, Z.; Xiang, J.; Wang, Y.; Song, B.; Gu, X.; et al. Clinical course and risk factors for mortality of adult inpatients with COVID-19 in Wuhan, China: A retrospective cohort study. Lancet 2020, 395, 1054–1062. [Google Scholar] [CrossRef]
  17. Mehta, P.; McAuley, D.F.; Brown, M.; Sanchez, E.; Tattersall, R.S.; Manson, J.J. HLH Across Speciality Collaboration, UK. COVID-19: Consider cytokine storm syndromes and immunosuppression. Lancet 2020, 395, 1033–1034. [Google Scholar] [CrossRef]
  18. Zhou, D.; Dai, S.M.; Tong, Q. COVID-19: A recommendation to examine the effect of hydroxychloroquine in preventing infection and progression. J. Antimicrob. Chemother. 2020, dkaa114. [Google Scholar] [CrossRef]
  19. Wu, D.; Yang, X.O. TH17 responses in cytokine storm of COVID-19: An emerging target of JAK2 inhibitor Fedratinib. J. Microbiol. Immunol. Infect. 2020. [Google Scholar] [CrossRef]
  20. Boor, P.P.C.; de Ruiter, P.E.; Asmawidjaja, P.S.; Lubberts, E.; van der Laan, L.J.W.; Kwekkeboom, J. JAK-inhibitor tofacitinib suppresses interferon alfa production by plasmacytoid dendritic cells and inhibits arthrogenic and antiviral effects of interferon alfa. Transl. Res. 2017, 188, 67–79. [Google Scholar] [CrossRef]
  21. Iwata-Yoshikawa, N.; Okamura, T.; Shimizu, Y.; Hasegawa, H.; Takeda, M.; Nagata, N. TMPRSS2 Contributes to Virus Spread and Immunopathology in the Airways of Murine Models after Coronavirus Infection. J. Virol. 2019, 93, e01815–e01818. [Google Scholar] [CrossRef] [Green Version]
  22. Hoffmann, M.; Kleine-Weber, H.; Schroeder, S.; Krüger, N.; Herrler, T.; Erichsen, S.; Schiergens, T.S.; Herrler, G.; Wu, N.H.; Nitsche, A.; et al. SARS-CoV-2 Cell Entry Depends on ACE2 and TMPRSS2 and Is Blocked by a Clinically Proven Protease Inhibitor. Cell 2020, 181, 71–280. [Google Scholar] [CrossRef] [PubMed]
  23. Zhou, Y.; Vedantham, P.; Lu, K.; Agudelo, J.; Carrion, R., Jr.; Nunneley, J.W.; Barnard, D.; Pöhlmann, S.; McKerrow, J.H.; Renslo, A.R.; et al. Protease inhibitors targeting coronavirus and filovirus entry. Antivir. Res. 2015, 116, 76–84. [Google Scholar] [CrossRef] [PubMed]
  24. Yamamoto, M.; Matsuyama, S.; Li, X.; Takeda, M.; Kawaguchi, Y.; Inoue, J.I.; Matsuda, Z. Identification of Nafamostat as a Potent Inhibitor of Middle East Respiratory Syndrome Coronavirus S Protein-Mediated Membrane Fusion Using the Split-Protein-Based Cell-Cell Fusion Assay. Antimicrob. Agents Chemother. 2016, 60, 6532–6539. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  25. Prussia, A.; Thepchatri, P.; Snyder, J.P.; Plemper, R.K. Systematic approaches towards the development of host-directed antiviral therapeutics. Int. J. Mol. Sci. 2011, 12, 4027–4052. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  26. Yan, R.; Zhang, Y.; Li, Y.; Xia, L.; Guo, Y.; Zhou, Q. Structural basis for the recognition of the SARS-CoV-2 by full-length human ACE2. Science 2020, 367, 1444–1448. [Google Scholar] [CrossRef] [Green Version]
  27. Walls, A.C.; Park, Y.J.; Tortorici, M.A.; Wall, A.; McGuire, A.T.; Veesler, D. Structure, Function, and Antigenicity of the SARS-CoV-2 Spike Glycoprotein. Cell 2020, 181, 281–292. [Google Scholar] [CrossRef]
  28. Sisk, J.M.; Frieman, M.B.; Machamer, C.E. Coronavirus S protein-induced fusion is blocked prior to hemifusion by Abl kinase inhibitors. J. Gen. Virol. 2018, 99, 619–630. [Google Scholar] [CrossRef]
  29. Coleman, C.M.; Sisk, J.M.; Mingo, R.M.; Nelson, E.A.; White, J.M.; Frieman, M.B. Abelson Kinase Inhibitors Are Potent Inhibitors of Severe Acute Respiratory Syndrome Coronavirus and Middle East Respiratory Syndrome Coronavirus Fusion. J. Virol. 2016, 90, 8924–8933. [Google Scholar] [CrossRef] [Green Version]
  30. Zhang, L.; Lin, D.; Sun, X.; Curth, U.; Drosten, C.; Sauerhering, L.; Becker, S.; Rox, K.; Hilgenfeld, R. Crystal structure of SARS-CoV-2 main protease provides a basis for design of improved α-ketoamide inhibitors. Science 2020, 368, 409–412. [Google Scholar] [CrossRef] [Green Version]
  31. Xia, S.; Zhu, Y.; Liu, M.; Lan, Q.; Xu, W.; Wu, Y.; Ying, T.; Liu, S.; Shi, Z.; Jiang, S.; et al. Fusion mechanism of 2019-nCoV and fusion inhibitors targeting HR1 domain in spike protein. Cell Mol. Immunol. 2020. [Google Scholar] [CrossRef]
  32. Báez-Santos, Y.M.; St John, S.E.; Mesecar, A.D. The SARS-coronavirus papain-like protease: Structure, function and inhibition by designed antiviral compounds. Antivir. Res. 2015, 115, 21–38. [Google Scholar] [CrossRef] [PubMed]
  33. Wu, C.; Liu, Y.; Yang, Y.; Zhang, P.; Zhong, W.; Wang, Y.; Wang, Q.; Xu, Y.; Li, M.; Li, X.; et al. Analysis of therapeutic targets for SARS-CoV-2 and discovery of potential drugs by computational methods. Acta Pharm. Sin. B 2020. [Google Scholar] [CrossRef] [PubMed]
  34. Zhou, Y.; Hou, Y.; Shen, J.; Huang, Y.; Martin, W.; Cheng, F. Network-based drug repurposing for novel coronavirus 2019-nCoV/SARS-CoV-2. Cell Discov. 2020, 6, 14. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  35. Harrison, C. Coronavirus puts drug repurposing on the fast track. Nat. Biotechnol. 2020. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  36. Nukoolkarn, V.; Lee, V.S.; Malaisree, M.; Aruksakulwong, O.; Hannongbua, S. Molecular dynamic simulations analysis of ritonavir and lopinavir as SARS-CoV 3CL(pro) inhibitors. J. Theor. Biol. 2008, 254, 861–867. [Google Scholar] [CrossRef] [PubMed]
  37. Stellbrink, H.J.; Arastéh, K.; Schürmann, D.; Stephan, C.; Dierynck, I.; Smyej, I.; Hoetelmans, R.M.; Truyers, C.; Meyvisch, P.; Jacquemyn, B.; et al. Antiviral activity, pharmacokinetics, and safety of the HIV-1 protease inhibitor TMC310911, coadministered with ritonavir, in treatment-naive HIV-1-infected patients. J. Acquir. Immune Defic. Syndr. 2014, 65, 283–289. [Google Scholar] [CrossRef]
  38. Chan, J.F.; Kok, K.H.; Zhu, Z.; Chu, H.; To, K.K.; Yuan, S.; Yuen, K.Y. Genomic characterization of the 2019 novel human-pathogenic coronavirus isolated from a patient with atypical pneumonia after visiting Wuhan. Emerg. Microbes Infect. 2020, 9, 221–236. [Google Scholar] [CrossRef] [Green Version]
  39. Fehr, A.R.; Perlman, S. Coronaviruses: An overview of their replication and pathogenesis. Methods Mol. Biol. 2015, 1282, 1–23. [Google Scholar]
  40. Daffis, S.; Szretter, K.J.; Schriewer, J.; Li, J.; Youn, S.; Errett, J.; Lin, T.Y.; Schneller, S.; Zust, R.; Dong, H.; et al. 2′-O methylation of the viral mRNA cap evades host restriction by IFIT family members. Nature 2010, 468, 452–456. [Google Scholar] [CrossRef]
  41. Züst, R.; Cervantes-Barragan, L.; Habjan, M.; Maier, R.; Neuman, B.W.; Ziebuhr, J.; Szretter, K.J.; Baker, S.C.; Barchet, W.; Diamond, M.S.; et al. Ribose 2′-O-methylation provides a molecular signature for the distinction of self and non-self mRNA dependent on the RNA sensor Mda5. Nat. Immunol. 2011, 12, 137–143. [Google Scholar] [CrossRef] [Green Version]
  42. García-Sastre, A. 2 methylate or not 2 methylate: Viral evasion of the type I interferon response. Nat. Immunol. 2011, 12, 114–115. [Google Scholar] [CrossRef] [PubMed]
  43. Bouvet, M.; Debarnot, C.; Imbert, I.; Selisko, B.; Snijder, E.J.; Canard, B.; Decroly, E. In vitro reconstitution of SARS-coronavirus mRNA cap methylation. PLoS Pathog. 2010, 6, e1000863. [Google Scholar] [CrossRef]
  44. Lugari, A.; Betzi, S.; Decroly, E.; Bonnaud, E.; Hermant, A.; Guillemot, J.C.; Debarnot, C.; Borg, J.P.; Bouvet, M.; Canard, B.; et al. Molecular mapping of the RNA Cap 2′-O-methyltransferase activation interface between severe acute respiratory syndrome coronavirus nsp10 and nsp16. J. Biol. Chem. 2010, 285, 33230–33241. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  45. Chen, Y.; Su, C.; Ke, M.; Jin, X.; Xu, L.; Zhang, Z.; Wu, A.; Sun, Y.; Yang, Z.; Tien, P.; et al. Biochemical and structural insights into the mechanisms of SARS coronavirus RNA ribose 2′-O-methylation by nsp16/nsp10 protein complex. PLoS Pathog. 2011, 7, e1002294. [Google Scholar] [CrossRef] [Green Version]
  46. Decroly, E.; Debarnot, C.; Ferron, F.; Bouvet, M.; Coutard, B.; Imbert, I.; Gluais, L.; Papageorgiou, N.; Sharff, A.; Bricogne, G.; et al. Crystal structure and functional analysis of the SARS-coronavirus RNA cap 2′-O-methyltransferase nsp10/nsp16 complex. PLoS Pathog. 2011, 7, e1002059. [Google Scholar] [CrossRef] [Green Version]
  47. Ke, M.; Chen, Y.; Wu, A.; Sun, Y.; Su, C.; Wu, H.; Jin, X.; Tao, J.; Wang, Y.; Ma, X.; et al. Short peptides derived from the interaction domain of SARS coronavirus nonstructural protein nsp10 can suppress the 2′-O-methyltransferase activity of nsp10/nsp16 complex. Virus Res. 2012, 167, 322–328. [Google Scholar] [CrossRef]
  48. Jin, X.; Chen, Y.; Sun, Y.; Zeng, C.; Wang, Y.; Tao, J.; Wu, A.; Yu, X.; Zhang, Z.; Tian, J.; et al. Characterization of the guanine-N7 methyltransferase activity of coronavirus nsp14 on nucleotide GTP. Virus Res. 2013, 176, 45–52. [Google Scholar] [CrossRef]
  49. Chen, Y.; Tao, J.; Sun, Y.; Wu, A.; Su, C.; Gao, G.; Cai, H.; Qiu, S.; Wu, Y.; Ahola, T.; et al. Structure-function analysis of severe acute respiratory syndrome coronavirus RNA cap guanine-N7-methyltransferase. J. Virol. 2013, 87, 6296–6305. [Google Scholar] [CrossRef] [Green Version]
  50. Sun, Y.; Wang, Z.; Tao, J.; Wang, Y.; Wu, A.; Yang, Z.; Wang, K.; Shi, L.; Chen, Y.; Guo, D. Yeast-based assays for the high-throughput screening of inhibitors of coronavirus RNA cap guanine-N7-methyltransferase. Antivir. Res. 2014, 104, 156–164. [Google Scholar] [CrossRef]
  51. Menachery, V.D.; Yount, B.L., Jr.; Josset, L.; Gralinski, L.E.; Scobey, T.; Agnihothram, S.; Katze, M.G.; Baric, R.S. Attenuation and restoration of severe acute respiratory syndrome coronavirus mutant lacking 2′-o-methyltransferase activity. J. Virol. 2014, 88, 4251–4264. [Google Scholar] [CrossRef] [Green Version]
  52. Wishart, D.S.; Feunang, Y.D.; Guo, A.C.; Lo, E.J.; Marcu, A.; Grant, J.R.; Sajed, T.; Johnson, D.; Li, C.; Sayeeda, Z.; et al. DrugBank 5.0: A major update to the DrugBank database for 2018. Nucleic Acids Res. 2018, 46, D1074–D1082. [Google Scholar] [CrossRef] [PubMed]
  53. Krieger, E.; Vriend, G. YASARA View—Molecular graphics for all devices—From smartphones to workstations. Bioinformatics 2014, 30, 2981–2982. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  54. Trott, O.; Olson, A.J. AutoDock Vina: Improving the speed and accuracy of docking with a new scoring function, efficient optimization, and multithreading. J. Comput. Chem. 2010, 31, 455–461. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  55. Ruiz-Torres, V.; Losada-Echeberría, M.; Herranz-López, M.; Barrajón-Catalán, E.; Galiano, V.; Micol, V.; Encinar, J.A. New Mammalian Target of Rapamycin (mTOR) Modulators Derived from Natural Product Databases and Marine Extracts by Using Molecular Docking Techniques. Mar. Drugs. 2018, 16, 385. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  56. Seeliger, D.; de Groot, B.L. Ligand docking and binding site analysis with PyMOL and Autodock/Vina. J. Comput. Aided Mol. Des. 2010, 24, 417–422. [Google Scholar] [CrossRef] [Green Version]
  57. Lionta, E.; Spyrou, G.; Vassilatis, D.K.; Cournia, Z. Structure-based virtual screening for drug discovery: Principles, applications and recent advances. Curr. Top Med. Chem. 2014, 14, 1923–1938. [Google Scholar] [CrossRef]
  58. Krieger, E.; Nielsen, J.E.; Spronk, C.A.; Vriend, G. Fast empirical pKa prediction by Ewald summation. J. Mol. Graph. Model. 2006, 25, 481–486. [Google Scholar] [CrossRef]
  59. Cuyàs, E.; Verdura, S.; Lozano-Sánchez, J.; Viciano, I.; Llorach-Parés, L.; Nonell-Canals, A.; Bosch-Barrera, J.; Brunet, J.; Segura-Carretero, A.; Sanchez-Martinez, M.; et al. The extra virgin olive oil phenolic oleacein is a dual substrate-inhibitor of catechol-O-methyltransferase. Food Chem. Toxicol. 2019, 128, 35–45. [Google Scholar] [CrossRef]
  60. Verdura, S.; Cuyàs, E.; Cortada, E.; Brunet, J.; Lopez-Bonet, E.; Martin-Castillo, B.; Bosch-Barrera, J.; Encinar, J.A.; Menendez, J.A. Resveratrol targets PD-L1 glycosylation and dimerization to enhance antitumor T-cell immunity. Aging 2020, 12, 8–34. [Google Scholar] [CrossRef] [Green Version]
  61. Cuyàs, E.; Gumuzio, J.; Verdura, S.; Brunet, J.; Bosch-Barrera, J.; Martin-Castillo, B.; Alarcón, T.; Encinar, J.A.; Martin, Á.G.; Menendez, J.A. The LSD1 inhibitor iadademstat (ORY-1001) targets SOX2-driven breast cancer stem cells: A potential epigenetic therapy in luminal-B and HER2-positive breast cancer subtypes. Aging 2020, 12, 4794. [Google Scholar] [CrossRef]
  62. Genheden, S.; Ryde, U. The MM/PBSA and MM/GBSA methods to estimate ligand-binding affinities. Expert Opin. Drug Discov. 2015, 10, 449–461. [Google Scholar] [CrossRef] [PubMed]
  63. Wang, C.; Nguyen, P.H.; Pham, K.; Huynh, D.; Le, T.B.; Wang, H.; Ren, P.; Luo, R. Calculating protein-ligand binding affinities with MMPBSA: Method and error analysis. J. Comput. Chem. 2016, 37, 2436–2446. [Google Scholar] [CrossRef] [Green Version]
  64. Hyde, J.L.; Diamond, M.S. Innate immune restriction and antagonism of viral RNA lacking 2′-O methylation. Virology 2015, 479, 66–74. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  65. Srinivasan, S.; Cui, H.; Gao, Z.; Liu, M.; Lu, S.; Mkandawire, W.; Narykov, O.; Sun, M.; Korkin, D. Structural Genomics of SARS-CoV-2 Indicates Evolutionary Conserved Functional Regions of Viral Proteins. Viruses 2020, 12, 360. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  66. Gordon, D.E.; Jang, G.M.; Bouhaddou, M.; Xu, J.; Obernier, K.; White, K.M.; O’Meara, M.J.; Rezelj, V.V.; Guo, J.Z.; Swaney, D.L.; et al. A SARS-CoV-2 protein interaction map reveals targets for drug repurposing. Nature 2020, 1–13. [Google Scholar] [CrossRef] [PubMed]
  67. Hebner, C.M.; Han, B.; Brendza, K.M.; Nash, M.; Sulfab, M.; Tian, Y.; Hung, M.; Fung, W.; Vivian, R.W.; Trenkle, J.; et al. The HCV non-nucleoside inhibitor Tegobuvir utilizes a novel mechanism of action to inhibit NS5B polymerase function. PLoS ONE 2012, 7, e39163. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  68. Wyles, D.L.; Rodriguez-Torres, M.; Lawitz, E.; Shiffman, M.L.; Pol, S.; Herring, R.W.; Massetto, B.; Kanwar, B.; Trenkle, J.D.; Pang, P.S.; et al. All-oral combination of ledipasvir, vedroprevir, tegobuvir, and ribavirin in treatment-naïve patients with genotype 1 HCV infection. Hepatology 2014, 60, 56–64. [Google Scholar] [CrossRef]
  69. Vliegen, I.; Paeshuyse, J.; Zhong, W.; Neyts, J. In vitro combinations containing Tegobuvir are highly efficient in curing cells from HCV replicon and in delaying/preventing the development of drug resistance. Antivir. Res. 2015, 120, 112–121. [Google Scholar] [CrossRef]
  70. Liu, M.; Xu, Q.; Guo, S.; Zuo, R.; Hong, Y.; Luo, Y.; Li, Y.; Gong, P.; Liu, Y. Design, synthesis, and structure-activity relationships of novel imidazo [4,5-c] pyridine derivatives as potent non-nucleoside inhibitors of hepatitis C virus NS5B. Bioorg Med Chem. 2018, 26, 2621–2631. [Google Scholar] [CrossRef]
  71. Ruan, Z.; Liu, C.; Guo, Y.; He, Z.; Huang, X.; Jia, X.; Yang, T. Potential Inhibitors Targeting RNA-Dependent RNA Polymerase Activity (NSP12) of SARS-CoV-2. Preprints 2020, 2020030024. [Google Scholar] [CrossRef]
  72. Khalili, J.S.; Zhu, H.; Mak, A.; Yan, Y.; Zhu, Y. Novel coronavirus treatment with ribavirin: Groundwork for evaluation concerning COVID-19. J. Med. Virol. 2020. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  73. Elfiky, A.A. Ribavirin, Remdesivir, Sofosbuvir, Galidesivir, and Tenofovir against SARS-CoV-2 RNA dependent RNA polymerase (RdRp): A molecular docking study. Life Sci. 2020, 117592. [Google Scholar] [CrossRef] [PubMed]
  74. Ferron, F.; Subissi, L.; Silveira De Morais, A.T.; Le, N.T.T.; Sevajol, M.; Gluais, L.; Decroly, E.; Vonrhein, C.; Bricogne, G.; Canard, B.; et al. Structural and molecular basis of mismatch correction and ribavirin excision from coronavirus RNA. Proc. Natl. Acad. Sci. USA 2018, 115, E162–E171. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  75. Ostenfeld, M.S.; Fehrenbacher, N.; Høyer-Hansen, M.; Thomsen, C.; Farkas, T.; Jäättelä, M. Effective tumor cell death by sigma-2 receptor ligand siramesine involves lysosomal leakage and oxidative stress. Cancer Res. 2005, 65, 8975–8983. [Google Scholar] [CrossRef] [Green Version]
  76. Ostenfeld, M.S.; Høyer-Hansen, M.; Bastholm, L.; Fehrenbacher, N.; Olsen, O.D.; Groth-Pedersen, L.; Puustinen, P.; Kirkegaard-Sørensen, T.; Nylandsted, J.; Farkas, T.; et al. Anti-cancer agent siramesine is a lysosomotropic detergent that induces cytoprotective autophagosome accumulation. Autophagy 2008, 4, 487–499. [Google Scholar] [CrossRef] [Green Version]
  77. Zhitomirsky, B.; Yunaev, A.; Kreiserman, R.; Kaplan, A.; Stark, M.; Assaraf, Y.G. Lysosomotropic drugs activate TFEB via lysosomal membrane fluidization and consequent inhibition of mTORC1 activity. Cell Death Dis. 2018, 9, 1191. [Google Scholar] [CrossRef]
  78. Ma, S.; Henson, E.S.; Chen, Y.; Gibson, S.B. Ferroptosis is induced following siramesine and lapatinib treatment of breast cancer cells. Cell Death Dis. 2016, 7, e2307. [Google Scholar] [CrossRef] [Green Version]
  79. Villalpando-Rodriguez, G.E.; Blankstein, A.R.; Konzelman, C.; Gibson, S.B. Lysosomal Destabilizing Drug Siramesine and the Dual Tyrosine Kinase Inhibitor Lapatinib Induce a Synergistic Ferroptosis through Reduced Heme Oxygenase-1 (HO-1) Levels. Oxid. Med. Cell. Longev. 2019, 2019, 9561281. [Google Scholar] [CrossRef]
  80. Gassen, N.C.; Niemeyer, D.; Muth, D.; Corman, V.M.; Martinelli, S.; Gassen, A.; Hafner, K.; Papies, J.; Mösbauer, K.; Zellner, A.; et al. SKP2 attenuates autophagy through Beclin1-ubiquitination and its inhibition reduces MERS-Coronavirus infection. Nat. Commun. 2019, 10, 5770. [Google Scholar] [CrossRef]
  81. Hirata, Y.; Yamamoto, H.; Atta, M.S.; Mahmoud, S.; Oh-hashi, K.; Kiuchi, K. Chloroquine inhibits glutamate-induced death of a neuronal cell line by reducing reactive oxygen species through sigma-1 receptor. J. Neurochem. 2011, 119, 839–847. [Google Scholar] [CrossRef]
  82. Tesei, A.; Cortesi, M.; Zamagni, A.; Arienti, C.; Pignatta, S.; Zanoni, M.; Paolillo, M.; Curti, D.; Rui, M.; Rossi, D.; et al. Sigma Receptors as Endoplasmic Reticulum Stress “Gatekeepers” and their Modulators as Emerging New Weapons in the Fight against Cancer. Front. Pharmacol. 2018, 9, 711. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  83. Wang, H.; Yang, P.; Liu, K.; Guo, F.; Zhang, Y.; Zhang, G.; Jiang, C. SARS coronavirus entry into host cells through a novel clathrin- and caveolae-independent endocytic pathway. Cell Res. 2008, 18, 290–301. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  84. Yang, N.; Shen, H.M. Targeting the Endocytic Pathway and Autophagy Process as a Novel Therapeutic Strategy in COVID-19. Int. J. Biol. Sci. 2020, 16, 1724–1731. [Google Scholar] [CrossRef] [PubMed]
  85. Fung, T.S.; Liu, D.X. Coronavirus infection, ER stress, apoptosis and innate immunity. Front. Microbiol. 2014, 5, 296. [Google Scholar] [CrossRef] [Green Version]
  86. Fung, T.S.; Liao, Y.; Liu, D.X. Regulation of Stress Responses and Translational Control by Coronavirus. Viruses 2016, 8, 184. [Google Scholar] [CrossRef]
  87. Shi, C.S.; Nabar, N.R.; Huang, N.N.; Kehrl, J.H. SARS-Coronavirus Open Reading Frame-8b triggers intracellular stress pathways and activates NLRP3 inflammasomes. Cell Death Discov. 2019, 5, 101. [Google Scholar] [CrossRef] [Green Version]
  88. Heading, C. Siramesine H Lundbeck. Curr. Opin. Investig. Drugs. 2001, 2, 266–270. [Google Scholar]
  89. Søby, K.K.; Mikkelsen, J.D.; Meier, E.; Thomsen, C. Lu 28-179 labels a sigma(2)-site in rat and human brain. Neuropharmacology 2002, 43, 95–100. [Google Scholar] [CrossRef]
  90. Schmidt, H.R.; Betz, R.M.; Dror, R.O.; Kruse, A.C. Structural basis for σ1 receptor ligand recognition. Nat. Struct. Mol. Biol. 2018, 25, 981–987. [Google Scholar] [CrossRef]
  91. Sahn, J.J.; Mejia, G.L.; Ray, P.R.; Martin, S.F.; Price, T.J. Sigma 2 Receptor/Tmem97 Agonists Produce Long Lasting Antineuropathic Pain Effects in Mice. ACS Chem. Neurosci. 2017, 8, 1801–1811. [Google Scholar] [CrossRef] [Green Version]
  92. Georgiadis, M.O.; Karoutzou, O.; Foscolos, A.S.; Papanastasiou, I. Sigma Receptor (σR) Ligands with Antiproliferative and Anticancer Activity. Molecules 2017, 22, 1408. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  93. Holland, S.J.; Pan, A.; Franci, C.; Hu, Y.; Chang, B.; Li, W.; Duan, M.; Torneros, A.; Yu, J.; Heckrodt, T.J.; et al. R428, a selective small molecule inhibitor of Axl kinase, blocks tumor spread and prolongs survival in models of metastatic breast cancer. Cancer Res. 2010, 70, 1544–1554. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  94. Sadahiro, H.; Kang, K.D.; Gibson, J.T.; Minata, M.; Yu, H.; Shi, J.; Chhipa, R.; Chen, Z.; Lu, S.; Simoni, Y.; et al. Activation of the Receptor Tyrosine Kinase AXL Regulates the Immune Microenvironment in Glioblastoma. Cancer Res. 2018, 78, 3002–3013. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  95. Ludwig, K.F.; Du, W.; Sorrelle, N.B.; Wnuk-Lipinska, K.; Topalovski, M.; Toombs, J.E.; Cruz, V.H.; Yabuuchi, S.; Rajeshkumar, N.V.; Maitra, A.; et al. Small-Molecule Inhibition of Axl Targets Tumor Immune Suppression and Enhances Chemotherapy in Pancreatic Cancer. Cancer Res. 2018, 78, 246–255. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  96. Landolt, L.; Furriol, J.; Babickova, J.; Ahmed, L.; Eikrem, Ø.; Skogstrand, T.; Scherer, A.; Suliman, S.; Leh, S.; Lorens, J.B.; et al. AXL targeting reduces fibrosis development in experimental unilateral ureteral obstruction. Physiol. Rep. 2019, 7, e14091. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  97. Tutusaus, A.; de Gregorio, E.; Cucarull, B.; Cristóbal, H.; Aresté, C.; Graupera, I.; Coll, M.; Colell, A.; Gausdal, G.; Lorens, J.B.; et al. A Functional Role of GAS6/TAM in Nonalcoholic Steatohepatitis Progression Implicates AXL as Therapeutic Target. Cell Mol. Gastroenterol. Hepatol. 2020, 9, 349–368. [Google Scholar] [CrossRef]
  98. Rothlin, C.V.; Ghosh, S.; Zuniga, E.I.; Oldstone, M.B.; Lemke, G. TAM receptors are pleiotropic inhibitors of the innate immune response. Cell 2007, 131, 1124–1136. [Google Scholar] [CrossRef] [Green Version]
  99. Meertens, L.; Labeau, A.; Dejarnac, O.; Cipriani, S.; Sinigaglia, L.; Bonnet-Madin, L.; Le Charpentier, T.; Hafirassou, M.L.; Zamborlini, A.; Cao-Lormeau, V.M.; et al. Axl Mediates ZIKA Virus Entry in Human Glial Cells and Modulates Innate Immune Responses. Cell Rep. 2017, 18, 324–333. [Google Scholar] [CrossRef]
  100. Persaud, M.; Martinez-Lopez, A.; Buffone, C.; Porcelli, S.A.; Diaz-Griffero, F. Infection by Zika viruses requires the transmembrane protein AXL, endocytosis and low pH. Virology 2018, 518, 301–312. [Google Scholar] [CrossRef]
  101. Chen, F.; Song, Q.; Yu, Q. Axl inhibitor R428 induces apoptosis of cancer cells by blocking lysosomal acidification and recycling independent of Axl inhibition. Am. J. Cancer Res. 2018, 8, 1466–1482. [Google Scholar]
  102. Lotsberg, M.L.; Wnuk-Lipinska, K.; Terry, S.; Tan, T.Z.; Lu, N.; Trachsel-Moncho, L.; Røsland, G.V.; Siraji, M.I.; Hellesøy, M.; Rayford, A.; et al. AXL Targeting Abrogates Autophagic Flux and Induces Immunogenic Cell Death in Drug-Resistant Cancer Cells. J. Thorac. Oncol. 2020. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  103. Burness, C.B. Sonidegib: First Global Approval. Drugs 2015, 75, 1559–1566. [Google Scholar] [CrossRef] [PubMed]
  104. Collier, N.J.; Ali, F.R.; Lear, J.T. The safety and efficacy of sonidegib for the treatment of locally advanced basal cell carcinoma. Expert Rev. Anticancer Ther. 2016, 16, 1011–1018. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  105. Jain, S.; Song, R.; Xie, J. Sonidegib: Mechanism of action, pharmacology, and clinical utility for advanced basal cell carcinomas. Onco Targets Ther. 2017, 10, 1645–1653. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  106. Smelkinson, M.G. The Hedgehog Signaling Pathway Emerges as a Pathogenic Target. J. Dev. Biol. 2017, 5, 14. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  107. Roskoski, R., Jr. Janus kinase (JAK) inhibitors in the treatment of inflammatory and neoplastic diseases. Pharmacol. Res. 2016, 111, 784–803. [Google Scholar] [CrossRef]
  108. Gong, X.; Darpo, B.; Xue, H.; Punwani, N.; He, K.; Barbour, A.M.; Epstein, N.; Landman, R.; Chen, X.; Yeleswaram, S. Evaluation of Clinical Cardiac Safety of Itacitinib, a JAK1 Inhibitor, in Healthy Participants. Clin. Pharmacol. Drug Dev. 2019. [Google Scholar] [CrossRef]
  109. Bissonnette, R.; Luchi, M.; Fidelus-Gort, R.; Jackson, S.; Zhang, H.; Flores, R.; Newton, R.; Scherle, P.; Yeleswaram, S.; Chen, X.; et al. A randomized, double-blind, placebo-controlled, dose-escalation study of the safety and efficacy of INCB039110, an oral janus kinase 1 inhibitor, in patients with stable, chronic plaque psoriasis. J. Dermatolog. Treat. 2016, 27, 332–338. [Google Scholar] [CrossRef]
  110. Kettle, J.G.; Åstrand, A.; Catley, M.; Grimster, N.P.; Nilsson, M.; Su, Q.; Woessner, R. Inhibitors of JAK-family kinases: An update on the patent literature 2013–2015, part 1. Expert Opin. Ther. Pat. 2017, 27, 127–143. [Google Scholar] [CrossRef]
  111. Schwartz, D.M.; Kanno, Y.; Villarino, A.; Ward, M.; Gadina, M.; O’Shea, J.J. JAK inhibition as a therapeutic strategy for immune and inflammatory diseases. Nat. Rev. Drug Discov. 2017, 16, 843–862. [Google Scholar] [CrossRef]
  112. Virtanen, A.T.; Haikarainen, T.; Raivola, J.; Silvennoinen, O. Selective JAKinibs: Prospects in Inflammatory and Autoimmune Diseases. BioDrugs 2019, 33, 15–32. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  113. Gadina, M.; Le, M.T.; Schwartz, D.M.; Silvennoinen, O.; Nakayamada, S.; Yamaoka, K.; O’Shea, J.J. Janus kinases to jakinibs: From basic insights to clinical practice. Rheumatology 2019, 58, i4–i16. [Google Scholar] [CrossRef] [PubMed]
  114. Park, J.H.; Frigault, M.J.; Maziarz, R.T.; Naim, A.; Burke, L.; Tian, C.; Porter, D.L. Trial in Progress: A Phase 2, Single-Arm, Open-Label Study of Itacitinib (ITA) for the Prevention of Chimeric Antigen Receptor (CAR) T-Cell–Induced Cytokine Release Syndrome (CRS). Biol. Blood Marrow Transplant. 2020, 16, S269. [Google Scholar] [CrossRef]
  115. Shimabukuro-Vornhagen, A.; Gödel, P.; Subklewe, M.; Stemmler, H.J.; Schlößer, H.A.; Schlaak, M.; Kochanek, M.; Böll, B.; von Bergwelt-Baildon, M.S. Cytokine release syndrome. J. Immunother. Cancer. 2018, 6, 56. [Google Scholar] [CrossRef] [Green Version]
  116. Chen, G.; Wu, D.; Guo, W.; Cao, Y.; Huang, D.; Wang, H.; Wang, T.; Zhang, X.; Chen, H.; Yu, H.; et al. Clinical and immunologic features in severe and moderate Coronavirus Disease 2019. J. Clin. Investig. 2020, 130, 2620–2629. [Google Scholar] [CrossRef] [Green Version]
  117. Pedersen, S.F.; Ho, Y.C. SARS-CoV-2: A Storm is Raging. J. Clin Investig. 2020, 130, 137647. [Google Scholar] [CrossRef]
  118. Favalli, E.G.; Ingegnoli, F.; De Lucia, O.; Cincinelli, G.; Cimaz, R.; Caporali, R. COVID-19 infection and rheumatoid arthritis: Faraway, so close! Autoimmun. Rev. 2020, 19, 102523. [Google Scholar] [CrossRef]
  119. Richardson, P.; Griffin, I.; Tucker, C.; Smith, D.; Oechsle, O.; Phelan, A.; Stebbing, J. Baricitinib as potential treatment for 2019-nCoV acute respiratory disease. Lancet 2020, 395, e30–e31. [Google Scholar] [CrossRef] [Green Version]
  120. Bell, I.M.; Gallicchio, S.N.; Wood, M.R.; Quigley, A.G.; Stump, C.A.; Zartman, C.B.; Fay, J.F.; Li, C.C.; Lynch, J.J.; Moore, E.L.; et al. Discovery of MK-3207: A Highly Potent, Orally Bioavailable CGRP Receptor Antagonist. ACS Med. Chem. Lett. 2010, 1, 24–29. [Google Scholar] [CrossRef] [Green Version]
  121. Hewitt, D.J.; Aurora, S.K.; Dodick, D.W.; Goadsby, P.J.; Ge, Y.J.; Bachman, R.; Taraborelli, D.; Fan, X.; Assaid, C.; Lines, C.; et al. Randomized controlled trial of the CGRP receptor antagonist MK-3207 in the acute treatment of migraine. Cephalalgia 2011, 31, 712–722. [Google Scholar] [CrossRef]
  122. Crowley, B.M.; Stump, C.A.; Nguyen, D.N.; Potteiger, C.M.; McWherter, M.A.; Paone, D.V.; Quigley, A.G.; Bruno, J.G.; Cui, D.; Culberson, J.C.; et al. Novel oxazolidinone calcitonin gene-related peptide (CGRP) receptor antagonists for the acute treatment of migraine. Bioorg. Med. Chem. Lett. 2015, 25, 4777–4781. [Google Scholar] [CrossRef] [PubMed]
  123. Negro, A.; Martelletti, P. Gepants for the treatment of migraine. Expert Opin. Investig. Drugs. 2019, 28, 555–567. [Google Scholar] [CrossRef] [PubMed]
  124. Tepper, S.J. History and Review of anti-Calcitonin Gene-Related Peptide (CGRP) Therapies: From Translational Research to Treatment. Headache 2018, 58 (Suppl. 3), 238–275. [Google Scholar] [CrossRef] [PubMed]
  125. De Vries, T.; Villalón, C.M.; MaassenVanDenBrink, A. Pharmacological treatment of migraine: CGRP and 5-HT beyond the triptans. Pharmacol. Ther. 2020, 107528. [Google Scholar] [CrossRef]
  126. Li, Y.C.; Bai, W.Z.; Hashikawa, T. The neuroinvasive potential of SARS-CoV2 may play a role in the respiratory failure of COVID-19 patients. J. Med. Virol. 2020, 92. [Google Scholar] [CrossRef]
  127. Baig, A.M.; Khaleeq, A.; Ali, U.; Syeda, H. Evidence of the COVID-19 Virus Targeting the CNS: Tissue Distribution, Host-Virus Interaction, and Proposed Neurotropic Mechanisms. ACS Chem. Neurosci. 2020, 11, 995–998. [Google Scholar] [CrossRef] [Green Version]
  128. Wu, Y.; Xu, X.; Chen, Z.; Duan, J.; Hashimoto, K.; Yang, L.; Liu, C.; Yang, C. Nervous system involvement after infection with COVID-19 and other coronaviruses. Brain Behav. Immun. 2020. [Google Scholar] [CrossRef]
  129. Han, H.; Xie, L.; Liu, R.; Yang, J.; Liu, F.; Wu, K.; Chen, L.; Hou, W.; Feng, Y.; Zhu, C. Analysis of heart injury laboratory parameters in 273 COVID-19 patients in one hospital in Wuhan, China. J. Med. Virol. 2020. [Google Scholar] [CrossRef]
  130. Fried, J.A.; Ramasubbu, K.; Bhatt, R.; Topkara, V.K.; Clerkin, K.J.; Horn, E.; Rabbani, L.; Brodie, D.; Jain, S.S.; Kirtane, A.; et al. The Variety of Cardiovascular Presentations of COVID-19. Circulation 2020. [Google Scholar] [CrossRef] [Green Version]
  131. Bansal, M. Cardiovascular disease and COVID-19. Diabetes Metab. Syndr. 2020, 14, 247–250. [Google Scholar] [CrossRef]
  132. Yang, J.; Zheng, Y.; Gou, X.; Pu, K.; Chen, Z.; Guo, Q.; Ji, R.; Wang, H.; Wang, Y.; Zhou, Y. Prevalence of comorbidities in the novel Wuhan coronavirus (COVID-19) infection: A systematic review and meta-analysis. Int. J. Infect. Dis. 2020, 94, 91–95. [Google Scholar] [CrossRef] [PubMed]
  133. Ding, Y.; He, L.; Zhang, Q.; Huang, Z.; Che, X.; Hou, J.; Wang, H.; Shen, H.; Qiu, L.; Li, Z.; et al. Organ distribution of severe acute respiratory syndrome (SARS) associated coronavirus (SARS-CoV) in SARS patients: Implications for pathogenesis and virus transmission pathways. J. Pathol. 2004, 203, 622–630. [Google Scholar] [CrossRef] [PubMed]
  134. Netland, J.; Meyerholz, D.K.; Moore, S.; Cassell, M.; Perlman, S. Severe acute respiratory syndrome coronavirus infection causes neuronal death in the absence of encephalitis in mice transgenic for human ACE2. J. Virol. 2008, 82, 7264–7275. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  135. Le Mével, J.C.; Lancien, F.; Mimassi, N.; Conlon, J.M. Brain neuropeptides in central ventilatory and cardiovascular regulation in trout. Front. Endocrinol. 2012, 3, 124. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  136. Dakhama, A.; Larsen, G.L.; Gelfand, E.W. Calcitonin gene-related peptide: Role in airway homeostasis. Curr. Opin. Pharmacol. 2004, 4, 215–220. [Google Scholar] [CrossRef] [PubMed]
  137. Rubio-Beltrán, E.; van den Brink, A.M. Understanding CGRP and Cardiovascular Risk. Handb. Exp. Pharmacol. 2019, 255, 131–140. [Google Scholar]
  138. Favoni, V.; Giani, L.; Al-Hassany, L.; Asioli, G.M.; Butera, C.; de Boer, I.; Guglielmetti, M.; Koniari, C.; Mavridis, T.; Vaikjärv, M.; et al. European Headache Federation School of Advanced Studies (EHF-SAS). CGRP and migraine from a cardiovascular point of view: What do we expect from blocking CGRP? J. Headache Pain. 2019, 20, 27. [Google Scholar] [CrossRef]
  139. Nahama, A.; Ramachandran, R.; Cisternas, A.F.; Ji, H. The role of afferent pulmonary innervation in poor prognosis of acute respiratory distress syndrome in COVID-19 patients and proposed use of resiniferatoxin (RTX) to improve patient outcomes in advanced disease state: A review. Med. Drug Discov. 2020. [Google Scholar] [CrossRef]
  140. Johnson, D.S.; Stiff, C.; Lazerwith, S.E.; Kesten, S.R.; Fay, L.K.; Morris, M.; Beidler, D.; Liimatta, M.B.; Smith, S.E.; Dudley, D.T.; et al. Discovery of PF-04457845: A Highly Potent, Orally Bioavailable, and Selective Urea FAAH Inhibitor. ACS Med. Chem. Lett. 2011, 2, 91–96. [Google Scholar] [CrossRef]
  141. Ahn, K.; Smith, S.E.; Liimatta, M.B.; Beidler, D.; Sadagopan, N.; Dudley, D.T.; Young, T.; Wren, P.; Zhang, Y.; Swaney, S.; et al. Mechanistic and pharmacological characterization of PF-04457845: A highly potent and selective fatty acid amide hydrolase inhibitor that reduces inflammatory and noninflammatory pain. J. Pharmacol. Exp. Ther. 2011, 338, 114–124. [Google Scholar] [CrossRef]
  142. D’Souza, D.C.; Cortes-Briones, J.; Creatura, G.; Bluez, G.; Thurnauer, H.; Deaso, E.; Bielen, K.; Surti, T.; Radhakrishnan, R.; Gupta, A.; et al. Efficacy and safety of a fatty acid amide hydrolase inhibitor (PF-04457845) in the treatment of cannabis withdrawal and dependence in men: A double-blind, placebo-controlled, parallel group, phase 2a single-site randomised controlled trial. Lancet Psychiatr. 2019, 6, 35–45. [Google Scholar] [CrossRef]
  143. Berry, H.; Coquelin, J.P.; Gordon, A.; Seymour, D. Antrafenine, naproxen and placebo in osteoarthritis: A comparative study. Br. J. Rheumatol. 1983, 22, 89–94. [Google Scholar] [CrossRef] [PubMed]
  144. Palanki, M.S.; Akiyama, H.; Campochiaro, P.; Cao, J.; Chow, C.P.; Dellamary, L.; Doukas, J.; Fine, R.; Gritzen, C.; Hood, J.D.; et al. Development of prodrug 4-chloro-3-(5-methyl-3-{[4-(2-pyrrolidin-1-ylethoxy)phenyl]amino}-1,2,4-benzotriazin-7-yl)phenyl benzoate (TG100801): A topically administered therapeutic candidate in clinical trials for the treatment of age-related macular degeneration. J. Med. Chem. 2008, 51, 1546–1559. [Google Scholar] [PubMed]
  145. Sharma, S.; Mulik, S.; Kumar, N.; Suryawanshi, A.; Rouse, B.T. An anti-inflammatory role of VEGFR2/Src kinase inhibitor in herpes simplex virus 1-induced immunopathology. J. Virol. 2011, 85, 5995–6007. [Google Scholar] [CrossRef] [Green Version]
  146. Abidi, A.; Shukla, P.; Ahmad, A. Lifitegrast: A novel drug for treatment of dry eye disease. J. Pharmacol. Pharmacother. 2016, 7, 194–198. [Google Scholar] [CrossRef] [Green Version]
  147. Paton, D.M. Lifitegrast: First LFA-1/ICAM-1 antagonist for treatment of dry eye disease. Drugs Today 2016, 52, 485–493. [Google Scholar] [CrossRef]
  148. Xia, J.; Tong, J.; Liu, M.; Shen, Y.; Guo, D. Evaluation of coronavirus in tears and conjunctival secretions of patients with SARS-CoV-2 infection. J. Med. Virol. 2020. [Google Scholar] [CrossRef] [Green Version]
  149. Guan, W.J.; Ni, Z.Y.; Hu, Y.; Liang, W.H.; Ou, C.Q.; He, J.X.; Liu, L.; Shan, H.; Lei, C.L.; Hui, D.S.C.; et al. China Medical Treatment Expert Group for Covid-19. Clinical Characteristics of Coronavirus Disease 2019 in China. N. Engl. J. Med. 2020. [Google Scholar] [CrossRef]
  150. Seah, I.; Agrawal, R. Can the Coronavirus Disease 2019 (COVID-19) Affect the Eyes? A Review of Coronaviruses and Ocular Implications in Humans and Animals. Ocul. Immunol. Inflamm. 2020, 28, 391–395. [Google Scholar] [CrossRef]
  151. Chen, L.; Deng, C.; Chen, X.; Zhang, X.; Chen, B.; Yu, H.; Qin, Y.; Xiao, K.; Zhnag, H.; Sun, X. Ocular manifestations and clinical characteristics of 534 cases of COVID-19 in China: A cross-sectional study. MedRxiv 2020. [Google Scholar] [CrossRef] [Green Version]
  152. Donnelly, R.F.; McCarron, P.A.; Tunney, M.M. Antifungal photodynamic therapy. Microbiol. Res. 2008, 163, 1–12. [Google Scholar] [CrossRef] [PubMed]
  153. Nikolaeva-Glomb, L.; Mukova, L.; Nikolova, N.; Kussovski, V.; Doumanova, L.; Mantareva, V.; Angelov, I.; Wöhrle, D.; Galabov, A.S. Photodynamic Effect of some Phthalocyanines on Enveloped and Naked Viruses. Acta Virol. 2017, 61, 341–346. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  154. Remichkova, M.; Mukova, L.; Nikolaeva-Glomb, L.; Nikolova, N.; Doumanova, L.; Mantareva, V.; Angelov, I.; Kussovski, V.; Galabov, A.S. Virus inactivation under the photodynamic effect of phthalocyanine zinc(II) complexes. Z. Naturforsch. C J. Biosci. 2017, 72, 123–128. [Google Scholar] [CrossRef] [PubMed]
  155. Korneev, D.; Kurskaya, O.; Sharshov, K.; Eastwood, J.; Strakhovskaya, M. Ultrastructural Aspects of Photodynamic Inactivation of Highly Pathogenic Avian H5N8 Influenza Virus. Viruses 2019, 11, 955. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  156. Kampf, G.; Todt, D.; Pfaender, S.; Steinmann, E. Persistence of coronaviruses on inanimate surfaces and their inactivation with biocidal agents. J. Hosp. Infect. 2020, 104, 246–251. [Google Scholar] [CrossRef] [Green Version]
  157. Van Doremalen, N.; Bushmaker, T.; Morris, D.H.; Holbrook, M.G.; Gamble, A.; Williamson, B.N.; Tamin, A.; Harcourt, J.L.; Thornburg, N.J.; Gerber, S.I.; et al. Aerosol and Surface Stability of SARS-CoV-2 as Compared with SARS-CoV-1. N. Engl. J. Med. 2020. [Google Scholar] [CrossRef]
  158. Desai, J.; Gan, H.; Barrow, C.; Jameson, M.; Atkinson, V.; Haydon, A.; Millward, M.; Begbie, S.; Brown, M.; Markman, B.; et al. Phase I, Open-Label, Dose-Escalation/Dose-Expansion Study of Lifirafenib (BGB-283), an RAF Family Kinase Inhibitor, in Patients with Solid Tumors. J. Clin. Oncol. 2020. [Google Scholar] [CrossRef]
  159. Drilon, A.; Siena, S.; Ou, S.I.; Patel, M.; Ahn, M.J.; Lee, J.; Bauer, T.M.; Farago, A.F.; Wheler, J.J.; Liu, S.V.; et al. Safety and Antitumor Activity of the Multitargeted Pan-TRK, ROS1, and ALK Inhibitor Entrectinib: Combined Results from Two Phase I Trials (ALKA-372-001 and STARTRK-1). Cancer Discov. 2017, 7, 400–409. [Google Scholar] [CrossRef] [Green Version]
  160. Drilon, A.; Siena, S.; Dziadziuszko, R.; Barlesi, F.; Krebs, M.G.; Shaw, A.T.; de Braud, F.; Rolfo, C.; Ahn, M.J.; Wolf, J.; et al. Entrectinib in ROS1 fusion-positive non-small-cell lung cancer: Integrated analysis of three phase 1-2 trials. Lancet Oncol. 2020, 21, 261–270. [Google Scholar] [CrossRef]
  161. Doebele, R.C.; Drilon, A.; Paz-Ares, L.; Siena, S.; Shaw, A.T.; Farago, A.F.; Blakely, C.M.; Seto, T.; Cho, B.C.; Tosi, D.; et al. Entrectinib in patients with advanced or metastatic NTRK fusion-positive solid tumours: Integrated analysis of three phase 1-2 trials. Lancet Oncol. 2020, 21, 271–282. [Google Scholar] [CrossRef]
  162. Facchinetti, F.; Friboulet, L. Profile of entrectinib and its potential in the treatment of ROS1-positive NSCLC: Evidence to date. Lung Cancer 2019, 10, 87–94. [Google Scholar] [CrossRef] [PubMed] [Green Version]
Figure 1. (A) Therapeutic relevance of targeting the nsp16/nsp10 2′-O-methyltransferase (2′-O-MTase) protein complex. Severe acute respiratory syndrome coronavirus-2 (SARS-CoV-2) encodes two S-adenosyl-l-methionine (SAM)-dependent methyltransferases (MTases), which sequentially methylate the RNA cap at the guanosine-7 and ribose-2′-O-positions, and are catalyzed by nsp14 N7-MTase and nsp16 2′-O-MTase, respectively. The molecular traits of the coronaviral 2′-O-MTase relative to other eukaryotic MTases make the nsp16/nsp10 protein complex an attractive target for inhibitor design. A unique feature of SARS-CoV-2 is that nsp16 obligatorily requires the non-structural protein nsp10 as a stimulatory factor for exerting its 2′-O-MTase activity by stabilizing the co-substrate SAM-binding pocket and extending the substrate RNA-binging groove. Drugs capable of inhibiting the formation of the nsp16/nsp10 complex (by targeting the nsp16/nsp10 activating interface) and the 2′-O-MTase activity (by targeting the SAM-binding site and the RNA-binding groove) can directly contribute to the suppression of genome replication. The pharmacological reduction of nsp16 2′-O-MTase activity can disrupt the 2′-O-methylation in the RNA cap structure and promote the accumulation of cap-0 structures with lower efficiency of translation, thereby suppressing the synthesis of viral proteins, while the viral RNA lacking a cap-1 structure can be recognized by innate RNA sensors and stimulate an early type I interferon response. (B) A computational approach for uncovering SARS-CoV2 2′-O-MTAse-targeting drugs. We performed a structure-based virtual screening (VS) procedure by employing molecular docking of almost 9000 U.S. Food and Drug Administration (FDA)-approved investigational and experimental (discovery-phase) drugs to predict candidates targeting one or many of the three differential traits of the intermolecular interactions of the SARS-CoV-2-specific nsp16/nsp10 heterodimer, namely: a SAM cofactor-binding site slightly different to that of main eukaryotic non-viral MTases; the interface between the stimulatory subunit nsp10 and the catalytic subunit nsp16, which is uniquely different from all other partner-independent 2′-O-MTases; and the (nsp10) allosterically regulated extension and narrowing of the substrate RNA-binding groove. First, we selected the top 20 high-scoring compounds (which therefore have a potentially higher affinity) for each site in a first filtering step; second, based on root mean square deviation (RMSD) values and the molecular mechanics Poisson–Boltzmann surface area (MM/PBSA) parameter, we selected a final list of 20 candidates that were predicted to strongly and stably interact with the abovementioned binding sites of the nsp16/nsp10 complex.
Figure 1. (A) Therapeutic relevance of targeting the nsp16/nsp10 2′-O-methyltransferase (2′-O-MTase) protein complex. Severe acute respiratory syndrome coronavirus-2 (SARS-CoV-2) encodes two S-adenosyl-l-methionine (SAM)-dependent methyltransferases (MTases), which sequentially methylate the RNA cap at the guanosine-7 and ribose-2′-O-positions, and are catalyzed by nsp14 N7-MTase and nsp16 2′-O-MTase, respectively. The molecular traits of the coronaviral 2′-O-MTase relative to other eukaryotic MTases make the nsp16/nsp10 protein complex an attractive target for inhibitor design. A unique feature of SARS-CoV-2 is that nsp16 obligatorily requires the non-structural protein nsp10 as a stimulatory factor for exerting its 2′-O-MTase activity by stabilizing the co-substrate SAM-binding pocket and extending the substrate RNA-binging groove. Drugs capable of inhibiting the formation of the nsp16/nsp10 complex (by targeting the nsp16/nsp10 activating interface) and the 2′-O-MTase activity (by targeting the SAM-binding site and the RNA-binding groove) can directly contribute to the suppression of genome replication. The pharmacological reduction of nsp16 2′-O-MTase activity can disrupt the 2′-O-methylation in the RNA cap structure and promote the accumulation of cap-0 structures with lower efficiency of translation, thereby suppressing the synthesis of viral proteins, while the viral RNA lacking a cap-1 structure can be recognized by innate RNA sensors and stimulate an early type I interferon response. (B) A computational approach for uncovering SARS-CoV2 2′-O-MTAse-targeting drugs. We performed a structure-based virtual screening (VS) procedure by employing molecular docking of almost 9000 U.S. Food and Drug Administration (FDA)-approved investigational and experimental (discovery-phase) drugs to predict candidates targeting one or many of the three differential traits of the intermolecular interactions of the SARS-CoV-2-specific nsp16/nsp10 heterodimer, namely: a SAM cofactor-binding site slightly different to that of main eukaryotic non-viral MTases; the interface between the stimulatory subunit nsp10 and the catalytic subunit nsp16, which is uniquely different from all other partner-independent 2′-O-MTases; and the (nsp10) allosterically regulated extension and narrowing of the substrate RNA-binding groove. First, we selected the top 20 high-scoring compounds (which therefore have a potentially higher affinity) for each site in a first filtering step; second, based on root mean square deviation (RMSD) values and the molecular mechanics Poisson–Boltzmann surface area (MM/PBSA) parameter, we selected a final list of 20 candidates that were predicted to strongly and stably interact with the abovementioned binding sites of the nsp16/nsp10 complex.
Viruses 12 00525 g001
Figure 2. Predicted potential inhibitors of the nsp16/nsp10 2′-O-MTAse complex. Structures of the potential drugs were predicted to target: (a) the (nsp10-stabilized) S-adenosyl-l-methionine (SAM)-binding pocket of nsp16 (green code), (b) the unique “activating surface” between nsp16 and nsp10 (brown code), and (c) the (nsp10-extended) RNA-binding groove of nsp16 (pink code).
Figure 2. Predicted potential inhibitors of the nsp16/nsp10 2′-O-MTAse complex. Structures of the potential drugs were predicted to target: (a) the (nsp10-stabilized) S-adenosyl-l-methionine (SAM)-binding pocket of nsp16 (green code), (b) the unique “activating surface” between nsp16 and nsp10 (brown code), and (c) the (nsp10-extended) RNA-binding groove of nsp16 (pink code).
Viruses 12 00525 g002
Figure 3. Incorporation models of the predicted SAM-binding-site-targeting drugs in SARS-CoV-2 2′-O-MTAse. The best poses of the predicted candidates coupled to the SAM-binding site of the nsp16/nsp10 complex before (0 ns) and after (100 ns) the molecular dynamics (MD) simulation are shown. The protein is represented as a function of the hydrophobicity of its surface amino acids and the Na+ and Cl ions have been eliminated to facilitate visualization. Each inset shows the detailed interactions of each drug candidate docked to the SAM-binding site of nsp16, indicating the participating amino acids involved in the interaction and the type of interaction (hydrogen bonds, hydrophilic interactions, salt bridges, Π-stacking, etc). The root mean square deviation (RMSD, Å) of each drug’s heavy atoms over the simulation time, measured after superposing the protein on its reference structure, is shown. In the context of the MD presented here, the RMSD incorporates traces of SAM in two independent control simulations (i.e., nsp16/nsp10 complex and nsp16 alone).
Figure 3. Incorporation models of the predicted SAM-binding-site-targeting drugs in SARS-CoV-2 2′-O-MTAse. The best poses of the predicted candidates coupled to the SAM-binding site of the nsp16/nsp10 complex before (0 ns) and after (100 ns) the molecular dynamics (MD) simulation are shown. The protein is represented as a function of the hydrophobicity of its surface amino acids and the Na+ and Cl ions have been eliminated to facilitate visualization. Each inset shows the detailed interactions of each drug candidate docked to the SAM-binding site of nsp16, indicating the participating amino acids involved in the interaction and the type of interaction (hydrogen bonds, hydrophilic interactions, salt bridges, Π-stacking, etc). The root mean square deviation (RMSD, Å) of each drug’s heavy atoms over the simulation time, measured after superposing the protein on its reference structure, is shown. In the context of the MD presented here, the RMSD incorporates traces of SAM in two independent control simulations (i.e., nsp16/nsp10 complex and nsp16 alone).
Viruses 12 00525 g003
Figure 4. Incorporation models of the predicted nsp16-activating-surface-targeting drugs in SARS-CoV-2 2′-O-MTAse. The best poses of the predicted candidates coupled to the nsp16 activating surface before (0 ns) and after (100 ns) the molecular dynamics (MD) simulation are shown. The protein has been represented as a function of the hydrophobicity of its surface amino acids and the Na+ and Cl ions have been eliminated to facilitate visualization. Each inset shows the detailed interactions of each drug candidate docked to the nsp16 activating surface, indicating the participating amino acids involved in the interaction and the type of interaction (hydrogen bonds, hydrophilic interactions, salt bridges, Π-stacking, etc). The root mean square deviation (RMSD, Å) of each drug’s heavy atoms over the simulation time, measured after superposing the protein on its reference structure, is shown. In the context of the MD presented here, the RMSD incorporates traces of SAM alone and in two additional control simulations (i.e., nsp16/nsp10 complex and nsp16 alone).
Figure 4. Incorporation models of the predicted nsp16-activating-surface-targeting drugs in SARS-CoV-2 2′-O-MTAse. The best poses of the predicted candidates coupled to the nsp16 activating surface before (0 ns) and after (100 ns) the molecular dynamics (MD) simulation are shown. The protein has been represented as a function of the hydrophobicity of its surface amino acids and the Na+ and Cl ions have been eliminated to facilitate visualization. Each inset shows the detailed interactions of each drug candidate docked to the nsp16 activating surface, indicating the participating amino acids involved in the interaction and the type of interaction (hydrogen bonds, hydrophilic interactions, salt bridges, Π-stacking, etc). The root mean square deviation (RMSD, Å) of each drug’s heavy atoms over the simulation time, measured after superposing the protein on its reference structure, is shown. In the context of the MD presented here, the RMSD incorporates traces of SAM alone and in two additional control simulations (i.e., nsp16/nsp10 complex and nsp16 alone).
Viruses 12 00525 g004
Figure 5. Incorporation models of predicted RNA-binding-groove-targeting drugs in SARS-CoV-2 2′-O-MTAse. The best poses of predicted candidates coupled to the RNA activating surface of nsp16 before (0 ns) and after (100 ns) the molecular dynamics (MD) simulation are shown. The protein has been represented as a function of the hydrophobicity of its surface amino acids and the Na+ and Cl- ions have been eliminated to facilitate visualization. Each inset shows the detailed interactions of each drug candidate docked to the RNA-binding groove of nsp16, indicating the participating amino acids involved in the interaction and the type of interaction (hydrogen bonds, hydrophilic interactions, salt bridges, Π-stacking, etc). The root mean square deviation (RMSD, Å) of each drug’s heavy atoms over the simulation time, measured after superposing the protein on its reference structure, is shown. In the context of the MD presented here, the RMSD incorporates traces of SAM alone and in two additional control simulations (i.e., nsp16/nsp10 complex and nsp16 alone).
Figure 5. Incorporation models of predicted RNA-binding-groove-targeting drugs in SARS-CoV-2 2′-O-MTAse. The best poses of predicted candidates coupled to the RNA activating surface of nsp16 before (0 ns) and after (100 ns) the molecular dynamics (MD) simulation are shown. The protein has been represented as a function of the hydrophobicity of its surface amino acids and the Na+ and Cl- ions have been eliminated to facilitate visualization. Each inset shows the detailed interactions of each drug candidate docked to the RNA-binding groove of nsp16, indicating the participating amino acids involved in the interaction and the type of interaction (hydrogen bonds, hydrophilic interactions, salt bridges, Π-stacking, etc). The root mean square deviation (RMSD, Å) of each drug’s heavy atoms over the simulation time, measured after superposing the protein on its reference structure, is shown. In the context of the MD presented here, the RMSD incorporates traces of SAM alone and in two additional control simulations (i.e., nsp16/nsp10 complex and nsp16 alone).
Viruses 12 00525 g005
Figure 6. Potential repurposed drug candidates against SARS-CoV2 2′-O-MTAse: a multi-layer network framework. Upon the selection of strong candidates based on the stability of high-scoring nsp16/nsp10–drug complexes (Table S2; Figure 3), we re-evaluated their pharmacological overlapping with functional modules of the virus–host interactome relevant for the viral lifecycle, as well as with clinical and laboratory features of COVID-19. The proposed multi-layer network framework might help to select the narrowest list of candidates that can be rapidly tested experimentally before evaluating their in vivo efficiency and side-effects.
Figure 6. Potential repurposed drug candidates against SARS-CoV2 2′-O-MTAse: a multi-layer network framework. Upon the selection of strong candidates based on the stability of high-scoring nsp16/nsp10–drug complexes (Table S2; Figure 3), we re-evaluated their pharmacological overlapping with functional modules of the virus–host interactome relevant for the viral lifecycle, as well as with clinical and laboratory features of COVID-19. The proposed multi-layer network framework might help to select the narrowest list of candidates that can be rapidly tested experimentally before evaluating their in vivo efficiency and side-effects.
Viruses 12 00525 g006

Share and Cite

MDPI and ACS Style

Encinar, J.A.; Menendez, J.A. Potential Drugs Targeting Early Innate Immune Evasion of SARS-Coronavirus 2 via 2’-O-Methylation of Viral RNA. Viruses 2020, 12, 525. https://doi.org/10.3390/v12050525

AMA Style

Encinar JA, Menendez JA. Potential Drugs Targeting Early Innate Immune Evasion of SARS-Coronavirus 2 via 2’-O-Methylation of Viral RNA. Viruses. 2020; 12(5):525. https://doi.org/10.3390/v12050525

Chicago/Turabian Style

Encinar, José Antonio, and Javier A. Menendez. 2020. "Potential Drugs Targeting Early Innate Immune Evasion of SARS-Coronavirus 2 via 2’-O-Methylation of Viral RNA" Viruses 12, no. 5: 525. https://doi.org/10.3390/v12050525

APA Style

Encinar, J. A., & Menendez, J. A. (2020). Potential Drugs Targeting Early Innate Immune Evasion of SARS-Coronavirus 2 via 2’-O-Methylation of Viral RNA. Viruses, 12(5), 525. https://doi.org/10.3390/v12050525

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop