Cystic fibrosis (CF) is a channelopathy caused by loss-of-function mutations of the cystic fibrosis transmembrane conductance regulator (CFTR) gene, which encodes a phosphorylation-activated and adenosine triphosphate (ATP)–gated chloride channel. In the past few years, high-throughput drug screening has successfully realized the first US Food and Drug Administration–approved therapy for CF, called ivacaftor (or VX-770). A more recent CFTR potentiator, GLPG1837 (N-(3-carbamoyl-5,5,7,7-tetramethyl-4,7-dihydro-5H-thieno[2,3-c]pyran-2-yl)-1H-pyrazole-3-carboxamide), has been shown to exhibit a higher efficacy than ivacaftor for the G551D mutation, yet the underlying mechanism of GLPG1837 remains unclear. Here we find that despite their differences in potency and efficacy, GLPG1837 and VX-770 potentiate CFTR gating in a remarkably similar manner. Specifically, they share similar effects on single-channel kinetics of wild-type CFTR. Their actions are independent of nucleotide-binding domain (NBD) dimerization and ATP hydrolysis, critical steps controlling CFTR’s gate opening and closing, respectively. By applying the two reagents together, we provide evidence that GLPG1837 and VX-770 likely compete for the same site, whereas GLPG1837 and the high-affinity ATP analogue 2′-deoxy-N6-(2-phenylethyl)-adenosine-5′-O-triphosphate (dPATP) work synergistically through two different sites. We also find that the apparent affinity for GLPG1837 is dependent on the open probability of the channel, suggesting a state-dependent binding of the drug to CFTR (higher binding affinity for the open state than the closed state), which is consistent with the classic mechanism for allosteric modulation. We propose a simple four-state kinetic model featuring an energetic coupling between CFTR gating and potentiator binding to explain our experimental results.
Introduction
CFTR, a member of the ATP-binding cassette (ABC) protein superfamily, is unique in that instead of being a transporter, CFTR is a chloride channel that plays a critical role in the regulation of water and salt movement across epithelium-lining tissues (Riordan et al., 1989; Quinton and Reddy, 1991; Bear et al., 1992). Loss-of-function mutations in the CFTR gene are the root cause of the lethal hereditary disease cystic fibrosis (CF), affecting one in every ∼2,500 newborns of Caucasian origin (Zielenski and Tsui, 1995; Rowe et al., 2005). Recently, the breakthrough in solving the atomic structure of human CFTR protein provided the first glimpse of the fundamental features of this channel at a molecular level (Liu et al., 2017). Like most ABC transporters, CFTR inherits two transmembrane domains (TMDs) that craft the chloride permeation pathway and two cytosolic nucleotide-binding domains (NBD1 and NBD2), where CFTR’s ligand ATP binds to gate the channel. In addition, CFTR possesses a unique regulatory (R) domain containing multiple serine residues for PKA-dependent phosphorylation (Ostedgaard et al., 2001). Once the R domain is phosphorylated, the gating cycle of CFTR is driven by ATP binding–induced NBD dimerization and hydrolysis-elicited partial separation of the NBD dimer (Vergani et al., 2003, 2005; Hwang and Sheppard, 2009; Tsai et al., 2009, 2010; Chaves and Gadsby, 2015).
Approximately 2,000 different mutations have been identified (http://www.genet.sickkids.on.ca/cftr/) since the first cloning of the CFTR gene in 1989 (Riordan et al., 1989). Currently, disease-associated mutations are grouped into six classes based on their defective mechanisms (Welsh and Smith, 1993; Zielenski and Tsui, 1995; Rowe et al., 2005; Veit et al., 2016). The most prevalent mutation, the deletion of phenylalanine at position 508 (ΔF508), results in both severe folding (Du et al., 2005; Thibodeau et al., 2005; Cui et al., 2007; Rosser et al., 2008; Du and Lukacs, 2009) and gating (Dalemans et al., 1991; Hwang et al., 1997; Cui et al., 2006) defects. On the other hand, the third most common disease-associated mutation, G551D, primarily impairs channel gating, reducing the open probability (Po) by more than 100-fold (Welsh and Smith, 1993; Cai et al., 2006; Bompadre et al., 2007). In the past few decades, tremendous efforts have been made to develop compounds targeting the CFTR protein itself. Those CFTR modulators that can increase the Po of CFTR are known as CFTR potentiators (Hwang and Sheppard, 1999; Rowe and Verkman, 2013; Barry et al., 2015; Jih et al., 2017). Among the plethora of CFTR potentiators, however, only a few have been extensively investigated. Based on their mechanisms of action, these well-studied reagents can be classified into three major groups: (1) ATP analogues that target the ATP-binding sites, such as N6-(2-phenylethyl)-ATP (PATP; Zhou et al., 2005; Bompadre et al., 2008; Tsai et al., 2010), 2′-deoxy-ATP (dATP; Aleksandrov et al., 2002; Cai et al., 2006), and 2′-deoxy-N6-(2-phenylethyl)-adenosine-5′-O-triphosphate (dPATP; Miki et al., 2010); (2) compounds that are proposed to work on the TMDs to stabilize the open channel configuration, including nitrate (Yeh et al., 2015) and N-(2,4-di-tert-butyl-5-hydroxyphenyl)-4-oxo-1,4-dihydroquinoline-3-carboxamide (VX-770; Fig. 1 A; Van Goor et al., 2009; Eckford et al., 2012; Jih and Hwang, 2013); and (3) 5-nitro-2-(3-phenylpropylamino)benzoate (NPPB; Fig. 1 B; Wang et al., 2005), which appears to promote gate opening with unsettled mechanisms (Csanády and Töröcsik, 2014; Lin et al., 2016).
The FDA’s approval of VX-770 for the treatment of CF in 2012 ushered in a new era of personalized medicine in CF. Although the original application of VX-770 was limited to patients carrying the G551D mutation, now it is used for a broad spectrum of mutations because in vitro studies showed that VX-770 potentiates the activity of a variety of CFTR mutants with gating abnormalities (Yu et al., 2012; Van Goor et al., 2014). One potential downside with VX-770 as a therapeutic reagent is that prolonged exposure of VX-770 may have negative impacts on the action of VX-809 (Cholon et al., 2014; Veit et al., 2014), a CFTR corrector designed to ameliorate the trafficking defect associated with the ΔF508 mutation (Van Goor et al., 2011). In addition, the VX-770–treated G551D channels have a Po still less than 10% of that of WT channels (Jih and Hwang, 2013; but compare Van Goor et al., 2009). Therefore, developing new CFTR potentiators with improved pharmacological properties remains an important goal.
Lately, a novel and potentially more effective CFTR potentiator, N-(3-carbamoyl-5,5,7,7-tetramethyl-4,7-dihydro-5H-thieno[2,3-c]pyran-2-yl)-1H-pyrazole-3-carboxamide (GLPG1837; Fig. 1 C), was identified by Conrath et al. (2016. The 30th Annual North American Cystic Fibrosis Conference. Poster 23) using a YFP-based high-throughput screening assay. In the current study, we demonstrate that although it has a chemical structure different from VX-770 (compare Fig. 1, A and C), GLPG1837 may share a common mechanism with VX-770 in modulating CFTR functions by binding to the same site. Our data also show a Po-dependent shift of the dose–response relationships for GLPG1837—the lower the Po, the higher the concentration required for reaching 50% of the maximal effect (K1/2), a phenomenon that can be explained by state-dependent binding. We therefore propose a classic allosteric model featuring an energetic coupling between binding of GLPG1837 and gating. Computer simulations using this simple model can indeed replicate most of our data.
Materials and methods
Cell culture and transfection
Chinese hamster ovary (CHO) cells were used for all patch-clamp experiments. CHO cells were grown at 37°C in Dulbecco’s modified Eagle’s medium supplemented with 10% (vol/vol) FBS. 1 d before the transfection, cells were trypsinized and transferred to 35-mm tissue culture dishes. We performed transfection with pcDNA plasmids carrying different CFTR constructs and green fluorescence protein encoding pEGFP-C3 (Takara Bio) by using PolyFect transfection reagent (Qiagen). After transfection, cells were incubated at 27°C for 2–3 d before experiments for microscopic current recordings and 3–6 d for macroscopic current recordings.
Mutagenesis
CFTR mutants were constructed using QuikChange XL kit (Agilent) according to the manufacture’s protocol. All DNA constructs were sequenced by the DNA Core Facility (University of Missouri, Columbia, MO) to confirm the mutation identity.
Electrophysiological recordings
In patch-clamp experiments, glass chips carrying the transfected cells were transferred to a chamber located on the stage of an inverted microscope (IX51; Olympus). A two-stage micropipette puller (PP-81; Narishige) was used to pull borosilicate capillary glasses into patch pipettes, which were then polished with a homemade microforge to a resistance of 2–4 MΩ in the bath solution. Membrane patches were excised to an inside-out configuration once a seal resistance >40 GΩ was reached. Subsequently, the pipette tip was placed to the outlet of a three-barrel perfusion system and perfused with 25 IU PKA and 2 mM ATP until the current reached a steady state. To maintain the phosphorylation level of the CFTR channels, 6 IU PKA was added to all other ATP-containing solutions applied thereafter. All electrophysiological data were recorded with a patch-clamp amplifier (EPC9; HEKA) at room temperature. An eight-pole Bessel filter (LPF-8; Warner Instruments) was used to filter the data at 100 Hz and digitized to a computer at a sampling rate of 500 Hz. For macroscopic and single-channel recordings, the membrane potential was held at −30 and −50 mV, respectively. Solution changes were effected with a fast solution change system (SF-77B; Warner Instruments) with a dead time of ∼30 ms (Tsai et al., 2009).
Of note, the effect of VX-770 cannot be washed out by a continuous perfusion of VX-770–free solution within the experimentally permissible time. As a result, we could not bracket our experiments with this reagent. Instead, all experiments with VX-770 were done after first obtaining a control in the same patch. All devices that had been in contact with VX-770 were repeatedly washed with 50% DMSO to minimize contamination by residual VX-770.
Chemicals and solution compositions
For all of the patch-clamp experiments, pipette solution contained (mM) 140 NMDG-Cl, 5 CaCl2, 2 MgCl2, and 10 HEPES, pH 7.4 with NMDG. Cells were perfused with a bath solution containing (mM) 145 NaCl, 2 MgCl2, 5 KCl, 1 CaCl2, 5 glucose, 5 HEPES, and 20 sucrose, pH 7.4 with NaOH, before patch excision. After an inside-out configuration was established, the patch was perfused with a standard perfusion solution containing (mM) 150 NMDG-Cl, 2 MgCl2, 1 CaCl2, 10 EGTA, and 8 Tris, pH 7.4 with NMDG.
PKA and MgATP were purchased from Sigma-Aldrich. MgATP was stored in 500 mM stock solutions at −20°C. The [MgATP] used in this study was 2 mM unless otherwise indicated in the figures. dPATP was custom-synthesized by Biolog Life Science Institute and stored in 10-mM stock solutions at −70°C. VX-770, provided by R. Bridges (Rosalind Franklin University, North Chicago, IL), was stored as a 100-µM stock in DMSO at −70°C. GLPG1837 was provided by Galapagos and stored as a 10-mM stock at −20°C. All chemicals were diluted to the concentrations indicated in each figure with the perfusion solution, and the pH was adjusted to 7.4 with NMDG.
Data analysis and statistics
Igor Pro program (Wave-Metrics) was used to measure the steady-state mean current amplitude and estimate the relaxation time constant. The current decay upon removal of ATP, dPATP, or GLPG1837 were fitted with built-in single or double exponential functions to obtain the relaxation time constants. Single-channel kinetic analysis was done with a program developed by Csanády (2000), which allows microscopic kinetic analysis for recordings containing fewer than eight openings. We took a more conservative approach to analyze only traces with less than four opening steps for all the microscopic experiments. The partition coefficients and logP values were calculated by using ACD/ChemSketch (www.acdlabs.com). Student’s t tests assuming equal variance were conducted with Excel for the comparisons showing statistical probability. Paired t test was conducted for analysis in Figs. 4 B and 7 D and Fig. S3 B. P < 0.05 was considered statistically significant. The error bars represent SEM in all figures.
Online supplemental material
Fig. S1 shows raw traces of locked-open events occasionally observed with WT-CFTR potentiated by GLPG1837. Figs. S2 and S3 present additional enhancement of G551D-CFTR currents by nitrate in the presence of dPATP plus GLPG1837. Fig. S4 demonstrates the effect of GLPG1837 on WT-CFTR in the absence of ATP. Figs. S5 and S6 are simulated results for data shown in Figs. 2 and 7, respectively. Table S1 lists the parameters used in the simulations. The online supplemental material also includes discussion on the limitations and possible alternatives for the allosteric modulation model proposed in Fig. 10.
Results
As a first step to characterize GLPG1837, we examined the effect of acute addition of this compound to macroscopic WT-CFTR currents preactivated with PKA and ATP to a steady state at a holding potential of −30 mV in excised inside-out patches. Fig. 2 A shows a representative real-time recording in which an application of 3 µM GLPG1837 in the continuous presence of ATP enhances the currents by 2.06 ± 0.08-fold (n = 19). This potentiation effect is reversible and concentration dependent. The dose responses at different concentrations of GLPG1837 were normalized to the currents at 3 µM GLPG1837 in the same patch, and the data can be fitted with the Hill equation, yielding a K1/2 of 0.23 ± 0.12 µM and a Hill coefficient of 0.70 ± 0.24 (Fig. 2 B).
Interestingly, when we closely examined the macroscopic current relaxations upon removal of ATP in the absence or presence of GLPG1837, a striking difference emerged. The current decline after washout of ATP in the absence of GLPG1837 can be well fitted with a double-exponential function, with 90% of the current decay within the first fast phase (Fig. 2 A, red line). However, in the presence of GLPG1837, removal of ATP shows a visually distinguishable slow phase in addition to the fast phase (Fig. 2 A, blue line). In addition, the residual currents in the continuous presence of GLPG1837 remain discernable for hundreds of seconds even in a solution without ATP (Fig. 2 A, boxed area). Fig. 2 C magnifies the boxed area in Fig. 2 A, highlighting that the residual currents long after ATP removal can be reversibly potentiated by GLPG1837. Fig. 2 (D and E) summarizes the parameters of double-exponential curve fitting for the current relaxations upon ATP washout. The fraction of second phase in the presence of GLPG1837 is significantly larger than that in the absence of GLPG1837 (22 ± 3% and 11 ± 2%, respectively; n = 7; P < 0.05). Of note, because the ATP-induced currents in the presence of GLPG1837 are approximately twofold larger than the control, the twofold increase of the second fraction by GLPG1837 indicates that the initial current after the completion of the fast-phase current drop is increased by fourfold with GLPG1837. In addition, the relaxation time constant of the fast current decay (τ1) is prolonged by GLPG1837 (τ1 = 1.51 ± 0.36 s with GLPG1837; τ1 = 0.62 ± 0.08 s without GLPG1837; n = 7; P < 0.05), whereas the time constants of the second phase (τ2) are not significantly different. These observations led us to conclude that GLPG1837 potentiates both ATP-dependent and ATP-independent gating of WT-CFTR.
We next tested the effects of GLPG1837 on WT-CFTR in patches containing fewer than four channels to understand the kinetic effects of GLPG1837. As shown in Fig. 3 A, single-channel activity of WT-CFTR was recorded in the same patch either with or without the application of GLPG1837. The open time is indeed visibly longer in the GLPG1837 than in its absence. Fig. 3 B summarizes results from single-channel kinetic analysis and confirms an increase of the Po by GLPG1837 (0.71 ± 0.01 and 0.44 ± 0.02 with and without GLPG1837, n = 9 and n = 7, respectively; P < 0.001). Although the mean closed time constant (τc) is slightly decreased, GLPG1837 predominantly increases the Po by prolonging the open time constant (τo) from 397 ± 56 ms (n = 7) to 851 ± 55 ms (n = 9; P < 0.001). Comparing the macroscopic and microscopic data, we noticed a discrepancy between the predicted Po estimated by macroscopic current fold increase (∼0.9) and the actual Po measured from microscopic kinetic analysis (0.71). Fig. S1 provides a possible explanation. In short, in the single-channel recording in the presence of GLPG1837, we occasionally observed long opening events that lasted for tens to hundreds of seconds even after ATP was removed from the solution (Fig. S1, A and B), also known as the locked-open phenomenon (Gunderson and Kopito, 1994). It is noteworthy that these locked-open events, rarely observed in WT channels gated with ATP alone, were seen more often in the presence of GLPG1837. Because we excluded the locked-open events in microscopic kinetic analysis, the Po was inevitably underestimated, resulting in a fold increase of microscopic Po lower than that of macroscopic currents.
Although a swath of reagents has been reported to enhance the gating function of CFTR, detailed mechanistic studies for these CFTR potentiators are limited (reviewed in Hwang and Sheppard, 1999; Rowe and Verkman, 2013; Jih et al., 2017). As described in the Introduction, we can divide the potentiators into three major categories based on their mechanisms of action: (1) ATP analogues, (2) VX-770 and nitrate, and (3) NPPB. Given that the chemical structure of GLPG1837 is distinct from ATP or its derivatives, it is unlikely that GLPG1837 works via the same mechanism as ATP analogues. However, whether the action of GLPG1837 is distinct from any other potentiators remains unclear. One unique pharmacological property distinguishing NPPB from VX-770 is that the current relaxation upon removal of ATP for hydrolysis-deficient mutants was shortened by NPPB (Lin et al., 2016) but prolonged by VX-770 (Kopeikin et al., 2014). We reasoned that this opposite effect on nonhydrolytic closing may serve as a tool to distinguish whether GLPG1837 works similarly to VX-770 or NPPB. Therefore, we tested the effects of GLPG1837 on the nonhydrolytic closing rate of E1371S-CFTR, a hydrolysis-deficient mutant whose ATPase activity was abolished by replacing the catalytic glutamate with a serine (Stratford et al., 2007). Fig. 4 A shows that the relaxation time constant upon removal of ATP in the absence of GLPG1837 can be fitted with a single-exponential function with a time constant of ∼50 s (red line), which was greatly prolonged to ∼103 s (blue line) in the presence of GLPG1837, suggesting that GLPG1837 decelerates nonhydrolytic closing of the channels. The relaxation time constants are summarized in Fig. 4 B.
This effect of GLPG1837 on nonhydrolytic closing was further confirmed in another hydrolysis-deficient mutant, D1370N/ΔR-CFTR (Csanády et al., 2010), which allows examinations of microscopic kinetics. As seen in Fig. 4 C, the Po of D1370N/ΔR-CFTR is ∼0.5 at 2 mM ATP; the opening bursts of D1370N/ΔR -CFTR are visibly prolonged by GLPG1837 in single-channel recordings. Of note, this set of experiments was performed in the ΔR background, which bears two distinct advantages: first, the experiments can be performed immediately after excising the patch without phosphorylation, and hence any dephosphorylation-induced effects are avoided; second, the reduced expression of the ΔR construct facilitates acquisition of single-channel data. In the presence of GLPG1837, the opening bursts for D1370N/ΔR -CFTR are lengthened to 1659 ± 295 ms (vs. 930 ± 143 ms in control; P < 0.05), and the interburst durations are considerably shortened from 694 ± 125 to 200 ± 39 ms, resulting in a Po of 0.87 ± 0.02 (n = 5, Fig. 4 D). Thus, data presented in Fig. 4 suggest that GLPG1837 may work differently than NPPB.
We next considered the possibility that GLPG1837 shares a similar mechanism of action with VX-770. Our current understanding for the action of VX-770 can be summarized in three major features: first, VX-770 increases both ATP-dependent and ATP-independent gating of CFTR (Eckford et al., 2012; Jih and Hwang, 2013); second, VX-770 slows down nonhydrolytic closing (Jih and Hwang, 2013; Kopeikin et al., 2014); third, VX-770 potentiates mutants with defective NBD dimerization such as G551D-CFTR (Van Goor et al., 2009; Yu et al., 2012) and ΔNBD2-CFTR (Yeh et al., 2015), presumably by perturbing the equilibrium between closed- and open-channel configurations in the TMDs (Jih and Hwang, 2013). Our experimental data of GLPG1837 presented so far match the first two features of VX-770. We therefore asked whether GLPG1837 acts on G551D-CFTR and ΔNBD2-CFTR as well.
Fig. 5 A shows that 3 µM GLPG1837 dramatically increases the currents of G551D-CFTR by ∼20-fold (19.85 ± 0.94-fold increase, n = 27). The maximal current increase at the saturating concentration of GLPG1837 (20 µM) reaches 35.62 ± 5.42-fold, which is larger than the maximal effect of VX-770 (∼6-fold increase of Po as shown in Van Goor et al., 2009, but ∼10-fold in Jih and Hwang, 2013). Interestingly, the biphasic current decay after ATP washout in the continuous presence of GLPG1837 was also seen with G551D-CFTR potentiated with VX-770, as reported in Lin et al. (2014). Subsequent removal of GLPG1837 results in further current decay, supporting the idea that GLPG1837 exerts its effects even in the absence of ATP. Inh-172 (Ma et al., 2002; Kopeikin et al., 2010) was applied at the end of the recording to ensure the consistence of the baseline. We further demonstrated that the response to GLPG1837 is independent of the formation of dimeric NBDs by treating ΔNBD2-CFTR with GLPG1837 (Fig. 5 B). The effectiveness of GLPG1837 on ΔNBD2-CFTR also reiterates the idea that GLPG1837 and NPPB work via different mechanisms, as deletion of NBD2 drastically reduces the efficacy of NPPB (Lin et al., 2016). Fig. 5 C depicts the dose–response relationships of GLPG1837 for G551D-CFTR and ΔNBD2-CFTR. Fitting of the Hill equation yields K1/2 values of 2.23 ± 0.25 and 1.19 ± 0.03 µM, respectively. It is interesting to note that the K1/2 values of G551D-CFTR and ΔNBD2-CFTR are ∼10-fold higher than that of WT-CFTR (0.23 ± 0.12 µM; black curve in Fig. 5 C). This finding prompted us to compare the K1/2 for GLPG1837 for mutants that exhibit similar gating defects, such as G1349D-CFTR (K1/2 = 1.49 ± 0.05 µM), and a mutation that does not alter the Po, D1370N-CFTR (K1/2 = 0.34 ± 0.07 µM). It appears that the dose–response curves are rightward-shifted in G551D-, ΔNBD2-, and G1349D-CFTR compared with WT- and D1370N-CFTR. The fundamental difference between the two groups (G551D-, ΔNBD2-, and G1349D-CFTR vs. WT- and D1370N-CFTR) is that the former has severe gating defects leading to a reduction of the Po (Cai et al., 2006; Bompadre et al., 2007; Miki et al., 2010), which is reflected by a large current fold increase upon application of GLPG1837, whereas D1370N-CFTR assumes a Po (Gunderson and Kopito, 1995) similar to that of WT channels, and the maximal current increase of both channels is only ∼2-fold (2.06 ± 0.08–fold for WT-CFTR and 1.59 ± 0.36–fold for D1370N-CFTR). As we plotted the relationship between maximal fold increase and K1/2 of GLPG1837, a positive correlation emerged (Fig. 6; R2 = 0.90). This finding raised the possibility that the apparent affinity of GLPG1837 may be a function of the Po of the channels: the lower the Po (i.e., the higher the maximal fold increase in Po), the lower the affinity.
One caveat in relating the Po and the apparent affinity of GLPG1837 (Fig. 6) is that different CFTR constructs were used; thus one cannot exclude the effects of mutation itself on the measured K1/2. Nonetheless, the hypothesis that the affinity of GLPG1837 is determined by the Po of the channels predicts two scenarios. First, if we performed the dose–response experiments on a channel whose Po has already been raised to a higher level by a different class of potentiators, the dose–response curve should be shifted to the left. Second, if we lower the Po of WT channels by decreasing the concentration of ATP, we should obtain a rightward-shifted dose–response relationship.
To test the first scenario, we conducted dose–response experiments on G551D-CFTR that had been exposed to a high-affinity ATP analogue, dPATP (Miki et al., 2010). Fig. 7 A shows that 20 µM dPATP potentiates G551D currents by 17.43 ± 1.28-fold (n = 17; Miki et al., 2010). The currents can be further enhanced by GLPG1837 in a dose-dependent manner, with a K1/2 of 0.47 ± 0.02 µM (Fig. 7 B). Thus, when the Po of G551D-CFTR is increased by dPATP, the apparent affinity for GLPG1837 also increases. It should be noted that the dose–response shift is not caused by a saturation of the currents by dPATP plus GLPG1837, as Fig. S2 shows that the G551D-CFTR currents in the presence of saturating concentrations of dPATP and GLPG1837 can be further enhanced by nitrate.
Further evidence that dPATP increases the apparent affinity of GLPG1837 comes from experiments that assess the apparent off-rate of GLPG1837 upon washout of the reagent. Fig. 7 C shows that the current relaxation upon removal of GLPG1837 in the presence of dPATP has a time constant of 18.4 s, which is shortened to 3.2 s when dPATP is replaced with ATP. The time constants in both conditions are summarized in Fig. 7 D (n = 4; P < 0.01). Of note, the dissociation of dPATP is decelerated by GLPG1837 as well. Fig. S3 A shows an experimental protocol similar to that of Fig. 7 C except that the current relaxation upon removal of dPATP was better fitted with a double-exponential function. Both the first (τ1) and second (τ2) time constants are prolonged to 10.8 and 178.4 s (blue line) in the presence of GLPG1837 (τ1 = 4.5 s and τ2 = 62.4 s in the absence of GLPG1837, red line). Comparisons of both time constants are summarized in Fig. S3 B.
So far the first scenario is confirmed: the apparent affinity of GLPG1837 is higher for G551D-CFTR when its Po has been elevated by other potentiators. We next used WT-CFTR to test the validity of the second scenario. The idea of Po-dependent affinity predicts that as we lower the concentration of ATP, which lowers the Po of WT channels, we should see a decrease in the apparent affinity of GLPG1837. Indeed, Fig. 8 A shows a representative trace of dose–response experiments of GLPG1837 conducted in the presence of 2.5 µM ATP. Although 3 µM is already the saturating concentration for GLPG1837 on WT-CFTR gated by 2 mM ATP (Fig. 2 B), the current at 2.5 µM ATP, in contrast, does not saturate at such a concentration of GLPG1837. In fact, we observe a further increase of the currents upon addition of 20 µM GLPG1837, with the K1/2 shifted to 1.77 ± 0.65 µM (Fig. 8 B, blue curve; c.f. K1/2 = 0.23 ± 0.12 µM at 2 mM ATP), supporting our idea that the apparent affinity of GLPG1837 decreases as the Po declines. In addition, when we completely removed ATP from the solution so that only residual currents were left, increasing [GLPG1837] from 3 to 20 µM enhanced “ATP-independent” gating of WT-CFTR (box in Fig. S4 A, expanded in Fig. S4 B). Although the current in the absence of ATP is too small to grant an accurate quantification for a complete dose–response relationship, these results echo the finding described in Fig. 8. Collectively, our data provide evidence that the apparent affinity of GLPG1837 is a function of the Po of the channels regardless of the strategies we used to manipulate the Po, including mutagenesis, combination of other potentiators, and changing ATP concentrations. The pharmacological and mechanistic implications will be elaborated in the Discussion.
Data presented in Fig. 7 also indicate that GLPG1837 and dPATP work through different mechanisms, as the former can greatly enhance the G551D-CFTR currents in the presence of a saturating concentration of dPATP. This idea is perhaps not surprising because the structure of GLPG1837 (Fig. 1 C) does not resemble that of ATP. GLPG1837 is also chemically distinct from NPPB (Fig. 1 B). They do not share a common mechanism, as a decrease of the nonhydrolytic closing rate observed on E1371S-CFTR treated with GLPG1837 (Fig. 4) contradicts the effect of NPPB. Interestingly, notwithstanding the unique Po-dependent changes in the apparent affinity of GLPG1837, the action of GLPG1837 has several features that resemble the characteristics of VX-770. They share the following similar effects as described previously: first, both ATP-dependent and -independent gating of WT-CFTR are increased (Fig. 2; Jih and Hwang, 2013); second, nonhydrolytic channel closing is decelerated (Fig. 4; Kopeikin et al., 2014); and third, gating potentiation occurs in mutants with defects in NBD dimerization (Fig. 5; Lin et al., 2014; Yeh et al., 2015). However, the fact that GLPG1837 is chemically distinguishable from VX-770 and its analogues (Hadida et al., 2014) leads us to ask how two structurally different compounds share such similarities in modulating CFTR gating.
To address this question, we tested the hypothesis that GLPG1837 and VX-770 may share a common binding site in CFTR. We reasoned that if GLPG1837 and VX-770 potentiate CFTR gating via the same binding site, occupancy of the binding site by one compound should preclude binding of the other. We chose G551D-CFTR to test this idea because the effects of both compounds have been carefully quantified on this mutant. Fig. 9 A shows a continuous recording of macroscopic G551D-CFTR currents preactivated with PKA and ATP. The application of 200 nM VX-770 increases the currents by ∼10-fold as expected, but more importantly, further addition of 3 µM GLPG1837 in the presence of VX-770 fails to affect the current amplitude despite the fact that GLPG1837, when applied alone, is more efficacious than VX-770 (∼20-fold for GLPG1837 vs. ∼10-fold for VX-770). On the other hand, when we reversed the protocol and applied GLPG1837 first, the currents were increased by 20-fold as expected, but the following addition of VX-770 actually lowered the current to the level equivalent to that by VX-770 alone (Fig. 9 B). Subsequent removal of VX-770 failed to bring the current back to the original level in GLPG1837. Although other more complex mechanisms can account for these results (see Supplement), the simplest explanation is that VX-770 and GLPG1837 are competing for the same site in CFTR. Of note, VX-770 assumes an extremely high affinity for CFTR (EC50 in the low nanomolar range, as reported in Hadida et al. [2014]); indeed, within the experimental time of tens of minutes, removal of VX-770 from the perfusate does not eliminate its effects (Jih and Hwang, 2013). Thus, if we consider the low-affinity GLPG1837 as a full agonist for the binding site, the high-affinity VX-770 is a partial agonist. The idea of competition for the same site predicts that once the binding site is occupied by the high-affinity VX-770, GLPG1837 will not be able to bind to exert its own effects on the channel because of a negligible dissociation of VX-770. In contrast, when the binding site is first taken by the low-affinity ligand GLPG1837, the high-affinity partial agonist VX-770 can replace GLPG1837. But, because of the lower efficacy of VX-770, the GLPG1837-potentiated currents are reduced gradually as more and more channels have their initial ligand substituted by VX-770.
Discussion
In this study, we investigated the detailed mechanism of an alternative CFTR potentiator in development, GLPG1837, which recently has entered phase II clinical trials in patients harboring the G551D or S1251N CFTR mutation. We provide evidence that GLPG1837 shares a common mechanism for gating potentiation with VX-770 (ivacaftor), an US Food and Drug Administration–approved drug for the treatment of CF (Van Goor et al., 2009). Because the high-affinity VX-770 obviates the effects of GLPG1837, the simplest explanation for this shared mechanism is that these two reagents bind to a common binding site. In addition, our results reveal a Po-dependent dose–response relationship for GLPG1837: the higher the Po, the smaller the K1/2. In this section, we first discuss the modi operandi of CFTR potentiators based on our current understanding of the gating mechanism of CFTR, and then use the idea of classic allostery to explain the Po-dependent apparent affinity of GLPG1837. Finally, we speculate on the structural implications of the binding site for GLPG1837 and VX-770 based on their chemical characteristics.
It is generally accepted that once phosphorylated, opening of CFTR’s gate is coupled to ATP binding-induced NBD dimerization, and gate closure is facilitated by ATP hydrolysis-catalyzed partial (or complete) separation of the NBDs (Vergani et al., 2003, 2005; Hwang and Sheppard, 2009; Tsai et al., 2009, 2010). Although the exact mechanism by which NBD dimerization is coupled to the movement of the gate in TMDs is still debated (Sohma and Hwang, 2015), there are three critical steps involved in the gating cycle of CFTR: (1) ATP binding, (2) formation of dimeric NBDs, and (3) opening of the gate in TMDs. Therefore, any reagents that provide a favorable condition for either step could act as a CFTR potentiator. For example, ATP analogues PATP, dATP, and dPATP increase the Po partly because they bind better and partly because they catalyze NBD dimerization more efficiently (Aleksandrov et al., 2002; Zhou et al., 2005; Miki et al., 2010). In contrast, VX-770 is proposed to enhance CFTR functions by modulating the gating transitions between the open and closed conformations in the TMDs. Although Csanády and Töröcsik (2014) provided evidence for modulation of the transition state energy by NPPB, a CFTR potentiator that also blocks the pore (Wang et al., 2005), Lin et al. (2016) suggested that NPPB may work by facilitating NBD dimerization. Furthermore, because of the energetic coupling between NBD dimerization and gate opening (Jih et al., 2012), a synergistic but dependent interaction between NPPB and VX-770 is observed (Lin et al., 2016).
In the current study, we demonstrated two types of pharmacological interactions. On one hand, GLPG1837 and VX-770 seem to work through a competitive mechanism by which binding of one reagent precludes the effects of the other, as though they are competing for the same site. On the other hand, GLPG1837 and dPATP show a clear synergism, as even at a maximally effective concentration of dPATP, GLPG1837 can further increase the Po of G551D-CFTR (Fig. 7 A). Moreover, the magnitude of potentiation by GLPG1837 depends on the presence or absence of dPATP. As depicted in Figs. 5 A and 7 A, 3 µM GLPG1837 and 20 µM dPATP cause ∼20- and ∼17-fold increases of G551D-CFTR currents, respectively, whereas the combined effect of both reagents is 94.55 ± 8.07-fold (n = 14), which is larger than the sum of the individual effects (∼37-fold), indicating a pharmacological synergy (Goodman, 2011), but smaller than the value predicted if these two potentiators worked in an independent manner. (See Yeh et al. [2015], for example.)
Despite the differences in the exact kinetic steps acted on by different potentiators, one fundamental feature all potentiators must possess is that binding of the potentiators promotes channel opening. To elucidate how two separate events, channel opening and potentiator binding, are coupled, we borrowed a general four-state allosteric modulation model used successfully to explain the mechanisms of action for classical ligand-gated channels (Fig. 10 A). Based on the assumptions that ligand can bind to both the closed and open states, and that channel opening can occur whether or not the ligand is bound, this model depicts binding and gate opening as separate but coupled transitions, where the thermodynamic laws demand that if ligand binding promotes opening, then opening must promote ligand binding. It should be noted that other models that allow ligand binding only to the closed or open state may also explain state-dependent binding (see Supplement). Nevertheless, we decided to use a more generalized allosteric model, as it is less constrained.
If we assign the equilibrium constants of gate opening/closing for the ligand-free state and ligand-bound state as KO and fKO, respectively, where f is greater than 1 for a ligand that increases the Po of the channel, the opening from the ligand-bound state is more favorable than opening from the ligand-free state by a factor of f. The free energy change (ΔΔG) needed to alter the equilibrium gating constant is –RTln(f) kcal/mol. In a system without other sources of energy input, this ΔΔG then comes from a stronger binding affinity of the ligand for the channel in the open state (KDO) than the closed state (KDC), where KDC and KDO also differ by a factor of f. In other words, binding of the ligand stabilizes the open state relative to the closed state by −RTln(f) kcal/mol.
Mathematically, this model also predicts that the apparent affinity (K1/2) measured in the dose–response experiments should lie between KDC and KDO (see Supplement and Figs. S5 and S6 for details). Of note, because the Po of G551D-CFTR is extremely low even in the presence of CFTR potentiators, the formula for Po calculation can be simplified to ko/kc, where ko and kc represent, respectively, the opening and closing rates. Thus, the maximal macroscopic current fold increase (or efficacy) of G551D-CFTR by GLPG1837 will be approximately the same as f, which is determined by the differences of free energy of binding between the open and closed states as described in the previous paragraph. In contrast, because the measured apparent affinity from the dose–response relationship always lies between the true binding affinity of the open state and that of the closed state, the potency of GLPG1837 is determined by how tight the compound binds to the open state as well as the closed state.
The idea of state-dependent binding implies that the structure of the binding site for GLPG1837 in an open state must differ from that in a closed state, leading to different affinities. However, a maximal fold increase of ∼35 for G551D-CFTR is equivalent to a gain of ∼2.1 kcal/mol of free energy for stabilizing the open state. This magnitude of ΔΔG can be accomplished by changes of just one hydrogen bond (Berg et al., 2002), suggesting that the binding site for GLPG1837 should not undergo a drastic structural alteration during gating. Although thermodynamic analysis only offers limited insights into the drug binding site itself, the finding that GLPG1837 and VX-770 may compete for the same site provides clues for the possible compositions of the amino acids participating in the formation of the binding pocket. Comparing the chemical structures of GLPG1837 and VX-770, we noticed that, in addition to a hydrophobic nature, they are both composed of aromatic rings linked by a peptide bond, a feature also shared by other known CFTR potentiators (Yang et al., 2003; Pedemonte et al., 2005). The exact length of the peptide bond may position the two opposing aromatic rings close to their respective binding partners that may possess aromatic and/or hydrophobic amino acids to stabilize binding. Furthermore, if, as proposed (Jih and Hwang, 2013; Yeh et al., 2015), the binding site indeed resides in the lipid–protein interface, membrane lipids should also contribute significantly to the hydrophobic interactions involved in drug binding.
Although GLPG1837 has an efficacy threefold larger than VX-770 on G551D-CFTR, its potency, reflected by the apparent affinity within a range of 0.2∼2 µM (see Results for details), is much less than that of VX-770 (EC50 at low nanomolar range; Hadida et al., 2014). The extremely tight binding (i.e., high affinity) of VX-770 to CFTR is consistent with the observation that its effect cannot be washed out by continuous perfusion of VX-770–free solution in the experimentally permissive time frame (Jih and Hwang, 2013), whereas removal of GLPG1837 from the solution readily brings the current back to the control level (Fig. 5, A and B). Technically, this tight binding of VX-770 makes it virtually impossible to obtain an accurate dose–response relationship to access its state-dependent binding. Looking into the physical properties of these two reagents, we noted a large difference in the partition coefficients (P) between VX-770 (logP = 6.34 ± 0.78) and GLPG1837 (logP = 3.22 ± 0.73). The hydrophobic nature of VX-770 explains why visible precipitates were seen at [VX-770] >1 µM in our perfusion solution. We also speculate that this hydrophobicity property may partly account for the high affinities of VX-770. Although more studies are required to examine how VX-770 possesses a higher affinity yet lower efficacy than GLPG1837, our results support the aforementioned thermodynamic argument that efficacy and potency are two separate physicochemical properties of a ligand.
Regardless of the detailed mechanism, the state-dependent binding of potentiators demonstrated in the current work lays the foundation for guiding future designs of therapeutic strategies. The observation that the apparent affinity of GLPG1837 is increased on dPATP-potentiated channels (and vice versa) raises the possibility of manipulating the effective concentration of one potentiator by first exposing the channels to another class of potentiators. In other words, combining two potentiators with different mechanisms of action not only results in pharmacological synergism in efficacy, but also mutually enhances their potency, which may be beneficial in a sense that such combination could lower the drug dosage used in clinics. In this aspect, experiments testing the idea of state-dependent binding for other CFTR potentiators are urgently needed. Moreover, we believe that a fundamental understanding of the mechanisms of action for CFTR potentiators, together with the exquisite molecular insights from the atomic structure of human CFTR (Liu et al., 2017), should pave the way for structure-based drug design targeting CFTR.
Acknowledgments
We thank Cindy Chu and Shenghui Hu for their technical assistance and Drs. Ashvani Singh and Chris Tse for their helpful discussion. We also thank Galapagos/AbbVie for providing GLPG1837.
This work is supported by the National Institutes of Health (grant R01DK55835) and the Cystic Fibrosis Foundation (grant Hwang11P0) to T.-C. Hwang.
Part of the work is supported by a service agreement with AbbVie Inc. The authors declare no additional competing financial interests.
Author contributions: conception and design of experiments: H.-I Yeh and T.-C. Hwang; collection and data analysis: H.-I Yeh; interpretation of data: H.-I Yeh and T.-C. Hwang; design of computer simulation program: Y. Sohma; drafting the article: H.-I Yeh, Y. Sohma, K. Conrath, and T.-C. Hwang; it is confirmed that all authors approved the final version of the manuscript. All experiments were performed in Dalton Cardiovascular Research Center, University of Missouri-Columbia.
Richard W. Aldrich served as editor.
References
- CF
cystic fibrosis
- dPATP
2′-deoxy-N6-(2-phenylethyl)-adenosine-5′-O-triphosphate
- GLPG1837
N-(3-carbamoyl-5,5,7,7-tetramethyl-4,7-dihydro-5H-thieno[2,3-c]pyran-2-yl)-1H-pyrazole-3-carboxamide
- NBD
nucleotide-binding domain
- NPPB
5-nitro-2-(3-phenylpropylamino)benzoate
- TMD
transmembrane domain
- VX-770
N-(2,4-di-tert-butyl-5-hydroxyphenyl)-4-oxo-1,4-dihydroquinoline-3-carboxamide)