Proteomic Analyses Reveal Divergent Ubiquitylation Site Patterns in Murine Tissues

Posttranslational modifications of proteins increase the complexity of the cellular proteome and enable rapid regulation of protein functions in response to environmental changes. Protein ubiquitylation is a central regulatory posttranslational modification that controls numerous biological processes including proteasomal degradation of proteins, DNA damage repair and innate immune responses. Here we combine high-resolution mass spectrometry with single-step immunoenrichment of di-glycine modified peptides for mapping of endogenous putative ubiquitylation sites in murine tissues. We identify more than 20,000 unique ubiquitylation sites on proteins involved in diverse biological processes. Our data reveals that ubiquitylation regulates core signaling pathways common for each of the studied tissues. In addition, we discover that ubiquitylation regulates tissue-specific signaling networks. Many tissue-specific ubiquitylation sites were obtained from brain highlighting the complexity and unique physiology of this organ. We further demonstrate that different di-glycine-lysine-specific monoclonal antibodies exhibit sequence preferences, and that their complementary use increases the depth of ubiquitylation site analysis, thereby providing a more unbiased view of protein ubiquitylation.Ubiquitin is a small 76-amino-acid protein that is conjugated to the ε-amino group of lysines in a highly orchestrated enzymatic cascade involving ubiquitin activating (E1), ubiquitin conjugating (E2), and ubiquitin ligase (E3) enzymes (1). Ubiquitylation is involved in the regulation of diverse cellular processes including protein degradation (2, 3, 4), DNA damage repair (5, 6), DNA replication (7), cell surface receptor endocytosis, and innate immune signaling (8, 9). Deregulation of protein ubiquitylation is implicated in the development of cancer and neurodegenerative diseases (10, 11). Inhibitors targeting the ubiquitin proteasome system are used in the treatment of hematologic malignancies such as multiple myeloma (12, 13).Recent developments in the mass spectrometry (MS)-based proteomics have greatly expedited proteome-wide analysis of posttranslational modifications (PTMs) (14–17). Large-scale mapping of ubiquitylation sites by mass spectrometry is based on the identification of the di-glycine remnant that results from trypsin digestion of ubiquitylated proteins and remains attached to ubiquitylated lysines (18). Recently, two monoclonal antibodies were developed that specifically recognize di-glycine remnant modified peptides enabling their efficient enrichment from complex peptide mixtures (19, 20). These antibodies have been used to identify thousands of endogenous ubiquitylation sites in human cells, and to quantify site-specific changes in ubiquitylation in response to different cellular perturbations (20–22). It should be noted that the di-glycine remnant is not specific for proteins modified by ubiquitin but also proteins modified by NEDD8 or ISG15 generate an identical di-glycine remnant on modified lysines making it impossible to distinguish between these modifications by mass spectrometry. However, expression of NEDD8 in mouse tissues was shown to be developmentally down-regulated (23), and ISG15 expression in bovine tissues is low in the absence of interferon stimulation (24). In cell culture experiments it was shown that a great majority of sites identified using di-glycine-lysine-specific antibodies stems from ubiquitylated peptides (20).The rates of cell proliferation and protein turnover in mammals vary dramatically between different tissues. Immortalized cell lines, often derived from cancer, are selected for high proliferation rates and fail to represent the complex conditions in tissues. Tissue proteomics can help to gain a more comprehensive understanding of physiological processes in multicellular organisms. Analysis of tissue proteome and PTMs can provide important insights into tissue-specific processes and signaling networks that regulate these processes (25–32). In addition, development of mass spectrometric methods for analysis of PTMs in diseased tissues might lead to the identification of biomarkers.In this study, we combined high-resolution mass spectrometry with immunoenrichment of di-glycine modified peptides to investigate endogenous ubiquitylation sites in murine tissues. We identified more than 20,000 ubiquitylation sites from five different murine tissues and report the largest ubiquitylation dataset obtained from mammalian tissues to date. Furthermore, we compared the performance of the two monoclonal di-glycine-lysine-specific antibodies available for enrichment of ubiquitylated peptides, and reveal their relative preferences for different amino acids flanking ubiquitylation sites.     

Proinsulin and the Genetics of Diabetes Mellitus

Insulin plays a central role in the regulation of vertebrate metabolism. The hormone, the post-translational product of a single-chain precursor, is a globular protein containing two chains, A (21 residues) and B (30 residues). Recent advances in human genetics have identified dominant mutations in the insulin gene causing permanent neonatal-onset DM² (1–4). The mutations are predicted to block folding of the precursor in the ER of pancreatic β-cells. Although expression of the wild-type allele would in other circumstances be sufficient to maintain homeostasis, studies of a corresponding mouse model (5–7) suggest that the misfolded variant perturbs wild-type biosynthesis (8, 9). Impaired β-cell secretion is associated with ER stress, distorted organelle architecture, and cell death (10). These findings have renewed interest in insulin biosynthesis (11–13) and the structural basis of disulfide pairing (14–19). Protein evolution is constrained not only by structure and function but also by susceptibility to toxic misfolding.Insulin plays a central role in the regulation of vertebrate metabolism. The hormone, the post-translational product of a single-chain precursor, is a globular protein containing two chains, A (21 residues) and B (30 residues). Recent advances in human genetics have identified dominant mutations in the insulin gene causing permanent neonatal-onset DM² (1–4). The mutations are predicted to block folding of the precursor in the ER of pancreatic β-cells. Although expression of the wild-type allele would in other circumstances be sufficient to maintain homeostasis, studies of a corresponding mouse model (5–7) suggest that the misfolded variant perturbs wild-type biosynthesis (8, 9). Impaired β-cell secretion is associated with ER stress, distorted organelle architecture, and cell death (10). These findings have renewed interest in insulin biosynthesis (11–13) and the structural basis of disulfide pairing (14–19). Protein evolution is constrained not only by structure and function but also by susceptibility to toxic misfolding.     

