Next Article in Journal
The Structure of an AAV5-AAVR Complex at 2.5 Å Resolution: Implications for Cellular Entry and Immune Neutralization of AAV Gene Therapy Vectors
Next Article in Special Issue
Fluoroquinolone Antibiotics Exhibit Low Antiviral Activity against SARS-CoV-2 and MERS-CoV
Previous Article in Journal
Application of Viral Vectors for Vaccine Development with a Special Emphasis on COVID-19
Previous Article in Special Issue
Drug Repositioning: New Approaches and Future Prospects for Life-Debilitating Diseases and the COVID-19 Pandemic Outbreak
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Target-Centered Drug Repurposing Predictions of Human Angiotensin-Converting Enzyme 2 (ACE2) and Transmembrane Protease Serine Subtype 2 (TMPRSS2) Interacting Approved Drugs for Coronavirus Disease 2019 (COVID-19) Treatment through a Drug-Target Interaction Deep Learning Model

1
Deargen, Inc., Daejeon 34051, Korea
2
Department of Microbiology, College of Natural Sciences, Dankook University, Cheonan 31116, Korea
*
Author to whom correspondence should be addressed.
Submission received: 5 September 2020 / Revised: 5 November 2020 / Accepted: 16 November 2020 / Published: 18 November 2020
(This article belongs to the Special Issue Drug-Repositioning Opportunities for Antiviral Therapy)

Abstract

:
Previously, our group predicted commercially available Food and Drug Administration (FDA) approved drugs that can inhibit each step of the replication of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) using a deep learning-based drug-target interaction model called Molecule Transformer-Drug Target Interaction (MT-DTI). Unfortunately, additional clinically significant treatment options since the approval of remdesivir are scarce. To overcome the current coronavirus disease 2019 (COVID-19) more efficiently, a treatment strategy that controls not only SARS-CoV-2 replication but also the host entry step should be considered. In this study, we used MT-DTI to predict FDA approved drugs that may have strong affinities for the angiotensin-converting enzyme 2 (ACE2) receptor and the transmembrane protease serine 2 (TMPRSS2) which are essential for viral entry to the host cell. Of the 460 drugs with Kd of less than 100 nM for the ACE2 receptor, 17 drugs overlapped with drugs that inhibit the interaction of ACE2 and SARS-CoV-2 spike reported in the NCATS OpenData portal. Among them, enalaprilat, an ACE inhibitor, showed a Kd value of 1.5 nM against the ACE2. Furthermore, three of the top 30 drugs with strong affinity prediction for the TMPRSS2 are anti-hepatitis C virus (HCV) drugs, including ombitasvir, daclatasvir, and paritaprevir. Notably, of the top 30 drugs, AT1R blocker eprosartan and neuropsychiatric drug lisuride showed similar gene expression profiles to potential TMPRSS2 inhibitors. Collectively, we suggest that drugs predicted to have strong inhibitory potencies to ACE2 and TMPRSS2 through the DTI model should be considered as potential drug repurposing candidates for COVID-19.

1. Introduction

Coronavirus disease 2019 (COVID-19) pandemic caused by the severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) [1] has become a serious threat to public health management. At least 55 million cases have been reported in more than 210 countries and territories for about eight months since December 2019, and the COVID-19 outbreak is currently in progress (Worldometer.info, 17 November 2020). While it is urgent to develop vaccines and treatments to prevent the spread of viral infections, unfortunately, no effective drug or therapy has been found for COVID-19 yet despite extensive efforts by researchers around the world except remdesivir (GS-5734). Remdesivir, developed as a treatment for Ebola virus disease, has been approved for emergency use in patients with severe COVID-19 by U.S. Food and Drug Administration (FDA) [2]. More recently, only remdesivir, now commercially named as Veklury, has been approved as the first treatment for COVID-19 by FDA (FDA press release, www.fda.gov, 22 October 2020). Multinational randomized clinical trials show that remdesivir does not make a significant difference in mortality, but is known to shorten the recovery period for patients [3]. However, more effective medications are still necessary for patients suffering from symptoms of SARS-CoV-2 infection.
Our group previously predicted several antiviral (atazanavir, remdesivir, lopinavir/ritonvir) and non-antiviral drugs (rapamycin, tiotropium) that might inhibit SARS-CoV-2 using a new deep drug–target interaction (DTI) prediction model called Molecule Transformer (MT)-DTI [4]. In the previous analysis, these drugs showed strong binding affinity values for viral proteins of SARS-CoV-2. Indeed, several antiviral drugs, such as atazanavir and remdesivir, predicted by MT-DTI, showed potent in vitro and in vivo inhibitory effects on SARS-CoV-2 infection [5,6]. Recently, we further examined that top-ranked non-antiviral drug, tiotropium bromide, has an inhibitory effect on some of the SARS-CoV-2-induced genes through transcriptome-based analysis [7]. While analyzing tiotropium bromide, it was suggested that there is a possibility of more therapeutic targets and options for COVID-19 that may be predicted by deep learning technology to fight against the on-going pandemic. To carry out the task, angiotensin-converting enzyme 2 (ACE2) and transmembrane protease serine subtype 2 (TMPRSS2) that the most studied host targets at the moment, were selected.
SARS-CoV-2 uses the ACE2 receptor and TMPRSS2 for infection into host cells [8,9]. ACE2, which serves as a binding receptor of the spike protein of SARS-CoV-2, involves in blood pressure regulation and cardiovascular regulation through enzymatic hydrolysis of angiotensin (Ang) II to Ang (1–7). Ang II cleavage by ACE2 counteracts the function of ACE that promotes Ang II formation [10,11]. Interestingly, since ACE2 is already known as a binding receptor for SARS-CoV [12], it supports that SARS-CoV-2, which maintains a very close evolutionary relationship with SARS-CoV, also utilizes ACE2 for host cell entry. Another proposed target protein, TMPRSS2, is a cell membrane-bound serine protease [13] which function is yet unclear since its identification, but the involvement of TMPRSS2 during SARS-CoV-2 entry is suggested as the host protease for the viral spike-ACE2 complex proteolytic activation [8,14].
Due to the role of ACE2 and TMPRSS2, which are essential for virus entry, they have been proposed as host targets to block SARS-CoV-2 entry [14]. Therefore, in this study, we predicted drug candidates that can control these two potential targets by using a deep-learning approach.

