Stabilization of a nucleotide-binding domain of the cystic fibrosis transmembrane conductance regulator yields insight into disease-causing mutations

The Folding Pathway of ABC Transporter CFTR: Effective and Robust

De novo protein folding into a native three-dimensional structure is indispensable for biological function, is instructed by its amino acid sequence, and occurs along a vectorial trajectory. The human proteome contains thousands of membrane-spanning proteins, whose biosynthesis begins on endoplasmic reticulum-associated ribosomes. Nearly half of all membrane proteins traverse the membrane more than once, including therapeutically important protein families such as solute carriers, G-protein-coupled receptors, and ABC transporters. These mediate a variety of functions like signal transduction and solute transport and are often of vital importance for cell function and tissue homeostasis. Missense mutations in multispan membrane proteins can lead to misfolding and cause disease; an example is the ABC transporter Cystic Fibrosis Transmembrane Conductance Regulator (CFTR). Even though our understanding of multispan membrane-protein folding still is rather rudimental, the cumulative knowledge of 20 years of basic research on CFTR folding has led to development of drugs that modulate the misfolded protein. This has provided the prospect of a life without CF to the vast majority of patients. In this review we describe our understanding of the folding pathway of CFTR in cells, which is modular and tolerates many defects, making it effective and robust. We address how modulator drugs affect folding and function of CFTR, and distinguish protein stability from its folding process. Since the domain architecture of (mammalian) ABC transporters are highly conserved, we anticipate that the insights we discuss here for folding of CFTR may lay the groundwork for understanding the general rules of ABC-transporter folding.

Transmembrane Helices 7 and 8 Confer Aggregation Sensitivity to the Cystic Fibrosis Transmembrane Conductance Regulator

The Cystic Fibrosis Transmembrane Conductance Regulator (CFTR) is a large multi-spanning membrane protein that is susceptible to misfolding and aggregation. We have identified here the region responsible for this instability. Temperature-induced aggregation of C-terminally truncated versions of CFTR demonstrated that all truncations up to the second transmembrane domain (TMD2), including the R region, largely resisted aggregation. Limited proteolysis identified a folded structure that was prone to aggregation and consisted of TMD2 and at least part of the Regulatory Region R. Only when both TM7 (TransMembrane helix 7) and TM8 were present, TMD2 fragments became as aggregation-sensitive as wild-type CFTR, in line with increased thermo-instability of late CFTR nascent chains and in silico prediction of aggregation propensity. In accord, isolated TMD2 was degraded faster in cells than isolated TMD1. We conclude that TMD2 extended at its N-terminus with part of the R region forms a protease-resistant structure that induces heat instability in CFTR and may be responsible for its limited intracellular stability.

ABC-transporter CFTR folds with high fidelity through a modular, stepwise pathway

The question how proteins fold is especially pointed for large multi-domain, multi-spanning membrane proteins with complex topologies. We have uncovered the sequence of events that encompass proper folding of the ABC transporter CFTR in live cells by combining kinetic radiolabeling with protease-susceptibility assays. We found that CFTR folds in two clearly distinct stages. The first, co-translational, stage involves folding of the 2 transmembrane domains TMD1 and TMD2, plus one nucleotide-binding domain, NBD1. The second stage is a simultaneous, post-translational increase in protease resistance for both TMDs and NBD2, caused by assembly of these domains onto NBD1. Our assays probe every 2-3 residues (on average) in CFTR. This in-depth analysis at amino-acid level allows detailed analysis of domain folding and importantly also the next level: assembly of the domains into native, folded CFTR. Defects and changes brought about by medicines, chaperones, or mutations also are amenable to analysis. We here show that the well-known disease-causing mutation F508del, which established cystic fibrosis as protein-folding disease, caused co-translational misfolding of NBD1 but not TMD1 nor TMD2 in stage 1, leading to absence of stage-2 folding. Corrector drugs rescued stage 2 without rescuing NBD1. Likewise, the DxD motif in NBD1 that was identified to be required for export of CFTR from the ER we found to be required already upstream of export as CFTR mutated in this motif phenocopies F508del CFTR. The highly modular and stepwise folding process of such a large, complex protein explains the relatively high fidelity and correctability of its folding.

Spatial covariance analysis reveals the residue-by-residue thermodynamic contribution of variation to the CFTR fold

