Structure-based discovery of CFTR potentiators and inhibitors

Abstract

The cystic fibrosis transmembrane conductance regulator (CFTR) is a crucial ion channel whose loss of function leads to cystic fibrosis, while its hyperactivation leads to secretory diarrhea. Small molecules that improve CFTR folding (correctors) or function (potentiators) are clinically available. However, the only potentiator, ivacaftor, has suboptimal pharmacokinetics and inhibitors have yet to be clinically developed. Here we combine molecular docking, electrophysiology, cryo-EM, and medicinal chemistry to identify novel CFTR modulators. We docked ~155 million molecules into the potentiator site on CFTR, synthesized 53 test ligands, and used structure-based optimization to identify candidate modulators. This approach uncovered novel mid-nanomolar potentiators as well as inhibitors that bind to the same allosteric site. These molecules represent potential leads for the development of more effective drugs for cystic fibrosis and secretory diarrhea, demonstrating the feasibility of large-scale docking for ion channel drug discovery.

Extended Data Figure 1.. Z2075279358 increased the ion flux rate in a bronchial epithelial cell line derived from a patient homozygous for ΔF508 (CFBE41o⁻ cells expressing ΔF508-CFTR and the fluorescent protein eYFP-H148Q/I152L/F46L).

Cells were treated with 20 μM forskolin +/− 1 μM test compound. The relative potentiation was calculated as the ratio of flux rates with and without compound. Data points represent a single measurement for each test compound. For ivacaftor and GLPG1837 data represent means and SEs for 9 or 2 measurements, respectively.

Extended Data Figure 2.. Cryo-EM reconstructions of the CFTR-’853 complex.

Summary of the image processing procedure.

Extended Data Figure 3.

SAR for optimization of compound ‘853.

Extended Data Figure 4.

Docked poses of I1421, I1408, and (S)-SX263.

Figure 1:. Ultra-large docking screen identifies novel CFTR potentiators.

(a) The workflow of this study. (b) Compound Z2075279358 (‘358) potentiates ΔF508 CFTR. CFBE41o⁻ cells homozygous for ΔF508 CFTR were treated with 1 μM lumacaftor and 20 μM forskolin. The relative potentiation was calculated as the ratio of flux rates with and without potentiator. Data points represent the means and standard errors (SEs) of 6 to 8 measurements (each shown as a dot). (c) Potentiation activity of 10 μM GLPG1837 or 5 μM compound against WT CFTR fused to a carboxy-terminal GFP tag. Inside-out membrane patches containing WT CFTR were excised from CHO cells and then fully phosphorylated by protein kinase A (PKA) in the presence of 3 mM ATP. The fold stimulation is defined as the ratio of the current in the presence and absence of added compound. Data represent means and SEs of 3–33 patches with individual measurements shown as dots. (d) The 2D structures of the potentiators GLPG1837, ivacaftor, and the 13 positive hits from the initial screen. (e) Representative macroscopic current trace and dose-response curve of WT CFTR in response to perfusion with ‘358. CFTR-containing membrane patches were fully phosphorylated by PKA. The current in the presence of 3 mM ATP before titration was used to normalize the current potentiated by different concentrations of ‘358. The EC₅₀ is estimated to be 2.2 ± 0.6 μM by fitting the dose-responses with the Hill equation. Data represent means and SEs from 3 patches.

Figure 2:. New modulators identified through an analog screen.

(a) Potentiation activity of ‘358 analogs. All compounds were tested at a single concentration of 10 M in inside-out membrane patches containing fully phosphorylated WT CFTR. The current stimulation levels of GLPG1837 and ‘358 are indicated as dashed lines. (b) Reported structures versus the NMR-determined structures. (c) Docked poses of the reported structure of ‘853 (left) versus the NMR-determined structure of ‘853 (right). (d) (S)-’853 potentiates WT CFTR currents, while (R)-’853 mildly inhibits them in inside-out patches. Both enantiomers were perfused at 100 M concentration. Data represent means and SEs of 20 ((S)-’853) or 18 ((R)-’853) patches. Statistical significance relative to no effect was tested by two-tailed Student’s t-test (****P = 2.3×10⁻⁹ for (S)-’853 and P = 1.6×10⁻⁷ for (R)-’853). (e) Competition assay showing that the presence of (R)-’853 (100 M) diminishes the potentiating effect of (S)-’853 (100 M). (f) Dose-response curve of (S)-’853 versus the racemic mixture (±)-’853, as described for Figure 1e. The EC₅₀ for (S)-’853 was estimated to be 2.1 ± 0.9 μM. Data represent means and SEs of 2–20 patches. 3 mM ATP was used in all panels.

Figure 3:. Z1834339853 binds to the same site as ivacaftor and GLPG1837.

(a) Cryo-EM structure of phosphorylated and ATP-bound CFTR (E1371Q) in complex with ‘853. (b) Zoomed-in views of the density of ‘853 (top) and a comparison between the docked pose (salmon) and the cryo-EM pose (magenta). (c) Representative macroscopic current traces and dose-response curves of fully phosphorylated WT, S308A, and Y304A CFTR in response to perfusion of (S)-’853 onto inside-out excised membrane patches. 3 mM ATP was used. Each data point represents the mean and SEs determined from 3 to 12 patches.

Figure 4:. Medicinal chemistry leads to novel CFTR potentiators and inhibitors.

(a) General formula of newly synthesized ‘853 analogs. (b) Effects of ‘853 analogs on currents measured in inside-out excised membrane patches containing fully phosphorylated WT CFTR. Measurements were made with 3 mM ATP. Data represent means and SEs of 2 to 46 patches. (c) The structures and dose-response curves of three of the most efficacious potentiators. Data represent means and SEs of 5–11 patches. (d) The structures and dose-response curves of the two most efficacious inhibitors. The dose-responses of the Y304A and S308A variants in response to I1412 perfusion were also shown (left). Data represent means and SEs of 2–8 patches.

Figure 5:. The activity of I1421 against 10 CF-causing mutations.

(a) The positions of the mutations mapped onto dephosphorylated and ATP-free CFTR (PDB 5UAK). (b) Representative macroscopic current traces in response to I1421 (10 μM) perfusion onto inside-out membrane patches excised from CHO cells. 3 mM ATP was used. (c) Potentiation activity of I1421 versus GLPG1837. The mean and SE values were determined from 2 to 7 patches. (d) Pharmacokinetic analysis of compound I1421. Plasma concentration-time profiles in male C57BL/6N mice following a single subcutaneous (SC), intraperitoneal (IP), per-oral (PO) (dose 10 mg/kg) or intravenous (IV) (3 mg/kg) administration. Data represent means and SDs. (e) Selected pharmacokinetic parameters of I1421. C_max: peak plasma concentration; T_max: the time when the peak plasma concentration was observed; AUC_last: the areas under the concentration time curve; T_1/2: terminal half-life; CL: clearance, V_ss: steady-state volume of distribution; %F: %bioavailability.