Esomeprazole attenuates inflammatory and fibrotic response in lung cells through the MAPK/Nrf2/HO1 pathway
Journal of Inflammation volume 18, Article number: 17 (2021)
Idiopathic pulmonary fibrosis (IPF) is an orphan disease characterized by progressive loss of lung function resulting in shortness of breath and often death within 3–4 years of diagnosis. Repetitive lung injury in susceptible individuals is believed to promote chronic oxidative stress, inflammation, and uncontrolled collagen deposition. Several preclinical and retrospective clinical studies in IPF have reported beneficial outcomes associated with the use of proton pump inhibitors (PPIs) such as esomeprazole. Accordingly, we sought to investigate molecular mechanism(s) by which PPIs favorably regulate the disease process.
We stimulated oxidative stress, pro-inflammatory and profibrotic phenotypes in primary human lung epithelial cells and fibroblasts upon treatment with bleomycin or transforming growth factor β (TGFβ) and assessed the effect of a prototype PPI, esomeprazole, in regulating these processes.
Our study shows that esomeprazole controls pro-inflammatory and profibrotic molecules through nuclear translocation of the transcription factor nuclear factor-like 2 (Nrf2) and induction of the cytoprotective molecule heme oxygenase 1 (HO1). Genetic deletion of Nrf2 or pharmacological inhibition of HO1 impaired esomeprazole-mediated regulation of proinflammatory and profibrotic molecules. Additional studies indicate that activation of Mitogen Activated Protein Kinase (MAPK) pathway is involved in the process. Our experimental data was corroborated by bioinformatics studies of an NIH chemical library which hosts gene expression profiles of IPF lung fibroblasts treated with over 20,000 compounds including esomeprazole. Intriguingly, we found 45 genes that are upregulated in IPF but downregulated by esomeprazole. Pathway analysis showed that these genes are enriched for profibrotic processes. Unbiased high throughput RNA-seq study supported antifibrotic effect of esomeprazole and revealed several novel targets.
Taken together, PPIs may play antifibrotic role in IPF through direct regulation of the MAPK/Nrf2/HO1 pathway to favorably influence the disease process in IPF.
Idiopathic pulmonary fibrosis (IPF) is a rare but deadly form of lung disease characterized by progressive loss of lung function that culminates in shortness of breath and often death within 3–4 years from the time of diagnosis. The disease has an incidence of 93.7 cases per 100,000  and its prevalence increases with age with a median age at diagnosis of 60 years [1, 2]. It is believed that repetitive subclinical injuries to the lungs of susceptible individuals favor a lung microenvironment that tips the balance towards pro-oxidant, pro-inflammatory, and profibrotic processes. A few years ago, the FDA approved two drugs for the treatment of IPF; pirfenidone and nintedanib [3, 4]. Although this is a significant landmark in our quest to treat IPF, these two drugs only slow the disease progression but are unable to halt the disease or reverse established fibrosis to cure IPF. Accordingly, there is an opportunity to search for and develop more effective therapies either from libraries of new chemical entities (NCEs) or through repurposing of existing drugs approved for other indications.
Among existing drugs, proton pump inhibitors (PPIs) have promising potential to be repurposed for the treatment of IPF [5, 6]. Originally approved to reduce gastric acidity, a number of studies have linked the use of PPIs with improvement in measures of lung function leading to significantly longer transplant-free survival time in patients with well-defined IPF [7,8,9,10,11]. A recent official clinical practice guideline representing leading Thoracic Societies also conditionally recommended the use of PPIs for IPF . However, mechanistic understanding of how PPIs regulate processes involved in lung remodeling is lacking. In this regard, we recently reported that PPIs directly inhibit an enzyme, dimethylarginine dimethylaminohydrolase (DDAH) , that is upregulated in lung tissues explanted from IPF patients [13, 14], and has been shown to promote experimental lung fibrosis .
DDAH is a cytosolic enzyme expressed in almost every mammalian cell to regulate levels of its endogenous substrate asymmetric dimethylarginine (ADMA) . ADMA is a competitive inhibitor of nitric oxide (NO) synthase (NOS) including the inducible isoform (iNOS). Inducible NOS has been shown to be significantly upregulated in animal models of lung injury  and its genetic or pharmacological inhibition has been demonstrated to be protective in models of bleomycin-induced lung injury [10, 16, 17]. It has also been reported that IPF patients have higher levels of circulating peroxynitrite (OONO−); a highly reactive oxidant and nitrating molecule that is generated by iNOS during inflammatory processes [18, 19]. Overall, the DDAH/iNOS pathway is pathologically upregulated in experimental and clinical lung fibrosis, and PPIs inhibit this pathway .
Since esomeprazole is not a selective DDAH/iNOS inhibitor, we hypothesized that the drug may target other key biological molecules to control lung inflammation and fibrosis. Accordingly, we performed cellular and molecular studies to screen known mediators of lung inflammation and fibrosis to determine if they are targeted by the drug, and whether the drug relies on these biological targets to control inflammatory and fibrotic responses in normal and IPF-derived lung cells. To validate our cell biological studies and identify gene networks, we performed bioinformatics analysis of the NIH Library of Integrated Network-Based Cellular Signatures (LINCS) database, which hosts gene expression profiles of IPF lung fibroblasts treated with over 20,000 compounds including esomeprazole. We also performed unbiased high throughput RNA-seq studies to verify the effect of esomeprazole on inflammation- and fibrosis- related genes/pathways, and to identify novel targets.
Materials and methods
Primary human IPF lung fibroblasts were purchased from Lonza (Walkersville, MD; cat # CC-7231) and cultured in DMEM (ThermoFisher; Waltham, MA; cat # 11995065) supplemented with 10% FBS and 1% Penicillin-Streptomycin (10,000 U/ml). Human primary Bronchial/Tracheal Epithelial Cells were purchased from Lonza (cat # CC-2541) and cultured in BEGM™ Bronchial Epithelial Cell Growth Medium BulletKit™ (Lonza; cat # CC-3170). Human lung endothelial cells were also purchased from Lonza (cat # CC-2527) and cultured in EGM2-MV medium (Lonza; cat # CC-3202).
Induction of oxidative stress and pro-inflammatory response by bleomycin
Lung epithelial cells (7 × 105 cells) were seeded in 25 cm2 flasks and incubated in a humidified 5% CO2 incubator at 37 °C. Once the cells reached about 70% confluency, they were induced with bleomycin (25 μg/ml final concentration in 20 μL) or control (equal volume of water), and cultured in 4 mL fully-supplemented media in the absence or presence of different concentrations of esomeprazole (1–100 μM) for 24 h. Subsequently, the cells were harvested for gene and protein expression studies as described below.
Induction of profibrotic response by TGFβ
Lung fibroblasts (7 × 105 cells) were seeded in 25 cm2 flasks and incubated at 37 °C/5%CO2. At about 70% confluency, the conditioned media was replaced with fresh DMEM and the cells were induced with recombinant human transforming growth factor beta (TGFβ; Peprotech; Rocky Hill, NJ; cat # 100–21) for 5 days at a final concentration of 10 ng/ml in the absence or presence of various concentrations of esomeprazole (1–100 μM). The culture medium was exchanged every 48 h with fresh media containing TGFβ with or without esomeprazole. The expression of extracellular matrix (ECM) components including collagen was assessed by quantitative RT-PCR as described below.
Effect of esomeprazole on Nrf2/HO1 antioxidant pathway
Lung epithelial cells or fibroblasts were cultured as described above. When the cells reached about 70% confluency, they were treated with various concentrations of esomeprazole (1–100 μM) for 6 to 24 h. Total RNA and protein were isolated for gene and protein expression studies, respectively. In some samples, nuclear and cytoplasmic proteins were separated for western blot studies to evaluate translocation of the transcription factor nuclear factor-like 2 (Nrf2).
Pharmacological inhibition of HO1, Nrf2, ERK1/2, and MEK1/2
To understand the dependence of esomeprazole on heme oxygenase 1 (HO1) to control inflammatory and fibrotic processes, as well as to elucidate key signaling molecules involved in the process, we assessed HO1 and its upstream signaling that includes Nrf2, extracellular signal-regulated kinase 1/2 (ERK1/2) and MAP-ERK kinase 1/2 (MEK1/2) [20, 21] in IPF lung fibroblasts. The cells were cultured as described above and treated with pharmacological inhibitors of HO1 (Tin protoporphyrin IX; SnPPIX; Cayman Chemical; Ann Arbor, MI; cat # 16375; 1 μM final concentration), Nrf2 (Trigonelline; United States Pharmacopeia; Rockville, MD; cat # 1686411; 10 μM final concentration), or MEK1/2 (U0126; Tocris Bioscience; Bristol, UK; cat # 1144; 10 μM final concentration) for 24 h. Quantitative RT-PCR and western blot were used to assess the effect of inhibiting these targets on the upregulation of HO1 by esomeprazole. In some studies, the effect of esomeprazole on these targets including translocation of Nrf2 and phosphorylation of ERK1/2 was assessed. In other studies, the effect of HO1 inhibition (by SnPPIX) on the regulation of pro-inflammatory (TNFα, IL-1β and IL-6) and profibrotic molecules (collagen types I, III and V) by esomeprazole was assessed. To complement the pharmacological approach, lung fibroblasts were isolated from Nrf2 knockout mice (8–10 weeks old), and the dependence of esomeprazole on Nrf2 to activate HO1 was assessed. To minimize age-dependent impairment of HO1 expression in the Nrf2 knockout cells, early passage (below passage 5) fibroblasts were used.
Gene expression analysis
Total RNA was isolated from lung fibroblasts, epithelial and endothelial cells using Direct-zol™ RNA MiniPrep Kit (Zymo Research; Irvine, CA; cat # R2050). The quality and quantity of RNA was evaluated by NanoQuant Plate™ using Tecan Spark 20 M spectrophotometer. A total of 1 μg RNA was reverse transcribed into cDNA using the SuperScript™ VILO™ Master Mix (ThermoFisher; cat # 11755) in a final reaction volume of 20 μl following the manufacturer’s protocol. The resulting cDNA was used for gene expression studies using TaqMan Gene Expression Assay (ThermoFisher; protocol # 4333458) and “best coverage” primers (ThermoFisher) on a Bio-Rad CFX96 RT-PCR system (Bio-Rad; Hercules, CA). Finally, the gene expression data was analyzed using the CFX Maestro Software (Bio-Rad) and fold changes in gene expression were calculated after normalizing to β-actin (ACTB).
Cytoplasmic and nuclear proteins were separated using NE-PER™ Nuclear and Cytoplasmic Extraction Kit (ThermoFisher; cat # 78833) following the manufacturer’s protocol. For the isolation of total cellular protein, cells were PBS-washed and lysed using cell lysis buffer (10 mM Na2HPO4, 0.1% Triton-X100; pH 7.4) containing 1x Halt™ Protease and Phosphatase Inhibitor Cocktail (ThermoFisher; cat # 78440). Subsequently, the cell lysate was vortexed for 1 min and incubated on ice for 30 min prior to recovering the protein by centrifugation at 16000 x g for 15 min at 4 °C. The protein concentration was determined spectrophotometrically, and equal amounts of protein (30 μg) were separated on Bis-Tris pre-cast SDS-PAGE mini-gel (ThermoFisher; cat # NP0322). The gel was transferred to polyvinylidene fluoride membrane using an iBlot™ 2 dry blotting system (ThermoFisher). Non-fat milk (5%) in Tris-buffered saline (TBS) containing 0.1% Tween-20 (TBST; pH 7.4) was used to block the membrane by incubating for 1 h at room temperature. Subsequently, the immunoblot membranes were incubated overnight with target-specific antibody including anti-HO1 (Enzo Life Sciences, Inc.; Farmingdale, NY; cat # BML-HC3001; 1:250), anti-Nrf2 (Abcam; Cambridge, MA; cat # ab62352; 1:250), anti-Keap1 (Abcam; cat # ab1194403; 1:1000), anti-phosho-ERK1/2 (Cell Signaling Technology; Danvers, MA; cat # 4370; 1:2000), and anti-phospho-MEK1/2 (Santa Cruz Biotechnology; Santa Cruz, CA; cat # sc-81,503; 1:200). For internal controls, anti-β-actin (Sigma; St. Louis, MO; cat # A2066; 1:1500) was used as control for cytoplasmic and total protein extracts, and anti-histone H3 (Abcam; cat # ab1791; 1:1000) was used as control for nuclear extracts. In some cases, total (phosphorylated plus unphosphorylated) protein was used as control to compare phosphorylation status of the same protein. Finally, the membranes were incubated for 1 h at room temperature with HRP-conjugated goat anti-rabbit (GE Healthcare; Chicago, IL; cat # NA 934 V; 1:5000), or goat anti-mouse (GE Healthcare; cat # NA 931 V; 1:5000) secondary antibody. Protein bands were visualized using Amersham ECL Prime Western Blotting Detection Reagent (Amersham Biosciences; Little Chalfont, UK; cat # RPN2236) using a ChemiDoc XRS imager system (Bio-Rad). Data was quantified using NIH’s Image J software.
Human lung epithelial cells were cultured in 6-well plates until 70% confluency and incubated further in the absence or presence of esomeprazole (1–100 μM) for 24 h. The cells were fixed with paraformaldehyde (4%) and permeabilized with 0.1% triton-X100 for 15 min prior to staining for HO1 using mouse anti-HO1 antibody (ThermoFisher; cat # MA1–112; 1:100). Secondary anti-mouse antibody (Santa Cruz; SC-362277; 1:5000) was used to detect signal. The cell membrane was stained with Alexa Fluor 488-conjugated phalloidin (ThermoFisher; cat # A12379; 1:400).
IPF gene signature containing lung transcriptomic dataset (GSE53845)  was queried against the LINCS database (http://www.lincscloud.org/) that contains over 1 million gene-expression profiles of cells treated with different small molecules, including esomeprazole . The comparison was performed using The Connectivity Map resource  that was developed to connect small molecules, genes and diseases using gene-expression signatures. Z-scored differential expressions were calculated for both upregulated and downregulated gene sets, and a list of hits were ranked from highest to lowest on the basis of their match strength. Functional enrichment analysis of the inversely correlated gene sets between IPF lungs and esomeprazole-treated cells from the LINCS database was generated using ToppFun  and visualized using Cytoscape software .
Total RNA was isolated from mouse lung fibroblasts treated with TGFβ (10 ng/mL) in the absence or presence of esomeprazole (100 μM) as described above, and integrity of the RNA was verified using Bioanalyzer 2100 (Agilent Technologies; Santa Clara, CA) at the Baylor College of Medicine, Genomic and RNA Profiling Core. Subsequently, messenger RNA (mRNA) library was prepared at Novogene Corporation (sacramento, CA) and subjected to sequencing using NovaSeq 6000 PE150 sequencing platform. Finally, the data was log-transformed, and clustered using KEGG and GO Enrichment Analyses of differentially expressed genes (DEGs). Statistical significance was set at p value below 0.05 (p < 0.05) and, fold-change threshold at 1.5.
Data are presented as mean value ± standard error mean (SEM) from at least triplicate experiments. For the RNA-seq studies, samples were run in duplicates for comparison. Multiple groups were compared using one-way analysis of variance (ANOVA) followed by Bonferroni posttest. Comparison between two groups was performed using student’s t-test (GraphPad Prism; La Jolla, CA). Densitometric analysis of western blot data was performed using Image J and compared for statistical significance using ANOVA.
Esomeprazole inhibits the DDAH/iNOS pathway
The DDAH/iNOS pathway is known to be involved in aberrant lung remodeling including lung inflammation and fibrosis. Genetic or pharmacological inhibition of DDAH or iNOS has been shown to significantly attenuate fibrosis and improve lung compliance [13, 16, 17]. In this study, we found that the expression of both DDAH and iNOS were significantly downregulated by esomeprazole (Fig. 1). For example, the gene expression of iNOS was increased by about 100-fold upon induction with bleomycin, and treatment with esomeprazole decreased the induction in iNOS expression by about 25-fold (Fig. 1).
Esomeprazole upregulates HO1 expression
The beneficial role of HO1 in lung homeostasis is well recognized [27,28,29,30,31]. As a rate-limiting enzyme in the breakdown of heme into three bioactive metabolites (i.e., carbon monoxide, bilirubin, and ferritin), HO1 is regarded as a critical mediator of lung physiology. The expression of HO1 was upregulated by esomeprazole in a dose-dependent manner (Fig. 2). As expected, the Western blot data showing induction of HO1 was reproduced by immunofluorescence (Fig. S1).
Esomeprazole translocates Nrf2 into the nucleus
Physiologically, Nrf2 is kept inactive by forming a complex with Kelch ECH associating protein 1 (Keap1) in the cytoplasm of mammalian cells . As a cytoprotective response to oxidative and electrophilic stress, Nrf2 is dissociated from Keap1 and translocates to the nucleus to activate antioxidant enzymes including HO1 and NADPH quinone oxidoreductase 1 (NQO1). In this study, we found that treatment of IPF lung fibroblasts with esomeprazole caused significant accumulation of Nrf2 protein in the nuclei (Fig. 3). To support this finding, we tested Nrf2 translocation in primary human lung endothelial cells and found similar results (Fig. S2). In addition to HO1, we also found that the expression of NQO1 was significantly and dose dependently upregulated by esomeprazole (Fig. S3). However, it is unlikely that the translocation of Nrf2 by esomeprazole is due to direct inhibition of Keap1 expression since we did not find changes in Keap1 protein expression (Fig. S4).
Esomeprazole activates MAPK via phosphorylation
The mitogen-activated protein kinase (MAPK) signaling cascade includes MEK1/2 and ERK1/2; kinase enzymes which activate several other proteins via phosphorylation. For example, the ERK signaling pathway, which represents ubiquitously expressed protein kinases such as ERK1 (44 kDa) and ERK2 (42 kDa), is involved in intracellular signaling functions to influence diverse biological processes including cell adhesion, proliferation, growth, and differentiation. Phosphorylation of ERKs by MEK1/2 is necessary for activation of their kinase enzymatic activity. Activated ERK1/2 regulates phosphorylation of hundreds of regulatory molecules and transcription factors including Nrf2 [33, 34]. In our study, we found that treatment of IPF lung fibroblasts with esomeprazole phosphorylated both ERK1/2 (Fig. 4a) and MEK1/2 (Fig. 4b). This observation was also reproduced in lung epithelial (Fig. S5) and endothelial cells (Fig. S6). As expected, there was no change in the expression of total ERK protein (Fig. S7).
Esomeprazole downregulates bleomycin-induced pro-inflammatory molecules through the HO1 pathway
The inflammatory molecule tumor necrosis factor α (TNFα) has been shown to activate latent TGFβ  while overexpression of interleukin 1β (IL-1β) sustains TGFβ expression and promotes lung fibrosis . Here we linked the dependence of esomeprazole on HO1 to suppress TNFα and interleukins such as IL-1β and IL-6 (Fig. 5). As shown, the HO1 inhibitor SnPPIX attenuated the anti-inflammatory effects of esomeprazole.
Esomeprazole downregulates TGFβ-induced collagen expression and the effect is blunted by pharmacological inhibition of the MAPK pathway
As major components of the ECM, several collagen types (e.g., type I, III, V) are pathologically involved in abnormal deposition of collagen and progression of lung fibrosis. Here, we report that esomeprazole significantly inhibits the expression of major collagen types and that this inhibition is mediated by the MAPK/ERK signaling (Fig. 6). As shown in the figure, pharmacological inhibition of MEK1/2 abrogated the antifibrotic effect of esomeprazole.
Inhibition of Nrf2 impairs esomeprazole-mediated HO1 activation
Activation of HO1 appears to be an important mechanism by which esomeprazole exerts antioxidant, anti-inflammatory and antifibrotic functions. In this study, we investigated whether activation of HO1 by esomeprazole is dependent on the transcription factor Nrf2. Pharmacological inhibition of Nrf2 with trigonelline (Fig. 7a) or genetic deletion of Nrf2 (Fig. 7b) significantly impaired the activation of HO1 by esomeprazole.
Inhibition of MEK1/2 impairs esomeprazole-mediated HO1 activation
The activation of HO1 by esomeprazole appears to be upstream of Nrf2/Keap1 towards the MAPK pathway. Accordingly, in addition to the phosphorylation of ERK1/2 by esomeprazole, pharmacological inhibition of MEK1/2 significantly impaired esomeprazole-mediated regulation of HO1 (Fig. 8).
Esomeprazole favorably regulates several myofibroblast- and lung development- related genes
Given that several retrospective clinical studies reported beneficial outcomes associated with the use of PPIs in IPF, we sought to gain mechanistic understanding of processes by which PPIs regulate lung remodeling. Accordingly, we queried the NIH’s Library of Integrated Network-Based Cellular Signatures (LINCS)  which hosts gene expression profiles of over 20,000 compounds including some of the most common drugs such as esomeprazole. Intriguingly, we found 45 genes that are significantly upregulated in IPF  but downregulated by esomeprazole (Fig. 9 and Table S1). These genes are enriched for profibrotic processes including ECM proteins such as collagen and matrix metalloproteinases (MMPs). We also found 34 lung development related genes that are downregulated in IPF but significantly upregulated by esomeprazole (Fig. 9 and Table S1). Interestingly, functional enrichment analysis using ToppFun application of the ToppGene Suite  showed that esomeprazole suppressed ECM, collagen metabolism and myofibroblast activation while favorably regulating genes involved in lung development, angiogenesis and wound healing (Fig. 10 and Table S1). Strikingly, RNA-seq data independently corroborated the effect of esomeprazole on ECM-related pathways including downregulation of several collagen types and other extracellular matrix components (Fig. 11, Fig. S8, and Table S2). By contrast, several lung development and wound healing related genes, including heme oxygenase 1, were upregulated by esomeprazole (Table S3).
Esomeprazole is a pleiotropic molecule: activation of HO1 in lung cells
Mounting evidence indicates that PPIs possess biological activities that extend beyond suppression of gastric acidity into antioxidant, anti-inflammatory and antifibrotic properties [37,38,39,40,41]. The antioxidant property of PPIs is reported to be due to direct scavenging of reactive oxygen species (ROS), restoration of detoxifying enzymes such as glutathione (GSH) and HO1 [37, 39, 41,42,43]. As a rate-limiting enzyme in the breakdown of heme into carbon monoxide (CO), bilirubin and ferritin, HO1 plays important role in the pathobiology of lung diseases including pulmonary fibrosis. For example, CO possesses multiple biological functions including antioxidant, anti-inflammatory and bronchodilator [31, 44], and ferritin acts as an antioxidant by sequestering free iron and suppressing iron-dependent redox reaction . Bilirubin is a cytoprotective molecule that has been shown to attenuate experimental lung fibrosis . HO1 itself, when overexpressed, plays protective role in animal models of lung fibrosis [31, 46]. By contrast, levels of HO1 have been found to be reduced in alveolar macrophages isolated from IPF patients and in areas of active fibrosis [47, 48]. Intriguingly, our study shows that a classic PPI, esomeprazole, significantly and dose dependently upregulates the expression of HO1 (Fig. 2). The upregulation of HO1 appears to involve nuclear translocation of Nrf2 (Figs. 3 and 7) and phosphorylation of MAPK family members (Fig. 4). Among several possibilities, dissociation of Nrf2 from Keap1 may be triggered by kinases that are upstream of the transcription factor including activated (i.e. phosphorylated) ERK1/2 and MEK1/2 [33, 34] (Fig. 12). Follow up studies should interrogate signaling cascades that are upstream of ERK and MEK including Raf, Ras and tyrosine receptor kinases. Intriguingly, pharmacological inhibition or genetic deletion of Nrf2 significantly impaired HO1 activation (Fig. 7) suggesting that activation of HO1 by esomeprazole is Nrf2-dependent. Our data also shows that the anti-inflammatory action of esomeprazole is mediated by HO1 (Fig. 5), and selective inhibition of HO1 impairs esomeprazole’s effect on key inflammatory molecules. In addition, the effect of esomeprazole on the DDAH/iNOS pathway may contribute to its overall anti-inflammatory effect (Fig. 1).
MAPK signaling pathway is involved in antifibrotic action of esomeprazole
In addition to the dependence of esomeprazole on the HO1 pathway to regulate processes involved in lung inflammation, the compound appears to depend on the MAPK signaling pathway to control TGFβ-induced expression of collagen (Fig. 6). In addition to impairing esomeprazole-mediated HO1 activation (Fig. 8), selective inhibition of MEK with U0126 appears to reduce TGFβ-induced expression of collagen compared to TGFβ only control (Fig. 6). This is consistent with an earlier report that demonstrated antifibrotic effect of MEK inhibition in experimental lung fibrosis . Taken together, the activation of HO1 through phosphorylation of ERK/MEK and nuclear translocation of Nrf2, as well as inhibition of the pro-inflammatory and profibrotic DDAH/iNOS pathway contributes to the overall anti-inflammatory and antifibrotic properties of esomeprazole.
Esomeprazole favorably regulates a network of genes involved in lung fibrosis
Consistent with our cell culture data, our bioinformatics analysis and RNA-seq studies comparing differentially expressed genes with the transcriptome of the PPI esomeprazole indicate that the antifibrotic action of esomeprazole may broadly involve suppression of ECM components including several types of collagen (e.g. collagen I, III, V) (Fig. 9 and Table S1). In fact, the RNA-seq study revealed that esomeprazole significantly downregulated over 800 extracellular matrix-related genes including collagens (Fig. 11, Fig. S8 and Table S2). The various collagen types that are upregulated in IPF fibroblasts (e.g. Figure 9) are expected to contribute to the abnormal accumulation of ECM proteins. Notably, other genes that are upregulated in IPF but found to be downregulated by esomeprazole include uncoupling protein 2 (UCP2), growth factors (e.g., FGF18), matrix metalloproteinases (e.g., MMP11, MMP13), Wilms’ Tumor 1 (WT1) and TWIST1 (Fig. 10 and Table S1). These proteins are known to play profibrotic role through increased ECM deposition, reduced collagen metabolism and/or transdifferentiation of fibroblasts into collagen-synthesizing myofibroblasts. For example, UCP2 is significantly induced by TGFβ and plays profibrotic role . The role of WT1 and TWIST1 as profibrotic molecules is also well established [51,52,53]. By contrast, esomeprazole upregulated several lung development- and wound healing- related genes that are found to be downregulated in the control IPF lung fibroblasts (Fig. 10). For example, upregulation of Nbeal2 is expected to suppress myofibroblast infiltration and accelerate normal wound healing .
In conclusion, PPIs in general and esomeprazole in particular may provide beneficial effect in IPF through upregulation of antioxidant, anti-inflammatory and antifibrotic molecules, as well as suppression of ECM proteins that are involved in collagen metabolism and myofibroblast activation. Given the premises of our data and the conditional recommendation of PPIs for the treatment of IPF , randomized controlled clinical trials evaluating the efficacy of PPIs for the treatment of IPF are warranted. However, such studies need to be designed with the intent to treat fibrosis and not gastric reflux per se. Accordingly, the dose of esomeprazole and other PPIs need to be adjusted to achieve plasma concentration of 50–100 μM in order to reliably regulate processes involved in lung remodeling. The standard antacid doses of PPIs only achieve plasma drug concentrations of 10–20 μM [55, 56]. However, higher doses of PPIs that can achieve antifibrotic concentrations can safely be administered to patients . Limitations of the study include the use of bleomycin and TGFβ to induce pro-inflammatory and profibrotic processes to understand an “idiopathic” disease. Nevertheless, there is no perfect model to study IPF at preclinical level and these stimuli are often used to study the disease process in vitro and in animal models.
Availability of data and materials
The datasets and materials used and analyzed during the current study are available from the corresponding author on reasonable request.
Raghu G, Chen SY, Yeh WS, Maroni B, Li Q, Lee YC, et al. Idiopathic pulmonary fibrosis in US Medicare beneficiaries aged 65 years and older: incidence, prevalence, and survival, 2001-11. Lancet Respir Med. 2014;2(7):566–72. https://doi.org/10.1016/S2213-2600(14)70101-8.
Hodgson U, Laitinen T, Tukiainen P. Nationwide prevalence of sporadic and familial idiopathic pulmonary fibrosis: evidence of founder effect among multiplex families in Finland. Thorax. 2002;57(4):338–42. https://doi.org/10.1136/thorax.57.4.338.
King TE Jr, Bradford WZ, Castro-Bernardini S, Fagan EA, Glaspole I, Glassberg MK, et al. A phase 3 trial of pirfenidone in patients with idiopathic pulmonary fibrosis. N Engl J Med. 2014;370(22):2083–92. https://doi.org/10.1056/NEJMoa1402582.
Richeldi L, du Bois RM, Raghu G, Azuma A, Brown KK, Costabel U, et al. Efficacy and safety of nintedanib in idiopathic pulmonary fibrosis. N Engl J Med. 2014;370(22):2071–82. https://doi.org/10.1056/NEJMoa1402584.
Ghebre YT. Proton pump inhibitors in IPF: a call for clinical trials. Front Pharmacol. 2018;9:499. https://doi.org/10.3389/fphar.2018.00499.
Ghebre YT, Raghu G. Idiopathic pulmonary fibrosis: novel concepts of proton pump inhibitors as Antifibrotic drugs. Am J Respir Crit Care Med. 2016;193(12):1345–52. https://doi.org/10.1164/rccm.201512-2316PP.
Raghu G, Freudenberger TD, Yang S, et al. High prevalence of abnormal acid gastro-oesophageal reflux in idiopathic pulmonary fibrosis. Eur Respir J. 2006;27(1):136–42. https://doi.org/10.1183/09031936.06.00037005.
Lee JS, Ryu JH, Elicker BM, Lydell CP, Jones KD, Wolters PJ, et al. Gastroesophageal reflux therapy is associated with longer survival in patients with idiopathic pulmonary fibrosis. Am J Respir Crit Care Med. 2011;184(12):1390–4. https://doi.org/10.1164/rccm.201101-0138OC.
Lee J, Collard HR, Anstrom KJ, Martinez FJ, Noth I, Roberts RS, et al. Anti-acid treatment and disease progression in idiopathic pulmonary fibrosis: an analysis of data from three randomised controlled trials. Lancet Respir Med. 2013;1(5):369–76. https://doi.org/10.1016/S2213-2600(13)70105-X.
Ghebremariam YT, Cooke JP, Gerhart W, Griego C, Brower JB, Doyle-Eisele M, et al. Pleiotropic effect of the proton pump inhibitor esomeprazole leading to suppression of lung inflammation and fibrosis. J Transl Med. 2015;13(1):249. https://doi.org/10.1186/s12967-015-0614-x.
Lee CM, Lee DH, Ahn BK, Hwang JJ, Yoon H, Shin CM, et al. Protective effect of proton pump inhibitor for survival in patients with Gastroesophageal reflux disease and idiopathic pulmonary fibrosis. J Neurogastroenterol Motil. 2016;22(3):444–51. https://doi.org/10.5056/jnm15192.
Raghu G, Rochwerg B, Zhang Y, Garcia CA, Azuma A, Behr J, et al. An official ATS/ERS/JRS/ALAT clinical practice guideline: treatment of idiopathic pulmonary fibrosis. An update of the 2011 clinical practice guideline. Am J Respir Crit Care Med. 2015;192(2):e3–e19. https://doi.org/10.1164/rccm.201506-1063ST.
Pullamsetti SS, Savai R, Dumitrascu R, et al. The role of dimethylarginine dimethylaminohydrolase in idiopathic pulmonary fibrosis. Sci Transl Med. 2011;3:87ra53.
Janssen W, Pullamsetti SS, Cooke J, Weissmann N, Guenther A, Schermuly RT. The role of dimethylarginine dimethylaminohydrolase (DDAH) in pulmonary fibrosis. J Pathol. 2013;229(2):242–9. https://doi.org/10.1002/path.4127.
Palm F, Onozato ML, Luo Z, Wilcox CS. Dimethylarginine dimethylaminohydrolase (DDAH): expression, regulation, and function in the cardiovascular and renal systems. Am J Physiol Heart Circ Physiol. 2007;293(6):H3227–45. https://doi.org/10.1152/ajpheart.00998.2007.
Seimetz M, Parajuli N, Pichl A, Veit F, Kwapiszewska G, Weisel FC, et al. Inducible NOS inhibition reverses tobacco-smoke-induced emphysema and pulmonary hypertension in mice. Cell. 2011;147(2):293–305. https://doi.org/10.1016/j.cell.2011.08.035.
Genovese T, Cuzzocrea S, Di Paola R, et al. Inhibition or knock out of inducible nitric oxide synthase result in resistance to bleomycin-induced lung injury. Respir Res. 2005;6(1):58. https://doi.org/10.1186/1465-9921-6-58.
Saleh D, Barnes PJ, Giaid A. Increased production of the potent oxidant peroxynitrite in the lungs of patients with idiopathic pulmonary fibrosis. Am J Respir Crit Care Med. 1997;155(5):1763–9. https://doi.org/10.1164/ajrccm.155.5.9154889.
Montaldo C, Cannas E, Ledda M, et al. Bronchoalveolar glutathione and nitrite/nitrate in idiopathic pulmonary fibrosis and sarcoidosis. Sarcoidosis Vasc Diffuse Lung Dis. 2002;19:54–8.
Zipper LM, Mulcahy RT. Erk activation is required for Nrf2 nuclear localization during pyrrolidine dithiocarbamate induction of glutamate cysteine ligase modulatory gene expression in HepG2 cells. Toxicol Sci. 2003;73(1):124–34. https://doi.org/10.1093/toxsci/kfg083.
Naidu S, Vijayan V, Santoso S, Kietzmann T, Immenschuh S. Inhibition and genetic deficiency of p38 MAPK up-regulates heme oxygenase-1 gene expression via Nrf2. J Immunol. 2009;182(11):7048–57. https://doi.org/10.4049/jimmunol.0900006.
DePianto DJ, Chandriani S, Abbas AR, et al. Heterogeneous gene expression signatures correspond to distinct lung pathologies and biomarkers of disease severity in idiopathic pulmonary fibrosis. Thorax. 2015;70(1):48–56. https://doi.org/10.1136/thoraxjnl-2013-204596.
Subramanian A, Narayan R, Corsello SM, Peck DD, Natoli TE, Lu X, et al. A next generation connectivity map: L1000 platform and the first 1,000,000 profiles. Cell. 2017;171(6):1437–52 e17. https://doi.org/10.1016/j.cell.2017.10.049.
Lamb J, Crawford ED, Peck D, Modell JW, Blat IC, Wrobel MJ, et al. The connectivity map: using gene-expression signatures to connect small molecules, genes, and disease. Science. 2006;313(5795):1929–35. https://doi.org/10.1126/science.1132939.
Chen J, Bardes EE, Aronow BJ, Jegga AG. ToppGene suite for gene list enrichment analysis and candidate gene prioritization. Nucleic Acids Res. 2009;37(Web Server):W305–11. https://doi.org/10.1093/nar/gkp427.
Shannon P, Markiel A, Ozier O, Baliga NS, Wang JT, Ramage D, et al. Cytoscape: a software environment for integrated models of biomolecular interaction networks. Genome Res. 2003;13(11):2498–504. https://doi.org/10.1101/gr.1239303.
Morse D, Choi AM. Heme oxygenase-1: the “emerging molecule” has arrived. Am J Respir Cell Mol Biol. 2002;27(1):8–16. https://doi.org/10.1165/ajrcmb.27.1.4862.
Morse D, Choi AM. Heme oxygenase-1: from bench to bedside. Am J Respir Crit Care Med. 2005;172(6):660–70. https://doi.org/10.1164/rccm.200404-465SO.
Slebos DJ, Ryter SW, Choi AM. Heme oxygenase-1 and carbon monoxide in pulmonary medicine. Respir Res. 2003;4(1):7. https://doi.org/10.1186/1465-9921-4-7.
Morse D, Lin L, Choi AM, et al. Heme oxygenase-1, a critical arbitrator of cell death pathways in lung injury and disease. Free Radic Biol Med. 2009;47(1):1–12. https://doi.org/10.1016/j.freeradbiomed.2009.04.007.
Constantin M, Choi AJ, Cloonan SM, et al. Therapeutic potential of heme oxygenase-1/carbon monoxide in lung disease. Int J Hypertens. 2012;2012:859235.
Kansanen E, Kuosmanen SM, Leinonen H, Levonen AL. The Keap1-Nrf2 pathway: mechanisms of activation and dysregulation in cancer. Redox Biol. 2013;1(1):45–9. https://doi.org/10.1016/j.redox.2012.10.001.
Busca R, Pouyssegur J, Lenormand P. ERK1 and ERK2 map kinases: specific roles or functional redundancy? Front Cell Dev Biol. 2016;4:53.
Cullinan SB, Zhang D, Hannink M, Arvisais E, Kaufman RJ, Diehl JA. Nrf2 is a direct PERK substrate and effector of PERK-dependent cell survival. Mol Cell Biol. 2003;23(20):7198–209. https://doi.org/10.1128/MCB.23.20.7198-7209.2003.
Sullivan DE, Ferris M, Pociask D, Brody AR. The latent form of TGFbeta (1) is induced by TNFalpha through an ERK specific pathway and is activated by asbestos-derived reactive oxygen species in vitro and in vivo. J Immunotoxicol. 2008;5(2):145–9. https://doi.org/10.1080/15476910802085822.
Kolb M, Margetts PJ, Anthony DC, Pitossi F, Gauldie J. Transient expression of IL-1beta induces acute lung injury and chronic repair leading to pulmonary fibrosis. J Clin Invest. 2001;107(12):1529–36. https://doi.org/10.1172/JCI12568.
Namazi MR, Jowkar F. A succinct review of the general and immunological pharmacologic effects of proton pump inhibitors. J Clin Pharm Ther. 2008;33(3):215–7. https://doi.org/10.1111/j.1365-2710.2008.00907.x.
Kedika RR, Souza RF, Spechler SJ. Potential anti-inflammatory effects of proton pump inhibitors: a review and discussion of the clinical implications. Dig Dis Sci. 2009;54(11):2312–7. https://doi.org/10.1007/s10620-009-0951-9.
Becker JC, Grosser N, Waltke C, Schulz S, Erdmann K, Domschke W, et al. Beyond gastric acid reduction: proton pump inhibitors induce heme oxygenase-1 in gastric and endothelial cells. Biochem Biophys Res Commun. 2006;345(3):1014–21. https://doi.org/10.1016/j.bbrc.2006.04.170.
Sasaki T, Nakayama K, Yasuda H, Yamaya M. A new strategy with proton pump inhibitors for the prevention of acute exacerbations in COPD. Ther Adv Respir Dis. 2011;5(2):91–103. https://doi.org/10.1177/1753465810392264.
Biswas K, Bandyopadhyay U, Chattopadhyay I, Varadaraj A, Ali E, Banerjee RK. A novel antioxidant and antiapoptotic role of omeprazole to block gastric ulcer through scavenging of hydroxyl radical. J Biol Chem. 2003;278(13):10993–1001. https://doi.org/10.1074/jbc.M210328200.
Lapenna D, de Gioia S, Ciofani G, et al. Antioxidant properties of omeprazole. FEBS Lett. 1996;382(1-2):189–92. https://doi.org/10.1016/0014-5793(96)00155-X.
Simon WA, Sturm E, Hartmann HJ, Weser U. Hydroxyl radical scavenging reactivity of proton pump inhibitors. Biochem Pharmacol. 2006;71(9):1337–41. https://doi.org/10.1016/j.bcp.2006.01.009.
Otterbein LE, Bach FH, Alam J, Soares M, Tao Lu H, Wysk M, et al. Carbon monoxide has anti-inflammatory effects involving the mitogen-activated protein kinase pathway. Nat Med. 2000;6(4):422–8. https://doi.org/10.1038/74680.
Wang HD, Yamaya M, Okinaga S, et al. Bilirubin ameliorates bleomycin-induced pulmonary fibrosis in rats. Am J Respir Crit Care Med. 2002;165(3):406–11. https://doi.org/10.1164/ajrccm.165.3.2003149.
Tsuburai T, Suzuki M, Nagashima Y, Suzuki S, Inoue S, Hashiba T, et al. Adenovirus-mediated transfer and overexpression of heme oxygenase 1 cDNA in lung prevents bleomycin-induced pulmonary fibrosis via a Fas-Fas ligand-independent pathway. Hum Gene Ther. 2002;13(16):1945–60. https://doi.org/10.1089/10430340260355356.
Lakari E, Pylkas P, Pietarinen-Runtti P, et al. Expression and regulation of hemeoxygenase 1 in healthy human lung and interstitial lung disorders. Hum Pathol. 2001;32(11):1257–63. https://doi.org/10.1053/hupa.2001.28937.
Ye Q, Dalavanga Y, Poulakis N, Sixt SU, Guzman J, Costabel U. Decreased expression of haem oxygenase-1 by alveolar macrophages in idiopathic pulmonary fibrosis. Eur Respir J. 2008;31(5):1030–6. https://doi.org/10.1183/09031936.00125407.
Madala SK, Schmidt S, Davidson C, Ikegami M, Wert S, Hardie WD. MEK-ERK pathway modulation ameliorates pulmonary fibrosis associated with epidermal growth factor receptor activation. Am J Respir Cell Mol Biol. 2012;46(3):380–8. https://doi.org/10.1165/rcmb.2011-0237OC.
Jiang L, Qiu W, Zhou Y, Wen P, Fang L, Cao H, et al. A microRNA-30e/mitochondrial uncoupling protein 2 axis mediates TGF-beta1-induced tubular epithelial cell extracellular matrix production and kidney fibrosis. Kidney Int. 2013;84(2):285–96. https://doi.org/10.1038/ki.2013.80.
Sontake V, Shanmukhappa SK, DiPasquale BA, et al. Fibrocytes regulate Wilms tumor 1-positive cell accumulation in severe fibrotic lung disease. J Immunol. 2015;195(8):3978–91. https://doi.org/10.4049/jimmunol.1500963.
Sontake V, Kasam RK, Sinner D, Korfhagen TR, Reddy GB, White ES, et al. Wilms' tumor 1 drives fibroproliferation and myofibroblast transformation in severe fibrotic lung disease. JCI Insight. 2018;3(16). https://doi.org/10.1172/jci.insight.121252.
Mammoto T, Jiang A, Jiang E, Mammoto A. Role of Twist1 phosphorylation in angiogenesis and pulmonary fibrosis. Am J Respir Cell Mol Biol. 2016;55(5):633–44. https://doi.org/10.1165/rcmb.2016-0012OC.
Deppermann C, Cherpokova D, Nurden P, Schulz JN, Thielmann I, Kraft P, et al. Gray platelet syndrome and defective thrombo-inflammation in Nbeal2-deficient mice. J Clin Invest. 2013;123(8):3331–42. https://doi.org/10.1172/JCI69210.
Shin JM, Sachs G. Pharmacology of proton pump inhibitors. Curr Gastroenterol Rep. 2008;10(6):528–34. https://doi.org/10.1007/s11894-008-0098-4.
Ghebremariam YT, LePendu P, Lee JC, Erlanson DA, Slaviero A, Shah NH, et al. Unexpected effect of proton pump inhibitors: elevation of the cardiovascular risk factor asymmetric dimethylarginine. Circulation. 2013;128(8):845–53. https://doi.org/10.1161/CIRCULATIONAHA.113.003602.
Brana I, Ocana A, Chen EX, Razak ARA, Haines C, Lee C, et al. A phase I trial of pantoprazole in combination with doxorubicin in patients with advanced solid tumors: evaluation of pharmacokinetics of both drugs and tissue penetration of doxorubicin. Investig New Drugs. 2014;32(6):1269–77. https://doi.org/10.1007/s10637-014-0159-5.
We would like to thank the Departments of Radiation Oncology and Medicine at Baylor College of Medicine for overall support. We would also like to thank Dr. Samana Shrestha (Baylor College of Medicine) for her technical support, and Prof. Dean Sheppard of the University of California, San Francisco for his editorial revision of the manuscript. We would also like to thank Novogene Corporation Inc. (Sacramento, CA) for the RNA-seq study and data analysis.
The study was supported (to YTG) in part by grants from the NHLBI (grant numbers K01HL118683; R01HL137703), NIAMS (grant number R56AR077445), American Heart Association (grant number 17GRNT33460159), the Cancer Prevention and Research Institute of Texas (grant number RP190497), and by intramural funding from Baylor College of Medicine (project ID 2690000104). AE is a recipient of T32 training grant from the NHLBI through Baylor College of Medicine (grant number T32HL007747–22). This content is solely the responsibility of the authors and does not necessarily represent the official views of the sponsors.
Ethics approval and consent to participate
Consent for publication
All authors consent to the publication of this manuscript.
YTG is an inventor on patents, owned by Stanford University and Baylor College of Medicine, that protect the use of agents, including proton pump inhibitors (PPIs), for therapeutic use of new indications. MDB is an inventor on the patent application owned by Baylor College of Medicine.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Comparison and overlap of differentially expressed genes in control and esomeprazole treated IPF lung fibroblasts from the Library of Integrated Network-Based Cellular Signatures (LINCS) database. The data shows that there are 2879 genes that are differentially expressed in the IPF cells and of which 141 genes overlap with differentially expressed genes following esomeprazole treatment. In the first tab of the spreadsheet, the 45 esomeprazole downregulated genes are indicated in blue and the 34 esomeprazole upregulated genes are shown in orange.
RNA-sequencing (RNA-seq) data comparing differential gene expression in mouse lung fibroblasts stimulated with the profibrotic cytokine TGFβ (control; 10 ng/mL) or treated with TGFβ (10 ng/mL) and esomeprazole (100 μM). The fold-change, gene symbol, and description of the significantly downregulated genes (p < 0.05) is shown. Experiment was run in duplicate for each condition.
RNA-sequencing (RNA-seq) data comparing differential gene expression in mouse lung fibroblasts stimulated with the profibrotic cytokine TGFβ (control; 10 ng/mL) or treated with TGFβ (10 ng/mL) and esomeprazole (100 μM). The fold-change, gene symbol, and description of the significantly upregulated genes (p < 0.05) is shown. Experiment was run in duplicate for each condition.
Immunofluorescence data showing increased expression of heme oxygenase 1 (HO1) protein by esomeprazole in primary human lung epithelial cells. The cells were treated with various concentrations of esomeprazole for 24 hours (1-100 μM) prior to staining with mouse anti-HO1 antibody (shown in red). The cell membrane was stained with Alexa Fluor 488-conjugated phalloidin and is shown in green. Figure S2. Western blot data showing nuclear translocation of nuclear factor-like 2 (Nrf2) in human lung endothelial cells. The cells were treated for 24 hours with esomeprazole (1-100 μM) or vehicle control (water) prior to isolation of nuclear protein. Data is representative of five independent experiments. Histone H3 is used as a loading control. Densitometric quantification of the protein bands relative to the housekeeping control protein histone H3 is shown in the lower panel. Figure S3. Quantitative RT-PCR (qRT-PCR) data showing dose-dependent upregulation of NADPH quinone oxidoreductase 1 (NQO1) by esomeprazole in human IPF lung fibroblasts. The cells were treated with various concentrations of esomeprazole for 6 hours prior to isolation of RNA for qRT-PCR. Data is Mean ± SEM from triplicate experiments. *p<0.05 compared to control. Figure S4. Western blot data showing no change in the protein expression of Kelch ECH associating protein 1 (Keap1) by esomeprazole in human lung epithelial cells. The cells were treated for 24 hours with various concentrations of esomeprazole (1-100 μM) or vehicle control prior to isolation of total protein. Data is representative of three independent experiments. Beta actin (ACTB) is used as a loading control. Densitometric quantification of the protein bands relative to ACTB is shown in the lower panel. Figure S5. Western blot data showing phosphorylation of ERK1 and ERK2 in human lung epithelial cells treated with vehicle control or various concentrations of esomeprazole (1-100 μM) for 24 hours. Data is representative of four independent experiments. Beta actin (ACTB) is used as a loading control. Densitometric quantification of the protein bands relative to ACTB is shown in the lower panel. Figure S6. Western blot data showing phosphorylation of ERK1 and ERK2 (pERK1/2) in human lung endothelial cells treated with vehicle control or various concentrations of esomeprazole (1-100 μM) for 24 hours. Data is representative of at least three independent experiments. Beta actin (ACTB) is used as a loading control. Densitometric quantification of the protein bands relative to ACTB is shown in the lower panel. Figure S7. Western blot data showing no change in the expression of total extracellular signal-regulated kinase 1/2 (ERK1/2) upon treatment of human IPF lung fibroblasts with vehicle control or esomeprazole (100 μM) for up to 2 hours. Data is representative of three independent experiments. Beta actin (ACTB) is used as a loading control. Densitometric quantification of the protein bands relative to ACTB is shown in the lower panel. Figure S8. Volcano plot of RNA-seq data from mouse lung fibroblasts stimulated with the profibrotic cytokine TGFβ (control; 10 ng/mL) or treated with TGFβ (10 ng/mL) and esomeprazole (100 μM). The plot shows the total number of significantly upregulated (1876; red) and downregulated (2035; green) genes by esomeprazole. The total number of unchanged genes (17,114) is shown in blue.
About this article
Cite this article
Ebrahimpour, A., Wang, M., Li, L. et al. Esomeprazole attenuates inflammatory and fibrotic response in lung cells through the MAPK/Nrf2/HO1 pathway. J Inflamm 18, 17 (2021). https://doi.org/10.1186/s12950-021-00284-6
- Proton pump inhibitors