In silico Drug Repurposing to combat COVID-19 based on Pharmacogenomics of Patient Transcriptomic Data

doi:10.21203/rs.3.rs-39128/v1

Download PDF

Research Article

In silico Drug Repurposing to combat COVID-19 based on Pharmacogenomics of Patient Transcriptomic Data

https://doi.org/10.21203/rs.3.rs-39128/v1

This work is licensed under a CC BY 4.0 License

Journal Publication

published 01 Jan, 2021

Read the published version in Clinical Microbiology and Infectious Diseases →

Version 1

posted

You are reading this latest preprint version

The ongoing global pandemic of coronavirus disease 2019 (COVID-19) continues to affect a growing number of populations in different parts of the world. In the current situation, drug repurposing is a viable strategy to combat COVID-19. The drugs targeting the host receptors that interact with SARS-CoV-2 are possible candidates. However, assessment of their effectiveness in COVID-19 patients is necessary before prioritizing them for further study. We attempted to shortlist the candidate drugs using an in-silico approach. First, we analysed two published transcriptomic data sets of COVID-19- and SARS-infected patients compared to healthy individuals to find the key pathways altered after infection. Then, using publicly available drug perturbational data sets in human cell lines from the Broad Institute Connectivity Map (CMAP), we assessed the effects of the approved drugs on the altered pathways. We also used the available pharmacogenomic data sets from the Genomics of Drug Sensitivity in Cancer (GDSC) portal to assess the effects of the altered pathways on resistance or sensitivity to the drugs in human cell lines. Our analysis identified many candidate drugs, some of which are already being investigated for treatment of COVID-19 and can serve as a basis for prioritizing additional viable candidate drugs for COVID-19.

Systems Biology

SARS-CoV-2

COVID-19

drug repurposing

patient transcriptomics

pharmacogenomics

in-silico

