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.     

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]     

Rearrangements in the Relative Orientation of Cytoplasmic Domains Induced by a Membrane-anchored Protein Mediate Modulations in Kv Channel Gating

Interdomain interactions between intracellular N and C termini have been described for various K⁺ channels, including the voltage-gated Kv2.1, and suggested to affect channel gating. However, no channel regulatory protein directly affecting N/C interactions has been demonstrated. Most Kv2.1 channel interactions with regulatory factors occur at its C terminus. The vesicular SNARE that is also present at a high concentration in the neuronal plasma membrane, VAMP2, is the only protein documented to affect Kv2.1 gating by binding to its N terminus. As its binding target has been mapped near a site implicated in Kv2.1 N/C interactions, we hypothesized that VAMP2 binding to the N terminus requires concomitant conformational changes in the C terminus, which wraps around the N terminus from the outside, to give VAMP2 access. Here, we first determined that the Kv2.1 N terminus, although crucial, is not sufficient to convey functional interaction with VAMP2, and that, concomitant to its binding to the “docking loop” at the Kv2.1 N terminus, VAMP2 binds to the proximal part of the Kv2.1 C terminus, C1a. Next, using computational biology approaches (ab initio modeling, docking, and molecular dynamics simulations) supported by molecular biology, biochemical, electrophysiological, and fluorescence resonance energy transfer analyses, we mapped the interaction sites on both VAMP2 and Kv2.1 and found that this interaction is accompanied by rearrangements in the relative orientation of Kv2.1 cytoplasmic domains. We propose that VAMP2 modulates Kv2.1 inactivation by interfering with the interaction between the docking loop and C1a, a mechanism for gating regulation that may pertain also to other Kv channels.Interdomain interactions between intracellularly located N and C termini have been described for various K⁺ channels, including inwardly rectifying Kir2.3 and Kir6.2 (1, 2), small conductance Ca²⁺-activated (hSK3) (3), and voltage-gated Kv2.1 (4) and Kv4.1 (5) channels. In the case of Kv2.1, two modes of interaction have been proposed: an association of the distal part of Kv2.1 C terminus (termed CTA domain; amino acids (aa) 741–853)⁴ with aa 67 and 75 of the Kv2.1 N terminus (4); or an association between the proximal part of the Kv2.1 C terminus (aa 444–477) and the predicted loop structure (aa 55–71) in the N-terminal T1 domain (6). In addition, involvement of the S4-S5 linker in this interaction has been suggested (7). Although these studies propose two different C-terminal sites, they indicate a specific loop in the N terminus of Kv2.1 (6, 8), which could be functionally related to the Shaker and Shal docking loops in the lateral part of their T1 domains (9, 10). These latter loops are responsible for the subfamily-specific association with β-subunits (Kvβ and KChIP, respectively). Further, the interaction between the N and C cytoplasmic termini (N/C interaction) of Kv2.1 has been shown to be dynamic and voltage-dependent and to involve structural rearrangements between these domains, which could affect both activation and inactivation gating of the channel (4, 6, 7). These rearrangements can be clearly detected with fluorescence resonance energy transfer (FRET) (11). A similar N/C interaction has been shown to affect gating of the closely related Kv4.1 channel (5, 12).It is conceivable that the specific packaging of Kv2.1 cytoplasmic termini (a relatively long C terminus (>400 aa) wrapping the N terminus (<190 aa) from the outside (4)) not only supports multiple interactions between the termini but also reflects the fact that most of the interactions of the channel with intracellular and membrane-bound regulatory factors occur at the C terminus, including channel phosphorylation (13–15), clustering through a unique proxinal restriction and clustering signal (16), and protein-protein interactions with both the plasma membrane SNAREs, syntaxin 1A and SNAP-25 (17–19), and the MiRP2 (KCNE3) peptide (20). For the Kv2.1 N terminus, on the other hand, there are only two examples of protein-protein interactions: a transient association with KChAP (21), which does not affect channel function; and an interaction with the vesicular SNARE partner VAMP2 (synaptobrevin 2), which is also present at a high concentration in the neuronal plasma membrane and enhances channel inactivation (8). Specifically, VAMP2 has been shown to associate with the extension of a docking loop in the lateral part of the T1 domain (8) near the site of interaction with the C terminus (4, 6). Thus, it is reasonable to hypothesize that interaction with VAMP2 will affect the N/C interaction, similar to proton-mediated Kir2.3 (1) and Kir1.1 (22) N/C interactions or the ATP-dependent Kir6.2 (2) N/C interaction. To date, no protein molecule that directly affects N/C interactions in a K⁺ channel has been demonstrated. Because VAMP2 was the first protein documented to affect Kv2.1 channel gating by binding to a specific N-terminal site, which is probably masked by the C terminus, we have put forward the idea that its interaction with the Kv2.1 N terminus requires conformational changes in the C terminus that will enable its access to the N terminus.Here we endeavored to gain a mechanistic and structural understanding of the Kv2.1-VAMP2 interaction. Based on our evidence, we propose that VAMP2 modulates Kv2.1 gating by interfering with the Kv2.1 cytoplasmic N/C interaction.     

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.     

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]     

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.     

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.     

