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]     

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.     

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]     

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]     

Ligand-specific Conformational Changes in the ��1 Glycine Receptor Ligand-binding Domain

Understanding the activation mechanism of Cys loop ion channel receptors is key to understanding their physiological and pharmacological properties under normal and pathological conditions. The ligand-binding domains of these receptors comprise inner and outer β-sheets and structural studies indicate that channel opening is accompanied by conformational rearrangements in both β-sheets. In an attempt to resolve ligand-dependent movements in the ligand-binding domain, we employed voltage-clamp fluorometry on α1 glycine receptors to compare changes mediated by the agonist, glycine, and by the antagonist, strychnine. Voltage-clamp fluorometry involves labeling introduced cysteines with environmentally sensitive fluorophores and inferring structural rearrangements from ligand-induced fluorescence changes. In the inner β-sheet, we labeled residues in loop 2 and in binding domain loops D and E. At each position, strychnine and glycine induced distinct maximal fluorescence responses. The pre-M1 domain responded similarly; at each of four labeled positions glycine produced a strong fluorescence signal, whereas strychnine did not. This suggests that glycine induces conformational changes in the inner β-sheet and pre-M1 domain that may be important for activation, desensitization, or both. In contrast, most labeled residues in loops C and F yielded fluorescence changes identical in magnitude for glycine and strychnine. A notable exception was H201C in loop C. This labeled residue responded differently to glycine and strychnine, thus underlining the importance of loop C in ligand discrimination. These results provide an important step toward mapping the domains crucial for ligand discrimination in the ligand-binding domain of glycine receptors and possibly other Cys loop receptors.Glycine receptor (GlyR)3 chloride channels are pentameric Cys loop receptors that mediate fast synaptic transmission in the nervous system (1, 2). This family also includes nicotinic acetylcholine receptors (nAChRs), γ-aminobutyric acid type A and type C receptors, and serotonin type 3 receptors. Individual subunits comprise a large ligand-binding domain (LBD) and a transmembrane domain consisting of four α-helices (M1–M4). The LBD consists of a 10-strand β-sandwich made of an inner β-sheet with six strands and an outer β-sheet with four strands (3). The ligand-binding site is situated at the interface of adjacent subunits and is formed by loops A–C from one subunit and loops D–F from the neighboring subunit (3).The activation mechanism of Cys loop receptors is currently the subject of intense investigation because it is key to understanding receptor function under normal and pathological conditions (4, 5). Based on structural analysis of Torpedo nAChRs, Unwin and colleagues (6, 7) originally proposed that agonist binding induced the inner β-sheet to rotate, whereas the outer β-sheet tilted slightly upwards with loop C clasping around the agonist. These movements were thought to be transmitted to the transmembrane domain via a differential movement of loop 2 (β1-β2) and loop 7 (β6-β7) (both part of the inner β-sheet) and the pre-M1 domain (which is linked via a β-strand to the loop C in the outer sheet). The idea of large loop C movements accompanying agonist binding is supported by structural and functional data (3, 8–13). However, a direct link between loop C movements and channel gating has proved more difficult to establish. Although computational modeling studies have suggested that this loop may be a major component of the channel opening mechanism (14–18), experimental support for this model is not definitive. Similarly, loop F is also thought to move upon ligand binding, although there is as yet no consensus as to whether these changes represent local or global conformational changes (11, 19–21). Recently, a comparison of crystal structures of bacterial Cys loop receptors in the closed and open states revealed that although both the inner and outer β-sheets exhibit different conformations in closed and open states, the pre-M1 domain remains virtually stationary (22, 23). It is therefore relevant to question whether loop C, loop F, and pre-M1 movements are essential for Cys loop receptor activation.Strychnine is a classical competitive antagonist of GlyRs (24, 25), and to date there is no evidence that it can produce LBD structural changes. In this study we use voltage-clamp fluorometry (VCF) to compare glycine- and strychnine-induced conformational changes in the GlyR loops 2, C, D, E, and F and the pre-M1 domain in an attempt to determine whether they signal ligand-binding events, local conformational changes, or conformational changes associated with receptor activation.In a typical VCF experiment, a domain of interest is labeled with an environmentally sensitive fluorophore, and current and fluorescence are monitored simultaneously during ligand application. VCF is ideally suited for identifying ligand-specific conformational changes because it can report on electrophysiologically silent conformational changes (26), such as those induced by antagonists. Indeed, VCF has recently provided valuable insights into the conformational rearrangements of various Cys loop receptors (19, 21, 27–33).     

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.     

Multipattern Consensus Regions in Multiple Aligned Protein Sequences and Their Segmentation