In December 2019, outbreak of a new virus causing severe respiratory disease was reported in Wuhan, China. The disease was attributed to infection by a novel member of the coronavirus family, named SARS-CoV-2¹. In contrast with the two other highly pathogenic coronaviruses, SARS-CoV (outbreak in 2002) and MERS-CoV (outbreak in 2012), the novel virus SARS-CoV-2 had a lower fatality rate, but it could spread from person to person more efficiently². The outbreak spread worldwide³, and by March 2020, the World Health Organization declared the disease, by then named coronavirus disease 2019 (COVID-19), a global pandemic. COVID-19 is currently a leading public health concern in countries worldwide as the hospitalization rate and death toll increase daily (https://www.who.int/emergencies/diseases/novel-coronavirus-2019). Apart from mitigation efforts, there is an urgent need for effective treatments for patients showing severe symptoms. After Chinese researchers published the first sequence of SARS-CoV-2 in January⁴, many groups worldwide began working to develop an effective vaccine. Some of the candidates are undergoing clinical trials, among them two promising candidates from Oxford University and Moderna, Inc., that are entering Phase II⁵. However, clinical approval of the vaccines will take some time; consequently, the best solution in the interim is to look for drugs that have already been approved for humans and that may be repurposed for COVID-19 treatment. The U.S. Food and Drug Administration recently approved remdesivir, an antiviral drug targeting the viral RNA-dependent RNA polymerase, for COVID-19 treatment after a clinical trial run by the National Institute of Allergy and Infectious Diseases found it to reduce the time to recovery in treated patients compared to the control group⁶. However, the drug did not show significant improvement in patient survival. Thus, there is still a persistent need to search for alternative candidate drugs that can be used as single or as combination therapy to benefit patients with severe illness.

Researchers are also exploring targeting of the host receptors of SARS-CoV-2 as another potential solution. Inhibitors of the two main receptors for SARS-CoV-2/SARS-CoV-1 entry in human cells, ACE2 and TMPRSS2, are being tested as potential treatment options⁷. Other host interactors have been curated from previously published literature in a recent study⁸, and a list of 332 high-confidence interactions between SARS-CoV-2 and human proteins was recently published by Gordon et al., who used affinity purification mass spectrometry on a human cell line expressing the viral proteins⁹.

We collected the drugs targeting these SARS-CoV-2/SARS-CoV-1-interacting proteins from DrugBank¹⁰ and the Drug Gene Interaction Database (DGIdb)¹¹ and assessed their effect in COVID-19- and SARS-infected patients in silico using published patient transcriptomic data and drug perturbational data sets. By combining recently published COVID-19-infected patient transcriptomic data¹² and previous SARS-infected patient transcriptomic data¹³, we found key pathways that are altered in infected patients compared to healthy individuals. The common altered pathways were related to androgen biosynthesis, steroid biosynthesis, gluconeogenesis, regulation of mitosis, the proteasome, and inflammatory pathways. Next, using the available drug perturbational data sets from the Broad Institute Connectivity map (CMAP project¹⁴, we assessed how the candidate drugs targeting SARS-CoV-2-interacting proteins affect the pathways that are altered after COVID-19 or SARS infection in a time-dependent manner. Using another available pharmacogenomic data resource, Genomics of Drug Sensitivity in Cancer (GDSC)¹⁵, we assessed whether the altered pathways affect sensitivity or resistance to our candidate drugs in human cell lines. Finally, we present a list of potential drugs for COVID-19 treatment that affect the COVID-19-altered pathways and incur less resistance in the presence of alteration in the pathways. Figure 1 summarizes our approach.

Analysis of transcriptomic data of COVID-19- and SARS-infected patients revealed common altered pathways

To get an understanding of the transcriptomic alterations in COVID-19-infected patients, we analysed a published RNA sequencing (RNA-Seq) data set of peripheral blood mononuclear cells (PBMC) of three COVID-19 patients from China compared to three healthy individuals¹². As there were not enough patient transcriptomic data available for COVID-19-infected people, we searched the literature for transcriptomic data sets of SARS-infected patients. We included a microarray-based transcriptomic data set of PBMCs of 40 SARS-infected and 10 heathy individuals from a paper published in 2007¹³. For both of these transcriptomic data sets, we performed single-sample Gene Set Enrichment Analysis (ssGSEA) for a set of curated pathways collected from the Broad Institute Molecular Signatures Database (MSigDB) and a set of immune function pathways collected from published resources^16-18. While checking the differential enrichment of immune functions between COVID-19 vs. healthy as well as SARS vs. healthy patients, we found that the macrophages, monocytes, dendritic cells (MoMaDC), B cells, pathogen defence, and interferon signalling were the common immune functions with high enrichment in both COVID-19- and SARS-infected patients compared to healthy patients (Figure 2). On the other hand, among other pathways, androgen biosynthesis, steroid biosynthesis, potassium channel, proteasomes, glycolysis, gluconeogenesis, tryptophan metabolism, and p53-dependent DNA damage response were the common functions that were highly enriched in both COVID-19- and SARS-infected patients compared to healthy patients (Figure 3). Some of the common pathways with less enrichment in the COVID-19- and SARS-infected patients are NK-cell- and T-cell-mediated responses (Figure 2c, 3c).

Assessing the effect of potential drug candidates on the pathways altered in COVID-19 and SARS infection

In an attempt at finding already approved drugs that could be repurposed for COVID-19 treatment, we explored the list of all human proteins interacting with SARS-CoV-2 that were collected from the affinity capture mass-spectrometry study by Gordon et al.⁹ and the list of SARS-CoV-1-interacting proteins curated from previous literature⁸. We collected the list of drugs targeting these host proteins from DrugBank and DGIdb databases^10,11. Additionally, we procured the list of drugs currently in trial for COVID-19 treatment from literature (https://www.excelra.com/COVID-19-drug-repurposing-database/data). The resulting list comprised 289 potential drugs for COVID-19. Using the drug perturbational data sets from CMAP¹⁴, we found the drugs that affected the gene expression of these targets (differences in gene expression at 24 hours vs. 6 hours of drug treatment). We mapped 221 drugs to the drug perturbational data set from CMAP. The network in Figure 4a shows the drugs and SARS-CoV targets, which get significantly downregulated at 24 hours compared to six hours of the respective drug treatments (76 drugs) along with the drugs in clinical trial for COVID-19 treatment (Figure 4a, Supplementary Dataset 1 shows the list of drugs and their targets). Many of these drugs have been already reported as potential treatments for COVID-19 and/or SARS (Table 1), including sirolimus, apicidin, mercaptopurine, midostaurin, camostat, metformin, fedratinib, ribavirin, entacapone, azacitidine, and chloroquine. Next, we sought to check if the drugs are effective in the context of transcriptomic changes in COVID-19 or SARS patients after infection. We checked the effect of these drugs on the perturbation of the pathways upregulated or downregulated after SARS and COVID-19 infection (Figure 4b). In CMAP, gene expression is profiled in cell lines treated with drugs at different dose/time points, and we see how drugs change the expression of a given set of genes. To get a pathway-based estimation instead of individual genes, we calculated ssGSEA scores for the differentially enriched pathways in COVID-19 or SARS infection for all drug-treated cell lines at different dose/time points. Then the ssGSEA scores of a pathway were compared between two time points (24 hours vs. 6 hours) in cells treated with a certain drug using student’s t-test. Among all drugs targeting SARS-CoV-2-interacting proteins, 60 drugs altered the COVID-19/SARS differentially enriched pathways in a time-dependent manner. In Figure 4b (see Supplementary Dataset 2 for related data), we show drugs that 1) significantly downregulate the pathways that get upregulated after COVID-19/SARS infection (blue cells on the left panel), 2) upregulate the pathways that get upregulated after COVID-19/SARS infection (red cells on the left panel), 3) upregulate the pathways that get downregulated after COVID-19/SARS infection (red cells on the right panel), and 4) downregulate the pathways that get downregulated after COVID-19/SARS infection (blue cells on the right panel). Among these drugs, vorinostat, panobinostat, tacrolimus, lovastatin, midostaurin, daunorubicin, cytarabine, digitoxin, digoxin, mercaptopurine, bortezomib, azacitidine, and clorafarabine can suppress the pathways that were upregulated due to COVID-19/SARS infection (p-value < 0.05, t-test with alternative hypothesis = less), while they can also boost the pathways that were downregulated due to COVID-19/SARS infection (p-value < 0.05, t-test with alternative hypothesis = greater). As examples, Figure 4c and 4d show the effect of the drug cytarabine on the proteasome (upregulated in COVID-19/SARS infection) and CTL (downregulated in COVID-19/SARS infection) pathways at 24 hours vs. 6 hours of drug treatment. In contrast, another class of drugs, including entacapone, ribavirin, decitabine, menadione, chloroquine, dacinostat, and belinostat, significantly upregulated the pathways upregulated in COVID-19 or SARS infection.

Assessing the effect of potential drug candidates on the pathways altered in COVID-19 and SARS infection

Drug treatment sensitivity or resistance also depends on the genomic and transcriptomic features of tissues. The upregulated or downregulated pathway signatures of virus-infected patients can impact the sensitivity or resistance to the drugs targeting SARS-Cov-2-interacting proteins. Considering this, we sought to check whether the pathways that get altered after SARS or COVID-19 infection have any effect on drug resistance. Drug sensitivity is calculated from GDSC¹⁵. If the upregulated (or downregulated) pathways are associated with reduced half maximal inhibitory concentration (IC50) of a drug targeting a SARS-CoV-interacting protein, then we can infer that the pathway alteration makes the cell lines less resistant to drug treatment. Figure 5a (see Supplementary Dataset 3 for related data) shows the network of drug and pathway relations from this analysis, where upregulation (or downregulation) of these pathways is associated with reduced resistance to the drugs in GDSC cell lines. The network includes drugs targeting SARS-CoV-interacting proteins and common pathways altered (upregulated or downregulated) in COVID-19/SARS infection. Figure 5b illustrates that the high enrichment of NLRP3 inflammosomes (upregulated in COVID-19- and/or SARS-infected patients) is associated with reduced resistance (low IC50) to the drug fedratinib (calculated from GDSC data). Similarly, Figure 5c and 5d show reduced resistance to the drugs cytarabine and temsirolimus in GDSC cell lines as affected by high enrichment of the pathways’ cell cycle and glycolysis.

Considering the time needed for developing a functional vaccine, the feasible way to fight the COVID-19 pandemic is repurposing drugs already approved for clinical use. We tried to address this issue by identifying drugs that target potential host proteins interacting with SARS-CoV-2 and predicting their effect on the patient transcriptomic signatures of COVID-19- and/or SARS-infected patients. We analyzed COVID-19- and SARS-infected patients’ and healthy individuals’ transcriptomic data sets^12,13 to identify the pathways commonly altered due to infection with these viruses. Thereafter, combining the potential human interactome of SARS-CoV-2 from a recently published study⁹ and SARS-CoV-1-interacting proteins curated in another publication⁸ with drug target databases^10,11, drug perturbational data sets¹⁴, and drug sensitivity screening data sets¹⁵, we propose a map of the drugs that can be effective in COVID-19 treatment. Our proposed framework for assessing drug efficacy is based on two approaches: 1) whether the drugs can regulate the pathways that get altered due to COVID-19/SARS infection and 2) whether the sensitivity or resistance to the drugs depends on the alteration of these pathways.

We acknowledge that there are some limitations to our study. First, we found only one publicly available data set for transcriptomic signatures of COVID-19-infected patients, and that study only included three patients and three healthy individuals¹². Since it would be difficult to draw meaningful conclusions on the patient transcriptomic changes from so few patient samples, we included a previously published data set of SARS-infected patient transcriptomes consisting of 40 SARS-infected patients and 10 healthy individuals¹³. We reasoned that the viruses SARS-CoV-1 (causing SARS) and SARS-CoV-2 (causing COVID-19) share many similarities, and the symptoms of the infected patients are also similar, though SARS infection causes more severe symptoms and a higher fatality rate¹⁹. Another limitation of our study is the availability of drug perturbational and drug sensitivity screening data in human cell lines. Many of the drugs targeting the potential SARS-CoV-2- or SARS-CoV-1-interacting proteins did not have drug perturbational data available in human cell lines from CMAP, and only a few drugs have the drug sensitivity screening data available from GDSC, which is limited to drugs commonly used in cancer treatment. Without the available drug perturbational data, we could not assess the effects of antiviral drugs like remdesivir, favipiravir, and oseltamivir in the human transcriptome and were limited to drugs targeting the potential human receptors of the viruses.

Our analysis of COVID-19- and SARS-infected patient transcriptomic signatures identified the common pathways upregulated or downregulated after infection with these viruses. Among the common pathways, androgen biosynthesis, steroid biosynthesis, proteasome, tryptophan metabolism, and mitotic cell cycle-related pathways were upregulated due to infection, while TCR, Interleukin 2, TGF-beta receptor signaling, and beta-cell development pathways were downregulated. Among the immune pathways, MoMaDC, B cells, and interferon responses were elevated in patients with COVID-19 and SARS infection, while T-cell and NK cell functions were downregulated. Many of these pathways have been previously implicated in viral infections. Androgens have already been linked to the severity of COVID-19 as possible mechanisms driving disease outcomes²⁰. Steroid, proteasome, and tryptophan metabolisms are linked to viral infections and subsequent host immune response^21-23. Moreover, monocyte and macrophage activation has been linked to the hyperinflammatory phenotypes in SARS-CoV-2 infection that have adverse clinical outcomes²⁴. By analysing the drug perturbational data sets from CMAP, we identified drugs that can upregulate or downregulate the pathways that get altered due to COVID-19/SARS infection. As at this point, it’s unclear which pathways will be clinically beneficial for COVID-19 treatment when inhibited or elevated. We provided the full map of all drugs that 1) downregulate the pathways that get upregulated due to COVID-19/SARS infection, 2) upregulate the pathways that get upregulated due to COVID-19/SARS infection, 3) upregulate the pathways that get downregulated due to COVID-19/SARS infection, and 4) downregulate the pathways that get downregulated due to COVID-19/SARS infection. Some drugs that could downregulate the androgen biosynthesis pathway (upregulated due to COVID-19/SARS infection), which is linked to severe outcomes, are ponatinib, amlexanox, lenalidomide, menadione, and decitabine. The drugs vorinostat, panobinostat, tacrolimus, lovastatin, midostaurin, daunorubicin, cytarabine, digitoxin, digoxin, mercaptopurine, bortezomib, azacitidine, and clofarabine could downregulate the mitotic cell cycle-related pathways and proteasome activity. The mitotic cell-cycle arrest and uniquitin-proteasome pathway systems play an important role in viral replication for coronaviruses^22,25. Drugs that could downregulate the MoMaDC cell pathway potentially linked to severe symptoms in COVID-19 patients²⁴ are menadione, decitabine, amlexanox, mercaptopurine, nicotinamide, and pevonedistat. The drugs that could upregulate the CTL or Interleukin 2 pathways, which are important for antiviral immune responses, are: icosapent, lovastatin, azacitidine, decitabine, vorinostat, entinostat, digoxin, ponatinib, apicidin, daunorubicin, and metformin. Many of these drugs have been recently reported as possible treatment options for COVID-19, including ponatinib, mercaptopurine, azacitidine, apicidin, daunorubicin, and metformin⁹.