Conformational Properties of ��-PrP

Prion propagation involves a conformational transition of the cellular form of prion protein (PrP^C) to a disease-specific isomer (PrP^Sc), shifting from a predominantly α-helical conformation to one dominated by β-sheet structure. This conformational transition is of critical importance in understanding the molecular basis for prion disease. Here, we elucidate the conformational properties of a disulfide-reduced fragment of human PrP spanning residues 91–231 under acidic conditions, using a combination of heteronuclear NMR, analytical ultracentrifugation, and circular dichroism. We find that this form of the protein, which similarly to PrP^Sc, is a potent inhibitor of the 26 S proteasome, assembles into soluble oligomers that have significant β-sheet content. The monomeric precursor to these oligomers exhibits many of the characteristics of a molten globule intermediate with some helical character in regions that form helices I and III in the PrP^C conformation, whereas helix II exhibits little evidence for adopting a helical conformation, suggesting that this region is a likely source of interaction within the initial phases of the transformation to a β-rich conformation. This precursor state is almost as compact as the folded PrP^C structure and, as it assembles, only residues 126–227 are immobilized within the oligomeric structure, leaving the remainder in a mobile, random-coil state.Prion diseases, such as Creutzfeldt-Jacob and Gerstmann-Sträussler-Scheinker in humans, scrapie in sheep, and bovine spongiform encephalopathy in cattle, are fatal neurological disorders associated with the deposition of an abnormally folded form of a host-encoded glycoprotein, prion (PrP)² (1). These diseases may be inherited, arise sporadically, or be acquired through the transmission of an infectious agent (2, 3). The disease-associated form of the protein, termed the scrapie form or PrP^Sc, differs from the normal cellular form (PrP^C) through a conformational change, resulting in a significant increase in the β-sheet content and protease resistance of the protein (3, 4). PrP^C, in contrast, consists of a predominantly α-helical structured domain and an unstructured N-terminal domain, which is capable of binding a number of divalent metals (5–12). A single disulfide bond links two of the main α-helices and forms an integral part of the core of the structured domain (13, 14).According to the protein-only hypothesis (15), the infectious agent is composed of a conformational isomer of PrP (16) that is able to convert other isoforms to the infectious isomer in an autocatalytic manner. Despite numerous studies, little is known about the mechanism of conversion of PrP^C to PrP^Sc. The most coherent and general model proposed thus far is that PrP^C fluctuates between the dominant native state and minor conformations, one or a set of which can self-associate in an ordered manner to produce a stable supramolecular structure composed of misfolded PrP monomers (3, 17). This stable, oligomeric species can then bind to, and stabilize, rare non-native monomer conformations that are structurally complementary. In this manner, new monomeric chains are recruited and the system can propagate.In view of the above model, considerable effort has been devoted to generating and characterizing alternative, possibly PrP^Sc-like, conformations in the hope of identifying common properties or features that facilitate the formation of amyloid oligomers. This has been accomplished either through PrP^Sc-dependent conversion reactions (18–20) or through conversion of PrP^C in the absence of a PrP^Sc template (21–25). The latter approach, using mainly disulfide-oxidized recombinant PrP, has generated a wide range of novel conformations formed under non-physiological conditions where the native state is relatively destabilized. These conformations have ranged from near-native (14, 26, 27), to those that display significant β-sheet content (21, 23, 28–33). The majority of these latter species have shown a high propensity for aggregation, although not all are on-pathway to the formation of amyloid. Many of these non-native states also display some of the characteristics of PrP^Sc, such as increased β-sheet content, protease resistance, and a propensity for oligomerization (28, 29, 31) and some have been claimed to be associated with the disease process (34).One such PrP folding intermediate, termed β-PrP, differs from the majority of studied PrP intermediate states in that it is formed by refolding the PrP molecule from the native α-helical conformation (here termed α-PrP), at acidic pH in a reduced state, with the disulfide bond broken (22, 35). Although no covalent differences between the PrP^C and PrP^Sc have been consistently identified to date, the role of the disulfide bond in prion propagation remains disputed (25, 36–39). β-PrP is rich in β-sheet structure (22, 35), and displays many of the characteristics of a PrP^Sc-like precursor molecule, such as partial resistance to proteinase K digestion, and the ability to form amyloid fibrils in the presence of physiological concentrations of salts (40).The β-PrP species previously characterized, spanning residues 91–231 of PrP, was soluble at low ionic strength buffers and monomeric, according to elution volume on gel filtration (22). NMR analysis showed that it displayed radically different spectra to those of α-PrP, with considerably fewer observable peaks and markedly reduced chemical shift dispersion. Data from circular dichroism experiments showed that fixed side chain (tertiary) interactions were lost, in contrast to the well defined β-sheet secondary structure, and thus in conjunction with the NMR data, indicated that β-PrP possessed a number of characteristics associated with a “molten globule” folding intermediate (22). Such states have been proposed to be important in amyloid and fibril formation (41). Indeed, antibodies raised against β-PrP (e.g. ICSM33) are capable of recognizing native PrP^Sc (but not PrP^C) (42–44). Subsequently, a related study examining the role of the disulfide bond in PrP folding confirmed that a monomeric molten globule-like form of PrP was formed on refolding the disulfide-reduced protein at acidic pH, but reported that, under their conditions, the circular dichroism response interpreted as β-sheet structure was associated with protein oligomerization (45). Indeed, atomic force microscopy on oligomeric full-length β-PrP (residues 23–231) shows small, round particles, showing that it is capable of formation of oligomers without forming fibrils (35). Notably, however, salt-induced oligomeric β-PrP has been shown to be a potent inhibitor of the 26 S proteasome, in a similar manner to PrP^Sc (46). Impairment of the ubiquitin-proteasome system in vivo has been linked to prion neuropathology in prion-infected mice (46).Although the global properties of several PrP intermediate states have been determined (30–32, 35), no information on their conformational properties on a sequence-specific basis has been obtained. Their conformational properties are considered important, as the elucidation of the chain conformation may provide information on the way in which these chains pack in the assembly process, and also potentially provide clues on the mechanism of amyloid assembly and the phenomenon of prion strains. As the conformational fluctuations and heterogeneity of molten globule states give rise to broad NMR spectra that preclude direct observation of their conformational properties by NMR (47–50), here we use denaturant titration experiments to determine the conformational properties of β-PrP, through the population of the unfolded state that is visible by NMR. In addition, we use circular dichroism and analytical ultracentrifugation to examine the global structural properties, and the distribution of multimeric species that are formed from β-PrP.     

