sign in sign up
<<
Home , Epicriptine, Pentostatin, Simeprevir, TMPRSS2, Venetoclax, Zofenopril

viruses

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

Yoonjung Choi 1, Bonggun Shin 1 , Keunsoo Kang 2 , Sungsoo Park 1 and Bo Ram Beck 1,*

1 Deargen, Inc., Daejeon 34051, Korea; [email protected] (Y.C.); [email protected] (B.S.); [email protected] (S.P.) 2 Department of Microbiology, College of Natural Sciences, Dankook University, Cheonan 31116, Korea; [email protected] * Correspondence: [email protected]

 Received: 5 September 2020; Accepted: 16 November 2020; Published: 18 November 2020 

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 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 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, , 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- C virus (HCV) drugs, including , , and . Notably, of the top 30 drugs, AT1R blocker and neuropsychiatric drug showed similar 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.

Keywords: COVID-19; SARS-CoV-2; coronavirus; deep learning; drug repurposing; ACE2; TMPRSS2

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,

Viruses 2020, 12, 1325; doi:10.3390/v12111325 www.mdpi.com/journal/viruses Viruses 2020, 12, 1325 2 of 11

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 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-, tiotropium bromide, has an inhibitory effect on some of the SARS-CoV-2-induced 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 Viruses 2020, 12, 1325 3 of 11 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, Table1 and Table S1; TMPRSS2, Table1 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 Viruses 2020, 12, 1325 4 of 11 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 (Table2). Although AC 50 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 -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, (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 virus (HCV) inhibitors, suggesting that these compounds may exhibit antiviral effects on HCV and SARS-CoV-2.

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).

ACE2 TMPRSS2

Small Molecules Predicted Kd in nM Small Molecules Predicted Kd in nM 0.02 0.37 Liothyronine 0.43 Pentostatin 0.47 Emtricitabine 0.45 Tazemetostat 2.88 Tiotropium 0.92 Tiotropium 4.03 Obeticholic acid 0.98 Eluxadoline 5.08 Mestranol 1.02 Pimecrolimus 5.79 1.10 Tacrolimus 5.82 Levorphanol 1.30 Ombitasvir 5.91 Levallorphan 1.37 6.06 Methscopolamine 1.40 Venetoclax 6.12 Enalaprilat 1.46 Daclatasvir 6.66 Bremelanotide 1.49 Tolvaptan 7.33 Norethynodrel 1.51 Aclidinium 8.46 Epicriptine 1.56 Paritaprevir 8.75 Rolitetracycline 1.57 Eprosartan 9.19 Diacetyl benzoyl lathyrol 1.77 10.59 Pizotifen 1.80 11.08 Dapagliflozin 1.80 Lisuride 11.20 Clobetasol 1.80 Erdafitinib 11.51 1.86 Letermovir 12.43 Viruses 2020, 12, 1325 5 of 11

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.