Considering that host transcriptomic features can affect the sensitivity or resistance to drugs¹⁵, we checked whether the pathways that get altered after COVID-19/SARS infection can reduce the resistance to the drugs under consideration using GDSC cell-line drug sensitivity screening data. We were limited to drugs approved for cancer treatment that have been screened in human cancer cell lines, so we could not test all potential drugs for COVID-19. However, among the drugs we tested, we found that there was reduced drug resistance (increased sensitivity) related to the pathway alteration signature due to COVID-19/SARS infection with the drugs navitoclax, bortezomib, fedratinib, lenalidomide, cytarabine, temsirolimus, paclitaxel, docetaxel, vorinostat, belinostat, entinostat, dacinostat, and vinblastine. Among these drugs, fedratinib is a potential treatment option for COVID-19²⁶. Fedratinib is a semi-selective janus kinase inhibitor with anti-inflammatory properties that can be beneficial for treatment of severe COVID-19 patients with hyperinflammatory symptoms that may lead to lung damage and death. Among the other cancer drugs, the histone deacetylase inhibitors vorinostat, belinostat, entinostat, and dacinostat have anti-inflammatory properties, and the drug lenalidomide has immunomodulatory properties that may have treatment benefits in COVID-19 patients with severe symptoms, possibly when treated with antiviral agents. The chemotherapeutic drugs paclitaxel, docetaxel, and cytarabine may also have potential antiviral properties^27,28; however, their suitability for treating virus infection is not well established, and therefore, these drugs should be carefully examined before being considered as treatment options.