Splitting the BLOSUM Score into Numbers of Biological Significance

Mathematical tools developed in the context of Shannon information theory were used to analyze the meaning of the BLOSUM score, which was split into three components termed as the BLOSUM spectrum (or BLOSpectrum). These relate respectively to the sequence convergence (the stochastic similarity of the two protein sequences), to the background frequency divergence (typicality of the amino acid probability distribution in each sequence), and to the target frequency divergence (compliance of the amino acid variations between the two sequences to the protein model implicit in the BLOCKS database). This treatment sharpens the protein sequence comparison, providing a rationale for the biological significance of the obtained score, and helps to identify weakly related sequences. Moreover, the BLOSpectrum can guide the choice of the most appropriate scoring matrix, tailoring it to the evolutionary divergence associated with the two sequences, or indicate if a compositionally adjusted matrix could perform better.[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]     

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]     

Regulation of Microtubule Dynamic Instability in Vitro by Differentially Phosphorylated Stathmin

Stathmin is an important regulator of microtubule polymerization and dynamics. When unphosphorylated it destabilizes microtubules in two ways, by reducing the microtubule polymer mass through sequestration of soluble tubulin into an assembly-incompetent T2S complex (two α:β tubulin dimers per molecule of stathmin), and by increasing the switching frequency (catastrophe frequency) from growth to shortening at plus and minus ends by binding directly to the microtubules. Phosphorylation of stathmin on one or more of its four serine residues (Ser¹⁶, Ser²⁵, Ser³⁸, and Ser⁶³) reduces its microtubule-destabilizing activity. However, the effects of phosphorylation of the individual serine residues of stathmin on microtubule dynamic instability have not been investigated systematically. Here we analyzed the effects of stathmin singly phosphorylated at Ser¹⁶ or Ser⁶³, and doubly phosphorylated at Ser²⁵ and Ser³⁸, on its ability to modulate microtubule dynamic instability at steady-state in vitro. Phosphorylation at either Ser¹⁶ or Ser⁶³ strongly reduced or abolished the ability of stathmin to bind to and sequester soluble tubulin and its ability to act as a catastrophe factor by directly binding to the microtubules. In contrast, double phosphorylation of Ser²⁵ and Ser³⁸ did not affect the binding of stathmin to tubulin or microtubules or its catastrophe-promoting activity. Our results indicate that the effects of stathmin on dynamic instability are strongly but differently attenuated by phosphorylation at Ser¹⁶ and Ser⁶³ and support the hypothesis that selective targeting by Ser¹⁶-specific or Ser⁶³-specific kinases provides complimentary mechanisms for regulating microtubule function.Stathmin is an 18-kDa ubiquitously expressed microtubule-destabilizing phosphoprotein whose activity is modulated by phosphorylation of its four serine residues, Ser¹⁶, Ser²⁵, Ser³⁸, and Ser⁶³ (1–7). Several classes of kinases have been identified that phosphorylate stathmin, including kinases associated with cell growth and differentiation such as members of the mitogen-activated protein kinase (MAPK)2 family, cAMP-dependent protein kinase (1–5, 8–11), and kinases associated with cell cycle regulation such as cyclin-dependent kinase 1 (3, 12–14). Phosphorylation of stathmin is required for cell cycle progression through mitosis and for proper assembly/function of the mitotic spindle (3, 13–16). Inhibition of stathmin phosphorylation produces strong mitotic phenotypes characterized by disassembly and disorganization of mitotic spindles and abnormal chromosome distributions (3, 13–14).Stathmin is known to destabilize microtubules in two ways. One is by binding to soluble tubulin and forming a stable complex that cannot polymerize into microtubules, consisting of one molecule of stathmin and two molecules of tubulin (T2S complex) (17–24). Addition of stathmin to microtubules in equilibrium with soluble tubulin results in sequestration of the tubulin and a reduction in the level of microtubule polymer (17–18, 22, 25–28). In addition to reducing the amount of assembled polymer, tubulin sequestration by stathmin has been shown to increase the switching frequency at microtubule plus ends from growth to shortening (called the catastrophe frequency) as the microtubules relax to a new steady state (17, 29). The second way is by binding directly to microtubules (27–30). The direct binding of stathmin to microtubules increases the catastrophe frequency at both ends of the microtubules and considerably more strongly at minus ends than at plus ends (27). Consistent with its strong catastrophe-promoting activity at minus ends, stathmin increases the treadmilling rate of steady-state microtubules in vitro (27). These results have led to the suggestion that stathmin might be an important cellular regulator of minus-end microtubule dynamics (27).Phosphorylation of stathmin diminishes its ability to regulate microtubule polymerization (3, 14, 25–26). Phosphorylation of Ser¹⁶ or Ser⁶³ appears to be more critical than phosphorylation of Ser²⁵ and Ser³⁸ for the ability of stathmin to bind to soluble tubulin and to inhibit microtubule assembly in vitro (3, 25). Inhibition of stathmin phosphorylation induces defects in spindle assembly and organization (3, 14) suggesting that not only soluble tubulin-microtubule levels are regulated by phosphorylation of stathmin, but the dynamics of microtubules could also be regulated in a phosphorylation-dependent manner.It is not known how phosphorylation at any of the four serine residues of stathmin affects its ability to regulate microtubule dynamics and, specifically, its ability to increase the catastrophe frequency at plus and minus ends due to its direct interaction with microtubules. Thus, we determined the effects of stathmin individually phosphorylated at either Ser¹⁶ or Ser⁶³ and doubly phosphorylated at both Ser²⁵ and Ser³⁸ on dynamic instability at plus and minus ends in vitro at microtubule polymer steady state and physiological pH (pH 7.2). We find that phosphorylation of Ser¹⁶ strongly reduces the direct catastrophe-promoting activity of stathmin at plus ends and abolishes it at minus ends, whereas phosphorylation of Ser⁶³ abolishes the activity at both ends. The effects of phosphorylation of individual serines correlated well with stathmin''s reduced abilities to form stable T2S complexes, to inhibit microtubule polymerization, and to bind to microtubules. In contrast, double phosphorylation of Ser²⁵ and Ser³⁸ did not alter the ability of stathmin to modulate dynamic instability at the microtubule ends, its ability to form a stable T2S complex, or its ability to bind to microtubules. The data further support the hypotheses that phosphorylation of stathmin on either Ser¹⁶ or Ser⁶³ plays a critical role in regulating microtubule polymerization and dynamics in cells.     

