key: cord-0885829-9l3x3owc authors: Creeden, Justin Fortune; Imami, Ali Sajid; Eby, Hunter M.; Gillman, Cassidy; Becker, Kathryn N.; Reigle, Jim; Andari, Elissar; Pan, Zhixing K; O’Donovan, Sinead M; McCullumsmith, Robert E; McCullumsmith, Cheryl B title: Fluoxetine as an Anti-Inflammatory Therapy in SARS-CoV-2 Infection date: 2021-02-25 journal: Biomed Pharmacother DOI: 10.1016/j.biopha.2021.111437 sha: ae9380f255867a42d8a67e2bfa077346005c7f05 doc_id: 885829 cord_uid: 9l3x3owc Hyperinflammatory response caused by infections such as Severe Acute Respiratory Syndrome Coronavirus 2 (SARS-CoV-2) increases organ failure, intensive care unit admission, and mortality. Cytokine storm in patients with Coronavirus Disease 2019 (COVID-19) drives this pattern of poor clinical outcomes and is dependent upon the activity of the transcription factor complex nuclear factor kappa-light-chain-enhancer of activated B cells (NF-kappaB) and its downstream target gene interleukin 6 (IL6) which interacts with IL6 receptor (IL6R) and the IL6 signal transduction protein (IL6ST or gp130) to regulate intracellular inflammatory pathways. In this study, we compare transcriptomic signatures from a variety of drug-treated or genetically suppressed (i.e. knockdown) cell lines in order to identify a mechanism by which antidepressants such as fluoxetine demonstrate non-serotonergic, anti-inflammatory effects. Our results demonstrate a critical role for IL6ST and NF-kappaB Subunit 1 (NFKB1) in fluoxetine’s ability to act as a potential therapy for hyperinflammatory states such as asthma, sepsis, and COVID-19. disease outcomes [1] . Clinical consequences of COVID-19 infection such as multisystem organ failure, intensive care unit admission, and death result from an exaggerated immune response whereby increased IL6 drives a destructive cytokine storm [2] [3] [4] [5] [6] . This pattern of poor clinical J o u r n a l P r e -p r o o f Page 4 outcomes parallels the Secondary Hemophagocytic Lymphohistiocytosis (SHLH) seen in another coronavirus-induced illness, Severe Acute Respiratory Syndrome (SARS) [7, 8] . Indeed, many infections result in an IL6-mediated cytokine storm that drives negative clinical outcomes. NF-kappaB-mediates increases in the transcriptional expression of the proinflammatory cytokine IL6 [9, 10] . Binding of IL6 to its receptor complex (IL6R) activates the IL6 signal transduction protein (IL6ST or gp130), a subunit of the IL6 receptor complex. IL6ST activation is integral to many intracellular cytokine-mediated inflammatory pathways, including the cytokine storm / cytokine release syndrome in COVID-19 patients [11] . Further, in vitro, in vivo, and clinical data demonstrate that many selective serotonin reuptake inhibitors (SSRIs) decrease IL6 signaling activity and reduce hyperinflammatory states [12] [13] [14] [15] [16] [17] . The SSRI fluoxetine causes pronounced inhibition of NF-kappaB signaling and IL6 expression correlating with decreased illness severity across a variety of animal disease models including cancer, pulmonary inflammation, and sepsis [17] [18] [19] [20] [21] [22] . Fluoxetine's well-established safety profile and its inhibition of IL6 signal transduction suggests that this drug may be advantageous when used early in COVID-19 infection to prevent or reduce the cytokine storm. Evidence supporting this hypothesis can be found across multiple animal models of disease whereby fluoxetine not only reduced IL6 signaling cascades but also prevented end organ damage [16, 18, [23] [24] [25] . In rat models of chronic obstructive pulmonary disease, fluoxetine decreased lung injury by inhibiting the transcription factor NF-kappaB, thus reducing downstream IL6 production [18] . In a lipopolysaccharide (LPS)-induced sepsis model of hyperinflammatory illness, pretreatment with fluoxetine outperformed pretreatment with corticosteroids and resulted in diminished end organ damage that included decreased incidence J o u r n a l P r e -p r o o f Page 5 of pulmonary arterial hypertension, decreased pulmonary arterial muscularization, and decreased extracellular matrix remodeling [23, 24] . Fluoxetine treatment in rats decreased bronchial asthma severity concurrent with reduced levels of inflammatory cytokines, including IL6 [24] . Fluoxetine has also demonstrated inhibition of NF-kappaB with corresponding decreases in inflammatory cytokines like IL6 in models of septic shock and allergic asthma, as well as in patients with depression [16, 25] . Recently, Zhou et al. performed phylogenetic analyses of multiple human coronavirus genomes to identify evolutionarily conserved targets for SARS-CoV-2 therapeutics. Their quantification of interactions between the human coronavirus-host interactome and subsequent gene set enrichment analyses comparing coronavirus-infected host cells and drug-treated human cells identified the SSRI paroxetine as a novel COVID-19 therapeutic candidate [26] . An observational study correlated antidepressant therapy with decreased risk of intubation or death for hospitalized patients with COVID-19 [27] . Preliminary reports also suggest that fluoxetine demonstrates direct inhibition of SARS-CoV-2 viral protein expression [28] , although the present study is focused on fluoxetine's anti-inflammatory activity in human pathophysiology. To investigate the potential clinical utility of using fluoxetine to reduce IL6 mediated inflammation in COVID-19 patients and to elucidate the mechanism by which fluoxetine acts as an anti-inflammatory agent, we used transcriptomic signatures from the Library of Integrated involved in fluoxetine's anti-inflammatory mechanism of action were compared to the transcriptomic changes caused by antidepressant drug therapy. Our resultant data highlights the significance of IL6ST and NF-kappaB in fluoxetine's anti-inflammatory mechanism of action and supports further exploration of fluoxetine as a candidate drug in the treatment of inflammatory pathologies such as those contributing to COVID-19. assay which measures the expression of 978 "hub" genes which represent ~82% of the information present in the transcriptome. The iLINCS Portal (www.ilincs.org) provides a convenient gateway to access these signatures. We used the iLINCS application programing interface (API) to download differential gene expression signatures. We acquired all available signatures for fluoxetine (33) , paroxetine (65), bupropion (29) , or dexamethasone (309) treated cell lines. In total, 436 signatures were obtained for all 4 drug treatments, with data generated in 45 different cell lines, at different timepoints (6h -24h) and different doses (0.01uM -80uM) (Table S1 ). To create high-confidence gene signatures, differential gene expression data were filtered such that only those genes with a log-fold change (LFC) ≥ 0.85 or ≤ -0.85 in drug-treated versus corresponding control cell lines were included in downstream gene signature analyses. This LFC threshold was found to effectively reduce excess noise from low-difference genes and Page 7 reduce signatures to their most significant components [29, 30] . The resulting signatures were searched against the iLINCS Consensus Gene Knockdown (CGS) database to generate similarity indices between (a) the drug treatment signatures within a given cell line and (b) individual inflammatory gene knockdown signatures for 27 inflammatory genes with a potential role in fluoxetine's anti-inflammatory mechanism of action (i.e. drug treatment gene signatures were compared to gene knockdown signatures generated only in the same cell line). A concordance score is an adjusted correlation coefficient that can range from -1 to 1. iLINCS only reports a match if the concordance score is ≤ -0.2 or ≥ 0.2. Concordance scores outside of those parameters are considered statistically insignificant. In this way, we identified 779 total signatures (concordance ≤ -0.2 or ≥ 0.2) derived from comparisons of the 168 drug treatment gene signatures to CGS for 27 inflammatory genes in 11 different cell lines (Table S2) . We again filtered our list to include only the highest absolute concordance score (range -0.7166 -0.7982) for each drug treatment, knockdown, and cell line combination. This resulted in 395 concordance scores comparing 27 inflammatory genes and 4 drugs of interest across 11 cell lines (Table S3) . Mean concordance scores are calculated by averaging the concordance scores for a given drug-gene comparison across cell lines. Statistically significant outliers (P < 0.5) within the list of drug-gene concordance scores among cell lines containing calculated concordance scores for less than 3 drugs of interest were identified using the extreme studentized deviate method (Grubbs' test) and subsequently excluded. The entire process was performed using R version 3.6.3. The entire bioinformatic pipeline including raw data, analytical scripts, and results can be reproduced using our publicly available J o u r n a l P r e -p r o o f To assess the potential clinical utility of fluoxetine in the prevention or mitigation of cytokine storm syndromes, and to investigate the mechanism by which fluoxetine realizes its anti-inflammatory effect, we generated differential gene expression signatures produced by fluoxetine, paroxetine, bupropion, or dexamethasone drug treatments. We also generated differential gene expression signatures produced by genetic knockdown of inflammatory genes related to the initiation and maintenance of IL6-mediated cytokine storm. We then compared the drug treatment signatures with the genetic knockdown signatures to generate concordance scores which measure the degree to which a given drug treatment parallels genetic knockdown of a given inflammatory gene. When a drug treatment signature is highly concordant with a genetic knockdown signature, we surmise that this drug and this gene knockdown induce a similar pattern of gene expression changes [31] . We compared gene signatures generated from cell lines that have been treated with fluoxetine, paroxetine, bupropion, and dexamethasone to matched cell lines in which IL6ST was genetically knocked down with siRNA ( Figure 2A ). In this way, the differential gene expression signatures for drug-treated A375 or HEPG2 cell lines were compared to IL6ST knockdown A375 or HEPG2 cell lines, respectively. This allowed us to quantify the degree of similarity between pharmacologically initiated mRNA changes and those resulting from targeted knockdown of a specific gene. In A375 cells, concordance scores quantifying the similarity between drug treated To further develop the relationship between our drugs of interest, IL6ST, and inflammatory pathways, we again analyzed differential gene expression signatures caused by drug treatments in A375 and HEPG2 cell lines. This time however, we compared drug treatment signatures to those produced by NFKB1 knockdown ( Figure 2B ). As before, comparisons of these signatures generated concordance scores quantifying the degree of similarity between a drug treatment and NFKB1 knockdown. In A375 cells, these concordance scores are as follows: These data are particularly interesting in contrast to dexamethasone-induced transcriptional changes. Surprisingly, treatments with fluoxetine approximate inflammatory IL6ST and NFBK1 gene knockdowns more significantly than dexamethasone-a steroid routinely used to treat cytokine induced hyperinflammation. Administration of dexamethasone in porcine respiratory coronavirus-infected pigs downregulated IL6 in the early disease stages suggesting that corticosteroids such as dexamethasone may be beneficial to COVID-19 patients if administered during the early acute phase of infection [34, 35] . A study comparing over 2,000 patients being treated with dexamethasone to over 4,000 control patients receiving usual care showed that hospitalized patients receiving respiratory support such as mechanical ventilation or Page 13 oxygen demonstrate lower 28-day mortality when treated with dexamethasone [36] . Concordance scores generated by comparisons of fluoxetine treatment and IL6ST knockdown or NFKB1 knockdown are also distinctly higher than those comparing bupropion or paroxetine treatment with the same inflammatory gene knockdowns. This is striking given that both bupropion and paroxetine share monoamine neurotransmission profiles with fluoxetine. Previous reports describe the NF-kappaB family of transcription factors as regulating many proinflammatory target genes including IL6ST [9, 10, 37] . Taken together, our data suggest that fluoxetine induces a reduction in IL6ST-mediated signal transduction via NF-kappaB-mediated transcriptional changes. These data also suggest that the anti-inflammatory action of fluoxetine is independent of its action on monoaminergic pathways. Instead, fluoxetine may manifest its antiinflammatory character through a mechanism of action related to the suppression of inflammatory genes such as IL6ST and NFKB1. IL6 signaling is enormously complex and capable of eliciting both pro-and antiinflammatory cellular behaviors [38, 39] . Classical IL6 pathways reliant upon membrane-bound IL6R tend to elicit anti-inflammatory responses. In contrast, trans-signaling-mediated responses are distinctly pro-inflammatory and proceed by way of IL6ST on cells that do not express membrane-bound IL6R [40] . These trans-signaling-mediated responses require soluble IL6R (sIL6R) to interact with IL6 and IL6ST proteins (Figure 3) . A fluoxetine-mediated reduction in IL6ST signal transduction could therefore specifically decrease hyper-inflammatory IL6 pathways while maintaining the classic anti-inflammatory IL6 pathways critical to the reestablishment of homeostasis in hyperinflammatory disorders [41] . Indeed, our analysis shows Page 14 insignificant similarity between fluoxetine treated cells and IL6R knockdown cells (concordance score = 0.2467; supplementary tables 2 and 3). Recent work from our group demonstrates antidepressant-mediated inhibition of IL6 and NF-KappaB signaling pathways limit organ damage, decrease pro-inflammatory cytokine production, decrease intracellular migration of early-stage inflammatory response, and improve animal survival after an overwhelming systemic inflammatory response [42] . The present study further supports the potential of antidepressants like fluoxetine to prevent cytokine storm severity and reduce the risk of severe organ dysfunction and death for COVID-19 patients. Indeed, preliminary reports from a multi-site fluoxetine clinical trial in France [27] , as well as our group's ongoing fluoxetine clinical trial in America (NCT04377308) provide encouraging results suggesting fluoxetine treatment significantly reduces morbidity and mortality in COVID-19 patients. Increased IL6 signal transduction causes cytokine storm which is associated with increased morbidity and mortality in COVID-19 patients. SARS-CoV-2 vaccine distribution is underway but immediate life-saving drug treatments are necessary to decrease case-fatality rates. NF-kappaB plays a pivotal role in cytokine storms through transcriptional activation of the proinflammatory cytokines IL6 and IL6ST. This study provides further evidence supporting the use of fluoxetine to decrease NF-kappaB signaling and thereby decrease the IL6ST signal transduction pathway driving cytokine storm in SARS-CoV-2 infection. Using fluoxetine to disrupt this NF-kappaB/IL6ST axis and mitigate the resultant cytokine storm may increase survival and decrease rates of hospitalization in COVID-19 patients. We generated differential gene expression signatures produced by fluoxetine, paroxetine, bupropion, or dexamethasone drug treatments across a variety of concentrations and treatment durations compared to control, untreated, wild-type cells. We also generated differential gene expression signatures produced by genetic knockdown of inflammatory genes related to the initiation and maintenance of IL6-mediated cytokine storm. Again, these signatures are in relation to control, untreated, wild-type cells. We then compared the drug treatment signatures with the genetic knockdown signatures to quantify similarity using concordance scores. When a drug treatment signature is highly concordant with a genetic knockdown signature, it generates a value approaching +1 and allows us to surmise that this drug and this genetic knockdown induce equivalent changes in gene expression. When no significant similarities exist, the concordance score approaches zero. When a drug treatment signature is highly discordant with a genetic Page 16 knockdown signature, it generates a value approaching -1 and allows us to surmise that this drug and this genetic knockdown induce inverted changes in gene expression. The transcription factor NF-kappaB binds the IL6 gene to induce transcriptional expression of the proinflammatory cytokine. Pro-inflammatory trans-signaling-mediated responses involve IL6 proteins binding to soluble IL6 receptors (sIL6R) which activates the IL6 signal transduction protein (IL6ST or gp130). IL6ST activation is integral to many intracellular cytokine-mediated inflammatory pathways including the cytokine storm / cytokine release syndrome driving many of the negative clinic outcomes associated with COVID-19. Immunotherapeutic implications of IL-6 blockade for cytokine storm Epidemiological and clinical characteristics of 99 cases of 2019 novel coronavirus pneumonia in Wuhan, China: a descriptive study Clinical features of patients infected with 2019 novel coronavirus in Wuhan COVID-19: consider cytokine storm syndromes and immunosuppression Clinical predictors of mortality due to COVID-19 based on an analysis of data of 150 patients from Wuhan, China Complex Immune Dysregulation in COVID-19 Patients with Severe Respiratory Failure Early enhanced expression of interferon-inducible protein-10 (CXCL-10) and other chemokines predicts adverse outcome in severe acute respiratory syndrome Dynamic changes in blood cytokine levels as clinical indicators in severe acute respiratory syndrome Activation of interleukin-6 gene expression through the NFkappa B transcription factor NF-kappaB signaling in inflammation Cytokine release syndrome in severe COVID-19: interleukin-6 receptor antagonist tocilizumab may be the key to reduce mortality Sertraline decreases serum level of interleukin-6 (IL-6) in hemodialysis patients with depression: results of a randomized double-blind, placebo-controlled clinical trial Selective serotonin reuptake inhibitors as a novel class of immunosuppressants IL-6 levels decrease with SSRI treatment in patients with major depression Paroxetine differentially modulates LPS-induced TNFalpha and IL-6 production in mouse macrophages Effects of SSRIs on peripheral inflammatory markers in patients with major depressive disorder: A systematic review and meta-analysis Antidepressants inhibit interferon-gamma-induced microglial production of IL-6 and nitric oxide Fluoxetine, an Antidepressant Drug, Inhibited Cigarette Smoke-Induced Pulmonary Inflammation and Apoptosis in Rats. Inflammation Fluoxetine is Neuroprotective in Early Brain Injury via its Anti-inflammatory and Anti-apoptotic Effects in a Rat Experimental Subarachnoid Hemorrhage Model Anti-inflammatory effects of fluoxetine in lipopolysaccharide(LPS)-stimulated microglial cells Disruption of NF-kappaB signaling by fluoxetine attenuates MGMT expression in glioma cells Fluoxetine inhibits NF-kappaB signaling in intestinal epithelial cells and ameliorates experimental colitis and colitis-associated colon cancer in mice Effects of escitalopram, R-citalopram, and reboxetine on serum levels of tumor necrosis factor-alpha, interleukin-10, and depression-like behavior in mice after lipopolysaccharide administration Fluoxetine protects against OVA induced bronchial asthma and depression in rats Anti-inflammatory properties of desipramine and fluoxetine Network-based drug repurposing for novel coronavirus 2019-nCoV/SARS-CoV-2 Association between antidepressant use and reduced risk of intubation or death in hospitalized patients with COVID-19: results from an observational study The serotonin reuptake inhibitor Fluoxetine inhibits SARS-CoV-2 Identification of candidate repurposable drugs to combat COVID-19 using a signature-based approach Oxytocin's anti-inflammatory and proimmune functions in COVID-19: a transcriptomic signature-based approach Accelerating Drug Discovery and Repurposing by Combining Transcriptional Signature Connectivity with Docking. bioRxiv Random-set methods identify distinct aspects of the enrichment signal in gene-set analysis. The Annals of Applied Statistics Connecting omics signatures of diseases, drugs, and mechanisms of actions with iLINCS. bioRxiv COVID-19 and treatment with NSAIDs and corticosteroids: should we be limiting their use in the clinical setting? Ecancermedicalscience Cytokine responses in porcine respiratory coronavirus-infected pigs treated with corticosteroids as a model for severe acute respiratory syndrome Dexamethasone in Hospitalized Patients with Covid-19 -Preliminary Report Coordinated binding of NF-kappaB family members in the response of human cells to lipopolysaccharide Trans-signaling is a dominant mechanism for the pathogenic actions of interleukin-6 in the brain The balance of interleukin (IL)-6, IL-6.soluble IL-6 receptor (sIL-6R), and IL-6.sIL-6R.sgp130 complexes allows simultaneous classic and trans-signaling Activation of the JAK/STAT3 and PI3K/AKT pathways are crucial for IL-6 trans-signaling-mediated pro-inflammatory response in human vascular endothelial cells Physiological effects of modulating the interleukin-6 axis Rolipram Protects Mice from Gram-negative Bacterium Escherichia coli-induced Inflammation and Septic Shock