Although the impact of genome variation on the thermodynamic properties of function on the protein fold has been studied in vitro, it remains a challenge to assign these relationships across the entire polypeptide sequence in vivo. Using the Gaussian process regression based principle of Spatial CoVariance, we globally assign on a residue-by-residue basis the biological thermodynamic properties that contribute to the functional fold of CFTR in the cell. We demonstrate the existence of a thermodynamically sensitive region of the CFTR fold involving the interface between NBD1 and ICL4 that contributes to its export from endoplasmic reticulum. At the cell surface a new set of residues contribute uniquely to the management of channel function. These results support a general ‘quality assurance’ view of global protein fold management as an SCV principle describing the differential pre- and post-ER residue interactions contributing to compartmentalization of the energetics of the protein fold for function. Our results set the stage for future analyses of the quality systems managing protein sequence-to-function-to-structure broadly encompassing genome design leading to protein function in complex cellular relationships responsible for diversity and fitness in biology in response to the environment.

Spatial covariance analysis shows how each residue in the critical fold region of chloride channel CFTR, defective in cystic fibrosis patients, contributes to CFTR’s export from the endoplasmic reticulum and function in the cell.

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.

Factoring in the Complexity of the Cystic Fibrosis Lung to Understand Aspergillus fumigatus and Pseudomonas aeruginosa Interactions

Pseudomonas aeruginosa has long been established as the most prevalent respiratory pathogen in Cystic Fibrosis (CF) patients, with opportunistic infection causing profound morbidity and mortality. Recently, Aspergillus fumigatus has also been recognised as a key contributor to CF lung deterioration, being consistently associated with decreased lung function and worsened prognosis in these patients. As clinical evidence for the common occurrence of combined infection with these two pathogens increases, research into the mechanism and consequences of their interaction is becoming more relevant. Clinical evidence suggests a synergistic effect of combined infection, which translates into a poorer prognosis for the patients. In vitro results from the laboratory have identified a variety of possible synergistic and antagonistic interactions between A. fumigatus and P. aeruginosa. Here, we present a comprehensive overview of the complex environment of the CF lung and discuss how it needs to be considered to determine the exact molecular interactions that A. fumigatus and P. aeruginosa undergo during combined infection and their effects on the host.

CFTR structure, stability, function and regulation

Cystic fibrosis transmembrane conductance regulator (CFTR) is a unique member of the ATP-binding cassette family of proteins because it has evolved into a channel. Mutations in CFTR cause cystic fibrosis, the most common genetic disease in people of European origin. The F508del mutation is found in about 90% of patients and here we present data that suggest its main effect is on CFTR stability rather than on the three-dimensional (3D) folded state. A survey of recent cryo-electron microscopy studies was carried out and this highlighted differences in terms of CFTR conformation despite similarities in experimental conditions. We further studied CFTR structure under various phosphorylation states and with the CFTR-interacting protein NHERF1. The coexistence of outward-facing and inward-facing conformations under a range of experimental conditions was suggested from these data. These results are discussed in terms of structural models for channel gating, and favour the model where the mostly disordered regulatory-region of the protein acts as a channel plug.

Disease-relevant mutations alter amino acid co-evolution networks in the second nucleotide binding domain of CFTR

Cystic Fibrosis (CF) is an inherited disease caused by mutations in the cystic fibrosis transmembrane conductance regulator (CFTR) ion channel. Mutations in CFTR cause impaired chloride ion transport in the epithelial tissues of patients leading to cardiopulmonary decline and pancreatic insufficiency in the most severely affected patients. CFTR is composed of twelve membrane-spanning domains, two nucleotide-binding domains (NBDs), and a regulatory domain. The most common mutation in CFTR is a deletion of phenylalanine at position 508 (ΔF508) in NBD1. Previous research has primarily concentrated on the structure and dynamics of the NBD1 domain; However numerous pathological mutations have also been found in the lesser-studied NBD2 domain. We have investigated the amino acid co-evolved network of interactions in NBD2, and the changes that occur in that network upon the introduction of CF and CF-related mutations (S1251N(T), S1235R, D1270N, N1303K(T)). Extensive coupling between the α- and β-subdomains were identified with residues in, or near Walker A, Walker B, H-loop and C-loop motifs. Alterations in the predicted residue network varied from moderate for the S1251T perturbation to more severe for N1303T. The S1235R and D1270N networks varied greatly compared to the wildtype, but these CF mutations only affect ion transport preference and do not severely disrupt CFTR function, suggesting dynamic flexibility in the network of interactions in NBD2. Our results also suggest that inappropriate interactions between the β-subdomain and Q-loop could be detrimental. We also identified mutations predicted to stabilize the NBD2 residue network upon introduction of the CF and CF-related mutations, and these predicted mutations are scored as benign by the MUTPRED2 algorithm. Our results suggest the level of disruption of the co-evolution predictions of the amino acid networks in NBD2 does not have a straightforward correlation with the severity of the CF phenotypes observed.

