Novel gain-of-function mutants identify a critical region within CFTR membrane-spanning domain 2 controlling cAMP-dependent and ATP-independent channel activation.

CFTR is an anion channel that has evolved from the mold of an ABC transporter. It possesses specific structural features, including a lateral portal between the cytoplasmic extensions of its transmembrane helices TM4 and TM6. This TM4-TM6 portal is lined by basic residues attracting anions from the cytosol towards the intracellular vestibule. Even though a symmetric, open portal is not observed at the level of the TM10/TM12 interface, basic amino acids are also present at this level, exposed to solvent in the vicinity of the regulatory R region, whose phosphorylation enables channel activation. Here, using all-atom molecular dynamics simulations in combination with functional and biochemical assays, we investigate the importance of these basic amino acids (R1158 and R1030), and of a neighboring aromatic amino acid (W846) in the regulation of CFTR activity. Results indicate that mutation of these amino acids globally increased channel activity and enabled channel opening by potentiators without the need to elevate cAMP levels. These effects (i) were observed even when the binding site of the potentiator VX-770 was mutated, revealing a probable independent mechanism, and (ii) were additive to one gain-of-function mutant within the selectivity filter. Taken together, our results indicate that the region of the membrane-spanning domain 2 (MSD2), symmetric to the lateral portal located between MSD1 TM4 and TM6, is a novel critical actor of CFTR regulation.

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.

Acting on the CFTR Membrane-Spanning Domains Interface Rescues Some Misfolded Mutants.

ABC transporters are large membrane proteins sharing a complex architecture, which comprises two nucleotide-binding domains (NBDs) and two membrane-spanning domains (MSDs). These domains are susceptible to mutations affecting their folding and assembly. In the CFTR (ABCC7) protein, a groove has been highlighted in the MSD1 at the level of the membrane inner leaflet, containing both multiple mutations affecting folding and a binding site for pharmaco-chaperones that stabilize this region. This groove is also present in ABCB proteins, however it is covered by a short elbow helix, while in ABCC proteins it remains unprotected, due to a lower position of the elbow helix in the presence of the ABCC-specific lasso motif. Here, we identified a MSD1 second-site mutation located in the vicinity of the CFTR MSD1 groove that partially rescued the folding defect of cystic fibrosis causing mutations located within MSD1, while having no effect on the most frequent mutation, F508del, located within NBD1. A model of the mutated protein 3D structure suggests additional interaction between MSD1 and MSD2, strengthening the assembly at the level of the MSD intracellular loops. Altogether, these results provide insightful information in understanding key features of the folding and function of the CFTR protein in particular, and more generally, of type IV ABC transporters.

Combining theoretical and experimental data to decipher CFTR 3D structures and functions.

Cryo-electron microscopy (cryo-EM) has recently provided invaluable experimental data about the full-length cystic fibrosis transmembrane conductance regulator (CFTR) 3D structure. However, this experimental information deals with inactive states of the channel, either in an apo, quiescent conformation, in which nucleotide-binding domains (NBDs) are widely separated or in an ATP-bound, yet closed conformation. Here, we show that 3D structure models of the open and closed forms of the channel, now further supported by metadynamics simulations and by comparison with the cryo-EM data, could be used to gain some insights into critical features of the conformational transition toward active CFTR forms. These critical elements lie within membrane-spanning domains but also within NBD1 and the N-terminal extension, in which conformational plasticity is predicted to occur to help the interaction with filamin, one of the CFTR cellular partners.

Order in Disorder as Observed by the "Hydrophobic Cluster Analysis" of Protein Sequences.

