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]     

Sip1, a Novel RS Domain-Containing Protein Essential for Pre-mRNA Splicing

Previous studies have shown that protein-protein interactions among splicing factors may play an important role in pre-mRNA splicing. We report here identification and functional characterization of a new splicing factor, Sip1 (SC35-interacting protein 1). Sip1 was initially identified by virtue of its interaction with SC35, a splicing factor of the SR family. Sip1 interacts with not only several SR proteins but also with U1-70K and U2AF65, proteins associated with 5′ and 3′ splice sites, respectively. The predicted Sip1 sequence contains an arginine-serine-rich (RS) domain but does not have any known RNA-binding motifs, indicating that it is not a member of the SR family. Sip1 also contains a region with weak sequence similarity to the Drosophila splicing regulator suppressor of white apricot (SWAP). An essential role for Sip1 in pre-mRNA splicing was suggested by the observation that anti-Sip1 antibodies depleted splicing activity from HeLa nuclear extract. Purified recombinant Sip1 protein, but not other RS domain-containing proteins such as SC35, ASF/SF2, and U2AF65, restored the splicing activity of the Sip1-immunodepleted extract. Addition of U2AF65 protein further enhanced the splicing reconstitution by the Sip1 protein. Deficiency in the formation of both A and B splicing complexes in the Sip1-depleted nuclear extract indicates an important role of Sip1 in spliceosome assembly. Together, these results demonstrate that Sip1 is a novel RS domain-containing protein required for pre-mRNA splicing and that the functional role of Sip1 in splicing is distinct from those of known RS domain-containing splicing factors.Pre-mRNA splicing takes place in spliceosomes, the large RNA-protein complexes containing pre-mRNA, U1, U2, U4/6, and U5 small nuclear ribonucleoprotein particles (snRNPs), and a large number of accessory protein factors (for reviews, see references 21, 22, 37, 44, and 48). It is increasingly clear that the protein factors are important for pre-mRNA splicing and that studies of these factors are essential for further understanding of molecular mechanisms of pre-mRNA splicing.Most mammalian splicing factors have been identified by biochemical fractionation and purification (3, 15, 19, 31–36, 45, 69–71, 73), by using antibodies recognizing splicing factors (8, 9, 16, 17, 61, 66, 67, 74), and by sequence homology (25, 52, 74).Splicing factors containing arginine-serine-rich (RS) domains have emerged as important players in pre-mRNA splicing. These include members of the SR family, both subunits of U2 auxiliary factor (U2AF), and the U1 snRNP protein U1-70K (for reviews, see references 18, 41, and 59). Drosophila alternative splicing regulators transformer (Tra), transformer 2 (Tra2), and suppressor of white apricot (SWAP) also contain RS domains (20, 40, 42). RS domains in these proteins play important roles in pre-mRNA splicing (7, 71, 75), in nuclear localization of these splicing proteins (23, 40), and in protein-RNA interactions (56, 60, 64). Previous studies by us and others have demonstrated that one mechanism whereby SR proteins function in splicing is to mediate specific protein-protein interactions among spliceosomal components and between general splicing factors and alternative splicing regulators (1, 1a, 6, 10, 27, 63, 74, 77). Such protein-protein interactions may play critical roles in splice site recognition and association (for reviews, see references 4, 18, 37, 41, 47 and 59). Specific interactions among the splicing factors also suggest that it is possible to identify new splicing factors by their interactions with known splicing factors.Here we report identification of a new splicing factor, Sip1, by its interaction with the essential splicing factor SC35. The predicted Sip1 protein sequence contains an RS domain and a region with sequence similarity to the Drosophila splicing regulator, SWAP. We have expressed and purified recombinant Sip1 protein and raised polyclonal antibodies against the recombinant Sip1 protein. The anti-Sip1 antibodies specifically recognize a protein migrating at a molecular mass of approximately 210 kDa in HeLa nuclear extract. The anti-Sip1 antibodies sufficiently deplete Sip1 protein from the nuclear extract, and the Sip1-depleted extract is inactive in pre-mRNA splicing. Addition of recombinant Sip1 protein can partially restore splicing activity to the Sip1-depleted nuclear extract, indicating an essential role of Sip1 in pre-mRNA splicing. Other RS domain-containing proteins, including SC35, ASF/SF2, and U2AF65, cannot substitute for Sip1 in reconstituting splicing activity of the Sip1-depleted nuclear extract. However, addition of U2AF65 further increases splicing activity of Sip1-reconstituted nuclear extract, suggesting that there may be a functional interaction between Sip1 and U2AF65 in nuclear extract.     