Our computational study provides an assessment of the potential effects of possible repurposable drugs targeting SARS-CoV-2- and SARS-CoV-1-interacting proteins in the context of transcriptomic signatures of COVID-19- or SARS-infected patients. Because we have limited understanding of the virus pathogenesis and the safety of targeting specific pathways for treatment of this disease, we have included a complete map of the drug-pathway relationships that can be important in the context of COVID-19 infection and treatment. Due to the urgency of identifying effective drugs for treatment of COVID-19 patients, many groups have been working on repurposing existing drugs that can be used for this disease⁸. Our study may be the first that links potentially repurposable drugs to the pathway signatures of infected patients. The drugs we presented in this study as potentially repurposable for COVID-19 treatment should be clinically tested for efficacy before use.

Data sets and pre-processing

The RNA-Seq transcriptomic data of the PBMCs from three COVID-19 infected patients and three healthy individuals were obtained from the study by Xiong et al.¹². The microarray transcriptomic data of the PBMCs from 40 SARS-infected and 10 healthy patients were obtained from the study by Cameron et al.¹³. The RNA-Seq data from the COVID-19 patients and healthy controls was normalized to fractions per kilobase per million (FPKM). The microarray data from SARS patients and healthy controls was normalized with locally weighted scatterplot smoothing (LOWESS). Transcriptomic data from both COVID-19 and SARS studies were log2-transformed before further processing.

