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]     

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]     

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.     

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]     

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]     

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]     

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.     

Dominant-negative Inhibitors of the Clostridium perfringens ?-Toxin

The Clostridium perfringens ϵ-toxin is responsible for a severe, often lethal intoxication. In this study, we characterized dominant-negative inhibitors of the ϵ-toxin. Site-specific mutations were introduced into the gene encoding ϵ-toxin, and recombinant proteins were expressed in Escherichia coli. Paired cysteine substitutions were introduced at locations predicted to form a disulfide bond. One cysteine in each mutant was introduced into the membrane insertion domain of the toxin; the second cysteine was introduced into the protein backbone. Mutant proteins with cysteine substitutions at amino acid positions I51/A114 and at V56/F118 lacked detectable cytotoxic activity in a MDCK cell assay. Cytotoxic activity could be reconstituted in both mutant proteins by incubation with dithiothreitol, indicating that the lack of cytotoxic activity was attributable to the formation of a disulfide bond. Fluorescent labeling of the cysteines also indicated that the introduced cysteines participated in a disulfide bond. When equimolar mixtures of wild-type ϵ-toxin and mutant proteins were added to MDCK cells, the I51C/A114C and V56C/F118C mutant proteins each inhibited the activity of wild-type ϵ-toxin. Further analysis of the inhibitory activity of the I51C/A114C and V56C/F118C mutant proteins indicated that these proteins inhibit the ability of the active toxin to form stable oligomeric complexes in the context of MDCK cells. These results provide further insight into the properties of dominant-negative inhibitors of oligomeric pore-forming toxins and provide the basis for developing new therapeutics for treating intoxication by ϵ-toxin.The Clostridium perfringens ϵ-toxin is one of the most potent bacterial toxins (1, 2). The ϵ-toxin can lead to a fatal enterotoxemia characterized by widespread vascular permeability and edema in the heart, lungs, brain, and kidneys (3–6). The disease most frequently affects livestock animals, though the toxin may also affect humans (7–9). Because of its extreme potency and the possibility of intoxicating humans, the C. perfringens ϵ-toxin is considered a select agent by the United States Department of Health and Human Services. A vaccine currently is approved for veterinary use, though multiple immunizations are required to provide long-term immunity (10–13). There also is an antitoxin approved for veterinary use. However, in the event that an animal exhibits symptoms of intoxication by ϵ-toxin, it is typically too late for the current antitoxin to be effective, and use of the antitoxin is typically limited to prophylactic treatment of unvaccinated animals within a herd (14). There is no treatment currently approved for use in humans. Thus, alternative countermeasures are needed that inhibit the activity of the toxin.One alternative method of countering the cytotoxic activity of bacterial toxins is through dominant-negative inhibitors. Dominant-negative inhibitors are non-cytotoxic mutant forms of active toxins that are able to inhibit the activity of wild-type toxin when the two proteins are mixed together. Such dominant-negative inhibitors have been described for a diverse set of toxins, including Helicobacter pylori VacA (15–19), Bacillus anthracis anthrax toxin protective antigen (20–25), Bacillus thuringiensis Cry1Ab (26), and Escherichia coli ClyA cytotoxin (27). Like VacA, protective antigen, Cry1Ab, and ClyA, the ϵ-toxin assembles into oligomeric complexes containing multiple toxin monomers (28–30). In the case of VacA and protective antigen, the most extensively studied examples of toxins inhibited by dominant-negative mutants, the number of mutations that inactivate the toxins is substantially greater than the number of mutations that lead to a dominant-negative phenotype (16, 17, 24, 31, 32). Although many of the mutations leading to dominant-negative toxins are located within regions of the toxins that are believed to form the membrane insertion domain, some mutations that inactivate the toxins (but are not dominant-negative) also map within the predicted membrane insertion domains (24, 32). Thus, a deeper understanding of the nature of the dominant-negative phenotype is needed.In this study, we sought to generate dominant-negative mutants of the ϵ-toxin. We hypothesized that mutations within the membrane insertion domain of ϵ-toxin, particularly mutations that are expected to restrict movement of this domain, would lead to dominant-negative inhibitors. We expressed wild-type and site-specific mutants of the ϵ-toxin as recombinant proteins in E. coli. The recombinant proteins were purified, and cytotoxicity was assessed using an established cell culture assay. Using this approach, we identified mutant proteins that inhibited the activity of wild-type ϵ-toxin in vitro and determined the mechanism of inhibition.     

MqsR, a Crucial Regulator for Quorum Sensing and Biofilm Formation, Is a GCU-specific mRNA Interferase in Escherichia coli

