Next Article in Journal
Diet, Obesity, and Depression: A Systematic Review
Next Article in Special Issue
Systems Approaches to Treatment Response to Imatinib in Severe Asthma: A Pilot Study
Previous Article in Journal
Mapping Brain Microstructure and Network Alterations in Depressive Patients with Suicide Attempts Using Generalized Q-Sampling MRI
Previous Article in Special Issue
The Role of SNP Interactions when Determining Independence of Novel Signals in Genetic Association Studies—An Application to ARG1 and Bronchodilator Response
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
JPM
Volume 11
Issue 3
10.3390/jpm11030175
jpm-logo
Submit to this Journal Review for this Journal Propose a Special Issue
► Article Menu

Article Menu

share Share announcement Help format_quote Cite question_answer Discuss in SciProfiles thumb_up
...
Endorse textsms
...
Comment
first_page
settings
Order Article Reprints
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Drug Repurposing to Treat Glucocorticoid Resistance in Asthma

by
Alberta L. Wang
1,2,*,
Ronald Panganiban
3,
Weiliang Qiu
1,
Alvin T. Kho
1,4,
Geoffrey Chupp
5,
Deborah A. Meyers
6,
Eugene R. Bleecker
6,
Scott T. Weiss
1,
Quan Lu
3 and
Kelan G. Tantisira
7
1
Channing Division of Network Medicine, Brigham and Women’s Hospital, Boston, MA 02115, USA
2
Division of Allergy and Clinical Immunology, Brigham and Women’s Hospital, Boston, MA 02115, USA
3
Program in Molecular and Integrative Physiological Sciences, Departments of Environmental Health and Molecular Metabolism, Harvard T.H. Chan School of Public Health, Boston, MA 02115, USA
4
Computational Health Informatics Program, Boston Children’s Hospital, Boston, MA 02115, USA
5
Division of Pulmonary, Critical Care, and Sleep Medicine, Yale University School of Medicine, New Haven, CT 06520, USA
6
Division of Genetics, Genomics and Precision Medicine, Department of Medicine, University of Arizona, Tucson, AZ 85724, USA
7
Division of Respiratory Medicine, Department of Pediatrics, University of California San Diego, San Diego, CA 92123, USA
*
Author to whom correspondence should be addressed.
J. Pers. Med. 2021, 11(3), 175; https://doi.org/10.3390/jpm11030175
Submission received: 1 January 2021 / Revised: 17 February 2021 / Accepted: 25 February 2021 / Published: 3 March 2021
(This article belongs to the Special Issue APAA: Asthma Pharmacogenetics across Ages)
Browse Figures
Versions Notes

Abstract

:
Corticosteroid resistance causes significant morbidity in asthma, and drug repurposing may identify timely and cost-effective adjunctive treatments for corticosteroid resistance. In 95 subjects from the Childhood Asthma Management Program (CAMP) and 19 subjects from the Severe Asthma Research Program (SARP), corticosteroid response was measured by the change in percent predicted forced expiratory volume in one second (FEV1). In each cohort, differential gene expression analysis was performed comparing poor (resistant) responders, defined as those with zero to negative change in FEV1, to good responders, followed by Connectivity Map (CMap) analysis to identify inversely associated (i.e., negatively connected) drugs that reversed the gene expression profile of poor responders to resemble that of good responders. Mean connectivity scores weighted by sample size were calculated. The top five drug compound candidates underwent in vitro validation in NF-κB-based luciferase reporter A549 cells stimulated by IL-1β ± dexamethasone. In CAMP and SARP, 134 and 178 respective genes were differentially expressed in poor responders. CMap analysis identified 46 compounds in common across both cohorts with connectivity scores < −50. γ-linolenic acid, ampicillin, exemestane, brinzolamide, and INCA-6 were selected for functional validation. γ-linolenic acid, brinzolamide, and INCA-6 significantly reduced IL-1β induced luciferase activity and potentiated the anti-inflammatory effect of dexamethasone in A549/NF-κB-luc reporter cells. These results demonstrate how existing drugs, including γ-linolenic acid, brinzolamide, and INCA-6, may be repurposed to improve corticosteroid response in asthmatics.

1. Introduction

Asthma is one of the most common chronic respiratory diseases worldwide affecting more than 339 million children and adults [1]. Inhaled corticosteroids (ICS) are considered one of the primary treatments for the control of asthma. However, up to 30–50% of asthmatics have poor response to corticosteroid treatment [2,3]. Corticosteroid refractory asthma represents about 50% of total costs for treating asthma [4]. While the new next generation biologics may help to address the issue of corticosteroid resistance, they remain very expensive and must be administered parenterally. Therefore, interrogations to identify alternatives for the treatment of steroid-refractory asthma are warranted. Drug repurposing studies in asthma to date have focused on identifying new drugs for the primary treatment of asthma and not for the treatment of resistance against corticosteroids or other classes of medications [5,6,7,8,9]. Gene expression studies have demonstrated an important role for transcriptomics in determining ICS treatment response, and drugs that reverse the gene expression profile of poor ICS responders to be more similar to that of ICS good responders may potentially restore treatment response [10,11,12].
The development of a novel drugs is expensive, time consuming, and associated with low success rates. On average, it takes 13–15 years and US $2–3 billion to get a new drug to market, and only 15% of drugs successfully progress to FDA approval [13,14]. Drug repositioning or repurposing reduces drug development time and costs 6.5 years and US $300 million by finding new indications for the 3000 existing approved drugs [13,14]. In fact, two-thirds of the applications approved by the FDA are new indications for existing drugs [15].
The Connectivity Map (CMap) at the Broad Institute is a genome-scale library of more than 1.5 million transcriptional responses to genetic and pharmacologic perturbations in cultured human cell lines that have been used to identify drug repurposing candidates in many fields and diseases, including psychiatry, cancer, diabetes, Alzheimer’s disease, and asthma [7,16,17,18,19,20,21,22,23]. The CMap dataset contains perturbational gene expression profiles from treatment of 54 cultivated cell lines with 5000 small-molecule compounds [16,17]. These reference gene expression profiles can be compared to gene expression profiles from a biological state of interest (i.e., poor corticosteroid response). Compounds that elicit a similar expression signature to the disease state of interest are termed positively connected, and compounds that elicit an expression signature inverse to the disease state of interest are termed negatively connected [16,17]. To date, CMap has predominantly been used to identify novel repurposed drugs for the primary treatment of disease, rather than focusing on compounds that may be used to overcome drug resistance to improve existing treatments for a disease. In this study, we used CMap to identify negatively connected drug compounds that reversed the gene expression signature of corticosteroid poor responders to resemble that of corticosteroid good responders. The top identified drug compounds were validated in in vitro cell line experiments providing novel proof-of-concept of repurposing existing compounds to alleviate corticosteroid resistance in asthmatics.

2. Materials and Methods

2.1. Subjects

