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

[This corrects the article DOI: 10.1371/journal.pone.0227668.].

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.

SubVis: an interactive R package for exploring the effects of multiple substitution matrices on pairwise sequence alignment.

Understanding how proteins mutate is critical to solving a host of biological problems. Mutations occur when an amino acid is substituted for another in a protein sequence. The set of likelihoods for amino acid substitutions is stored in a matrix and input to alignment algorithms. The quality of the resulting alignment is used to assess the similarity of two or more sequences and can vary according to assumptions modeled by the substitution matrix. Substitution strategies with minor parameter variations are often grouped together in families. For example, the BLOSUM and PAM matrix families are commonly used because they provide a standard, predefined way of modeling substitutions. However, researchers often do not know if a given matrix family or any individual matrix within a family is the most suitable. Furthermore, predefined matrix families may inaccurately reflect a particular hypothesis that a researcher wishes to model or otherwise result in unsatisfactory alignments. In these cases, the ability to compare the effects of one or more custom matrices may be needed. This laborious process is often performed manually because the ability to simultaneously load multiple matrices and then compare their effects on alignments is not readily available in current software tools. This paper presents SubVis, an interactive R package for loading and applying multiple substitution matrices to pairwise alignments. Users can simultaneously explore alignments resulting from multiple predefined and custom substitution matrices. SubVis utilizes several of the alignment functions found in R, a common language among protein scientists. Functions are tied together with the Shiny platform which allows the modification of input parameters. Information regarding alignment quality and individual amino acid substitutions is displayed with the JavaScript language which provides interactive visualizations for revealing both high-level and low-level alignment information.

Measuring protein dynamics in live cells: protocols and practical considerations for fluorescence fluctuation microscopy.

Quantitative analysis of protein complex stoichiometries and mobilities are critical for elucidating the mechanisms that regulate cellular pathways. Fluorescence fluctuation spectroscopy (FFS) techniques can measure protein dynamics, such as diffusion coefficients and formation of complexes, with extraordinary precision and sensitivity. Complete calibration and characterization of the microscope instrument is necessary in order to avoid artifacts during data acquisition and to capitalize on the full capabilities of FFS techniques. We provide an overview of the theory behind FFS techniques, discuss calibration procedures, provide protocols, and give practical considerations for performing FFS experiments. One important parameter recovered from FFS measurements is the relative molecular brightness that can correlate with oligomerization. Three methods for measuring molecular brightness (fluorescence correlation spectroscopy, photon-counting histogram, and number and brightness analysis) recover similar values when measuring samples under ideal conditions in vitro. However, examples are given illustrating that these different methods used for calculating molecular brightness of fluorescent molecules in cells are not always equivalent. Methods relying on spot measurements are more prone to bleaching and movement artifacts that can lead to underestimation of brightness values. We advocate for the use of multiple FFS techniques to study molecular brightnesses to overcome and compliment limitations of individual techniques.

Altered dynamics of a lipid raft associated protein in a kidney model of Fabry disease.

Accumulation of globotriaosylceramide (Gb3) and other neutral glycosphingolipids with galactosyl residues is the hallmark of Fabry disease, a lysosomal storage disorder caused by deficiency of the enzyme alpha-galactosidase A (α-gal A). These lipids are incorporated into the plasma membrane and intracellular membranes, with a preference for lipid rafts. Disruption of raft mediated cell processes is implicated in the pathogenesis of several human diseases, but little is known about the effects of the accumulation of glycosphingolipids on raft dynamics in the context of Fabry disease. Using siRNA technology, we have generated a polarized renal epithelial cell model of Fabry disease in Madin-Darby canine kidney cells. These cells present increased levels of Gb3 and enlarged lysosomes, and progressively accumulate zebra bodies. The polarized delivery of both raft-associated and raft-independent proteins was unaffected by α-gal A knockdown, suggesting that accumulation of Gb3 does not disrupt biosynthetic trafficking pathways. To assess the effect of α-gal A silencing on lipid raft dynamics, we employed number and brightness (N&B) analysis to measure the oligomeric status and mobility of the model glycosylphosphatidylinositol (GPI)-anchored protein GFP-GPI. We observed a significant increase in the oligomeric size of antibody-induced clusters of GFP-GPI at the plasma membrane of α-gal A silenced cells compared with control cells. Our results suggest that the interaction of GFP-GPI with lipid rafts may be altered in the presence of accumulated Gb3. The implications of our results with respect to the pathogenesis of Fabry disease are discussed.

Multiple motifs regulate apical sorting of p75 via a mechanism that involves dimerization and higher-order oligomerization.