Printor, a Novel TorsinA-interacting Protein Implicated in Dystonia Pathogenesis

Early onset generalized dystonia (DYT1) is an autosomal dominant neurological disorder caused by deletion of a single glutamate residue (torsinA ΔE) in the C-terminal region of the AAA⁺ (ATPases associated with a variety of cellular activities) protein torsinA. The pathogenic mechanism by which torsinA ΔE mutation leads to dystonia remains unknown. Here we report the identification and characterization of a 628-amino acid novel protein, printor, that interacts with torsinA. Printor co-distributes with torsinA in multiple brain regions and co-localizes with torsinA in the endoplasmic reticulum. Interestingly, printor selectively binds to the ATP-free form but not to the ATP-bound form of torsinA, supporting a role for printor as a cofactor rather than a substrate of torsinA. The interaction of printor with torsinA is completely abolished by the dystonia-associated torsinA ΔE mutation. Our findings suggest that printor is a new component of the DYT1 pathogenic pathway and provide a potential molecular target for therapeutic intervention in dystonia.Early onset generalized torsion dystonia (DYT1) is the most common and severe form of hereditary dystonia, a movement disorder characterized by involuntary movements and sustained muscle spasms (1). This autosomal dominant disease has childhood onset and its dystonic symptoms are thought to result from neuronal dysfunction rather than neurodegeneration (2, 3). Most DYT1 cases are caused by deletion of a single glutamate residue at positions 302 or 303 (torsinA ΔE) of the 332-amino acid protein torsinA (4). In addition, a different torsinA mutation that deletes amino acids Phe³²³–Tyr³²⁸ (torsinA Δ323–328) was identified in a single family with dystonia (5), although the pathogenic significance of this torsinA mutation is unclear because these patients contain a concomitant mutation in another dystonia-related protein, ϵ-sarcoglycan (6). Recently, genetic association studies have implicated polymorphisms in the torsinA gene as a genetic risk factor in the development of adult-onset idiopathic dystonia (7, 8).TorsinA contains an N-terminal endoplasmic reticulum (ER)³ signal sequence and a 20-amino acid hydrophobic region followed by a conserved AAA⁺ (ATPases associated with a variety of cellular activities) domain (9, 10). Because members of the AAA⁺ family are known to facilitate conformational changes in target proteins (11, 12), it has been proposed that torsinA may function as a molecular chaperone (13, 14). TorsinA is widely expressed in brain and multiple other tissues (15) and is primarily associated with the ER and nuclear envelope (NE) compartments in cells (16–20). TorsinA is believed to mainly reside in the lumen of the ER and NE (17–19) and has been shown to bind lamina-associated polypeptide 1 (LAP1) (21), lumenal domain-like LAP1 (LULL1) (21), and nesprins (22). In addition, recent evidence indicates that a significant pool of torsinA exhibits a topology in which the AAA⁺ domain faces the cytoplasm (20). In support of this topology, torsinA is found in the cytoplasm, neuronal processes, and synaptic terminals (2, 3, 15, 23–26) and has been shown to bind cytosolic proteins snapin (27) and kinesin light chain 1 (20). TorsinA has been proposed to play a role in several cellular processes, including dopaminergic neurotransmission (28–31), NE organization and dynamics (17, 22, 32), and protein trafficking (27, 33). However, the precise biological function of torsinA and its regulation remain unknown.To gain insights into torsinA function, we performed yeast two-hybrid screens to search for torsinA-interacting proteins in the brain. We report here the isolation and characterization of a novel protein named printor (protein interactor of torsinA) that interacts selectively with wild-type (WT) torsinA but not the dystonia-associated torsinA ΔE mutant. Our data suggest that printor may serve as a cofactor of torsinA and provide a new molecular target for understanding and treating dystonia.     

Analysis of Free Energy Signals Arising from Nucleotide Hybridization Between rRNA and mRNA Sequences during Translation in Eubacteria

A decoding algorithm is tested that mechanistically models the progressive alignments that arise as the mRNA moves past the rRNA tail during translation elongation. Each of these alignments provides an opportunity for hybridization between the single-stranded, -terminal nucleotides of the 16S rRNA and the spatially accessible window of mRNA sequence, from which a free energy value can be calculated. Using this algorithm we show that a periodic, energetic pattern of frequency 1/3 is revealed. This periodic signal exists in the majority of coding regions of eubacterial genes, but not in the non-coding regions encoding the 16S and 23S rRNAs. Signal analysis reveals that the population of coding regions of each bacterial species has a mean phase that is correlated in a statistically significant way with species () content. These results suggest that the periodic signal could function as a synchronization signal for the maintenance of reading frame and that codon usage provides a mechanism for manipulation of signal phase.[1,2,3,4,5,6,7,8,9,10,11,12,13,14,15,16,17,18,19,20,21,22,23,24,25,26,27,28,29,30,31,32]     

Semaphorin3A Signaling Mediated by Fyn-dependent Tyrosine Phosphorylation of Collapsin Response Mediator Protein 2 at Tyrosine 32