Genetic Analysis of the nuo Locus,Which Encodes the Proton-Translocating NADH Dehydrogenase in Escherichia coli

Complex I (EC 1.6.99.3) of the bacterium Escherichia coli is considered to be the minimal form of the type I NADH dehydrogenase, the first enzyme complex in the respiratory chain. Because of its small size and relative simplicity, the E. coli enzyme has become a model used to identify and characterize the mechanism(s) by which cells regulate the synthesis and assembly of this large respiratory complex. To begin dissecting the processes by which E. coli cells regulate the expression of nuo and the assembly of complex I, we undertook a genetic analysis of the nuo locus, which encodes the 14 Nuo subunits comprising E. coli complex I. Here we present the results of studies, performed on an isogenic collection of nuo mutants, that focus on the physiological, biochemical, and molecular consequences caused by the lack of or defects in several Nuo subunits. In particular, we present evidence that NuoG, a peripheral subunit, is essential for complex I function and that it plays a role in the regulation of nuo expression and/or the assembly of complex I.

Complex I (NADH:ubiquinone oxidoreductase; EC 1.6.99.3), a type I NADH dehydrogenase that couples the oxidation of NADH to the generation of a proton motive force, is the first enzyme complex of the respiratory chain (2, 35, 47). The Escherichia coli enzyme, considered to be the minimal form of complex I, consists of 14 subunits instead of the 40 to 50 subunits associated with the homologous eukaryotic mitochondrial enzyme (17, 29, 30, 48–50). E. coli also possesses a second NADH dehydrogenase, NDH-II, which does not generate a proton motive force (31). E. coli complex I resembles eukaryotic complex I in many ways (16, 17, 30, 49): it performs the same enzymatic reaction and is sensitive to a number of the same inhibitors, it consists of subunits homologous to those found in all proton-translocating NADH:ubiquinone oxidoreductases studied thus far, it is comprised of a large number of subunits relative to the number that comprise other respiratory enzymes, and it contains flavin mononucleotide and FeS center prosthetic groups. Additionally, it possesses an L-shaped topology (14, 22) like that of its Neurospora crassa homolog (27), and it consists of distinct fragments or subcomplexes. Whereas eukaryotic complex I can be dissected into a peripheral arm and a membrane arm, the E. coli enzyme consists of three subcomplexes referred to as the peripheral, connecting, and membrane fragments (29) (Fig. ​(Fig.1A).1A). The subunit composition of these three fragments correlates approximately with the organization of the 14 structural genes (nuoA to nuoN) (49) of the nuo (for NADH:ubiquinone oxidoreductase) locus (Fig. ​(Fig.1B),1B), an organization that is conserved in several other bacteria, including Salmonella typhimurium (3), Paracoccus denitrificans (53), Rhodobacter capsulatus (12), and Thermus thermophilus (54). The 5′ half of the locus contains a promoter (nuoP), previously identified and located upstream of nuoA (8, 49), and the majority of genes that encode subunits homologous to the nucleus-encoded subunits of eukaryotic complex I and to subunits of the Alcaligenes eutrophus NAD-reducing hydrogenase (17, 29, 30, 49). In contrast, the 3′ half contains the majority of the genes that encode subunits homologous to the mitochondrion-encoded subunits of eukaryotic complex I and to subunits of the E. coli formate-hydrogen lyase complex (17, 29, 30, 49). Whereas the nuclear homologs NuoE, NuoF, and NuoG constitute the peripheral fragment (also referred to as the NADH dehydrogenase fragment [NDF]), the nuclear homologs NuoB, NuoC, NuoD, and NuoI constitute the connecting fragment. The mitochondrial homologs NuoA, NuoH, NuoJ, NuoK, NuoL, NuoM, and NuoN constitute the membrane fragment (29). E. coli complex I likely evolved by fusion of preexisting protein assemblies constituting modules for electron transfer and proton translocation (17–19, 30). Open in a separate windowFIG. 1Schematic of E. coli complex I and the nuo locus. Adapted with permission of the publisher (17, 29, 30, 49). (A) E. coli complex I is comprised of three distinct fragments: the peripheral (light gray), connecting (white), and membrane (dark gray) fragments (17, 29). The peripheral fragment (NDF) is comprised of the nuclear homologs NuoE, -F, and -G and exhibits NADH dehydrogenase activity that oxidizes NADH to NAD⁺; the connecting fragment is comprised of the nuclear homologs NuoB, -C, -D, and -I; and the membrane fragment is comprised of the mitochondrial homologs NuoA, -H, and -J to -N and catalyzes ubiquinone (Q) to its reduced form (QH₂). FMN, flavin mononucleotide. (B) The E. coli nuo locus encodes the 14 Nuo subunits that constitute complex I. The 5′ half of the locus contains a previously identified promoter (nuoP) and the majority of genes that encode the peripheral and connecting subunits (light gray and white, respectively). The 3′ half of the locus contains the majority of the genes encoding the membrane subunits (dark gray). The 3′ end of nuoG encodes a C-Terminal region (CTR) of the NuoG subunit (hatched).Because of its smaller size and relative simplicity, researchers recently have begun to utilize complex I of E. coli, and that of its close relative S. typhimurium, to identify and characterize the mechanism(s) by which cells regulate the synthesis and assembly of this large respiratory complex (3, 8, 46) and to investigate the diverse physiological consequences caused by defects in this enzyme (4, 6, 10, 40, 59). Such defects affect the ability of cells to perform chemotaxis (40), to grow on certain carbon sources (4, 6, 10, 40, 57), to survive stationary phase (59), to perform energy-dependent proteolysis (4), to regulate the expression of at least one gene (32), and to maintain virulence (5).To begin dissecting the processes by which E. coli cells regulate the expression of nuo and the assembly of complex I, we undertook a genetic analysis of the nuo locus. Here, we present the results of studies, performed on an isogenic collection of nuo mutants, that focus on the physiological, biochemical, and molecular consequences caused by the lack of or defects in several Nuo subunits. In particular, we present evidence that NuoG, a peripheral subunit, is essential for complex I function and that it plays a role in the regulation of nuo expression and/or the assembly of complex I.     

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]     