Hsp104 facilitates the endoplasmic-reticulum-associated degradation of disease-associated and aggregation-prone substrates

Misfolded proteins in the endoplasmic reticulum (ER) are selected for ER-associated degradation (ERAD). More than 60 disease-associated proteins are substrates for the ERAD pathway due to the presence of missense or nonsense mutations. In yeast, the Hsp104 molecular chaperone disaggregates detergent-insoluble ERAD substrates, but the spectrum of disease-associated ERAD substrates that may be aggregation prone is unknown. To determine if Hsp104 recognizes aggregation-prone ERAD substrates associated with human diseases, we developed yeast expression systems for a hydrophobic lipid-binding protein, apolipoprotein B (ApoB), along with a chimeric protein harboring a nucleotide-binding domain from the cystic fibrosis transmembrane conductance regulator (CFTR) into which disease-causing mutations were introduced. We discovered that Hsp104 facilitates the degradation of ER-associated ApoB as well as a truncated CFTR chimera in which a premature stop codon corresponds to a disease-causing mutation. Chimeras containing a wild-type version of the CFTR domain or a different mutation were stable and thus Hsp104 independent. We also discovered that the detergent solubility of the unstable chimera was lower than the stable chimeras, and Hsp104 helped retrotranslocate the unstable chimera from the ER, consistent with disaggregase activity. To determine why the truncated chimera was unstable, we next performed molecular dynamics simulations and noted significant unraveling of the CFTR nucleotide-binding domain. Because human cells lack Hsp104, these data indicate that an alternate disaggregase or mechanism facilitates the removal of aggregation-prone, disease-causing ERAD substrates in their native environments.

Different SUMO paralogues determine the fate of wild-type and mutant CFTRs: biogenesis versus degradation

A pathway for cystic fibrosis transmembrane conductance regulator (CFTR) degradation is initiated by Hsp27, which cooperates with Ubc9 and binds to the common F508del mutant to modify it with SUMO-2/3. These SUMO paralogues form polychains, which are recognized by the ubiquitin ligase, RNF4, for proteosomal degradation. Here, protein array analysis identified the SUMO E3, protein inhibitor of activated STAT 4 (PIAS4), which increased wild-type (WT) and F508del CFTR biogenesis in CFBE airway cells. PIAS4 increased immature CFTR threefold and doubled expression of mature CFTR, detected by biochemical and functional assays. In cycloheximide chase assays, PIAS4 slowed immature F508del degradation threefold and stabilized mature WT CFTR at the plasma membrance. PIAS4 knockdown reduced WT and F508del CFTR expression by 40–50%, suggesting a physiological role in CFTR biogenesis. PIAS4 modified F508del CFTR with SUMO-1 in vivo and reduced its conjugation to SUMO-2/3. These SUMO paralogue-specific effects of PIAS4 were reproduced in vitro using purified F508del nucleotide-binding domain 1 and SUMOylation reaction components. PIAS4 reduced endogenous ubiquitin conjugation to F508del CFTR by ∼50% and blocked the impact of RNF4 on mutant CFTR disposal. These findings indicate that different SUMO paralogues determine the fates of WT and mutant CFTRs, and they suggest that a paralogue switch during biogenesis can direct these proteins to different outcomes: biogenesis versus degradation.

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.

ABC transportome inventory of human pathogenic yeast Candida glabrata: Phylogenetic and expression analysis