Collapsin response mediator protein 2 (CRMP2) is an intracellular protein that mediates signaling of Semaphorin3A (Sema3A), a repulsive axon guidance molecule. Fyn, a Src-type tyrosine kinase, is involved in the Sema3A signaling. However, the relationship between CRMP2 and Fyn in this signaling pathway is still unknown. In our research, we demonstrated that Fyn phosphorylated CRMP2 at Tyr³² residues in HEK293T cells. Immunohistochemical analysis using a phospho-specific antibody at Tyr³² of CRMP showed that Tyr³²-phosphorylated CRMP was abundant in the nervous system, including dorsal root ganglion neurons, the molecular and Purkinje cell layer of adult cerebellum, and hippocampal fimbria. Overexpression of a nonphosphorylated mutant (Tyr³² to Phe³²) of CRMP2 in dorsal root ganglion neurons interfered with Sema3A-induced growth cone collapse response. These results suggest that Fyn-dependent phosphorylation of CRMP2 at Tyr³² is involved in Sema3A signaling.Collapsin response mediator proteins (CRMPs)⁴ have been identified as intracellular proteins that mediate Semaphorin3A (Sema3A) signaling in the nervous system (1). CRMP2 is one of the five members of the CRMP family. CRMPs also mediate signal transduction of NT3, Ephrin, and Reelin (2–4). CRMPs interact with several intracellular molecules, including tubulin, Numb, kinesin1, and Sra1 (5–8). CRMPs are involved in axon guidance, axonal elongation, cell migration, synapse maturation, and the generation of neuronal polarity (1, 2, 4, 5).CRMP family proteins are known to be the major phosphoproteins in the developing brain (1, 9). CRMP2 is phosphorylated by several Ser/Thr kinases, such as Rho kinase, cyclin-dependent kinase 5 (Cdk5), and glycogen synthase kinase 3β (GSK3β) (2, 10–13). The phosphorylation sites of CRMP2 by these kinases are clustered in the C terminus and have already been identified. Rho kinase phosphorylates CRMP2 at Thr⁵⁵⁵ (10). Cdk5 phosphorylates CRMP2 at Ser⁵²², and this phosphorylation is essential for sequential phosphorylations by GSK3β at Ser⁵¹⁸, Thr⁵¹⁴, and Thr⁵⁰⁹ (2, 11–13). These phosphorylations disrupt the interaction of CRMP2 with tubulin or Numb (2, 3, 13). The sequential phosphorylation of CRMP2 by Cdk5 and GSK3β is an essential step in Sema3A signaling (11, 13). Furthermore, the neurofibrillary tangles in the brains of people with Alzheimer disease contain hyperphosphorylated CRMP2 at Thr⁵⁰⁹, Ser⁵¹⁸, and Ser⁵²² (14, 15).CRMPs are also substrates of several tyrosine kinases. The phosphorylation of CRMP2 by Fes/Fps and Fer has been shown to be involved in Sema3A signaling (16, 17). Phosphorylation of CRMP2 at Tyr⁴⁷⁹ by a Src family tyrosine kinase Yes regulates CXCL12-induced T lymphocyte migration (18). We reported previously that Fyn is involved in Sema3A signaling (19). Fyn associates with PlexinA2, one of the components of the Sema3A receptor complex. Fyn also activates Cdk5 through the phosphorylation at Tyr¹⁵ of Cdk5 (19). In dorsal root ganglion (DRG) neurons from fyn-deficient mice, Sema3A-induced growth cone collapse response is attenuated compared with control mice (19). Furthermore, we recently found that Fyn phosphorylates CRMP1 and that this phosphorylation is involved in Reelin signaling (4). Although it has been shown that CRMP2 is involved in Sema3A signaling (1, 11, 13), the relationship between Fyn and CRMP2 in Sema3A signaling and the tyrosine phosphorylation site(s) of CRMPs remain unknown.Here, we show that Fyn phosphorylates CRMP2 at Tyr³². Using a phospho-specific antibody against Tyr³², we determined that the residue is phosphorylated in vivo. A nonphosphorylated mutant CRMP2Y32F inhibits Sema3A-induced growth cone collapse. These results indicate that tyrosine phosphorylation by Fyn at Tyr³² is involved in Sema3A signaling.     

Deep,Quantitative Coverage of the Lysine Acetylome Using Novel Anti-acetyl-lysine Antibodies and an Optimized Proteomic Workflow

