Unique Shine-Dalgarno Sequences in Cyanobacteria and Chloroplasts Reveal Evolutionary Differences in Their Translation Initiation.

Microorganisms require efficient translation to grow and replicate rapidly, and translation is often rate-limited by initiation. A prominent feature that facilitates translation initiation in bacteria is the Shine-Dalgarno (SD) sequence. However, there is much debate over its conservation in Cyanobacteria and in chloroplasts which presumably originated from endosymbiosis of ancient Cyanobacteria. Elucidating the utilization of SD sequences in Cyanobacteria and in chloroplasts is therefore important to understand whether 1) SD role in Cyanobacterial translation has been reduced prior to chloroplast endosymbiosis or 2) translation in Cyanobacteria and in plastid has been subjected to different evolutionary pressures. To test these alternatives, we employed genomic, proteomic, and transcriptomic data to trace differences in SD usage among Synechocystis species, Microcystis aeruginosa, cyanophages, Nicotiana tabacum chloroplast, and Arabidopsis thaliana chloroplast. We corrected their mis-annotated 16S rRNA 3' terminus using an RNA-Seq-based approach to determine their SD/anti-SD locational constraints using an improved measurement DtoStart. We found that cyanophages well-mimic Cyanobacteria in SD usage because both have been under the same selection pressure for SD-mediated initiation. Whereas chloroplasts lost this similarity because the need for SD-facilitated initiation has been reduced in plastids having much reduced genome size and different ribosomal proteins as a result of host-symbiont coevolution. Consequently, SD sequence significantly increases protein expression in Cyanobacteria but not in chloroplasts, and only Cyanobacterial genes compensate for a lack of SD sequence by having weaker secondary structures at the 5' UTR. Our results suggest different evolutionary pressures operate on translation initiation in Cyanobacteria and in chloroplast.

Staphylococcus aureus responds to allicin by global S-thioallylation - Role of the Brx/BSH/YpdA pathway and the disulfide reductase MerA to overcome allicin stress.

The prevalence of methicillin-resitant Staphylococcus aureus (MRSA) in hospitals and the community poses an increasing health burden, which requires the discovery of alternative antimicrobials. Allicin (diallyl thiosulfinate) from garlic exhibits broad-spectrum antimicrobial activity against many multidrug resistant bacteria. The thiol-reactive mode of action of allicin involves its S-thioallylations of low molecular weight (LMW) thiols and protein thiols. To investigate the mode of action and stress response caused by allicin in S. aureus, we analyzed the transcriptome signature, the targets for S-thioallylation in the proteome and the changes in the bacillithiol (BSH) redox potential (EBSH) under allicin stress. Allicin caused a strong thiol-specific oxidative and sulfur stress response and protein damage as revealed by the induction of the PerR, HypR, QsrR, MhqR, CstR, CtsR, HrcA and CymR regulons in the RNA-seq transcriptome. Allicin also interfered with metal and cell wall homeostasis and caused induction of the Zur, CsoR and GraRS regulons. Brx-roGFP2 biosensor measurements revealed a strongly increased EBSH under allicin stress. In the proteome, 57 proteins were identified with S-thioallylations under allicin treatment, including translation factors (EF-Tu, EF-Ts), metabolic and redox enzymes (AldA, GuaB, Tpx, KatA, BrxA, MsrB) as well as redox-sensitive MarR/SarA-family regulators (MgrA, SarA, SarH1, SarS). Phenotype and biochemical analyses revealed that BSH and the HypR-controlled disulfide reductase MerA are involved in allicin detoxification in S. aureus. The reversal of protein S-thioallylation was catalyzed by the Brx/BSH/YpdA pathway. Finally, the BSSB reductase YpdA was shown to use S-allylmercaptobacillithiol (BSSA) as substrate to regenerate BSH in S. aureus. In conclusion, allicin results in an oxidative shift of EBSH and protein S-thioallylation, which can be reversed by YpdA and the Brx/BSH/YpdA electron pathways in S. aureus to regenerate thiol homeostasis.

Development of Assay Systems for Amber Codon Decoding at the Steps of Initiation and Elongation in Mycobacteria.

Genetic analysis of the mechanism of protein synthesis in Gram-positive bacteria has remained largely unexplored because of the unavailability of appropriate in vivo assay systems. We developed chloramphenicol acetyltransferase (CAT)-based in vivo reporter systems to study translation initiation and elongation in Mycobacterium smegmatis The CAT reporters utilize specific decoding of amber codons by mutant initiator tRNA (i-tRNA, metU) molecules containing a CUA anticodon (metUCUA). The assay systems allow structure-function analyses of tRNAs without interfering with the cellular protein synthesis and function with or without the expression of heterologous GlnRS from Escherichia coli We show that despite their naturally occurring slow-growth phenotypes, the step of i-tRNA formylation is vital in translation initiation in mycobacteria and that formylation-deficient i-tRNA mutants (metUCUA/A1, metUCUA/G72, and metUCUA/G72G73) with a Watson-Crick base pair at the 1·72 position participate in elongation. In the absence of heterologous GlnRS expression, the mutant tRNAs are predominantly aminoacylated (glutamylated) by nondiscriminating GluRS. Acid urea gels show complete transamidation of the glutamylated metUCUA/G72G73 tRNA to its glutaminylated form (by GatCAB) in M. smegmatis In contrast, the glutamylated metUCUA/G72 tRNA did not show a detectable level of transamidation. Interestingly, the metUCUA/A1 mutant showed an intermediate activity of transamidation and accumulated in both glutamylated and glutaminylated forms. These observations suggest important roles for the discriminator base position and/or a weak Watson-Crick base pair at 1·72 for in vivo recognition of the glutamylated tRNAs by M. smegmatis GatCAB.IMPORTANCE Genetic analysis of the translational apparatus in Gram-positive bacteria has remained largely unexplored because of the unavailability of appropriate in vivo assay systems. We developed chloramphenicol acetyltransferase (CAT)-based reporters which utilize specific decoding of amber codons by mutant tRNAs at the steps of initiation and/or elongation to allow structure-function analysis of the translational machinery. We show that formylation of the initiator tRNA (i-tRNA) is crucial even for slow-growing bacteria and that i-tRNA mutants with a CUA anticodon are aminoacylated by nondiscriminating GluRS. The discriminator base position, and/or a weak Watson-Crick base pair at the top of the acceptor stem, provides important determinants for transamidation of the i-tRNA-attached Glu to Gln by the mycobacterial GatCAB.

Structure and dynamics study of translation initiation factor 1 from Staphylococcus aureus suggests its RNA binding mode.

Translation initiation, the rate-limiting step in the protein synthesis, is tightly regulated. As one of the translation initiation factors, translation initiation factor 1 (IF1) plays crucial roles not only in translation but also in many cellular processes that are important for genomic stability, such as the activity of RNA chaperones. Here, we characterize the RNA interactions and dynamics of IF1 from Staphylococcus aureus Mu50 (IF1Sa) by NMR spectroscopy. First, the NMR-derived solution structure of IF1Sa revealed that IF1Sa adopts an oligonucleotide/oligosaccharide binding (OB)-fold. Structural comparisons showed large deviations in the α-helix and the following loop, which are potential RNA-binding regions of the OB-fold, as well as differences in the electrostatic potential surface among bacterial IF1s. Second, the 15N NMR relaxation data for IF1Sa indicated the flexible nature of the α-helix and the following loop region of IF1Sa. Third, RNA-binding properties were studied using FP assays and NMR titrations. FP binding assays revealed that IF1Sa binds to RNAs with moderate affinity. In combination with the structural analysis, the NMR titration results revealed the RNA binding sites. Taken together, these results show that IF1Sa binds RNAs with moderate binding affinity via the residues that occupy the large surface area of its ß-barrel. These findings suggest that IF1Sa is likely to bind RNA in various conformations rather than only at a specific site and indicate that the flexible RNA binding mode of IF1Sa is necessary for its interaction with various RNA substrates.

The conserved GTPase LepA contributes mainly to translation initiation in Escherichia coli.

LepA is a paralog of EF-G found in all bacteria. Deletion of lepA confers no obvious growth defect in Escherichia coli, and the physiological role of LepA remains unknown. Here, we identify nine strains (ΔdksA, ΔmolR1, ΔrsgA, ΔtatB, ΔtonB, ΔtolR, ΔubiF, ΔubiG or ΔubiH) in which ΔlepA confers a synthetic growth phenotype. These strains are compromised for gene regulation, ribosome assembly, transport and/or respiration, indicating that LepA contributes to these functions in some way. We also use ribosome profiling to deduce the effects of LepA on translation. We find that loss of LepA alters the average ribosome density (ARD) for hundreds of mRNA coding regions in the cell, substantially reducing ARD in many cases. By contrast, only subtle and codon-specific changes in ribosome distribution along mRNA are seen. These data suggest that LepA contributes mainly to the initiation phase of translation. Consistent with this interpretation, the effect of LepA on ARD is related to the sequence of the Shine-Dalgarno region. Global perturbation of gene expression in the ΔlepA mutant likely explains most of its phenotypes.

A comprehensive analysis of the importance of translation initiation factors for Haloferax volcanii applying deletion and conditional depletion mutants.

Translation is an important step in gene expression. The initiation of translation is phylogenetically diverse, since currently five different initiation mechanisms are known. For bacteria the three initiation factors IF1 - IF3 are described in contrast to archaea and eukaryotes, which contain a considerably higher number of initiation factor genes. As eukaryotes and archaea use a non-overlapping set of initiation mechanisms, orthologous proteins of both domains do not necessarily fulfill the same function. The genome of Haloferax volcanii contains 14 annotated genes that encode (subunits of) initiation factors. To gain a comprehensive overview of the importance of these genes, it was attempted to construct single gene deletion mutants of all genes. In 9 cases single deletion mutants were successfully constructed, showing that the respective genes are not essential. In contrast, the genes encoding initiation factors aIF1, aIF2γ, aIF5A, aIF5B, and aIF6 were found to be essential. Factors aIF1A and aIF2ß are encoded by two orthologous genes in H. volcanii. Attempts to generate double mutants failed in both cases, indicating that also these factors are essential. A translatome analysis of one of the single aIF2ß deletion mutants revealed that the translational efficiency of the second ortholog was enhanced tenfold and thus the two proteins can replace one another. The phenotypes of the single deletion mutants also revealed that the two aIF1As and aIF2ßs have redundant but not identical functions. Remarkably, the gene encoding aIF2α, a subunit of aIF2 involved in initiator tRNA binding, could be deleted. However, the mutant had a severe growth defect under all tested conditions. Conditional depletion mutants were generated for the five essential genes. The phenotypes of deletion mutants and conditional depletion mutants were compared to that of the wild-type under various conditions, and growth characteristics are discussed.

Precise and reliable gene expression via standard transcription and translation initiation elements.

An inability to reliably predict quantitative behaviors for novel combinations of genetic elements limits the rational engineering of biological systems. We developed an expression cassette architecture for genetic elements controlling transcription and translation initiation in Escherichia coli: transcription elements encode a common mRNA start, and translation elements use an overlapping genetic motif found in many natural systems. We engineered libraries of constitutive and repressor-regulated promoters along with translation initiation elements following these definitions. We measured activity distributions for each library and selected elements that collectively resulted in expression across a 1,000-fold observed dynamic range. We studied all combinations of curated elements, demonstrating that arbitrary genes are reliably expressed to within twofold relative target expression windows with ∼93% reliability. We expect the genetic element definitions validated here can be collectively expanded to create collections of public-domain standard biological parts that support reliable forward engineering of gene expression at genome scales.

Bacteroides mobilizable and conjugative genetic elements: antibiotic resistance among clinical isolates

La conjugación es uno de los mecanismos más importantes de la transferencia horizontal de genes en procariotas, lo que lleva a la variación genética dentro de una especie y la adquisición de nuevos rasgos, como la resistencia a los antibióticos. Bacteroides es un anaerobio obligado y un patógeno oportunista importante. La resistencia a los antibióticos entre especies de Bacteroides está aumentando rápidamente, debido en gran parte a la difusión de los factores de transferencia de ADN (plásmidos y transposones) albergado por los miembros de este género. Los factores de transferencia se pueden dividir en dos clases, conjugativos y movilizables. Las especies de Bacteroides han presentado plásmidos de resistencia a los antibióticos, todos los cuales han sido intensamente estudiados. Estos plásmidos codifican de alto nivel de resistencia MLS conferida por un gen erm conservado. Se ha informado una observación interesante asociada a la transferencia de varios de estos tipos de elementos, todo lo cual confiere mayor resistencia y que aparecen en gran medida la eficiencia de transferencia después de la exposición a la tetraciclina. Muchos de los transposones conjugativos (CTn) en Bacteroides están relacionados con varios elementos genéticos (como CTnDOT, CnTnERL, NBU y otros). CTnDOT lleva un gen de resistencia a la tetraciclina, tetQ, y un gen de resistencia a la eritromicina, ermF. La resistencia a los medicamentos utilizados para tratar infecciones por Bacteroides, tales como la clindamicina, también ha ido en aumento. Estos elementos conjugativos han sido encontrados en los aislados clínicos de Bacteroides. Por lo tanto, la transferencia horizontal de genes posiblemente podría jugar un papel importante en la creciente incidencia de la resistencia bacteriana en este grupo(AU)

The conjugation is one of the most important mechanisms of horizontal gene transfer in prokaryotes, leading to genetic variation within a species and the acquisition of new traits, such as antibiotic resistance. Bacteroides is an obligate anaerobe of the colon and a significant opportunistic pathogen. Antibiotic resistance among Bacteroides spp. is rapidly increasing, largely due to the dissemination of DNA transfer factors (plasmids and transposons) harbored by members of this genus. Transfer factors can be divided into two classes, conjugative and mobilizable. Species of the intestinal Bacteroides have yielded different resistance plasmids, all of which have been intensely studied, the plasmids encode high-level MLS resistance conferred by a conserved erm gene. It has been reported an interesting observation associated with the transfer of several of these types of elements, all of which conferred Tcr and displayed greatly increased transfer efficiency following exposure to tetracycline. Many of the conjugative transposons (CTns) in Bacteroides are related to various genetic elements (such as CTnDOT, CTnERL, NBU and others). CTnDOT carries a tetracycline resistance gene, tetQ, and an erythromycin resistance gene, ermF. Resistance to drugs used to treat Bacteroides infections, such as clindamycin, has also been increasing. These conjugal elements have been found in Bacteroides clinical isolates. Thus, horizontal gene transfer could conceivably have played a role in the rising incidence of resistance in this bacterial group(AU)

Crystal structure of TTHA0061, an uncharacterized protein from Thermus thermophilus HB8, reveals a novel fold.

The crystal structure of an uncharacterized protein TTHA0061 from Thermus thermophilus HB8, was determined and refined to 1.8 A by a single wavelength anomalous dispersion (SAD) method. The structural analysis and comparison of TTHA0061 with other existing structures in the Protein Data Bank (PDB) revealed a novel fold, suggesting that this protein may belong to a translation initiation factor or ribosomal protein family. Differential scanning calorimetry analysis suggested that the thermostability of TTHA0061 increased at pH ranges of 5.8-6.2, perhaps due to the abundance of glutamic acid residues.

Translational control in cancer.

Remarkable progress has been made in defining a new understanding of the role of mRNA translation and protein synthesis in human cancer. Translational control is a crucial component of cancer development and progression, directing both global control of protein synthesis and selective translation of specific mRNAs that promote tumour cell survival, angiogenesis, transformation, invasion and metastasis. Translational control of cancer is multifaceted, involving alterations in translation factor levels and activities unique to different types of cancers, disease stages and the tumour microenvironment. Several clinical efforts are underway to target specific components of the translation apparatus or unique mRNA translation elements for cancer therapeutics.

Transient kinetics, fluorescence, and FRET in studies of initiation of translation in bacteria.

Initiation of mRNA translation in prokaryotes requires the small ribosomal subunit (30S), initiator fMet-tRNA(fMet), three initiation factors, IF1, IF2, and IF3, and the large ribosomal subunit (50S). During initiation, the 30S subunit, in a complex with IF3, binds mRNA, IF1, IF2.GTP, and fMet-tRNA(fMet) to form a 30S initiation complex which then recruits the 50S subunit to yield a 70S initiation complex, while the initiation factors are released. Here we describe a transient kinetic approach to study the timing of elemental steps of 30S initiation complex formation, 50S subunit joining, and the dissociation of the initiation factors from the 70S initiation complex. Labeling of ribosomal subunits, fMet-tRNA(fMet), mRNA, and initiation factors with fluorescent reporter groups allows for the direct observation of the formation or dissociation of complexes by monitoring changes in the fluorescence of single dyes or fluorescence resonance energy transfer (FRET) between two fluorophores. Subunit joining was monitored by light scattering or by FRET between dyes attached to the ribosomal subunits. The kinetics of chemical steps, that is, GTP hydrolysis by IF2 and peptide bond formation following the binding of aminoacyl-tRNA to the 70S initiation complex, were measured by the quench-flow technique. The methods described here are based on results obtained with initiation components from Escherichia coli but can be adopted for mechanistic studies of initiation in other prokaryotic or eukaryotic systems.

MED: a new non-supervised gene prediction algorithm for bacterial and archaeal genomes.

BACKGROUND: Despite a remarkable success in the computational prediction of genes in Bacteria and Archaea, a lack of comprehensive understanding of prokaryotic gene structures prevents from further elucidation of differences among genomes. It continues to be interesting to develop new ab initio algorithms which not only accurately predict genes, but also facilitate comparative studies of prokaryotic genomes. RESULTS: This paper describes a new prokaryotic genefinding algorithm based on a comprehensive statistical model of protein coding Open Reading Frames (ORFs) and Translation Initiation Sites (TISs). The former is based on a linguistic "Entropy Density Profile" (EDP) model of coding DNA sequence and the latter comprises several relevant features related to the translation initiation. They are combined to form a so-called Multivariate Entropy Distance (MED) algorithm, MED 2.0, that incorporates several strategies in the iterative program. The iterations enable us to develop a non-supervised learning process and to obtain a set of genome-specific parameters for the gene structure, before making the prediction of genes. CONCLUSION: Results of extensive tests show that MED 2.0 achieves a competitive high performance in the gene prediction for both 5' and 3' end matches, compared to the current best prokaryotic gene finders. The advantage of the MED 2.0 is particularly evident for GC-rich genomes and archaeal genomes. Furthermore, the genome-specific parameters given by MED 2.0 match with the current understanding of prokaryotic genomes and may serve as tools for comparative genomic studies. In particular, MED 2.0 is shown to reveal divergent translation initiation mechanisms in archaeal genomes while making a more accurate prediction of TISs compared to the existing gene finders and the current GenBank annotation.

RNA polymerase pausing regulates translation initiation by providing additional time for TRAP-RNA interaction.

RNA polymerase (RNAP) pause sites have been identified in several prokaryotic genes. Although the presumed biological function of RNAP pausing is to allow synchronization of RNAP position with regulatory factor binding and/or RNA folding, a direct causal link between pausing and changes in gene expression has been difficult to establish. RNAP pauses at two sites in the Bacillus subtilis trpEDCFBA operon leader. Pausing at U107 and U144 participates in transcription attenuation and trpE translation control mechanisms, respectively. Substitution of U144 caused a substantial pausing defect in vitro and in vivo. These mutations led to increased trp operon expression that was suppressed by overproduction of TRAP, indicating that pausing at U144 provides additional time for TRAP to bind to the nascent transcript and promote formation of an RNA structure that blocks translation of trpE. These results establish that pausing is capable of playing a role in regulating translation in bacteria.

Detection of prokaryotic promoters from the genomic distribution of hexanucleotide pairs.

BACKGROUND: In bacteria, sigma factors and other transcriptional regulatory proteins recognize DNA patterns upstream of their target genes and interact with RNA polymerase to control transcription. As a consequence of evolution, DNA sequences recognized by transcription factors are thought to be enriched in intergenic regions (IRs) and depleted from coding regions of prokaryotic genomes. RESULTS: In this work, we report that genomic distribution of transcription factors binding sites is biased towards IRs, and that this bias is conserved amongst bacterial species. We further take advantage of this observation to develop an algorithm that can efficiently identify promoter boxes by a distribution-dependent approach rather than a direct sequence comparison approach. This strategy, which can easily be combined with other methodologies, allowed the identification of promoter sequences in ten species and can be used with any annotated bacterial genome, with results that rival with current methodologies. Experimental validations of predicted promoters also support our approach. CONCLUSION: Considering that complete genomic sequences of over 1000 bacteria will soon be available and that little transcriptional information is available for most of them, our algorithm constitutes a promising tool for the prediction of promoter sequences. Importantly, our methodology could also be adapted to identify DNA sequences recognized by other regulatory proteins.

An unsupervised classification scheme for improving predictions of prokaryotic TIS.