2. Materials and Methods

2.1. Drug Screening Database Design

Approved small molecule drugs were collected from DrugBank database [15] and listed for the drug screening database of the present study. Then, all stereoisomers of approved drugs were collected from ZINC database [16] and merged to the drug screening database by using simplified molecular-input line-entry system (SMILES) strings as queries. SMILES strings identical in both DrugBank and ZINC database were treated as approved isomers, and drugs with no SMILES match were treated as racemates (e.g., thalidomide, lenalidomide, pomalidomide, etc.). A total of 1400 drugs were selected for the screening process.

2.2. Affinity Prediction of Each Drug to ACE2 and TMPRSS2

A drug-target interaction deep learning model previously described as MT-DTI was used for the affinity prediction [17]. The model is pretrained with approximately 100,000,000 SMILES of chemicals in PubChem database [18] to produce a better chemical representation that considers global contexts of a given SMILES sequence through a self-attention mechanism. This self-attention mechanism proved its effectiveness in many natural language process tasks compared to local context-based models. With the pretrained MT chemical representation and a convolutional neural network (CNN) protein amino sequence representation, the MT-DTI model was trained with a combined and curated chemical-protein pairs from the Drug Target Common database [19] and BindingDB database [20] to predict affinity in Kd value of given chemical-protein pairs. The performance evaluation (Concordance index, 0.882 and 0.887 in KIBA and DAVIS datasets; mean squared error, 0.152 and 0.245 in KIBA and DAVIS dataset) and comparison to other models of MT-DTI is reported in Shin et al. (2019). Finally, the screening database was processed to the model with amino acid sequences of human ACE2 (UniProt entry Q9BYF1, NCBI Reference Sequence: NP_001358344.1) and TMPRSS2 (UniProt entry O15393, NCBI Reference Sequence: NP_001128571.1) as targets.
After the prediction, results were screened for drugs with affinity Kd < 100 nM to each target protein. Drugs approved as racemic mixture were excluded as they cannot be expressed in one SMILES string and may have differential ratio of isomers and biological outcomes per each isomer in the human. DTI prediction results without any filters are provided as supplementary materials.

2.3. Cross-Prediction of MT-DTI Results through AutoDock Vina Docking Tool