Hydrophobic cluster analysis (HCA) is an original approach for protein sequence analysis, which provides access to the foldable repertoire of the protein universe, including yet unannotated protein segments ("dark proteome"). Foldable segments correspond to ordered regions, as well as to intrinsically disordered regions (IDRs) undergoing disorder to order transitions. In this review, how HCA can be used to give insight into this last category of foldable segments is illustrated, with examples matching known 3D structures. After reviewing the HCA principles, examples of short foldable segments are given, which often contain short linear motifs, typically matching hydrophobic clusters. These segments become ordered upon contact with partners, with secondary structure preferences generally corresponding to those observed in the 3D structures within the complexes. Such small foldable segments are sometimes larger than the segments of known 3D structures, including flanking hydrophobic clusters that may be critical for interaction specificity or regulation, as well as intervening sequences allowing fuzziness. Cases of larger conditionally disordered domains are also presented, with lower density in hydrophobic clusters than well-folded globular domains or with exposed hydrophobic patches, which are stabilized by interaction with partners.

Cis variants identified in F508del complex alleles modulate CFTR channel rescue by small molecules.

Molecules correcting the trafficking (correctors) and gating defects (potentiators) of the cystic fibrosis causing mutation c.1521_1523delCTT (p.Phe508del) begin to be a useful treatment for CF patients bearing p.Phe508del. This mutation has been identified in different genetic contexts, alone or in combination with variants in cis. Until now, 21 exonic variants in cis of p.Phe508del have been identified, albeit at a low frequency. The aim of this study was to evaluate their impact on the efficacy of CFTR-directed corrector/potentiator therapy (Orkambi). The analysis by minigene showed that two out of 15 cis variants tested increased exon skipping (c.609C > T and c.2770G > A). Four cis variants were studied functionally in the absence of p.Phe508del, one of which was found to be deleterious for protein maturation c.1399C > T (p.Leu467Phe). In the presence of p.Phe508del, this variant was the only to prevent the response to Orkambi treatment. This study showed that some patients carrying p.Phe508del complex alleles are predicted to poorly respond to corrector/potentiator treatments. Our results underline the importance to validate treatment efficacy in the context of complex alleles.

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.

The respective roles of polar/nonpolar binary patterns and amino acid composition in protein regular secondary structures explored exhaustively using hydrophobic cluster analysis.

Several studies have highlighted the leading role of the sequence periodicity of polar and nonpolar amino acids (binary patterns) in the formation of regular secondary structures (RSS). However, these were based on the analysis of only a few simple cases, with no direct mean to correlate binary patterns with the limits of RSS. Here, HCA-derived hydrophobic clusters (HC) which are conditioned binary patterns whose positions fit well those of RSS, were considered. All the HC types, defined by unique binary patterns, which were commonly observed in three-dimensional (3D) structures of globular domains, were analyzed. The 180 HC types with preferences for either α-helices or ß-strands distinctly contain basic binary units typical of these RSS. Therefore a general trend supporting the "binary pattern preference" assumption was observed. HC for which observed RSS are in disagreement with their expected behavior (discordant HC) were also examined. They were separated in HC types with moderate preferences for RSS, having "weak" binary patterns and versatile RSS and HC types with high preferences for RSS, having "strong" binary patterns and then displaying nonpolar amino acids at the protein surface. It was shown that in both cases, discordant HC could be distinguished from concordant ones by well-differentiated amino acid compositions. The obtained results could, thus, help to complement the currently available methods for the accurate prediction of secondary structures in proteins from the only information of a single amino acid sequence. This can be especially useful for characterizing orphan sequences and for assisting protein engineering and design.

Full-open and closed CFTR channels, with lateral tunnels from the cytoplasm and an alternative position of the F508 region, as revealed by molecular dynamics.

