Binding screen for cystic fibrosis transmembrane conductance regulator correctors finds new chemical matter and yields insights into cystic fibrosis therapeutic strategy

Inhibitory targeting cGAS-STING-TBK1 axis: Emerging strategies for autoimmune diseases therapy

The cGAS-STING signaling plays an integral role in the host immune response, and the abnormal activation of cGAS-STING is highly related to various autoimmune diseases. Therefore, targeting the cGAS-STING-TBK1 axis has become a promising strategy in therapy of autoimmune diseases. Herein, we summarized the key pathways mediated by the cGAS-STING-TBK1 axis and various cGAS-STING-TBK1 related autoimmune diseases, as well as the recent development of cGAS, STING, or TBK1 selective inhibitors and their potential application in therapy of cGAS-STING-TBK1 related autoimmune diseases. Overall, the review highlights that inhibiting cGAS-STING-TBK1 signaling is an attractive strategy for autoimmune disease therapy.

Comprehensive Analysis of Chemical Structures That Have Been Tested as CFTR Activating Substances in a Publicly Available Database CandActCFTR

Background: Cystic fibrosis (CF) is a genetic disease caused by mutations in CFTR, which encodes a chloride and bicarbonate transporter expressed in exocrine epithelia throughout the body. Recently, some therapeutics became available that directly target dysfunctional CFTR, yet research for more effective substances is ongoing. The database CandActCFTR aims to provide detailed and comprehensive information on candidate therapeutics for the activation of CFTR-mediated ion conductance aiding systems-biology approaches to identify substances that will synergistically activate CFTR-mediated ion conductance based on published data. Results: Until 10/2020, we derived data from 108 publications on 3,109 CFTR-relevant substances via the literature database PubMed and further 666 substances via ChEMBL; only 19 substances were shared between these sources. One hundred and forty-five molecules do not have a corresponding entry in PubChem or ChemSpider, which indicates that there currently is no single comprehensive database on chemical substances in the public domain. Apart from basic data on all compounds, we have visualized the chemical space derived from their chemical descriptors via a principal component analysis annotated for CFTR-relevant biological categories. Our online query tools enable the search for most similar compounds and provide the relevant annotations in a structured way. The integration of the KNIME software environment in the back-end facilitates a fast and user-friendly maintenance of the provided data sets and a quick extension with new functionalities, e.g., new analysis routines. CandActBase automatically integrates information from other online sources, such as synonyms from PubChem and provides links to other resources like ChEMBL or the source publications. Conclusion: CandActCFTR aims to establish a database model of candidate cystic fibrosis therapeutics for the activation of CFTR-mediated ion conductance to merge data from publicly available sources. Using CandActBase, our strategy to represent data from several internet resources in a merged and organized form can also be applied to other use cases. For substances tested as CFTR activating compounds, the search function allows users to check if a specific compound or a closely related substance was already tested in the CF field. The acquired information on tested substances will assist in the identification of the most promising candidates for future therapeutics.

Pharmacological chaperones improve intra-domain stability and inter-domain assembly via distinct binding sites to rescue misfolded CFTR

Protein misfolding is involved in a large number of diseases, among which cystic fibrosis. Complex intra- and inter-domain folding defects associated with mutations in the cystic fibrosis transmembrane regulator (CFTR) gene, among which p.Phe508del (F508del), have recently become a therapeutical target. Clinically approved correctors such as VX-809, VX-661, and VX-445, rescue mutant protein. However, their binding sites and mechanisms of action are still incompletely understood. Blind docking onto the 3D structures of both the first membrane-spanning domain (MSD1) and the first nucleotide-binding domain (NBD1), followed by molecular dynamics simulations, revealed the presence of two potential VX-809 corrector binding sites which, when mutated, abrogated rescue. Network of amino acids in the lasso helix 2 and the intracellular loops ICL1 and ICL4 allosterically coupled MSD1 and NBD1. Corrector VX-445 also occupied two potential binding sites on MSD1 and NBD1, the latter being shared with VX-809. Binding of both correctors on MSD1 enhanced the allostery between MSD1 and NBD1, hence the increased efficacy of the corrector combination. These correctors improve both intra-domain folding by stabilizing fragile protein-lipid interfaces and inter-domain assembly via distant allosteric couplings. These results provide novel mechanistic insights into the rescue of misfolded proteins by small molecules.

Differential Effects of Oleuropein and Hydroxytyrosol on Aggregation and Stability of CFTR NBD1-ΔF508 Domain