ATP-binding cassette (ABC) is one of the two major superfamilies of transporters present across the evolutionary scale. ABC superfamily members came to prominence due to their ability to extrude broad spectrum of substrates and to confer multi drug resistance (MDR). Overexpression of some ABC transporters in clinical isolates of Candida species was attributed to the development of MDR phenotypes. Among Candida species, Candida glabrata is an emerging drug resistant species in human fungal infections. A comprehensive analysis of such proteins in C. glabrata is required to untangle their role not only in MDR but also in other biological processes. Bioinformatic analysis of proteins encoded by genome of human pathogenic yeast C. glabrata identified 25 putative ABC protein coding genes. On the basis of phylogenetic analysis, domain organization and nomenclature adopted by the Human Genome Organization (HUGO) scheme, these proteins were categorized into six subfamilies such as Pleiotropic Drug Resistance (PDR)/ABCG, Multi Drug Resistance (MDR)/ABCB, Multi Drug Resistance associated Protein (MRP)/ABCC, Adrenoleukodystrophy protein (ALDp)/ABCD, RNase L Inhibitor (RLI)/ABCE and Elongation Factor 3 (EF3)/ABCF. Among these, only 18 ABC proteins contained transmembrane domains (TMDs) and were grouped as membrane proteins, predominantly belonging to PDR, MDR, MRP, and ALDp subfamilies. A comparative phylogenetic analysis of these ABC proteins with other yeast species revealed their orthologous relationship and pointed towards their conserved functions. Quantitative real time PCR (qRT-PCR) analysis of putative membrane localized ABC protein encoding genes of C. glabrata confirmed their basal expression and showed variable transcriptional response towards antimycotic drugs. This study presents first comprehensive overview of ABC superfamily proteins of a human fungal pathogen C. glabrata, which is expected to provide an important platform for in depth analysis of their physiological relevance in cellular processes and drug resistance.

Structural analysis reveals pathomechanisms associated with pseudoxanthoma elasticum-causing mutations in the ABCC6 transporter

Mutations in ABC subfamily C member 6 (ABCC6) transporter are associated with pseudoxanthoma elasticum (PXE), a disease resulting in ectopic mineralization and affecting multiple tissues. A growing number of mutations have been identified in individuals with PXE. For most of these variants, no mechanistic information is available regarding their role in normal and pathophysiologies. To assess how PXE-associated mutations alter ABCC6 biosynthesis and structure, we biophysically and biochemically evaluated the N-terminal nucleotide-binding domain. A high-resolution X-ray structure of nucleotide-binding domain 1 (NBD1) of human ABCC6 was obtained at 2.3 Å that provided a template on which to evaluate PXE-causing mutations. Biochemical analysis of mutations in this domain indicated that multiple PXE-causing mutations altered its structural properties. Analyses of the full-length protein revealed a strong correlation between the alterations in NBD properties and the processing and expression of ABCC6. These results suggest that a significant fraction of PXE-associated mutations located in NBD1 causes changes in its structural properties and that these mutation-induced alterations directly affect the maturation of the full-length ABCC6 protein.

Structural stability of purified human CFTR is systematically improved by mutations in nucleotide binding domain 1

The Cystic Fibrosis Transmembrane Conductance Regulator (CFTR) is an ABC transporter containing two transmembrane domains forming a chloride ion channel, and two nucleotide binding domains (NBD1 and NBD2). CFTR has presented a formidable challenge to obtain monodisperse, biophysically stable protein. Here we report a comprehensive study comparing effects of single and multiple NBD1 mutations on stability of both the NBD1 domain alone and on purified full length human CFTR. Single mutations S492P, A534P, I539T acted additively, and when combined with M470V, S495P, and R555K cumulatively yielded an NBD1 with highly improved structural stability. Strategic combinations of these mutations strongly stabilized the domain to attain a calorimetric T_m > 70 °C. Replica exchange molecular dynamics simulations on the most stable 6SS-NBD1 variant implicated fluctuations, electrostatic interactions and side chain packing as potential contributors to improved stability. Progressive stabilization of NBD1 directly correlated with enhanced structural stability of full-length CFTR protein. Thermal unfolding of the stabilized CFTR mutants, monitored by changes in intrinsic fluorescence, demonstrated that Tm could be shifted as high as 67.4 °C in 6SS-CFTR, more than 20 °C higher than wild-type. H1402S, an NBD2 mutation, conferred CFTR with additional thermal stability, possibly by stabilizing an NBD-dimerized conformation. CFTR variants with NBD1-stabilizing mutations were expressed at the cell surface in mammalian cells, exhibited ATPase and channel activity, and retained these functions to higher temperatures. The capability to produce enzymatically active CFTR with improved structural stability amenable to biophysical and structural studies will advance mechanistic investigations and future cystic fibrosis drug development.

CFTR structure

Structural studies of the cystic fibrosis transmembrane conductance regulator (CFTR) protein are critical to understand molecular mechanisms involved in gating of the apical anion channel as well as the way in which the gating is regulated, especially by the regulatory region (R region). They are also instrumental for understanding the root cause of cystic fibrosis (CF) and supporting the development of therapeutic strategies. In this short review, we summarize recent progress in the knowledge of the CFTR 3D structure and briefly discuss implications for CF drug development.