The Childhood Asthma Management Program (CAMP) and Severe Asthma Research Program (SARP) III comprise the two study cohorts. CAMP was a multicenter randomized control trial of inhaled budesonide, nedocromil, or placebo for mild-to-moderate persistent asthma in North America [24]. Immortalized B cells were cultured from 145 subjects randomized to inhaled budesonide, and the subjects were divided into tertiles based on ICS response measured by the change in percent predicted forced expiratory volume in one second (FEV1) from baseline at enrollment to month two on inhaled budesonide (Figure S1a) [25]. Parental informed consent with subject assent was obtained, and the study was approved by each CAMP study center institutional review board (IRB) and by the CAMP Coordinating Center (Mass General Brigham IRB 1999P0015492). The 48 children in the lowest tertile of ICS response were considered poor ICS responders (median change in FEV1 % predicted of −2.7 (±6.7), and the 47 children in the highest tertile were considered good responders (median change in FEV1 % predicted of 15.8 (±16.7), Kruskal–Wallis Test p-value < 0.001). For both poor and good responders, the immortalized B cells for each subject were split into two equal halves. One-half was treated with dexamethasone, and the other half was treated with sham. After 6 h, expression levels of 22,184 gene probes were measured using the Illumina HumanRef8 v2 BeadChip (San Diego, CA, USA). Each pair of dexamethasone-treated and sham-treated arrays was in the same batch. Details of the gene expression data quality control have been published [25]. The data were variance-stabilized then quantile normalized and log2-transformed, and 20,917 gene probes passed quality control.
SARP is a multisite longitudinal observational cohort study of children and adults that investigates the mechanisms of severe asthma [26]. From the Wake Forest University site, data were available on 19 adult subjects who underwent extensive phenotyping at baseline and at two research visits prior to and two weeks after an intramuscular injection of 40mg triamcinolone acetonide [27]. Sputum induction was performed before and after the corticosteroid evoked phenotype, and RNA sequencing was obtained on the sputum cells (Appendix A). The IRB at each center approved the study (Wake Forest University IRB 00021507). The data were normalized using the DESeq2 R package then log2-transformed, and 26,422 transcripts present at both visits passed quality control [27,28]. The 19 subjects were divided into good and poor responders based on the change in percent predicted FEV1 before and after steroid administration, which was the same measurement used in CAMP. Good responders (n = 8) had a change in percent predicted FEV1 greater than zero (median 2.6 (±3.9)), and poor responders (n = 11) had a change in percent predicted FEV1 less than or equal to zero (median −4.0 (±5.6), Kruskal–Wallis Test p-value < 0.001) (Figure S1b).

2.2. Statistical Analysis

All statistical analyses were performed in R (Vienna, Austria). Univariate T-tests, Wilcoxon rank-sum tests, χ2 tests, and Fisher’s exact tests were used to evaluate the demographics of good and poor responders within CAMP and SARP. Differential gene expression analysis was performed separately for CAMP and SARP to detect genes differentially expressed between poor and good responders. The outcome variable was the difference in gene expression pre- and post-corticosteroid administration. Multivariable linear regression controlling for age, sex, and ancestry was performed using the limma and iCheck R packages [29,30]. Due to small sample sizes, no genes were differentially expressed after false discovery rate adjustment for multiple testing, and genes with a nominal p-value < 0.01 were taken forward to CMap analysis. The differentially expressed genes were categorized into two groups, genes over-expressed in poor responders and genes under-expressed in poor responders, for input into CMap. The org.Hs.eg.db R package was used to annotate Entrez gene IDs to HGNC gene symbols and chromosome locations [31].

2.3. Connectivity Map Analysis

CMap analysis was performed using CLUE (CMap and LINCS Unified Environment), a cloud computing platform at the Broad Institute that facilitates access and manipulation of CMap data [32]. CMap was queried separately using the gene expression profiles of corticosteroid poor responders from CAMP and SARP to identify similar and dissimilar reference perturbagen signatures. CMap input requires Entrez gene IDs, and therefore, differentially expressed gene without Entrez IDs cannot be queried. Out of all available cell lines in CMap, the A549 cell line was selected for analysis because it is a lung-derived airway epithelial cell, and asthma is, at least in part, a disease of the airway epithelium. Furthermore, the anti-inflammatory effects of inhaled corticosteroid treatment are largely targeted at the airway epithelium, and for consistency, A549 cells were used for functional validation studies as described below. The CMap query outputs a connectivity score for each drug compound perturbagen, ranging from −100 to +100. A negative connectivity score of −100 indicates that the compound completely reversed the gene expression profile of poor responders to resemble that of good responders by inducing the expression of genes that are down-regulated in poor responders and suppressing the expression of genes that are up-regulated in poor responders. Hence, drug compounds with negative connectivity scores less than −50 were prioritized as these drugs had the highest likelihood of reversing corticosteroid resistance. Next, weighted mean connectivity scores were calculated based on sample size for the drug compounds with connectivity scores less than −50 in both cohorts.

2.4. In Vitro Validation

The top 10 available drug compounds with the highest weighted mean negative connectivity scores were screened, and five were chosen for in vitro validation based on safety data in humans, mechanistic plausibility, and drug availability. Functionally validation was performed using a previously described glucocorticoid-mediated tethered transrepression of NF-κB activity assay [33]. A549/NF-κB-luc reporter cells were treated with the selected compounds at indicated concentrations for 1 h. Standard drug concentrations were used (1 µM or 1µg/mL) and increased up to the limit where no cell death occurred. These A549/NF-κB-luc reporter cells are responsive to interleukin (IL)-1β stimulation, which is reduced by dexamethasone, a glucocorticoid receptor agonist, indicating that these cells exhibit glucocorticoid-mediated tethered transrepression of NF-κB [33]. One hour post-treatment, the cells were stimulated with 5 ng/mL IL-1β ± 2 nM dexamethasone. Luciferase assays were performed after 18 hours of treatment using the Luciferase Assay System (Promega, Madison, WI, USA), and luciferase activity was measured using the Synergy 2 plate reader (BioTek, Winooski, VT, USA). The luciferase activity in untreated (no drug compounds) cells exposed to 5 ng/mL IL-1β was normalized to 1. The luciferase activities in cells treated with the drug compounds were compared to that of untreated cells. Data were expressed as means ± standard errors. Statistical significance was evaluated using Student’s t-test (n = 4, p < 0.05).

3. Results

3.1. Demographics

The characteristics of the CAMP and SARP subjects included in this study are presented in Table 1. CAMP was composed of children (mean age 8.7 (±2.2) years) with mostly moderate persistent asthma. The subjects from SARP were all adults (mean age 46.5 (±8.4) years) with mostly severe persistent asthma. Both cohorts had a female sex and European ancestry predominance. In CAMP, the pretreatment percent predicted FEV1 and FEV1 to forced vital capacity ratio were lower in good responders compared to poor responders, and overall, lung function measurements were preserved in CAMP. The demographics of SARP did not vary across responder status. To mitigate the heterogeneity across populations, differential gene expression and CMap analyses were performed independently in each cohort, and then, the independent sets of drug compounds identified on CMap analysis were consolidated.

3.2. Differential Gene Expression Analysis

Differential gene expression analysis in CAMP revealed 72 under-expressed and 62 over-expressed genes in corticosteroid poor responders compared to good responders (Table S1). In SARP, 66 genes were under-expressed and 112 genes were over-expressed in corticosteroid poor responders compared to good responders (Table S2).

3.3. Connectivity Map Analysis

Next, the differentially expressed genes in CAMP and SARP were separately analyzed in CMap. Out of 2395 perturbagen compounds in the A549 cell line in CMap, 673 compounds had a negative connectivity score in CAMP, of which 277 compounds had a connectivity score less than −50 (Table S3). In SARP, 684 compounds had a negative connectivity score, and 261 compounds had a connectivity score less than −50 (Table S4). CAMP and SARP had 200 compounds in common with connectivity scores less than zero and 46 compounds in common with connectivity scores less than −50 (Table 2). Among the identified candidate compounds, eight were glucocorticoid receptor agonists, which is consistent with the direction of effect expected and a logical verification of our results. The compounds were prioritized based on weighted mean negative connectivity scores, and out of the top 10 drug compounds, γ-linolenic acid, ampicillin, exemestane, brinzolamide, and INCA-6 were selected for functional validation based on safety data in humans, mechanistic plausibility, and drug availability. Enrofloxacin, the compound with the highest weighted mean negative connectivity score, was not selected because it is a veterinary antibiotic. Due to limited data in humans and their conflicting mechanisms, serotonin receptor antagonist SB-203186 and serotonin receptor agonist dipropyl-5ct were not selected for functional validation.

3.4. In Vitro Validation

The five selected drug compounds were tested for their ability to modulate the anti-inflammatory effect of glucocorticoids using an assay based on the tethered transrepression of NF-κB activity by dexamethasone [33]. γ-linolenic acid, brinzolamide, and INCA-6 significantly reduced IL-1β induced luciferase activity in A549/NF-κB-luc reporter cells (Figure 1). Treatment with 30 μM of γ-linolenic acid led to a 20% relative reduction in luciferase activity after IL-1β stimulation compared to untreated cells (p = 0.004), and in the presence of dexamethasone, γ-linolenic acid treated cells had a significantly greater absolute reduction in luciferase activity (p = 0.004). Treatment with 30 μM of brinzolamide led to a 25% relative reduction in luciferase activity after IL-1β stimulation compared to untreated cells (p = 0.003), and the addition of dexamethasone led to a significantly greater absolute reduction in luciferase activity in brinzolamide treated cells (p = 0.013). Treatment with 30 μM of INCA-6 led to a 35% relative reduction in luciferase activity after IL-1β stimulation compared to untreated cells (p = 0.0001), and the addition of dexamethasone led to significantly greater absolute reduction in luciferase activity in INCA-6 treated cells (p = 0.002). Ampicillin and exemestane did not change IL-1β induced luciferase activity, and in the presence of dexamethasone, individual treatment with the highest concentrations of ampicillin and exemestane increased IL-1β induced luciferase activity compared to untreated samples (p = 0.0143 and p = 0.00003, respectively). Together, these data suggest that γ-linolenic acid, brinzolamide, and INCA-6 can each act as independent anti-inflammatory agents and synergistically augment the anti-inflammatory effects of corticosteroids.

4. Discussion

Poor response to corticosteroid treatment causes significant morbidity and mortality in asthma, and no treatments exist to directly alleviate corticosteroid resistance [34]. Few drug repurposing studies have been performed in asthma, and none have studied the repurposing of existing drugs to alleviate corticosteroid resistance [5,6,7,8,9]. This study is unique in its investigation of drugs that can be used to improve existing asthma therapies. Most prior studies used a drug-based approach by selecting specific compounds, including antithyroid and antihelminth agents, to repurpose as a new indication for the treatment of asthma, whereas this study used a disease-based drug repurposing approach to investigate the novel indication of glucocorticoid resistance in asthma [5,6,7,8,9].
In asthmatics treated with corticosteroids, we determined the intraindividual gene expression changes differentiating poor from good corticosteroid response as measured by FEV1 and identified 200 existing compounds that reversed the gene expression profile of poor responders to resemble that of good responders. These identified drug compounds have the potential to be repurposed as adjunctive treatments to improve corticosteroid response. Moreover, three of the top five candidate compounds were validated in in vitro human airway epithelial cells. γ-linolenic acid, brinzolamide, and INCA-6 each possessed anti-inflammatory activity independent of dexamethasone administration and further reduced IL-1β induced NF-κB activation in the presence dexamethasone.
γ-linolenic acid is an omega-6 fatty acid with anti-inflammatory properties found in plant-based oils such as borage oil, blackcurrant seed oil, and evening primrose oil [35,36]. As part of the polyunsaturated fatty acid biosynthesis pathway, γ-linolenic acid is elongated to form dihomo-γ-linolenic acid, which in turn is metabolized into anti-inflammatory series 1 prostaglandins, in particular prostaglandin E1 (PGE1), and 15-hydroxyeicosatrienoic acid (15-HETrE) by cyclooxygenase-1 and -2 and 15-lipoxygenase, respectively [35,36,37]. Paradoxically, dihomo-γ-linolenic acid can also undergo desaturation to form arachidonic acid, a precursor of many pro-inflammatory lipid mediators [36,37]. However, γ-linolenic acid is predominantly anti-inflammatory as both PGE1 and 15-HETrE metabolites inhibit the biosynthesis of arachidonic acid-derived eicosanoids [37]. In addition, dietary supplementation with γ-linolenic acid in asthmatics has been shown to reduce inflammation by attenuating leukotriene production and improve both quality of life and symptom control, and in agreement with these prior studies our results demonstrate the ability of γ-linolenic acid to potentiate the anti-inflammatory effect of glucocorticoids in asthma [38,39,40,41].
Brinzolamide is a carbonic anhydrase inhibitor used to treat elevated intraocular pressure in ocular hypertension and glaucoma. Carbonic anhydrase catalyzes the conversion of carbon dioxide and water to carbonic acid and also plays an important role in regulating type 2 inflammation. Carbonic anhydrase expression is induced by type 2 cytokines and elevated in allergen-induced pulmonary eosinophils [42]. Moreover, carbonic anhydrase positively regulates mast cell development, and inhibition of carbonic anhydrase selectively reduces human mast cell development and allergen-induced mast cell responses [43]. In addition, the inhibition of carbonic anhydrase by acetazolamide has been demonstrated to protect against bronchoconstriction induced by sodium metabisulfite, cold air, and deep inspiration in asthma [44,45,46]. These anti-inflammatory and protective effects of carbonic anhydrase inhibitors in asthma support our findings of its ability improve glucocorticoid response.
INCA-6 competitively binds to and inhibits calcineurin to prevent the transcriptional activation of cytokines, including IL-2, IL-4, IL-5, and granulocyte-macrophage colony-stimulating factor, leading to a reduction in both T lymphocyte proliferation and eosinophil proliferation and survival [47]. Calcineurin inhibitors also inhibit the release of histamine and the generation of arachadonic-acid derived eicosanoids from mast cells and basophils [48]. Resistance to glucocorticoids is a class effect caused by reduced inhibition of T lymphocyte proliferation by glucocorticoids and is not due to altered drug pharmacokinetics or metabolism [49]. Treatment with anti-T-lymphocyte agents including the calcineurin inhibitor cyclosporine has been shown to inhibit T cell proliferation in both glucocorticoid sensitive and resistance asthmatics [49]. Cyclosporine has also been found to improve lung function and reduce chronic oral corticosteroid dose and exacerbation rate [50,51]. However, the routine use of calcineurin inhibitors such as cyclosporine and tacrolimus in the treatment of asthma has been limited by adverse effects such as renal impairment, hypertension, hypertrichosis, headache, and peripheral neuropathy [50,51]. Calcineurin inhibitors are favorable candidates for drug repurposing to improve glucocorticoid sensitivity, but the selection of patients and optimal dosage needs to be carefully determined.
Included among the top 46 compounds identified on CMap analysis (Table 2) were both serotonin receptor agonists and antagonists, calcium channel activators and blockers, and adrenergic receptor agonists and antagonists. These antagonistic results suggest that serotonin, calcium, and adrenergic signaling have an important yet complex role in regulating corticosteroid response in asthmatics. Further research is needed to confirm the anti-inflammatory effects of these additionally identified drug classes. Contrary to expectations based on CMap, ampicillin and exemestane at the highest concentrations increased luciferase activity in A549/NF-κB-luc reporter cells in in vitro experiments. Ampicillin is a semisynthetic antibiotic derived from penicillin, and its routine use to treat asthma exacerbations has not been shown to accelerate improvement [52]. Exemestane is an aromatase inhibitor that blocks the conversion of androgens to estrogen. Estrogen promotes type 2 airway inflammation in asthma, which supports a plausible anti-inflammatory role for exemestane [53]. These pro-inflammatory effects of ampicillin and exemestane were concentration-dependent, and a different effect may have been detected at lower concentrations that were not tested. In addition, in vitro cell lines may not replicate true biological conditions and a different effect of these drugs may also be found in vivo. Further in vivo experiments are needed to determine the safety, efficacy, and optimal dosing of the identified compounds as adjunctive therapies to improve glucocorticoid response in asthma.
Genetic mutations affecting glucocorticoid resistance are only found in a small minority of asthmatics and do not fully explain the steroid resistance phenotype seen in asthma. Therefore, we focused on gene expression changes, which capture a wide variety of upstream processes, including environmental, genetics, and epigenetic causes that affect corticosteroid response. In addition, by examining relative instead of absolute corticosteroid resistance, our results have greater applicability across the spectrum of corticosteroid resistance in asthmatics. Strengths of this study include the inclusion of diverse cohorts, allowing for broad generalizability of our results across age, sex, ancestry, and asthma severity. Because of the heterogeneity between CAMP and SARP, differential gene expression and CMap analyses were performed separately in each cohort instead of a replication analysis, which produced two independent sets of negatively connected drug compounds. On direct comparison of these two independent sets of negatively connected drug compounds, 200 identical drug compounds were identified indicating the potential for these drugs to reverse poor corticosteroid response in asthmatics across a spectrum of age ranges, asthma severity, cell types, gene expression platforms, and corticosteroid response measurements. Furthermore, both the CMap analysis and in vitro validation were performed in A549 cells, which are representative of type II alveolar epithelial cells and suitable for studying drug delivery and metabolism in the pulmonary epithelium [54].
Important limitations to this study exist. CMap contains more than 2000 small-molecule compound perturbational signatures in A549 cells, and drug compounds not in the database were not investigated. CAMP used genotyping microarray data, and SARP used RNA sequencing data. It is possible that the use of RNA sequencing in both populations in addition to larger sample sizes would have increased the discovery of genes associated with steroid treatment response and hence, CMap hits. In addition, we only considered FEV1 response as an index of steroid nonresponsiveness. Previous work from our group has shown that steroid nonresponsiveness can present as other phenotypes not evaluated here [55]. Finally, although a glucocorticoid-resistant in vitro model was not used, the three compounds (γ-linolenic acid, brinzolamide, and INCA-6) that synergistically augmented the anti-inflammatory effect of dexamethasone also demonstrated anti-inflammatory activity at high concentrations independent of glucocorticoid administration and response. None of these theoretical limitations invalidate our overall approach and study design.

5. Conclusions

In summary, we identified several existing drugs with anti-inflammatory properties that may potentially alleviate glucocorticoid resistance in asthma. Experimental validation confirmed our results by demonstrating the potentiation of the effect of dexamethasone in pulmonary epithelial cells by γ-linolenic acid, brinzolamide, and INCA-6. Future in vivo studies are needed before these drugs can be approved as adjunctive therapies for asthma. Drug repurposing is timely and cost effective and has tremendous promise in treating poor glucocorticoid response. This novel approach can be generalized to treat resistance to standard therapies in diseases beyond asthma.

Supplementary Materials

The following are available online at https://www.mdpi.com/2075-4426/11/3/175/s1, Figure S1: Change in FEV1 percent predicted in CAMP and SARP, Table S1: Genes in CAMP differentially expressed in corticosteroid poor responders compared to corticosteroid good responders, Table S2: Genes in SARP differentially expressed in corticosteroid poor responders compared to corticosteroid good responders, Table S3: Connectivity Map compounds with negative connectivity scores less than −50 in CAMP, Table S4: Connectivity Map compounds with negative connectivity scores less than −50 in SARP.

Author Contributions

Conceptualization, K.G.T.; methodology, A.L.W., R.P., W.Q., A.T.K., Q.L. and K.G.T.; validation, R.P. and Q.L.; formal analysis, A.L.W., W.Q., A.T.K., G.C. and K.G.T.; resources, S.T.W., K.G.T., D.A.M. and E.R.B.; writing—original draft preparation, A.L.W. and W.Q.; writing—review and editing, all authors; funding acquisition, A.L.W., S.T.W., K.G.T., D.A.M., E.R.B. and Q.L. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the NIH, grant numbers R01 HL127332, R01 HL129935, U01 HL65899, R01 NR13391, P01 HL132825, R01 HL139496, R01 HL125734, and the Thrasher Research Fund Award #15115.

Institutional Review Board Statement

The study was conducted according to the guidelines of the Declaration of Helsinki, and approved by the Institutional Review Board of Mass General Brigham (1999P0015492) and Wake Forest University (00021507).

Informed Consent Statement

Informed consent was obtained from all subjects involved in the study.

Data Availability Statement

The data will be made publicly available on NCBI GEO upon publication.

Conflicts of Interest

The authors declare no conflict of interest. The funders had no role in the design of the study; in the collection, analyses, or interpretation of data; in the writing of the manuscript, or in the decision to publish the results.

Appendix A

SARP III RNA Sequencing: cDNA library was prepared with the Illumina TruSeq Stranded Total RNA (San Diego, CA, USA), quantified using qPCR and the PicoGreen dsDNA Assay Kit (Thermo Fisher Scientific, Waltham, MA, USA), and assessed for size and concentration using the 2100 Bioanalyzer System and 2200 TapeStation (Agilent, Santa Clara, CA, USA). Sequencing was performed on the Illumina HiSeq 2500 System with an average depth of 20 million reads. STAR was used to align the raw reads in sample FASTQ files to reference human genome version GRCh38 using default parameters after adapter sequence trimming by Skewer [56,57]. Raw counts for each mapped gene transcript were generated using the Bioconductor package Rsubread [58]. The counts were normalized using the Bioconductor package DESeq2 and log2-transformed [28]. A total of 65,988 gene transcripts in 38 samples (19 subjects) were generated. Gene transcripts with zero expression in all subjects were removed; 26,422 gene transcripts across both timepoints were retained. The data will be made publicly available on NCBI GEO upon publication.

References

  1. GBD 2016 Disease and Injury Incidence and Prevalence Collaborators. Global, regional, and national incidence, prevalence, and years lived with disability for 328 diseases and injuries for 195 countries, 1990–2016: A systematic analysis for the Global Burden of Disease Study 2016. Lancet 2017, 390, 1211–1259. [Google Scholar] [CrossRef] [Green Version]
  2. Szefler, S.J.; Martin, R.J.; King, T.S.; Boushey, H.A.; Cherniack, R.M.; Chinchilli, V.M.; Craig, T.J.; Dolovich, M.; Drazen, J.M.; Fagan, J.K.; et al. Significant variability in response to inhaled corticosteroids for persistent asthma. J. Allergy Clin. Immunol. 2002, 109, 410–418. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  3. Tantisira, K.G.; Lake, S.; Silverman, E.S.; Palmer, L.J.; Lazarus, R.; Silverman, E.K.; Liggett, S.B.; Gelfand, E.W.; Rosenwasser, L.J.; Richter, B.; et al. Corticosteroid pharmacogenetics: Association of sequence variants in CRHR1 with improved lung function in asthmatics treated with inhaled corticosteroids. Hum. Mol. Genet. 2004, 13, 1353–1359. [Google Scholar] [CrossRef] [Green Version]
  4. Jang, A.S. Steroid response in refractory asthmatics. Korean J. Intern. Med. 2012, 27, 143–148. [Google Scholar] [CrossRef]
  5. Suzuki, S.; Ogawa, M.; Ohta, S.; Arima, K.; Nunomura, S.; Nanri, Y.; Mitamura, Y.; Yoshihara, T.; Nakamura, Y.; Yamauchi, K.; et al. The potential for repositioning antithyroid agents as antiasthma drugs. J. Allergy Clin. Immunol. 2016, 138, 1458–1461.e1458. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  6. Miner, K.; Labitzke, K.; Liu, B.; Wang, P.; Henckels, K.; Gaida, K.; Elliott, R.; Chen, J.J.; Liu, L.; Leith, A.; et al. Drug Repurposing: The Anthelmintics Niclosamide and Nitazoxanide Are Potent TMEM16A Antagonists That Fully Bronchodilate Airways. Front. Pharm. 2019, 10, 51. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  7. Shin, E.; Lee, Y.C.; Kim, S.R.; Kim, S.H.; Park, J. Drug Signature-based Finding of Additional Clinical Use of LC28-0126 for Neutrophilic Bronchial Asthma. Sci. Rep. 2015, 5, 17784. [Google Scholar] [CrossRef] [Green Version]
  8. Al-Khami, A.A.; Ghonim, M.A.; Del Valle, L.; Ibba, S.V.; Zheng, L.; Pyakurel, K.; Okpechi, S.C.; Garay, J.; Wyczechowska, D.; Sanchez-Pino, M.D.; et al. Fuelling the mechanisms of asthma: Increased fatty acid oxidation in inflammatory immune cells may represent a novel therapeutic target. Clin. Exp. Allergy 2017, 47, 1170–1184. [Google Scholar] [CrossRef]
  9. Howell, C.; Smith, J.R.; Shute, J.K. Targeting matrix metalloproteinase-13 in bronchial epithelial repair. Clin. Exp. Allergy 2018, 48, 1214–1221. [Google Scholar] [CrossRef] [Green Version]
  10. Keskin, O.; Farzan, N.; Birben, E.; Akel, H.; Karaaslan, C.; Maitland-van der Zee, A.H.; Wechsler, M.E.; Vijverberg, S.J.; Kalayci, O. Genetic associations of the response to inhaled corticosteroids in asthma: A systematic review. Clin. Transl. Allergy 2019, 9, 2. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  11. Himes, B.E.; Jiang, X.; Wagner, P.; Hu, R.; Wang, Q.; Klanderman, B.; Whitaker, R.M.; Duan, Q.; Lasky-Su, J.; Nikolos, C.; et al. RNA-Seq transcriptome profiling identifies CRISPLD2 as a glucocorticoid responsive gene that modulates cytokine function in airway smooth muscle cells. PLoS ONE 2014, 9, e99625. [Google Scholar] [CrossRef]
  12. Qiu, W.; Guo, F.; Glass, K.; Yuan, G.C.; Quackenbush, J.; Zhou, X.; Tantisira, K.G. Differential connectivity of gene regulatory networks distinguishes corticosteroid response in asthma. J. Allergy Clin. Immunol. 2018, 141, 1250–1258. [Google Scholar] [CrossRef] [Green Version]
  13. Nosengo, N. Can you teach old drugs new tricks? Nature 2016, 534, 314–316. [Google Scholar] [CrossRef] [PubMed]
  14. Hay, M.; Thomas, D.W.; Craighead, J.L.; Economides, C.; Rosenthal, J. Clinical development success rates for investigational drugs. Nat. Biotechnol. 2014, 32, 40–51. [Google Scholar] [CrossRef] [PubMed]
  15. Austin, D.H. A Congressional Budget Office Study: Research and Development in the Pharmaceutical Industry; No. 2589; Congress of the United States: Washington, DC, USA, 2006; p. 55.
  16. Lamb, J.; Crawford, E.D.; Peck, D.; Modell, J.W.; Blat, I.C.; Wrobel, M.J.; Lerner, J.; Brunet, J.P.; Subramanian, A.; Ross, K.N.; et al. The Connectivity Map: Using gene-expression signatures to connect small molecules, genes, and disease. Science 2006, 313, 1929–1935. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  17. Subramanian, A.; Narayan, R.; Corsello, S.M.; Peck, D.D.; Natoli, T.E.; Lu, X.; Gould, J.; Davis, J.F.; Tubelli, A.A.; Asiedu, J.K.; et al. A Next Generation Connectivity Map: L1000 Platform and the First 1,000,000 Profiles. Cell 2017, 171, 1437–1452.e1417. [Google Scholar] [CrossRef]
  18. So, H.C.; Chau, C.K.; Chiu, W.T.; Ho, K.S.; Lo, C.P.; Yim, S.H.; Sham, P.C. Analysis of genome-wide association data highlights candidates for drug repositioning in psychiatry. Nat. Neurosci. 2017, 20, 1342–1349. [Google Scholar] [CrossRef]
  19. Nygren, P.; Fryknäs, M.; Agerup, B.; Larsson, R. Repositioning of the anthelmintic drug mebendazole for the treatment for colon cancer. J. Cancer Res. Clin. Oncol. 2013, 139, 2133–2140. [Google Scholar] [CrossRef] [Green Version]
  20. Xiao, Z.X.; Chen, R.Q.; Hu, D.X.; Xie, X.Q.; Yu, S.B.; Chen, X.Q. Identification of repaglinide as a therapeutic drug for glioblastoma multiforme. Biochem. Biophys. Res. Commun. 2017, 488, 33–39. [Google Scholar] [CrossRef]
  21. Kwon, O.S.; Lee, H.; Kong, H.J.; Kwon, E.J.; Park, J.E.; Lee, W.; Kang, S.; Kim, M.; Kim, W.; Cha, H.J. Connectivity map-based drug repositioning of bortezomib to reverse the metastatic effect of GALNT14 in lung cancer. Oncogene 2020. [Google Scholar] [CrossRef]
  22. Zhang, M.; Luo, H.; Xi, Z.; Rogaeva, E. Drug repositioning for diabetes based on ’omics’ data mining. PLoS ONE 2015, 10, e0126082. [Google Scholar] [CrossRef]
  23. Vargas, D.M.; De Bastiani, M.A.; Zimmer, E.R.; Klamt, F. Alzheimer’s disease master regulators analysis: Search for potential molecular targets and drug repositioning candidates. Alzheimers Res. 2018, 10, 59. [Google Scholar] [CrossRef]
  24. Childhood Asthma Management Program Research Group. The Childhood Asthma Management Program (CAMP): Design, rationale, and methods. Control Clin. Trials 1999, 20, 91–120. [Google Scholar] [CrossRef]
  25. Qiu, W.; Rogers, A.J.; Damask, A.; Raby, B.A.; Klanderman, B.J.; Duan, Q.L.; Tyagi, S.; Niu, S.; Anderson, C.; Cahir-Mcfarland, E.; et al. Pharmacogenomics: Novel loci identification via integrating gene differential analysis and eQTL analysis. Hum. Mol. Genet. 2014, 23, 5017–5024. [Google Scholar] [CrossRef] [Green Version]
  26. Teague, W.G.; Phillips, B.R.; Fahy, J.V.; Wenzel, S.E.; Fitzpatrick, A.M.; Moore, W.C.; Hastie, A.T.; Bleecker, E.R.; Meyers, D.A.; Peters, S.P.; et al. Baseline Features of the Severe Asthma Research Program (SARP III) Cohort: Differences with Age. J. Allergy Clin. Immunol. Pract. 2018, 6, 545–554.e544. [Google Scholar] [CrossRef]
  27. Phipatanakul, W.; Mauger, D.T.; Sorkness, R.L.; Gaffin, J.M.; Holguin, F.; Woodruff, P.G.; Ly, N.P.; Bacharier, L.B.; Bhakta, N.R.; Moore, W.C.; et al. Effects of Age and Disease Severity on Systemic Corticosteroid Responses in Asthma. Am. J. Respir. Crit. Care Med. 2017, 195, 1439–1448. [Google Scholar] [CrossRef] [PubMed]
  28. Love, M.I.; Huber, W.; Anders, S. Moderated estimation of fold change and dispersion for RNA-seq data with DESeq2. Genome Biol. 2014, 15, 550. [Google Scholar] [CrossRef] [Green Version]
  29. Ritchie, M.E.; Phipson, B.; Wu, D.; Hu, Y.; Law, C.W.; Shi, W.; Smyth, G.K. limma powers differential expression analyses for RNA-sequencing and microarray studies. Nucleic Acids Res. 2015, 43, e47. [Google Scholar] [CrossRef] [PubMed]
  30. Qiu, W.; Guo, B.; Anderson, C.; Klanderman, B.; Carey, V.; Raby, B. iCheck: QC Pipeline and Data Analysis Tools for High-Dimensional Illumina mRNA Expression Data; R Package Version 1.18.0; 2020. [Google Scholar]
  31. Carlson, M. org.Hs.eg.db: Genome Wide Annotation for Human; R Package Version 3.8.2; 2019. [Google Scholar]
  32. Iordanidou, M.; Paraskakis, E.; Tavridou, A.; Paschou, P.; Chatzimichael, A.; Manolopoulos, V.G. G894T polymorphism of eNOS gene is a predictor of response to combination of inhaled corticosteroids with long-lasting beta2-agonists in asthmatic children. Pharmacogenomics 2012, 13, 1363–1372. [Google Scholar] [CrossRef]
  33. Jiang, X.; Dahlin, A.; Weiss, S.T.; Tantisira, K.; Lu, Q. A high-throughput chemical screen identifies novel inhibitors and enhancers of anti-inflammatory functions of the glucocorticoid receptor. Sci. Rep. 2017, 7, 7405. [Google Scholar] [CrossRef] [Green Version]
  34. O’Byrne, P.; Fabbri, L.M.; Pavord, I.D.; Papi, A.; Petruzzelli, S.; Lange, P. Asthma progression and mortality: The role of inhaled corticosteroids. Eur. Respir. J. 2019, 54. [Google Scholar] [CrossRef] [Green Version]
  35. Innes, J.K.; Calder, P.C. Omega-6 fatty acids and inflammation. Prostaglandins Leukot. Essent. Fat. Acids 2018, 132, 41–48. [Google Scholar] [CrossRef] [Green Version]
  36. Timoszuk, M.; Bielawska, K.; Skrzydlewska, E. Evening Primrose (Oenothera biennis) Biological Activity Dependent on Chemical Composition. Antioxidants 2018, 7, 108. [Google Scholar] [CrossRef] [Green Version]
  37. Sergeant, S.; Rahbar, E.; Chilton, F.H. Gamma-linolenic acid, Dihommo-gamma linolenic, Eicosanoids and Inflammatory Processes. Eur. J. Pharm. 2016, 785, 77–86. [Google Scholar] [CrossRef] [Green Version]
  38. Arm, J.P.; Boyce, J.A.; Wang, L.; Chhay, H.; Zahid, M.; Patil, V.; Govindarajulu, U.; Ivester, P.; Weaver, K.L.; Sergeant, S.; et al. Impact of botanical oils on polyunsaturated fatty acid metabolism and leukotriene generation in mild asthmatics. Lipids Health Dis. 2013, 12, 141. [Google Scholar] [CrossRef] [Green Version]
  39. Surette, M.E.; Koumenis, I.L.; Edens, M.B.; Tramposch, K.M.; Clayton, B.; Bowton, D.; Chilton, F.H. Inhibition of leukotriene biosynthesis by a novel dietary fatty acid formulation in patients with atopic asthma: A randomized, placebo-controlled, parallel-group, prospective trial. Clin. Ther. 2003, 25, 972–979. [Google Scholar] [CrossRef]
  40. Surette, M.E.; Stull, D.; Lindemann, J. The impact of a medical food containing gammalinolenic and eicosapentaenoic acids on asthma management and the quality of life of adult asthma patients. Curr. Med. Res. Opin. 2008, 24, 559–567. [Google Scholar] [CrossRef] [PubMed]
  41. Lindemann, J.; David Pampe, E.; Peterkin, J.J.; Orozco-Cronin, P.; Belofsky, G.; Stull, D. Clinical study of the effects on asthma-related QOL and asthma management of a medical food in adult asthma patients. Curr. Med. Res. Opin. 2009, 25, 2865–2875. [Google Scholar] [CrossRef] [PubMed]
  42. Teh, A.L.; Pan, H.; Chen, L.; Ong, M.L.; Dogra, S.; Wong, J.; MacIsaac, J.L.; Mah, S.M.; McEwen, L.M.; Saw, S.M.; et al. The effect of genotype and in utero environment on interindividual variation in neonate DNA methylomes. Genome Res. 2014, 24, 1064–1074. [Google Scholar] [CrossRef] [Green Version]
  43. Henry, E.K.; Sy, C.B.; Inclan-Rico, J.M.; Espinosa, V.; Ghanny, S.S.; Dwyer, D.F.; Soteropoulos, P.; Rivera, A.; Siracusa, M.C. Carbonic anhydrase enzymes regulate mast cell-mediated inflammation. J. Exp. Med. 2016, 213, 1663–1673. [Google Scholar] [CrossRef] [Green Version]
  44. O’Connor, B.J.; Yeo, C.T.; Chen-Worsdell, Y.M.; Barnes, P.J.; Chung, K.F. Effect of acetazolamide and amiloride against sodium metabisulphite-induced bronchoconstriction in mild asthma. Thorax 1994, 49, 1096–1098. [Google Scholar] [CrossRef] [Green Version]
  45. O’Donnell, W.J.; Rosenberg, M.; Niven, R.W.; Drazen, J.M.; Israel, E. Acetazolamide and furosemide attenuate asthma induced by hyperventilation of cold, dry air. Am. Rev. Respir. Dis. 1992, 146, 1518–1523. [Google Scholar] [CrossRef]
  46. Spicuzza, L.; Ciancio, N.; Pellegrino, R.; Bellofiore, S.; Polosa, R.; Ricciardolo, F.L.; Brusasco, V.; Di Maria, G.U. The effect of inhaled furosemide and acetazolamide on bronchoconstriction induced by deep inspiration in asthma. Monaldi. Arch. Chest Dis. 2003, 59, 150–154. [Google Scholar] [PubMed]
  47. Sano, T.; Nakamura, Y.; Matsunaga, Y.; Takahashi, T.; Azuma, M.; Okano, Y.; Shimizu, E.; Ogushi, F.; Sone, S.; Ogura, T. FK506 and cyclosporin A inhibit granulocyte/macrophage colony-stimulating factor production by mononuclear cells in asthma. Eur. Respir. J. 1995, 8, 1473–1478. [Google Scholar] [PubMed]
  48. Harrison, C.A.; Bastan, R.; Peirce, M.J.; Munday, M.R.; Peachell, P.T. Role of calcineurin in the regulation of human lung mast cell and basophil function by cyclosporine and FK506. Br. J. Pharm. 2007, 150, 509–518. [Google Scholar] [CrossRef] [PubMed]
  49. Corrigan, C.J.; Bungre, J.K.; Assoufi, B.; Cooper, A.E.; Seddon, H.; Kay, A.B. Glucocorticoid resistant asthma: T-lymphocyte steroid metabolism and sensitivity to glucocorticoids and immunosuppressive agents. Eur. Respir. J. 1996, 9, 2077–2086. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  50. Alexander, A.G.; Kay, A.B.; Barnes, N.C. Trial of cyclosporin in corticosteroid-dependent chronic severe asthma. Lancet 1992, 339, 324–328. [Google Scholar] [CrossRef]
  51. Lock, S.H.; Kay, A.B.; Barnes, N.C. Double-blind, placebo-controlled study of cyclosporin A as a corticosteroid-sparing agent in corticosteroid-dependent asthma. Am. J. Respir. Crit. Care Med. 1996, 153, 509–514. [Google Scholar] [CrossRef]
  52. Normansell, R.; Sayer, B.; Waterson, S.; Dennett, E.J.; Del Forno, M.; Dunleavy, A. Antibiotics for exacerbations of asthma. Cochrane Database Syst. Rev. 2018, 6, Cd002741. [Google Scholar] [CrossRef]
  53. Yung, J.A.; Fuseini, H.; Newcomb, D.C. Hormones, sex, and asthma. Ann. Allergy Asthma. Immunol. 2018, 120, 488–494. [Google Scholar] [CrossRef] [Green Version]
  54. Foster, K.A.; Oster, C.G.; Mayer, M.M.; Avery, M.L.; Audus, K.L. Characterization of the A549 cell line as a type II pulmonary epithelial cell model for drug metabolism. Exp. Cell Res. 1998, 243, 359–366. [Google Scholar] [CrossRef] [PubMed]
  55. Rogers, A.J.; Tantisira, K.G.; Fuhlbrigge, A.L.; Litonjua, A.A.; Lasky-Su, J.A.; Szefler, S.J.; Strunk, R.C.; Zeiger, R.S.; Weiss, S.T. Predictors of poor response during asthma therapy differ with definition of outcome. Pharmacogenomics 2009, 10, 1231–1242. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  56. Dobin, A.; Davis, C.A.; Schlesinger, F.; Drenkow, J.; Zaleski, C.; Jha, S.; Batut, P.; Chaisson, M.; Gingeras, T.R. STAR: Ultrafast universal RNA-seq aligner. Bioinformatics 2013, 29, 15–21. [Google Scholar] [CrossRef] [PubMed]
  57. Jiang, H.; Lei, R.; Ding, S.-W.; Zhu, S. Skewer: A fast and accurate adapter trimmer for next-generation sequencing paired-end reads. BMC Bioinform. 2014, 15, 182. [Google Scholar] [CrossRef]
  58. Liao, Y.; Smyth, G.K.; Shi, W. The R package Rsubread is easier, faster, cheaper and better for alignment and quantification of RNA sequencing reads. Nucleic Acids Res. 2019, 47, e47. [Google Scholar] [CrossRef] [Green Version]
Jpm 11 00175 g001a 550Jpm 11 00175 g001b 550
Figure 1. Effect of the selected drug compounds on the response of NF-κB reporter cells to IL-1β ± dexamethasone (dex). In vitro validation of the top five drug compounds from the Connectivity Map analysis in the Childhood Asthma Management Program and Severe Asthma Research Program III: (a) γ-linolenic-acid, (b) ampicillin, (c) exemestane, (d) brinzolamide, and (e) INCA-6. The luciferase activities of reporter cells treated with the drug compounds were compared to that of untreated cells. Data represent means ± standard errors compared using t-tests. * p < 0.05 (drug compounds vs. untreated in cells exposed to IL-1β alone); # p < 0.05 (drug compounds vs. untreated in cells exposed to IL-1β ± dex).
Figure 1. Effect of the selected drug compounds on the response of NF-κB reporter cells to IL-1β ± dexamethasone (dex). In vitro validation of the top five drug compounds from the Connectivity Map analysis in the Childhood Asthma Management Program and Severe Asthma Research Program III: (a) γ-linolenic-acid, (b) ampicillin, (c) exemestane, (d) brinzolamide, and (e) INCA-6. The luciferase activities of reporter cells treated with the drug compounds were compared to that of untreated cells. Data represent means ± standard errors compared using t-tests. * p < 0.05 (drug compounds vs. untreated in cells exposed to IL-1β alone); # p < 0.05 (drug compounds vs. untreated in cells exposed to IL-1β ± dex).
Jpm 11 00175 g001aJpm 11 00175 g001b
Table
Table 1. Demographics of Childhood Asthma Management Program (CAMP) and Severe Asthma Research Program (SARP).
Table 1. Demographics of Childhood Asthma Management Program (CAMP) and Severe Asthma Research Program (SARP).
CAMPSARP
Good Responders
(n = 47)		Poor Responders
(n = 48)		p-ValueGood Responders
(n = 8)		Poor Responders
(n = 11)		p-Value
Age (mean ± SD)9.0 (±2.2)8.5 (±2.1)0.350.9 (±9.1)43.3 (±6.6)0.05
Sex 0.61 1.00
  Female25 (53.2%)23 (47.9%)7 (87.5%)9 (81.8%)
  Male22 (46.8%)25 (52.1%)1 (12.5%)2 (18.2%)
Ancestry 0.6 1.00
  European39 (83.0%)42 (87.5%)5 (62.5%)8 (72.7%)
  African5 (10.6%)4 (8.3%)3 (37.5%)3 (27.3%)
  Hispanic1 (2.1%)2 (4.2%)--
  Other2 (4.3%)0 (0.0%)--
Asthma severity 0.61 0.10
  Mild23 (48.9%)21 (43.8%)3 (37.5%)1 (9.1%)
  Moderate24 (51.1%)27 (56.2%)2 (25.0%)1 (9.1%)
  Severe--3 (37.5%)9 (81.8%)
Body mass index (mean ± SD)18.1 (±2.8)18.0 (±2.9)0.934.9 (±8.2)35.3 (±8.3)0.90
Baseline serum IgE level
(kU/L, median ± IQR)		230.0 (±310.0)150.0 (±236.0)0.3205.0 (±396.0)36.0 (±195.0)0.09
Baseline serum eosinophil count
(cells/µL, median ± IQR)		493.0 (±481.0)390.0 (±400.0)0.1244.0 (±220.0)172.0 (±77.0)0.08
Baseline PC20 (median ± IQR)0.6 (±1.0)1.1 (±1.7)0.050.4 (±0.0)1.5 (±0.8)0.50
Pre-treatment FEV1 % predicted
(mean ± SD)		86.0% (±15.2%)99.8% (±14.1%)<0.00172.5% (±22.5%)82.7% (±11.9%)0.20
Post-treatment FEV1 % predicted
(mean ± SD)		101.0% (±11.5%)100.5% (±15.1%)0.876.5% (±23.6%)73.6% (±17.9%)0.80
Pre-treatment FEV1/FVC
(mean ± SD)		75.1% (±9.9%)80.8% (±7.4%)0.00270.4% (±11.0%)75.8% (±8.6%)0.20
Post-treatment FEV1/FVC
(mean ± SD)		82.0% (±6.5%)81.5% (±8.3%)0.771.6% (±9.2%)73.0% (±9.2%)0.80
Table
Table 2. Mean negative connectivity scores in CAMP and SARP weighted by sample sizes.
Table 2. Mean negative connectivity scores in CAMP and SARP weighted by sample sizes.
NameDescriptionCAMP (n = 95)
Connectivity Score 		SARP (n = 19)
Connectivity Score 		Weighted Mean
Connectivity Score
enrofloxacinBacterial DNA gyrase inhibitor−96.51−99.79−97.06
SB-203186Serotonin receptor antagonist−93.70−91.12−93.27
γ-linolenic-acidAnti-inflammatory omega-6 fatty acid−94.78−80.70−92.43
dipropyl-5ct (3-(N,N-Dipropylaminoethyl)-1H-indole-5-carboxamide maleate)Serotonin receptor agonist−92.65−83.03−91.05
ampicillinBacterial cell wall synthesis inhibitor−97.60−56.59−90.76
exemestaneAromatase inhibitor−86.59−98.24−88.53
brinzolamideCarbonic anhydrase inhibitor−86.81−89.74−87.30
INCA-6Calcineurin inhibitor−87.45−78.14−85.90
SCH-23390Dopamine receptor antagonist−88.17−72.69−85.59
brazilinNitric oxide production inhibitor−86.00−80.72−85.12
pyrazinamideFatty acid synthase inhibitor−82.85−90.52−84.13
vinburnineAdrenergic receptor antagonist−79.35−99.25−82.67
NS-1619Calcium channel activator−79.38−98.59−82.58
caffeineAdenosine receptor antagonist−78.83−81.33−79.25
masitinibKIT inhibitor−79.32−76.71−78.88
isoflupredoneGlucocorticoid receptor agonist−82.06−62.39−78.78
fluocinonideGlucocorticoid receptor agonist−80.62−66.33−78.24
fluoxetineSelective serotonin reuptake inhibitor−74.30−96.36−77.98
quinpiroleDopamine receptor agonist−75.20−90.40−77.73
saquinavirHIV protease inhibitor−72.31−94.39−75.99
retinolRetinoid receptor ligand−72.33−80.84−73.75
doxapramPotassium channel blocker−74.67−61.12−72.41
somatostatinSomatostatin receptor agonist−68.97−81.93−71.13
thalidomideTNF production inhibitor−73.74−57.45−71.02
L-690488Inositol monophosphatase inhibitor−65.92−86.80−69.40
fludroxycortideGlucocorticoid receptor agonist−69.60−58.68−67.78
dexamethasoneGlucocorticoid receptor agonist−63.63−86.38−67.42
nitrendipineCalcium channel blocker−70.41−52.10−67.36
acyclovirDNA polymerase inhibitor−64.90−71.25−65.96
DMAB-anabaseineAdrenergic receptor agonist−68.48−51.34−65.62
m-chlorophenylbiguanideSerotonin receptor agonist−61.73−82.67−65.22
SB-590885RAF inhibitor−65.27−59.71−64.34
trimebutineOpioid receptor agonist−58.98−84.10−63.17
teicoplaninBacterial cell wall synthesis inhibitor−64.10−56.12−62.77
PD-102807Acetylcholine receptor antagonist−61.37−61.17−61.34
benzatropineAcetylcholine receptor antagonist−55.99−81.83−60.30
GBR-13069Dopamine uptake inhibitor−57.34−70.82−59.59
vinblastineMicrotubule inhibitor−55.40−76.83−58.97
sulindacCyclooxygenase inhibitor−58.98−55.76−58.44
dexketoprofenCyclooxygenase inhibitor−59.67−50.92−58.21
hydrocortisoneGlucocorticoid receptor agonist−50.67−91.73−57.51
fludrocortisoneGlucocorticoid receptor agonist−50.42−92.94−57.51
OMDM-2FAAH inhibitor−57.04−58.07−57.21
diflorasoneCorticosteroid agonist−54.44−70.40−57.10
beclometasoneGlucocorticoid receptor agonist−52.13−72.93−55.60
marbofloxacinBacterial DNA gyrase inhibitor−55.04−51.70−54.48
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Wang, A.L.; Panganiban, R.; Qiu, W.; Kho, A.T.; Chupp, G.; Meyers, D.A.; Bleecker, E.R.; Weiss, S.T.; Lu, Q.; Tantisira, K.G. Drug Repurposing to Treat Glucocorticoid Resistance in Asthma. J. Pers. Med. 2021, 11, 175. https://doi.org/10.3390/jpm11030175

AMA Style

Wang AL, Panganiban R, Qiu W, Kho AT, Chupp G, Meyers DA, Bleecker ER, Weiss ST, Lu Q, Tantisira KG. Drug Repurposing to Treat Glucocorticoid Resistance in Asthma. Journal of Personalized Medicine. 2021; 11(3):175. https://doi.org/10.3390/jpm11030175

Chicago/Turabian Style

Wang, Alberta L., Ronald Panganiban, Weiliang Qiu, Alvin T. Kho, Geoffrey Chupp, Deborah A. Meyers, Eugene R. Bleecker, Scott T. Weiss, Quan Lu, and Kelan G. Tantisira. 2021. "Drug Repurposing to Treat Glucocorticoid Resistance in Asthma" Journal of Personalized Medicine 11, no. 3: 175. https://doi.org/10.3390/jpm11030175

APA Style

Wang, A. L., Panganiban, R., Qiu, W., Kho, A. T., Chupp, G., Meyers, D. A., Bleecker, E. R., Weiss, S. T., Lu, Q., & Tantisira, K. G. (2021). Drug Repurposing to Treat Glucocorticoid Resistance in Asthma. Journal of Personalized Medicine, 11(3), 175. https://doi.org/10.3390/jpm11030175

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop