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]     

Characterization of Mycobacterium leprae RecA Intein, a LAGLIDADG Homing Endonuclease, Reveals a Unique Mode of DNA Binding, Helical Distortion, and Cleavage Compared with a Canonical LAGLIDADG Homing Endonuclease

Mycobacterium leprae, which has undergone reductive evolution leaving behind a minimal set of essential genes, has retained intervening sequences in four of its genes implicating a vital role for them in the survival of the leprosy bacillus. A single in-frame intervening sequence has been found embedded within its recA gene. Comparison of the M. leprae recA intervening sequence with the known intervening sequences indicated that it has the consensus amino acid sequence necessary for being a LAGLIDADG-type homing endonuclease. In light of massive gene decay and function loss in the leprosy bacillus, we sought to investigate whether its recA intervening sequence encodes a catalytically active homing endonuclease. Here we show that the purified M. leprae RecA intein (PI-MleI) binds to cognate DNA and displays endonuclease activity in the presence of alternative divalent cations, Mg²⁺ or Mn²⁺. A combination of approaches, including four complementary footprinting assays such as DNase I, copper-phenanthroline, methylation protection, and KMnO₄, enhancement of 2-aminopurine fluorescence, and mapping of the cleavage site revealed that PI-MleI binds to cognate DNA flanking its insertion site, induces helical distortion at the cleavage site, and generates two staggered double strand breaks. Taken together, these results implicate that PI-MleI possesses a modular structure with separate domains for DNA target recognition and cleavage, each with distinct sequence preferences. From a biological standpoint, it is tempting to speculate that our findings have implications for understanding the evolution of the LAGLIDADG family of homing endonucleases.Mycobacterium leprae, a Gram-positive rod-shaped bacillus, mostly found in warm tropical countries, is the bacterium that causes leprosy in humans (1). The lack of understanding of the basic biology of M. leprae is believed to be the key factor for the failure of leprosy research to advance. The genome sequence of M. leprae contains 3.27 Mb and has an average G + C content of 57.8%, values much lower than the corresponding values for Mycobacterium tuberculosis, which are ∼4.41 Mb and 65.6% G + C, respectively (2). There are some 1500 genes that are common to both M. leprae and M. tuberculosis. The comparative genome analysis suggests that both species of mycobacteria are derived from a common ancestor and, at one stage, had gene pools of similar size. The downsizing of the M. tuberculosis genome from ∼4.41 to 3.27 Mb of M. leprae would account for the loss of some 1200 protein-coding sequences (1, 3). There is evidence that many of the genes that were present in the genome of M. leprae have truly been lost (1, 3). Comparative genomics of M. leprae with that of M. tuberculosis indicate that the former has undergone substantial downsizing, losing more than 2000 genes, thus suggesting an extreme case of reductive evolution in a microbial pathogen (1). With the availability of the M. leprae genome sequence, using functional genomics approaches, it is possible to identify the gene products, elucidate the mechanism of their action, and identify novel drug targets for rational design of new therapeutic regimens and drugs to treat leprosy.Eubacterial RecA proteins catalyze a set of biochemical reactions that are essential for homologous recombination, DNA repair, restoration of stalled replication forks, and SOS response (4–7). RecA protein and the process of homologous recombination, which is the main mechanism of genetic exchange, are evolutionarily conserved among a range of organisms (4, 7). Perhaps the most striking development in the field of RecA protein biology was the discovery of an in-frame insertion of an intein-coding sequence in the recA genes of M. tuberculosis and M. leprae (8, 9). In these organisms, RecA is synthesized as a large precursor, which undergoes protein splicing to excise the intein, and the two flanking domains called exteins are ligated together to generate a functionally active RecA protein (9, 10). The milieu in which RecA precursor undergoes splicing differs substantially between M. tuberculosis and M. leprae. M. leprae RecA precursor (79 kDa) undergoes splicing only in mycobacterial species, whereas M. tuberculosis RecA precursor (85 kDa) is spliced efficiently in Escherichia coli as well (9–11). Intriguingly, M. tuberculosis and M. leprae RecA inteins differ greatly in their size, primary sequence, and location within the recA gene, thereby suggesting two independent origins during evolution (9). The occurrence of inteins in the obligate mycobacterial pathogens, M. tuberculosis, M. leprae, and Mycobacterium microti, suggested that RecA inteins might play a role in mycobacterial functions related to pathogenesis or virulence (9). Previously, we have shown that M. tuberculosis RecA intein (PI-MtuI),² which contains Walker A motif, displays dual target specificity in the presence of alternative cofactors in an ATP-dependent manner (12, 13).Since their discovery in Saccharomyces cerevisiae (14, 15), a large number of putative homing endonucleases have been found in a diverse range of proteins in all the three domains of life (16–19). The majority of inteins possess the protein splicing and homing endonuclease activities (18, 19). Homing endonucleases are a class of diverse rare-cutting enzymes that promote site-specific transposition of their encoding genetic elements by inflicting double-stranded DNA breaks via different cleavage mechanisms in alleles lacking these elements (18–23). In addition, these are characterized by their ability to bind long DNA target sites (14–40 bp), and their tolerance of minor sequence changes in their binding region. These have been divided into highly divergent subfamilies on the basis of conserved sequence and structural motifs as follows: LAGLIDADG, GIY-YIG, HNH, His-Cys box, and the more recently identified PD(D/E)XK families (18–24). LAGLIDADG homing enzymes, which include the largest family, contain one or two copies of the conserved dodecapeptide motif and utilize an extended protein-DNA interface covering up to 40 bp to acquire their necessary specificity (18–22). The LAGLIDADG sequence is a part of the conserved 10- or 12-residue sequence motif defining the family of LAGLIDADG-type homing endonucleases; therefore, it is designated as deca- or dodecapeptide motif (19).Comparison of the M. leprae recA intervening sequence with known intervening sequences indicated that it has the consensus amino acid sequence necessary for being a LAGLIDADG-type homing endonuclease (25, 26). In light of massive gene decay and function loss in the leprosy bacillus, and dissimilarities in size and primary structures among mycobacterial inteins, we sought to investigate whether M. leprae recA intervening sequence encodes a catalytically active homing endonuclease. In this study, we show that the purified M. leprae RecA intein (PI-MleI) binds to cognate DNA and displays endonuclease activity in the presence of alternative divalent cations Mg²⁺ or Mn²⁺. Furthermore, using a variety of approaches, we have mapped the positions of PI-MleI binding as well as cleavage in the cognate DNA, thus providing the most comprehensive analysis of PI-MleI. Taken together, these results suggest that PI-MleI possesses a modular structure with functionally separable domains for DNA target recognition and cleavage, each with distinct sequence preferences. These results provide insights into understanding the function and evolution of the family of LAGLIDADG homing endonucleases.     

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]     

Acetylome Analysis Reveals Diverse Functions of Lysine Acetylation in Mycobacterium tuberculosis

The lysine acetylation of proteins is a reversible post-translational modification that plays a critical regulatory role in both eukaryotes and prokaryotes. Mycobacterium tuberculosis is a facultative intracellular pathogen and the causative agent of tuberculosis. Increasing evidence shows that lysine acetylation may play an important role in the pathogenesis of M. tuberculosis. However, only a few acetylated proteins of M. tuberculosis are known, presenting a major obstacle to understanding the functional roles of reversible lysine acetylation in this pathogen. We performed a global acetylome analysis of M. tuberculosis H37Ra by combining protein/peptide prefractionation, antibody enrichment, and LC-MS/MS. In total, we identified 226 acetylation sites in 137 proteins of M. tuberculosis H37Ra. The identified acetylated proteins were functionally categorized into an interaction map and shown to be involved in various biological processes. Consistent with previous reports, a large proportion of the acetylation sites were present on proteins involved in glycolysis/gluconeogenesis, the citrate cycle, and fatty acid metabolism. A NAD⁺-dependent deacetylase (MRA_1161) deletion mutant of M. tuberculosis H37Ra was constructed and its characterization showed a different colony morphology, reduced biofilm formation, and increased tolerance of heat stress. Interestingly, lysine acetylation was found, for the first time, to block the immunogenicity of a peptide derived from a known immunogen, HspX, suggesting that lysine acetylation plays a regulatory role in immunogenicity. Our data provide the first global survey of lysine acetylation in M. tuberculosis. The dataset should be an important resource for the functional analysis of lysine acetylation in M. tuberculosis and facilitate the clarification of the entire metabolic networks of this life-threatening pathogen.Mycobacterium tuberculosis was responsible for 1.3 million deaths and 8.6 million new cases of tuberculosis (TB)¹ worldwide in 2012 (1). This global public health crisis remains a serious problem, with the emergence of drug-resistant M. tuberculosis, especially multidrug-resistant and extensively drug-resistant M. tuberculosis, and also the emergence of coinfections of TB and human immunodeficiency virus (2, 3). To counter the increasing threat of TB, it is critical to understand fundamental aspects of TB-related biology. Such studies will not only provide new drug targets for the design of novel therapeutic agents, but also facilitate the development of novel diagnostic tools and new vaccines.Acetylation is one of the important protein modifications and occurs both co- and post-translationally on the α-amino group at the N terminus of the protein, so-called “N-terminal acetylation,” or on the ε-amino group on the side chain of lysine (4). Lysine acetylation is one of the most common post-translational modifications to proteins in both eukaryotes and prokaryotes. As a dynamic and reversible process, protein acetylation plays important roles in many cellular physiological processes, including cell-cycle regulation and apoptosis, cell morphology (5), metabolic pathways (6–8), protein interactions (9), and enzymatic activity (8, 10). In recent years, great advances have been made in proteomic studies, and a large number of lysine-acetylated proteins have been identified in many eukaryotes, including human (5, 11, 12), rat (13), mouse (11), Drosophila (14), Arabidopsis (15, 16), Saccharomyces cerevisiae (17), and protozoans (18, 19). The global analysis of lysine acetylation has also been reported in bacteria, including Escherichia coli (20–22), Erwinia amylovora (23), Bacillus subtilis (24), and Salmonella enterica (6). These acetylome studies have generated large datasets of bacterial proteins acetylated on lysine residues and have demonstrated the diverse cellular functions of lysine acetylation in bacteria.Increasing evidence shows that protein acetylation occurs and plays an important regulatory role in mycobacteria (8, 25–31). For example, Lange et al. reported the N-terminal acetylation of early secreted antigenic target 6 (ESAT-6) protein (31). Rv1151c is reported to be an NAD⁺-dependent protein deacetylase in M. tuberculosis that deacetylates and thus regulates the activity of acetyl-CoA synthase (25, 32). Two cyclic adenosine monophosphate (cAMP)-binding proteins in M. smegmatis and M. tuberculosis (MSMEG_5458 and Rv0998, respectively) show similarity to the GNAT family of acetyltransferases and could acetylate a universal stress protein (USP, MSMEG_4207) (30). Subsequent structural studies revealed the fine mechanisms of how cAMP regulates the protein lysine acetyltransferase in mycobacteria (27, 28). Very recently, reversible lysine acetylation was shown to regulate the activity of several fatty acyl-CoA synthetases in M. tuberculosis (8, 26), and also to regulate acetate and propionate metabolism in M. smegmatis (8, 26). However, to the best of our knowledge, only a few acetylated proteins in M. tuberculosis have been identified, presenting a major obstacle to further understanding the regulatory roles of reversible lysine acetylation in this life-threatening pathogen.To fill this gap in our knowledge, we undertook a systematic study of the functional roles of lysine acetylation in M. tuberculosis. We performed an acetylomic analysis of M. tuberculosis H37Ra using high-accuracy MS combined with the identification of 226 unique lysine acetylation sites on 137 proteins. This set of M. tuberculosis proteins acetylated on lysine residues supports the emerging view that lysine acetylation is a general and fundamental regulatory process, and is not restricted to eukaryotes. It also opens the way for its detailed functional and evolutionary analysis of lysine acetylation in M. tuberculosis. The identified acetylated proteins that are involved in several important biological processes were functionally categorized into an interaction map. This is the first time that an interaction network of acetylated proteins in M. tuberculosis has been constructed, and should allow us to better understand the significance of acetylation in key cellular mechanisms in M. tuberculosis. To further explore the effects of lysine acetylation on the physiology of M. tuberculosis H37Ra, MRA_1161, the gene encoding the only known protein deacetylase in this bacterium, was deleted. The roles of MRA_1161 in the colony morphology, carbon source utilization, heat stress tolerance, and biofilm formation of M. tuberculosis were analyzed. The effect of lysine acetylation on the immunogenicity of a known immunogen, HspX, was also tested.     

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]     

In Vitro and In Vivo Biology of Recombinant Adenovirus Vectors with E1, E1/E2A,or E1/E4 Deleted

Isogenic, E3-deleted adenovirus vectors defective in E1, E1 and E2A, or E1 and E4 were generated in complementation cell lines expressing E1, E1 and E2A, or E1 and E4 and characterized in vitro and in vivo. In the absence of complementation, deletion of both E1 and E2A completely abolished expression of early and late viral genes, while deletion of E1 and E4 impaired expression of viral genes, although at a lower level than the E1/E2A deletion. The in vivo persistence of these three types of vectors was monitored in selected strains of mice with viral genomes devoid of transgenes to exclude any interference by immunogenic transgene-encoded products. Our studies showed no significant differences among the vectors in the short-term maintenance and long-term (4-month) persistence of viral DNA in liver and lung cells of immunocompetent and immunodeficient mice. Furthermore, all vectors induced similar antibody responses and comparable levels of adenovirus-specific cytotoxic T lymphocytes. These results suggest that in the absence of transgenes, the progressive deletion of the adenovirus genome does not extend the in vivo persistence of the transduced cells and does not reduce the antivirus immune response. In addition, our data confirm that, in the absence of transgene expression, mouse cellular immunity to viral antigens plays a minor role in the progressive elimination of the virus genome.Replication-deficient human adenoviruses (Ad) have been widely investigated as ex vivo and in vivo gene delivery systems for human gene therapy. The ability of these vectors to mediate the efficient expression of candidate therapeutic or vaccine genes in a variety of cell types, including postmitotic cells, is considered an advantage over other gene transfer vectors (3, 28, 49). However, the successful application of currently available E1-defective Ad vectors in human gene therapy has been hampered by the fact that transgene expression is only transient in vivo (2, 15, 16, 33, 36, 46). This short-lived in vivo expression of the transgene has been explained, at least in part, by the induction in vivo of cytotoxic immune responses to cells infected with the Ad vector. Studies with rodent systems have suggested that cytotoxic T lymphocytes (CTLs) directed against virus antigens synthesized de novo in the transduced tissues play a major role in eliminating cells containing the E1-deleted viral genome (56–58, 61). Consistent with the concept of cellular antiviral immunity, expression of transgenes is significantly extended in experimental rodent systems that are deficient in various components of the cellular immune system or that have been rendered immunocompromised by administration of pharmacological agents (2, 33, 37, 48, 60, 64).Based on the assumption that further reduction of viral antigen expression may lower the immune response and thus extend persistence of transgene expression, previous studies have investigated the consequences of deleting both E1 and an additional viral regulatory region, such as E2A or E4. The E2A region encodes a DNA binding protein (DBP) with specific affinity for single-stranded Ad DNA. The DNA binding function is essential for the initiation and elongation of viral DNA synthesis during the early phase of Ad infection. During the late phase of infection, DBP plays a central role in the activation of the major late promoter (MLP) (for a recent review, see reference 44). The E4 region, located at the right end of the viral genome, encodes several regulatory proteins with pleiotropic functions which are involved in the accumulation, splicing, and transport of early and late viral mRNAs, in DNA replication, and in virus particle assembly (reviewed in reference 44). The simultaneous deletion of E1 and E2A or of E1 and E4 should therefore further reduce the replication of the virus genome and the expression of early and late viral genes. Such multidefective vectors have been generated and tested in vitro and in vivo (9, 12, 17, 19–21, 23, 24, 26, 34, 40, 52, 53, 59, 62, 63). Recombinant vectors with E1 deleted and carrying an E2A temperature-sensitive mutation (E2Ats) have been shown in vitro to express much smaller amounts of virus proteins, leading to extended transgene expression in cotton rats and mice (19, 20, 24, 59). To eliminate the risks of reversion of the E2Ats point mutation to a wild-type phenotype, improved vectors with both E1 and E2A deleted were subsequently generated in complementation cell lines coexpressing E1 and E2A genes (26, 40, 63). In vitro analysis of human cells infected by these viruses demonstrated that the double deletion completely abolished viral DNA replication and late protein synthesis (26). Similarly, E1/E4-deleted vectors have been generated in various in vitro complementation systems and tested in vitro and in vivo (9, 17, 23, 45, 52, 53, 62). These studies showed that deletion of both E1 and E4 did indeed reduce significantly the expression of early and late virus proteins (17, 23), leading to a decreased anti-Ad host immune response (23), reduced hepatotoxicity (17, 23, 52), and improved in vivo persistence of the transduced liver cells (17, 23, 52).Interpretation of these results is difficult, however, since all tested E1- and E1/E4-deleted vectors encoded the bacterial β-galactosidase (βgal) marker, whose strong immunogenicity is known to influence the in vivo persistence of Ad-transduced cells (32, 37). Moreover, the results described above are not consistent with the conclusions from other studies showing, in various immunocompetent mouse models, that cellular immunity to Ad antigens has no detectable impact on the persistence of the transduced cells (37, 40, 50, 51). Furthermore, in contrast to results of earlier studies (19, 20, 59), Fang et al. (21) demonstrated that injection of E1-deleted/E2Ats vectors into immunocompetent mice and hemophilia B dogs did not lead to an improvement of the persistence of transgene expression compared to that with isogenic E1-deleted vectors. Similarly, Morral et al. (40) did not observe any difference in persistence of transgene expression in mice injected with either vectors deleted in E1 only or vectors deleted in both E1 and E2A. Finally, the demonstration that some E4-encoded products can modulate transgene expression (1, 17, 36a) makes the evaluation of E1- and E1/E4-deleted vectors even more complex when persistence of transgene expression is used for direct comparison of the in vivo persistence of cells transduced by the two types of vectors.The precise influence of the host immune response to viral antigens on the in vivo persistence of the transduced cells, and hence the impact of further deletions in the virus genome, therefore still remains unclear. To investigate these questions, we generated a set of isogenic vectors with single deletions (AdE1°) and double deletions (AdE1°E2A° and AdE1°E4°) and their corresponding complementation cell lines and compared the biologies and immunogenicities of these vectors in vitro and in vivo. To eliminate any possible influence of transgene-encoded products on the interpretation of the in vivo results, we used E1-, E1/E2A-, and E1/E4-deleted vectors with no transgenes.     

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]