In absence of experimental 3D structures, several homology models, based on ABC exporter 3D structures, have provided significant insights into the molecular mechanisms underlying the function of the cystic fibrosis transmembrane conductance regulator (CFTR) protein, a chloride channel whose defects are associated with cystic fibrosis (CF). Until now, these models, however, did not furnished much insights into the continuous way that ions could follow from the cytosol to the extracellular milieu in the open form of the channel. Here, we have built a refined model of CFTR, based on the outward-facing Sav1866 experimental 3D structure and integrating the evolutionary and structural information available today. Molecular dynamics simulations revealed significant conformational changes, resulting in a full-open channel, accessible from the cytosol through lateral tunnels displayed in the long intracellular loops (ICLs). At the same time, the region of nucleotide-binding domain 1 in contact with one of the ICLs and carrying amino acid F508, the deletion of which is the most common CF-causing mutation, was found to adopt an alternative but stable position. Then, in a second step, this first stable full-open conformation evolved toward another stable state, in which only a limited displacement of the upper part of the transmembrane helices leads to a closure of the channel, in a conformation very close to that adopted by the Atm1 ABC exporter, in an inward-facing conformation. These models, supported by experimental data, provide significant new insights into the CFTR structure-function relationships and into the possible impact of CF-causing mutations.

Structure-function analysis of the human ferroportin iron exporter (SLC40A1): effect of hemochromatosis type 4 disease mutations and identification of critical residues.

Ferroportin (SLC40A1) is the only known iron exporter in mammals and is considered a key coordinator of the iron balance between intracellular and systemic iron homeostasis. However, the structural organization of ferroportin in the lipid bilayer remains controversial and very little is known about the mechanism underlying iron egress. In the present study, we have developed an approach based on comparative modeling, which has led to the construction of a model of the three-dimensional (3D) structure of ferroportin by homology to the crystal structure of a Major Facilitator Superfamily member (EmrD). This model predicts atomic details for the organization of ferroportin transmembrane helices and is in agreement with our current understanding of the ferroportin function and its interaction with hepcidin. Using in vitro experiments, we demonstrate that this model can be used to identify novel critical amino acids. In particular, we show that the tryptophan residue 42 (p.Trp42), which is localized within the extracellular end of the ferroportin pore, is likely involved in both the iron export function and in the mechanism of inhibition by hepcidin. Thus, our 3D model provides a new perspective for understanding the molecular basis of ferroportin functions and dysfunctions.

The Hedgehog-inducible ubiquitin ligase subunit WSB-1 modulates thyroid hormone activation and PTHrP secretion in the developing growth plate.

WSB-1 is a SOCS-box-containing WD-40 protein of unknown function that is induced by Hedgehog signalling in embryonic structures during chicken development. Here we show that WSB-1 is part of an E3 ubiquitin ligase for the thyroid-hormone-activating type 2 iodothyronine deiodinase (D2). The WD-40 propeller of WSB-1 recognizes an 18-amino-acid loop in D2 that confers metabolic instability, whereas the SOCS-box domain mediates its interaction with a ubiquitinating catalytic core complex, modelled as Elongin BC-Cul5-Rbx1 (ECS(WSB-1)). In the developing tibial growth plate, Hedgehog-stimulated D2 ubiquitination via ECS(WSB-1) induces parathyroid hormone-related peptide (PTHrP), thereby regulating chondrocyte differentiation. Thus, ECS(WSB-1) mediates a mechanism by which 'systemic' thyroid hormone can effect local control of the Hedgehog-PTHrP negative feedback loop and thus skeletogenesis.

Digging into the 3D Structure Predictions of AlphaFold2 with Low Confidence: Disorder and Beyond.

AlphaFold2 (AF2) has created a breakthrough in biology by providing three-dimensional structure models for whole-proteome sequences, with unprecedented levels of accuracy. In addition, the AF2 pLDDT score, related to the model confidence, has been shown to provide a good measure of residue-wise disorder. Here, we combined AF2 predictions with pyHCA, a tool we previously developed to identify foldable segments and estimate their order/disorder ratio, from a single protein sequence. We focused our analysis on the AF2 predictions available for 21 reference proteomes (AFDB v1), in particular on their long foldable segments (>30 amino acids) that exhibit characteristics of soluble domains, as estimated by pyHCA. Among these segments, we provided a global analysis of those with very low pLDDT values along their entire length and compared their characteristics to those of segments with very high pLDDT values. We highlighted cases containing conditional order, as well as cases that could form well-folded structures but escape the AF2 prediction due to a shallow multiple sequence alignment and/or undocumented structure or fold. AF2 and pyHCA can therefore be advantageously combined to unravel cryptic structural features in whole proteomes and to refine predictions for different flavors of disorder.