Introduction of antibodies specific for acetylated lysine has significantly improved the detection of endogenous acetylation sites by mass spectrometry. Here, we describe a new, commercially available mixture of anti-lysine acetylation (Kac) antibodies and show its utility for in-depth profiling of the acetylome. Specifically, seven complementary monoclones with high specificity for Kac were combined into a final anti-Kac reagent which results in at least a twofold increase in identification of Kac peptides over a commonly used Kac antibody. We outline optimal antibody usage conditions, effective offline basic reversed phase separation, and use of state-of-the-art LC-MS technology for achieving unprecedented coverage of the acetylome. The methods were applied to quantify acetylation sites in suberoylanilide hydroxamic acid-treated Jurkat cells. Over 10,000 Kac peptides from over 3000 Kac proteins were quantified from a single stable isotope labeling by amino acids in cell culture labeled sample using 7.5 mg of peptide input per state. This constitutes the deepest coverage of acetylation sites in quantitative experiments obtained to-date. The approach was also applied to breast tumor xenograft samples using isobaric mass tag labeling of peptides (iTRAQ4, TMT6 and TMT10-plex reagents) for quantification. Greater than 6700 Kac peptides from over 2300 Kac proteins were quantified using 1 mg of tumor protein per iTRAQ 4-plex channel. The novel reagents and methods we describe here enable quantitative, global acetylome analyses with depth and sensitivity approaching that obtained for other well-studied post-translational modifications such as phosphorylation and ubiquitylation, and should have widespread application in biological and clinical studies employing mass spectrometry-based proteomics.Lysine acetylation (Kac)¹ is a well conserved, reversible post-translational modification (PTM) involved in multiple cellular processes (1). Acetylation is regulated by two classes of enzymes: lysine acetyltransferases (KATs) and histone deacetylases (HDACs) (2–4). This modification was originally identified as a nuclear event on histone proteins and has been long appreciated for its role in epigenetic and DNA-dependent processes. With the help of a growing number of large-scale acetylation studies, it has become evident that lysine acetylation is ubiquitous, also occurring on cytoplasmic and mitochondrial proteins and has a role in signaling, metabolism, and immunity (1, 4–6). Therefore, the examination of lysine acetylation on nonhistone proteins has gained a prominent role in PTM analysis.To date, the identification of large numbers of acetylation sites has been challenging because of the substoichiometric nature of this modification (7, 8). Additionally, global acetylation is generally less abundant than phosphorylation and ubiquitylation (1). The introduction of antibodies specific for lysine acetylation has significantly improved the ability to enrich and identify thousands of sites (9–14). A landmark study by Choudhary et al. used anti-Kac antibodies to globally map 3600 lysine acetylation sites on 1750 proteins, thereby demonstrating the feasibility of profiling the acetylome (10). A more recent study by Lundby et al. investigated the function and distribution of acetylation sites in 16 different rat tissues, and identified, in aggregate, 15,474 acetylation sites from 4541 proteins (12).Although anti-acetyl lysine antibodies have been a breakthrough for globally mapping acetylation sites (9–12), it remains a challenge to identify large numbers of lysine acetylation sites from a single sample, as is now routinely possible for phosphorylation and ubiquitylation (13, 15–18). To improve the depth-of-coverage in acetylation profiling experiments there is a clear need for (1) alternative anti-acetyl lysine antibodies with higher specificity, (2) optimized antibody usage parameters, and (3) robust proteomic workflows that permit low to moderate protein input. In this study, we describe a newly commercialized mixture of anti-Kac antibodies and detail a complete proteomic workflow for achieving unprecedented coverage of the acetylome from a single stable isotope labeling by amino acids in cell culture (SILAC) labeled sample as well as isobaric tags for relative and absolute quantitation (iTRAQ)- and tandem mass tag (TMT)-labeled samples.     

Algorithms for Finding Small Attractors in Boolean Networks

A Boolean network is a model used to study the interactions between different genes in genetic regulatory networks. In this paper, we present several algorithms using gene ordering and feedback vertex sets to identify singleton attractors and small attractors in Boolean networks. We analyze the average case time complexities of some of the proposed algorithms. For instance, it is shown that the outdegree-based ordering algorithm for finding singleton attractors works in time for , which is much faster than the naive time algorithm, where is the number of genes and is the maximum indegree. We performed extensive computational experiments on these algorithms, which resulted in good agreement with theoretical results. In contrast, we give a simple and complete proof for showing that finding an attractor with the shortest period is NP-hard.[1,2,3,4,5,6,7,8,9,10,11,12,13,14,15,16,17,18,19,20,21,22,23,24,25,26,27,28,29,30,31,32]     

Conserved Glutamate Residues Glu-343 and Glu-519 Provide Mechanistic Insights into Cation/Nucleoside Cotransport by Human Concentrative Nucleoside Transporter hCNT3

Human concentrative nucleoside transporter 3 (hCNT3) utilizes electrochemical gradients of both Na⁺ and H⁺ to accumulate pyrimidine and purine nucleosides within cells. We have employed radioisotope flux and electrophysiological techniques in combination with site-directed mutagenesis and heterologous expression in Xenopus oocytes to identify two conserved pore-lining glutamate residues (Glu-343 and Glu-519) with essential roles in hCNT3 Na⁺/nucleoside and H⁺/nucleoside cotransport. Mutation of Glu-343 and Glu-519 to aspartate, glutamine, and cysteine severely compromised hCNT3 transport function, and changes included altered nucleoside and cation activation kinetics (all mutants), loss or impairment of H⁺ dependence (all mutants), shift in Na⁺:nucleoside stoichiometry from 2:1 to 1:1 (E519C), complete loss of catalytic activity (E519Q) and, similar to the corresponding mutant in Na⁺-specific hCNT1, uncoupled Na⁺ currents (E343Q). Consistent with close-proximity integration of cation/solute-binding sites within a common cation/permeant translocation pore, mutation of Glu-343 and Glu-519 also altered hCNT3 nucleoside transport selectivity. Both residues were accessible to the external medium and inhibited by p-chloromercuribenzene sulfonate when converted to cysteine.Physiologic nucleosides and the majority of synthetic nucleoside analogs with antineoplastic and/or antiviral activity are hydrophilic molecules that require specialized plasma membrane nucleoside transporter (NT)³ proteins for transport into or out of cells (1–4). NT-mediated transport is required for nucleoside metabolism by salvage pathways and is a critical determinant of the pharmacologic actions of nucleoside drugs (3–6). By regulating adenosine availability to purinoreceptors, NTs also modulate a diverse array of physiological processes, including neurotransmission, immune responses, platelet aggregation, renal function, and coronary vasodilation (4, 6, 7). Two structurally unrelated NT families of integral membrane proteins exist in human and other mammalian cells and tissues as follows: the SLC28 concentrative nucleoside transporter (CNT) family and the SLC29 equilibrative nucleoside transporter (ENT) family (3, 4, 6, 8, 9). ENTs are normally present in most, possibly all, cell types (4, 6, 8). CNTs, in contrast, are found predominantly in intestinal and renal epithelia and other specialized cell types, where they have important roles in absorption, secretion, distribution, and elimination of nucleosides and nucleoside drugs (1–3, 5, 6, 9).The CNT protein family in humans is represented by three members, hCNT1, hCNT2, and hCNT3. Belonging to a CNT subfamily phylogenetically distinct from hCNT1/2, hCNT3 utilizes electrochemical gradients of both Na⁺ and H⁺ to accumulate a broad range of pyrimidine and purine nucleosides and nucleoside drugs within cells (10, 11). hCNT1 and hCNT2, in contrast, are Na⁺-specific and transport pyrimidine and purine nucleosides, respectively (11–13). Together, hCNT1–3 account for the three major concentrative nucleoside transport processes of human and other mammalian cells. Nonmammalian members of the CNT protein family that have been characterized functionally include hfCNT, a second member of the CNT3 subfamily from the ancient marine prevertebrate the Pacific hagfish Eptatretus stouti (14), CeCNT3 from Caenorhabditis elegans (15), CaCNT from Candida albicans (16), and the bacterial nucleoside transporter NupC from Escherichia coli (17). hfCNT is Na⁺- but not H⁺-coupled, whereas CeCNT3, CaCNT, and NupC are exclusively H⁺-coupled. Na⁺:nucleoside coupling stoichiometries are 1:1 for hCNT1 and hCNT2 and 2:1 for hCNT3 and hfCNT3 (11, 14). H⁺:nucleoside coupling ratios for hCNT3 and CaCNT are 1:1 (11, 16).Although much progress has been made in molecular studies of ENT proteins (4, 6, 8), studies of structurally and functionally important regions and residues within the CNT protein family are still at an early stage. Topological investigations suggest that hCNT1–3 and other eukaryote CNT family members have a 13 (or possibly 15)-transmembrane helix (TM) architecture, and multiple alignments reveal strong sequence similarities within the C-terminal half of the proteins (18). Prokaryotic CNTs lack the first three TMs of their eukaryotic counterparts, and functional expression of N-terminally truncated human and rat CNT1 in Xenopus oocytes has established that these three TMs are not required for Na⁺-dependent uridine transport activity (18). Consistent with this finding, chimeric studies involving hCNT1 and hfCNT (14) and hCNT1 and hCNT3 (19) have demonstrated that residues involved in Na⁺- and H⁺-coupling reside in the C-terminal half of the protein. Present in this region of the transporter, but of unknown function, is a highly conserved (G/A)XKX₃NEFVA(Y/M/F) motif common to all eukaryote and prokaryote CNTs.By virtue of their negative charge and consequent ability to interact directly with coupling cations and/or participate in cation-induced and other protein conformational transitions, glutamate and aspartate residues play key functional and structural roles in a broad spectrum of mammalian and bacterial cation-coupled transporters (20–30). Little, however, is known about their role in CNTs. This study builds upon a recent mutagenesis study of conserved glutamate and aspartate residues in hCNT1 (31) to undertake a parallel in depth investigation of corresponding residues in hCNT3. By employing the multifunctional capability of hCNT3 as a template for these studies, this study provides novel mechanistic insights into the molecular mechanism(s) of CNT-mediated cation/nucleoside cotransport, including the role of the (G/A)XKX₃NEFVA(Y/M/F) motif.     

The Wavelet-Based Cluster Analysis for Temporal Gene Expression Data

A variety of high-throughput methods have made it possible to generate detailed temporal expression data for a single gene or large numbers of genes. Common methods for analysis of these large data sets can be problematic. One challenge is the comparison of temporal expression data obtained from different growth conditions where the patterns of expression may be shifted in time. We propose the use of wavelet analysis to transform the data obtained under different growth conditions to permit comparison of expression patterns from experiments that have time shifts or delays. We demonstrate this approach using detailed temporal data for a single bacterial gene obtained under 72 different growth conditions. This general strategy can be applied in the analysis of data sets of thousands of genes under different conditions.[1,2,3,4,5,6,7,8,9,10,11,12,13,14,15,16,17,18,19,20,21,22,23,24,25,26,27,28,29]     

Protein-tyrosine Phosphatase-�� and Src Functionally Link Focal Adhesions to the Endoplasmic Reticulum to Mediate Interleukin-1-induced Ca2+ Signaling

Calcium (Ca²⁺) signaling by the pro-inflammatory cytokine interleukin-1 (IL-1) is dependent on focal adhesions, which contain diverse structural and signaling proteins including protein phosphatases. We examined here the role of protein-tyrosine phosphatase (PTP) α in regulating IL-1-induced Ca²⁺ signaling in fibroblasts. IL-1 promoted recruitment of PTPα to focal adhesions and endoplasmic reticulum (ER) fractions, as well as tyrosine phosphorylation of the ER Ca²⁺ release channel IP₃R. In response to IL-1, catalytically active PTPα was required for Ca²⁺ release from the ER, Src-dependent phosphorylation of IP₃R1 and accumulation of IP₃R1 in focal adhesions. In pulldown assays and immunoprecipitations PTPα was required for the association of PTPα with IP₃R1 and c-Src, and this association was increased by IL-1. Collectively, these data indicate that PTPα acts as an adaptor to mediate functional links between focal adhesions and the ER that enable IL-1-induced Ca²⁺ signaling.The interleukin-1 (IL-1)³ family of pro-inflammatory cytokines mediates host responses to infection and injury. Impaired control of IL-1 signaling leads to chronic inflammation and destruction of extracellular matrices (1, 2), as seen in pathological conditions such as pulmonary fibrosis (3), rheumatoid arthritis (4, 5), and periodontitis (6). IL-1 elicits multiple signaling programs, some of which trigger Ca²⁺ release from the endoplasmic reticulum (ER) as well as expression of multiple cytokines and inflammatory factors including c-Fos and c-Jun (7, 8), and matrix metalloproteinases (9, 10), which mediate extracellular matrix degradation via mitogen-activated protein kinase-regulated pathways (11).In anchorage-dependent cells including fibroblasts and chondrocytes, focal adhesions (FAs) are required for IL-1-induced Ca²⁺ release from the ER and activation of ERK (12–14). FAs are actin-enriched adhesive domains composed of numerous (>50) scaffolding and signaling proteins (15–17). Many FA proteins are tyrosine-phosphorylated, including paxillin, focal adhesion kinase, and src family kinases, all of which are crucial for the assembly and disassembly of FAs (18–21). Protein-tyrosine phosphorylation plays a central role in regulating many cellular processes including adhesion (22, 23), motility (24), survival (25), and signal transduction (26–29). Phosphorylation of proteins by kinases is balanced by protein-tyrosine phosphatases (PTP), which can enhance or attenuate downstream signaling by dephosphorylation of tyrosine residues (30–32).PTPs can be divided into two main categories: receptor-like and intracellular PTPs (33). Two receptor-like PTPs have been localized to FA (leukocyte common antigen-related molecule and PTPα). Leukocyte common antigen-related molecule can dephosphorylate and mediate degradation of p130^cas, which ultimately leads to cell death (34, 35). PTPα contains a heavily glycosylated extracellular domain, a transmembrane domain, and two intracellular phosphatase domains (33, 36). The amino-terminal domain predominantly mediates catalytic activity, whereas the carboxyl-terminal domain serves a regulatory function (37, 38). PTPα is enriched in FA (23) and is instrumental in regulating FA dynamics (39) via activation of c-Src/Fyn kinases by dephosphorylating the inhibitory carboxyl tyrosine residue, namely Tyr⁵²⁹ (22, 40–42) and facilitation of integrin-dependent assembly of Src-FAK and Fyn-FAK complexes that regulate cell motility (43). Although PTPα has been implicated in formation and remodeling of FAs (44, 45), the role of PTPα in FA-dependent signaling is not defined.Ca²⁺ release from the ER is a critical step in integrin-dependent IL-1 signal transduction and is required for downstream activation of ERK (13, 46). The release of Ca²⁺ from the ER depends on the inositol 1,4,5-triphosphate receptor (IP₃R), which is an IP₃-gated Ca²⁺ channel (47). All of the IP₃R subtypes (subtypes 1–3) have been localized to the ER, as well as other the plasma membrane and other endomembranes (48–50). Further, IP₃R may associate with FAs, enabling the anchorage of the ER to FAs (51, 52). However, the molecule(s) that provide the structural link for this association has not been defined.FA-restricted, IL-1-triggered signal transduction in anchorage-dependent cells may rely on interacting proteins that are enriched in FAs and the ER (53). Here, we examined the possibility that PTPα associates with c-Src and IP₃R to functionally link FAs to the ER, thereby enabling IL-1 signal transduction.     

Deletion of the Chloride Transporter Slc26a7 Causes Distal Renal Tubular Acidosis and Impairs Gastric Acid Secretion

SLC26A7 (human)/Slc26a7 (mouse) is a recently identified chloride-base exchanger and/or chloride transporter that is expressed on the basolateral membrane of acid-secreting cells in the renal outer medullary collecting duct (OMCD) and in gastric parietal cells. Here, we show that mice with genetic deletion of Slc26a7 expression develop distal renal tubular acidosis, as manifested by metabolic acidosis and alkaline urine pH. In the kidney, basolateral Cl⁻/HCO₃⁻ exchange activity in acid-secreting intercalated cells in the OMCD was significantly decreased in hypertonic medium (a normal milieu for the medulla) but was reduced only mildly in isotonic medium. Changing from a hypertonic to isotonic medium (relative hypotonicity) decreased the membrane abundance of Slc26a7 in kidney cells in vivo and in vitro. In the stomach, stimulated acid secretion was significantly impaired in isolated gastric mucosa and in the intact organ. We propose that SLC26A7 dysfunction should be investigated as a potential cause of unexplained distal renal tubular acidosis or decreased gastric acid secretion in humans.The collecting duct segment of the distal kidney nephron plays a major role in systemic acid base homeostasis by acid secretion and bicarbonate absorption. The acid secretion occurs via H⁺-ATPase and H-K-ATPase into the lumen and bicarbonate is absorbed via basolateral Cl⁻/HCO₃⁻ exchangers (1–4). The tubules, which are located within the outer medullary region of the kidney collecting duct (OMCD),² have the highest rate of acid secretion among the distal tubule segments and are therefore essential to the maintenance of acid base balance (2).The gastric parietal cell is the site of generation of acid and bicarbonate through the action of cytosolic carbonic anhydrase II (5, 6). The intracellular acid is secreted into the lumen via gastric H-K-ATPase, which works in conjunction with a chloride channel and a K⁺ recycling pathway (7–10). The intracellular bicarbonate is transported to the blood via basolateral Cl⁻/HCO₃⁻ exchangers (11–14).SLC26 (human)/Slc26 (mouse) isoforms are members of a conserved family of anion transporters that display tissue-specific patterns of expression in epithelial cells (15–24). Several SLC26 members can function as chloride/bicarbonate exchangers. These include SLC26A3 (DRA), SLC26A4 (pendrin), SLC26A6 (PAT1 or CFEX), SLC26A7, and SLC26A9 (25–31). SLC26A7 and SLC26A9 can also function as chloride channels (32–34).SLC26A7/Slc26a7 is predominantly expressed in the kidney and stomach (28, 29). In the kidney, Slc26a7 co-localizes with AE1, a well-known Cl⁻/HCO₃⁻ exchanger, on the basolateral membrane of (acid-secreting) A-intercalated cells in OMCD cells (29, 35, 36) (supplemental Fig. 1). In the stomach, Slc26a7 co-localizes with AE2, a major Cl⁻/HCO₃⁻ exchanger, on the basolateral membrane of acid secreting parietal cells (28). To address the physiological function of Slc26a7 in the intact mouse, we have generated Slc26a7 ko mice. We report here that Slc26a7 ko mice exhibit distal renal tubular acidosis and impaired gastric acidification in the absence of morphological abnormalities in kidney or stomach.     

