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]     

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]     

F-Actin Dynamics in Neurospora crassa

This study demonstrates the utility of Lifeact for the investigation of actin dynamics in Neurospora crassa and also represents the first report of simultaneous live-cell imaging of the actin and microtubule cytoskeletons in filamentous fungi. Lifeact is a 17-amino-acid peptide derived from the nonessential Saccharomyces cerevisiae actin-binding protein Abp140p. Fused to green fluorescent protein (GFP) or red fluorescent protein (TagRFP), Lifeact allowed live-cell imaging of actin patches, cables, and rings in N. crassa without interfering with cellular functions. Actin cables and patches localized to sites of active growth during the establishment and maintenance of cell polarity in germ tubes and conidial anastomosis tubes (CATs). Recurrent phases of formation and retrograde movement of complex arrays of actin cables were observed at growing tips of germ tubes and CATs. Two populations of actin patches exhibiting slow and fast movement were distinguished, and rapid (1.2 μm/s) saltatory transport of patches along cables was observed. Actin cables accumulated and subsequently condensed into actin rings associated with septum formation. F-actin organization was markedly different in the tip regions of mature hyphae and in germ tubes. Only mature hyphae displayed a subapical collar of actin patches and a concentration of F-actin within the core of the Spitzenkörper. Coexpression of Lifeact-TagRFP and β-tubulin–GFP revealed distinct but interrelated localization patterns of F-actin and microtubules during the initiation and maintenance of tip growth.Actins are highly conserved proteins found in all eukaryotes and have an enormous variety of cellular roles. The monomeric form (globular actin, or G-actin) can self-assemble, with the aid of numerous actin-binding proteins (ABPs), into microfilaments (filamentous actin, or F-actin), which, together with microtubules, form the two major components of the fungal cytoskeleton. Numerous pharmacological and genetic studies of fungi have demonstrated crucial roles for F-actin in cell polarity, exocytosis, endocytosis, cytokinesis, and organelle movement (6, 7, 20, 34, 35, 51, 52, 59). Phalloidin staining, immunofluorescent labeling, and fluorescent-protein (FP)-based live-cell imaging have revealed three distinct subpopulations of F-actin-containing structures in fungi: patches, cables, and rings (1, 14, 28, 34, 60, 63, 64). Actin patches are associated with the plasma membrane and represent an accumulation of F-actin around endocytic vesicles (3, 26, 57). Actin cables are bundles of actin filaments stabilized with cross-linking proteins, such as tropomyosins and fimbrin, and are assembled by formins at sites of active growth, where they form tracks for myosin V-dependent polarized secretion and organelle transport (10, 16, 17, 27, 38, 47, 48). Cables, unlike patches, are absolutely required for polarized growth in the budding yeast Saccharomyces cerevisiae (34, 38). Contractile actomyosin rings are essential for cytokinesis in budding yeast, whereas in filamentous fungi, actin rings are less well studied but are known to be involved in septum formation (20, 28, 34, 39, 40).Actin cables and patches have been particularly well studied in budding yeast. However, there are likely to be important differences between F-actin architecture and dynamics in budding yeast and those in filamentous fungi, as budding yeasts display only a short period of polarized growth during bud formation, which is followed by isotropic growth over the bud surface (10). Sustained polarized growth during hyphal morphogenesis is a defining feature of filamentous fungi (21), making them attractive models for studying the roles of the actin cytoskeleton in cell polarization, tip growth, and organelle transport.In Neurospora crassa and other filamentous fungi, disruption of the actin cytoskeleton leads to rapid tip swelling, which indicates perturbation of polarized tip growth, demonstrating a critical role for F-actin in targeted secretion to particular sites on the plasma membrane (7, 22, 29, 56). Immunofluorescence studies of N. crassa have shown that F-actin localizes to hyphal tips as “clouds” and “plaques” (7, 54, 59). However, immunolabeling has failed to reveal actin cables in N. crassa and offers limited insights into F-actin dynamics. Live-cell imaging of F-actin architecture and dynamics has not been accomplished in N. crassa, yet it is expected to yield key insights into cell polarization, tip growth, and intracellular transport.We took advantage of a recently developed live-cell imaging probe for F-actin called Lifeact (43). Lifeact is a 17-amino-acid peptide derived from the N terminus of the budding yeast actin-binding protein Abp140 (5, 63) and has recently been demonstrated to be a universal live-cell imaging marker for F-actin in eukaryotes (43). Here, we report the successful application of fluorescent Lifeact fusion constructs for live-cell imaging of F-actin in N. crassa. We constructed two synthetic genes consisting of Lifeact fused to “synthetic” green fluorescent protein (sGFP) (S65T) (henceforth termed GFP) (12) or red fluorescent protein (TagRFP) (33) and expressed these constructs in various N. crassa strains. In all strain backgrounds, fluorescent Lifeact constructs clearly labeled actin patches, cables, and rings and revealed a direct association of F-actin structures with sites of cell polarization and active tip growth. Our results demonstrate the efficacy of Lifeact as a nontoxic live-cell imaging probe in N. crassa.     

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.     

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]     

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.