BACKGROUND: Although it is not difficult for state-of-the-art gene finders to identify coding regions in prokaryotic genomes, exact prediction of the corresponding translation initiation sites (TIS) is still a challenging problem. Recently a number of post-processing tools have been proposed for improving the annotation of prokaryotic TIS. However, inherent difficulties of these approaches arise from the considerable variation of TIS characteristics across different species. Therefore prior assumptions about the properties of prokaryotic gene starts may cause suboptimal predictions for newly sequenced genomes with TIS signals differing from those of well-investigated genomes. RESULTS: We introduce a clustering algorithm for completely unsupervised scoring of potential TIS, based on positionally smoothed probability matrices. The algorithm requires an initial gene prediction and the genomic sequence of the organism to perform the reannotation. As compared with other methods for improving predictions of gene starts in bacterial genomes, our approach is not based on any specific assumptions about prokaryotic TIS. Despite the generality of the underlying algorithm, the prediction rate of our method is competitive on experimentally verified test data from E. coli and B. subtilis. Regarding genomes with high G+C content, in contrast to some previously proposed methods, our algorithm also provides good performance on P. aeruginosa, B. pseudomallei and R. solanacearum. CONCLUSION: On reliable test data we showed that our method provides good results in post-processing the predictions of the widely-used program GLIMMER. The underlying clustering algorithm is robust with respect to variations in the initial TIS annotation and does not require specific assumptions about prokaryotic gene starts. These features are particularly useful on genomes with high G+C content. The algorithm has been implemented in the tool "TICO" (TIs COrrector) which is publicly available from our web site.

The GTPase, CpgA(YloQ), a putative translation factor, is implicated in morphogenesis in Bacillus subtilis.

YloQ, from Bacillus subtilis, was identified previously as an essential nucleotide-binding protein of unknown function. YloQ was successfully over-expressed in Escherichia coli in soluble form. The purified protein displayed a low GTPase activity similar to that of other small bacterial GTPases such as Bex/Era. Based on the demonstrated GTPase activity and the unusual order of the yloQ G motifs, we now designate this protein as CpgA (circularly permuted GTPase). An unexpected property of this low abundance GTPase was the demonstration, using gel filtration and ultracentrifugation analysis, that the protein formed stable dimers, dependent upon the concentration of YloQ(CpgA), but independent of GTP. In order to investigate function, cpgA was placed under the control of the pspac promotor in the B. subtilis chromosome. When grown in E or Spizizen medium in the absence of IPTG, the rate of growth was significantly reduced. A large proportion of the cells exhibited a markedly perturbed morphology, with the formation of swollen, bent or 'curly' shapes. To confirm that this was specifically due to depleted CpgA a plasmid-borne cpgA under pxyl control was introduced. This restored normal cell shape and growth rate, even in the absence of IPTG, provided xylose was present. The crystal structure of CpgA(YloQ) suggests a role as a translation initiation factor and we discuss the possibility that CpgA is involved in the translation of a subset of proteins, including some required for shape maintenance.

Selenocysteine tRNA-specific elongation factor SelB is a structural chimaera of elongation and initiation factors.

In all three kingdoms of life, SelB is a specialized translation elongation factor responsible for the cotranslational incorporation of selenocysteine into proteins by recoding of a UGA stop codon in the presence of a downstream mRNA hairpin loop. Here, we present the X-ray structures of SelB from the archaeon Methanococcus maripaludis in the apo-, GDP- and GppNHp-bound form and use mutational analysis to investigate the role of individual amino acids in its aminoacyl-binding pocket. All three SelB structures reveal an EF-Tu:GTP-like domain arrangement. Upon binding of the GTP analogue GppNHp, a conformational change of the Switch 2 region in the GTPase domain leads to the exposure of SelB residues involved in clamping the 5' phosphate of the tRNA. A conserved extended loop in domain III of SelB may be responsible for specific interactions with tRNA(Sec) and act as a ruler for measuring the extra long acceptor arm. Domain IV of SelB adopts a beta barrel fold and is flexibly tethered to domain III. The overall domain arrangement of SelB resembles a 'chalice' observed so far only for initiation factor IF2/eIF5B. In our model of SelB bound to the ribosome, domain IV points towards the 3' mRNA entrance cleft ready to interact with the downstream secondary structure element.