The Cell Tropism of Human Immunodeficiency Virus Type 1 Determines the Kinetics of Plasma Viremia in SCID Mice Reconstituted with Human Peripheral Blood Leukocytes

Most individuals infected with human immunodeficiency virus type 1 (HIV-1) initially harbor macrophage-tropic, non-syncytium-inducing (M-tropic, NSI) viruses that may evolve into T-cell-tropic, syncytium-inducing viruses (T-tropic, SI) after several years. The reasons for the more efficient transmission of M-tropic, NSI viruses and the slow evolution of T-tropic, SI viruses remain unclear, although they may be linked to expression of appropriate chemokine coreceptors for virus entry. We have examined plasma viral RNA levels and the extent of CD4⁺ T-cell depletion in SCID mice reconstituted with human peripheral blood leukocytes following infection with M-tropic, dual-tropic, or T-tropic HIV-1 isolates. The cell tropism was found to determine the course of viremia, with M-tropic viruses producing sustained high viral RNA levels and sparing some CD4⁺ T cells, dual-tropic viruses producing a transient and lower viral RNA spike and extremely rapid depletion of CD4⁺ T cells, and T-tropic viruses causing similarly lower viral RNA levels and rapid-intermediate rates of CD4⁺ T-cell depletion. A single amino acid change in the V3 region of gp120 was sufficient to cause one isolate to switch from M-tropic to dual-tropic and acquire the ability to rapidly deplete all CD4⁺ T cells.The envelope gene of human immunodeficiency virus type 1 (HIV-1) determines the cell tropism of the virus (11, 32, 47, 62), the use of chemokine receptors as cofactors for viral entry (4, 17), and the ability of the virus to induce syncytia in infected cells (55, 60). Cell tropism is closely linked to but probably not exclusively determined by the ability of different HIV-1 envelopes to bind CD4 and the CC or the CXC chemokine receptors and initiate viral fusion with the target cell. Macrophage-tropic (M-tropic) viruses infect primary cultures of macrophages and CD4⁺ T cells and use CCR5 as the preferred coreceptor (2, 5, 15, 23, 26, 31). T-cell-tropic (T-tropic) viruses can infect primary cultures of CD4⁺ T cells and established T-cell lines, but not primary macrophages. T-tropic viruses use CXCR4 as a coreceptor for viral entry (27). Dual-tropic viruses have both of these properties and can use either CCR5 or CXCR4 (and infrequently other chemokine receptors [25]) for viral entry (24, 37, 57). M-tropic viruses are most frequently transmitted during primary infection of humans and persist throughout the duration of the infection (63). Many, but not all, infected individuals show an evolution of virus cell tropism from M-tropic to dual-tropic and finally to T-tropic with increasing time after infection (21, 38, 57). Increases in replicative capacity of viruses from patients with long-term infection have also been noted (22), and the switch to the syncytium-inducing (SI) phenotype in T-tropic or dual-tropic isolates is associated with more rapid disease progression (10, 20, 60). Primary infection with dual-tropic or T-tropic HIV, although infrequent, often leads to rapid disease progression (16, 51). The viral and host factors that determine the higher transmission rate of M-tropic HIV-1 and the slow evolution of dual- or T-tropic variants remain to be elucidated (4).These observations suggest that infection with T-tropic, SI virus isolates in animal model systems with SCID mice grafted with human lymphoid cells or tissue should lead to a rapid course of disease (1, 8, 44–46). While some studies in SCID mice grafted with fetal thymus and liver are in agreement with this concept (33, 34), our previous studies with the human peripheral blood leukocyte-SCID (hu-PBL-SCID) mouse model have shown that infection with M-tropic isolates (e.g., SF162) causes more rapid CD4⁺ T-cell depletion than infection with T-tropic, SI isolates (e.g., SF33), despite similar proviral copy numbers, and that this property mapped to envelope (28, 41, 43). However, the dual-tropic 89.6 isolate (19) caused extremely rapid CD4⁺ T-cell depletion in infected hu-PBL-SCID mice that was associated with an early and transient increase in HIV-1 plasma viral RNA (29). The relationship between cell tropism of the virus isolate and the pattern of disease in hu-PBL-SCID mice is thus uncertain. We have extended these studies by determining the kinetics of HIV-1 RNA levels in serial plasma samples of hu-PBL-SCID mice infected with primary patient isolates or laboratory stocks that differ in cell tropism and SI properties. The results showed significant differences in the kinetics of HIV-1 replication and CD4⁺ T-cell depletion that are determined by the cell tropism of the virus isolate.     