Predicted MT-DTI Rank Alphalisa-AC50 Small Molecules TruHit-AC50 (µM) ACE2 Kd (nM) Out of 460 (µM) Enalaprilat 1.46 11 7.52 - Daclatasvir 5.81 79 5.97-8.00 4.41–6.22 Mifepristone 10.3 124 5.97 2.78 Estradiol cypionate 11.79 136 9.47 5.55 19.34 196 6.70 8.79 Ombitasvir 20.91 208 4.50–5.97 2.78–3.50 Norgestimate 21.06 210 9.47 4.41 d-alpha-Tocopherol acetate 21.18 212 10.62 2.21 25.85 243 0.40–0.60 0.25–0.31 Pasireotide 28.29 257 3.67 2.69 33.64 281 2.67 0.44 38.26 303 8.44 4.94 Venetoclax 44.77 324 8.44 3.50 Anthralin 61.92 380 2.38 6.22 Ciclopirox 88.25 439 2.67 1.56 88.98 440 5.66–6.70 3.12 Posaconazole 94.7 449 2.01–6.70 2.78–3.50

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 , 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 (Table3, 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 (Table3, 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. Viruses 2020, 12, 1325 6 of 11

Table 3. Cross-prediction results of ACE2 and TMPRSS2 interacting drugs through MT-DTI and AutoDock Vina.

ACE2 TMPRSS2 Predicted K AutoDock Vina Predicted K AutoDock Vina Small Molecules d Small Molecules d (nM) ∆G (kcal/mol) (nM) ∆G (kcal/mol) Bremelanotide 1.49 8.1 Tazemetostat 2.88 7.1 − − 7.04 8 Eluxadoline 5.08 7.5 − − 14.68 8 Entrectinib 11.08 7.9 − − 14.94 9.2 Erdafitinib 11.51 7.1 − − Tezacaftor 15.90 8 Aprepitant 13.52 7.5 − − Dutasteride 16.35 8.8 Canagliflozin 27.12 7 − − Rifapentine 17.81 8 Naldemedine 46.74 7.2 − − Acetyldigitoxin 18.93 8.6 Adapalene 52.49 7.5 − − Alatrofloxacin 24.46 8 68.38 7.7 − − Deslanoside 24.75 9.1 71.78 7.4 − − 25.99 8.4 75.18 7.8 − − 28.37 8.2 − Naldemedine 29.94 8.8 − Ciclesonide 31.43 8.1 − Ubrogepant 33.50 8.1 − Ergotamine 33.64 9.1 − Lumacaftor 36.26 8 − Venetoclax 44.77 9.2 − Adapalene 45.11 8 − Letermovir 50.73 8 − Paritaprevir 53.21 9.2 − Entrectinib 54.40 8.5 − Glycyrrhizic acid 61.94 8.1 − Simeprevir 67.03 8.7 − 75.81 8.5 − Lifitegrast 94.67 8.3 − Posaconazole 94.70 8 −

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, , 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- (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 (Table1). 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]) Viruses 2020, 12, 1325 7 of 11 and the top 30 drugs of MT-DTI through CMap Linked User Environment (CLUE, https://clue.io) (Table4). 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].

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. − Predicted TMPRSS2 K Bromhexine Probucol Small Molecules d (nM) Connectivity Score Connectivity Score Dasatinib 0.37 17.58 15.38 Tacrolimus 5.82 10.40 70.67 Eprosartan 9.19 91.13 51.17 Lisuride 11.20 34.70 92.18 Aprepitant 13.52 54.69 12.72 15.31 33.30 41.78 15.95 0.44 5.95

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 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, , , , , , , , , , 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 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. Viruses 2020, 12, 1325 8 of 11

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 http://www.mdpi.com/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.[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- (GS-5734) for the Treatment of Ebola and Emerging Viruses. J. Med. Chem. 2017, 60, 1648–1661. [CrossRef] 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. [CrossRef] Viruses 2020, 12, 1325 9 of 11

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. [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. [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 , inhibits SARS-CoV-2 replication and pro-inflammatory cytokine production. Antimicrob. Agents Chemother. 2020, 64.[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. [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. [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. [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. [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. [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. [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. [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. [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. [CrossRef][PubMed] 16. Sterling, T.; Irwin, J.J. ZINC 15 - Ligand Discovery for Everyone. J. Chem. Inf. Model. 2015, 55, 2324–2337. [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. 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. [CrossRef] 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. [CrossRef] 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. [CrossRef] 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. 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. [CrossRef][PubMed] 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. [CrossRef][PubMed] Viruses 2020, 12, 1325 10 of 11

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. [CrossRef] 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. [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.[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. [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. [CrossRef] 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. [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.[CrossRef] 31. Uno, Y. Camostat mesilate therapy for COVID-19. Intern. Emerg. Med. 2020, 15, 1577–1578. [CrossRef] [PubMed] 32. Elfiky, A.A. Anti-HCV, nucleotide inhibitors, repurposing against COVID-19. Life Sci. 2020, 248, 117477. [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. [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. [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. [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. [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. [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.[CrossRef] 39. Ledford, H. Coronavirus breakthrough: Dexamethasone is first drug shown to save lives. Nature 2020, 582, 469. [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.[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. [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. [CrossRef][PubMed] 43. Anighoro, A.; Bajorath, J.; Rastelli, G. Polypharmacology: Challenges and opportunities in drug discovery. J. Med. Chem. 2014, 57, 7874–7887. [CrossRef][PubMed] Viruses 2020, 12, 1325 11 of 11

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. [CrossRef][PubMed]

Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

© 2020 by the authors. Licensee MDPI, Basel, Switzerland. This article is an open access article distributed under the terms and conditions of the Creative Commons Attribution (CC BY) license (http://creativecommons.org/licenses/by/4.0/).