C terminus of nucleotide binding domain 1 contains critical features for cystic fibrosis transmembrane conductance regulator trafficking and activation.

The cystic fibrosis transmembrane conductance regulator (CFTR) is a Cl(-) channel physiologically important in fluid-transporting epithelia and pathologically relevant in several human diseases. Here, we show that mutations in the C terminus of the first nucleotide binding domain comprising the latest beta strands (beta(c)5 and beta(c)6) influence the trafficking, channel activity, and pharmacology of CFTR. We mutated CFTR amino acids located in the beta(c)5-beta(c)6 hairpin, within the beta(c)5 strand (H620Q), within the beta-turn linking the two beta strands (E621G, G622D), as well as within (S623A, S624A) and at the extremity (G628R) of the beta(c)6 strand. Functional analysis reveals that the current density was largely reduced for G622D and G628R channels compared with wt CFTR, similar for E621G and S624A, but increased for H620Q and S623A. For G622D and G628R, the abnormal activity is likely due to a defective maturation process, as assessed by the augmented activity and mature C-band observed in the presence of the trafficking corrector miglustat. In addition, in presence of the CFTR activator benzo[c]quinolizinium, the CFTR current density compared with that of wt CFTR was abolished for G622D and G628R channels, but similar for H620Q, S623A, and S624A or slightly increased for E621G. Finally, G622D and G628R were activated by the CFTR agonists genistein, RP-107, and isobutylmethylxanthine. Our results identify the C terminus of the CFTR first nucleotide binding domain as an important molecular site for the trafficking of CFTR protein, for the control of CFTR channel gating, and for the pharmacological effect of a dual activity agent.

Delineation of the Xrcc4-interacting region in the globular head domain of cernunnos/XLF.

In mammals, the majority of DNA double-strand breaks are processed by the nonhomologous end-joining (NHEJ) pathway, composed of seven factors: Ku70, Ku80, DNA-PKcs, Artemis, Xrcc4 (X4), DNA-ligase IV (L4), and Cernunnos/XLF. Cernunnos is part of the ligation complex, constituted by X4 and L4. To improve our knowledge on the structure and function of Cernunnos, we performed a systematic mutagenesis study on positions selected from an analysis of the recent three-dimensional structures of this factor. Ten of 27 screened mutants were nonfunctional in several DNA repair assays. Outside amino acids critical for the expression and stability of Cernunnos, we identified three amino acids (Arg(64), Leu(65), and Leu(115)) essential for the interaction with X4 and the proper function of Cernunnos. Docking the crystal structures of the two factors further validated this probable interaction surface of Cernunnos with X4.

LOTUS, a new domain associated with small RNA pathways in the germline.

We describe here LOTUS, a hitherto uncharacterized small globular domain, which was identified using sensitive sequence profile analysis. The LOTUS domain is found in germline-specific proteins that are present in the nuage/polar granules of germ cells. TDRD5 and TDRD7, two mammalian members of the germline Tudor group, possess three copies of the LOTUS domain in their extreme N-termini. The Tudor domains of these proteins bind symmetric dimethyl arginines present on the germ cell-specific Piwi proteins, which form a particular clade of Argonaute proteins. Piwi proteins interact with a specific class of non-coding RNAs [piwi-interacting RNAs (piRNAs)] and play a key role in the repression (silencing) of transposons and possibly other germline-specific functions. A LOTUS domain is also present in the Oskar protein, a critical component of the pole plasm in the Drosophila oocyte, which is required for germ cell formation. LOTUS domains are found in various proteins from metazoans and plants, are often associated with RNA-specific modules and are likely to adopt a winged helix fold. This suggests a germline-specific role in the mRNA localization and/or translation or a specific function toward piRNAs.

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.

A histidine in the beta-CASP domain of Artemis is critical for its full in vitro and in vivo functions.