Cystic Fibrosis (CF) is caused by loss of function mutations in the Cystic Fibrosis transmembrane conductance regulator (CFTR). The folding and assembly of CFTR is inefficient. Deletion of F508 in the first nucleotide binding domain (NBD1-ΔF508) further disrupts protein stability leading to endoplasmic reticulum retention and proteasomal degradation. Stabilization and prevention of NBD1-ΔF508 aggregation is critical to rescuing the folding and function of the entire CFTR channel. We report that the phenolic compounds Oleuropein and Hydroxytryosol reduce aggregation of NBD1-ΔF508. The NBD1-ΔF508 aggregate size was smaller in the presence of Hydroxytryosol as determined by dynamic light scattering. Neither phenolic compound increased the thermal stability of NBD1-ΔF508 as measured by differential scanning fluorimetry. Interestingly, Hydroxytyrosol inhibited the stabilizing effect of the indole compound BIA, a known stabilizer, on NBD1-ΔF508. Molecular docking studies predicted that Oleuropein preferred to bind in the F1-type core ATP-binding subdomain in NBD1. In contrast, Hydroxytyrosol preferred to bind in the α4/α5/α6 helical bundle of the ABCα subdomain of NBD1 next to the putative binding site for BIA. This result suggests that Hydroxytyrosol interferes with BIA binding, thus providing an explanation for the antagonistic effect on NBD1 stability upon incubation with both compounds. To our knowledge, these studies are the first to explore the effects of these two phenolic compounds on the aggregation and stability of NBD1-ΔF508 domain of CFTR.

Stability Prediction for Mutations in the Cytosolic Domains of Cystic Fibrosis Transmembrane Conductance Regulator

Cystic Fibrosis (CF) is caused by mutations to the Cystic Fibrosis Transmembrane Conductance Regulator (CFTR) chloride channel. CFTR is composed of two membrane spanning domains, two cytosolic nucleotide-binding domains (NBD1 and NBD2) and a largely unstructured R-domain. Multiple CF-causing mutations reside in the NBDs and some are known to compromise the stability of these domains. The ability to predict the effect of mutations on the stability of the cytosolic domains of CFTR and to shed light on the mechanisms by which they exert their effect is therefore important in CF research. With this in mind, we have predicted the effect on domain stability of 59 mutations in NBD1 and NBD2 using 15 different algorithms and evaluated their performances via comparison to experimental data using several metrics including the correct classification rate (CCR), and the squared Pearson correlation (R²) and Spearman's correlation (ρ) calculated between the experimental ΔT_m values and the computationally predicted ΔΔG values. Overall, the best results were obtained with FoldX and Rosetta. For NBD1 (35 mutations), FoldX provided R² and ρ values of 0.64 and -0.71, respectively, with an 86% correct classification rate (CCR). For NBD2 (24 mutations), FoldX R², ρ, and CCR were 0.51, -0.73, and 75%, respectively. Application of the Rosetta high-resolution protocol (Rosetta_hrp) to NBD1 yielded R², ρ, and CCR of 0.64, -0.75, and 69%, respectively, and for NBD2 yielded R², ρ, and CCR of 0.29, -0.27, and 50%, respectively. The corresponding numbers for the Rosetta's low-resolution protocol (Rosetta_lrp) were R² = 0.47, ρ = -0.69, and CCR = 69% for NBD1 and R² = 0.27, ρ = -0.24, and CCR = 63% for NBD2. For NBD1, both algorithms suggest that destabilizing mutations suffer from destabilizing vdW clashes, whereas stabilizing mutations benefit from favorable H-bond interactions. Two triple consensus approaches based on FoldX, Rosetta_lpr, and Rosetta_hpr were attempted using either "majority-voting" or "all-voting". The all-voting consensus outperformed the individual predictors, albeit on a smaller data set. In summary, our results suggest that the effect of mutations on the stability of CFTR's NBDs could be largely predicted. Since NBDs are common to all ABC transporters, these results may find use in predicting the effect and mechanism of the action of multiple disease-causing mutations in other proteins.

New insights into structure and function of bis-phosphinic acid derivatives and implications for CFTR modulation