A Ca2+ Release-activated Ca2+ (CRAC) Modulatory Domain (CMD) within STIM1 Mediates Fast Ca2+-dependent Inactivation of ORAI1 Channels

STIM1 and ORAI1, the two limiting components in the Ca²⁺ release-activated Ca²⁺ (CRAC) signaling cascade, have been reported to interact upon store depletion, culminating in CRAC current activation. We have recently identified a modulatory domain between amino acids 474 and 485 in the cytosolic part of STIM1 that comprises 7 negatively charged residues. A STIM1 C-terminal fragment lacking this domain exhibits enhanced interaction with ORAI1 and 2–3-fold higher ORAI1/CRAC current densities. Here we focused on the role of this CRAC modulatory domain (CMD) in the fast inactivation of ORAI1/CRAC channels, utilizing the whole-cell patch clamp technique. STIM1 mutants either with C-terminal deletions including CMD or with 7 alanines replacing the negative amino acids within CMD gave rise to ORAI1 currents that displayed significantly reduced or even abolished inactivation when compared with STIM1 mutants with preserved CMD. Consistent results were obtained with cytosolic C-terminal fragments of STIM1, both in ORAI1-expressing HEK 293 cells and in RBL-2H3 mast cells containing endogenous CRAC channels. Inactivation of the latter, however, was much more pronounced than that of ORAI1. The extent of inactivation of ORAI3 channels, which is also considerably more prominent than that of ORAI1, was also substantially reduced by co-expression of STIM1 constructs missing CMD. Regarding the dependence of inactivation on Ca²⁺, a decrease in intracellular Ca²⁺ chelator concentrations promoted ORAI1 current fast inactivation, whereas Ba²⁺ substitution for extracellular Ca²⁺ completely abrogated it. In summary, CMD within the STIM1 cytosolic part provides a negative feedback signal to Ca²⁺ entry by triggering fast Ca²⁺-dependent inactivation of ORAI/CRAC channels.The Ca²⁺ release-activated Ca²⁺ (CRAC)⁵ channel is one of the best characterized store-operated entry pathways (1–7). Substantial efforts have led to identification of two key components of the CRAC channel machinery: the stromal interaction molecule 1 (STIM1), which is located in the endoplasmic reticulum and acts as a Ca²⁺ sensor (8–10), and ORAI1/CRACM1, the pore-forming subunit of the CRAC channel (11–13). Besides ORAI1, two further homologues named ORAI2 and ORAI3 belong to the ORAI channel family (12, 14).STIM1 senses endoplasmic reticulum store depletion primarily by its luminal EF-hand in its N terminus (8, 15), redistributes close to the plasma membrane, where it forms puncta-like structures, and co-clusters with ORAI1, leading to inward Ca²⁺ currents (12, 16–19). The STIM1 C terminus, located in the cytosol, contains two coiled-coil regions overlapping with an ezrin-radixin-moesin (ERM)-like domain followed by a serine/proline- and a lysine-rich region (2, 8, 20–22). Three recent studies have described the essential ORAI-activating region within the ERM domain, termed SOAR (Stim ORAI-activating region) (23), OASF (ORAI-activating small fragment) (24), and CAD (CRAC-activating domain) (25), including the second coiled coil domain and the following ∼55 amino acids. We and others have provided evidence that store depletion leads to a dynamic coupling of STIM1 to ORAI1 (26–28) that is mediated by a direct interaction of the STIM1 C terminus with ORAI1 C terminus probably involving the putative coiled-coil domain in the latter (27).Furthermore, different groups have proven that the C terminus of STIM1 is sufficient to activate CRAC as well as ORAI1 channels independent of store depletion (22–25, 27, 29). We have identified that OASF-(233–474) or shorter fragments exhibit further enhanced coupling to ORAI1 resulting in 3-fold increased constitutive Ca²⁺ currents. A STIM1 fragment containing an additional cluster of anionic amino acids C-terminal to position 474 displays weaker interaction with ORAI1 as well as reduced Ca²⁺ current comparable with that mediated by wild-type STIM1 C terminus. Hence, we have suggested that these 11 amino acids (474–485) act in a modulatory manner onto ORAI1; however, their detailed mechanistic impact within the STIM1/ORAI1 signaling machinery has remained so far unclear.In this study, we focused on the impact of this negative cluster on fast inactivation of STIM1-mediated ORAI Ca²⁺ currents. Lis et al. (30) have shown that all three ORAI homologues display distinct inactivation profiles, where ORAI2 and ORAI3 show a much more pronounced fast inactivation than ORAI1. Moreover, it has been reported (31) that different expression levels of STIM1 to ORAI1 affect the properties of CRAC current inactivation. Yamashita et al. (32) have demonstrated a linkage between the selectivity filter of ORAI1 and its Ca²⁺-dependent fast inactivation. Here we provide evidence that a cluster of acidic residues within the C terminus of STIM1 is involved in the fast inactivation of ORAI1 and further promotes that of ORAI3 and native CRAC currents.