The list of human proteins interacting with SARS-CoV-2 was collected from the study by Gordon et al.⁹. In addition, the human proteins interacting with SARS-CoV are collected from the article by Zhou et al.⁸. The drugs targeting the SARS-CoV-2-interacting proteins were collected from DrugBank and DGIdb^10,11. The drug perturbational data set containing dose- and time-dependent transcriptomic profiles upon drug treatments in human cell lines was collected from CMAP¹⁴. The drug-screening data measuring IC50 upon drug treatments in human cell lines and cell-line transcriptomic profiles are collected from GDSC¹⁵.

Computational analysis and metadata

The general pathway gene sets were obtained from the C2 collection of MSigDB²⁹. We curated immunologic signatures specific to immune functions from four published resources: the LM22 immune infiltration signature¹⁶, the LM7 immune infiltration signature¹⁷, the ImSig signature of solid tumor immune infiltration¹⁸, and the NanoString immune signature panel (https://www.nanostring.com). The combined 3,129 curated gene sets and 61 immune gene sets were used for a single-sample gene-set enrichment analysis on the RNA-seq and microarray transcriptomic tumor data sets for COVID-19- and SARS-infected patient transcriptomes as well as for cell-line transcriptomic data for the CMAP drug perturbational studies and GDSC data sets. R-package gene-set variation analysis was used³⁰. Differential enrichments of pathways between COVID-19- and SARS-infected patients vs. healthy controls were calculated with limma R package using ssGSEA scores. The pathways with adjusted p-value < 0.05 are identified as differentially enriched.

A one-sided student’s t-test was used for calculating the p-value of differences in pathway ssGSEA scores between time-dependent drug treatment responses (from CMAP) at 24 hours vs. 6 hours. For pathways that are differentially regulated due to COVID-19/SARS infection, we checked whether a drug can significantly downregulate the pathway at 24 hours compared to 6 hours of transfection in human cell lines from CMAP by performing t-tests on the pathway ssGSEA scores between these two time points using an alternative hypothesis equal to less. Similarly, to check whether a drug can significantly upregulate the pathway at 24 hours compared to 6 hours of transfection, we performed t-tests on the pathway ssGSEA scores between these two time points using an alternative hypothesis equal to greater.

The pathway-dependent drug sensitivity/resistance for each pathway and drug combination from GDSC is calculated using the following approach. First, ssGSEA scores for the pathways are calculated for the GDSC cell lines. Then, the cell lines are stratified into two sets based on the ssGSEA scores of the pathways. For upregulated pathways, cell lines with pathway ssGSEA z-score > 2 have high pathway enrichment, and for downregulated pathways, cell lines with pathway ssGSEA z-score < -2 have low pathway enrichment. The difference between drug IC50 in GDSC cell lines with high ssGSEA scores (z-score > 2) for COVID-19/SARS upregulated pathways and all other GDSC cell lines was calculated using a one-sided t-test with an alternative hypothesis equal to greater. The difference between drug IC50 in GDSC cell lines with low ssGSEA scores (z-score < -2) for COVID-19/SARS downregulated pathways and all other GDSC cell lines was calculated using a one-sided t-test with an alternative hypothesis equal to less.

Acknowledgements

This research was supported by the Intramural Research Program of the National Institutes of Health, National Cancer Institute (NCI).

Author Contributions

SD and US conceived the study design. SD performed all analysis. SD, US and KC performed revision and manuscript writing.

Additional Information

The authors have no competing interests.

Funding

This research was funded by the Intramural Research Program of the National Institutes of Health, National Cancer Institute (NCI).

Data Availability

All analysis data generated during this study are included in this published article as supplementary datasets.

1 Wu, F. et al. A new coronavirus associated with human respiratory disease in China. Nature 579, 265-269, doi:10.1038/s41586-020-2008-3 (2020).

2 Gates, B. Responding to Covid-19 - A Once-in-a-Century Pandemic? N Engl J Med 382, 1677-1679, doi:10.1056/NEJMp2003762 (2020).

3 Wang, C., Horby, P. W., Hayden, F. G. & Gao, G. F. A novel coronavirus outbreak of global health concern. Lancet 395, 470-473, doi:10.1016/S0140-6736(20)30185-9 (2020).

4 Chan, J. F. et al. Genomic characterization of the 2019 novel human-pathogenic coronavirus isolated from a patient with atypical pneumonia after visiting Wuhan. Emerg Microbes Infect 9, 221-236, doi:10.1080/22221751.2020.1719902 (2020).

5 Lurie, N., Saville, M., Hatchett, R. & Halton, J. Developing Covid-19 Vaccines at Pandemic Speed. N Engl J Med 382, 1969-1973, doi:10.1056/NEJMp2005630 (2020).

6 Grein, J. et al. Compassionate Use of Remdesivir for Patients with Severe Covid-19. N Engl J Med, doi:10.1056/NEJMoa2007016 (2020).

7 Hoffmann, M. et al. SARS-CoV-2 Cell Entry Depends on ACE2 and TMPRSS2 and Is Blocked by a Clinically Proven Protease Inhibitor. Cell 181, 271-280 e278, doi:10.1016/j.cell.2020.02.052 (2020).

8 Zhou, Y. et al. Network-based drug repurposing for novel coronavirus 2019-nCoV/SARS-CoV-2. Cell Discov 6, 14, doi:10.1038/s41421-020-0153-3 (2020).

9 Gordon, D. E. et al. A SARS-CoV-2 protein interaction map reveals targets for drug repurposing. Nature, doi:10.1038/s41586-020-2286-9 (2020).

10 Wishart, D. S. et al. DrugBank 5.0: a major update to the DrugBank database for 2018. Nucleic Acids Res 46, D1074-D1082, doi:10.1093/nar/gkx1037 (2018).

11 Cotto, K. C. et al. DGIdb 3.0: a redesign and expansion of the drug-gene interaction database. Nucleic Acids Res 46, D1068-D1073, doi:10.1093/nar/gkx1143 (2018).

12 Xiong, Y. et al. Transcriptomic characteristics of bronchoalveolar lavage fluid and peripheral blood mononuclear cells in COVID-19 patients. Emerg Microbes Infect 9, 761-770, doi:10.1080/22221751.2020.1747363 (2020).

13 Cameron, M. J. et al. Interferon-mediated immunopathological events are associated with atypical innate and adaptive immune responses in patients with severe acute respiratory syndrome. J Virol 81, 8692-8706, doi:10.1128/JVI.00527-07 (2007).

14 Lamb, J. et al. The Connectivity Map: using gene-expression signatures to connect small molecules, genes, and disease. Science 313, 1929-1935, doi:10.1126/science.1132939 (2006).

15 Yang, W. et al. Genomics of Drug Sensitivity in Cancer (GDSC): a resource for therapeutic biomarker discovery in cancer cells. Nucleic Acids Res 41, D955-961, doi:10.1093/nar/gks1111 (2013).

16 Newman, A. M. et al. Robust enumeration of cell subsets from tissue expression profiles. Nat Methods 12, 453-457, doi:10.1038/nmeth.3337 (2015).

17 Tosolini, M. et al. Assessment of tumor-infiltrating TCRVgamma9Vdelta2 gammadelta lymphocyte abundance by deconvolution of human cancers microarrays. Oncoimmunology 6, e1284723, doi:10.1080/2162402X.2017.1284723 (2017).

18 Nirmal, A. J. et al. Immune Cell Gene Signatures for Profiling the Microenvironment of Solid Tumors. Cancer Immunol Res 6, 1388-1400, doi:10.1158/2326-6066.CIR-18-0342 (2018).

19 El Zowalaty, M. E. & Jarhult, J. D. From SARS to COVID-19: A previously unknown SARS- related coronavirus (SARS-CoV-2) of pandemic potential infecting humans - Call for a One Health approach. One Health 9, 100124, doi:10.1016/j.onehlt.2020.100124 (2020).

20 Sharifi, N. & Ryan, C. J. Androgen hazards with COVID-19. Endocr Relat Cancer 27, E1-E3, doi:10.1530/ERC-20-0133 (2020).

21 Reading, P. C., Moore, J. B. & Smith, G. L. Steroid hormone synthesis by vaccinia virus suppresses the inflammatory response to infection. J Exp Med 197, 1269-1278, doi:10.1084/jem.20022201 (2003).

22 Raaben, M. et al. The ubiquitin-proteasome system plays an important role during various stages of the coronavirus infection cycle. J Virol 84, 7869-7879, doi:10.1128/JVI.00485-10 (2010).

23 Meier, M. A. et al. Activation of the tryptophan/serotonin pathway is associated with severity and predicts outcomes in pneumonia: results of a long-term cohort study. Clin Chem Lab Med 55, 1060-1069, doi:10.1515/cclm-2016-0912 (2017).

24 Merad, M. & Martin, J. C. Pathological inflammation in patients with COVID-19: a key role for monocytes and macrophages. Nat Rev Immunol, doi:10.1038/s41577-020-0331-4 (2020).

25 Dove, B., Brooks, G., Bicknell, K., Wurm, T. & Hiscox, J. A. Cell cycle perturbations induced by infection with the coronavirus infectious bronchitis virus and their effect on virus replication. J Virol 80, 4147-4156, doi:10.1128/JVI.80.8.4147-4156.2006 (2006).

26 Wu, D. & Yang, X. O. TH17 responses in cytokine storm of COVID-19: An emerging target of JAK2 inhibitor Fedratinib. J Microbiol Immunol Infect, doi:10.1016/j.jmii.2020.03.005 (2020).

27 Ryang, J., Yan, Y., Song, Y., Liu, F. & Ng, T. B. Anti-HIV, antitumor and immunomodulatory activities of paclitaxel from fermentation broth using molecular imprinting technique. AMB Express 9, 194, doi:10.1186/s13568-019-0915-1 (2019).

28 Renis, H. E. Antiviral activity of cytarabine in herpesvirus-infected rats. Antimicrob Agents Chemother 4, 439-444, doi:10.1128/aac.4.4.439 (1973).

29 Subramanian, A. et al. Gene set enrichment analysis: a knowledge-based approach for interpreting genome-wide expression profiles. Proc Natl Acad Sci U S A 102, 15545-15550, doi:10.1073/pnas.0506580102 (2005).

30 Hanzelmann, S., Castelo, R. & Guinney, J. GSVA: gene set variation analysis for microarray and RNA-seq data. BMC Bioinformatics 14, 7, doi:10.1186/1471-2105-14-7 (2013).

Table 1. Literature evidence of drugs that can be repurposed to target SARS-CoV-2-interacting proteins.

Drug	Targets interacting with SARS-CoV-2	Literature/clinical trial evidence
sirolimus	PSMA2,ACSL3,SCARB1,PLOD2,POLA1,XPO1,COMT,BCL2L2,IKBKB,CHEK2,PPIA,GLA,SIGMAR1,IDE	clinical trial: NCT04341675
irbesartan	BCL2,GFER,IKBKB,SIRT5,NARS2	PMID: 32129518
ribavirin	ATP6V1A,PLOD2,BCL2,GGCX,BCL2L2,IMPDH2	PMID: 32227493
apicidin	PSMA2,DNMT1,SCARB1,PLAT,PLOD2,ATF5,PCSK6,COMT,CHEK2,RTN4,SIGMAR1,HDAC2	PMID: 32353859
azacitidine	PSMA2,JUN,SCARB1,ABCC1,PRKAR2B,PCSK6,COMT,CHEK2,F2RL1,BRD2	PMID: 32353859
entacapone	JUN,DNMT1,PRKACA,ABCC1,COMT,F2RL1,TBK1	PMID: 32353859
mercaptopurine	DNMT1,SCARB1,HDAC2,PRKAR2B,XPO1,COMT,PPIA,IDE,NARS2	PMID: 18313035,PMID: 25542975,PMID: 19374142
ruxolitinib	MARK2	PMID: 32353859
midostaurin	MARK2	PMID: 32353859
oseltamivir	NEU1	DOI: https://doi.org/10.1016/S2214-109X(20)30114-5
valproic acid	HDAC2	PMID: 32353859
lisinopril	ACE2	doi: 10.1001/jamacardio.2020.1624
mycophenolic acid	IMPDH2	approved drug
camostat	PRSS1	clinical trial: NCT04353284
silmitasertib	CSNK2A2	PMID: 32353859
metformin hydrochloride	NDUFAF2,NDUFAF1,NDUFB9	PMID: 32360697, PMID: 32347974
fedratinib	BRD4	PMID: 32205092
nafamostat	PRSS1	DOI: 10.1128/AAC.00754-20
gemcitabine hydrochloride	POLA1	https://www.cebm.net/COVID-19/registered-trials-and-analysis/
migalastat	GLA	PMID: 32353859
chloroquine	ACE2	PMID: 32074550
daunorubicin	ABCC1	Chemrxiv: https://orcid.org/0000-0002-4745-366X
ponatinib	RIPK1	PMID: 32353859

Supplementary Dataset 1. Data related to Figure 4a. The list of potential drugs targeting SARS-CoV-2 or SARS-CoV-1 interacting proteins.

Supplementary Dataset 2. Data related to Figure 4b. Drugs that significantly alter the pathways that get affected after COVID-19/SARS infection.

Supplementary Dataset 3. Data related to Figure 5a. Drug and pathway relations, where upregulation (or downregulation) of the pathway is associated with reduced resistance to the drugs in GDSC cell lines.

SupplementaryDataset1.xlsx
The list of potential drugs targeting SARS-CoV-2 or SARS-CoV-1 interacting proteins
SupplementaryDataset2.xlsx
Drugs that significantly alter the pathways that get affected after COVID-19/SARS infection
SupplementaryDataset3.xlsx
Drug and pathway relations, where upregulation (or downregulation) of the pathway is associated with reduced resistance to the drugs in GDSC cell lines.

Download PDF

Journal Publication

published 01 Jan, 2021

Read the published version in Clinical Microbiology and Infectious Diseases →

Version 1

posted

You are reading this latest preprint version

In silico Drug Repurposing to combat COVID-19 based on Pharmacogenomics of Patient Transcriptomic Data

Status:

Journal Publication

Version 1

Abstract

Figures

Introduction

Results

Analysis of transcriptomic data of COVID-19- and SARS-infected patients revealed common altered pathways

Assessing the effect of potential drug candidates on the pathways altered in COVID-19 and SARS infection

Assessing the effect of potential drug candidates on the pathways altered in COVID-19 and SARS infection

Discussion

Methods

Data sets and pre-processing

Computational analysis and metadata

Declarations

References

Tables

Supplementary Datasets

Supplementary Files

Status:

Journal Publication

Version 1