Identification of Protein Domains That Control Proton and Calcium Sensitivity of ASIC1a

The acid-sensing ion channels (ASICs) open in response to extracellular acidic pH, and individual subunits display differential sensitivity to protons and calcium. ASIC1a acts as a high affinity proton sensor, whereas ASIC2a requires substantially greater proton concentrations to activate. Using chimeras composed of ASIC1a and ASIC2a, we determined that two regions of the extracellular domain (residues 87–197 and 323–431) specify the high affinity proton response of ASIC1a. These two regions appear to undergo intersubunit interactions within the multimeric channel to specify proton sensitivity. Single amino acid mutations revealed that amino acids around Asp³⁵⁷ play a prominent role in determining the pH dose response of ASIC1a. Within the same region, mutation F352L abolished PcTx1 modulation of ASIC1a. Surprisingly, we determined that another area of the extracellular domain was required for calcium-dependent regulation of ASIC1a activation, and this region functioned independently of high affinity proton sensing. These results indicate that specific regions play overlapping roles in pH-dependent gating and PcTx1-dependent modulation of ASIC1a activity, whereas a distinct region determines the calcium dependence of ASIC1a activation.The acid-sensing ion channels (ASICs)³ are proton-gated ion channels expressed in neurons throughout the central and peripheral nervous system (1–3). ASICs are activated by extracellular acidosis, and protons act as ligands triggering channel opening (4). Disruption of the accn2 gene (which encodes ASIC1) dramatically reduces proton-gated currents in central neurons and alters a variety of behaviors, including fear, learning, and memory (5, 6). ASIC1 also contributes to neuronal damage and death during the prolonged acidosis accompanying cerebral ischemia (7). Specifically, mice with disruptions in the accn2 gene display 60% smaller lesion size compared with normal mice in models of stroke (8). Application of PcTx1, a venom peptide that prevents ASIC1a activation, is similarly neuroprotective, even when administered hours after injury (8, 9). Thus, ASIC1a represents a novel pharmacological target for the prevention of neuronal death following stroke.Mammals have four genes encoding ASICs (accn1 to -4) that encode at least six different ASIC subunits (1–3, 10). Like all members of the DEG/ENaC family, individual ASIC subunits have two transmembrane regions separated by a large cysteine-rich extracellular region. Three ASIC subunits associate to form homomeric or heteromeric channels with distinct biophysical characteristics (11–14). Specifically, ASIC1a homomeric channels activate at pH values much closer to neutral pH compared with ASIC2a homomeric channels. The high affinity proton sensitivity of ASIC1a plays a prominent role in acidosis-induced neuronal death, and modulators that alter the pH dose response of ASIC1a affect neuronal sensitivity to prolonged acidosis (8, 9, 15). For example, the neuroprotective venom peptide PcTx1 increases the proton sensitivity of the ASIC1a channel, allowing the channel to desensitize at neutral pH and become unresponsive to subsequent acidic shifts in pH (16, 17). The large extracellular region of ASIC1a is thought to be the site of proton/modulator interaction and governs the characteristics of channel gating (10, 11, 18). However, the exact molecular mechanisms defining ASIC1a activation and the protein domains that are responsible for the apparent proton sensitivity of ASIC1a remain unclear. Here, we used chimeras containing specific regions from both ASIC1a and ASIC2a to identify the specific protein regions that confer high affinity proton sensing, PcTx1 sensitivity, and calcium modulation to ASIC1a.     

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.