The sorting signals that direct proteins to the apical surface of polarized epithelial cells are complex and can include posttranslational modifications, such as N- and O-linked glycosylation. Efficient apical sorting of the neurotrophin receptor p75 is dependent on its O-glycosylated membrane proximal stalk, but how this domain mediates targeting is unknown. Protein oligomerization or clustering has been suggested as a common step in the segregation of all apical proteins. Like many apical proteins, p75 forms dimers, and we hypothesized that formation of higher-order clusters mediated by p75 dimerization and interactions of the stalk facilitate its apical sorting. Using fluorescence fluctuation techniques (photon-counting histogram and number and brightness analyses) to study p75 oligomerization status in vivo, we found that wild-type p75-green fluorescent protein forms clusters in the trans-Golgi network (TGN) but not at the plasma membrane. Disruption of either the dimerization motif or the stalk domain impaired both clustering and polarized delivery. Manipulation of O-glycan processing or depletion of multiple galectins expressed in Madin-Darby canine kidney cells had no effect on p75 sorting, suggesting that the stalk domain functions as a structural prop to position other determinants in the lumenal domain of p75 for oligomerization. Additionally, a p75 mutant with intact dimerization and stalk motifs but with a dominant basolateral sorting determinant (Δ250 mutant) did not form oligomers, consistent with a requirement for clustering in apical sorting. Artificially enhancing dimerization restored clustering to the Δ250 mutant but was insufficient to reroute this mutant to the apical surface. Together these studies demonstrate that clustering in the TGN is required for normal biosynthetic apical sorting of p75 but is not by itself sufficient to reroute a protein to the apical surface in the presence of a strong basolateral sorting determinant. Our studies shed new light on the hierarchy of polarized sorting signals and on the mechanisms by which newly synthesized proteins are segregated in the TGN for eventual apical delivery.

Sialylation of N-linked glycans mediates apical delivery of endolyn in MDCK cells via a galectin-9-dependent mechanism.

The sialomucin endolyn is implicated in adhesion, migration, and differentiation of various cell types. Along rat kidney tubules, endolyn is variously localized to the apical surface and endosomal/lysosomal compartments. Apical delivery of newly synthesized rat endolyn predominates over direct lysosomal delivery in polarized Madin-Darby canine kidney cells. Apical sorting depends on terminal processing of a subset of lumenal N-glycans. Here we dissect the requirements of N-glycan processing for apical targeting and investigate the underlying mechanism. Modulation of glycan branching and subsequent polylactosamine elongation by knockdown of N-acetylglucosaminyltransferase III or V had no effect on apical delivery of endolyn. In contrast, combined but not individual knockdown of sialyltransferases ST3Gal-III, ST3Gal-IV, and ST6Gal-I, which together are responsible for addition of α2,3- and α2,6-linked sialic acids on N-glycans, dramatically decreased endolyn surface polarity. Endolyn synthesized in the presence of kifunensine, which blocks terminal N-glycan processing, reduced its interaction with several recombinant canine galectins, and knockdown of galectin-9 (but not galectin-3, -4, or -8) selectively disrupted endolyn polarity. Our data suggest that sialylation enables recognition of endolyn by galectin-9 to mediate efficient apical sorting. They raise the intriguing possibility that changes in glycosyltransferase expression patterns and/or galectin-9 distribution may acutely modulate endolyn trafficking in the kidney.

FK506 binding protein 8 peptidylprolyl isomerase activity manages a late stage of cystic fibrosis transmembrane conductance regulator (CFTR) folding and stability.

Cystic fibrosis (CF) is caused by mutations in the apical chloride channel cystic fibrosis transmembrane conductance regulator (CFTR) with 90% of patients carrying at least one deletion of the F508 (ΔF508) allele. This mutant form of CFTR is characterized by a folding and trafficking defect that prevents exit from the endoplasmic reticulum. We previously reported that ΔF508 CFTR can be recovered in a complex with Hsp90 and its co-chaperones as an on-pathway folding intermediate, suggesting that Δ508 CF disease arises due to a failure of the proteostasis network (PN), which manages protein folding and degradation in the cell. We have now examined the role of FK506-binding protein 8 (FKBP8), a component of the CFTR interactome, during the biogenesis of wild-type and ΔF508 CFTR. FKBP8 is a member of the peptidylprolyl isomerase family that mediates the cis/trans interconversion of peptidyl prolyl bonds. Our results suggest that FKBP8 is a key PN factor required at a post-Hsp90 step in CFTR biogenesis. In addition, changes in its expression level or alteration of its activity by a peptidylprolyl isomerase inhibitor alter CFTR stability and transport. We propose that CF is caused by the sequential failure of the prevailing PN pathway to stabilize ΔF508-CFTR for endoplasmic reticulum export, a pathway that can be therapeutically managed.

Multiple biosynthetic trafficking routes for apically secreted proteins in MDCK cells.