The mqsR gene has been shown to be positively regulated by the quorum-sensing autoinducer AI-2, which in turn activates a two-component system, the qseB-qseC operon. This operon plays an important role in biofilm formation in Escherichia coli. However, its cellular function has remained unknown. Here, we found that 1 base downstream of mqsR there is a gene, ygiT, that is co-transcribed with mqsR. Induction of mqsR caused cell growth arrest, whereas ygiT co-induction recovered cell growth. We demonstrate that MqsR (98 amino acid residues), which has no homology to the well characterized mRNA interferase MazF, is a potent inhibitor of protein synthesis that functions by degrading cellular mRNAs. In vivo and in vitro primer extension experiments showed that MqsR is an mRNA interferase specifically cleaving mRNAs at GCU. The mRNA interferase activity of purified MqsR was inhibited by purified YgiT (131 residues). MqsR forms a stable 2:1 complex with YgiT, and the complex likely functions as a repressor for the mqsR-ygiT operon by specifically binding to two different palindromic sequences present in the 5′-untranslated region of this operon.It has been reported that quorum sensing is involved in biofilm formation (1–4). mqsR expression was found to be induced by 8-fold in biofilms (5) and also by the quorum-sensing signal autoinducer AI-2, which is a species-nonspecific signaling molecule produced by both Gram-negative and Gram-positive bacteria, including Escherichia coli (6). It has been reported that induction of mqsR activates a two-component system, the qseB-qseC operon, which is known to play an important role in biofilm formation (6). Thus, it has been proposed that MqsR (98 amino acid residues) is a regulator of biofilm formation because it activates qseB, which controls the flhDC expression required for motility and biofilm formation in E. coli (6). However, the cellular function of MqsR has remained unknown.Interestingly, all free-living bacteria examined to date contain a number of suicide or toxin genes in their genomes (7, 8). Many of these toxins are co-transcribed with their cognate antitoxins in an operon (termed toxin-antitoxin (TA)² operon) and form a stable complex in the cell, so their toxicity is subdued under normal growth conditions (9–11). However, the stability of antitoxins is substantially lower than that of their cognate toxins, so any stress causing cellular damage or growth inhibition that induces proteases alters the balance between toxin and antitoxin, leading to toxin release in the cell.To date, 16 (12) TA systems have been reported on the E. coli genome, including relB-relE (13, 14), chpBI-chpBK (15), mazEF (16–18), yefM-yoeB (19, 20), dinJ-yafQ (21, 22), hipBA and hicAB (23, 24), prlF-yhaV (25), and ybaJ-hha (26). Interestingly, all of these TA operons appear to use similar modes of regulation: the formation of complexes between antitoxins and their cognate toxins to neutralize toxin activity and the ability of TA complexes to autoregulate their expression. The cellular targets of some toxins have been identified. CcdB directly interacts with gyrase A and blocks DNA replication (27, 28). RelE, which by itself has no endoribonuclease activity, appears to act as a ribosome-associating factor that promotes mRNA cleavage at the ribosome A-site (13, 29, 30). PemK (31), ChpBK (15), and MazF (32) are unique among toxins because they target cellular mRNAs for degradation by functioning as sequence-specific endoribonucleases to effectively inhibit protein synthesis and thereby cell growth.MazF, ChpBK, and PemK have been characterized as sequence-specific endoribonucleases that cleave mRNA at the ACA, ACY (Y is U, A, or G), and UAH (H is C, A, or U) sequences, respectively. They are completely different from other known endoribonucleases such as RNases E, A, and T1, as these toxins function as protein synthesis inhibitors by interfering with the function of cellular mRNAs. It is well known that small RNAs, such as mRNA-interfering cRNA (33), microRNA (34), and small interfering RNA (35), interfere with the function of specific RNAs. These small RNAs bind to specific mRNAs to inhibit their expression. Ribozymes also act on their target RNAs specifically and interfere with their function (36). Therefore, MazF, ChpBK, and PemK homologs form a novel endoribonuclease family that exhibits a new mRNA-interfering mechanism by cleaving mRNAs at specific sequences. Thus, they have been termed “mRNA interferases” (2).During our search for TA systems on the E. coli genome, we found that the mqsR gene is co-transcribed with a downstream gene, ygiT. These two genes appear to function as a TA system, as their size is small (98 residues for MqsR and 131 residues for YgiT) and their respective open reading frames are separated by 1 bp. In this study, we demonstrate that MqsR-YgiT is a new E. coli TA system consisting of a toxin, MqsR, and an antitoxin, YgiT. Moreover, we identify MqsR as a novel mRNA interferase that does not exhibit homology to MazF. This toxin cleaves RNA at GCU sequences in vivo and in vitro. The implication of this finding as to how this mRNA interferase is involved in cell physiology and biofilm formation will be discussed.