Oligo kernels for datamining on biological sequences: a case study on prokaryotic translation initiation sites.

BACKGROUND: Kernel-based learning algorithms are among the most advanced machine learning methods and have been successfully applied to a variety of sequence classification tasks within the field of bioinformatics. Conventional kernels utilized so far do not provide an easy interpretation of the learnt representations in terms of positional and compositional variability of the underlying biological signals. RESULTS: We propose a kernel-based approach to datamining on biological sequences. With our method it is possible to model and analyze positional variability of oligomers of any length in a natural way. On one hand this is achieved by mapping the sequences to an intuitive but high-dimensional feature space, well-suited for interpretation of the learnt models. On the other hand, by means of the kernel trick we can provide a general learning algorithm for that high-dimensional representation because all required statistics can be computed without performing an explicit feature space mapping of the sequences. By introducing a kernel parameter that controls the degree of position-dependency, our feature space representation can be tailored to the characteristics of the biological problem at hand. A regularized learning scheme enables application even to biological problems for which only small sets of example sequences are available. Our approach includes a visualization method for transparent representation of characteristic sequence features. Thereby importance of features can be measured in terms of discriminative strength with respect to classification of the underlying sequences. To demonstrate and validate our concept on a biochemically well-defined case, we analyze E. coli translation initiation sites in order to show that we can find biologically relevant signals. For that case, our results clearly show that the Shine-Dalgarno sequence is the most important signal upstream a start codon. The variability in position and composition we found for that signal is in accordance with previous biological knowledge. We also find evidence for signals downstream of the start codon, previously introduced as transcriptional enhancers. These signals are mainly characterized by occurrences of adenine in a region of about 4 nucleotides next to the start codon. CONCLUSIONS: We showed that the oligo kernel can provide a valuable tool for the analysis of relevant signals in biological sequences. In the case of translation initiation sites we could clearly deduce the most discriminative motifs and their positional variation from example sequences. Attractive features of our approach are its flexibility with respect to oligomer length and position conservation. By means of these two parameters oligo kernels can easily be adapted to different biological problems.

Isolation and characterization of ribosomes and translation initiation factors from the gram-positive soil bacterium Streptomyces lividans.

A primer extension inhibition (toeprint) assay was developed using ribosomes and ribosomal subunits from Streptomyces lividans. This assay allowed the study of ribosome binding to streptomycete leaderless and leadered mRNA. Purified 30S subunits were unable to form a ternary complex on aph leaderless mRNA, whereas 70S ribosomes could form ternary complexes on this mRNA. 30S subunits formed ternary complexes on leadered aph and malE mRNA. The translation initiation factors (IF1, IF2, and IF3) from S. lividans were isolated and included in toeprint and filter binding assays with leadered and leaderless mRNA. Generally, the IFs reduced the toeprint signal on leadered mRNA; however, incubation of IF1 and IF2 with 30S subunits that had been washed under high-salt conditions promoted the formation of a ternary complex on aph leaderless mRNA. Our data suggest that, as reported for Escherichia coli, initiation complexes with leaderless mRNAs might use a novel pathway involving 70S ribosomes or 30S subunits bound by IF1 and IF2 but not IF3. Some mRNA-ribosome-initiator tRNA reactions that yielded weak or no toeprint signals still formed complexes in filter binding assays, suggesting the occurrence of interactions that are not stable in the toeprint assay.

The genetic core of the universal ancestor.

Molecular analysis of conserved sequences in the ribosomal RNAs of modern organisms reveals a three-domain phylogeny that converges in a universal ancestor for all life. We used the Clusters of Orthologous Groups database and information from published genomes to search for other universally conserved genes that have the same phylogenetic pattern as ribosomal RNA, and therefore constitute the ancestral genetic core of cells. Our analyses identified a small set of genes that can be traced back to the universal ancestor and have coevolved since that time. As indicated by earlier studies, almost all of these genes are involved with the transfer of genetic information, and most of them directly interact with the ribosome. Other universal genes have either undergone lateral transfer in the past, or have diverged so much in sequence that their distant past could not be resolved. The nature of the conserved genes suggests innovations that may have been essential to the divergence of the three domains of life. The analysis also identified several genes of unknown function with phylogenies that track with the ribosomal RNA genes. The products of these genes are likely to play fundamental roles in cellular processes.