Many newly synthesized membrane proteins traverse endocytic intermediates en route to the surface in polarized epithelial cells; however, the biosynthetic itinerary of secreted proteins has not been elucidated. We monitored the trafficking route of two secreted proteins with different apical sorting signals: the N-glycan-dependent cargo glycosylated growth hormone (gGH) and Ensol, a soluble version of endolyn whose apical sorting is independent of N-glycans. Both proteins were observed to colocalize in part with apical recycling endosome (ARE) markers. Cargo that lacks an apical targeting signal and is secreted in a nonpolarized manner did not localize to the ARE. Expression of a dominant-negative mutant of myosin Vb, which disrupts ARE export of glycan-dependent membrane proteins, selectively inhibited apical release of gGH but not Ensol. Fluorescence recovery after photobleaching (FRAP) measurements revealed that gGH in the ARE was less mobile than Ensol, consistent with tethering to a sorting receptor. However, knockdown of galectin-3 or galectin-4, lectins implicated in apical sorting, had no effect on the rate or polarity of gGH secretion. Together, our results suggest that apically secreted cargoes selectively access the ARE and are exported via differentially regulated pathways.

At the crossroads of homoeostasis and disease: roles of the PACS proteins in membrane traffic and apoptosis.

The endomembrane system in mammalian cells has evolved over the past two billion years from a simple endocytic pathway in a single-celled primordial ancestor to complex networks supporting multicellular structures that form metazoan tissue and organ systems. The increased organellar complexity of metazoan cells requires additional trafficking machinery absent in yeast or other unicellular organisms to maintain organ homoeostasis and to process the signals that control proliferation, differentiation or the execution of cell death programmes. The PACS (phosphofurin acidic cluster sorting) proteins are one such family of multifunctional membrane traffic regulators that mediate organ homoeostasis and have important roles in diverse pathologies and disease states. This review summarizes our current knowledge of the PACS proteins, including their structure and regulation in cargo binding, their genetics, their roles in secretory and endocytic pathway traffic, interorganellar communication and how cell-death signals reprogramme the PACS proteins to regulate apoptosis. We also summarize our current understanding of how PACS genes are dysregulated in cancer and how viral pathogens ranging from HIV-1 to herpesviruses have evolved to usurp the PACS sorting machinery to promote virus assembly, viral spread and immunoevasion.

Akt and 14-3-3 control a PACS-2 homeostatic switch that integrates membrane traffic with TRAIL-induced apoptosis.

TRAIL selectively kills diseased cells in vivo, spurring interest in this death ligand as a potential therapeutic. However, many cancer cells are resistant to TRAIL, suggesting the mechanism mediating TRAIL-induced apoptosis is complex. Here we identify PACS-2 as an essential TRAIL effector, required for killing tumor cells in vitro and virally infected hepatocytes in vivo. PACS-2 is phosphorylated at Ser437 in vivo, and pharmacologic and genetic studies demonstrate Akt is an in vivo Ser437 kinase. Akt cooperates with 14-3-3 to regulate the homeostatic and apoptotic properties of PACS-2 that mediate TRAIL action. Phosphorylated Ser437 binds 14-3-3 with high affinity, which represses PACS-2 apoptotic activity and is required for PACS-2 to mediate trafficking of membrane cargo. TRAIL triggers dephosphorylation of Ser437, reprogramming PACS-2 to promote apoptosis. Together, these studies identify the phosphorylation state of PACS-2 Ser437 as a molecular switch that integrates cellular homeostasis with TRAIL-induced apoptosis.

HIV-1 Nef binds PACS-2 to assemble a multikinase cascade that triggers major histocompatibility complex class I (MHC-I) down-regulation: analysis using short interfering RNA and knock-out mice.

Human immunodeficiency virus, type 1, negative factor (Nef) initiates down-regulation of cell-surface major histocompatibility complex-I (MHC-I) by assembling an Src family kinase (SFK)-ZAP70/Syk-phosphoinositide 3-kinase (PI3K) cascade through the sequential actions of two sites, Nef EEEE(65) and PXXP(75). The internalized MHC-I molecules are then sequestered in endosomal compartments by a process requiring Nef Met(20). How Nef assembles the multikinase cascade to trigger the MHC-I down-regulation pathway is unknown. Here we report that EEEE(65)-dependent binding to the sorting protein PACS-2 targets Nef to the paranuclear region, enabling PXXP(75) to bind and activate a trans-Golgi network (TGN)-localized SFK. This SFK then phosphorylates ZAP-70 to recruit class I PI3K by interaction with the p85 C-terminal Src homology 2 domain. Using splenocytes and embryonic fibroblasts from PACS-2(-/-) mice, we confirm genetically that Nef requires PACS-2 to localize to the paranuclear region and assemble the multikinase cascade. Moreover, genetic loss of PACS-2 or inhibition of class I PI3K prevents Nef-mediated MHC-I down-regulation, demonstrating that short interfering RNA knockdown of PACS-2 phenocopies the gene knock-out. This PACS-2-dependent targeting pathway is not restricted to Nef, because PACS-2 is also required for trafficking of an endocytosed cation-independent mannose 6-phosphate receptor reporter from early endosomes to the TGN. Together, these results demonstrate PACS-2 is required for Nef action and sorting of itinerant membrane cargo in the TGN/endosomal system.