C407 is a compound that corrects the Cystic Fibrosis Transmembrane Conductance Regulator (CFTR) protein carrying the p.Phe508del (F508del) mutation. We investigated the corrector effect of c407 and its derivatives on F508del-CFTR protein. Molecular docking and dynamics simulations combined with site-directed mutagenesis suggested that c407 stabilizes the F508del-Nucleotide Binding Domain 1 (NBD1) during the co-translational folding process by occupying the position of the p.Phe1068 side chain located at the fourth intracellular loop (ICL4). After CFTR domains assembly, c407 occupies the position of the missing p.Phe508 side chain. C407 alone or in combination with the F508del-CFTR corrector VX-809, increased CFTR activity in cell lines but not in primary respiratory cells carrying the F508del mutation. A structure-based approach resulted in the synthesis of an extended c407 analog G1, designed to improve the interaction with ICL4. G1 significantly increased CFTR activity and response to VX-809 in primary nasal cells of F508del homozygous patients. Our data demonstrate that in-silico optimized c407 derivative G1 acts by a mechanism different from the reference VX-809 corrector and provide insights into its possible molecular mode of action. These results pave the way for novel strategies aiming to optimize the flawed ICL4-NBD1 interface.

Molecular Docking and QSAR Studies as Computational Tools Exploring the Rescue Ability of F508del CFTR Correctors

Cystic fibrosis (CF) is the autosomal recessive disorder most recurrent in Caucasian populations. Different mutations involving the cystic fibrosis transmembrane regulator protein (CFTR) gene, which encodes the CFTR channel, are involved in CF. A number of life-prolonging therapies have been conceived and deeply investigated to combat this disease. Among them, the administration of the so-called CFTR modulators, such as correctors and potentiators, have led to quite beneficial effects. Recently, based on QSAR (quantitative structure activity relationship) studies, we reported the rational design and synthesis of compound 2, an aminoarylthiazole-VX-809 hybrid derivative exhibiting promising F508del-CFTR corrector ability. Herein, we explored the docking mode of the prototype VX-809 as well as of the aforementioned correctors in order to derive useful guidelines for the rational design of further analogues. In addition, we refined our previous QSAR analysis taking into account our first series of in-house hybrids. This allowed us to optimize the QSAR model based on the chemical structure and the potency profile of hybrids as F508del-CFTR correctors, identifying novel molecular descriptors explaining the SAR of the dataset. This study is expected to speed up the discovery process of novel potent CFTR modulators.

Recent Strategic Advances in CFTR Drug Discovery: An Overview

Cystic fibrosis transmembrane conductance regulator (CFTR)-rescuing drugs have already transformed cystic fibrosis (CF) from a fatal disease to a treatable chronic condition. However, new-generation drugs able to bind CFTR with higher specificity/affinity and to exert stronger therapeutic benefits and fewer side effects are still awaited. Computational methods and biosensors have become indispensable tools in the process of drug discovery for many important human pathologies. Instead, they have been used only piecemeal in CF so far, calling for their appropriate integration with well-tried CF biochemical and cell-based models to speed up the discovery of new CFTR-rescuing drugs. This review will give an overview of the available structures and computational models of CFTR and of the biosensors, biochemical and cell-based assays already used in CF-oriented studies. It will also give the reader some insights about how to integrate these tools as to improve the efficiency of the drug discovery process targeted to CFTR.

A simple model for determining affinity from irreversible thermal shifts

Thermal denaturation (Tm) data are easy to obtain; it is a technique that is used by both small labs and large‐scale industrial organizations. The link between ligand affinity (K_D) and ΔTm is understood for reversible denaturation; however, there is a gap in our understanding of how to quantitatively interpret ΔTm for the many proteins that irreversibly denature. To better understand the origin, and extent of applicability, of a K_D to ΔTm correlate, we define equations relating K_D and ΔTm for irreversible protein unfolding, which we test with computational models and experimental data. These results suggest a general relationship exists between K_D and ΔTm for irreversible denaturation.

Domain-interface dynamics of CFTR revealed by stabilizing nanobodies

The leading cause of cystic fibrosis (CF) is the deletion of phenylalanine 508 (F508del) in the first nucleotide-binding domain (NBD1) of the cystic fibrosis transmembrane conductance regulator (CFTR). The mutation affects the thermodynamic stability of the domain and the integrity of the interface between NBD1 and the transmembrane domain leading to its clearance by the quality control system. Here, we develop nanobodies targeting NBD1 of human CFTR and demonstrate their ability to stabilize both isolated NBD1 and full-length protein. Crystal structures of NBD1-nanobody complexes provide an atomic description of the epitopes and reveal the molecular basis for stabilization. Furthermore, our data uncover a conformation of CFTR, involving detachment of NBD1 from the transmembrane domain, which contrast with the compact assembly observed in cryo-EM structures. This unexpected interface rearrangement is likely to have major relevance for CF pathogenesis but also for the normal function of CFTR and other ABC proteins.