Artemis is a key factor of the nonhomologous end-joining (NHEJ) pathway, which is critical for DNA double-strand break (DSB) repair in eukaryotic cells. It belongs to the beta-CASP family of nucleases, forming a distinct group within the metallo-beta-lactamase superfamily. Proteins of this group are specific for nucleic acids and contain an original domain, the beta-CASP domain, which serves as a cap covering the active site displayed by the metallo-beta-lactamase domain.Here, we have identified in the highly divergent sequences of the beta-CASP domains from DNA-specific nucleases two conserved residues (Artemis E213 and H254), which are not present in RNA-specific enzymes, and shown that H254 plays a key role in the Artemis function, as it is critical for its full activity in vitro. Moreover, inherited mutation of H254 results in radiosensitive severe combined immune deficiency (RS-SCID) in humans. This residue might play a key role in specificity towards DNA, if not directly in zinc binding.

A concanavalin A-like lectin domain in the CHS1/LYST protein, shared by members of the BEACH family.

CHS1/LYST, the causative protein of the Chediak-Higashi syndrome (CHS), belongs to the BEACH (named after BEige And Chediak-Higashi) family, which includes various large proteins sharing the same C-terminal domain architecture [a PH (Pleckstrin homology)-BEACH domain followed by WD repeats). Members of the BEACH family are generally defined as vesicle-trafficking regulatory proteins, but their functions remain to be determined at the molecular level. Here, using a panel of sensitive methods of sequence analysis, we show that the N-terminal regions of BEACH proteins contain an as yet not described domain, which shares striking similarities with clostridial neurotoxins and defines a novel family within the concanavalin A (ConA)-like lectin superfamily. These results suggest that the BEACH ConA-like lectin domain could be involved in oligosaccharide binding associated with protein traffic and sorting along the secretory pathway, especially in relation with components of the vesicle fusion machinery.

Ubiquitination-induced conformational change within the deiodinase dimer is a switch regulating enzyme activity.

Ubiquitination is a critical posttranslational regulator of protein stability and/or subcellular localization. Here we show that ubiquitination can also regulate proteins by transiently inactivating enzymatic function through conformational change in a dimeric enzyme, which can be reversed upon deubiquitination. Our model system is the thyroid hormone-activating type 2 deiodinase (D2), an endoplasmic reticulum-resident type 1 integral membrane enzyme. D2 exists as a homodimer maintained by interacting surfaces at its transmembrane and globular cytosolic domains. The D2 dimer associates with the Hedgehog-inducible ubiquitin ligase WSB-1, the ubiquitin conjugase UBC-7, and VDU-1, a D2-specific deubiquitinase. Upon binding of T4, its natural substrate, D2 is ubiquitinated, which inactivates the enzyme by interfering with D2's globular interacting surfaces that are critical for dimerization and catalytic activity. This state of transient inactivity and change in dimer conformation persists until deubiquitination. The continuous association of D2 with this regulatory protein complex supports rapid cycles of deiodination, conjugation to ubiquitin, and enzyme reactivation by deubiquitination, allowing tight control of thyroid hormone action.

Molecular models of the open and closed states of the whole human CFTR protein.

Cystic fibrosis transmembrane conductance regulator (CFTR), involved in cystic fibrosis (CF), is a chloride channel belonging to the ATP-binding cassette (ABC) superfamily. Using the experimental structure of Sav1866 as template, we previously modeled the human CFTR structure, including membrane-spanning domains (MSD) and nucleotide-binding domains (NBD), in an outward-facing conformation (open channel state). Here, we constructed a model of the CFTR inward-facing conformation (closed channel) on the basis of the recent corrected structures of MsbA and compared the structural features of those two states of the channel. Interestingly, the MSD:NBD coupling interfaces including F508 (DeltaF508 being the most common CF mutation) are mainly left unchanged. This prediction, completed by the modeling of the regulatory R domain, is supported by experimental data and provides a molecular basis for a better understanding of the functioning of CFTR, especially of the structural features that make CFTR the unique channel among the ABC transporters.