Impact Factor 3.845 | CiteScore 3.92
More on impact ›

Original Research ARTICLE

Front. Pharmacol., 14 February 2019 | https://doi.org/10.3389/fphar.2019.00051

Drug Repurposing: The Anthelmintics Niclosamide and Nitazoxanide Are Potent TMEM16A Antagonists That Fully Bronchodilate Airways

Kent Miner1, Katja Labitzke2†, Benxian Liu1, Paul Wang3, Kathryn Henckels1, Kevin Gaida1, Robin Elliott1, Jian Jeffrey Chen4, Longbin Liu4†, Anh Leith5†, Esther Trueblood6,7,8†, Kelly Hensley6,7,8†, Xing-Zhong Xia1, Oliver Homann9, Brian Bennett1, Mike Fiorino1, John Whoriskey1, Gang Yu1, Sabine Escobar5, Min Wong1, Teresa L. Born5†, Alison Budelsky5†, Mike Comeau5, Dirk Smith5†, Jonathan Phillips1, James A. Johnston1, Joseph G. McGivern3, Kerstin Weikl2†, David Powers3, Karl Kunzelmann10, Deanna Mohn1*, Andreas Hochheimer2*† and John K. Sullivan1*†
  • 1Department of Inflammation Research, Amgen Inc., Thousand Oaks, CA, United States
  • 2Department of Therapeutic Discovery, Amgen Inc., Regensburg, Germany
  • 3Department of Therapeutic Discovery, Amgen Inc., Thousand Oaks, CA, United States
  • 4Department of Medicinal Chemistry, Amgen Inc., Thousand Oaks, CA, United States
  • 5Department of Inflammation Research, Amgen Inc., Seattle, WA, United States
  • 6Department of Comparative Biology and Safety Sciences, Amgen Inc., Seattle, WA, United States
  • 7Department of Comparative Biology and Safety Sciences, Amgen Inc., Thousand Oaks, CA, United States
  • 8Department of Comparative Biology and Safety Sciences, Amgen Inc., South San Francisco, CA, United States
  • 9Genome Analysis Unit, Amgen Inc., South San Francisco, CA, United States
  • 10Institut für Physiologie, Universität Regensburg, Regensburg, Germany

There is an unmet need in severe asthma where approximately 40% of patients exhibit poor β-agonist responsiveness, suffer daily symptoms and show frequent exacerbations. Antagonists of the Ca2+-activated Cl channel, TMEM16A, offers a new mechanism to bronchodilate airways and block the multiple contractiles operating in severe disease. To identify TMEM16A antagonists we screened a library of ∼580,000 compounds. The anthelmintics niclosamide, nitazoxanide, and related compounds were identified as potent TMEM16A antagonists that blocked airway smooth muscle depolarization and contraction. To evaluate whether TMEM16A antagonists resist use- and inflammatory-desensitization pathways limiting β-agonist action, we tested their efficacy under harsh conditions using maximally contracted airways or airways pretreated with a cytokine cocktail. Stunningly, TMEM16A antagonists fully bronchodilated airways, while the β-agonist isoproterenol showed only partial effects. Thus, antagonists of TMEM16A and repositioning of niclosamide and nitazoxanide represent an important additional treatment for patients with severe asthma and COPD that is poorly controlled with existing therapies. It is of note that drug repurposing has also attracted wide interest in niclosamide and nitazoxanide as a new treatment for cancer and infectious disease. For the first time we identify TMEM16A as a molecular target for these drugs and thus provide fresh insights into their mechanism for the treatment of these disorders in addition to respiratory disease.

Introduction

Asthma is a heterogenous disease with multiple clinical phenotypes and newly described endotypes (Wenzel, 2012; Holgate et al., 2015), however, core features of airway constriction and hyper-responsiveness are likely common to multiple endotypes. Over 300 million people worldwide are afflicted with asthma, with 5–10% of patients suffering from severe disease which remains poorly controlled despite maximal treatment with an inhaled corticosteroid (ICS) and long-acting β-agonist (LABA). Inflammation is poorly controlled in severe asthmatics, which often necessitates addition of oral corticosteroids, and ∼40% of patients show poor β-agonist responsiveness (Moore et al., 2007, 2010). Agonists of the β2-adrenergic receptor (β2AR) signal through the stimulatory G-protein, Gsα, to activate adenylate cyclase (AC), increase cAMP levels and activate PKA. While this mechanism offers the advantage of delivering a negative intracellular signal (cAMP/PKA) to airway smooth muscle (ASM) cells that blocks multiple contractants, repeated β-agonist use and poorly controlled inflammation can cause desensitization of this pathway. Phosphodiesterases (PDE) can degrade cAMP, and increased levels of PDE4D have been detected in biopsies of asthmatic ASM cells leading to reduced β-agonist cAMP responses (Trian et al., 2011). Proinflammatory cytokines or allergic inflammation has been found to increase expression of the inhibitory G-protein, Giα (Hakonarson et al., 1996), uncouple the β2AR from Gsα-induced activation of adenylate cyclase (Shore et al., 1997) and upregulate COX-2 and PGE2 production inducing heterologous β2AR desensitization of ASM cells and reduced β-agonist responsiveness (Moore et al., 2001).

Ion channels represent an alternative target class to modulate contraction which may resist these desensitization pathways, but it has been poorly understood what channels control excitation-contraction coupling in ASM cells. In vascular and airway smooth muscle, bronchoconstrictors have been recognized for years to induce calcium-activated chloride channels (CaCCs) promoting chloride efflux and spontaneous transient inward currents (STICs) contributing to depolarization and contraction, but the molecular identity of the channel has remained elusive (Janssen and Sims, 1995; Large and Wang, 1996; Wellman and Nelson, 2003).

In 2008, three separate groups employing distinct methods identified TMEM16A as the long sought-after CaCC recognized in secretory epithelial cells, smooth muscle cells and sensory neurons (Caputo et al., 2008; Schroeder et al., 2008; Yang et al., 2008). While these initial studies solidified TMEM16A’s role in the bronchial epithelium as responsible for the elevated Ca2+-dependent chloride secretion induced by Th2 cytokines, subsequent studies by Huang et al. (2012) expanded this perspective suggesting TMEM16A may have a dual role in controlling both ASM contraction and epithelial mucin secretion, a hallmark of allergic asthma. Although this research, and more recent findings by Danielsson et al. (2015), introduced TMEM16A antagonists as a novel approach to bronchodilate airways, only a limited number of compounds were explored and details of their relative efficacy compared to β-agonists remained unclear.

To identify additional TMEM16A antagonists and characterize their potential as novel bronchodilators, we screened a library of over half a million compounds. For the first time, we describe the approved drug niclosamide, nitazoxanide, and related compounds as potent TMEM16A antagonists that provided robust bronchodilation of airways and resist use- and inflammatory-desensitization pathways limiting β-agonist action. We also discuss separate research from drug-repurposing screens finding the anthelmintics niclosamide and nitazoxanide have efficacy for treating a variety of other disorders, including cancer, and discuss this in the context of TMEM16A, which could represent a molecular target explaining some of these activities.

Materials and Methods

Materials

CaCCinh-A01, MONNA, and T16Ainh-A01 were purchased from Sigma-Aldrich and tizoxanide was purchased from J&W Pharmlab LLC. Niclosamide (Cpd 1) and niclosamide-related compounds Cpd 2, Cpd 3, Cpd 4, and Cpd 6 – 18, nitazoxanide and CID# 2806957 were from the Amgen small molecule compound collection, as were the benchmark antagonists niflumic acid, dichlorophen, benzbromarone, 1PBC and NTTP. Cpd 5 was synthesized according to a literature procedure (Macielag et al., 1998). The source for other compounds and reagents is provided in each experimental section.

Cell Line Generation and High-Throughput Screen