Differential Scanning Fluorimetry and Hydrogen Deuterium Exchange Mass Spectrometry to Monitor the Conformational Dynamics of NBD1 in Cystic Fibrosis

Cystic fibrosis (CF) is one of the most common, lethal autosomal recessive diseases in Caucasians with a life expectancy of 37-47 years. The CF transmembrane conductance regulator (CFTR) is a plasma membrane ion channel, confined to apical membrane of epithelia, and ensures transepithelial water and solute movement across secretory epithelia in several organs. Numerous CF mutations, including the most prevalent deletion of F508 (ΔF508) in the nucleotide binding domain 1 (NBD1) leads to CFTR global misfolding and premature intracellular degradation at the endoplasmic reticulum (ER). To better understand the misfolding mechanism caused by CF-causing point mutations in the NBD1, which is poorly understood, differential scanning fluorimetry (DSF) and hydrogen deuterium exchange coupled with mass spectrometry (HDX-MS) are the choice of techniques. These established methods can measure the conformational dynamics of the NBD1 globally and at peptide resolution level by monitoring backbone amide HDX, respectively, and will be instrumental to evaluate the mechanism of action of CF mutations and folding correctors that rescue CFTR folding defects via stabilizing the mutant NBD1.

Ligand binding to a remote site thermodynamically corrects the F508del mutation in the human cystic fibrosis transmembrane conductance regulator

Many disease-causing mutations impair protein stability. Here, we explore a thermodynamic strategy to correct the disease-causing F508del mutation in the human cystic fibrosis transmembrane conductance regulator (hCFTR). F508del destabilizes nucleotide-binding domain 1 (hNBD1) in hCFTR relative to an aggregation-prone intermediate. We developed a fluorescence self-quenching assay for compounds that prevent aggregation of hNBD1 by stabilizing its native conformation. Unexpectedly, we found that dTTP and nucleotide analogs with exocyclic methyl groups bind to hNBD1 more strongly than ATP and preserve electrophysiological function of full-length F508del-hCFTR channels at temperatures up to 37 °C. Furthermore, nucleotides that increase open-channel probability, which reflects stabilization of an interdomain interface to hNBD1, thermally protect full-length F508del-hCFTR even when they do not stabilize isolated hNBD1. Therefore, stabilization of hNBD1 itself or of one of its interdomain interfaces by a small molecule indirectly offsets the destabilizing effect of the F508del mutation on full-length hCFTR. These results indicate that high-affinity binding of a small molecule to a remote site can correct a disease-causing mutation. We propose that the strategies described here should be applicable to identifying small molecules to help manage other human diseases caused by mutations that destabilize native protein conformation.

Structure-guided combination therapy to potently improve the function of mutant CFTRs

Available corrector drugs are unable to effectively rescue the folding defects of CFTR-ΔF508 (or CFTR-F508del), the most common disease-causing mutation of the cystic fibrosis transmembrane conductance regulator, a plasma membrane (PM) anion channel, and thus to substantially ameliorate clinical phenotypes of cystic fibrosis (CF). To overcome the corrector efficacy ceiling, here we show that compounds targeting distinct structural defects of CFTR can synergistically rescue mutant expression and function at the PM. High-throughput cell-based screens and mechanistic analysis identified three small-molecule series that target defects at nucleotide-binding domain (NBD1), NBD2 and their membrane-spanning domain (MSD) interfaces. Although individually these compounds marginally improve ΔF508-CFTR folding efficiency, function and stability, their combinations lead to ~50-100% of wild-type-level correction in immortalized and primary human airway epithelia and in mouse nasal epithelia. Likewise, corrector combinations were effective against rare missense mutations in various CFTR domains, probably acting via structural allostery, suggesting a mechanistic framework for their broad application.

ΔF508-CFTR Modulator Screen Based on Cell Surface Targeting of a Chimeric Nucleotide Binding Domain 1 Reporter

The most common cystic fibrosis-causing mutation in the cystic fibrosis transmembrane conductance regulator (CFTR) is deletion of phenylalanine at residue 508 (∆F508). The ∆F508 mutation impairs folding of nucleotide binding domain 1 (NBD1) and interfacial interactions of NBD1 and the membrane spanning domains. Here, we report a domain-targeted screen to identify ∆F508-CFTR modulators that act on NBD1. A biochemical screen for ΔF508-NBD1 cell surface expression was done in Madin-Darby canine kidney cells expressing a chimeric reporter consisting of ΔF508-NBD1, the CD4 transmembrane domain, and an extracellular horseradish peroxidase (HRP) reporter. Using a luminescence readout of HRP activity, the screen was robust with a Z' factor of 0.7. The screening of ~20,000 synthetic small molecules allowed the identification of compounds from four chemical classes that increased ∆F508-NBD1 cell surface expression by up to 4-fold; for comparison, a 12-fold increased cell surface expression was found for a wild-type NBD1 chimera. While the compounds were inactive as correctors of full-length ΔF508-CFTR, several carboxamide-benzothiophenes had potentiator activity with low micromolar EC₅₀. Interestingly, the potentiators did not activate G551D or wild-type CFTR. Our results provide a proof of concept for a cell-based NBD1 domain screen to identify ∆F508-CFTR modulators that target the NBD1 domain.