AutoDock Vina (version 1.1.2), which is a molecular docking and virtual screening application [21], was used to predict binding affinities (kcal/mol) between target proteins (ACE2 [PDB ID: 1R42] or TMPRSS2) and given chemical compounds. SMILES of given chemical compounds were converted to the PDBQT format using Open Babel (version 2.3.2) [22] with the following options: --gen3d and -p 7.4. A homology model of TMPRSS2 was built using the SWISS-MODEL server (https://swissmodel.expasy.org/) due to the unavailability of the 3D structure in the protein data bank (PDB) [23]. AutoDockTools was used to convert coordinates of the given target proteins into the form needed for docking calculations by adding charges and hydrogen atoms, and merging nonpolar hydrogens (version 1.5.6) [24]. Then, binding affinities between the target proteins and given chemical compounds were calculated using AutoDock Vina. The following parameters were set for the search space: —center_x 90.654 —center_y 66.551 —center_z 39.925 —size_x 38 —size_y 92 —size_z 80 for ACE2 and —center_x 10.250 —center_y -1.306 —center_z 13.472 —size_x 28 —size_y 28 —size_z 28 for TMPRSS2. The exhaustiveness parameter was set to 10.

2.4. Connectivity Map Analysis

The molecular signatures of protease inhibitors proposed as TMPRSS2 inhibitors, including bromhexine, probucol, were used to compare the gene expression signatures of the top 30 drugs in MT-DTI results for TMPRSS2 using the Touchstone database of the CLUE platform (https://clue.io) [25]. Each compound was characterized according to the connectivity score ranging from −100 to 100. Positive scores indicate the similarity between the signatures of the two compounds, and negative scores indicate the difference between the two signatures. The score corresponds to the degree of similarity or dissimilarity.

3. Results

To identify potent FDA-approved drugs capable of inhibiting the entry of SARS-CoV-2, we performed an in silico screening method using the MT-DTI deep learning-based model [4]. The MT-DTI model can predict the binding affinities of Kd values from on chemical sequences (SMILES) and amino acid sequences (FASTA) of a target protein. The sequence-based drug-target affinity prediction approach predicted drugs with a strong binding affinity (Kd < 100 nM) against ACE2 and TMPRSS2 from 1400 FDA-approved drugs (ACE2, Table 1 and Table S1; TMPRSS2, Table 1 and Table S2). The NCATS OpenData portal, which was developed to provide insight into drug repurposing for the treatment of COVID-19 (https://opendata.ncats.nih.gov/covid19) [26], provides drug results that affect protein–protein interaction between ACE2 and RBD (Receptor-binding domain) of the SARS-CoV-2 S protein through AlphaLISA and TruHit (counterscreen) ACE2-Spike protein–protein interaction proximity assay. We first compared the results obtained from the MT-DTI model with the drugs provided by the OpenData portal. Of the 460 drugs that may bind ACE2 with strong affinity, 17 drugs showed the interference of the interaction of ACE2 with SARS-CoV-2 spike protein (Table 2). Although AC50 values of the proximity assay are not affinity values, submicromolar to one-digit micromolar AC50 results supports the ACE2 affinity prediction of 17 selected drugs. Interestingly, enalaprilat, an ACE inhibitor, was identified as a chemical compound that has a strong affinity for ACE2 and is likely to inhibit the entry of SARS-CoV-2. According to previous results, the angiotensin-converting enzyme (ACE) and its homologous ACE2 both belong to the ACE family of dipeptidyl carboxypeptidase but provide the opposite effect for the renin-angiotensin system. ACE produces Ang II, which induces enhanced inflammation, elevated blood pressure, and increased coagulation. However, ACE2 is responsible for the breakdown and inactivation of Ang II [10,11]. Based on our predictions, enalaprilat suggests that it is a promising drug candidate that can reduce inflammation and blood pressure by inhibiting ACE activity and simultaneously inhibit the entry of SARS-Cov-2 through interaction with ACE2. Interestingly, venetoclax (Kd 6.12 nM), posaconazole (Kd 17.11 nM), daclatasvir (Kd 6.65 nM), and ombitasvir (Kd 5.91 nM) also predicted to have strong affinity for TMPRSS2, as shown in Table S2. In particular, daclatasvir and ombitasvir are hepatitis C virus (HCV) inhibitors, suggesting that these compounds may exhibit antiviral effects on HCV and SARS-CoV-2.
Another strategy to block the entry of SARS-CoV-2 by targeting ACE2 is to increase the human recombinant soluble ACE2 (hrsACE2), thereby inhibiting the membrane association of host cells with SARS-CoV-2 [27]. Indeed, hrsACE2 has already undergone phase 1 and phase 2 clinical trials to treat acute respiratory distress syndrome and is now considered a candidate for treatment for COVID-19. In the cell, the soluble form of ACE2 is endogenously generated through ACE2 ectodomain shedding by the disintegrin metalloproteinase 17 (ADAM-17) [28]. On the other hand, calmodulin inhibits ectodomain shedding by interacting with ACE2 [29]. Therefore, inhibition of calmodulin-ACE2 interaction through the calmodulin antagonists may be an effective strategy to reduce the infectivity of the virus. We predicted that the calmodulin antagonist, the oxide of the amitriptyline, binds to ACE2 with a Kd of 40.38 nM (Table S1). However, it requires experimental verification that the drug’s interaction with ACE2 further promotes ACE2 shedding.
To compare MT-DTI affinity prediction results to widely used docking study tool AutoDock Vina, Gibbs free energy (∆G) values of each drug-target pair were predicted. For the ACE2 cross-prediction result, ∆G −8.0 kcal/mol is given as threshold, and 27 drugs show correlative prediction results out of 460 DTI predicted results (Table 3, left). The TMPRSS2 cross-prediction result was screened with ∆G −7.0 kcal/mol threshold as the numbers of predicted drugs are less than 100, and 11 out of 75 drugs showed correlation to the MT-DTI prediction (Table 3, right). However, TMPRSS2 docking results through a homology modeled TMPRSS2 may not be precise as ACE2 due to the absence of an actual 3D structure, therefore, it should be considered again when crystallography results are available. All ∆G values predicted through AutoDock Vina are annotated in Tables S1 and S2.
We have reported in the previous MT-DTI studies that tiotropium bromide has a strong binding affinity to the SARS-CoV-2 viral proteins [4]. A recent RNA-seq analysis of normal bronchial epithelium (NHBE) cells infected with SARS-CoV-2 has shown that the treatment of tiotropium can have an advantageous prescription effect [7]. Interestingly, tiotropium bromide was also predicted to strongly bind with ACE2 (Kd 0.92 nM) in the present study, suggesting that the drug is a promising candidate for COVID-19.
We further investigated 75 drugs predicted to have an affinity of less than 100 nM for TMPRSS2 that is a serine protease for S protein priming of SARS-CoV-2. Clinically proven protease inhibitors, including bromhexine, aprotinin, camostat, and nafamostat, have been suggested as potential treatment options for COVID-19 [30,31]. However, the MT-DTI results for TMPRSS2 showed affinity of more than 100 nM for these protease inhibitors. Intriguingly, we found that five of the top 30 drugs that strongly bind to TMPRSS2 are antiviral drugs with a Kd value of less than 20 nM (Table S2). Three out of five (ombitasvir, daclatasvir, and paritaprevir) are anti-hepatitis C virus (HCV) drugs, and recently, there are reports that antiviral agents used to treat HCV could be another option for SARS-CoV-2 treatment [32]. Although HCV and coronavirus are not closely related viruses, it can be suggested that those drugs may act on life cycles of both viruses in a mode of action perspective as they are single-stranded RNA viruses. These three anti-HCV drugs predicted to have Kd values of less than 10 nM for TMPRSS2, implicating a strong inhibitory potency (Table 1).
To screen for more effective drugs that can inhibit TMPRSS2, we used the connectivity map (CMap) database [33] that allows gene-expression signature matching between the drugs proposed as TMPRSS2 inhibitors and the top 30 drugs of the MT-DTI result. CMap is a database that provides transcriptome profiling information caused through genetic and chemical perturbations [33]. We compared gene expression changes by the previously proposed TMPRSS2 inhibitors (Bromhexine [34] and Probucol [35]) and the top 30 drugs of MT-DTI through CMap Linked User Environment (CLUE, https://clue.io) (Table 4). Scores above +90 indicate a high positive correlation between the gene-expression signature of a given perturbagen and the gene-expression signature by the query. Among the top 30 drugs, the gene expression by eprosartan and lisuride was strongly correlated with the gene expression signatures of bromhexine and probucol, respectively. Interestingly, eprosartan is a type I angiotensin II receptor (AT1R) blocker, preventing Ang II from binding to the AT1R receptor in the renin-angiotensin system [36]. As the binding of Ang II and AT1R causes fibrosis and damage of lung tissues, AT1R blocker, eprosartan, is considered a suitable therapeutic target for COVID-19 [35,36]. In addition to eprosartan, neuropsychiatric agents, lisuride, have recently been predicted as a candidate drug for treating COVID-19 through coexpression-based drug enrichment analysis of COVID-19 induced genes [37].

4. Discussion

An AI-integrated drug repurposing approach, such as MT-DTI, is a promising solution for the development of potential drugs that can overcome life-threatening diseases caused by COVID-19. For example, atazanavir, predicted by MT-DTI in the previous study [4], is registered in the REVOLUTIOn trial (NCT04468087) for the treatment of COVID-19. As a case study for ACE2 affinity prediction, a recent report on a clinical trial RECOVERY (NCT04381936) suggested dexamethasone, an anti-inflammatory corticosteroid, as a possible treatment option for COVID-19 [38,39]. Interestingly, dexamethasone was predicted to have an affinity of Kd 9.50 nM to ACE2 (Table S1). In addition, approved drugs sharing a similar polycyclic core structure of steroids such as obeticholic acid (Kd 0.98 nM), mestranol (Kd 1.02 nM), norethynodrel (Kd 1.51 nM), clobetasol (Kd 1.80 nM), norethisterone (Kd 1.92 nM), fluoxymesterone (Kd 2.29 nM), cholic acid (Kd 2.86 nM), and more were listed on top in the ACE2 affinity prediction. However, original indications and side effects must be considered since steroids are involved in hormonal regulations.
Furthermore, the MT-DTI prediction resulted that ACE inhibitors, including enalaprilat, zofenopril, lisinopril, benazepril, trandolapril, cilazapril, perindopril, ramipril, fosinopril, moexipril, spirapril, have a strong binding affinity of less than Kd 70 nM for ACE2 (Table S1). Although there is concern that inhibition of ACE2 activity may promote an excessive inflammatory response resulting from an increase in Ang II [40], ACE inhibitors can be considered as a treatment option for COVID-19 in two respects. First of all, inhibition of ACE function that elicits an immune response via Ang II may alleviate the side effects of loss of ACE2 function by ACE inhibitors. Another is that ACE inhibitors do not inhibit the activity of ACE2, but rather interfere with the interaction between ACE2 and the viral S protein. Indeed, ACE2 is not only known to be insensitive to ACE inhibitors [41], but the NCATS OpenData portal provides experimental results that the enalapril maleate disrupts the protein–protein interaction between ACE2 and S protein with an AC50 of 7.5 µM. Therefore, our results suggest that the ACE inhibitor is a promising therapeutic option.
We further identified three anti-HCV drugs, including ombitasvir, daclatasvir, and paritaprevir, in the top 30 predictive candidates for TMPRSS2. Another group recently reported that anti-HCV reagents were potential drug repurposing candidates for COVID-19 [42], supporting the results. Above all, we were able to present more effective drug candidates by merging AI-based prediction results and CMap-based drug repurposing approach in this study. Eprosartan and lisuride, drugs predicted to have overall strong binding affinity values for TMPRSS2, demonstrated that they could effectively control TMPRSS2 by showing gene expression signatures similar to the previously proposed TMPRSS2 inhibitors (bromhexine, a mucolytic drug for respiratory diseases; probucol, an anti-hyperlipidemic drug). Bromhexine and probucol are known to block TMPRSS2, but both drugs are currently used as treatments for different indications and have different mode of actions other than TMPRSS2 inhibition. Consistent with this, the connectivity score between the two previously proposed TMPRSS2 inhibitors is not high. Therefore, the drugs predicted by MT-DTI do not seem to show a strong correlation for both drugs simultaneously. Interestingly, the dasatinib, which was predicted to have the strongest affinity for TMPRSS2, had little correlation with bromhexine and probucol implicating that chemical and transcriptional influence of dasatinib may be different compared to both formerly proposed TMPRSS2 inhibitors. It suggests that CMap analysis is a useful approach for filtering out possible false-positive candidates that have opposite effect or little biological effect from DTI prediction results. However, at least, any positive connectivity score, better with higher positive connectivity score, should be considered and may guide researchers to pick drugs with priority for drug repurposing.
Finally, to go deeper to target-centered drug repurposing, drugs appearing both in ACE2 and TMPRSS2 such as aclidinium, buprenorphine, emtricitabine, lurasidone, and tiotropium may be investigated as dual targeting drugs in polypharmacological perspective [43,44]. Therefore, we hope that our results can contribute to the therapeutic management of COVID-19 through evaluation efficacy and safety in vitro and in vivo studies.
Taken together, the drugs resulted through the MT-DTI affinity prediction for ACE2 and TMPRSS2 are suggested promising drug candidates for an effective drug repurposing strategy to treat COVID-19 in a target-wise perspective.

Supplementary Materials

The following are available online at https://0-www-mdpi-com.brum.beds.ac.uk/1999-4915/12/11/1325/s1, Table S1: Predicted candidates identified by MT-DTI against ACE2 (NCBI Reference Sequence: NP_001358344.1). Table S2: Predicted candidates identified by MT-DTI against TMPRSS2 (NCBI Reference Sequence: NP_001128571.1).

Author Contributions

Conceptualization, Y.C. and B.R.B.; Data curation, B.S., S.P. and B.R.B.; Methodology, B.S., K.K., S.P. and B.R.B.; Software, B.S. and S.P.; Writing—original draft, Y.C. and B.R.B.; Writing—review and editing, K.K. and B.R.B. All authors have read and agreed to the published version of the manuscript.

Funding

The present study was supported by the National Research Foundation of Korea (NRF, grant number NRF-2020M3A9I2081698) funded by the Ministry of Science and ICT (MIST).

Conflicts of Interest

Y.C., B.S., S.P. and B.R.B. are employed by Deargen Inc., and K.K. is one of the cofounders of, and a shareholder in, Deargen Inc.

References

  1. Zhu, N.; Zhang, D.; Wang, W.; Li, X.; Yang, B.; Song, J.; Zhao, X.; Huang, B.; Shi, W.; Lu, R.; et al. A novel coronavirus from patients with pneumonia in China, 2019. N. Engl. J. Med. 2020. [Google Scholar] [CrossRef]
  2. Siegel, D.; Hui, H.C.; Doerffler, E.; Clarke, M.O.; Chun, K.; Zhang, L.; Neville, S.; Carra, E.; Lew, W.; Ross, B.; et al. Discovery and Synthesis of a Phosphoramidate Prodrug of a Pyrrolo[2,1-f][triazin-4-amino] Adenine C-Nucleoside (GS-5734) for the Treatment of Ebola and Emerging Viruses. J. Med. Chem. 2017, 60, 1648–1661. [Google Scholar] [CrossRef] [Green Version]
  3. Singh, A.K.; Singh, A.; Singh, R.; Misra, A. Remdesivir in COVID-19: A critical review of pharmacology, pre-clinical and clinical studies. Diabetes Metab. Syndr. Clin. Res. Rev. 2020, 14, 641–648. [Google Scholar] [CrossRef]
  4. Beck, B.R.; Shin, B.; Choi, Y.; Park, S.; Kang, K. Predicting commercially available antiviral drugs that may act on the novel coronavirus (SARS-CoV-2) through a drug-target interaction deep learning model. Comput. Struct. Biotechnol. J. 2020, 18, 784–790. [Google Scholar] [CrossRef] [PubMed]
  5. Choy, K.T.; Wong, A.Y.L.; Kaewpreedee, P.; Sia, S.F.; Chen, D.; Hui, K.P.Y.; Chu, D.K.W.; Chan, M.C.W.; Cheung, P.P.H.; Huang, X.; et al. Remdesivir, lopinavir, emetine, and homoharringtonine inhibit SARS-CoV-2 replication in vitro. Antiviral Res. 2020, 178, 104786. [Google Scholar] [CrossRef] [PubMed]
  6. Fintelman-Rodrigues, N.; Sacramento, C.Q.; Ribeiro Lima, C.; Souza da Silva, F.; Ferreira, A.C.; Mattos, M.; de Freitas, C.S.; Cardoso Soares, V.; da Silva Gomes Dias, S.; Temerozo, J.R.; et al. Atazanavir, alone or in combination with ritonavir, inhibits SARS-CoV-2 replication and pro-inflammatory cytokine production. Antimicrob. Agents Chemother. 2020, 64. [Google Scholar] [CrossRef] [PubMed]
  7. Kang, K.; Kim, H.H.; Choi, Y. Tiotropium is predicted to be a promising drug for COVID-19 through transcriptome-based comprehensive molecular pathway analysis. Viruses 2020, 12, 776. [Google Scholar] [CrossRef] [PubMed]
  8. 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, 271–280.e8. [Google Scholar] [CrossRef] [PubMed]
  9. Xiao, L.; Sakagami, H.; Miwa, N. ACE2: The key molecule for understanding the pathophysiology of severe and critical conditions of COVID-19: Demon or angel? Viruses 2020, 12, 491. [Google Scholar] [CrossRef]
  10. Donoghue, M.; Hsieh, F.; Baronas, E.; Godbout, K.; Gosselin, M.; Stagliano, N.; Donovan, M.; Woolf, B.; Robison, K.; Jeyaseelan, R.; et al. A novel angiotensin-converting enzyme-related carboxypeptidase (ACE2) converts angiotensin I to angiotensin 1-9. Circ. Res. 2000, 87, e1–e9. [Google Scholar] [CrossRef]
  11. Keidar, S.; Kaplan, M.; Gamliel-Lazarovich, A. ACE2 of the heart: From angiotensin I to angiotensin (1-7). Cardiovasc. Res. 2007, 73, 463–469. [Google Scholar] [CrossRef] [PubMed]
  12. Hamming, I.; Timens, W.; Bulthuis, M.L.C.; Lely, A.T.; Navis, G.J.; van Goor, H. Tissue distribution of ACE2 protein, the functional receptor for SARS coronavirus. A first step in understanding SARS pathogenesis. J. Pathol. 2004, 203, 631–637. [Google Scholar] [CrossRef] [PubMed]
  13. Paoloni-Giacobino, A.; Chen, H.; Peitsch, M.C.; Rossier, C.; Antonarakis, S.E. Cloning of the TMPRSS2 gene, which encodes a novel serine protease with transmembrane, LDLRA, and SRCR domains and maps to 21q22.3. Genomics 1997, 44, 309–320. [Google Scholar] [CrossRef] [PubMed]
  14. Tay, M.Z.; Poh, C.M.; Rénia, L.; MacAry, P.A.; Ng, L.F.P. The trinity of COVID-19: Immunity, inflammation and intervention. Nat. Rev. Immunol. 2020, 20, 363–374. [Google Scholar] [CrossRef]
  15. 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]
  16. Sterling, T.; Irwin, J.J. ZINC 15 - Ligand Discovery for Everyone. J. Chem. Inf. Model. 2015, 55, 2324–2337. [Google Scholar] [CrossRef]
  17. Shin, B.; Park, S.; Kang, K.; Ho, J.C. Self-Attention Based Molecule Representation for Predicting Drug-Target Interaction. arXiv 2019, arXiv:1908.06760. [Google Scholar]
  18. Kim, S.; Chen, J.; Cheng, T.; Gindulyte, A.; He, J.; He, S.; Li, Q.; Shoemaker, B.A.; Thiessen, P.A.; Yu, B.; et al. PubChem 2019 update: Improved access to chemical data. Nucleic Acids Res. 2019, 47, D1102–D1109. [Google Scholar] [CrossRef] [Green Version]
  19. Tanoli, Z.R.; Alam, Z.; Vähä-Koskela, M.; Ravikumar, B.; Malyutina, A.; Jaiswal, A.; Tang, J.; Wennerberg, K.; Aittokallio, T. Drug Target Commons 2.0: A community platform for systematic analysis of drug-target interaction profiles. Database 2018, 2018, bay083. [Google Scholar] [CrossRef] [Green Version]
  20. Liu, T.; Lin, Y.; Wen, X.; Jorissen, R.N.; Gilson, M.K. BindingDB: A web-accessible database of experimentally determined protein-ligand binding affinities. Nucleic Acids Res. 2007, 35, D198–D201. [Google Scholar] [CrossRef] [Green Version]
  21. Oleg, T.; Arthur, J.O. 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]
  22. O’Boyle, N.M.; Banck, M.; James, C.A.; Morley, C.; Vandermeersch, T.; Hutchison, G.R. Open Babel: An Open chemical toolbox. J. Cheminform. 2011, 3, 33. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  23. Biasini, M.; Bienert, S.; Waterhouse, A.; Arnold, K.; Studer, G.; Schmidt, T.; Kiefer, F.; Cassarino, T.G.; Bertoni, M.; Bordoli, L.; et al. SWISS-MODEL: Modelling protein tertiary and quaternary structure using evolutionary information. Nucleic Acids Res. 2014, 42, W252–W258. [Google Scholar] [CrossRef] [PubMed]
  24. Morris, G.; Huey, R.; Linkstrom, W.; Sanner, M.; Belew, R.; Goodsell, D. Olson AutoDock4 and AutoDockTools4: Automated Docking with Selective Receptor Flexibility. J. Comput. Chem. 2010, 30, 2785–2791. [Google Scholar] [CrossRef] [Green Version]
  25. Duan, Q.; Flynn, C.; Niepel, M.; Hafner, M.; Muhlich, J.L.; Fernandez, N.F.; Rouillard, A.D.; Tan, C.M.; Chen, E.Y.; Golub, T.R.; et al. LINCS Canvas Browser: Interactive web app to query, browse and interrogate LINCS L1000 gene expression signatures. Nucleic Acids Res. 2014, 42, W449–W460. [Google Scholar] [CrossRef]
  26. Brimacombe, K.R.; Zhao, T.; Eastman, R.T.; Hu, X.; Wang, K.; Backus, M.; Baljinnyam, B.; Chen, C.Z.; Chen, L.; Eicher, T.; et al. An OpenData portal to share COVID-19 drug repurposing data in real time. bioRxiv Prepr. Serv. Biol. 2020. [Google Scholar] [CrossRef]
  27. Monteil, V.; Kwon, H.; Prado, P.; Hagelkrüys, A.; Wimmer, R.A.; Stahl, M.; Leopoldi, A.; Garreta, E.; Hurtado del Pozo, C.; Prosper, F.; et al. Inhibition of SARS-CoV-2 Infections in Engineered Human Tissues Using Clinical-Grade Soluble Human ACE2. Cell 2020, 181, 905–913.e7. [Google Scholar] [CrossRef]
  28. Lambert, D.W.; Yarski, M.; Warner, F.J.; Thornhill, P.; Parkin, E.T.; Smith, A.I.; Hooper, N.M.; Turner, A.J. Tumor necrosis factor-α convertase (ADAM17) mediates regulated ectodomain shedding of the severe-acute respiratory syndrome-coronavirus (SARS-CoV) receptor, angiotensin-converting enzyme-2 (ACE2). J. Biol. Chem. 2005, 280, 30113–30119. [Google Scholar] [CrossRef] [Green Version]
  29. Lambert, D.W.; Clarke, N.E.; Hooper, N.M.; Turner, A.J. Calmodulin interacts with angiotensin-converting enzyme-2 (ACE2) and inhibits shedding of its ectodomain. FEBS Lett. 2008, 582, 385–390. [Google Scholar] [CrossRef]
  30. Sonawane, K.; Barale, S.; Dhanavade, M.; Waghmare, S.; Nadaf, N.; Kamble, S.; Mohammed, A.A.; Makandar, A.; Fandilolu, P.; Dound, A.; et al. Homology Modeling and Docking Studies of TMPRSS2 with Experimentally Known Inhibitors Camostat Mesylate, Nafamostat and Bromhexine Hydrochloride to Control SARS-Coronavirus-2. ChemRxiv 2020. [Google Scholar] [CrossRef]
  31. Uno, Y. Camostat mesilate therapy for COVID-19. Intern. Emerg. Med. 2020, 15, 1577–1578. [Google Scholar] [CrossRef] [PubMed]
  32. Elfiky, A.A. Anti-HCV, nucleotide inhibitors, repurposing against COVID-19. Life Sci. 2020, 248, 117477. [Google Scholar] [CrossRef] [PubMed]
  33. Qu, X.A.; Rajpal, D.K. Applications of Connectivity Map in drug discovery and development. Drug Discov. Today 2012, 17, 1289–1298. [Google Scholar] [CrossRef] [PubMed]
  34. Maggio, R.; Corsini, G.U. Repurposing the mucolytic cough suppressant and TMPRSS2 protease inhibitor bromhexine for the prevention and management of SARS-CoV-2 infection. Pharmacol. Res. 2020, 157, 104837. [Google Scholar] [CrossRef] [PubMed]
  35. Balmeh, N.; Mahmoudi, S.; Mohammadi, N.; Karabedianhajiabadi, A. Predicted therapeutic targets for COVID-19 disease by inhibiting SARS-CoV-2 and its related receptors. Inform. Med. Unlocked 2020, 20, 100407. [Google Scholar] [CrossRef]
  36. Sun, M.L.; Yang, J.M.; Sun, Y.P.; Su, G.H. Inhibitors of RAS Might Be a Good Choice for the Therapy of COVID-19 Pneumonia. Zhonghua Jie He He Hu Xi Za Zhi 2020, 43, E014. [Google Scholar] [CrossRef]
  37. Jia, Z.; Song, X.; Shi, J.; Wang, W.; He, K. Transcriptome-based drug repositioning for coronavirus disease 2019 (COVID-19). Pathog. Dis. 2020, 78, 1–7. [Google Scholar] [CrossRef]
  38. Horby, P.; Lim, W.S.; Emberson, J.R.; Mafham, M.; Bell, J.L.; Linsell, L.; Staplin, N.; Brightling, C.; Ustianowski, A.; Elmahi, E.; et al. Dexamethasone in Hospitalized Patients with Covid-19—Preliminary Report. N. Engl. J. Med. 2020. [Google Scholar] [CrossRef]
  39. Ledford, H. Coronavirus breakthrough: Dexamethasone is first drug shown to save lives. Nature 2020, 582, 469. [Google Scholar] [CrossRef]
  40. Du, F.; Liu, B.; Zhang, S. COVID-19: The role of excessive cytokine release and potential ACE2 down-regulation in promoting hypercoagulable state associated with severe illness. J. Thromb. Thrombolysis 2020. [Google Scholar] [CrossRef]
  41. Turner, A.J.; Tipnis, S.R.; Guy, J.L.; Rice, G.I.; Hooper, N.M. ACEH/ACE2 is a novel mammalian metallocarboxypeptidase and a homologue of angiotensin-converting enzyme insensitive to ACE inhibitors. Can. J. Physiol. Pharmacol. 2002, 80, 346–353. [Google Scholar] [CrossRef] [PubMed]
  42. Hosseini, F.S.; Amanlou, M. Anti-HCV and anti-malaria agent, potential candidates to repurpose for coronavirus infection: Virtual screening, molecular docking, and molecular dynamics simulation study. Life Sci. 2020, 258, 118205. [Google Scholar] [CrossRef] [PubMed]
  43. Anighoro, A.; Bajorath, J.; Rastelli, G. Polypharmacology: Challenges and opportunities in drug discovery. J. Med. Chem. 2014, 57, 7874–7887. [Google Scholar] [CrossRef] [PubMed]
  44. Zhang, W.; Bai, Y.; Wang, Y.; Xiao, W. Polypharmacology in Drug Discovery: A Review from Systems Pharmacology Perspective. Curr. Pharm. Des. 2016, 22, 3171–3181. [Google Scholar] [CrossRef] [PubMed] [Green Version]
Table 1. Top 20 predicted drugs identified by MT-DTI against ACE2 (NCBI Reference Sequence: NP_001358344.1) and TMPRSS2 (NCBI Reference Sequence: NP_001128571.1).
Table 1. Top 20 predicted drugs identified by MT-DTI against ACE2 (NCBI Reference Sequence: NP_001358344.1) and TMPRSS2 (NCBI Reference Sequence: NP_001128571.1).
ACE2TMPRSS2
Small MoleculesPredicted Kd in nMSmall MoleculesPredicted Kd in nM
Pentostatin0.02Dasatinib0.37
Liothyronine0.43Pentostatin0.47
Emtricitabine0.45Tazemetostat2.88
Tiotropium0.92Tiotropium4.03
Obeticholic acid0.98Eluxadoline5.08
Mestranol1.02Pimecrolimus5.79
Gemcitabine1.10Tacrolimus5.82
Levorphanol1.30Ombitasvir5.91
Levallorphan1.37Brexpiprazole6.06
Methscopolamine1.40Venetoclax6.12
Enalaprilat1.46Daclatasvir6.66
Bremelanotide1.49Tolvaptan7.33
Norethynodrel1.51Aclidinium8.46
Epicriptine1.56Paritaprevir8.75
Rolitetracycline1.57Eprosartan9.19
Diacetyl benzoyl lathyrol1.77Cobimetinib10.59
Pizotifen1.80Entrectinib11.08
Dapagliflozin1.80Lisuride11.20
Clobetasol1.80Erdafitinib11.51
Lurasidone1.86Letermovir12.43
Table 2. The list of drugs those may affect protein–protein interaction between ACE2 and SARS-CoV-2 S protein with ACE2 Kd < 100 nM and comparison to AlphaLISA and TruHit proximity assay results provided by NCATS OpenData COVID-19.
Table 2. The list of drugs those may affect protein–protein interaction between ACE2 and SARS-CoV-2 S protein with ACE2 Kd < 100 nM and comparison to AlphaLISA and TruHit proximity assay results provided by NCATS OpenData COVID-19.
Small MoleculesPredicted ACE2 Kd (nM)MT-DTI Rank Out of 460Alphalisa-AC50 (μM)TruHit-AC50 (μM)
Enalaprilat1.46117.52-
Daclatasvir5.81795.97-8.004.41–6.22
Mifepristone10.31245.972.78
Estradiol cypionate11.791369.475.55
Estramustine phosphate19.341966.708.79
Ombitasvir20.912084.50–5.972.78–3.50
Norgestimate21.062109.474.41
d-alpha-Tocopherol acetate21.1821210.622.21
Mitoxantrone25.852430.40–0.600.25–0.31
Pasireotide28.292573.672.69
Ergotamine33.642812.670.44
Flupentixol38.263038.444.94
Venetoclax44.773248.443.50
Anthralin61.923802.386.22
Ciclopirox88.254392.671.56
Thiethylperazine88.984405.66–6.703.12
Posaconazole94.74492.01–6.702.78–3.50
Table 3. Cross-prediction results of ACE2 and TMPRSS2 interacting drugs through MT-DTI and AutoDock Vina.
Table 3. Cross-prediction results of ACE2 and TMPRSS2 interacting drugs through MT-DTI and AutoDock Vina.
ACE2TMPRSS2
Small MoleculesPredicted Kd (nM)AutoDock Vina ∆G (kcal/mol)Small MoleculesPredicted Kd (nM)AutoDock Vina ∆G (kcal/mol)
Bremelanotide1.49−8.1Tazemetostat2.88−7.1
Talazoparib7.04−8Eluxadoline5.08−7.5
Avapritinib14.68−8Entrectinib11.08−7.9
Dihydroergocristine14.94−9.2Erdafitinib11.51−7.1
Tezacaftor15.90−8Aprepitant13.52−7.5
Dutasteride16.35−8.8Canagliflozin27.12−7
Rifapentine17.81−8Naldemedine46.74−7.2
Acetyldigitoxin18.93−8.6Adapalene52.49−7.5
Alatrofloxacin24.46−8Droperidol68.38−7.7
Deslanoside24.75−9.1Larotrectinib71.78−7.4
Dihydroergocornine25.99−8.4Zanubrutinib75.18−7.8
Irinotecan28.37−8.2
Naldemedine29.94−8.8
Ciclesonide31.43−8.1
Ubrogepant33.50−8.1
Ergotamine33.64−9.1
Lumacaftor36.26−8
Venetoclax44.77−9.2
Adapalene45.11−8
Letermovir50.73−8
Paritaprevir53.21−9.2
Entrectinib54.40−8.5
Glycyrrhizic acid61.94−8.1
Simeprevir67.03−8.7
Glecaprevir75.81−8.5
Lifitegrast94.67−8.3
Posaconazole94.70−8
Table 4. Correlation results of seven drugs with similar gene-expression patterns to TMPRSS2 inhibitors (bromhexine and probucol) among top 30 predicted drugs with TMPRSS2 Kd < 100 nM. There was no gene-expression pattern available for those drugs not listed. Bold numbers indicate connectivity score higher than 90 in the −100–100 correlation scale.
Table 4. Correlation results of seven drugs with similar gene-expression patterns to TMPRSS2 inhibitors (bromhexine and probucol) among top 30 predicted drugs with TMPRSS2 Kd < 100 nM. There was no gene-expression pattern available for those drugs not listed. Bold numbers indicate connectivity score higher than 90 in the −100–100 correlation scale.
Small MoleculesPredicted TMPRSS2 Kd (nM)Bromhexine Connectivity ScoreProbucol Connectivity Score
Dasatinib0.3717.5815.38
Tacrolimus5.8210.4070.67
Eprosartan9.1991.1351.17
Lisuride11.2034.7092.18
Aprepitant13.5254.6912.72
Panobinostat15.3133.3041.78
Bosutinib15.950.445.95
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Choi, Y.; Shin, B.; Kang, K.; Park, S.; Beck, B.R. Target-Centered Drug Repurposing Predictions of Human Angiotensin-Converting Enzyme 2 (ACE2) and Transmembrane Protease Serine Subtype 2 (TMPRSS2) Interacting Approved Drugs for Coronavirus Disease 2019 (COVID-19) Treatment through a Drug-Target Interaction Deep Learning Model. Viruses 2020, 12, 1325. https://0-doi-org.brum.beds.ac.uk/10.3390/v12111325

AMA Style

Choi Y, Shin B, Kang K, Park S, Beck BR. Target-Centered Drug Repurposing Predictions of Human Angiotensin-Converting Enzyme 2 (ACE2) and Transmembrane Protease Serine Subtype 2 (TMPRSS2) Interacting Approved Drugs for Coronavirus Disease 2019 (COVID-19) Treatment through a Drug-Target Interaction Deep Learning Model. Viruses. 2020; 12(11):1325. https://0-doi-org.brum.beds.ac.uk/10.3390/v12111325

Chicago/Turabian Style

Choi, Yoonjung, Bonggun Shin, Keunsoo Kang, Sungsoo Park, and Bo Ram Beck. 2020. "Target-Centered Drug Repurposing Predictions of Human Angiotensin-Converting Enzyme 2 (ACE2) and Transmembrane Protease Serine Subtype 2 (TMPRSS2) Interacting Approved Drugs for Coronavirus Disease 2019 (COVID-19) Treatment through a Drug-Target Interaction Deep Learning Model" Viruses 12, no. 11: 1325. https://0-doi-org.brum.beds.ac.uk/10.3390/v12111325

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