Cysteine string protein monitors late steps in cystic fibrosis transmembrane conductance regulator biogenesis.

We examined the role of the cysteine string protein (Csp) in cystic fibrosis transmembrane conductance regulator (CFTR) biogenesis in relation to another J-domain protein, Hdj-2, a recognized CFTR cochaperone. Increased expression of Csp produced a dose-dependent reduction in mature (band C) CFTR and an increase in immature (band B) CFTR. Exogenous expression of Hdj-2 also increased CFTR band B, but unlike Csp, Hdj-2 increased band C as well. The Csp-induced block of CFTR maturation required Hsp70, because a J-domain mutant (H43Q) that interferes with the ability of Csp to stimulate Hsp70 ATPase activity relieved the Csp-induced block of CFTR maturation. Nevertheless, Csp H43Q still increased immature CFTR. Csp-induced band B CFTR was found adjacent to the nucleus, co-localizing with calnexin, and it remained detergent-soluble. These data indicate that Csp did not block CFTR maturation by promoting the aggregation or degradation of immature CFTR. Csp knockdown by RNA interference produced a 5-fold increase in mature CFTR and augmented cAMP-stimulated CFTR currents. Thus, the production of mature CFTR is inversely related to the expression level of Csp. Both Csp and Hdj-2 associated with the CFTR R-domain in vitro, and Hdj-2 binding was displaced by Csp, suggesting common interaction sites. Combined expression of Csp and Hdj-2 mimicked the effect of Csp alone, a block of CFTR maturation. But together, Csp and Hdj-2 produced additive increases in CFTR band B, and this did not depend on their interactions with Hsp70, consistent with direct chaperone actions of these proteins. Like Hdj-2, Csp reduced the aggregation of NBD1 in vitro in the absence of Hsp70. Our data suggest that both Csp and Hdj-2 facilitate the biosynthesis of immature CFTR, acting as direct CFTR chaperones, but in addition, Csp is positioned later in the CFTR biogenesis cascade where it regulates the production of mature CFTR by limiting its exit from the endoplasmic reticulum.

Distinct roles for the Hsp40 and Hsp90 molecular chaperones during cystic fibrosis transmembrane conductance regulator degradation in yeast.

Aberrant secreted proteins can be destroyed by ER-associated protein degradation (ERAD), and a prominent, medically relevant ERAD substrate is the cystic fibrosis transmembrane conductance regulator (CFTR). To better define the chaperone requirements during CFTR maturation, the protein was expressed in yeast. Because Hsp70 function impacts CFTR biogenesis in yeast and mammals, we first sought ER-associated Hsp40 cochaperones involved in CFTR maturation. Ydj1p and Hlj1p enhanced Hsp70 ATP hydrolysis but CFTR degradation was slowed only in yeast mutated for both YDJ1 and HLJ1, suggesting functional redundancy. In contrast, CFTR degradation was accelerated in an Hsp90 mutant strain, suggesting that Hsp90 preserves CFTR in a folded state, and consistent with this hypothesis, Hsp90 maintained the solubility of an aggregation-prone domain (NBD1) in CFTR. Soluble ERAD substrate degradation was unaffected in the Hsp90 or the Ydj1p/Hlj1p mutants, and surprisingly CFTR degradation was unaffected in yeast mutated for Hsp90 cochaperones. These results indicate that Hsp90, but not the Hsp90 complex, maintains CFTR structural integrity, whereas Ydj1p/Hlj1p catalyze CFTR degradation.

Localization of the BiP molecular chaperone with respect to endoplasmic reticulum foci containing the cystic fibrosis transmembrane conductance regulator in yeast.

Almost all secreted proteins pass through the endoplasmic reticulum (ER), an organelle that is equipped to tolerate and/or degrade misfolded proteins. We report here that yeast expressing the cystic fibrosis transmembrane conductance regulator (CFTR) concentrate the protein at defined sites in the ER membrane that are not necessarily enriched for the ER molecular chaperone BiP. We propose that these sites are Russell bodies, an ER subcompartment in which misfolded proteins are stored and can be targeted for degradation.