The High-Throughput-Screening clone HEK293T:eYFP: TMEM16A(abc) stably expressing eYFP (H149Q,I153L) and human TMEM16A(abc) was generated by stable transfection followed by functional clone selection. HEK293T:eYFP (H149Q,I153L) was cultivated in propagation medium DMEM (PAN #P04-03550) containing 10% FBS (PAN #3302-P281402), 1x Pen/Strep (PAN #P06-07100), 1x L-glutamine (PAN #P04-80100), 10 mM HEPES (PAN #P05-01100) and 0.5 μg/ml puromycin (Sigma #P7255) and transfected with pIREShyg3:huTMEM16A(abc). Stable cell pool was generated by cultivation in propagation medium containing 0.5 μg/ml puromycin and 200 μg/ml hygromycin B for 20 days. Cell pool was then diluted in the same medium and seeded in 384-well microtiter plates. Grown clones were replicated into daughter plates using a CyBio Felix pipetting platform (Analytik Jena AG). Daughter plates were analyzed using the eYFP assay essentially as described in the HTS protocol below using a FLIPR-Tetra. Clones showing maximal quenching of eYFP fluorescence intensity after addition of trigger solution (assay buffer containing NaI and ionomycin) were selected for propagation from the mother plate and further characterized with regard to HTS assay development performance parameters including inhibition by the benchmark compound benzbromarone, S/B and relevant HTS quality control statistics (Z-factor, Z′-factor, RZ-factor, RZ′-factor). One clone was finally selected as HTS clone and used throughout the entire HTS campaign to identify TMEM16A antagonists.

For the HTS screen, the screening cell line HEK293T:eYFP(H149Q,I153L):TMEM16A(abc) stably expressing halide sensitive YFP (eYFP (H149Q,I153L)) and human TMEM16A(abc) were propagated in DMEM (PAN #P04-03550) containing 10% FBS (PAN #3302-P281402), 1x Pen/Strep (PAN #P06-07100), 1x L-glutamine (PAN #P04-80100), 10 mM HEPES (PAN #P05-01100), 0.5 μg/ml puromycin (Sigma #P7255), and 150 μg/ml hygromycin B (Life Technologies #10687-010) in Corning TC-treated flasks. Twenty four hours prior to the assay 11,000 cells were seeded per well in 30 μl assay plating medium [propagation medium containing 0.01% Pluronic F-68 (PAN #P08-2100)] in Corning CellBIND 384-Well Flat Clear Bottom Black Polystyrene Microplates (Corning #3683) and cultivated for 24 h at 37°C. Before initiating the fully automated and scheduled HTS sequence, cell plates were transferred to an incubator and stored at RT and 0% CO2. Library compounds were delivered in Greiner 384 PP plates (Greiner #784201) containing 1mM compound dissolved in DMSO in columns 1–22 and DMSO in columns 23 and 24 and were stored in a stacker at RT. Further dilutions of compounds were generated by fully automated liquid handling using a CyBio 384-well/25 μl pipetting head (Analytik Jena AG) by diluting 0.5 μl compound (1 mM in DMSO) in 20 μl assay buffer [1xHBSS (Life Technologies #14025050) containing 10 mM HEPES (PAN #P05-01100) and 0.01% Pluronic F-68 (PAN #P08-2100), adjusted to pH 7.4] in Greiner 384 PP plates (Greiner #784201). Intermediate compound plates contained 25 μM library compounds in assay buffer containing 2.5% DMSO in columns 1–22, assay buffer containing 2.5% DMSO as neutral control in column 23 and 100 μM Benzbromarone (Sigma #B5774) in assay buffer containing 2.5% DMSO in column 24 as ligand control (benchmark antagonist) and were stored in a stacker at RT. Simultaneous to compound dilution medium was aspirated from cell plates using a Biotek 405 microplate washer followed by washing the cells twice with 65 μl/well assay buffer. With the final aspiration of assay buffer, the residual volume was adjusted to 15 μl. Ten μl compound solution were then added to the washed cell plates using a CyBio 384-well/25 μl pipetting head (Analytik Jena AG) resulting in a final volume of 25 μl assay buffer containing 10 μM compound and 1% DMSO. Cell plates were incubated for 30 min at RT and 0% CO2. After incubation cell plates were transferred to a FLIPR-Tetra and 25 μl trigger solution [assay buffer containing 20 mM sodium iodide (Sigma #383112) and 4 μM ionomycin (Sigma #I0634; stock solution 10 mM in DMSO)] was added before reading iodide-quenching of eYFP fluorescence intensity. Final concentration of library compounds in columns 1–22 at read-out was 5 μM in the presence of 10 mM iodide, 2 μM ionomycin and 0.5% DMSO. Neutral control column 23 received no ligand and the ligand control column received a final concentration 20 μM benzbromarone in the presence of 10 mM iodide, 2 μM ionomycin, which corresponded to EC80 in this assay setup, and 0.5% DMSO. FLIPR-Tetra settings included excitation wavelength (470–495 nm), emission wavelength (515–575 nm), 7 reads before dispensing the trigger solution for baseline recording followed by 40 reads after trigger solution dispense. Data analysis was performed using Genedata Screener with aggregation rule “Max-Min/Max” (Min: baseline from 0 to 10 s prior to trigger solution dispense) and normalization based on neutral control (0% inhibition corresponding to 100% TMEM16A activity) and inhibitor control Benzbromarone (100% inhibition corresponding to 0% TMEM16A activity). HTS quality control parameters and statistics (S/B, Z-factor, RZ-factor, Z′-factor and RZ′-factor) were automatically calculated and recorded using Genedata Screener. For unbiased primary hit selection, the POC cut-off (hit threshold) was calculated as Median + 3xIQR + 20% of all test wells. Confirmation Screening primary was run in triplicates (three independent, consecutive runs) with confirmed hits defined as median (replicate POC values) > POC cut-off (POC cut-off primary screen). Compounds were selected for dose-response screening after database mining to eliminate frequent hitters and fluorescent primary hit compounds, which were detected by calculating baseline fluorescence intensity measured by FLIPR-Tetra before addition of trigger solution for each compound well. Cut-off was set as Baseline FI < 10x STDV of the neutral control. Non-fluorescent compounds were selected and dose response testing was performed with 22-step, 1:2 dilution following the eYFP assay format procedure described above. Threshold was set at IC50 < 5 μM.

The 22 dose response dilutions were also tested in Calcium flux assays using FLIPR Calcium 6 Assay Kit (Molecular Devices # R8190 and R8191) to eliminate compounds, which inhibit ionomycin-dependent increase of intracellular calcium concentration. HEK293T cells were propagated in Corning TC-treated flasks in propagation medium consisting of DMEM (PAN #P04-03550) containing 10% FBS (PAN #3302-P281402), 1x Pen/Strep (PAN #P06-07100), 1x L-glutamine (PAN #P04-80100), 10 mM HEPES (PAN #P05-01100), 1x sodium-pyruvate (PAN #P04-43100). Prior to the assay 13,000 cells were plated per well in 30 μl propagation medium in Corning Cellbind 384 well microtiter plates (Corning #3683) and cultivated for 24 h. Cells were stained with labeling dye solution, which was assembled following a protocol provided by manufacturer followed by addition of 8 mM probenecid [Sigma-Aldrich #P8761, 250 mM stock solution in 1xHBSS (PAN #P04-49505) mixed with an equal volume of 1N NaOH and pH adjusted to 7.4]. Twenty-five μl labeling dye were transferred per well using Agilent Bravo liquid handling platform followed by incubation for 2 h at 37°C and 5% CO2. Simultaneously, intermediate compound plates were prepared by transferring 0.64 μl of compound solution from compound plates to 15 μl assay buffer [1xHBSS (Life Technologies #14025050), 20 mM HEPES (PAN #P05-01100), adjusted to pH 7.4] in columns 1–22 of intermediate plates using a CyBio 384/25 μl pipetting head (Analytik Jena AG) followed by storage in a stacker. In order to normalize the assay to identify compounds decreasing calcium signaling final assay concentration of 1 μM ionomycin was tested as neutral control in column 23 and buffer only as blank control in column 24. When cell labeling was complete, 10 μl of diluted compound were added per well using Agilent Bravo Platform and were incubated for 45 min at room temperature. Cell plates were then transferred to FLIPR-Tetra and 20 μl trigger solution [assay buffer containing 4.25 μM ionomycin (Sigma #I0634; stock solution 10 mM in DMSO)] was added to achieve a final Ionomycin concentration of 1 μM, which corresponded to the EC80 in this assay setup, in the presence of 0.5% DMSO (final). FLIPR-Tetra settings included excitation wavelength (470–495 nm), emission wavelength (515–575nm), baseline was recorded before dispensing the trigger solution followed by reading the ionomycin response. Data analysis was performed using Genedata Screener with aggregation rule “Max-Min/Max” (Min: baseline prior to trigger solution dispense) and normalization based on neutral control (0% inhibition) and blank control without ionomycin (100% inhibition). Dose response cut-off was set at IC50 < 10 μM to eliminate compounds, which reduce ionomycin-triggered calcium signaling.

Cell Line Generation and Medicinal Chemistry Support

The TMEM16A HEK293T eYFP (H149Q, I153L) stable cell line for medchem support was generated by transfecting 10 μg of linearized hTMEM16A(abc) DNA constructs into HEK293T eYFP (H149Q,I153L) stable cells. The transfected cells were cultured under 200 μg/ml of hygromycin B selection for 20 days and the stable pool was single cell seeded in 96 well plate and duplicated plate were made after single cell clones were formed. YFP function assay were performed on duplicated plate and clones producing higher changes in fluorescence from the YFP assay were chosen for further expansion to form a stable cell line.

For routine tests of compound activity, HEK293T cells stably co-expressing hTMEM16A (abc) and the halide-sensitive YFP were cultured with DMEM/F-12, HEPES (Life Technologies, Catalog Number: 11330-032) plus 10% (heat inactive) FBS (Life Technologies, Catalog Number: 10082-147), 1X Pen Strep Glutamine (Gibco 10378-016), 0.5 μg/ml of puromycin and 200 μg/ml Hygromycin B. The day before the assay, 30 μl/well of the cell’s suspension was seeded to achieve 15000 cells/well in 384 well plate (Corning® BioCoatTM Poly-D-Lysine 384 Well black/clear plate, Cat# 354663). Plates were incubated for 24 h at 37°C with 5% CO2. On the assay day, medium was removed and replaced with 20 μl/well of assay buffer (Hank’s Balanced Salt Solution with 10 mM HEPES pH 7.4), 10 μl/well of 40x diluted compounds was added to assay plate and it was incubated at room temperature for 30 min, after which baseline fluorescence was read in the FLIPR instrument for 10 s. Then, 5 μl/well of a 3x trigger solution (24 mM iodide, 3 μM ionomycin in assay buffer) was added and the fluorescence kinetic trace was recorded for 2 min.

Q-Patch Electrophysiology for Measuring Compound Effects on TMEM16A

HEK293 cells stably expressing the TMEM16A (acd) variant were purchased from SB Drug Discovery. The HEK293 TMEM16A (abc) stable cell line was from ChanTest and COLO-205 cells were purchased from ATCC. The same standard buffers and recording conditions listed below were used for QPatch studies on all three cell lines.

Recording Solutions: were made weekly as below and stored at room temperature until used. It is recommended not to add ATP to internal solution until immediately before using on the QPatch. The external buffer (in mM) was 140 NaCl, 4 KCl, 2 CaCl2, 1 MgCl2, 10 HEPES, 10 Glucose, pH7.4, while the internal buffer (in mM) was 110 CsCl, 20 TEA-Cl, 5.374 CaCl2, 1.75 MgCl2, 10 EGTA, 10 HEPES, 4 Na2ATP, pH 7.2.

Once established in whole-cell configuration, cells were clamped to a holding potential (Vhold) of 0 mV on 48 well single hole patch plate. A standardized IV protocol was used to elicit ionic current through the TMEM16A chloride channel at 20 s intervals. Steady-state voltage pulses were from −100 mV to +100 mV in +20 mV steps for a duration of 500 ms. After each pulse the voltage was returned to −100 mV for 50 ms to obtain tail currents, and then returned to the holding potential of 0 mV until the beginning of the next sweep. Block of TMEM16A chloride current due to test article was measured at the +100 mV depolarization sweep of the third and final IV protocol run per concentration addition. This ensures a minimum of 60 s per concentration incubation before measuring current block. Currents were acquired at 10 kHz and filtered at 2 kHz, leak subtraction and Rseries were disabled. A 1 MΩ minimum resistance was set because cells are held in an open channel state and therefore resistance is low.

A cumulative concentration response was measured where each cell is exposed to five concentrations of test article with a dilution factor of 1:5 (e.g., 0.048, 0.24, 1.2, 6, and 30 μM). Cells were recorded for ∼60 s per solution/compound addition. Initially, external solution is applied twice to allow currents to stabilize. Then vehicle (from the saline reservoir) is added twice to monitor any effect 0.3% DMSO might have on currents. This leads directly into the concentration runs by applying Concentration 1 (single addition) for 1 min. Concentration 2 is then applied for the subsequent minute, etc. Recovery/washout was monitored for a final minute with a final external solution addition.

Peak outward current magnitude at 540 ms through 550 ms of the IV step was measured for each sweep at 10 s intervals. The final measurement at +100 mV is used for calculating the effect of each compound application. Each cell’s current magnitude was then normalized to itself at the initial 0.3% vehicle control period prior to compound application. This step is to account for differences in current size for each cell. Normalized responses to test article were plotted against their concentrations to reveal a concentration-inhibition plot.

A Normalized Group Hill fit was then performed on the plotted results to yield a pooled IC50 value, reconstructed from a minimum of two cells/concentration. The baseline response as well as the full response were constrained to 1 and 0, respectively.

IonWorks Barracuda Electrophysiology Studies

HEK293 cells stably expressing the human TMEM16A (abc) variant and CHO cells expressing the human CFTR gene were from ChanTest. The IonWorks Barracuda (IWB) procedure and data provided are from outsourced studies performed at ChanTest (Charles River Discovery, Cleveland, OH, United States) who were blinded as to the identity of the test articles or compounds provided by Amgen. As a control, benchmark inhibitors CFTRinh-172, GlyH-101 and benzbromarone were included amongst the blinded compounds submitted and performed as expected.

For IonWorks Barracuda studies on TMEM16A compounds, in brief eight test article concentrations were applied to naïve cells (n = 4, where n = the number of replicate wells/concentration) via steel needles of a 384-channel pipettor. Each application consists of addition of 20 μl of 2X concentrated test article solution to the total 40 μl final volume of the extracellular well of the Population Patch ClampTM (PPC) planar electrode. This addition is followed by mixing (one time) of the PPC well content. Duration of exposure to each test article concentration was at least 5 min. The electrophysiology procedure used: (a) Intracellular solution containing 50 mM CsCl, 90 mM CsF, 5 mM MgCl2, 1 mM EGTA, 10 mM HEPES, adjusted to pH 7.2 with CsOH; (b) Amphotericin B for patch perforation, where 30 mg/ml stock solution of amphotericin B in DMSO is added to internal solution to final concentration of 33.3 μg/ml; (c) Extracellular solution containing HEPES-buffered physiological saline (HBPS): 137 mM NaCl, 4 mM KCl, 1.8 mM CaCl2, 1 mM MgCl2, 10 mM HEPES, 10 mM glucose, adjusted to pH 7.4 with NaOH; and (d) Ionomycin stimulation of chloride currents where 10 μM ionomycin is added to all test solutions including vehicle and positive controls. The current was elicited by a 500-ms step pulse to 0 mV followed 1000-ms step pulse to −100 mV from holding potential, −30 mV, with stimulation frequency 0.05 Hz. The specific recording procedure was as follows: extracellular buffer is loaded into the PPC plate wells (11 μl per well). Cell suspension is then pipetted into the wells (9 μl per well) of the PPC planar electrode. After establishment of a whole-cell configuration via patch perforation (∼7–10 min exposure to amphotericin B), membrane currents were recorded using the on-board patch clamp amplifiers. Recordings (scans) were performed as follows: three scans before and fifteen scans during the 5-min interval after ionomycin and test article application. A full dose-response of benzbromarone was included on each plate as a positive control, while multiple replicates of DMSO were included as negative control. Final DMSO concentration for test and control articles was 0.3%.

For measuring compound effects on CFTR chloride currents, compounds were serially diluted in HEPES-buffered physiological saline to 2X final concentration allowing for an 8-point dose-response analysis. Test article concentrations were applied to naïve cells (n = 4, where n = the number of replicate wells/concentration) via steel needles, where each application will consist of addition of 20 μl of 2X concentrated test article solution to a final 40 μl volume in the extracellular well of the Population Patch ClampTM (PPC) planar electrode. After mixing (three times), duration of exposure to compound is at least 5 min. Final solutions contain 0.3% DMSO. The electrophysiology procedure used: (a) Intracellular solution (mM): CsCl, 50; CsF 90; MgCl2, 5; EGTA, 1; HEPES, 10; adjusted to pH 7.2 with KOH, (b) Extracellular, HB PS solution (composition in mM): NaCl, 137.0; KCl, 4.0; CaCl2, 1.8; MgCl2, 1; HEPES, 10; adjusted to pH 7.4 with NaOH; and (c) Stimulation, where CFTR current is activated with 20 μM forskolin added to all test solutions including vehicle and positive controls. The current is measured using a pulse pattern consisting of a voltage step to +60 mV, 100 ms duration; voltage ramp from +60 to −120 mV, 1000 ms; step to −120 mV, 200 ms; and step to 0 mV, 200 ms; from holding potential, −30 mV. Current amplitudes were measured at the voltage step to 0 mV. The recording procedure was as follows: extracellular buffer is loaded into the PPC plate wells (11 μl per well). Cell suspension was pipetted into the wells (9 μl per well) of the PPC planar electrode. After establishment of a whole-cell configuration via amphotericin patch perforation similar to that described earlier, membrane currents were recorded using the patch clamp amplifier in the IonWorksTM Barracuda system. Two recordings (scans) were then performed, one scan before and a second scan 5 min after test article (TA) application. Every plate included a full dose-response of CFTRinh-172 as a positive control. Data acquisition and analyses were performed using the IonWorks BarracudaTM system operation software (version 2.0.2). The decrease in current amplitude after test article application was used to calculate the percent block relative to control. Results for each test article concentration (n ≥ 2) were averaged; the mean and standard deviation values were calculated and used to generate dose-response curves. Inhibition effect was calculated as: % Inhibition = (1 – ITA/IBaseline) × 100%, where IBaseline and ITA were the currents measured in control (before addition of a test article) and in the presence of a test article, respectively.

Histamine-Induced Depolarization of Human Bronchial Smooth Muscle Cells

Cultured human bronchial smooth muscle cells (BSMC) from Lonza were resuspended in Smooth Muscle Growth Medium-2 (SmGM-2; Lonza) at a concentration of 4 × 105 cells per ml. One hundred microliters of cells per well were plated in a black wall clear bottom polystyrene 96-well tissue culture plate and incubated in a 37°C humidified incubator with 5% CO2 overnight. Cells were serum starved by removing the SmGM-2 and replaced with one hundred microliters per well of Smooth Muscle Cell Basal Medium 2 (SmBM-2), phenol red-free (PromoCell), and incubated in a 37°C humidified incubator with 5% CO2 for 24 h. SmBM-2 was replaced with 100 microliters of fresh SmBM-2. Compounds were dissolved in 100% DMSO and serially diluted ½ log in a polypropylene 96-well microtiter plate (drug plate). Columns 6 and 12 were reserved as controls (HI control and LO control respectively) and contained only DMSO. Serially diluted compounds were diluted in SmBM-2 to 10X the final concentration. Twenty-five microliters of 10X compound titrations, and one hundred microliters of Blue dye-loading buffer (FLIPR Membrane Potential Assay Kit; Molecular Devices) were added to the cells. Cells were pre-incubated at room temperature with compound for 0.5 h. Five micromolar (10X) histamine was prepared in SmBM-2. Using the FLIPR-Tetra (Molecular Devices), three measurements of the baseline fluorescence were taken over the span of 1 min. Twenty-five microliters per well of 10X histamine were added to the first 11 columns of the plate containing the compound treated cells. Twenty-five microliters of SmBM-2 were added to column 12 for the LO control. Fifty-five fluorescence measurements were taken over the span of 14 min. The area under the curve (AUC) from the fluorescence kinetic traces (normalized to baseline fluorescence) were calculated. The amount of fluorescence in the presence of compound compared with that in the presence of DMSO vehicle alone (HI control) was calculated using the formula: % control (POC) = (compd – average LO)/(average HI – average LO)100.

While the method above using histamine represents the standard protocol we used to evaluate compound effects in blocking pro-contractile depolarization of human BSMCs, for some experiments provided as Supplementary Figures 144 we employed the cholinergics methacholine or carbachol instead of histamine, or evaluated if the TMEM16A opener, Eact, could itself act like contractants to induce BSMC membrane depolarization.

Measuring Bronchodilation of Mouse Tracheal Rings by Wire Myograph

Instrument and Reagents

A series of three Danish Myograph Technologies (DMT) 620M Multi Wire Myograph Systems instruments were typically used for each experiment, with each instrument containing four chambers and allowing tests on four tracheal rings. The instruments were interfaced with PowerLab 4/35 or PowerLab 8/35 data acquisition systems and a computer running DMT Device Enabler and LabChart Pro v8 (AD Instruments) for automatic recognition of devices and simultaneous recording of data. Carbachol (Cat# C4382) and isoproterenol (Cat# I2760) were from Sigma. Mouse IL-13 and mouse IL-1β were from R&D Systems. Histamine, theophylline and mouse TNFα were from Amgen.

Mice

Mice used for wire myograph studies were housed in groups at an AAALAC, International accredited facility. Animals were cared for in accordance with the Guide for the Care and Use of Laboratory Animals, 8th Edition. All research protocols were reviewed and approved by the Amgen Institutional Animal Care and Use Committee. Female C57 BL/6 (CRL, >12 weeks of age) were housed in individual ventilated caging (IVC) system on an irradiated corncob bedding (Envigo Teklad 7097). Lighting in animal holding rooms was maintained on 12:12 h light:dark cycle, and the ambient temperature and humidity range was at 68–79°F and 30–70%, respectively. Animals had ad libitum access to irradiated pelleted feed2 (Envigo Teklad Global Rodent Diet- soy protein free extruded 2020X) and reverse-osmosis (RO) chlorinated (0.3–0.5 ppm) water via an automatic watering system. Cages were changed biweekly inside an engineered cage changing station.

Standard Wire Myograph Studies for Measuring Bronchodilation

The trachea from C57 BL/6 mice was dissected and collected in PBS with Mg2+ and Ca2+. After further trimming, two 2 mm sections (rings) per trachea were then placed in DMEM media containing PSG/HEPES/AA/NaPyr for 1 h and used immediately or maintained overnight at 37°C in a humidified incubator with 5% CO2. Rings were then mounted into chambers of DMT wire myograph containing L-shaped mounting pins with 2–3 mN of force applied to secure rings on the pins and the rings were allowed to equilibrate for 25–30 min in a physiological saline buffer (PSS, 130 mM NaCl, 4.7 mM KCl, 1.2 mM MgS04, 14.9 mM NaHCO3, 1.2 mM KH2PO4, 0.026 mM EDTA, 1.6 mM CaCl2, 5.5 mM Dextrose). Throughout the experiment, rings were maintained at physiological temperature and gas conditions by heating buffer reservoir and chamber to 37°C and bubbling a 95% O2, 5% CO2 gas mixture into chamber.

After equilibration, a standard tissue “Wake-Up” procedure was applied by altering tension and then treatment with elevated potassium as follows: 3 mN tension was applied for 5–10 min, then another 2–3 mN tension to total of 5.5 mN and then to 7–8 mN for another 5 min; the tissue pre-tension was then set to 5 mN and allowed to equilibrate for 5–10 min followed by the removal of the PSS buffer and addition of KPSS (60 mM KCl, 74.7 mM NaCl in PSS solution), with force changes monitored until they reached a plateau. The tissue was then washed four times with PSS over 5 min and the KPSS treatment and washes repeated; the tissue was then allowed to sit in PSS for 10 min and the tension was set to 5 mN.

To qualify every tissue ring and determine its sensitivity to contractant and reference bronchodilator, 6 ml of PSS was added to each chamber and after 5 min ascending doses of the carbachol (CCh) contractant were added to determine the EC25, EC50, EC75, and EC95 for CCh for each ring (see Supplementary Figure 33 for example of raw traces and dose-response), where sufficient time (at least 5 min) was allowed for response to plateau after each dose. At the end of the CCh dose-response study when no further increases in force were observed, the β-agonist and reference bronchodilator isoproterenol was added to 10 μM final concentration to monitor tissue relaxation. The tissue was then washed three times and then three times more with PSS for 5 min until tension returned to baseline.

Only tissue that passed the qualification tests above advanced to studies on test article compounds. Fresh PSS (6 ml) was then added to each tissue bath, tension was reset to 5 mN and the rings were treated with the calculated EC25, EC50, EC75, or EC95 concentration of CCh from GraphPad Prism analysis and allowed to incubate for at least 25 min and a plateau was achieved. For standard ascending dose-response studies on TMEM16A antagonists, an EC75 or EC95 concentration of CCh was used to pre-contract tracheal rings. Ascending doses of the isoproterenol positive control, the vehicle negative controls or the test compound were then added to appropriate rings and incubated for 10 min or until plateau was achieved after each addition. At the end of the experiment when no further changes were observed to test article or positive control, theophylline was added to 2.2 mM final to fully relax airways. The dose-response curved and EC50 values for bronchodilation were derived from GraphPad Prism analysis of the data normalized to theophylline as 100% relaxation.

Bronchodilation of Tracheal Rings as Function of Differing Levels of Pre-contraction

Standard procedures as listed above were used to prepare, mount and Wake-Up tracheal rings and the concentrations of carbachol providing EC25, EC50, EC75, or EC95 levels of contraction for each ring was determined as described above. Rings with differing levels of CCh pre-contraction were then treated with ascending doses of isoproterenol, benzbromarone, niclosamide or Compound 4 to determine the efficacy of bronchodilation, which was normalized to theophylline added at the end as control of 100% relaxation.

Effects of Cytokines on Compound Efficacy in Relaxing Mouse Tracheal Rings

Airway rings were treated overnight in 2 ml DMEM + PSG/HEPES/AA/NaPyr buffer alone or with the same solution supplemented with 100 ng/ml of mTNFα, mIL-13 and mIL-1β. Unlike like the typical Wake-Up protocol and ring qualifying tests, mouse tracheal rings were not exposed to KPSS and rings were not observed to relax with a β-agonist following discovery of CCh dose-responses. After several washes with PSS, EC75 concentrations of CCh were used to pre-contract rings for >30 min alone or in the presence of cytokines. Ascending doses of isoproterenol, benzbromarone, or niclosamide were given to determine efficacy of bronchodilation normalized to theophylline as control of 100% relaxation.

Evaluating β-Agonist Short-Term Use-Dependent Desensitization

Standard procedures were used to prepare and mount tissues. However, similar to the experiments assessing the effects of cytokines in relaxing tracheal rings, no KPSS or qualifying tests with isoproterenol to observe relaxation following discovery of CCh dose-responses. To assess use-dependent desensitization, EC50 concentrations of CCh were used to pre-contract rings followed by a dose-response with isoproterenol. Rings were then washed several times with PSS, pre-contracted again with EC50 CCh, and given another round of isoproterenol dose-response normalized to theophylline as control of 100% relaxation.

Assessing Compound Duration of Action in Bronchodilating

Mouse airway rings were exposed to a typical Wake-Up procedure including the use of KPSS and discovery of EC75 CCh dose responses. However, no isoproterenol was used to observe tissue relaxation as part of the qualification procedure. With pre-contraction with EC75 CCh for >30 min, rings were subjected to a compound dose-response, followed by a complete washout of compound, and re-contracted with EC75 CCh. Traces were recorded for several hours thereafter and normalized to theophylline as control of 100% relaxation.

Measuring Bronchodilation of Human Bronchial Rings by Wire Myograph

Outsourced Bronchodilator Studies Performed at Biopta (Glasgow, United Kingdom)

Details of bronchodilator studies shown in Figures 5H, 6 and Supplementary Figure 40 are provided below as part of the end assay report on Study CUR011 by Lee Christie (Biopta). Donors with any of the following conditions were excluded from the study: asthma, COPD, emphysema, lung cancer, cystic fibrosis, pulmonary fibrosis, pulmonary hypertension, pneumonia. Donors were also excluded if they had smoked in the past 12 months. Non-transplantable lung tissue was acquired from organ donors post mortem. Any macroscopically diseased/necrotic tissue was rejected. Furthermore, any tissues that did not respond to functional checks were rejected.

Human airway rings were set up under isometric conditions on a wire myograph in order to examine the influence of the test articles on bronchodilation. In order to assess tissue viability, the airways were challenged with carbachol (10 μM) and then isoprenaline (10 μM) to assess their constriction and relaxation responses, respectively. Airways that did not respond to these initial checks were not used.

Viable airways were pre-constricted with histamine (10 μM), then exposed to one of the following cumulative concentration response curves (CCRCs): positive control, vehicle (water, volume matched to positive control), test article, vehicle (DMSO concentration matched to test articles). At the end of each CCRC, theophylline (2.2 mM) was added to induce maximal relaxation of the tissue.

Specific methodology was as follows: (a) Quaternary branches of human airway rings were dissected free from surrounding parenchyma, cut into 2 mm rings and mounted in 5 ml organ baths containing physiological saline solution (composition: 119.0 mM NaCl, 4.7 mM KCl, 1.2 mM MgS04, 24.9 mM NaHCO3, 1.2 mM KH2PO4, 2.5 mM CaCl2, 11.1 mM glucose, and 5 μM indomethacin), aerated with 95% O2/5% CO2, and warmed to 37°C. The tissues were allowed to equilibrate for approximately 30 min, with washes approximately every 10 min, before being set to a tension of approximately 1.0 – 1.5 g and then allowed to equilibrate for a further 90 min with washes approximately every 15 min. Airways were re-tensioned to 1.0 – 1.5 g if the tension had dropped below 1.0 g during the first 30 min of equilibration. (b) Airways were exposed to carbachol (bath concentration 10 μM) in order to measure their contractile responses and then isoprenaline (10 μM) in order to assess their relaxation. Airways were washed out and allowed to return to baseline. (c) Airways were pre-constricted with histamine (10 μM) before conducting a CCRC to the test article, the positive control (isoprenaline), test article vehicle or positive control vehicle. It should be noted there were no significant effect by any of the vehicles in the studies reported.

Bronchodilator Studies Performed Internally at Amgen (Thousand Oaks, CA, United States)

Bronchodilator studies using human bronchial rings (Figures 5E–G and Supplementary Figure 38) used the same instruments and reagents described earlier for studies on mouse tracheal rings.

Non-transplantable human lungs were obtained through the International Institute for the Advancement of Medicine, IIAM (Edison, NJ, United States) from non-smoking donors who had been ventilated for <3 days with acceptable blood gasses and used within 24 h of cross-clamp time.

Human 4th order bronchial rings were isolated, sliced into 2 mm long rings and then mounted individually in chambers of wire myographs containing 6 ml room temperature PSS per chamber. Chambers are aerated with 95% O2/5% CO2 throughout the experiment. Bronchial rings are allowed to equilibrate in PSS while chambers are warmed to 37°C. Tension on the airway rings was gradually increased until reaching a steady state passive tension of 9.8 mN. The rings were then equilibrated for additional 40–60 min, changing buffer every 15–20 min. Tension was adjusted if it dropped below 9.8 mN. Airway rings were “woken up” by exposing to pre-warmed, aerated KPSS and allowing rings to reach plateau of constriction, followed by washing rings four times with PSS. This procedure was repeated two additional times.

Airway contraction response was then assessed by treating with increasing doses of the contractant, carbachol (10 nM to 10 μM), waiting for response plateaus (at least 5 min between additions). Airways are then relaxed by addition of 10 μM isoproterenol. Airways that demonstrated expected contraction and relaxation were used to evaluate test compounds after a washout and re-equilibration period.

Airway rings were contracted with EC75 of carbachol and then exposed to one of the following CCRCs: vehicle (DMSO concentration matched to test article), isoproterenol (positive control), test compound. At the end of each CCRC, rings were fully relaxed by addition of 2.2 mM theophylline. One hundred percent contraction was calculated by subtracting the stable tension remaining in the airway ring after theophylline addition from the tension achieved after addition of EC75 carbachol while 100% relaxation is the net tension remaining after theophylline addition. Dose response curves and EC50 values for bronchodilation were determined in GraphPad Prism.

Evaluating the TMEM16A Expression in Airway Smooth Muscle Cells and the Bronchial Epithelium

RNA Sequencing

The RNA-Seq Amgen lung cell dataset containing primary human bronchial epithelial cells and airway smooth muscle cells has been described earlier, as has the method for RNA sequencing (Aisenberg et al., 2016), with modifications and additional details provided below. Similar RNA sequencing and data analysis methods were applied to generate the new RNA-Seq dataset described here for untreated or IL-13 treated mature bronchial epithelial ALI cultures. Information on patient demographics for each donor sample are provided in Supplementary Tables 24.

Data were analyzed using the Array Studio (Omicsoft, Cary, NC, United States) platform, as previously described (Aisenberg et al., 2016), but using v9.0 of the Oshell software and the GENCODE (Harrow et al., 2012) gene model (release 24; Comprehensive version was used for alignment, and Basic version was used for quantification). FPKM values were normalized using a modified version of upper-quartile normalization (Mortazavi et al., 2008; Robinson and Oshlack, 2010) in which the 70th percentile FPKM among genes was fixed at a value of 10 (excluding those genes belonging to families with high homology or with maximal transcript length < 500 bp). Data from the Cancer Cell Line Encyclopedia (Barretina et al., 2012) were processed by OmicSoft using the same analysis pipeline and normalized using the same methods.

Generating Human Bronchial Epithelial ALI Cultures and Effects of Th2 Cytokines

To determine the effects of IL-13 on TMEM16A versus Muc5AC mRNA expression over time in the human bronchial epithelium and evaluate TMEM16A alternative splicing, normal and COPD human bronchial epithelial cells from Lonza were grown on 6.5 mm permeable supports (Corning Transwell 3470) submerged in apical and basolateral growth media (BEGM, Lonza, CC-3170), until confluent at ∼5 days. Once confluent, apical media was removed, and basolateral growth media was replaced with ALI maintenance media (PneumaCult, STEMCELL Technologies #05001) to initiate the air-liquid interface (ALI). Cultures were then maintained at ALI for ∼21 days, with basolateral maintenance media being replenished on Monday–Wednesday–Friday schedule, to generate fully differentiated human bronchial epithelial cells. Cultures were then left untreated or treated with 20 ng/ml IL-13 for 1, 3, 5, or 7 days. For each time point, RNA was prepared from the ALI cultures for the NextGen RNA sequencing following the Qiagen RNeasy Mini kit (cat # 74104) protocol. Each ALI culture was solubilized with 350 μl of Qiagen lysis buffer RLT on the apical transwell and immediately spun through a Qiashredder (Qiagen # 79654) at 15,000 rpm for 2 min. The samples were either stored at −80°C and processed at a later date or immediately prepared following the Qiagen RNeasy Mini kit protocol including Part 2 which contains the Qiagen RNase-Free DNase set (Cat # 79254) treatment. Total RNA was eluted in 45 μl of RNase free water supplied from the kit and quantified on a Nanodrop ND-1000 spectrophotometer.

TMEM16A protein expression was determined by immunohistochemistry of untreated or IL-4 or IL-13 treated human bronchial epithelial ALI cultures. Fully differentiated human bronchial epithelial ALI cultures were generated using methods similar to above, except cultures were expanded on permeable supports for 3 days and maintained at ALI for 31 days prior to treatment. At 31 days post-airlift, cultures were left untreated or treated with either 20 ng/ml IL-4 or 20 ng/ml IL-13 for 1, 2, 3, or 5 days. For each timepoint, Transwells were fixed both apically and basolaterally with 2% paraformaldehyde for 2 h at room temperature, then washed with 70% ethanol. Permeable membranes containing the ALI cultures were then cut away from the plastic insert using a scalpel and processed for immunohistochemistry as described below.

TMEM16A Expression in Tissue Isolated From Naïve or Asthmatic Cynomolgus Monkeys

Tissue from naïve or asthmatic cynomolgus monkeys was obtained from Charles River Laboratories. Naïve adult cynos were never challenged with Ascaris suum aerosol, were from the standard CRL colony and were negative by intradermal screening for A. suum. Adult cynos exhibiting sensitivity to A. suum antigen and characterized as asthmatic, presumably due to early exposure and allergy to A. suum or similar parasite, were maintained in a separate colony and confirmed over time as reproducible sensitivity to an A. suum aerosol challenge. Naïve, unchallenged asthmatic (>2 weeks), asthmatic acute challenged (4 h post aerosol A. suum) and asthmatic subacute challenged (24 h post A. suum aerosol) cynos were euthanized and lung vasculature infused with cold HypoThermosol biopreservation media. Lungs and trachea were shipped on ice packs overnight and necropsy was immediately performed on arrival to prepare specimens for immunohistochemistry and bronchodilation studies.

Immunohistochemistry

The routinely formalin-fixed, paraffin embedded tissue blocks were sectioned at a 4 μm thickness and processed for IHC. Paraffin was removed from the tissue sections with xylene and the sections were rehydrated with graded ethanol and immersed in distilled water. Antigen retrieval was performed using Diva Decloaker pretreatment reagent (Biocare Medical, Concord, CA, United States) in a Biocare Decloaking Chamber (Biocare Medical, Concord, CA, United States) set to reach 125°C for 30 s, then 90°C for 10 s. Tissue sections were processed at room temperature in a Lab Vision 720 automated staining instrument (Thermo Scientific, Waltham, MA, United States). Endogenous peroxidase was blocked using Hydrogen Peroxide Block (cat# TA-125-HP, Thermo Scientific, Waltham, MA, United States) for 10 min. Protein Block (cat# X0909, Dako North America) was applied to the tissue sections for 10 min. Tissue sections were incubated with a rabbit anti-human TMEM16A polyclonal antibody (cat# Ab53212; 2 μg/ml) or rabbit anti-human TMEM16A monoclonal antibody (SP31, cat# Ab64085; 1:50) from Abcam, San Francisco, CA, United States. Envision+ System HRP labeled Polymer (cat# K4003, Dako North America) was used to detect the primary antibody. Sections were incubated with DAB plus chromogen substrate (cat#K3468) for 5 min.

The specificity of the two anti-TMEM16A antibodies was confirmed by western blot and immunohistochemical (IHC) or immunofluorescence (IF) studies that showed selective staining of TMEM16A transfected versus untransfected HEK293 cells. The performance and specificity of the antibodies was also validated by IHC studies against a panel of cynomolgus monkey and human tissues where they showed specific staining of salivary gland acinar cells as a positive control tissue (Romanenko et al., 2010; Chenevert et al., 2012). In our hands, both antibodies performed similarly in IHC studies and selectively stained goblet cells and submucosal glands in the lung. The utility of these two antibodies for monitoring TMEM16A expression has also been described in several earlier reports (Davis et al., 2010; Scudieri et al., 2012; Britschgi et al., 2013; Gallos et al., 2013; Caci et al., 2015).

Curve Fitting, Statistical Analysis, and Small Molecule Informatics

Data were analyzed using GraphPad Prism v7.0 (GraphPad Software, San Diego, CA, United States). All data is expressed as mean ± SEM, unless otherwise indicated. Concentration-response curves were fit by non-linear regression. Data were compared using the Student’s t-test for unpaired samples or by two-way ANOVA using Dunnett’s or Sidak’s multiple comparison tests, as appropriate. P-values < 0.05 were considered statistically significant.

Compound lipophilicity was determined using the SwissADME free web tool (Daina et al., 2017) and the XLogP3 additive model for logP calculation (Cheng et al., 2007).

Results

High Throughput Screening of Small Molecule Library Identifies Niclosamide as a Novel TMEM16A Antagonist

To identify novel inhibitors of the CaCC TMEM16A, we generated a stable cell line HEK293T:eYFP:TMEM16A(abc) co-expressing the ‘abc’ splice variant of human TMEM16A along with the halide-sensitive YFP mutant [H148Q,I152L] referred to as eYFP. To this end, HEK293T:eYFP cells, which show no significant response to ionomycin and iodide, were stably transfected with a TMEM16A(abc) expression cassette. Clones suitable for high-throughput screening were identified by functional clone selection using the eYFP assay (data not shown). For the HTS eYFP assay, cells were incubated with the test compounds for 30 min and then treated with 10 mM iodide and 2 μM ionomycin. Ionomycin triggers the Ca2+-dependent activation of TMEM16A, which allows iodide ions to enter the cell and quench eYFP fluorescence. eYFP assay results were normalized based on the benchmark TMEM16A inhibitor benzbromarone providing 100 percent of inhibition (POI) of TMEM16A activity. This corresponds to 0 percent of control (POC) response of Ca2+-dependent TMEM16A activity.

A library of ∼580,000 compounds was screened at a 5 μM final concentration in the TMEM16A halide-sensitive YFP assay. A summary of the screen and the hit triaging process is provided in Supplementary Table 1 and Supplementary Figure 1, and described further below. A total of 1,445 primary hits were identified using a hit cut-off of POC < 72.6%, which reflects compounds providing > 27.4% inhibition. All primary hits were tested again in three independent, consecutive runs and confirmed 673 hits. In order to eliminate false-positive TMEM16A inhibitors, auto-fluorescent compounds were removed revealing 491 non-fluorescent hits. A Ca2+ flux counterscreen was run to eliminate compounds non-specifically interfering with ionomycin-triggered calcium signaling further reducing the number of hits to 328 compounds. Dose response testing of these hits in the TMEM16A eYFP assay yielded 145 hits with an IC50 < 5 μM and finally 130 hits that passed quality control by mass spec analysis. Figure 1A shows the distribution of hit potency versus average POC response in the HTS assay. Niclosamide (Compound 1) and a related analog, Compound 4, were identified for the first time from our screen as potent TMEM16A inhibitors. Niclosamide with an IC50 of 140 nM in the eYFP HTS assay (Figure 1A) represents one of the most potent TMEM16A antagonists described to date. Compound 4, a structurally related analog (Figure 2), showed lesser activity with an IC50 of 970 nM, but both this compound and niclosamide showed only partial block of the halide-sensitive YFP response with a POC of 61.1 ± 7.8 and 69.0 ± 3.3, respectively, reflecting just 38.9 and 31.0% inhibition. In fact, many of the hits from the screen provided less than 50% inhibition of the iodide eYFP quenching after ionomycin-activation of TMEM16A (Figure 1A). This contrasts with the benchmark antagonist benzbromarone, which fully blocked the halide-sensitive response and thus served as a positive control for 100% inhibition (POC = 0). Interestingly, the TMEM16A inhibitor 1PBC was also included in our small molecule library and showed an activity of 74.8 POC (Figure 1A) corresponding to 25.2% inhibition. To characterize this compound further we performed 22-point dose-response analysis and found that 1PBC inhibited TMEM16A with IC50 of 1.05 μM but only partially blocked the halide-sensitive response just like niclosamide.

FIGURE 1
www.frontiersin.org

Figure 1. High-throughput screen using a halide-sensitive YFP assay identifies niclosamide as a highly potent, nanomolar inhibitor of TMEM16A. (A) Potency range of the 145 hits from the TMEM16A high-throughput screen and triplicate hit confirmation at 5 μM (n of 4) that reduced the YFP quenching by less than 72.6 percent of control (POC) response to ionomycin (Iono)/iodide alone. The dotted horizontal and vertical lines demarcate the cut-off values (POC < 72.6, IC50 < 5 μM) for validated hits. Representative data from the medchem TMEM16A YFP assay showing FLIPR kinetic traces for the DMSO vehicle control or compounds at 55.6 μM (B), and dose-response curves indicating niclosamide (C), Cpd 4 (D) and 1PBC (E) only partially blocked the ionomycin-induced YFP response, while benzbromarone (F) and CaCCinh-A01 (H) gave full block. Niclosamide and related compounds, Cpd 2 – Cpd 5, were more potent TMEM16A antagonists than benchmarks (G,H), but only partially blocked the iodide/YFP response (I). Nitazoxanide and tizoxanide were identified as additional inhibitors of TMEM16A (H). Mean ± SEM, n of 3–58. ∗∗∗∗P < 0.0001; significant difference from the benzbromarone benchmark antagonist; ns, not significant (unpaired t-test).

FIGURE 2
www.frontiersin.org

Figure 2. Structures of small molecule antagonists and their activity in the TMEM16A halide-sensitive YFP assay. Molecular structure of the hits niclosamide and Cpd 4 identified from our high-throughput screen and related compounds discovered during our hit-to-lead process, which includes the approved drug nitazoxanide and its metabolic product tizoxanide (A). The structures of TMEM16A benchmark antagonists described earlier are shown in panel (B) on the right. The average potency of each compound in inhibiting the TMEM16A halide-sensitive YFP response was determined by 22-point dose-response studies, with the average IC50 (μM) listed beneath each compound (n = 3–58; with exception of compounds labeled with which were single determinations).

Side by side comparison of TMEM16A antagonists using HEK293T:eYFP:TMEM16A(abc) cells in medchem follow-up assays revealed that treatment with 1PBC, Cpd 4 and niclosamide caused only partial inhibition of TMEM16A activity whereas benzbromarone blocks almost completely (Figure 1B). This was further corroborated in 22-point dose-response analyses showing that the benchmark antagonist benzbromarone fully inhibited the iodide eYFP quenching after ionomycin-activation of TMEM16A (Figure 1F), while niclosamide, Cpd 4 and 1PBC showed only partial inhibition even at the highest concentration of 55.6 μM (Figures 1C–E). The scope of their block in this follow-up characterization was more pronounced compared to inhibition during HTS. For instance, niclosamide provided about 31% inhibition during HTS but achieved up to 67% inhibition in the dose-response analysis (Figure 1C). However, its potency was unchanged. Remarkably, niclosamide with an average IC50 of 132 ± 67 nM (n = 42) was strikingly more potent than the benchmark antagonist benzbromarone having an average IC50 of 1.508 ± 0.504 (n = 58) μM (Figure 1G and Supplementary Figure 2). Compound 4 revealed an average IC50 of 0.901 ± 0.524 μM (n = 53) in comparison to 1PBC with an average IC50 of 1.523 ± 1.360 μM (n = 17) and CaCCinh-A01 with an average IC50 of 8.125 ± 2.720 μM (n = 17).

Activity of Additional Niclosamide-Related Compounds and Identification of the Approved Drug Nitazoxanide as a New TMEM16A Antagonist

The structure and activity of additional niclosamide-related compounds are provided in Figure 2A. Compounds 2, 3, and 5 provided sub-μM block of the TMEM16A eYFP response and improved activity compared to benzbromarone (Figure 1G), but like niclosamide only partially inhibited the ionomycin-induced eYFP quenching by iodide (Figures 1H,I). Besides benzbromarone, we also evaluated the TMEM16A antagonists dichlorophen, CaCCinh-A01, MONNA, niflumic acid and T16Ainh-A01. The TMEM16A antagonist CaCCinh-A01 fully inhibited the ionomycin-induced iodide eYFP quenching (Figures 1H,I) like benzbromarone, as did the benchmark antagonists MONNA and dichlorophen which had an average IC50 of 6.295 ± 1.470 (n = 19) and 4.365 ± 1.210 μM (n = 58), respectively, in the halide-sensitive YFP assay (Figure 2 and Supplementary Figure 3). Niflumic acid showed lesser activity with an average IC50 of 26.570 ± 11.800 μM (n = 55) while the compound T16Ainh-A01 was inactive despite repeat tests, IC50 > 50 μM (n = 16) (Supplementary Figure 3).

The niclosamide-related Compounds 2–5 represent only a small subset of compounds we’ve tested in this series to develop a structure-activity relationship (Supplementary Figure 4). Compounds 2 and 5 are instructive in showing the nitro group is unnecessary for TMEM16A bioactivity and can be replaced with a chlorine atom (Figure 2A). The World Health Organization includes niclosamide on its list of Essential Medicines based on its efficacy, safety and cost-effectiveness. Another approved drug, nitazoxanide, is synthesized using the scaffold of niclosamide. Since nitazoxanide and its metabolic product tizoxanide appear structurally similar to niclosamide (Figure 2A), we tested these drugs for activity in blocking TMEM16A. Both nitazoxanide and tizoxanide were found to be antagonists of TMEM16A (Figure 1H) and exhibit a pharmacology similar to niclosamide in partially inhibiting the iodide/eYFP response.

The TMEM16A inhibitor, NTTP, is noteworthy as it shares the mechanism of 1PBC in blocking TMEM16A by binding four basic residues in the TMEM16A selectivity filter. We have confirmed NTTP is a TMEM16A antagonist (Supplementary Figure 3). Interestingly, an examination of the chemical structure of NTTP reveals it is a highly similar analog of tizoxanide (Figure 2). This suggests nitazoxanide and its metabolic product tizoxanide may share the mechanism of NTTP in binding the pore region of TMEM16A to block channel conductivity.

In summary, the direct comparison of niclosamide with eight distinct TMEM16A benchmark antagonists revealed that niclosamide with an IC50 of 132 nM in the halide-sensitive YFP assay is 10–200 times more potent than other antagonists. Our results further suggest niclosamide and related compounds may utilize a unique mechanism for channel block. Therefore, we performed additional experiments to characterize their impact on chloride currents and effects in blocking airway smooth muscle cell depolarization and contraction of airways.

Analysis of Niclosamide, Nitazoxanide and Related Compounds for Inhibition of TMEM16A Calcium-Activated Chloride Currents

While our screen identified potent TMEM16A antagonists, we wondered why many just partially inhibited the ionomycin-induced quenching of the eYFP halide sensor by iodide in the YFP assay (Figures 1H,I). Therefore, we asked how these TMEM16A antagonists perform when measuring their efficacy in blocking calcium-activated chloride currents. Electrophysiology studies on this channel, however, can be notoriously difficult as TMEM16A exhibits rapid calcium-dependent inactivation (Wang and Kotlikoff, 1997; Tian et al., 2012). An additional challenge is that the ephys assay must have enough throughput to test numerous compounds and quickly return data to the chemist to guide compound synthesis and the hit-to-lead (HTL) process. Because no such assay was yet available, we developed an automated QPatch electrophysiology assay that enables prolonged recordings and potency determinations on up to 100 compounds per week. Whole-cell currents from HEK293 cells stably transfected with TMEM16A (acd) were recorded using 170 nM free intracellular calcium and 20 mV steps from −100 to +100 mV from a holding potential of 0 mV (Figure 3A). Three voltage protocols were run per sample addition allowing at least 60 s incubation time following each vehicle or compound addition. The currents were voltage-dependent and outward rectifying as expected for TMEM16A recordings in low intracellular calcium. Importantly, the currents were stable for about 10 min and insensitive to the DMSO vehicle control as shown in the representative QPatch instrument recording in Figure 3B where the dots reflect the measured current after each 20 mV step from −100 to +100 mV and three IV protocols are shown for each sample addition. While the optimized assay typically showed < 20% rundown over this time frame, it should be noted significant rundown becomes an issue with recordings that are longer than 12–15 min. The benchmark TMEM16A antagonist benzbromarone fully inhibited both the outward and inward currents (Figures 3C,D). Typical QPatch recording for dose-response analysis included two saline and DMSO additions to validate current stability and voltage-dependence, followed by five concentrations of antagonists as shown in Figure 3C. The sustained TMEM16A current following the second DMSO addition and third IV protocol served as the baseline current for calculating percent inhibition. Routine calculations of compound potency measured their effect in reducing the outward current at +100 mV.

FIGURE 3
www.frontiersin.org

Figure 3. A medium-throughput, optimized planar patch clamp electrophysiology assay for stable recordings and routine measures of compound activity in blocking TMEM16A. (A) Representative whole-cell currents of a HEK293 cell stably transfected with TMEM16A (acd) using an intracellular solution containing 170 nM free calcium. The standard QPatch voltage protocol had cells clamped from a holding potential of 0 mV (50 ms) to voltages between –100 and +100 mV in 20 mV steps (500 ms) followed by a step to –100 mV for 50 ms. There were 1.2 s between sweeps, 20 s between protocols, and three protocols run per sample addition which allowed at least 60 s incubation time following each vehicle or compound addition. (B) The TMEM16A currents were stable for at least 10 min and after repeat additions of the vehicle control (0.3% DMSO) in the optimized assay as shown in the representative recording, where the dots reflect the measured current from 540 to 550 ms after each 20 mV step from –100 to +100 mV, and three IV protocols are shown for each sample addition. Benzbromarone caused a dose-dependent inhibition of TMEM16A, which exhibited the typical outward rectification as shown in the representative instrument recording in (C) and the IV curves plotted in (D) of data from a separate recording showing currents measured at each voltage step during the third IV protocol.

Niclosamide, nitazoxanide and related compounds were found by QPatch electrophysiology to be potent inhibitors of the TMEM16A Ca2+-activated Cl current (Figure 4A). Significantly, while these compounds and 1PBC only partially inhibited the iodide/eYFP response (Figures 1H,I), they provided nearly complete inhibition of the chloride current as shown in the bar plots in Figure 4B and the dose-response curves of Figures 4C,D. Over 80% of the calcium-activated chloride current was inhibited by niclosamide, Compounds 2–5, nitazoxanide and 1PBC, which was a similar scope of block to that observed with the benchmark antagonist, CaCCinh-A01. This stands in stark contrast to the results from the TMEM16A/eYFP assay using iodide as the permeant anion where the maximal percent inhibition by niclosamide of 62.4 ± 12.0 (n = 39), as well as related compounds, was considerably less than that of CaCCinh-A01 which provided 102.1 ± 9.5 (n = 53) percent inhibition (Figure 1I). This would suggest their efficacy depends on the permeate anion, with there being much greater impact in blocking chloride versus iodide permeation. It’s also possible niclosamide could compete for chloride/iodide binding sites in the pore. The pharmacology of niclosamide and related compounds, like 1PBC, also contrasts with the benchmark antagonist benzbromarone, which fully inhibited both the iodide eYFP response and the TMEM16A chloride current. Importantly, however, 1PBC, niclosamide and related compounds, Cpd 2 and Cpd 4, fully inhibited the TMEM16A calcium-activated chloride current yet were inactive in blocking the cAMP-induced CFTR chloride current as measured by IonWorks Barracuda (IWB) electrophysiology (Figure 4E and Supplementary Figure 26). In contrast, the benchmark inhibitors benzbromarone, CaCCinh-A01 and MONNA were less selective as they reduced both TMEM16A and CFTR chloride currents (Figure 4E and Supplementary Figure 27).

FIGURE 4
www.frontiersin.org

Figure 4. Niclosamide and related compounds provide potent and nearly full inhibition of TMEM16A chloride currents. Whole cell patch clamp studies using 170 nM free intracellular calcium revealed the average potency of TMEM16A antagonists in blocking the TMEM16A outward current at +100 mV (A), the compounds maximum antagonist activity (B) and concentration-response relationship for inhibition (C,D). Niclosamide, Cpd 2 and Cpd4 showed potent inhibition of the TMEM16A calcium-activated chloride current (A–D) but spared the cAMP-induced CFTR chloride current (E) exhibiting < 20 percent of inhibition (POI) of the forskolin-induced current. (F) Shows representative current–voltage traces from TMEM16A whole cell patch clamp electrophysiology studies using 170 nM free intracellular calcium, where the y-axis is current in nA and the x-axis is membrane potential changes from –100 mV to +100 mV in 20 mV steps from Vhold of 0mV. While benzbromarone inhibited both the outward and inward currents (Figure 3), Cpd 4, 1PBC and niclosamide suppressed the outward current but showed concentration-dependent effects on the inward current ranging from inhibition to stimulation (F). The stimulation was especially noticeable for niclosamide and 1PBC at –100 mV and higher concentrations. In contrast, 1PBC (H) and Cpd 4 (G) were found to provide dose dependent inhibition (not stimulation) at –100 mV when tested by perforated patch clamp electrophysiology, where the ionomycin stim (10 μM) and compounds (from 0.01 to 30 μM) were added at the 0 s time point. The DMSO Vehicle Control (gray, open squares; mean ± SD, n = 8) included on each patch plate shows the rapid calcium-dependent activation and inactivation of TMEM16A described earlier. ∗∗P < 0.01, ∗∗∗P < 0.001; significantly different from benzbromarone (A,B) or niclosamide (B); ns, not significant (unpaired t-test). Panels (G,H), P < 0.05, ∗∗P < 0.01, #P < 0.001, $P < 0.0001 from two-way ANOVA for differences from the Vehicle Control; all other data points showed no significant difference (p > 0.05) compared to Vehicle. Number of replicates: n = 3–11 (A–D), n = 3–4 (E), mean ± SEM. Average of quadruplicate measures (G,H).

Compounds 2–4 were more potent inhibitors of the calcium-activated chloride current than the benchmark antagonist benzbromarone (Figure 4A). This finding is consistent with results from the TMEM16A eYFP assay (Figure 1G). Niclosamide, however, which stood out as the most potent compound in the halide-sensitive YFP assay with average IC50 of 0.132 ± 0.067 μM (n = 42), showed significantly less activity by QPatch electrophysiology with an average IC50 of 2.37 ± 2.11 μM (n = 7) (Figure 4A). A consideration of the two TMEM16A assays may provide an explanation for these differences. In the high-throughput eYFP assay there is a 30 min pre-incubation with compounds prior to ionomycin activation of the channel, while the automated QPatch electrophysiology assay has only about 1 min incubation time following each compound addition in the cumulative concentration response analyses (Figure 3). If there is insufficient time for the compound to block, one might expect to see currents still decreasing between the first and third I/V traces after each compound addition indicating steady-state block was not achieved. Supplementary Figures 621 provide representative QPatch instrument current-time recordings for niclosamide, nitazoxanide and related compounds, as well as, seven benchmark antagonists. While benzbromarone acts quickly (Figure 3C and Supplementary Figure 6) to block TMEM16A currents, we routinely observe currents continue to decline between the 1st and 3rd IV protocol after niclosamide addition (Supplementary Figure 12) indicating steady-state block is not achieved and its potency may be underestimated. Indeed, if we modify the protocol to look at single concentrations of niclosamide incubated for the maximum 12–15 min where we see stable recordings, we find for 1.11 μM niclosamide it took about 10 min to achieve steady state block (Supplementary Figure 13A) while approximately 4 min are needed when using the higher concentration of 3.3 μM (Supplementary Figure 13B). Similar studies on Compounds 2–4 (Supplementary Figure 1418) indicate they also act slowly, and their potency may be underestimated by our automated QPatch assay. Unfortunately, extending the incubation time after each compound or vehicle addition resulted in too much current run-down and use of just single concentrations per recording chamber was too limiting to assay throughput. Despite these technical limitations, the electrophysiology data in Figures 4A–D provided important evidence that niclosamide and related compounds were potent inhibitors of the TMEM16A Ca2+-activated Cl current.

As mentioned earlier, electrophysiology studies on TMEM16A are hampered by a strong calcium-dependent inactivation so we employed low intracellular calcium (170 nM free) and measured compound impact in blocking the outward current at +100 mV as little inward current is generated (Figure 3). Under physiological conditions, however, pro-contractile agonists would significantly increase intracellular calcium to the micromolar range resulting in a robust and biologically relevant TMEM16A inward current leading to depolarization and contraction of ASM cells. Interestingly, while only minor inward currents are observed when using low intracellular calcium and a holding potential of 0 mV, we found 1PBC and niclosamide caused a paradoxical increase in the inward current at −100 mV while robustly inhibiting the outward current in whole cell patch clamp electrophysiology studies (Figure 4F and Supplementary Figures 1113). The benchmark TMEM16A antagonist, niflumic acid, was also reported earlier to show such voltage-dependent effects on TMEM16A (Bradley et al., 2014; Liu et al., 2015). We have confirmed these results, showing that niflumic acid, just like niclosamide, is blocking the outward yet stimulating the inward current (Supplementary Figure 10). The benchmark antagonist NTTP showed a similar pharmacology (Supplementary Figure 19), while the niclosamide-related compounds, Cpd 3 and Cpd 4 inhibited the outward current but showed no significant effect in either stimulating or inhibiting the inward current (Figure 4F and Supplementary Figures 16, 18). In contrast, the benchmark antagonists benzbromarone, CaCCinh-A01, and MONNA inhibited both the outward and inward current (Figure 3D and Supplementary Figures 68).

While our recording conditions were similar to earlier whole cell patch clamp studies on TMEM16A, we were concerned that the non-physiological conditions of symmetrical chloride, low intracellular calcium and large hyperpolarization steps from 0 mV to −100 mV could be contributing to these paradoxical effects. Perforated patch clamp electrophysiology offers an advantage over whole cell patch clamp as it preserves many of the intracellular components of the cell. We therefore examined a subset of compounds for effects in blocking the inward current at −100 mV following ionomycin-activation of the TMEM16A chloride current. IonWorks Barracuda perforated patch clamp electrophysiology studies revealed a large TMEM16A inward current following ionomycin-activation of the channel, which rapidly inactivated (Figures 4G,H). 1PBC (30 μM) showed robust inhibition of this inward current (Figure 4H) which was in stark contrast to the effects observed by whole cell patch clamp studies (Figure 4F) where instead a slight stimulation was observed. Similarly, niflumic acid, which caused a dose-dependent stimulation of the inward current by whole cell patch clamp analysis (Supplementary Figure 10) instead by perforated patch clamp electrophysiology showed inhibition (not stimulation) of the inward current (Supplementary Figure 25C). These findings are consistent with the physiological effect of niflumic acid in blocking agonist-evoked and spontaneous-transient inward currents in airway smooth muscle (Liu and Farley, 1996; Gallos et al., 2013). In contrast, the benchmark TMEM16A antagonist benzbromarone was effective in blocking the inward currents in both assays (Figure 3D and Supplementary Figure 25B).

We would suggest the paradoxical increase in the inward current is not physiologically relevant, but instead an artifact of the whole cell patch clamp recording conditions that often involve huge depolarization steps. It is worth noting TMEM16A antagonists anthracene-9-carboxylic acid (A-9-C) and niflumic acid were shown earlier to inhibit the outward current but to cause an anomalous stimulation of the inward current that was dependent on a large prior depolarization step (Piper et al., 2002; Piper and Greenwood, 2003). The absence of us seeing stimulation of inward currents by NFA or 1PBC in perforated patch clamp may be due to our use of a differing holding potential (−30 mV) in these studies compared to whole cell patch clamp (0 mV), and thus there being lesser prior depolarization before stepping to −100 mV. In addition, as voltage dependence of TMEM16A is due to the pushing of cytosolic Ca2+ ions into the channel (by depolarized clamp voltage) as a result of lipid changes that facilitate better access of Ca2+ to intramembrane binding sites, this ‘artifact’ will be seen more readily in fast whole-cell patch clamp using pipette filling solutions with calcium compared to perforated patch clamp studies where cytosolic calcium is not compromised by pipette calcium. While further research will be required to clarify this, it seems clear that 1PBC, niclosamide and related compounds (including several hundred not described here; see Supplementary Figure 4) exhibit a distinct pharmacology from the benchmark antagonist benzbromarone which may reflect a different mechanism in blocking the channel. To further assess their physiological effects, we evaluated their efficacy in modulating airway smooth muscle depolarization and contraction following treatment with pro-contractiles.

TABLE 1
www.frontiersin.org

Table 1. Niclosamide and related compounds are potent inhibitors of TMEM16A conductivity that block histamine depolarization of human airway smooth muscle cells and relax mouse tracheal rings pre-contracted with carbachol.

TMEM16A Antagonists Block ASM Depolarization and Bronchodilate Airways

Pro-contractiles such as histamine and methacholine signal through Gq-coupled GPCRs on the cell surface to cause membrane depolarization and contraction of airway smooth muscle cells. To enable high-throughput measures of the impact of TMEM16A antagonists on ASM responses, we developed a membrane potential assay using primary human ASM cells. As demonstrated in the FLIPR-Tetra kinetic traces shown in Figure 5A, histamine causes a rapid increase in fluorescence in the membrane potential assay related to ASM depolarization that is fully suppressed by niclosamide and the TMEM16A antagonist, benzbromarone. Figure 5B provides representative dose-response results on six TMEM16A antagonists indicating they fully inhibited the depolarization induced by an EC90 amount of histamine (Supplementary Figure 28) and provided similar maximum inhibition compared to benchmark antagonists (Figure 5C). The niclosamide-related compounds, Cpd2 – Cpd4, exhibited significantly greater potency than the benchmark antagonists CaCCinh-A01 and niflumic acid in blocking ASM depolarization (Figure 5B). Heartened by these results we advanced compounds to tissue studies using mouse tracheal rings to evaluate their effect in bronchodilating airways. As shown in Figure 5D, niclosamide, Cpd 2 and Cpd 4 were highly potent (sub-μM) in providing full relaxation (bronchodilation) of carbachol precontracted airways. Tizoxanide, the metabolic product of nitazoxanide, also provided efficient bronchodilation of mouse airways. Table 1 provides the full results from studies on these compounds and additional TMEM16A antagonists, where it can be seen the enhanced efficacy of niclosamide and related compounds in blocking TMEM16A corresponds to improved potency in bronchodilating airways. For the first time, we also provide data exploring the effects of 1PBC on ASM physiology. While 1PBC to date has only been explored for effects on TMEM16A, our data would suggest it is highly effective in blocking depolarization and contraction of ASM cells (Figures 5B,D). NTTP, which is structurally related to tizoxanide (Figure 2), was also efficacious in relaxing mouse tracheal rings (Supplementary Figure 36). For over 50 years there has been no new bronchodilators, with β-agonists remaining the only agent and mechanism for blocking the multiple contractiles operating in disease. Antagonists of TMEM16A offer a new mechanism, thus it’s important to point out our studies with ASM cells using histamine and our studies with mouse trachea using the cholinergic carbachol support the idea that TMEM16A antagonists block the multiple contractiles that can operate in disease. To extend our bronchodilation studies to human tissue and smaller airways, we’ve also tested antagonists for relaxation of human 4th order bronchi. Niclosamide and Compound 4 fully relaxed human bronchial rings pre-contracted with carbachol (Figures 5E,F) or histamine (Figure 5H). Of the TMEM16A antagonists, niclosamide routinely provided the most potent bronchodilation. In general, its efficacy in relaxing human bronchi compared favorably to its efficacy in relaxing mouse tracheal rings where it averaged an EC50 of 0.379 μM (Table 1). Of the benchmark TMEM16A antagonists, 1PBC had the greatest potency in relaxing carbachol precontracted human bronchial rings with EC50 of 0.91 ± 0.20 μM, followed by benzbromarone, CaCCinh-A01, dichlorophen and niflumic acid with EC50 values of 2.08 ± 1.95, 9.24, 10.64 ± 12.0, and 14.40 ± 0.76 μM, respectively (Figures 5E–G and Supplementary Figure 38). The short-acting β-agonist salbutamol had an EC50 of 119 nM for bronchodilating histamine precontracted human airways that was significantly lower than the EC50 of niclosamide (967 nM) and Cpd 4 (3.45 μM), but a similar concentration of ∼10 μM was needed to fully relax the tissue (Figure 5H).

FIGURE 5
www.frontiersin.org

Figure 5. Niclosamide and related compounds block histamine depolarization of airway smooth muscle cells and provided robust bronchodilation of carbachol or histamine pre-contracted mouse and human airways. Representative kinetic traces of fluorescence recordings from the FLIPR membrane potential assay using primary human ASM cells (A). Histamine-induced a robust increase in fluorescence corresponding to depolarization which was fully blocked by the TMEM16A antagonists benzbromarone and niclosamide (A). Niclosamide and related compounds provided potent inhibition of the histamine-induced ASM depolarization (B), with the average maximum inhibition from dose-response studies exceeding 80% or <20 percent of control (POC) response to vehicle (Control); n = 2–15 (C). In wire myograph assays using carbachol pre-contracted mouse tracheal rings, niclosamide and related analogs, Cpd 2 and Cpd 4, provided sub-μM bronchodilation and improved potency over the two TMEM16A benchmark antagonists run in parallel (D). Niclosamide and Compound 4 also provided robust relaxation of human bronchial rings pre-contracted with carbachol, with results in panels (E,F) being wire myograph studies using rings from the same donor and involving simultaneous tests in parallel of the four compounds shown. Using a separate donor, effective bronchodilation of carbachol pre-contracted human airway rings was also observed with five benchmark TMEM16A antagonists of distinct chemotypes (G). Niclosamide and Compound 4, additionally provided robust relaxation of histamine pre-contracted human bronchial rings (n of 2–3) as shown in the wire myograph results in panel (H) on airway rings from different donors. The results from separate myograph studies on the β-agonist salbutamol are overlaid for comparison (n of 3 donors). As a precaution, all wire myograph studies using human bronchial rings included 5 μM indomethacin in the bath solution to eliminate any possible indirect effects due to prostaglandins. Mean ± SEM. Significant differences from benzbromarone or niclosamide (C) was determined with unpaired t-test (ns, not significant).

TMEM16A Antagonists Have a Slow Onset of Action, but Could Offer Sustained Effects in Bronchodilating Airways

The onset and duration of bronchodilation are two additional variables to explore when developing new asthma therapies. These can be examined with wire myograph by assessing force changes over time. The raw wire myograph traces of force changes over time for isoproterenol, benzbromarone, and niclosamide relaxation of human bronchial rings and the average concentration-response relationships are shown in Figure 6. The hydrophilic β-agonist isoproterenol (XLogP3 = −3.49) had a quick onset of action, relaxing the histamine pre-contracted human bronchial tissue in less than 7 min at 10 or 100 nM doses (Figure 6A). In contrast, the lipophilic TMEM16A antagonists benzbromarone and niclosamide (XLogP3 of 5.66 and 5.32, respectively) had a slow onset of action in relaxing human bronchial rings pre-contracted with histamine. Benzbromarone took 80–140 min at a 10 μM dose (Figure 6B), while niclosamide required ∼50 min at 10 μM and ∼92 min at a 1 μM before force changes reduced to a plateau (Figure 6C). Compound 4 is also lipophilic (XLogP3 = 4.63) and took between 50 and 125 min at 10 μM in three separate wire myograph recordings before fully relaxing the tissue (Supplementary Figure 40). Despite the significant difference in onset of action between the β-agonist isoproterenol and the TMEM16A antagonists, all compounds fully relaxed human bronchial rings pre-contracted with histamine (Figure 6D).

FIGURE 6
www.frontiersin.org

Figure 6. The hydrophilic β-agonist, isoproterenol, acts quickly to relax human bronchial rings while the lipophilic TMEM16A antagonists, benzbromarone and niclosamide, have a slow onset of action. Raw wire myograph recordings of force changes over time are shown for isoproterenol (A), benzbromarone (B), and niclosamide (C) relaxation of human bronchial rings pre-contracted with histamine and the average concentration-response relationship (D). Force changes after histamine addition are stable over time and unaffected by repeat additions of Vehicle control (A - top panel and C - bottom panel; vertical black dotted lines). Isoproterenol (ISO) is hydrophilic (XLogP3 = –3.49) and has a quick onset of action, taking less than ∼ 7 min at 10 or 100 nM concentrations to achieve steady-state relaxation of the human bronchial rings (A), with results from four donors (d1–d4) being shown. In contrast, the lipophilic TMEM16A antagonist benzbromarone (Benzbro) (XLogP3 = 5.66) has a slow onset of action (B) in relaxing bronchial rings from donors d1 and d2, taking 80–140 min at 10 μM before force changes have reached a plateau. Niclosamide which is also lipophilic, XLogP3 = 5.32, similarly acted slowly taking ∼50 min at 10 μM and ∼92 min at 1 μM to achieve steady-state relaxation of the bronchial tissue from donors d3 and d4 (C). To confirm responsiveness of each donor, isoproterenol is routinely run as a positive control. Mean ± SEM, N = 2–5.

Similar to the effects on human bronchial rings, benzbromarone, niclosamide and Cpd 4 had a slow onset of action in relaxing mouse tracheal rings pre-contracted with carbachol (Figure 7). Benzbromarone had the slowest onset of action taking 40 – 60 min at 1 μM to achieve steady-state relaxation of the tissue (Figure 7D), while niclosamide took 17 – 35 min at the same dose (Figure 7E) and Cpd 4 required 25 – 45 min at 316 nM (Figure 7F) before relaxation had reached a plateau. As was the case with human bronchial rings, the hydrophilic β-agonist isoproterenol acted quickly to relax mouse airways taking just 2 – 4 min to relax the tissue (Supplementary Figure 39). For the lipophilic TMEM16A antagonists, a technical limitation in wire myograph experiments due to the slow onset of action is a tendency to overlook activity at lower concentrations where insufficient time is allowed to relax the tissue. We think this can lead to compression of the dose-response curves and noticeably large Hill Slopes, such as those observed in Figures 7A–C.

FIGURE 7
www.frontiersin.org

Figure 7. The TMEM16A antagonists, benzbromarone, niclosamide, and Compound 4 fully relax mouse tracheal rings but have a slow onset of action. The average concentration-response relationship (A–C) and raw wire myograph recordings of force changes over time are shown for benzbromarone (A,D), niclosamide (B,E), and Compound 4 (C,F) for relaxation of mouse tracheal rings pre-contracted with an EC75 concentration of carbachol. The lipophilic compounds benzbromarone (XLogP3 = 5.66), niclosamide (XLogP3 = 5.32), and Cpd 4 (XLogP3 = 4.63) took an extended period of time to achieve steady-state relaxation of the tissue but had efficacy similar to theophylline in fully relaxing the mouse airways. At a 1 μM dose, benzbromarone (D) exhibited a slow onset of action taking between 40–60 min to achieve steady-state relaxation of the tissue, while niclosamide took between 17–35 min at the same dose (E). Cpd 4 also had a slow onset of action taking between 25–45 min at 316 nM to achieve steady-state relaxation (F). In contrast the hydrophilic β-agonist isoproterenol (XLogP3 = –3.49), while only partially relaxing airway rings pre-contracted with an EC75 level of carbachol, acted quickly taking just 2–4 min to bronchodilate the tissue (Supplementary Figure 39). To enable direct comparison of the activity of the molecules, compounds were run in parallel on the same day in separate chambers of the wire myograph. Theophylline was routinely added at the end of the recording to define full bronchodilation and calculate percent relaxation of the tissue. The results shown are from 3 to 5 separate experiments. Mean ± SEM.

To explore compound duration of action, we treated tissue with compounds, washed and then re-challenged the tissue with pro-contractiles to see if bronchodilation was sustained or transient. Preliminary tests of the β-agonist isoproterenol indicate its effects were short-acting as expected, while bronchodilation with niclosamide and benzbromarone were long lasting (Supplementary Figure 37).

Antagonists Resist Use- and Inflammatory-Desensitization Mechanisms

TMEM16A antagonists offer a new mechanism to bronchodilate airways that we’d expect to resist use- and inflammatory-desensitization mechanisms that limit β-agonist action. To test this hypothesis, we employed mouse tracheal rings which can be readily attained and compared TMEM16A antagonists versus the β-agonist isoproterenol for bronchodilation after various challenges meant to reflect the situation of severe asthma with poorly controlled inflammation and airway hyperresponsiveness (AHR), which reflects increased sensitivity to contractants. While β-agonists fully relax partially precontracted airways, it’s been shown earlier they only partially relax maximally precontracted airways (Lemoine and Overlack, 1992), which may be relevant to drug therapy and disease control. We have confirmed this desensitization phenomena using the β-agonist isoproterenol, but importantly find TMEM16A antagonists fully bronchodilate airways irrespective of the level of carbachol precontraction (Figure 8). To accomplish these studies, the EC25, EC50, EC75, and EC95 for carbachol bronchoconstriction is first determined for each airway ring (Supplementary Figure 33), the tissue is then washed to return tissue to baseline tension and then treated with various levels of contractant. As shown in Figure 8A, the β-agonist isoproterenol fully (100%) relaxed airway rings treated with an EC25 level of carbachol but caused only a partial 42% relaxation of airways treated with an EC95 level of carbachol. In contrast, the TMEM16A antagonists benzbromarone, niclosamide and Compound 4 fully relaxed maximally (EC95) contracted airways, provided roughly 100% bronchodilation irrespective of contractile force (Figures 8B–D). Thus, TMEM16A antagonists as an add-on therapy operating with a distinct MOA may help patients retain disease control when β-agonists are functioning sub-optimally. Indeed, an airway ring maximally contracted with an EC95 level of carbachol showing just 38% bronchodilation with 1 μM isoproterenol, achieved roughly 100% bronchodilation after subsequent treatment with the TMEM16A antagonist, Cpd 4 (Figure 8A).

FIGURE 8
www.frontiersin.org

Figure 8. TMEM16A antagonists fully relax airways in the presence of low and high concentrations of contractants, while the β-agonist isoproterenol showed only limited efficacy on maximally contracted airways. Mouse tracheal rings were pre-contracted with EC25 – EC95 concentrations of carbachol (CCh) and then treated with the β-agonist isoproterenol (A), or the TMEM16A antagonists benzbromarone (B), niclosamide (C) or Compound 4 (D) to measure bronchodilation using a DMT wire myograph (n = 3–5 rings per condition). A maximally contracted airway ring using an EC95 level of carbachol that showed only partial bronchodilation after the highest isoproterenol dose, was fully relaxed by add-on of the TMEM16A antagonist Cpd 4 (% relaxation shifts from 37.7% → 101.7% at EC95); see arrow to open blue square, cross-hatched (A). Mean ± SEM. P < 0.05, ∗∗P < 0.01, #P < 0.001, $P < 0.0001, two-way ANOVA with Dunnett’ s multiple comparison test for differences from the EC25 of CCh in each group. Only the isoproterenol group showed a significant difference in the dose-response relationship depending on the level of CCh pre-contraction, while no significant effects were observed for the benzbromarone, niclosamide or Cpd 4 groups.

Severe asthma and COPD is associated with poorly controlled inflammation. Using a cocktail of cytokines (IL-1β, TNFα, IL-13) to mimic this process, we show the efficacy of isoproterenol in bronchodilating airways is dramatically reduced when airways are pretreated with cytokines (Figure 9A). In contrast benzbromarone and niclosamide fully bronchodilated cytokine-treated airways (Figures 9B,C), suggesting TMEM16A antagonists offer a new mechanism to resist the inflammatory desensitization pathways that limit β-agonist action. Repeat treatment with β-agonists are also known to induce tachyphylaxis or use-dependent desensitization, which involves β2-adrenergic receptor phosphorylation, arrestin recruitment and receptor internalization and degradation resulting in reduced bronchodilation. To ascertain whether these pathways leading to β-agonist desensitization also interfere with bronchodilation by TMEM16A antagonists, we measured the degree of bronchodilation after repeated treatment with isoproterenol. As shown in Figure 9D, there was less bronchodilation and a rightward shift in potency after the 2nd round of β-agonist treatment which reflected use-dependent desensitization. Subsequent addition of Compound 4 to this same tissue resulted in full bronchodilation indicating TMEM16A antagonists restore bronchodilation of airways undergoing β-agonist use dependent desensitization (Figure 9D). The combined data indicates TMEM16A antagonists provide a new mechanism for robust relaxation of airway smooth muscle, even when they encounter high levels of contractants or are exposed to inflammatory conditions that limit β-agonist action.

FIGURE 9
www.frontiersin.org

Figure 9. TMEM16A antagonists fully bronchodilate airways and resist use- and inflammatory-desensitization pathways that limit β-agonist efficacy. Treatment of mouse tracheal rings overnight with a cytokine cocktail of IL-β, IL-13, and TNFα strongly reduced isoproterenol bronchodilation of carbachol pre-contracted airways (A), while the TMEM16A antagonists benzbromarone (B) and niclosamide (C) remained effective in fully relaxing the cytokine treated airways (n = 4; 2 rings/group, 2 experiments). As a model of β-adrenergic receptor use-dependent desensitization, carbachol pre-contracted (EC50 or EC75 level of CCh) mouse tracheal rings were treated with the β-agonist isoproterenol (round 1), washed, re-contracted with carbachol and tested again (round 2) for bronchodilation. Isoproterenol showed a rightward shift in potency after the second treatment with the β-agonist and a reduction in maximal efficacy in bronchodilation (D). Airway rings showing desensitization and only partial bronchodilation after this second treatment and highest isoproterenol dose, were fully relaxed by add-on of the TMEM16A antagonist Compound 4 (% relaxation shifts from 78.6% 98.9% at EC50 and from 42.7% 95.4% at EC75); see arrows to open cross-hatched circle and square (D). P < 0.05, ∗∗P < 0.01, #P < 0.001, $P < 0.0001 from two-way ANOVA with Sidak’ s multiple comparison test for differences between untreated and cytokine treated airway rings for each compound treatment group; all other data points showed no significant (ns) difference (p > 0.05). Mean ± SEM.

TMEM16A Expression in Human ASM and Bronchial Epithelium, Upregulation by Th2 Cytokines and Overexpression in Lungs of Monkeys With Allergic Asthma

Using an RNAseq dataset described earlier (Aisenberg et al., 2016), we evaluated TMEM16A expression in primary human airway smooth muscle cells and bronchial epithelial cells derived from healthy donors or patients with asthma or COPD. TMEM16A expression was higher in bronchial epithelial cells compared to ASM cells (Figure 10A), with the serial cultures of epithelial cells mostly comprised of basal cells. Expression was detected in both normal and diseased samples and confirmed at the sequence level by evaluating RNAseq read coverage (Supplementary Figure 41). Inflammatory cytokines can also act on the bronchial epithelium to induce goblet cell hyperplasia and mucin hypersecretion leading to airway narrowing and loss of lung function. To evaluate the effects of the Th2 cytokines, IL-4 and IL-13, on TMEM16A expression in the mature bronchial epithelium, we cultured normal and COPD human bronchial epithelial (HBE) cells at an air-liquid interface (ALI) to generate a polarized, pseudostratified bronchial epithelium composed of basal cells, secretory cells, and ciliated cells. The mature ALI cultures were then untreated or treated with IL-13 for 1, 3, 5, or 7 days and processed for RNAseq whole transcriptome sequencing to evaluate gene expression and alternative splicing. TMEM16A was dramatically upregulated after IL-13 treatment of HBE ALI cultures quickly reaching a peak after just 1 day (Figure 10B, left panel), while elevated expression of the Muc5AC mucin occurred more slowly and progressively from days 3 to 7 (Figure 10B, right panel). Importantly, we show human bronchial epithelial ALIs cultures treated with the cytokine IL-13 express the same TMEM16A ‘abc’ splice variant as untreated cultures suggesting Th2 inflammation does not induce alternative splicing (Supplementary Figure 41). Immunohistochemical analysis of separate ALI cultures untreated or treated with IL-4 or IL-13 for 2 days confirmed TMEM16A expression is also increased at the protein level (Figure 10C). Increased staining was observed at the apical surface (air interface, right side of ALI) and at the basolateral surface near the transwell membrane insert (liquid interface, left side of ALI) of the mature bronchial epithelial ALI cultures after treatment with the Th2 cytokines (Figure 10C). In the latter case elevated TMEM16A staining appeared to be of basal cells, while in the former case it appeared there was increased staining of mucin secreting goblet cells following treatment with the Th2 cytokines. To evaluate TMEM16A expression in the lung of naive or asthmatic cynomolgus (cyno) monkeys, tissue was obtained from naïve adult cynos that were non-allergic to Ascaris suum allergen and from adult cynos exhibiting repeated sensitivity to A. suum allergen that were characterized as asthmatic. Significant increases in TMEM16A expression was detected at the apical surface of the bronchial epithelium and in submucosal glands of cynomolgus monkeys with allergic asthma compared to naïve non-allergic animals (Figure 10D). Immunohistochemistry of lung tissue from a separate monkey with allergic asthma localized the increased TMEM16A staining to the apical surface of goblet cells in the bronchial epithelium (Figure 10E).

FIGURE 10
www.frontiersin.org

Figure 10. TMEM16A is expressed in cultured human airway smooth muscle and bronchial epithelial cells, quickly upregulated in mature human bronchial epithelial air-liquid interface (ALI) cultures after IL-4 or IL-13 treatment and overexpressed in the lung of asthmatic cynomolgus monkeys on the apical surface of goblet cells and in submucosal glands. TMEM16A (ANO1) expression in cultured human bronchial epithelial cells (HBE, submerged) and human airway smooth muscle cells from normal, asthmatic or COPD donors (A) was determined from an RNAseq dataset described earlier (Aisenberg et al., 2016). TMEM16A expression was higher in bronchial epithelial cells compared to airway smooth muscle cells. The results for ASM cells (normal or diseased) are from 3 separate donors each, while the results for HBE cells are the average response from 8 normal, 5 asthma and 3 COPD donors. The Th2 cytokine, IL-13, dramatically increased TMEM16A (ANO1) gene expression in mature bronchial epithelial ALI cultures reaching a maximum after 1 day (B, left panel), while elevated expression of the MUC5AC gene (B, right panel) occurred more slowly from day 1 through day 5 or 7. The RNAseq data shown are from ALI cultures from two normal (green; C1,C4) and three COPD (purple; C2, C3, C5) donors, untreated (Ctrl, open circles) or treated with IL-13 (closed circles). The Th2 cytokines, IL-4 and IL-13, also upregulated TMEM16A protein expression in human bronchial epithelial air-liquid interface cultures (C). Fully differentiated human bronchial epithelial ALI cultures were untreated (C, left panel) or treated for 48 h with IL-4 (C, middle panel) or IL-13 (C, right panel), fixed and stained with a rabbit anti-TMEM16A monoclonal antibody from Abcam (Ab64085). Increased TMEM16A staining (brown) was detected in basal cells located near the membrane insert and basolateral surface (left side of ALI culture) and on goblet cells located at the apical surface (right side) of IL-4 or IL-13 treated ALI cultures. Elevated TMEM16A protein expression was also detected in the bronchiolar epithelium and submucosal glands of cynomolgus monkeys with allergic asthma (D). Cynomolgus monkeys repeatedly challenged with Ascaris suum were compared to naïve non-allergic monkeys for TMEM16A expression using the polyclonal antibody Ab53212 (Abcam). Increased TMEM16A staining (brown) was observed on the apical surface of the bronchial epithelium and in submucosal glands of monkey with allergic asthma (D, right panel; animal #384630) compared to a naïve non-allergic monkey (D, left panel; animal #384606). Lung specimens from a separate asthmatic cynomolgus monkey (#384750) were collected 24 h post A. suum challenge. Immunohistochemical staining with the anti-TMEM16A antibody indicates TMEM16A is expressed selectively on the apical surface of goblet cells within the bronchial epithelium (E).

Antagonists Effects on Endogenous Calcium-Activated Chloride Currents and Identification of TMEM16A as Molecular Target for Niclosamide in Cancer

TMEM16A is located within the 11q13 amplicon, one of the most frequently amplified chromosomal regions in human cancer that correlates with a poor prognosis. To validate the pharmacology of antagonists in modulating native calcium-activated chloride currents, we profiled RNAseq data on large panel of cancer cell lines to identify high expressers. Of the 110 cell lines showing the most elevated TMEM16A expression (Supplementary Figure 43), roughly 18% were colorectal cancer cell lines and 6% were gastric cancer cell lines. We chose one of the colorectal cancer (CRC) cell lines, COLO205, for further evaluation by QPatch electrophysiology. COLO205 cells gave good seals and a robust TMEM16A current with the characteristic calcium and voltage-dependence. Application of three benchmark TMEM16A antagonists (benzbromarone, CaCCinh-A01 and 1PBC) fully blocked this native chloride current (Figures 11A–C and Supplementary Figures 22, 23). Importantly, the paradoxical pharmacology observed earlier in electrophysiology studies on 1PBC using cloned TMEM16A in HEK293 (Figure 4F), was also observed in studies with COLO205 cells expressing native TMEM16A currents. To our knowledge this represents the first description that COLO205 cells express a functional TMEM16A channel, but interestingly the colorectal cell line, HT-29, was exploited in one of the first screens to identify antagonists of the CaCC (De La Fuente et al., 2008) before its molecular identity was defined. While the confirmation of antagonist activity in blocking endogenous TMEM16A may seem trivial, notable differences have been observed between TMEM16A overexpressed in HEK293 cells and cell lines expressing TMEM16A endogenously (Schreiber et al., 2018), making it important to validate responses in both settings.

FIGURE 11
www.frontiersin.org

Figure 11. TMEM16A activity of additional niclosamide analogs and the impact of benchmark antagonists on endogenous calcium-activated chloride currents in the colorectal cancer cell line, COLO205. RNAseq identified several colorectal cancer cell lines as overexpressing TMEM16A (Supplementary Figure 43). Shown are results from QPatch electrophysiology studies indicating COLO205 cells exhibit the characteristic current–voltage relationship of TMEM16A which was sensitive to the TMEM16A antagonists benzbromarone (A), CaCCinh-A01 (B) and 1PBC (C), the later showing the same anomalous effects on the inward current observed with cloned TMEM16A in HEK293 cells. The structure activity relationship of several other niclosamide analogs is provided in (D). Drug repurposing has identified niclosamide as a promising new treatment for cancer, but its mechanism is not fully understood. Mook et al. (2015) have found niclosamide and related molecules suppress signaling and proliferation of cancer cells, including from the colon cancer cell line HCT-116. We find these compounds are also potent TMEM16A antagonists. Our compound numbering (Cpd#), the compound numbering by Mook et al. (2015) (RM#), and the calculated potency (IC50, μM) of the molecules in our TMEM16A halide-sensitive eYFP assay are shown (D).

Drug repurposing is an attractive approach to reduce the time and cost of drug development as it exploits preexisting clinical data on the pharmacokinetics, pharmacodynamics and safety of approved compounds. From such efforts niclosamide has emerged as a promising anticancer agent (Li et al., 2014). It inhibits proliferation of the colon cancer cell lines HCT116, COLO205 and HT29 with an IC50 of 0.37, 0.98, and 4.17 μM, respectively, which falls within the range of potency we find for blocking TMEM16A conductivity (Table 1). In a drug repositioning screen of 1600 clinical compounds using HCT116 spheroid cultures, Senkowski et al. (2015) identified niclosamide, as well as nitazoxanide, as promising new anticancer agents. Wnt/β-catenin signaling is a principle driver of colorectal carcinomas (Ahmed et al., 2016) and niclosamide was found to downregulate Wnt signaling and elicit colorectal antitumor responses (Osada et al., 2011). To optimize this chemotype further, Mook et al. (2015) synthesized over 35 niclosamide analogs to understand its SAR for inhibition of Wnt and to identify compounds with improved systemic exposure. A comparison of the 304 niclosamide analogs we’ve characterized for TMEM16A activity (Supplementary Figure 4) found 15 were identical to compounds they had characterized, were potent antagonists of both TMEM16A and Wnt signaling and had activity that correlated (Supplementary Figure 44). Niclosamide had an IC50 of 0.34 μM for inhibition of Wnt/β-catenin transcription using a TOPFlash reporter assay (Mook et al., 2015), which matches closely the IC50 of 0.13 μM we measured for inhibition of the TMEM16A, a molecular target also implicated in cancer. The molecular structure and activity of the fifteen niclosamide analogs is shown in Figure 11D. While niclosamide has a safe track record in man, compounds with a nitro group are often avoided due to potential toxicity concerns (Patterson and Wyllie, 2014). Removal of the nitro group of niclosamide resulted in only minor reduction in TMEM16A activity (Cpd 11), while equivalent potency was observed when it was replaced with a chlorine group (Cpd 12). Several additional derivatives were identified having sub-μM potency on TMEM16A. While none showed improved activity compared to niclosamide, understanding what structural substitutions are tolerated is equally important to identify derivates with improved solubility, pharmacokinetics and stability in vivo.

Discussion

For the first time, we identify the approved drugs niclosamide and nitazoxanide as TMEM16A antagonists. With an IC50 of 132 nM in blocking the TMEM16A halide sensitive YFP response, niclosamide is one of the most potent antagonists described to date. To provide perspective on their activity, niclosamide and related compounds were compared to eight benchmark TMEM16A antagonists described earlier. The pharmacology of niclosamide and related compounds in providing potent but partial block of the iodide eYFP response yet full block of the TMEM16A chloride current appears most similar to the benchmark antagonists 1PBC and NTTP. Thus, niclosamide and related molecules may share their mechanism as pore blockers. While we describe the activity of seventeen niclosamide analogs, over 300 niclosamide analogs have been tested and the partial block of the iodide eYFP response is characteristic of this chemotype. Further investigation is needed to understand why these compounds show partial inhibition of the iodide/eYFP response, but our findings are consistent with those of Peters et al. (2015) who similarly found 1PBC and NTTP only partially inhibited the TMEM16A iodide/eYFP response but fully inhibited the calcium-activated chloride current (see supplementary section of their online article). Our detailed evaluation of niclosamide SAR and the background literature on this anthelmintic motivated us to also analyze the approved drug nitazoxanide for TMEM16A activity since nitazoxanide was synthesized in early 1970s using the niclosamide scaffold, had improved oral bioavailability and exhibited no significant impact on cytochrome P450 enzymes (White, 2004; Anderson and Curran, 2007; Rossignol, 2014). Remarkably, nitazoxanide and its metabolic product tizoxanide also blocked TMEM16A with low single digit μM potency and showed a similar pharmacology to niclosamide in the eYFP and ephys assays. Upon further scrutiny of these compounds, we were surprised when we realized tizoxanide is a highly similar structural analog of NTTP, a known TMEM16A antagonist. We suggest this provides further evidence they may share this mechanism as pore blockers and propose they may interact with the same charged residues, Arg 511 and Lys 599, which Peters et al. (2015) find are important for NTTP and 1PBC inhibition of TMEM16A. These two basic residues along with Arg 531 which upon mutation to alanine alter anion selectivity of TMEM16A (Peters et al., 2015), have been confirmed recently in the cryoEM structures of mouse TMEM16A to reside near the external narrow neck of the channel pore where the positive charges facilitate anion entrance and permeability (Dang et al., 2017; Paulino et al., 2017a,b). Whether these compounds also compete for anion binding to this pore region, warrants further study. To our knowledge the mechanism of action for the benchmark TMEM16A antagonists CaCCinh-A01, benzbromarone, dichlorophen and MONNA remains unknown, but because these compounds show a distinct pharmacology from niclosamide and related compounds it is likely they have different site/s of action. Niflumic acid which has a pharmacology similar to niclosamide, 1PBC and NTTP in partially blocking the halide-sensitive YFP response, was proposed over a decade earlier to be a pore blocker (Piper et al., 2002) and to preferentially block smaller anions such as chloride compared to larger anions like iodide (Ni et al., 2014). At this time, we’ve run only limited counterscreens for selectivity over other chloride channels, but promisingly find niclosamide and related compounds showed little impact on the CFTR chloride current, as did the pore blockers 1PBC and niflumic acid.

TMEM16A is expressed on airway smooth muscle cells which display the classical electrophysiology properties of CaCC (Gallos et al., 2013; Zhang et al., 2013), including a rapid calcium-dependent inactivation (Wang and Kotlikoff, 1997) which we also observe after ionomycin-activation of cloned TMEM16A in HEK293 cells (Figures 4G,H and Supplementary Figure 25). Alternative splicing has been reported to be an important mechanism to regulate the voltage and Ca2+ dependence of TMEM16A chloride currents (Ferrera et al., 2009). We have confirmed TMEM16A expression using an RNAseq dataset (Aisenberg et al., 2016) of primary human ASM and bronchial epithelial (HBE) cells (Figure 10), detecting almost exclusively the ‘abc’ splice variant lacking segment ‘d’ in untreated or IL-13 treated HBE cells, but find ASM cells express a mixture of splice variants with variable amounts of segment ‘b’ and ‘d’ (Supplementary Figure 41). In this later case, this may include the ‘acd’ variant we used for QPatch electrophysiology studies, while in the former it would include the ‘abc’ variant originally identified in bronchial epithelial cells (Caputo et al., 2008) that was used in our high-throughput TMEM16A eYFP screen. While we’ve primarily used the ‘acd’ variant for patch clamp electrophysiology studies on antagonists, it’s worth noting the pore blocker 1PBC showed similar efficacy in blocking TMEM16A ‘acd’ and ‘abc’ Ca2+-activated Cl currents (Figure 4C and Supplementary Figures 11, 24). Irrespective, since pro-contractiles such as histamine and methacholine activate CaCCs leading to membrane depolarization, calcium mobilization and contraction of smooth muscle cells, we advanced TMEM16A antagonists for direct interrogation of impacts on these processes.

Niclosamide and related compounds blocked histamine depolarization of ASM cells at sub-micromolar concentrations. These studies confirm earlier findings that TMEM16A antagonists hyperpolarize human ASM cells (Yim et al., 2013; Danielsson et al., 2014; Danielsson et al., 2015) and extend them to demonstrate antagonists also reverse the depolarization induced by contractants. To provide perspective on the activity of niclosamide and related compounds, we compared them to seven small molecule antagonists described earlier of distinct chemotypes (Figure 2). Consistent with TMEM16A’s role in controlling ASM excitation-contraction coupling, we find antagonists characteristically blocked pro-contractile depolarization and importantly provided robust bronchodilation of mouse tracheal rings and human small airways (Figure 5). In alignment with their enhanced potency on TMEM16A, niclosamide and related compounds showed improved activity in relaxing airways. Our findings extend earlier studies demonstrating niflumic acid relaxes TEA or K+ gluconate pre-contracted guinea pig tracheal rings, rat small airways and human tracheal smooth muscle strips (Yim et al., 2013; Danielsson et al., 2014), benzbromarone and tannic acid suppress Substance P contraction of guinea pig tracheal rings (Gallos et al., 2013) and benzbromarone relaxes methacholine-contracted mouse peripheral airways or human airway smooth muscle strips pre-contracted with acetylcholine or LTD4 (Danielsson et al., 2015). Importantly, these findings and our results using carbachol and histamine indicate TMEM16A antagonists offer a new mechanism to block the multiple contractiles operating in disease. More recently using smooth-muscle-specific TMEM16A knockout mice, Wang et al. (2018) demonstrate agonist-induced inward currents and contractility of airway smooth muscle is profoundly reduced in the absence of TMEM16A, cementing its role as a key regulator of ASM excitation and contraction.

Asthma is a chronic inflammatory disorder characterized by airway hyperresponsiveness (AHR). Significantly, Zhang et al. (2013) find in mouse model of asthma, ovalbumin-sensitization upregulates TMEM16A expression in ASM cells and benzbromarone and niflumic acid prevented AHR and contractions evoked by methacholine. We also have devised schemes to evaluate the efficacy of TMEM16A antagonist under extreme conditions of contraction and inflammation that reduce β-agonist action (Lemoine and Overlack, 1992; Shore et al., 1997). Remarkably, maximally contracted mouse airways were fully relaxed by the TMEM16A antagonists niclosamide, Cpd 4 and benzbromarone, while the β-agonist isoproterenol showed only partial bronchodilation which importantly converted to full bronchodilation with add-on of a TMEM16A antagonist (Figure 8). Severe asthma and COPD are also associated with poorly controlled inflammation, steroid insensitivity and reduced β-agonist responsiveness. Thus, there is a need for new bronchodilators with a distinct MOA that can operate amidst this inflammation. The cytokines IL-1β and TNFα were shown earlier to reduce β-agonist relaxation of airway smooth muscle (Shore et al., 1997; Moore et al., 2001). More recently bitter taste receptors were shown to provide a new mechanism to bronchodilate airways which resists IL-13 desensitization events limiting β-agonist action (Robinett et al., 2014). Using mouse tracheal rings treated with a more robust new cocktail of three cytokines, IL-1β, IL-13, and TNFα, we demonstrate β-agonist bronchodilation is profoundly reduced, yet TMEM16A antagonists remain remarkably effective in fully relaxing airways (Figure 9). Finally, because asthma treatment guidelines include β-agonists as a standard of care but allow for add-on therapies for poorly controlled disease, we investigated if TMEM16A antagonists can operate with β-agonists, including under situations of tachyphylaxis or use-dependent desensitization. Promisingly, airways exhibiting hyporesponsiveness after repeat β-agonist treatment were fully relaxed upon addition of the TMEM16A antagonist, Cpd 4. Future studies to explore whether niclosamide and related compounds attenuate allergen-challenged and sensitized airways are warranted given the findings above. While the cytokines IL-1β, TNFα, and IL-13 were chosen as a surrogate of inflammation, it should be noted they do not fully reflect the cytokine networks operant in asthmatic and/or COPD airways. In addition, while TMEM16A antagonism or deficiency was found earlier to ameliorate bronchial hyperresponsiveness in the Ova mouse model of allergic asthma (Zhang et al., 2013; Wang et al., 2018) that is dependent on eosinophils and Th2 inflammation, it remains unexplored in newer models of asthma exacerbations and COPD that are steroid-resistant (Ito et al., 2008; Fricker et al., 2014; Jones et al., 2017) in which neutrophils and/or macrophages serve a greater role. Further pre-clinical studies in these models would help provide deeper insights into the utility of TMEM16A antagonists for these disorders, as well as, severe non-eosinophilic/neutrophilic asthma (Moore et al., 2014; Esteban-Gorgojo et al., 2018) for which there is an urgent need for new therapies.

Secretory epithelial cells, smooth muscle cells and sensory neurons express the CaCC TMEM16A. In the bronchial epithelium, TMEM16A is strongly upregulated by Th2 cytokines which correlates with goblet cell hyperplasia and mucin hypersecretion, as well as, increased calcium-dependent chloride secretion (Caputo et al., 2008; Huang et al., 2012; Scudieri et al., 2012). We have confirmed these results finding TMEM16A is dramatically upregulated in the bronchial epithelium after IL-13 treatment and expressed on the apical surface of goblet cells, but also show that cynomolgus monkeys with allergic asthma have elevated TMEM16A staining within the submucosal glands (Figure 10), an upregulation also strikingly noted in submucosal glands of cystic fibrosis patients and cultured bronchial epithelial cells treated with bacterial pyocyanin or with a P. aeruginosa culture supernatant (Caci et al., 2015). Interestingly, like our findings with IL-13, Gorrieri et al. (2016) similarly find TMEM16A expression peaks quickly, 24 h after IL-4 treatment, a time frame we find precedes maximum Muc5AC expression occurring 5–7 days later (Figure 10). Notably, using an IL-13 mouse model of asthma, Nakano et al. (2006) reported that niflumic acid inhibited goblet cell hyperplasia and airway hyperresponsiveness. Huang et al. (2012) extended these findings demonstrating TMEM16A is upregulated in the asthmatic bronchial epithelium and its inhibition with benzbromarone and dichlorophen blocked mucin secretion and ASM contraction. More recently, TMEM16A inhibition with antagonists or siRNA has been reported to reduce IL-13 induced Muc5AC production and goblet cell hyperplasia (Lin et al., 2015; Zhang et al., 2015; Qin et al., 2016), while TMEM16A overexpression had opposite effects in increasing Muc5AC expression (Lin et al., 2015; Qin et al., 2016). Thus, TMEM16A inhibitors may repress two of the key features leading to airway obstruction in severe asthma and COPD, bronchoconstriction and mucin hypersecretion.

In contrast for cystic fibrosis, activation of TMEM16A on the apical membrane of epithelial cells has been proposed as one strategy to bypass the CFTR chloride secretory defect to improve hydration and mucociliary clearance, but bronchoconstriction and pain have been pointed out as possible unwanted side effects (Sondo et al., 2014; Liu et al., 2016). To evaluate the effects of direct channel activation, we employed the TMEM16A opener, Eact (Namkung et al., 2011), and investigated its effects on airway smooth muscle and the bronchial epithelium. Interestingly, Eact induced robust ASM cell depolarization and calcium flux like contractants, which was blocked by antagonists (Supplementary Figures 31, 32). The average EC50 of 1.6 μM for Eact induced ASM depolarization, closely matched its EC50 of 3.0 μM in activating TMEM16A (Namkung et al., 2011). These findings are consistent with those from Danielsson et al. (2015), who found Eact depolarized ASM cells like contractants, while antagonists hyperpolarized ASM cells. To investigate its effects on bronchial epithelial cells, we incubated mature bronchial epithelial air-liquid interface (ALI) cultures for 2–3 days with the Eact opener. Muc5AC and leukotriene (LTC4, LTB4) secretion from bronchial epithelial ALI cultures was significantly increased (2–3 fold and 6–12 fold, respectively, p < 0.05) following treatment with the TMEM16A opener, Eact (Supplementary Figure 42). While these findings supported the concept that agonists may have reciprocal effects to antagonists, our enthusiasm with these results was tempered by recent findings Eact was also a direct activator of TrpV1 (Liu et al., 2016), which can physically associate with TMEM16A in sensory neurons to enhance pain sensation (Takayama et al., 2015). As TrpV1 is also expressed by bronchial epithelial and ASM cells (McGarvey et al., 2014; Yocum et al., 2017), it’s activation alone, or with TMEM16A, could also explain some of our findings with Eact. Interestingly, Benedetto et al. (2017) also recently report TMEM16A and CFTR physically interact and cross-regulate each other in differentiated epithelial cells. Tissue specific knockout of TMEM16A in mouse intestine and airways abolished not only Ca2+-activated Cl secretion, but also abrogated CFTR-mediated Cl secretion. Somewhat surprising, there was no overt phenotype despite the complete absence of chloride currents in knockout tissues and mucociliary clearance was not compromised, but unexpectedly enhanced in tracheas from knockout mice. While this may support airway Na+ transport is physiologically more relevant for hydration and TMEM16A knock-down may enhance mucociliary clearance, further studies are needed to understand TMEM16A’s role under disease conditions and the effects of pharmacological agents in modulating goblet cell hyperplasia and mucin hypersecretion. Niclosamide probably inhibits TMEM16A in mucus producing cells which inhibits mucus production and/or secretion, however, there remains a need for TMEM16A openers to explore the consequence of direct channel activation.

Increased airway smooth muscle mass, goblet cell hyperplasia, mucus plugs, and subepithelial matrix deposition are a hallmark of the remodeled airway of severe asthmatics which is thought to underlie airway hyperresponsiveness and lead to obstruction (Holgate et al., 2015). Bronchial thermoplasty is a successful interventional approach for severe refractory asthma that uses thermal energy to reduce the excess airway smooth muscle mass, which decreases exacerbations and improves asthma-related quality of life (Laxmanan and Hogarth, 2015). It is interesting to speculate if TMEM16A antagonists by targeting airway smooth muscle might have similar effects and reduce the need for thermoplasty? Indeed, the identification of disease-modifying asthma therapies (DMAT) was recommended in 2015 by an NHLBI Working Group as a key area for future research directions in asthma (Levy et al., 2015). The authors indicate such treatments should target fundamental pathobiological mechanisms involved in asthma and/or alterations in lung structure, but they note it will be important to determine the critical window for disease modification and clarify what end-effector mechanisms (e.g., smooth muscle structure or function, mucus production, or airway fibrosis) are reversible or modifiable. Whether airway remodeling in severe disease is reversible or not is still actively debated, however, Hassan et al. (2010) found increased levels of proliferating ASM cells in tissue of severe asthmatics versus healthy controls and proposed airway smooth muscle remodeling is an active process in long-standing severe disease and not just an end-stage phenomenon. In a clinical study in mild asthmatics, Grainge et al. (2011) provided evidence methacholine-induced airway smooth muscle bronchoconstriction alone in the absence of inflammation led to the same remodeling as that observed following an eosinophilic inflammatory response to dust mite allergen. These findings were conspicuous as they illustrated the importance of effective control of bronchoconstriction to prevent pathological remodeling (Tschumperlin, 2011). We’d expect TMEM16A antagonists to act quickly to relieve bronchoconstriction and suppress mucin hypersecretion leading to improved lung function, but it would be worth investigating if they also reduce airway smooth muscle mass (or mucus plugs) over time and are disease modifying. The advent of several new imaging approaches (e.g., QCT, OR-OCT, MDCT) makes the identification of such disease-modifying asthma therapies, that modify or reverse disease rather than simply relieving symptoms, an exciting new possibility (Adams et al., 2016; Berair et al., 2017; Dunican et al., 2018).

Beyond respiratory disease, antagonists of TMEM16A have been proposed to have utility in treating a wide variety of other diseases including, pulmonary hypertension (Namkung et al., 2010; Forrest et al., 2012; Heinze et al., 2014), secretory diarrhea (Ousingsawat et al., 2011; Thiagarajah et al., 2015), polycystic kidney disease (Buchholz et al., 2014), pain (Cho et al., 2012; Lee et al., 2014; Pineda-Farias et al., 2015) and cancer. TMEM16A has been implicated in numerous cancers, including head and neck, esophageal, gastric, colorectal, lung, breast, prostate, pancreatic and glioma [reviewed recently by Wang et al. (2017)]. TMEM16A overexpression contributes directly to tumorigenesis and cancer progression and its knockdown with shRNA or small molecule antagonists ameliorates disease (Duvvuri et al., 2012; Liu et al., 2012; Britschgi et al., 2013; Wang et al., 2017). It associates with EGFR to regulate cancer cell proliferation and activates several pathways important for tumor growth (Britschgi et al., 2013; Bill et al., 2015). Interestingly, drug repurposing efforts have similarly found the anthelmintics niclosamide and nitazoxanide hold great promise for the treatment of cancer, as well as numerous other disorders, including hypertension, secretory diarrhea, and as a broad-spectrum anti-infective agent (Rossignol et al., 2006; Li et al., 2014; Rossignol, 2014; Chen et al., 2017). Underlying mechanisms proposed for its action have included uncoupling of oxidative phosphorylation, modulation of Wnt/β-catenin, mTORC1, STAT3, NF-κB and Notch signaling, but identifying a unifying mechanism for niclosamide’s action has remained elusive (Chen et al., 2017). We’d suggest antagonism of TMEM16A should now be considered as a possible contributing factor. This is especially so for secretory diarrhea and cancer where the role of TMEM16A has been well established. Activation or upregulation of the CaCC TMEM16A, an upstream target on the cell surface, can initiate a variety downstream signaling pathways affecting various pathophysiological processes [e.g., as shown for Ras-Raf-MEK-ERK in cancer (Duvvuri et al., 2012)]. Receptor-mediated increases intracellular calcium have also been found to be Cl dependent and TMEM16A antagonists have been found to modulate calcium homeostasis (Danielsson et al., 2015; Cabrita et al., 2017). Channels can also associate with other cell surface molecules to form large macromolecular signaling complexes. Indeed, TMEM16A has been shown to interact with the ezrin-radixin-moesin network (Perez-Cornejo et al., 2012) and with TrpV1 (Takayama et al., 2015), TrpC6 (Wang et al., 2016) and the IP3 receptor (Jin et al., 2013) for coupled GPCR and ion channel signaling. Interestingly, the key signaling molecule β-catenin involved in numerous cancers was also found recently to associate with the ion channels KCNQ1 (Rapetti-Mauss et al., 2017) and BKCa (Bian et al., 2011) to modulate Wnt signaling.

It is important to mention some limitations of our research. Our work was focused in broadly evaluating many compounds to develop a structure-activity relationship and advance a lead, as opposed to a more narrow but deep understanding of the pharmacology and mechanism of any single compound. Future studies to better understand niclosamide and nitazoxanide’s impact on permeation of different anions, their activity on native chloride channels, specificity and mechanism of block, and activity in animal models of disease, are warranted. This is especially so given the wide interest in these compounds from a drug repositioning standpoint and our new discovery they are TMEM16A antagonists. The bioavailability and suitability of these compounds for respiratory disease also needs to be explored in greater detail. While the short-acting β-agonists (SABA) salbutamol and isoproterenol are hydrophilic (XLogP3 of 0.31 and −3.49, respectively) and act quickly (<7 min) to open airways which is necessary for a rescue treatment, the TMEM16A antagonists niclosamide, Cpd 4 and benzbromarone are quite lipophilic (XLogP3 of 5.32, 4.63, and 5.66, respectively), have a slow onset of action (17–60 min) but may show sustained effects in bronchodilating mouse airways. The long-acting β-agonist (LABA) salmeterol is also quite lipophilic (XLogP3 of 3.90) and has been similarly observed to have a slow onset of action (90 min) but sustained effects in relaxing airways (Stocks et al., 2014). Such compounds, however, can be challenging from a development standpoint as they tend to have poor solubility, but they have ideal characteristics as an inhalation maintenance therapy since they offer long-lasting effects and are retained in the lung which is a safety advantage. Whether niclosamide and related compounds have similar advantages should be examined further, including if they share the long duration of action of the LABA salmeterol, which is believed principally to be the result of a physicochemical interaction of the lipophilic compound with membrane lipid bilayers acting as a depot to allow drug rebinding, the so-called ‘plasmalemma diffusion microkinetic model’ (Anderson et al., 1994). In addition, although these compounds operate with a distinct mechanism (TMEM16A) from β-agonists or other respiratory medicines (e.g., anti-cholinergics, leukotriene receptor antagonists, corticosteroids), future development should include formal drug-interaction studies. An TMEM16A antagonist would be best positioned initially as an add-on maintenance therapy for severe asthmatics whose disease is poorly controlled with their current asthma medicines. Based on GINA guidelines this would be at Step 5, where add-on therapy with the anti-muscarinic tiotropium and the biologics targeting IgE, IL-5 and more recently IL-4/IL-13 are allowed. Optimizing the formulation and delivery of these compounds to the lung will be another important activity (Borghardt et al., 2018; Strong et al., 2018), but the successful development of other lipophilic drugs (e.g., salmeterol, indacaterol, vilanterol) for respiratory disease, suggests effective formulation and inhalation delivery is possible for such compounds, even for diseases of the small airways, such as COPD. While the majority of TMEM16A antagonists were lipophilic, nitazoxanide and tizoxanide with XLogP3 values of 2.04 and 3.15, respectively, were unique in being somewhat more hydrophilic. These compounds merit further explorations as they may offer some development advantages. Additional research to better understand the molecular interaction of antagonists with TMEM16A is also warranted to improve compound activity and clarify their mechanism in blocking the channel.

TMEM16A antagonists offer an exciting new strategy to treat a wide variety of disorders. We demonstrate antagonists block airway smooth muscle depolarization and contraction, offering a promising new mechanism to bronchodilate airways that resists use- and inflammatory-desensitization pathways limiting β-agonist efficacy in severe asthma and COPD. We identify the approved drugs niclosamide, nitazoxanide and related compounds as some of the most potent TMEM16A antagonist described to date and provide a molecular target to compounds of high interest from drug repurposing efforts. Further studies will be needed to assess their efficacy in these new disease indications, but niclosamide and nitazoxanide’s safe track record in man position these chemical series as a good starting point for additional studies.

Ethics Statement

Animal experimentation and tissue: All animal procedures and tissue collections were conducted in an Association for Assessment and Accreditation of Laboratory Animal Care accredited facility in accordance with the requirements and guidelines of the US National Research Council and in compliance with the protocols approved by the Institutional Animal Care and Use Committee of the following institutions: Amgen Inc., Thousand Oaks, CA, United States; Charles River, Shrewsbury, MA, United States. All human lung specimens were collected under Institutional Review Board approval. Written informed consent was obtained from all participants in this study. In all cases, materials obtained were surplus to standard clincal practice. Patient identity and PHI/identifying information were redacted from tissues and clinical data.

Author Contributions

AH, KL, and KW developed and conducted the high-throughput screen and performed follow-up experiments to identify validated hits. PW and DP performed medchem support in determining compound potency in eYFP assay. BL, KH, JM, and JS developed TMEM16A electrophysiology assays, performed experiments, and/or analyzed the data. KG and X-ZX developed and analyzed compounds for effects on ASM membrane potential or calcium flux. GY characterized compounds for effects on HBE ALI cultures. KM, RE, DM, and AL developed and characterized compounds for bronchodilation of mouse, cynomolgus monkey, or human tissue. KM ran routine bronchodilation studies to support medchem and performed studies evaluating compound ability to resist use- and inflammatory-desensitization. ET, KH, SE, OH, KM, BB, JW, and MF conducted experiments evaluating TMEM16A expression and splice variants. JS, DS, MC, TB, AB, JJ, JP, MW, and KK were involved in conception, data analysis or enablement. JC and LL synthesized compounds and analyzed the data. JS initiated the project, supervised the program, and wrote the first draft of the manuscript; while AH and KK made critical edits on the final draft. All authors discussed and commented on the manuscript.

Funding

The studies were funded by Amgen Inc., but the corporate funders had no role in study design, data collection, or the decision to submit the work for publication.

Conflict of Interest Statement

KM, KL, BL, PW, KH, KG, RE, JC, LL, AL, ET, KH, X-ZX, OH, BB, MF, JW, GY, SE, MW, TB, AB, MC, DS, JP, JJ, JM, KW, DP, DM, AH, and JS are current or past Amgen employees who may own shares in the company. TB is employed by Sartorius, AB by Lilly, ET and KH by Seattle Genetics, KL by Clariant, DS by Insight Bioconsulting and KW by Assay.Works. AH is affiliated with ISAR Bioscience GmbH and JS with PolestarBio LLC but declare no commercial or financial interests in the subject material disclosed in this manuscript.

The remaining author declares that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Acknowledgments

The authors wish to thank the Amgen Protein Technology group for TMEM16A expression constructs, the scientific team at Sophion and ChanTest (CRL) for help in developing the QPatch and IonWorks Barracuda electrophysiology assays or counterscreens against CFTR, scientists at Biopta (ReproCELL) for bronchodilation studies, and past members of the Airway Smooth Muscle Working Group, Asthma, COPD and IPF Teams at Amgen who contributed to building a foundation in the respiratory efforts which enabled these studies, including Heather Thomas, Ken Schooley, Ryan Brown, Deb Hopkins, Anna Pirrone, Leanne Peiser, Cindy Willis, Daniela Metz, and Heather Arnett. This manuscript has been released as a pre-print at bioRxiv (Miner et al., 2018).

Supplementary Material

The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fphar.2019.00051/full#supplementary-material

References

Adams, D. C., Hariri, L. P., Miller, A. J., Wang, Y., Cho, J. L., Villiger, M., et al. (2016). Birefringence microscopy platform for assessing airway smooth muscle structure and function in vivo. Sci. Transl. Med. 8:359ra131. doi: 10.1126/scitranslmed.aag1424

PubMed Abstract | CrossRef Full Text | Google Scholar

Ahmed, K., Shaw, H. V., Koval, A., and Katanaev, V. L. (2016). A second WNT for old drugs: drug repositioning against WNT-Dependent Cancers. Cancers (Basel) 8:66. doi: 10.3390/cancers8070066

PubMed Abstract | CrossRef Full Text | Google Scholar

Aisenberg, W. H., Huang, J., Zhu, W., Rajkumar, P., Cruz, R., Santhanam, L., et al. (2016). Defining an olfactory receptor function in airway smooth muscle cells. Sci. Rep. 6:38231. doi: 10.1038/srep38231

PubMed Abstract | CrossRef Full Text | Google Scholar

Anderson, G. P., Linden, A., and Rabe, K. F. (1994). Why are long-acting beta-adrenoceptor agonists long-acting? Eur. Respir. J. 7, 569–578.

PubMed Abstract | Google Scholar

Anderson, V. R., and Curran, M. P. (2007). Nitazoxanide: a review of its use in the treatment of gastrointestinal infections. Drugs 67, 1947–1967.

Google Scholar

Barretina, J., Caponigro, G., Stransky, N., Venkatesan, K., Margolin, A. A., Kim, S., et al. (2012). The cancer cell line encyclopedia enables predictive modelling of anticancer drug sensitivity. Nature 483, 603–607. doi: 10.1038/nature11003

PubMed Abstract | CrossRef Full Text | Google Scholar

Benedetto, R., Ousingsawat, J., Wanitchakool, P., Zhang, Y., Holtzman, M. J., Amaral, M., et al. (2017). Epithelial chloride transport by CFTR requires TMEM16A. Sci. Rep. 7:12397. doi: 10.1038/s41598-017-10910-0

PubMed Abstract | CrossRef Full Text | Google Scholar

Berair, R., Hartley, R., Mistry, V., Sheshadri, A., Gupta, S., Singapuri, A., et al. (2017). Associations in asthma between quantitative computed tomography and bronchial biopsy-derived airway remodelling. Eur. Respir. J. 49:1601507. doi: 10.1183/13993003.01507-2016

PubMed Abstract | CrossRef Full Text | Google Scholar

Bian, S., Bai, J. P., Chapin, H., Le Moellic, C., Dong, H., Caplan, M., et al. (2011). Interactions between beta-catenin and the HSlo potassium channel regulates HSlo surface expression. PLoS One 6:e28264. doi: 10.1371/journal.pone.0028264

PubMed Abstract | CrossRef Full Text | Google Scholar

Bill, A., Gutierrez, A., Kulkarni, S., Kemp, C., Bonenfant, D., Voshol, H., et al. (2015). ANO1/TMEM16A interacts with EGFR and correlates with sensitivity to EGFR-targeting therapy in head and neck cancer. Oncotarget 6, 9173–9188. doi: 10.18632/oncotarget.3277

PubMed Abstract | CrossRef Full Text | Google Scholar

Borghardt, J. M., Kloft, C., and Sharma, A. (2018). Inhaled therapy in respiratory disease: the complex interplay of pulmonary kinetic processes. Can. Respir. J. 2018, 2732017. doi: 10.1155/2018/2732017

PubMed Abstract | CrossRef Full Text | Google Scholar

Bradley, E., Fedigan, S., Webb, T., Hollywood, M. A., Thornbury, K. D., McHale, N. G., et al. (2014). Pharmacological characterization of TMEM16A currents. Channels (Austin) 8, 308–320.

Google Scholar

Britschgi, A., Bill, A., Brinkhaus, H., Rothwell, C., Clay, I., Duss, S., et al. (2013). Calcium-activated chloride channel ANO1 promotes breast cancer progression by activating EGFR and CAMK signaling. Proc. Natl. Acad. Sci. U.S.A. 110, E1026–E1034. doi: 10.1073/pnas.1217072110

PubMed Abstract | CrossRef Full Text | Google Scholar

Buchholz, B., Faria, D., Schley, G., Schreiber, R., Eckardt, K. U., and Kunzelmann, K. (2014). Anoctamin 1 induces calcium-activated chloride secretion and proliferation of renal cyst-forming epithelial cells. Kidney Int. 85, 1058–1067. doi: 10.1038/ki.2013.418

PubMed Abstract | CrossRef Full Text | Google Scholar

Cabrita, I., Benedetto, R., Fonseca, A., Wanitchakool, P., Sirianant, L., Skryabin, B. V., et al. (2017). Differential effects of anoctamins on intracellular calcium signals. FASEB J. 31, 2123–2134. doi: 10.1096/fj.201600797RR

PubMed Abstract | CrossRef Full Text | Google Scholar

Caci, E., Scudieri, P., Di Carlo, E., Morelli, P., Bruno, S., De Fino, I., et al. (2015). Upregulation of TMEM16A protein in bronchial epithelial cells by bacterial pyocyanin. PLoS One 10:e0131775. doi: 10.1371/journal.pone.0131775

PubMed Abstract | CrossRef Full Text | Google Scholar

Caputo, A., Caci, E., Ferrera, L., Pedemonte, N., Barsanti, C., Sondo, E., et al. (2008). TMEM16A, a membrane protein associated with calcium-dependent chloride channel activity. Science 322, 590–594. doi: 10.1126/science.1163518

PubMed Abstract | CrossRef Full Text | Google Scholar

Chen, W., Mook, R. A. Jr., Premont, R. T., and Wang, J. (2017). Niclosamide: beyond an antihelminthic drug. Cell Signal. 41, 89–96. doi: 10.1016/j.cellsig.2017.04.001

PubMed Abstract | CrossRef Full Text | Google Scholar

Chenevert, J., Duvvuri, U., Chiosea, S., Dacic, S., Cieply, K., Kim, J., et al. (2012). DOG1: a novel marker of salivary acinar and intercalated duct differentiation. Mod. Pathol. 25, 919–929. doi: 10.1038/modpathol.2012.57

PubMed Abstract | CrossRef Full Text | Google Scholar

Cheng, T., Zhao, Y., Li, X., Lin, F., Xu, Y., Zhang, X., et al. (2007). Computation of octanol-water partition coefficients by guiding an additive model with knowledge. J. Chem. Inform. Model. 47, 2140–2148. doi: 10.1021/ci700257y

PubMed Abstract | CrossRef Full Text | Google Scholar

Cho, H., Yang, Y. D., Lee, J., Lee, B., Kim, T., Jang, Y., et al. (2012). The calcium-activated chloride channel anoctamin 1 acts as a heat sensor in nociceptive neurons. Nat. Neurosci. 15, 1015–1021. doi: 10.1038/nn.3111

PubMed Abstract | CrossRef Full Text | Google Scholar

Daina, A., Michielin, O., and Zoete, V. (2017). SwissADME: a free web tool to evaluate pharmacokinetics, drug-likeness and medicinal chemistry friendliness of small molecules. Sci. Rep. 7:42717. doi: 10.1038/srep42717

PubMed Abstract | CrossRef Full Text | Google Scholar

Dang, S., Feng, S., Tien, J., Peters, C. J., Bulkley, D., Lolicato, M., et al. (2017). Cryo-EM structures of the TMEM16A calcium-activated chloride channel. Nature 552, 426–429. doi: 10.1038/nature25024

PubMed Abstract | CrossRef Full Text | Google Scholar

Danielsson, J., Perez-Zoghbi, J., Bernstein, K., Barajas, M. B., Zhang, Y., Kumar, S., et al. (2015). Antagonists of the TMEM16A calcium-activated chloride channel modulate airway smooth muscle tone and intracellular calcium. Anesthesiology 123, 569–581. doi: 10.1097/ALN.0000000000000769

PubMed Abstract | CrossRef Full Text | Google Scholar

Danielsson, J., Yim, P., Rinderspacher, A., Fu, X. W., Zhang, Y., Landry, D. W., et al. (2014). Chloride channel blockade relaxes airway smooth muscle and potentiates relaxation by beta-agonists. Am. J. Physiol. Lung. Cell Mol. Physiol. 307, L273–L282. doi: 10.1152/ajplung.00351.2013

PubMed Abstract | CrossRef Full Text | Google Scholar

Davis, A. J., Forrest, A. S., Jepps, T. A., Valencik, M. L., Wiwchar, M., Singer, C. A., et al. (2010). Expression profile and protein translation of TMEM16A in murine smooth muscle. Am. J. Physiol. Cell Physiol. 299, C948–C959. doi: 10.1152/ajpcell.00018.2010

PubMed Abstract | CrossRef Full Text | Google Scholar

De La Fuente, R., Namkung, W., Mills, A., and Verkman, A. S. (2008). Small-molecule screen identifies inhibitors of a human intestinal calcium-activated chloride channel. Mol. Pharmacol. 73, 758–768. doi: 10.1124/mol.107.043208

PubMed Abstract | CrossRef Full Text | Google Scholar

Dunican, E. M., Elicker, B. M., Gierada, D. S., Nagle, S. K., Schiebler, M. L., Newell, J. D., et al. (2018). Mucus plugs in patients with asthma linked to eosinophilia and airflow obstruction. J. Clin. Invest. 128, 997–1009. doi: 10.1172/JCI95693

PubMed Abstract | CrossRef Full Text | Google Scholar

Duvvuri, U., Shiwarski, D. J., Xiao, D., Bertrand, C., Huang, X., Edinger, R. S., et al. (2012). TMEM16A induces MAPK and contributes directly to tumorigenesis and cancer progression. Cancer Res. 72, 3270–3281. doi: 10.1158/0008-5472.CAN-12-0475-T

PubMed Abstract | CrossRef Full Text | Google Scholar

Esteban-Gorgojo, I., Antolin-Amerigo, D., Dominguez-Ortega, J., and Quirce, S. (2018). Non-eosinophilic asthma: current perspectives. J. Asthma Allergy 11, 267–281. doi: 10.2147/JAA.S153097

PubMed Abstract | CrossRef Full Text | Google Scholar

Ferrera, L., Caputo, A., Ubby, I., Bussani, E., Zegarra-Moran, O., Ravazzolo, R., et al. (2009). Regulation of TMEM16A chloride channel properties by alternative splicing. J. Biol. Chem. 284, 33360–33368. doi: 10.1074/jbc.M109.046607

PubMed Abstract | CrossRef Full Text | Google Scholar

Forrest, A. S., Joyce, T. C., Huebner, M. L., Ayon, R. J., Wiwchar, M., Joyce, J., et al. (2012). Increased TMEM16A-encoded calcium-activated chloride channel activity is associated with pulmonary hypertension. Am. J. Physiol. Cell Physiol. 303, C1229–C1243. doi: 10.1152/ajpcell.00044.2012

PubMed Abstract | CrossRef Full Text | Google Scholar

Fricker, M., Deane, A., and Hansbro, P. M. (2014). Animal models of chronic obstructive pulmonary disease. Expert Opin Drug Discov. 9, 629–645. doi: 10.1517/17460441.2014.909805

PubMed Abstract | CrossRef Full Text | Google Scholar

Gallos, G., Remy, K. E., Danielsson, J., Funayama, H., Fu, X. W., Chang, H. Y., et al. (2013). Functional expression of the TMEM16 family of calcium-activated chloride channels in airway smooth muscle. Am. J. Physiol. Lung. Cell Mol. Physiol. 305, L625–L634. doi: 10.1152/ajplung.00068.2013

PubMed Abstract | CrossRef Full Text | Google Scholar

Gorrieri, G., Scudieri, P., Caci, E., Schiavon, M., Tomati, V., Sirci, F., et al. (2016). Goblet cell hyperplasia requires high bicarbonate transport to support mucin release. Sci. Rep. 6:36016. doi: 10.1038/srep36016

PubMed Abstract | CrossRef Full Text | Google Scholar

Grainge, C. L., Lau, L. C., Ward, J. A., Dulay, V., Lahiff, G., Wilson, S., et al. (2011). Effect of bronchoconstriction on airway remodeling in asthma. N. Engl. J. Med. 364, 2006–2015. doi: 10.1056/NEJMoa1014350

PubMed Abstract | CrossRef Full Text | Google Scholar

Hakonarson, H., Herrick, D. J., Serrano, P. G., and Grunstein, M. M. (1996). Mechanism of cytokine-induced modulation of beta-adrenoceptor responsiveness in airway smooth muscle. J. Clin. Invest. 97, 2593–2600. doi: 10.1172/JCI118708

PubMed Abstract | CrossRef Full Text | Google Scholar

Harrow, J., Frankish, A., Gonzalez, J. M., Tapanari, E., Diekhans, M., Kokocinski, F., et al. (2012). GENCODE: the reference human genome annotation for The ENCODE Project. Genome Res. 22, 1760–1774. doi: 10.1101/gr.135350.111

PubMed Abstract | CrossRef Full Text | Google Scholar

Hassan, M., Jo, T., Risse, P. A., Tolloczko, B., Lemiere, C., Olivenstein, R., et al. (2010). Airway smooth muscle remodeling is a dynamic process in severe long-standing asthma. J. Allergy Clin. Immunol. 125, 1037.e3–1045.e3. doi: 10.1016/j.jaci.2010.02.031

PubMed Abstract | CrossRef Full Text | Google Scholar

Heinze, C., Seniuk, A., Sokolov, M. V., Huebner, A. K., Klementowicz, A. E., Szijarto, I. A., et al. (2014). Disruption of vascular Ca2+-activated chloride currents lowers blood pressure. J. Clin. Invest. 124, 675–686. doi: 10.1172/JCI70025

PubMed Abstract | CrossRef Full Text | Google Scholar

Holgate, S. T., Wenzel, S., Postma, D. S., Weiss, S. T., Renz, H., and Sly, P. D. (2015). Asthma. Nat. Rev. Dis. Primers 1:15025. doi: 10.1038/nrdp.2015.25

PubMed Abstract | CrossRef Full Text | Google Scholar

Huang, F., Zhang, H., Wu, M., Yang, H., Kudo, M., Peters, C. J., et al. (2012). Calcium-activated chloride channel TMEM16A modulates mucin secretion and airway smooth muscle contraction. Proc. Natl. Acad. Sci. U.S.A. 109, 16354–16359. doi: 10.1073/pnas.1214596109

PubMed Abstract | CrossRef Full Text | Google Scholar

Ito, K., Herbert, C., Siegle, J. S., Vuppusetty, C., Hansbro, N., Thomas, P. S., et al. (2008). Steroid-resistant neutrophilic inflammation in a mouse model of an acute exacerbation of asthma. Am. J. Respir. Cell Mol. Biol. 39, 543–550. doi: 10.1165/rcmb.2008-0028OC

PubMed Abstract | CrossRef Full Text | Google Scholar

Janssen, L. J., and Sims, S. M. (1995). Ca(2+)-dependent Cl- current in canine tracheal smooth muscle cells. Am. J. Physiol. 269(1 Pt 1),C163–C169.

Google Scholar

Jin, X., Shah, S., Liu, Y., Zhang, H., Lees, M., Fu, Z., et al. (2013). Activation of the Cl channel ANO1 by localized calcium signals in nociceptive sensory neurons requires coupling with the IP3 receptor. Sci. Signal. 6:ra73. doi: 10.1126/scisignal.2004184

CrossRef Full Text | Google Scholar

Jones, B., Donovan, C., Liu, G., Gomez, H. M., Chimankar, V., Harrison, C. L., et al. (2017). Animal models of COPD: what do they tell us? Respirology 22, 21–32. doi: 10.1111/resp.12908

PubMed Abstract | CrossRef Full Text | Google Scholar

Large, W. A., and Wang, Q. (1996). Characteristics and physiological role of the Ca(2+)-activated Cl conductance in smooth muscle. Am. J. Physiol. 271(2 Pt 1), C435–C454.

Google Scholar

Laxmanan, B., and Hogarth, D. K. (2015). Bronchial thermoplasty in asthma: current perspectives. J. Asthma Allergy 8, 39–49. doi: 10.2147/JAA.S49306

PubMed Abstract | CrossRef Full Text | Google Scholar

Lee, B., Cho, H., Jung, J., Yang, Y. D., Yang, D. J., and Oh, U. (2014). Anoctamin 1 contributes to inflammatory and nerve-injury induced hypersensitivity. Mol. Pain 10:5. doi: 10.1186/1744-8069-10-5

PubMed Abstract | CrossRef Full Text | Google Scholar

Lemoine, H., and Overlack, C. (1992). Highly potent beta-2 sympathomimetics convert to less potent partial agonists as relaxants of guinea pig tracheae maximally contracted by carbachol. Comparison of relaxation with receptor binding and adenylate cyclase stimulation. J. Pharmacol. Exp. Ther. 261,258–270.

PubMed Abstract | Google Scholar

Levy, B. D., Noel, P. J., Freemer, M. M., Cloutier, M. M., Georas, S. N., Jarjour, N. N., et al. (2015). Future research directions in asthma. an nhlbi working group report. Am. J. Respir. Crit. Care Med. 192, 1366–1372. doi: 10.1164/rccm.201505-0963WS

PubMed Abstract | CrossRef Full Text | Google Scholar

Li, Y., Li, P. K., Roberts, M. J., Arend, R. C., Samant, R. S., and Buchsbaum, D. J. (2014). Multi-targeted therapy of cancer by niclosamide: a new application for an old drug. Cancer Lett. 349, 8–14. doi: 10.1016/j.canlet.2014.04.003

PubMed Abstract | CrossRef Full Text | Google Scholar

Lin, J., Jiang, Y., Li, L., Liu, Y., Tang, H., and Jiang, D. (2015). TMEM16A mediates the hypersecretion of mucus induced by Interleukin-13. Exp. Cell Res. 334, 260–269. doi: 10.1016/j.yexcr.2015.02.026

PubMed Abstract | CrossRef Full Text | Google Scholar

Liu, S., Feng, J., Luo, J., Yang, P., Brett, T. J., and Hu, H. (2016). Eact, a small molecule activator of TMEM16A, activates TRPV1 and elicits pain- and itch-related behaviours. Br. J. Pharmacol. 173, 1208–1218. doi: 10.1111/bph.13420

PubMed Abstract | CrossRef Full Text | Google Scholar

Liu, W., Lu, M., Liu, B., Huang, Y., and Wang, K. (2012). Inhibition of Ca(2+)-activated Cl(-) channel ANO1/TMEM16A expression suppresses tumor growth and invasiveness in human prostate carcinoma. Cancer Lett. 326, 41–51. doi: 10.1016/j.canlet.2012.07.015

PubMed Abstract | CrossRef Full Text | Google Scholar

Liu, X., and Farley, J. M. (1996). Acetylcholine-induced chloride current oscillations in swine tracheal smooth muscle cells. J. Pharmacol. Exp. Ther. 276, 178–186.

Google Scholar

Liu, Y., Zhang, H., Huang, D., Qi, J., Xu, J., Gao, H., et al. (2015). Characterization of the effects of Cl(-) channel modulators on TMEM16A and bestrophin-1 Ca(2)(+) activated Cl(-) channels. Pflugers Arch. 467, 1417–1430. doi: 10.1007/s00424-014-1572-5

PubMed Abstract | CrossRef Full Text | Google Scholar

Macielag, M. J., Demers, J. P., Fraga-Spano, S. A., Hlasta, D. J., Johnson, S. G., Kanojia, R. M., et al. (1998). Substituted salicylanilides as inhibitors of two-component regulatory systems in bacteria. J. Med. Chem. 41, 2939–2945. doi: 10.1021/jm9803572

PubMed Abstract | CrossRef Full Text | Google Scholar

McGarvey, L. P., Butler, C. A., Stokesberry, S., Polley, L., McQuaid, S., Abdullah, H., et al. (2014). Increased expression of bronchial epithelial transient receptor potential vanilloid 1 channels in patients with severe asthma. J. Allergy Clin. Immunol. 133, 704.e4–712.e4. doi: 10.1016/j.jaci.2013.09.016

PubMed Abstract | CrossRef Full Text | Google Scholar

Miner, K., Labitzke, K., Liu, B., Wang, P., Henckels, K., Gaida, K., et al. (2018). Drug repurposing: the anthelmintics niclosamide and nitazoxanide are potent TMEM16A antagonists that fully bronchodilate airways. bioRxiv (Preprint) doi: 10.1101/254888

CrossRef Full Text | Google Scholar

Mook, R. A. Jr., Wang, J., Ren, X. R., Chen, M., Spasojevic, I., Barak, L. S., et al. (2015). Structure-activity studies of Wnt/beta-catenin inhibition in the Niclosamide chemotype: identification of derivatives with improved drug exposure. Bioorg. Med. Chem. 23, 5829–5838. doi: 10.1016/j.bmc.2015.07.001

PubMed Abstract | CrossRef Full Text | Google Scholar

Moore, P. E., Lahiri, T., Laporte, J. D., Church, T., Panettieri, R. A. Jr., and Shore, S. A. (2001). Selected contribution: synergism between TNF-alpha and IL-1 beta in airway smooth muscle cells: implications for beta-adrenergic responsiveness. J. Appl. Physiol. 91, 1467–1474.

PubMed Abstract | Google Scholar

Moore, W. C., Bleecker, E. R., Curran-Everett, D., Erzurum, S. C., Ameredes, B. T., Bacharier, L., et al. (2007). Characterization of the severe asthma phenotype by the national heart, lung, and blood institute’s severe asthma research program. J. Allergy Clin. Immunol. 119, 405–413. doi: 10.1016/j.jaci.2006.11.639

PubMed Abstract | CrossRef Full Text | Google Scholar

Moore, W. C., Hastie, A. T., Li, X., Li, H., Busse, W. W., Jarjour, N. N., et al. (2014). Sputum neutrophil counts are associated with more severe asthma phenotypes using cluster analysis. J. Allergy Clin. Immunol. 133, 1557.e5–1563.e5. doi: 10.1016/j.jaci.2013.10.011

PubMed Abstract | CrossRef Full Text | Google Scholar

Moore, W. C., Meyers, D. A., Wenzel, S. E., Teague, W. G., Li, H., Li, X., et al. (2010). Identification of asthma phenotypes using cluster analysis in the Severe Asthma Research Program. Am. J. Respir. Crit. Care Med. 181, 315–323. doi: 10.1164/rccm.200906-0896OC

PubMed Abstract | CrossRef Full Text | Google Scholar

Mortazavi, A., Williams, B. A., McCue, K., Schaeffer, L., and Wold, B. (2008). Mapping and quantifying mammalian transcriptomes by RNA-Seq. Nat. Methods 5, 621–628. doi: 10.1038/nmeth.1226

PubMed Abstract | CrossRef Full Text | Google Scholar

Nakano, T., Inoue, H., Fukuyama, S., Matsumoto, K., Matsumura, M., Tsuda, M., et al. (2006). Niflumic acid suppresses interleukin-13-induced asthma phenotypes. Am. J. Respir. Crit. Care Med. 173, 1216–1221. doi: 10.1164/rccm.200410-1420OC

PubMed Abstract | CrossRef Full Text | Google Scholar

Namkung, W., Thiagarajah, J. R., Phuan, P. W., and Verkman, A. S. (2010). Inhibition of Ca2+-activated Cl- channels by gallotannins as a possible molecular basis for health benefits of red wine and green tea. FASEB J. 24, 4178–4186. doi: 10.1096/fj.10-160648

PubMed Abstract | CrossRef Full Text | Google Scholar

Namkung, W., Yao, Z., Finkbeiner, W. E., and Verkman, A. S. (2011). Small-molecule activators of TMEM16A, a calcium-activated chloride channel, stimulate epithelial chloride secretion and intestinal contraction. FASEB J. 25, 4048–4062. doi: 10.1096/fj.11-191627

PubMed Abstract | CrossRef Full Text | Google Scholar

Ni, Y. L., Kuan, A. S., and Chen, T. Y. (2014). Activation and inhibition of TMEM16A calcium-activated chloride channels. PLoS One 9:e86734. doi: 10.1371/journal.pone.0086734

PubMed Abstract | CrossRef Full Text | Google Scholar

Osada, T., Chen, M., Yang, X. Y., Spasojevic, I., Vandeusen, J. B., Hsu, D., et al. (2011). Antihelminth compound niclosamide downregulates Wnt signaling and elicits antitumor responses in tumors with activating APC mutations. Cancer Res. 71, 4172–4182. doi: 10.1158/0008-5472.CAN-10-3978

PubMed Abstract | CrossRef Full Text | Google Scholar

Ousingsawat, J., Mirza, M., Tian, Y., Roussa, E., Schreiber, R., Cook, D. I., et al. (2011). Rotavirus toxin NSP4 induces diarrhea by activation of TMEM16A and inhibition of Na+ absorption. Pflugers Arch. 461, 579–589. doi: 10.1007/s00424-011-0947-0

PubMed Abstract | CrossRef Full Text | Google Scholar

Patterson, S., and Wyllie, S. (2014). Nitro drugs for the treatment of trypanosomatid diseases: past, present, and future prospects. Trends Parasitol. 30, 289–298. doi: 10.1016/j.pt.2014.04.003

PubMed Abstract | CrossRef Full Text | Google Scholar

Paulino, C., Kalienkova, V., Lam, A. K. M., Neldner, Y., and Dutzler, R. (2017a). Activation mechanism of the calcium-activated chloride channel TMEM16A revealed by cryo-EM. Nature 552, 421–425. doi: 10.1038/nature24652

PubMed Abstract | CrossRef Full Text | Google Scholar

Paulino, C., Neldner, Y., Lam, A. K., Kalienkova, V., Brunner, J. D., Schenck, S., et al. (2017b). Structural basis for anion conduction in the calcium-activated chloride channel TMEM16A. Elife 6:e26232. doi: 10.7554/eLife.26232

PubMed Abstract | CrossRef Full Text | Google Scholar

Perez-Cornejo, P., Gokhale, A., Duran, C., Cui, Y., Xiao, Q., Hartzell, H. C., et al. (2012). Anoctamin 1 (Tmem16A) Ca2+-activated chloride channel stoichiometrically interacts with an ezrin-radixin-moesin network. Proc. Natl. Acad. Sci. U.S.A. 109, 10376–10381. doi: 10.1073/pnas.1200174109

PubMed Abstract | CrossRef Full Text | Google Scholar

Peters, C. J., Yu, H., Tien, J., Jan, Y. N., Li, M., and Jan, L. Y. (2015). Four basic residues critical for the ion selectivity and pore blocker sensitivity of TMEM16A calcium-activated chloride channels. Proc. Natl. Acad. Sci. U.S.A. 112, 3547–3552. doi: 10.1073/pnas.1502291112

PubMed Abstract | CrossRef Full Text | Google Scholar

Pineda-Farias, J. B., Barragan-Iglesias, P., Loeza-Alcocer, E., Torres-Lopez, J. E., Rocha-Gonzalez, H. I., Perez-Severiano, F., et al. (2015). Role of anoctamin-1 and bestrophin-1 in spinal nerve ligation-induced neuropathic pain in rats. Mol. Pain 11:41. doi: 10.1186/s12990-015-0042-1

PubMed Abstract | CrossRef Full Text | Google Scholar

Piper, A. S., and Greenwood, I. A. (2003). Anomalous effect of anthracene-9-carboxylic acid on calcium-activated chloride currents in rabbit pulmonary artery smooth muscle cells. Br. J. Pharmacol. 138, 31–38. doi: 10.1038/sj.bjp.0705000

PubMed Abstract | CrossRef Full Text | Google Scholar

Piper, A. S., Greenwood, I. A., and Large, W. A. (2002). Dual effect of blocking agents on Ca2+-activated Cl(-) currents in rabbit pulmonary artery smooth muscle cells. J. Physiol. 539(Pt 1), 119–131.

PubMed Abstract | Google Scholar

Qin, Y., Jiang, Y., Sheikh, A. S., Shen, S., Liu, J., and Jiang, D. (2016). Interleukin-13 stimulates MUC5AC expression via a STAT6-TMEM16A-ERK1/2 pathway in human airway epithelial cells. Int. Immunopharmacol. 40, 106–114. doi: 10.1016/j.intimp.2016.08.033

PubMed Abstract | CrossRef Full Text | Google Scholar

Rapetti-Mauss, R., Bustos, V., Thomas, W., McBryan, J., Harvey, H., Lajczak, N., et al. (2017). Bidirectional KCNQ1:beta-catenin interaction drives colorectal cancer cell differentiation. Proc. Natl. Acad. Sci. U.S.A. 114, 4159–4164. doi: 10.1073/pnas.1702913114

PubMed Abstract | CrossRef Full Text | Google Scholar

Robinett, K. S., White, C. J., Akoluk, A., An, S. S., Panettieri, R. A. Jr., and Liggett, S. B. (2014). Bitter taste receptor function in asthmatic and nonasthmatic human airway smooth muscle cells. Am. J. Respir. Cell Mol. Biol. 50, 678–683. doi: 10.1165/rcmb.2013-0439RC

PubMed Abstract | CrossRef Full Text | Google Scholar

Robinson, M. D., and Oshlack, A. (2010). A scaling normalization method for differential expression analysis of RNA-seq data. Genome Biol. 11:R25. doi: 10.1186/gb-2010-11-3-r25

PubMed Abstract | CrossRef Full Text | Google Scholar

Romanenko, V. G., Catalan, M. A., Brown, D. A., Putzier, I., Hartzell, H. C., Marmorstein, A. D., et al. (2010). Tmem16A encodes the Ca2+-activated Cl- channel in mouse submandibular salivary gland acinar cells. J. Biol. Chem. 285, 12990–13001. doi: 10.1074/jbc.M109.068544

PubMed Abstract | CrossRef Full Text | Google Scholar

Rossignol, J. F. (2014). Nitazoxanide: a first-in-class broad-spectrum antiviral agent. Antiviral Res. 110, 94–103. doi: 10.1016/j.antiviral.2014.07.014

PubMed Abstract | CrossRef Full Text | Google Scholar

Rossignol, J. F., Abu-Zekry, M., Hussein, A., and Santoro, M. G. (2006). Effect of nitazoxanide for treatment of severe rotavirus diarrhoea: randomised double-blind placebo-controlled trial. Lancet 368, 124–129. doi: 10.1016/S0140-6736(06)68852-1

PubMed Abstract | CrossRef Full Text | Google Scholar

Schreiber, R., Ousingsawat, J., Wanitchakool, P., Sirianant, L., Benedetto, R., Reiss, K., et al. (2018). Regulation of TMEM16A/ANO1 and TMEM16F/ANO6 ion currents and phospholipid scrambling by Ca(2+) and plasma membrane lipid. J. Physiol. 596, 217–229. doi: 10.1113/JP275175

PubMed Abstract | CrossRef Full Text | Google Scholar

Schroeder, B. C., Cheng, T., Jan, Y. N., and Jan, L. Y. (2008). Expression cloning of TMEM16A as a calcium-activated chloride channel subunit. Cell 134,1019–1029. doi: 10.1016/j.cell.2008.09.003

PubMed Abstract | CrossRef Full Text | Google Scholar

Scudieri, P., Caci, E., Bruno, S., Ferrera, L., Schiavon, M., Sondo, E., et al. (2012). Association of TMEM16A chloride channel overexpression with airway goblet cell metaplasia. J. Physiol. 590, 6141–6155. doi: 10.1113/jphysiol.2012.240838

PubMed Abstract | CrossRef Full Text | Google Scholar

Senkowski, W., Zhang, X., Olofsson, M. H., Isacson, R., Hoglund, U., Gustafsson, M., et al. (2015). Three-dimensional cell culture-based screening identifies the anthelmintic drug nitazoxanide as a candidate for treatment of colorectal cancer. Mol. Cancer Ther. 14, 1504–1516. doi: 10.1158/1535-7163.MCT-14-0792

PubMed Abstract | CrossRef Full Text | Google Scholar

Shore, S. A., Laporte, J., Hall, I. P., Hardy, E., and Panettieri, R. A. Jr. (1997). Effect of IL-1 beta on responses of cultured human airway smooth muscle cells to bronchodilator agonists. Am. J. Respir. Cell Mol. Biol. 16, 702–712. doi: 10.1165/ajrcmb.16.6.9191472

PubMed Abstract | CrossRef Full Text | Google Scholar

Sondo, E., Caci, E., and Galietta, L. J. (2014). The TMEM16A chloride channel as an alternative therapeutic target in cystic fibrosis. Int. J. Biochem. Cell Biol. 52, 73–76. doi: 10.1016/j.biocel.2014.03.022

PubMed Abstract | CrossRef Full Text | Google Scholar

Stocks, M. J., Alcaraz, L., Bailey, A., Bonnert, R., Cadogan, E., Christie, J., et al. (2014). Discovery of AZD3199, an inhaled ultralong acting beta2 receptor agonist with rapid onset of action. ACS Med. Chem. Lett. 5, 416–421. doi: 10.1021/ml4005232

PubMed Abstract | CrossRef Full Text | Google Scholar

Strong, P., Ito, K., Murray, J., and Rapeport, G. (2018). Current approaches to the discovery of novel inhaled medicines. Drug Discov. Today 23, 1705–1717. doi: 10.1016/j.drudis.2018.05.017

PubMed Abstract | CrossRef Full Text | Google Scholar

Takayama, Y., Uta, D., Furue, H., and Tominaga, M. (2015). Pain-enhancing mechanism through interaction between TRPV1 and anoctamin 1 in sensory neurons. Proc. Natl. Acad. Sci. U.S.A. 112, 5213–5218. doi: 10.1073/pnas.1421507112

PubMed Abstract | CrossRef Full Text | Google Scholar

Thiagarajah, J. R., Donowitz, M., and Verkman, A. S. (2015). Secretory diarrhoea: mechanisms and emerging therapies. Nat. Rev. Gastroenterol. Hepatol. 12, 446–457. doi: 10.1038/nrgastro.2015.111

PubMed Abstract | CrossRef Full Text | Google Scholar

Tian, Y., Schreiber, R., and Kunzelmann, K. (2012). Anoctamins are a family of Ca2+-activated Cl- channels. J. Cell Sci. 125(Pt 21), 4991–4998. doi: 10.1242/jcs.109553

PubMed Abstract | CrossRef Full Text | Google Scholar

Trian, T., Burgess, J. K., Niimi, K., Moir, L. M., Ge, Q., Berger, P., et al. (2011). beta2-Agonist induced cAMP is decreased in asthmatic airway smooth muscle due to increased PDE4D. PLoS One 6:e20000. doi: 10.1371/journal.pone.0020000

PubMed Abstract | CrossRef Full Text | Google Scholar

Tschumperlin, D. J. (2011). Physical forces and airway remodeling in asthma. N. Engl. J. Med. 364, 2058–2059. doi: 10.1056/NEJMe1103121

PubMed Abstract | CrossRef Full Text | Google Scholar

Wang, H., Zou, L., Ma, K., Yu, J., Wu, H., Wei, M., et al. (2017). Cell-specific mechanisms of TMEM16A Ca2+-activated chloride channel in cancer. Mol. Cancer 16:152. doi: 10.1186/s12943-017-0720-x

PubMed Abstract | CrossRef Full Text | Google Scholar

Wang, P., Zhao, W., Sun, J., Tao, T., Chen, X., Zheng, Y. Y., et al. (2018). Inflammatory mediators mediate airway smooth muscle contraction through a G protein-coupled receptor-transmembrane protein 16A-voltage-dependent Ca(2+) channel axis and contribute to bronchial hyperresponsiveness in asthma. J. Allergy Clin. Immunol. 141, 1259.11–1268.e11. doi: 10.1016/j.jaci.2017.05.053

PubMed Abstract | CrossRef Full Text | Google Scholar

Wang, Q., Leo, M. D., Narayanan, D., Kuruvilla, K. P., and Jaggar, J. H. (2016). Local coupling of TRPC6 to ANO1/TMEM16A channels in smooth muscle cells amplifies vasoconstriction in cerebral arteries. Am. J. Physiol. Cell Physiol. 310, C1001–C1009. doi: 10.1152/ajpcell.00092.2016

PubMed Abstract | CrossRef Full Text | Google Scholar

Wang, Y. X., and Kotlikoff, M. I. (1997). Inactivation of calcium-activated chloride channels in smooth muscle by calcium/calmodulin-dependent protein kinase. Proc. Natl. Acad. Sci. U.S.A. 94, 14918–14923.

PubMed Abstract | Google Scholar

Wellman, G. C., and Nelson, M. T. (2003). Signaling between SR and plasmalemma in smooth muscle: sparks and the activation of Ca2+-sensitive ion channels. Cell Calcium 34, 211–229.

Google Scholar

Wenzel, S. E. (2012). Asthma phenotypes: the evolution from clinical to molecular approaches. Nat. Med. 18, 716–725. doi: 10.1038/nm.2678

PubMed Abstract | CrossRef Full Text | Google Scholar

White, C. A. Jr. (2004). Nitazoxanide: a new broad spectrum antiparasitic agent. Expert Rev. Anti Infect. Ther. 2, 43–49.

Google Scholar

Yang, Y. D., Cho, H., Koo, J. Y., Tak, M. H., Cho, Y., Shim, W. S., et al. (2008). TMEM16A confers receptor-activated calcium-dependent chloride conductance. Nature 455, 1210–1215. doi: 10.1038/nature07313

PubMed Abstract | CrossRef Full Text | Google Scholar

Yim, P. D., Gallos, G., Perez-Zoghbi, J. F., Trice, J., Zhang, Y., Siviski, M., et al. (2013). Chloride channel blockers promote relaxation of TEA-induced contraction in airway smooth muscle. J. Smooth Muscle Res. 49,112–124.

PubMed Abstract | Google Scholar

Yocum, G. T., Chen, J., Choi, C. H., Townsend, E. A., Zhang, Y., Xu, D., et al. (2017). Role of transient receptor potential vanilloid 1 in the modulation of airway smooth muscle tone and calcium handling. Am. J. Physiol. Lung. Cell Mol. Physiol. 312, L812–L821. doi: 10.1152/ajplung.00064.2017

PubMed Abstract | CrossRef Full Text | Google Scholar

Zhang, C. H., Li, Y., Zhao, W., Lifshitz, L. M., Li, H., Harfe, B. D., et al. (2013). The transmembrane protein 16A Ca(2+)-activated Cl- channel in airway smooth muscle contributes to airway hyperresponsiveness. Am. J. Respir. Crit. Care Med. 187, 374–381. doi: 10.1164/rccm.201207-1303OC

PubMed Abstract | CrossRef Full Text | Google Scholar

Zhang, Y., Wang, X., Wang, H., Jiao, J., Li, Y., Fan, E., et al. (2015). TMEM16A-mediated mucin secretion in IL-13-induced nasal epithelial cells from chronic rhinosinusitis patients. Allergy Asthma Immunol. Res. 7, 367–375. doi: 10.4168/aair.2015.7.4.367

PubMed Abstract | CrossRef Full Text | Google Scholar

Keywords: TMEM16A antagonist, niclosamide, nitazoxanide, bronchodilator, desensitization, airway smooth muscle (ASM), calcium-activated chloride channel, drug repositioning

Citation: Miner K, Labitzke K, Liu B, Wang P, Henckels K, Gaida K, Elliott R, Chen JJ, Liu L, Leith A, Trueblood E, Hensley K, Xia X-Z, Homann O, Bennett B, Fiorino M, Whoriskey J, Yu G, Escobar S, Wong M, Born TL, Budelsky A, Comeau M, Smith D, Phillips J, Johnston JA, McGivern JG, Weikl K, Powers D, Kunzelmann K, Mohn D, Hochheimer A and Sullivan JK (2019) Drug Repurposing: The Anthelmintics Niclosamide and Nitazoxanide Are Potent TMEM16A Antagonists That Fully Bronchodilate Airways. Front. Pharmacol. 10:51. doi: 10.3389/fphar.2019.00051

Received: 23 November 2018; Accepted: 18 January 2019;
Published: 14 February 2019.

Edited by:

Heike Wulff, University of California, Davis, United States

Reviewed by:

George Chandy, University of California, Irvine, United States
Peter Bradding, University of Leicester, United Kingdom

Copyright © 2019 Miner, Labitzke, Liu, Wang, Henckels, Gaida, Elliott, Chen, Liu, Leith, Trueblood, Hensley, Xia, Homann, Bennett, Fiorino, Whoriskey, Yu, Escobar, Wong, Born, Budelsky, Comeau, Smith, Phillips, Johnston, McGivern, Weikl, Powers, Kunzelmann, Mohn, Hochheimer and Sullivan. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.

*Correspondence: Deanna Mohn, dmohn@amgen.com orcid.org/0000-0001-5748-7591 Andreas Hochheimer, andreas.hochheimer@isarbioscience.de orcid.org/0000-0002-1572-9334 John K. Sullivan, jksullivan@polestarbio.com orcid.org/0000-0001-9840-852X

Present address: Katja Labitzke, Clariant Produkte (Deutschland) GmbH, Planegg, Germany, Longbin Liu, CHDI Management, CHDI Foundation, Los Angeles, CA, United States, Anh Leith, Department of Genome Sciences, University of Washington, Seattle, WA, United States, Esther Trueblood, Seattle Genetics, Bothell, WA, United States, Kelly Hensley, Seattle Genetics, Bothell, WA, United States, Teresa L. Born, Sartorius Stedim BioOutsource, Cambridge, MA, United States, Alison Budelsky, Immunology Research, Lilly Research Laboratories, San Diego, CA, United States, Dirk Smith, Insight BioConsulting, Bainbridge Island, WA, United States, Kerstin Weikl, Assay.Works GmbH, Regensburg, Germany, Andreas Hochheimer, ISAR Bioscience GmbH, Planegg, Germany, John K. Sullivan, PolestarBio LLC, Newbury Park, CA, United States