��-Defensins in Enteric Innate Immunity: FUNCTIONAL PANETH CELL ��-DEFENSINS IN MOUSE COLONIC LUMEN*

Paneth cells are a secretory epithelial lineage that release dense core granules rich in host defense peptides and proteins from the base of small intestinal crypts. Enteric α-defensins, termed cryptdins (Crps) in mice, are highly abundant in Paneth cell secretions and inherently resistant to proteolysis. Accordingly, we tested the hypothesis that enteric α-defensins of Paneth cell origin persist in a functional state in the mouse large bowel lumen. To test this idea, putative Crps purified from mouse distal colonic lumen were characterized biochemically and assayed in vitro for bactericidal peptide activities. The peptides comigrated with cryptdin control peptides in acid-urea-PAGE and SDS-PAGE, providing identification as putative Crps. Matrix-assisted laser desorption ionization time-of-flight mass spectrometry experiments showed that the molecular masses of the putative α-defensins matched those of the six most abundant known Crps, as well as N-terminally truncated forms of each, and that the peptides contain six Cys residues, consistent with identities as α-defensins. N-terminal sequencing definitively revealed peptides with N termini corresponding to full-length, (des-Leu)-truncated, and (des-Leu-Arg)-truncated N termini of Crps 1–4 and 6. Crps from mouse large bowel lumen were bactericidal in the low micromolar range. Thus, Paneth cell α-defensins secreted into the small intestinal lumen persist as intact and functional forms throughout the intestinal tract, suggesting that the peptides may mediate enteric innate immunity in the colonic lumen, far from their upstream point of secretion in small intestinal crypts.Antimicrobial peptides (AMPs)² are released by epithelial cells onto mucosal surfaces as effectors of innate immunity (1–5). In mammals, most AMPs derive from two major families, the cathelicidins and defensins (6). The defensins comprise the α-, β-, and θ-defensin subfamilies, which are defined by the presence of six cysteine residues paired in characteristic tridisulfide arrays (7). α-Defensins are highly abundant in two primary cell lineages: phagocytic leukocytes, primarily neutrophils, of myeloid origin and Paneth cells, which are secretory epithelial cells located at the base of the crypts of Lieberkühn in the small intestine (8–10). Neutrophil α-defensins are stored in azurophilic granules and contribute to non-oxidative microbial cell killing in phagolysosomes (11, 12), except in mice whose neutrophils lack defensins (13). In the small bowel, α-defensins and other host defense proteins (14–18) are released apically as components of Paneth cell secretory granules in response to cholinergic stimulation and after exposure to bacterial antigens (19). Therefore, the release of Paneth cell products into the crypt lumen is inferred to protect mitotically active crypt cells from colonization by potential pathogens and confer protection against enteric infection (7, 20, 21).Under normal, homeostatic conditions, Paneth cells are not found outside the small bowel, although they may appear ectopically in response to local inflammation throughout the gastrointestinal tract (22, 23). Paneth cell numbers increase progressively throughout the small intestine, occurring at highest numbers in the distal ileum (24). Mouse Paneth cells express numerous α-defensin isoforms, termed cryptdins (Crps) (25), that have broad spectrum antimicrobial activities (6, 26). Collectively, α-defensins constitute approximately seventy percent of the bactericidal peptide activity in mouse Paneth cell secretions (19), selectively killing bacteria by membrane-disruptive mechanisms (27–30). The role of Paneth cell α-defensins in gastrointestinal mucosal immunity is evident from studies of mice transgenic for human enteric α-defensin-5, HD-5, which are immune to infection by orally administered Salmonella enterica sv. typhimurium (S. typhimurium) (31).The biosynthesis of mature, bactericidal α-defensins from their inactive precursors requires activation by lineage-specific proteolytic convertases. In mouse Paneth cells, inactive ∼8.4-kDa Crp precursors are processed intracellularly into microbicidal ∼4-kDa Crps by specific cleavage events mediated by matrix metalloproteinase-7 (MMP-7) (32, 33). MMP-7 null mice exhibit increased susceptibility to systemic S. typhimurium infection and decreased clearance of orally administered non-invasive Escherichia coli (19, 32). Although the α-defensin proregions are sensitive to proteolysis, the mature, disulfide-stabilized peptides resist digestion by their converting enzymes in vitro, whether the convertase is MMP-7 (32), trypsin (34), or neutrophil serine proteinases (35). Because α-defensins resist proteolysis in vitro, we hypothesized that Paneth cell α-defensins resist degradation and remain in a functional state in the large bowel, a complex, hostile environment containing varied proteases of both host and microbial origin.Here, we report on the isolation and characterization of a population of enteric α-defensins from the mouse colonic lumen. Full-length and N-terminally truncated Paneth cell α-defensins were identified and are abundant in the distal large bowel lumen.