key: cord-0715319-gy9gthag authors: Wu, Qihui; Fan, Xiude; Hong, Honghai; Gu, Yong; Liu, Zhihong; Fang, Shuhuan; Wang, Qi; Cai, Chuipu; Fang, Jiansong title: Comprehensive Assessment of Side Effects in COVID-19 Drug Pipeline from a Network Perspective date: 2020-09-21 journal: Food Chem Toxicol DOI: 10.1016/j.fct.2020.111767 sha: b42de94780f8f85cf240c0b99d15b0a7e24d7744 doc_id: 715319 cord_uid: gy9gthag Currently, coronavirus disease 2019 (COVID-19), has posed an imminent threat to global public health. Although some current therapeutic agents have showed potential prevention or treatment, a growing number of associated adverse events have occurred on patients with COVID-19 in the course of medical treatment. Therefore, a comprehensive assessment of the safety profile of therapeutic agents against COVID-19 is urgently needed. In this study, we proposed a network-based framework to identify the potential side effects of current COVID-19 drugs in clinical trials. We established the associations between 116 COVID-19 drugs and 30 kinds of human tissues based on network proximity and gene-set enrichment analysis (GSEA) approaches. Additionally, we focused on four types of drug-induced toxicities targeting four tissues, including hepatotoxicity, renal toxicity, lung toxicity, and neurotoxicity, and validated our network-based predictions by preclinical and clinical evidence available. Finally, we further performed pharmacovigilance analysis to validate several drug-tissue toxicities via data mining adverse event reporting data, and we identified several new drug-induced side effects without labeling in Food and Drug Administration (FDA) drug instructions. Overall, this study provides forceful approaches to assess potential side effects on COVID-19 drugs, which will be helpful for their safe use in clinical practice and promoting the discovery of antiviral therapeutics against SARS-CoV-2. The current Coronavirus disease 2019 has presented a serious threat to global public health (Del Rio and Malani, 2020) . As of July 18th, 2020, the confirmed cases have reached up to 14 million and caused more than 600,000 deaths . However, many aspects of transmission, infection, and treatment of COVID-19 remain unclear (Wiersinga et al., 2020) . Nowadays, a large number of clinical trials are being conducted on repurposed or experimental drugs to identify prophylactic and therapeutic treatments against COVID-19. Although some current therapeutic agents have showed potential in ameliorating symptoms and reducing mortality in patients (Ledford, 2020) , a growing number of associated adverse events have occurred on patients with COVID-19 in the course of medical treatment. For example, the malaria medication chloroquine, has been brought into multiple clinical trials (e.g. ClinicalTrials.gov number, NCT04353336) to determine its effect on inhibiting severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2). Unfortunately, a latest research has demonstrated that hydroxychloroquine or chloroquine used alone or with a macrolide, do not show a benefit on in-hospital outcomes for COVID-19 and on the contrary, were relevant to decrease in-hospital survival as well as increase frequency of ventricular arrhythmias (Mehra et al., 2020) . Therefore, a comprehensive evaluation of the safety profile of potential therapeutic agents against COVID-19 is of J o u r n a l P r e -p r o o f paramount importance (Javorac et al., 2020) . Generally, conventional experimentally approaches including in vivo and in vitro assays for toxicity determination are time-consuming and costly (Merz et al., 2014) . Thus, computational strategies to assess side effects of therapeutic drugs are propitious to improve efficiency and accelerate the development of drug discovery (Cai et al., 2018; Cai et al., 2019; Wu et al., 2019) . Recent technological and computational advances in genomics and systems biology have offered possibilities for identifying side effects or therapeutic agents by integrating disease proteins and drug-target interactions in human protein-protein interactome Cheng et al., 2019; Fang et al., 2019) . Since drug targets do not work in isolation from the complex system of proteins that comprise the molecular environment of the cell, each drug-target interaction must be examined in an appropriate integrative context (Menche et al., 2015) . Several network-based methodologies, such as network proximity, offer powerful tools for efficient identification of potentially adverse effects or new indications for existing drugs (Fang et al., 2020) . In order to prioritize these drug-side effect pairs identified by network-based approaches, rigorous validation is mandatory . Since approved drugs are already in clinical practice, their potential side effect can be validated via pharmacovigilance analysis that utilizes adverse drug event data available (Böhm et al., 2012) . The Genotype-Tissue Expression (GTEx) is a functional genomics J o u r n a l P r e -p r o o f database that provides genotype-tissue expression to study the relationship between genetic variation and gene expression in human tissues (GTEx project, 2013) . In this study, we developed a network-based framework to identify side effects of the ongoing drugs in COVID-19 pipeline (Figure 1) . We posited that systematic examination of the relationship between drug targets (or drug regulatory genes) and specific expression genes in each tissue from GTEx, will serve as a foundation for generating predictive models to identify potential side effects in tissue. Specifically, we first collected the ongoing drugs in COVID-19 through searching for clinical trials and built their drug-target network from six data sources available as well as drug-gene signature from the Connectivity Map (CMap). To determine the potential adverse events related to COVID-19 drugs, we established the associations between COVID-19 drugs and 30 kinds of human tissues based on network proximity and gene-set enrichment analysis (GSEA) approaches. Moreover, we concentrated on side effects of four tissues including liver, kidney, lung and nerve, and validated several drug-tissue associations utilizing large-scale adverse drug event data and more than half of them supported our network-based prediction. Our overall framework is summarized in The ongoing drugs in COVID-19 pipeline were integrated from current clinical J o u r n a l P r e -p r o o f trials. We first searched for clinical trials (accessed in 21st April, 2020) related to COVID-19 medication by using terms OR SARS-COV-2' in ClinicalTrials.gov database (https://clinicaltrials.gov/). Only trials labeled with "Drug" in "Interventions" column were preserved. In total, we collected 116 drugs from 308 COVID-19 clinical trials after removing duplicates (Table S1 ). Each drug was subsequently matched by the unified DrugBank ID and classified by the Anatomical Therapeutic Chemical Classification (ATC) code. The host proteins of coronavirus (CoV) were collected from four recent public references (Gordon et al., 2020; Maryam Al-Motawa, 2020a; Shen et al., 2020; Zhou et al., 2020) . Zhou et al. has integrated 119 high-quality experimental validated host proteins for drug repurposing for human CoV (Zhou et al., 2020) . Additionally, a latest research has identified 332 high-confidence SARS-CoV-2-human protein-protein interactions (PPIs) which were connected to multiple biological processes (Gordon et al., 2020) . Moreover, 204 and 50 host proteins of SARS-CoV-2 have been separately acquired from two recent proteomics publication (Maryam Al-Motawa, 2020a; Shen et al., 2020) . In total, 705 CoV-related host proteins were consequently integrated (Table S2 ). We assembled high-quality physical drug target interactions (DTIs) for ongoing COVID-19 drugs from six authoritative databases, including DrugBank J o u r n a l P r e -p r o o f database (v4.3) (Wishart et al., 2018) , the Therapeutic Target Database (TTD, v4.3.02) (Yan et al., 2018) , the PharmGKB database (Barbarino et al., 2018) , ChEMBL (v20) (Gaulton et al., 2017) , BindingDB (Gilson et al., 2016) , and IUPHAR/BPS Guide to PHARMACOLOGY (Harding et al., 2018) The RNA-seq data (RPKM value) of 30 tissues were extracted from the GTEx V6 release (https://gtexportal.org/home/). We asserted that a gene is expressed in a specific tissue if its RPKM value is larger or equal to 1 in over 80 % of samples, and vice versa. To quantify the expression significance of a tissue-expressed gene i in a tissue t, the average expression E i and the standard deviation δ_E i of a gene's expression across all considered tissues were calculated. The significance of gene expression in tissue t is defined as below: Z_E i, t = E i, t − E i /δ_E i J o u r n a l P r e -p r o o f The gene was considered as high expressed in a specific tissue if its Z-expression was higher than 2 in this tissue. Given A, the module formed by specific high expressed proteins in a given tissue (tissue module), and B, the module formed by target proteins of a given drug (drug module), the closest distance represents the average shortest path length of all the nodes in tissue module (A) to drug module (B) in the human protein-protein interactome, which can be described as: ", # indicates the shortest path length between tissue protein " and drug target #. To assess the significance of the network distance between each drug module and tissue module, we constructed a reference distance distribution corresponded to the expected distance of each two randomly selected groups of proteins from human protein-protein interactome which have the same size and degree distribution as the original protein modules. This process was repeated 1000 times. According to the mean and standard deviation (% & ) of the reference distribution, the Z-score (' & ) was computed by converting an observed (non-Euclidean) distance to a normalized distance: The gene set enrichment analysis (GSEA) was performed as an additional prioritization method. Drug-gene signatures were collected from the Connectivity Map (CMap2.0), a pharmacogenomics database inducing over 7000 gene expression profiles from cultured human cell lines treated with 1,309 bioactive molecules at different concentrations, covering 6100 individual instances (Lamb et al., 2006) . CMap offers a measure to evaluate the extent of differential expression for a given probe set. The amplitude (") is defined as follows: where ) refers to the scaled and threshold average difference value for the drug treatment group, and * represents the threshold average difference value for the control group. Therefore, " > 0 and " < 0 indicate increased expression (upregulation) and decreased expression (downregulation) upon treatment respectively, while " = 0 means no differential expression. In this study, drug gene signatures with absolute value of " greater than two-fold induction (|"| > 0.67) were defined as drug-gene pairs. Finally, we obtained the drug regulatory gene network for 38 COVID-19 ongoing drugs, covering 91,164 drug-gene pairs. (Table S4) . We assumed that a drug would have a correlation with a specific tissue if its regulatory genes tend to be the specific high expressed proteins in this tissue. For each drug-tissue pair, we calculated the number of proteins in each tissue-specific protein network, drug regulatory gene network, as well as their overlap. The statistical significance (P-values) of the gene-set enrichment was calculated by Fisher's exact test and corrected by Benjamini-Hochberg J o u r n a l P r e -p r o o f method (Benjamini and Hochberg, 1995) . We set corrected P-value (q) < 0.05 as a threshold to define significantly predicted drug-tissue pairs. We retrieved adult adverse drug event data from the US Food and Drug Administration's (FDA) Adverse Event Reporting System (AERS) compiled from the first quarter of 2010 to the last quarter of 2019 using OpenVigil 2.1 (Böhm et al., 2012) (http://openvigil.sourceforge.net), which is a platform for data mining and analysis of adverse drug event data. In the AERS, adverse events can be specified at different levels of the medical dictionary for regulatory activities (MedDRA) terminology. We employed 9 preferred terms (PTs) for collecting relevant cases associated with "drug-induced liver injury", 7 PTs for cases associated with "drug-induced kidney injury", 4 PTs for cases associated with "drug-induced lung injury", and 4 PTs for cases associated with "drug-induced neuropathy". Detailed PTs and their number for each drug-induced organ injury were listed in the Table 1 . Moreover, the major advantage of the OpenVigil 2.1 is that the platform could perform Yates' chi-square test and calculate disproportionality (like reporting odds ratio (ROR), proportional reporting ratios (PRRs), and relative reporting ratio (RRR)) from adverse drug event data to determine the relationship between the drug and adverse event. In this study, number of adverse events of drugs >4, Chi square values >4, and the lower limit of the 95% confidence interval (CI) of the ROR>1 indicate that the adverse event and the drug is related and can be considered as "likely an adverse reaction". The RORs and 95% CIs for drugs and adverse events of interest were shown using the forest plots. J o u r n a l P r e -p r o o f The statistical analysis in this study was performed by Python (v3.2, http://www.python.org/) and R platforms (v3.01, http://www.r-project.org/). Networks were generated by Cytoscape (v3.2.1, http://www.cytoscape.org/) and Gephi (v0.9.2, https://gephi.org/). In the graphical network, drugs or genes were presented by nodes and interactions were encoded by edges. The degree of each node was calculated to evaluate its topological property, as this characterizes the most important nodes in a network. As of April 21, 2020, there were 116 agents in 308 clinical trials of anti-SARS-CoV-2 therapies. Here, a comprehensive overview of available information on potential drugs for treatment of COVID-19 was presented ( Figure 2) . Most of ongoing drugs in COVID-19 clinical trials are approved drugs since existing drug could shorten the time and reduce the cost in comparison with developing novel compounds from scratch (Kupferschmidt and Cohen, 2020) . In addition, there are still 16 investigational agents (e.g. remdesivir) and 4 experimental agents (e.g. tetrandrine) enrolled into the clinical trials for their potential therapeutical effects on COVID-19. As shown in and inhibiting the inflammatory processes) (Wiersinga et al., 2020) . As recent study demonstrated that chemotherapy could produce selective toxicity to SARS-CoV-2 via improving the levels of endogenous reactive metabolites (Maryam Al-Motawa, 2020b), scientists have turned their attention to antineoplastic agents. Besides, coronavirus was able to decrease immune capacity since it primarily exerts cytotoxic action on the hematopoietic and immune systems with reducing the number of neutrophils (Remuzzi and Remuzzi, 2020) . Thus, immunomodulating agents were also brought into the clinical trials in view of its positive regulation on immunity of the organism. Furthermore, Figure 2 shows that 16 drugs intend to be classified as antiinfectives for systemic use. Since SARS-CoV-2 is a single-stranded RNA-enveloped virus that can be infected from person to person, drugs for antiinfectives for systemic use such as azithromycin (Parra-Lara et al., 2020) are also considered to own the potential of suppressing the virus. Among the 116 ongoing agents in pipeline, 76% of them (88 agents) are being tested in phase 2 and phase 3 trials with uncertain effects rather than confirmed therapeutic effects on COVID-19. Although a few agents in COVID-19 trials (Figure 2 ) seem to show promising results based on current researches and clinical trials, they may play a limited role in COVID-19 treatment or even remain controversial (Ko et al., 2020) . Meanwhile, it also needs to be alerted for the potential toxicity of these drugs. We constructed a comprehensive drug-target network of ongoing drugs in COVID-19 pipeline through integrating curated physical DTIs and CoV host proteins. As depicted in Figure 3 In this study, the expression level of genes across 30 human tissues is quantified by Z-expression, including colon, esophagus, fallopian tube, heart, kidney, liver, lung, muscle, nerve, ovary, adipose tissue, pancreas, pituitary, prostate, salivary, skin, small intestine, spleen, stomach, testis, thyroid, adrenal gland, uterus, vagina, bladder, blood vessel, blood, brain, breast, and cervix uteri. Genes in a tissue-protein network with Z-expression greater than 2 are considered as high expression. As shown in Figure 4A , testis tissue possesses the largest number of high expressed genes (n=3,521), followed by brain (n=733) and liver (n=465). We further prioritized the high-expressed genes in each tissue by Z-expression and highlighted the top 10 specific genes of the 30 human tissues ( Figure 4B) . Interestingly, we find that among these genes, there are no overlap between any two tissues, implying they tend to have a high multiplicity and tissue specificity. These high-expressed genes seem to possess high association with the occurrence of drug-induced tissue toxicity. To identify the potential tissue toxicities of COVID-19 ongoing drugs, we applied two complementary approaches to compute the significance of the association between each tissue and drug, including: (i) network proximity approach which quantifies the network distance between each drug module and tissue module in human protein-protein interactome; and (ii) GSEA approach that calculates the statistical significance of the gene-set enrichment between each tissue-specific protein network and drug regulatory gene network. We herein assumed that a drug may induce potential toxicity on specific tissue if they have significant associations (P network proximity or q GESA < 0.05). As illustrated in Figure 5A , there are in total 42 drugs significantly associated with multiple tissues via network proximity approach (P network proximity < 0.05), covering 140 significant drug-tissue pairs. It suggests that those drugs may cause corresponding predicted tissue toxicities. For example, deferoxamine, which targets kidney organ and exhibits great tissue specificity (Z-score = -2.01), has been demonstrated the potential risk of inducing acute renal failure (Prasannan et al., 2003) . Moreover, deferoxamine also has significant association with liver tissue (Z-score = -1.97). In vitro studies have suggested that deferoxamine was toxic to hepatoma cells in the absence of J o u r n a l P r e -p r o o f with different human tissues, including 132 significant drug-tissue pairs. Taking sirolimus (also known as rapamycin) as an example, GSEA gave a q GESA value of 0.003 with lung tissue, indicating its potential lung toxicity. An in vivo study has revealed that sirolimus could augment lipopolysaccharide-induced lung injury and apoptosis in C57BL/6J mice through involving signal transducer and activator of transcription 1 (STAT1) and the induction of STAT1-dependent apoptosis genes (Fielhaber et al., 2012) . Notably, we found that 9 drug-tissue pairs were simultaneously identified by network proximity and GSEA approaches, alerting their potential tissue toxicities, including captopril-liver, sirolimus-blood vessel, atovaquone-adrenal gland, deferoxamine-blood, doxycycline-blood vessel, deferoxamine-testis, deferoxamine-spleen, deferoxamine-blood vessel, and captopril-ovary. The detailed information of the network-based prediction results are provided in Table S5 . J o u r n a l P r e -p r o o f To present specifically how drug candidate manifests itself and targets tissue toxicity within the drug-target networks and gene expression profiles, we further focus on adverse events on four tissues including lung, nerve, liver and kidney. As exhibited in Figure 6 , there are 42 drug-tissue pairs covering 32 agents predicted by network proximity or GSEA approaches. Among them, 64.3% (27/42) of the predicted tissue toxicities were validated by drug labels retrieved from Drug@FDA database, clinical evidences including clinical studies and case reports, as well as preclinical evidence available (Table S6) . Given the significance of clinical data, we further assessed whether these drug-tissue associations mentioned above were of statistical difference regarding the side effects of drug in specific tissue via pharmacovigilance analysis. Among 42 predicted significant drug-tissue pairs above, we selected J o u r n a l P r e -p r o o f drug-tissue pairs for further validation using subject matter expertise based on a combination of factors: (i) not being labeled by FDA drug instructions; (ii) strength of the network-based prediction; and (iv) availability of supportive experimental evidence. Applying these criteria resulted in identifying 15 potential drug-tissue toxicities across four tissues for further validation ( Table 1 ). Specifically, we first retrieved a total of 3,004,759 adverse events from the AERS. We identified 25,443, 60,024, 16,808, and 114,264 cases who developed drug-induced liver injury, drug-induced kidney injury, drug-induced lung injury, and drug-induced neuropathy while receiving any drugs, respectively. The statistical significant drug-induced toxicity across four tissues together with their RORs and 95% CIs were shown in Figure 7 and the detailed information on all drug-tissue associations can be found in Figure S1 . According to the criteria, drug-induced liver injury was significantly related to ibuprofen, enoxaparin, cholecalciferol, folic acid, and sirolimus. Indeed, recent literatures have reported that these drugs do cause liver injury, such as sirolimus (Neff et al., 2004; Umemura et al., 2014) . Surprisingly, we found some drug-tissue associations which had not been labeled by FDA drug instructions were of statistical difference, such as the potential hepatotoxicity of enoxaparin (Hahn et al., 2015; Pivarnik et al., 2016) . According to the pharmacovigilance analysis of adverse drug event data, enoxaparin seems to possess significant association with hepatotoxicity (ROR=1.312, 95% The current 2019-nCoV outbreak is still moving fast, and a rebound of victims of COVID-19 in some areas become a more serious concern . Athough numerous drugs are under therapeutic evaluation in clinical trials toward SARS-CoV-2, several medications such as hydroxychloroquine and azithromycin have been reported to exacerbate cardiotoxicity and liver injury in COVID-19 patients (Boeckmans et al., 2020) . Additionally, the pharmacological and toxicity properties of some investigational or experimental drugs in clinical trials have not been fully explored, meaning that there is greater safety risk within these drugs. Therefore, it is necessary to J o u r n a l P r e -p r o o f comprehensively assess side effects of ongoing COVID-19 drugs, which could benefit for the discovery and clinical application of pharmacotherapy for COVID-19. In this study, we utilized network-based approaches for rapid identification of potential side effects for ongoing COVID-19 drugs. Specifically, after integration of drug-target network, drug-gene signatures, tissue-specific gene expression profiles from GTEx, and human protein-protein interactome, we established the associations between 116 COVID-19 drugs and 30 kinds of human tissues based on network proximity and GSEA approaches. Our network proximity and GSEA approaches achieved a high accuracy to identify the associations between COVID-19 drugs and human tissues. For example, ibuprofen was deemed to own potential toxicity on liver based on our network-based prediction (Figure 6) . Emerging evidence shows that ibuprofen may elicit hepatotoxicity via altering several major pathways including energy metabolism, protein degradation, fatty acid metabolism as well as antioxidant system in mice (Tiwari et al., 2020) . Additionally, an in vivo study (Wen et al., 2018) also revealed that ibuprofen treatment on rats increased the levels of alanine transaminase (ALT) and aspartate aminotransferase (AST), causing more serious liver injury than did aspirin. Moreover, we found several drugs could induce potential toxicities in multiple tissues. Taking and acute renal failure (Tan et al., 2011) induced by captopril to some extent verified our predictions. Interestingly, we found that captopril-liver pair was simultaneously identified by network proximity and GSEA methods, further indicating its high probability of liver injury. In fact, Drug@FDA files on captopril have clearly labled its potential hepatoxicity. Moreover, we found several potential tissue toxicities of drugs which hadn't been labeled by FDA drug However, several limitations should be acknowledged in present study. First, both network proximity and GSEA approaches were applied to predict tissue toxicity of drugs by integrating physical drug-target interactions and transcriptome data on CMap. Given lacking of large scale published available experimental data, whether combining the two methods (network proximity and GSEA) synergistically enhance the prediction performance remain further evaluation. Moreover, this study mainly focuses on the prediction and validation of drug-tissue associations, the specific mechanism of actions (MOAs) related with tissue toxicity need to be explored in future. In addition, adverse drug events obtained through the adverse reaction reporting system should be treated with caution, since there are inherent limitations, such as selective reporting, under-reporting, and lack of information about total drug consumption. The results of disproportionality analyses using AERS data should not be considered as the quantification of the true risk, but rather as the assessment of signal detection. Finally, as current COVID-19 medications are being renewed very quickly, we can not cover all the ongoing drugs in COVID-19 pipeline. In summary, this study provides a powerful network-based methodology J o u r n a l P r e -p r o o f for efficient identification of drug-induced potential side effects in COVID-19 drug pipeline, followed by pharmacovigilance analysis that utilized large scale clinical adverse drug event data. If broadly applied, our computational framework can be used to assess potential tissue toxicity for all drugs, such as cancer or cardiovascular drugs. In this manner, we can build bridges between basic medical research and clinical outcomes, which could shorten the time to prioritize side effects associated with drugs and alert the potential adverse effects of drugs during pharmacotherapy course of COVID-19. The authors declare that they have no known competing financial interests or personal relationships that could have appeared to influence the work reported in this paper. OpenVigil--free eyeballs on AERS J o u r n a l P r e -p r o o f pharmacovigilance data PharmGKB: A worldwide resource for pharmacogenomic information Controlling the False Discovery Rate: A Practical and Powerful Approach to Multiple Testing Imatinib might constitute a treatment option for lung involvement in COVID-19 COVID-19 and drug-induced liver injury: a problem of plenty or a petty point? In Silico Pharmacoepidemiologic Evaluation of Drug-Induced Cardiovascular Complications Using Combined Classifiers Deep Learning-Based Prediction of Drug-Induced Cardiotoxicity Network-based approach to prediction and population-based validation of in silico A genome-wide positioning systems network algorithm for in silico drug repurposing Deferoxamine toxicity in hepatoma and primary rat cortical brain cultures Novel Coronavirus-Important Information for Clinicians Case reports study of the first five patients COVID-19 treated with remdesivir in France Quantitative and systems pharmacology 4. Network-based analysis of drug pleiotropy on coronary artery disease Harnessing endophenotypes and network medicine for Alzheimer's drug repurposing Inhibition of mammalian target of rapamycin augments lipopolysaccharide-induced lung injury and apoptosis BindingDB in 2015: A public database for medicinal chemistry, computational chemistry and systems pharmacology The potential of JAK/STAT pathway inhibition by ruxolitinib in the treatment of COVID-19 A SARS-CoV-2 protein interaction map reveals targets for drug repurposing Enoxaparin-Induced Liver Injury: Case Report and Review of the Literature and FDA Adverse Event Reporting System (FAERS). Drug Saf Case Rep The IUPHAR/BPS Guide to PHARMACOLOGY in 2018: updates and expansion to encompass the new guide to IMMUNOPHARMACOLOGY An overview of safety assessment of the medicines currently used in the treatment of COVID-19 disease Uromodulin p.Cys147Trp mutation drives kidney disease by activating ER stress and apoptosis Arguments in favour of remdesivir for treating SARS-CoV-2 infections Race to find COVID-19 treatments accelerates The Connectivity Map: using gene-expression signatures to connect small molecules, genes, and disease Coronavirus breakthrough: dexamethasone is first drug shown to save lives Captopril. LiverTox: Clinical and Research Information on Drug-Induced Liver Injury Fluvoxamine Virus-host interactome and proteomic survey of PBMCs from COVID-19 patients reveal potential virulence factors influencing SARS-CoV-2 pathogenesis Vulnerabilities of the SARS-CoV-2 virus to proteotoxicityopportunity for repurposed chemotherapy of COVID-19 infection RETRACTED: of COVID-19: a multinational registry analysis Disease networks. Uncovering disease-disease relationships through the incomplete interactome Methodology to Assess Clinical Liver Safety Data Sirolimus-associated hepatotoxicity in liver transplantation Azithromycin and SARS-CoV-2 infection: where we are now and where we are going Enoxaparin-induced hepatotoxicity: an under-recognised complication of enoxaparin therapy Acute renal failure following deferoxamine overdose The Genotype-Tissue Expression (GTEx) project COVID-19 and Italy: what next? Molecular characterization of an aquaporin-2 mutation causing a severe form of nephrogenic diabetes insipidus Proteomic and Metabolomic Characterization of COVID-19 Patient Sera Serum uromodulin is associated with the severity of clinicopathological findings in ANCA-associated glomerulonephritis Captopril induced reversible acute renal failure in a premature neonate with double outlet right ventricle and congestive heart failure Gender-specific changes in energy metabolism and protein degradation as major pathways affected in livers of mice treated with ibuprofen Early involvement of cellular stress and inflammatory signals in the pathogenesis of tubulointerstitial kidney disease due to UMOD mutations Liver damage, inflammation, and enhanced tumorigenesis after persistent mTORC1 inhibition Metabolism of liver CYP450 and ultrastructural changes after long-term administration of J o u r n a l P r e -p r o o f aspirin and ibuprofen Pathophysiology, Transmission, Diagnosis, and Treatment of Coronavirus Disease Identification and Mechanism Exploration of Hepatotoxic Ingredients in Traditional Chinese Medicine Glycycoumarin protects mice against acetaminophen-induced liver injury predominantly via activating sustained autophagy Network-based drug repurposing for novel coronavirus 2019-nCoV/SARS-CoV-2 This work was supported by the National Natural Science Foundation of China