The implications of CFTR structural studies for cystic fibrosis drug development

Development of Cystic Fibrosis Transmembrane conductance Regulator (CFTR) modulators, targeting the root cause of cystic fibrosis (CF), represents a challenge in the era of personalized medicine, as CFTR mutations lead to a variety of phenotypes, which likely require different, specific treatments. CF drug development is also complicated by the need to preserve the right balance between stability and flexibility, required for optimal function of the CFTR protein. In this review, we highlight how structural data can be exploited in this context to understand the molecular mechanisms of disease-associated mutations, to characterize the mechanisms of action of known modulators and to rationalize the search for novel, specific compounds.

The catalytic mechanism of cyclic GMP-AMP synthase (cGAS) and implications for innate immunity and inhibition

Cyclic GMP-AMP synthase (cGAS) is activated by ds-DNA binding to produce the secondary messenger 2',3'-cGAMP. cGAS is an important control point in the innate immune response; dysregulation of the cGAS pathway is linked to autoimmune diseases while targeted stimulation may be of benefit in immunoncology. We report here the structure of cGAS with dinucleotides and small molecule inhibitors, and kinetic studies of the cGAS mechanism. Our structural work supports the understanding of how ds-DNA activates cGAS, suggesting a site for small molecule binders that may cause cGAS activation at physiological ATP concentrations, and an apparent hotspot for inhibitor binding. Mechanistic studies of cGAS provide the first kinetic constants for 2',3'-cGAMP formation, and interestingly, describe a catalytic mechanism where 2',3'-cGAMP may be a minor product of cGAS compared with linear nucleotides.

Discovery of PF-06928215 as a high affinity inhibitor of cGAS enabled by a novel fluorescence polarization assay

Cyclic GMP-AMP synthase (cGAS) initiates the innate immune system in response to cytosolic dsDNA. After binding and activation from dsDNA, cGAS uses ATP and GTP to synthesize 2', 3' -cGAMP (cGAMP), a cyclic dinucleotide second messenger with mixed 2'-5' and 3'-5' phosphodiester bonds. Inappropriate stimulation of cGAS has been implicated in autoimmune disease such as systemic lupus erythematosus, thus inhibition of cGAS may be of therapeutic benefit in some diseases; however, the size and polarity of the cGAS active site makes it a challenging target for the development of conventional substrate-competitive inhibitors. We report here the development of a high affinity (KD = 200 nM) inhibitor from a low affinity fragment hit with supporting biochemical and structural data showing these molecules bind to the cGAS active site. We also report a new high throughput cGAS fluorescence polarization (FP)-based assay to enable the rapid identification and optimization of cGAS inhibitors. This FP assay uses Cy5-labelled cGAMP in combination with a novel high affinity monoclonal antibody that specifically recognizes cGAMP with no cross reactivity to cAMP, cGMP, ATP, or GTP. Given its role in the innate immune response, cGAS is a promising therapeutic target for autoinflammatory disease. Our results demonstrate its druggability, provide a high affinity tool compound, and establish a high throughput assay for the identification of next generation cGAS inhibitors.

Molecular modelling and molecular dynamics of CFTR

The cystic fibrosis transmembrane conductance regulator (CFTR) protein is a member of the ATP-binding cassette (ABC) transporter superfamily that functions as an ATP-gated channel. Considerable progress has been made over the last years in the understanding of the molecular basis of the CFTR functions, as well as dysfunctions causing the common genetic disease cystic fibrosis (CF). This review provides a global overview of the theoretical studies that have been performed so far, especially molecular modelling and molecular dynamics (MD) simulations. A special emphasis is placed on the CFTR-specific evolution of an ABC transporter framework towards a channel function, as well as on the understanding of the effects of disease-causing mutations and their specific modulation. This in silico work should help structure-based drug discovery and design, with a view to develop CFTR-specific pharmacotherapeutic approaches for the treatment of CF in the context of precision medicine.