Decomposing a biological sequence into its functional regions is an important prerequisite to understand the molecule. Using the multiple alignments of the sequences, we evaluate a segmentation based on the type of statistical variation pattern from each of the aligned sites. To describe such a more general pattern, we introduce multipattern consensus regions as segmented regions based on conserved as well as interdependent patterns. Thus the proposed consensus region considers patterns that are statistically significant and extends a local neighborhood. To show its relevance in protein sequence analysis, a cancer suppressor gene called p53 is examined. The results show significant associations between the detected regions and tendency of mutations, location on the 3D structure, and cancer hereditable factors that can be inferred from human twin studies.[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]     

In Vivo Identification of Nuclear Factors Interacting with the Conserved Elements of Box C/D Small Nucleolar RNAs

The U16 small nucleolar RNA (snoRNA) is encoded by the third intron of the L1 (L4, according to the novel nomenclature) ribosomal protein gene of Xenopus laevis and originates from processing of the pre-mRNA in which it resides. The U16 snoRNA belongs to the box C/D snoRNA family, whose members are known to assemble in ribonucleoprotein particles (snoRNPs) containing the protein fibrillarin. We have utilized U16 snoRNA in order to characterize the factors that interact with the conserved elements common to the other members of the box C/D class. In this study, we have analyzed the in vivo assembly of U16 snoRNP particles in X. laevis oocytes and identified the proteins which interact with the RNA by label transfer after UV cross-linking. This analysis revealed two proteins, of 40- and 68-kDa apparent molecular size, which require intact boxes C and D together with the conserved 5′,3′-terminal stem for binding. Immunoprecipitation experiments showed that the p40 protein corresponds to fibrillarin, indicating that this protein is intimately associated with the RNA. We propose that fibrillarin and p68 represent the RNA-binding factors common to box C/D snoRNPs and that both proteins are essential for the assembly of snoRNP particles and the stabilization of the snoRNA.One of the most interesting recent findings related to ribosome biogenesis has been the identification of a large number of small RNAs localized in the nucleolus (snoRNAs). So far, more than 60 snoRNAs have been identified in vertebrates (17), and more than 30 have been identified in yeast (2). The total number of snoRNAs is not known, but it is likely to be close to 200 (33, 38). These snoRNAs, with the exception of the mitochondrial RNA processing (MRP) species (38), can be grouped into two major families on the basis of conserved structural and sequence elements. The first group includes molecules referred to as box C/D snoRNAs, whereas the second one comprises the species belonging to the box H/ACA family (2, 15).The two families differ in many aspects. The box C/D snoRNAs are functionally heterogeneous. Most of them function as antisense RNAs in site-specific ribose methylation of the pre-rRNA (1, 10, 17, 26); a minority have been shown to play a direct role in pre-rRNA processing in both yeast and metazoan cells (11, 21). The box C/D snoRNAs play their role by means of unusually long (up to 21 contiguous nucleotides) regions of complementarity to highly conserved sequences of 28S and 18S rRNAs (1). In contrast, several members of the H/ACA RNA family have been shown to direct site-specific isomerization of uridines into pseudouridines and to display shorter regions of complementarity to rRNA (14, 24). Mutational analysis suggests that H/ACA snoRNAs can also play a role as antisense RNAs by base pairing with complementary regions on rRNA (15, 24).Another difference between the two families can be seen by comparison of secondary structures. A Y-shaped motif, where a 5′,3′-terminal stem adjoins the C and D conserved elements, has been proposed for many box C/D snoRNAs (16, 26, 40, 42), whereas box H/ACA snoRNAs have been proposed to fold into two conserved hairpin structures connected by a single-stranded hinge region, followed by a short 3′ tail (15).Despite these differences, analogies have been found in the roles played by the conserved box elements. Mutational analysis and competition experiments indicated that C/D and H/ACA boxes are required both for processing and stable accumulation of the mature snoRNA, suggesting that they represent binding sites for specific trans-acting factors (2, 3, 8, 15, 16, 28, 36, 41).All snoRNAs are associated with proteins to form specific ribonucleoparticles (snoRNPs). The study of these particles began only recently, and so far, very few aspects of their structure and biosynthesis have been clarified. The only detailed analysis performed was on the mammalian U3 (19) and the yeast snR30 (20) snoRNPs. Of the identified components, a few appear to be more general factors: fibrillarin, which was shown to be associated with C/D snoRNPs (3, 4, 8, 13, 28, 31, 39), and the nucleolar protein GAR1, which was found associated with H/ACA snoRNAs in yeast (20). Just as the study of small nuclear RNP (snRNP) particles was crucial to the understanding of the splicing process, a detailed structural and functional analysis of snoRNP particles will be essential to elucidate the complex process of ribosome biosynthesis.In this study, we have analyzed the snoRNP assembly of wild-type and mutant U16 snoRNAs by following the kinetics of complex formation in the in vivo system of the Xenopus laevis oocyte. By a UV cross-linking technique, we have identified two proteins, of 40- and 68-kDa apparent molecular mass, which require intact boxes C and D together with the terminal stem for their binding. The 40-kDa species is specifically recognized by fibrillarin antibodies, indicating that this protein is intimately associated with the RNA.     

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]     

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.