key: cord-0953440-fhmr1og4 authors: Ke, Yi-Yu; Peng, Tzu-Ting; Yeh, Teng-Kuang; Huang, Wen-Zheng; Chang, Shao-En; Wu, Szu-Huei; Hung, Hui-Chen; Hsu, Tsu-An; Lee, Shiow-Ju; Song, Jeng-Shin; Lin, Wen-Hsing; Chiang, Tung-Jung; Lin, Jiunn-Horng; Sytwu, Huey-Kang; Chen, Chiung-Tong title: Artificial intelligence approach fighting COVID-19 with repurposing drugs date: 2020-05-15 journal: Biomed J DOI: 10.1016/j.bj.2020.05.001 sha: 6e8e38f33035fd67b4e603ff02814ef1c112375d doc_id: 953440 cord_uid: fhmr1og4 BACKGROUND: The ongoing COVID-19 pandemic has caused more than 193,825 deaths during the past few months. A quick-to-be-identified cure for the disease will be a therapeutic medicine that has prior use experiences in patients in order to resolve the current pandemic situation before it could become worsening. Artificial intelligence (AI) technology is hereby applied to identify the marketed drugs with potential for treating COVID-19. MATERIAL AND METHODS: An AI platform was established to identify potential old drugs with anti-coronavirus activities by using two different learning databases; one consisted of the compounds reported or proven active against SARS-CoV, SARS-CoV-2, human immunodeficiency virus, influenza virus, and the other one containing the known 3C-like protease inhibitors. All AI predicted drugs were then tested for activities against a feline coronavirus in in vitro cell-based assay. These assay results were feedbacks to the AI system for relearning and thus to generate a modified AI model to search for old drugs again. RESULTS: After a few runs of AI learning and prediction processes, the AI system identified 80 marketed drugs with potential. Among them, 8 drugs (bedaquiline, brequinar, celecoxib, clofazimine, conivaptan, gemcitabine, tolcapone, and vismodegib) showed in vitro activities against the proliferation of a feline infectious peritonitis (FIP) virus in Fcwf-4 cells. In addition, 5 other drugs (boceprevir, chloroquine, homoharringtonine, tilorone, and salinomycin) were also found active during the exercises of AI approaches. CONCLUSION: Having taken advantages of AI, we identified old drugs with activities against FIP coronavirus. Further studies are underway to demonstrate their activities against SARS-CoV-2 in vitro and in vivo at clinically achievable concentrations and doses. With prior use experiences in patients, these old drugs if proven active against SARS-CoV-2 can readily be applied for fighting COVID-19 pandemic. use experiences in patients, these old drugs if proven active against SARS-CoV-2 can readily be applied for fighting COVID-19 pandemic. The outbreak of a new coronavirus SARS-CoV-2, previously known as the 2019-nCoV, has caused 2.81 million confirmed cases of the infection. There are more than 193,825 patients died of the coronavirus disease 2019 (COVID-19) affecting 210 countries and territories and the death toll is still climbing up sharp. SARS-CoV-2 was first reported found in Wuhan, China and, soon after, a human-to-human transmission is confirmed followed by a global COVID-19 pandemic as WHO declared on the 11th of March. Majority of the patients of COVID-19 present with fever and respiratory symptoms and death in some serious complications [1] . There are a number of clinical studies ongoing for treating COVID-19 with potential old drugs and investigational new drugs. In order to quickly resolve the current pandemic with more medicinal cures, repurpose of old drug is one of the strategies highly appreciated. We have hereby practiced with an artificial intelligence (AI) system for searching marketed drugs potentially with antiviral activities against coronaviruses. An AI system could be an efficient tool to quickly screen large number of compounds with assigned learning datasets and an intention to find drugs for specific purposes such as for treating COVID-19. The established AI system analyzed the most important descriptors, generated relevant prediction models, performed screening for the marketed drugs to quickly identify those with potential activities inhibiting SARS-CoV-2. An in vitro cell model for feline coronavirus replication was set up to evaluate the AI-identified drugs for antiviral activity verification. Feline coronavirus is an a-coronavirus and the virus causing enteritis in domestic and wild cats. Approximately 5e15% of infected cats develop feline infectious peritonitis (FIP) which is fatal to cats [2] . The infection by FIP virus in cats presented similar features to the severe acute respiratory syndrome (SARS) infection such as pulmonary lesions in humans [3] . It was documented that both the nucleoside analog GS-441524 and 3C-like protease inhibitor GC376 exhibited antiviral activities against FIP virus in vitro and were effective for treating FIP in cats [4, 5] . Remdesivir (GS-5734) is the prodrug of GS-441524 that inhibits viral RNA-dependent RNA polymerase and prevents viral replication in vitro including the middle east respiratory syndrome (MERS) virus, Ebola virus, Lassa fever virus, Junin virus and respiratory syncytial virus [6] . Remdesivir is now an investigational new drug that shows promising outcomes in compassionate uses. The use of FIP virus in vitro replication cell model in the Table 1 Learning datasets for Al training. present study is therefore relevant and shall provide a useful screening tool to identify marketed drugs with a broadspectrum antiviral activity. Two different learning datasets were generated in which one is consisted of the compounds reported or proven active against SARS-CoV, SARS-CoV-2, human immunodeficiency virus (HIV), and influenza virus and the other contains the known 3C-like protease inhibitors. An AI-system was established and, based on the learning datasets to predict drugs potentially active against coronavirus out of the marketed drugs. The predicted drugs were then tested for activities against a feline coronavirus in in vitro cell-based assay for a verification. During the AI practices, these assay results were served as feedbacks to the AI system for relearning and thus to generate a modified AI model to search for old drugs. Three types of molecular descriptors, extended connectivity fingerprint (ECFP), functional-class fingerprints (FCFPs), and octanolewater partition coefficient (ALogP_count), were performed to AI learning. The extended connectivity fingerprints (ECFPs) are generated in a molecule-directed manner by systematically recording the neighborhood of each non-hydrogen atom into multiple circular layers up to a given diameter of that molecule. The functional-class fingerprints (FCFPs) [7] detail circular fingerprints via the pharmacophore identification of atoms, which reports topological pharmacophore fingerprints. The ALogP_count is an array of 120 numbers that correspond to the 120 Ghose and Crippen atom types [8] . The 613 descriptors in total were used for the AI practices and prediction of the promising drugs. All these descriptors were generated by Discovery Studio v18.1/Calculates ligand properties program (BIOVIA Inc., San Diego, CA, USA), including ALogP_count (101 descriptors), ECFP_4 (256 descriptors), and Fig. 1 Flow scheme of AI approaches along with antiviral activity verification assays. Two independent datasets were compiled as the learning inputs to generate two AI prediction models of different approaches. Firstly, Model 1 was generated form the known drugs with antiviral activities. Model 2 was generated form the 3C-like protease inhibitors. A database for the marketapproved drugs were adopted in which the AI system screens for potential drugs with antiviral activities. The AI-predicted drugs were verified with antiviral activities by a cell-based FIP virus replication assay. Finally, these assay results served as feedbacks for the AI relearning and evolution progress. The Modified AI model was established to screen further and again verified by the FIP virus replication assay. The processing flows are indicated with the arrows. 3 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60 61 62 63 64 65 66 67 68 69 70 71 72 73 74 75 76 77 78 79 80 81 82 83 84 85 86 87 88 89 90 91 92 93 94 95 96 97 98 99 100 101 102 103 104 105 106 107 108 109 110 111 112 113 114 115 116 117 118 119 120 121 122 123 124 125 126 127 128 129 130 BJ291_proof ■ 15 May 2020 ■ 3/8 FCFP_4 (256 descriptors) as shown in the supplemental information. The system used a Deep Neural Network (DNN) [9] algorithm to identify the most important descriptors and gave different weightings to generate AI prediction models. We collected information of those compounds and drugs that had been reported with promising activities for treating COVID-19, such as anti-influenza drugs [10] , drugs and compounds inhibiting SARS-CoV [11] , HIV drugs [12] , and some showing activities against SARS-CoV-2 as recently reported [13, 14] . The other AI prediction model is generated from a learning dataset of 210 drugs and compounds reported with inhibitory activities on the 3C-like protease of SARS-CoV listed in Table 1 [15e35]. The database adopted to search for marketed drugs was the DrugBank at website: https://www.drugbank.ca. The serotype II FIP virus of Taiwan isolate NTU156 strain, was a kind gift from National Taiwan University [36] . The virus was propagated and titrated in the feline catus whole fetus-4 (Fcwf-4) cells (ATCC®CRL-2787). Fcwf-4 cells were amplified and maintained in Dulbecco's modified Eagle's medium (DMEM) containing 10% fetal bovine serum (FBS) with 1% penicillin/streptomycin at 37 C with 5% CO 2 . Confluent Fcwf-4 cells were seeded in 96-well plates and were treated with various concentrations of testing compounds at up to 50 mM at 37 C with 5% CO 2 for 48 h. Cell cytotoxicity was visualized and measured by crystal violet staining. In brief, 3 Â 10 4 Fcwf-4 cells per well were seeded in 96-well plates containing DMEM supplemented with 10% FBS and incubated at 37 C with 5% CO 2 . Sixteen hr post inoculation, the cells were infected with NTU156 at 300 TCID 50 per well and incubated at 37 C for 1 h and the supernatant were discarded followed by addition of 0.4, 2, 10, and 50 mM of testing compounds in DMEM medium containing 2% FBS. The culture plates were incubated for additional 48 h at 37 C with 5% CO 2 followed by fixing with 10% formalin and stained with 0.1% crystal violet 20% ethanol (w/v). The cytopathic effect (CPE) caused by FIP virus were measured visually for intensity of crystal violet staining. Back titration of FIP virus were performed parallelly and the TCID 50 values were calculated according a Reed and Muench method [37] . Results from AI model 1 approach AI prediction model 1 was exercised based on the learning dataset listed in Table 1 , containing some of the drugs for influenza, from previous reports showing activities against SARS-CoV, HIV drugs, and drugs and compounds of potential inhibiting SARS-CoV-2 as recently reported as shown in Fig. 1 . The AI model identified 22 drugs to be having the potential of antiviral effects. Followed by evaluations on these 22 drugs with the in vitro cell-based FIP virus replication assay, 2 promising drugs clofazimine and gemcitabine were verified with activity against the FIP coronavirus (Table 2) Q2 (see Fig. 2 ). An independent approach from the AI model 1, AI prediction model 2 was practiced based on the learning dataset of 210 inhibitors of the 3C-like protease of SARS-CoV listed in Table 1 . The AI model predicted 38 drugs with potential inhibiting the 3C-like protease of SARS-CoV and they were subject to tests in the FIP virus replication assay. Two out of the 38 4 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60 61 62 63 64 AI-identified drugs, celecoxib and tolcapone, showed antiviral activities against the FIP coronavirus (Table 3) . The results from the two AI models 1 and 2 were shown in Tables 2 and 3 , respectively, and were served together as a feedback dataset to the AI system for relearning and therefore a modified AI prediction system was established and used for the follow-up prediction exercise. During our integrated efforts for repurposing drugs to tackle COVID-19, there were additional drugs found active in our in vitro assay against the FIP virus (Table 4) , such as boceprevir [38] , chloroquine [39] , homoharringtonine [40] , tilorone [41] , and salinomycin [42] as they were also reported previously with potential to inhibit SARS-CoV or SARS-CoV-2. These 5 drugs were included as a sub-dataset for the relearning and thus to generate the modified AI prediction system which was used to screen and identify those 20 drugs out of the approved 2684 drugs for tests in the FIP virus replication assay (Table 5 ). With the AI evolution, the system showed a higher fidelity and 4 drugs (bedaquiline, brequinar, conivaptan, and vismodegib) out of the 20 AI-identified were verified with antiviral activities against the FIP coronavirus. Predicted by the AI system, 8 drugs were identified with antiviral activities against FIP virus proliferation in the Fcwf-4 cells. Together with the 5 additional drugs identified during the current AI approach and drug repurposing practices, there are 13 drugs with great potentials for further development toward treating COVID-19. Chloroquine showed inhibition activity against FIP virus at 10 mM, which is in agreement with the recent report about its broad-spectrum antiviral effects and it is one of the drugs used to treat COVID-19 patients in China and authorized for emergency use in United States [14] . Remdesivir is a Table 3 Assay results of the drugs predicted by AI Model 2. Cytotoxicity Viral Inhibition >0.4 NA a NA, not available to detect due to cytotoxicity. Table 2 Assay results of the drugs predicted by AI Model 1. Cytotoxicity Viral Inhibition a NA, not available to detect due to cytotoxicity. b GC376, a reference compound. b i o m e d i c a l j o u r n a l x x x ( x x x x ) x x x 5 SARS-CoV-2 RNA-dependent RNA polymerase inhibitor and a promising drug candidate for treating COVID-19 [4] ; now used in COVID-19 patients in clinical trials in a number of countries including Taiwan. GS-441524 is the active metabolite of remdesivir and showed activity against FIP virus at around 2 mM similar to the previously reported for GS-441524 being developed as a veterinary medicine treating feline coronavirus infections [4] . The results therefore provide solid evidence to support the validity of using cell-based FIP virus infection in vitro model as an antiviral activity assay in the present study. To have the promising drugs further being tested in a SARS-CoV-2 in vitro assay and in vivo animal models is planned to be ongoing. Among these AI-identified active drugs, four drugs (brequinar, clofazimine, gemcitabine and conivaptan) are of potential in their clinical uses and hereby elaborated further. Brequinar is a synthetic quinolinecarboxylic acid analogue and a dihydroorotate dehydrogenase (DHODH) inhibitor with experiences of patients uses for antitumor activities [43] . It was also used as an immunosuppressant in combination with other drugs such as cyclosporine for synergism [44] . It is of high water solubility and given by IV infusion. The biodistribution studies in rats of brequinar showed that liver and small intestines were the most exposed organs. The elimination rate of brequinar from all tissues were approximately equal to that from the blood circulation indicating no tissue accumulation of the intact drug molecules [45] . The present finding demonstrated that brequinar has an inhibiting activity of 2e10 mM. The systemically exposed plasma concentration range of brequinar in patients are around 73e261 mM for doses of 135e300 m 2 /kg in 5-day infusion period suggesting that brequinar may be effective against coronaviruses in appropriate dose regimens given alone or in combination with others [44] . Clofazimine is a water-soluble dye developed for the treatment of tuberculosis and used in combination for leprosy. It has a long circulation half-life of 25 days and the serum concentrations remain relatively low with evident tissue accumulation to extremely high levels approximately as 1.5 mg/g in jejunum. After 20 weeks treatment, the peak concentrations in lungs were 106 mg/g and 9020 mg/g in the lungs and spleen, respectively, in mice [46] . Clofazimine exhibited a dose-related relationship of tissue and serum concentrations (1.56e2.23 mM) from an oral dose of 100e300 mg in patients. A long onset of at least two weeks was observed in patients. Taking the advantages of the very slow clearance from body long half-life, and high drug concentration in the lungs [47] , clofazimine with minimum adverse effects could be further investigated for treating COVID-19. Conivaptan is a dual V1A/V2 vasopressin antagonist and given intravenously for the treatment of euvolemic and hypervolemic hyponatremia. The peak concentration C 0.5hr of the loading dose at 40 mg/day was about 1.4 mM and the maintenance dose at 40 mg/day for 4-day infusion was around 0.4 mM in patients [48] . Although these concentrations were lower than the inhibitory concentrations at approximately 50 mM, a clinical regimen may be manageable to seek possible efficacy and tolerability in COVID-19 patients with proper study design. Gemcitabine is an antimetabolite and first-line anticancer drug for pancreatic cancer. It also exhibits activities against breast cancer, ovarian cancer, non-small cell lung cancer, and bladder cancer. Gemcitabine is currently used alone or combined in cancer patients received 600e1200 mg/m 2 gemcitabine of 60e120 min i.v. infusion patients in treatment cycles with achievable plasma concentrations of 25e139 mM [49] . Side effects of gemcitabine include pneumonitis, pulmonary fibrosis, pulmonary edema, and adult respiratory distress syndrome. These pulmonary events could lead to fatal respiratory failure in some cases despite discontinuation of gemcitabine therapy 2 weeks ahead. A treatment regimen of low dose and frequency may be explored for single use against COVID-19 in animal models and for combined uses with the other drugs [50] . The present study revealed that the inhibitory activity of gemcitabine on coronavirus FIPV was approximately at 2 mM. A proper combined use with the other drugs in a lower dose of gemcitabine may overcome the pulmonary adverse effects while active inhibiting SARS-CoV-2 for cure of COVID-19 in further preclinical and clinical studies. AI approach is an efficient way to help quickly identify potential drugs as antiviral therapeutics. Some of the AI-identified drugs show activities against the proliferation of feline coronavirus in vitro. At the time of writing, there are still thousands of people dying of COVID-19 pandemic every day and we hope any of these AI-identified old drugs could soon make it to an effective cure for COVID-19. A special government fund to tackle the COVID-19 pandemic from the Q3 Ministry of Health and Welfare, Taiwan. 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60 61 62 63 64 65 66 67 68 69 70 71 72 73 74 75 76 77 78 79 80 81 82 83 84 85 86 87 88 89 90 91 92 93 94 95 96 97 98 99 100 101 102 103 104 105 106 107 108 109 110 111 112 113 114 115 116 117 118 119 120 121 122 123 124 125 126 127 128 129 130 BJ291_proof ■ 15 May 2020 ■ 7/8 SARS-CoV active drugs -Promazine compounds) Results from AI model 2 The continuing 2019-nCoV epidemic threat of novel coronaviruses to global health -the latest 2019 novel coronavirus outbreak in Wuhan, China Serologic studies of naturally occurring feline infectious peritonitis Human severe acute respiratory syndrome (SARS) and feline coronaviruses Efficacy and safety of the nucleoside analog GS-441524 for treatment of cats with naturally occurring feline infectious peritonitis Efficacy of a 3C-like protease inhibitor in treating various forms of acquired feline infectious peritonitis Broad-spectrum antiviral GS-5734 inhibits both epidemic and zoonotic coronaviruses Extended-connectivity fingerprints Atomic physicochemical parameters for three-dimensional-structure-directed quantitative structure-activity relationships. 2. Modeling dispersive and hydrophobic interactions Deep neural nets as a method for quantitative structure-activity relationships New treatments for influenza Pharmacokinetics of anti-SARS-CoV agent niclosamide and its analogs in rats Anti-HIV drugs for cancer therapeutics: back to the future? Traditional Chinese medicine is a resource for drug discovery against 2019 novel coronavirus (SARS-CoV-2) Clinical trials on drug repositioning for COVID-19 treatment. Rev Panam Salud Publ High-throughput screening identifies inhibitors of the SARS coronavirus main proteinase Inhibition of SARS-CoV 3C-like protease activity by theaflavin-3,3'-digallate (TF3). Evid Based Complement Discovering severe acute respiratory syndrome coronavirus 3CL protease inhibitors: virtual screening, surface plasmon resonance, and fluorescence resonance energy transfer assays Synthesis and evaluation of isatin derivatives as effective SARS coronavirus 3CL protease inhibitors Design, synthesis and antiviral efficacy of a series of potent chloropyridyl ester-derived SARS-CoV 3CLpro inhibitors Structure-based design, synthesis, and biological evaluation of peptidomimetic SARS-CoV 3CLpro inhibitors Discovery, synthesis, and structure-based optimization of a series of N-(tert-butyl)-2-(N-arylamido)-2-(pyridin-3-yl) acetamides (ML188) as potent noncovalent small molecule inhibitors of the severe acute respiratory syndrome coronavirus (SARS-CoV) 3CL protease Synthesis and evaluation of keto-glutamine analogues as potent inhibitors of severe acute respiratory syndrome 3CLpro Identification of novel small-molecule inhibitors of severe acute respiratory syndrome-associated coronavirus by chemical genetics Individual and common inhibitors of coronavirus and picornavirus main proteases Structure-based drug design and structural biology study of novel nonpeptide inhibitors of severe acute respiratory syndrome coronavirus main protease Structurebased virtual screening against SARS-3CL(pro) to identify novel non-peptidic hits Design, synthesis, and biological evaluation of novel dipeptide-type SARS-CoV 3CL protease inhibitors: structure-activity relationship study Development of potent dipeptide-type SARS-CoV 3CL protease inhibitors with novel P3 scaffolds: design, synthesis, biological evaluation, and docking studies Discovery of a novel family of SARS-CoV protease inhibitors by virtual screening and 3D-QSAR studies Specific plant terpenoids and lignoids possess potent antiviral activities against severe acute respiratory syndrome coronavirus Small molecules targeting severe acute respiratory syndrome human coronavirus Design and synthesis of cinanserin analogs as severe acute respiratory syndrome coronavirus 3CL protease inhibitors Aryl methylene ketones and fluorinated methylene ketones as reversible inhibitors for severe acute respiratory syndrome (SARS) 3C-like proteinase Design, synthesis, and evaluation of inhibitors for severe acute respiratory syndrome 3C-like protease based on phthalhydrazide ketones or heteroaromatic esters Full genome analysis of a novel type II feline coronavirus NTU156 A simple method of estimating fifty per cent endpoints A new lead for nonpeptidic active-sitedirected inhibitors of the severe acute respiratory syndrome coronavirus main protease discovered by a combination of screening and docking methods Remdesivir and chloroquine effectively inhibit the recently emerged novel coronavirus (2019-nCoV) in vitro Remdesivir, lopinavir, emetine, and homoharringtonine inhibit SARS-CoV-2 replication in vitro Tilorone: a broad-spectrum antiviral invented in the USA and commercialized in Russia and beyond Development of small-molecule MERS-CoV inhibitors The emergence of dihydroorotate dehydrogenase (DHODH) as a therapeutic target in acute myeloid leukemia Phase-I clinical and pharmacokinetic trial of brequinar sodium (Dup 785 Nsc 368390) Distribution of the novel anticancer drug candidate brequinar sodium (Dup-785, Nsc-368390) into Normal and tumor-tissues of nude-mice bearing humancolon carcinoma xenografts Tissue concentration, systemic distribution and toxicity of clofazimine-an autopsy study Pharmacokinetics and pharmacodynamics of clofazimine in a mouse model of tuberculosis Pharmacokinetics of conivaptan use in patients with severe hepatic impairment Pharmacokinetic study of gemcitabine, given as prolonged infusion at fixed dose rate, in combination with cisplatin in patients with advanced non-small-cell lung cancer Capillary leak syndrome: etiologies, pathophysiology, and management We thank Dr. Kung-Yee Liang for his support and helpful editing on the study and manuscript writing. The authors declare that they have no competing interests. Supplementary data to this article can be found online at https://doi.org/10.1016/j.bj.2020.05.001. r e f e r e n c e s