The sequence comparison index: A novel method for comparing proteins and proteomes

Historically morphological features were used as the primary means to classify organisms. However, the age of molecular genetics has allowed us to approach this field from the perspective of the organism's genetic code. Early work used highly conserved sequences, such as ribosomal RNA. The increasing number of complete genomes in the public data repositories provides the opportunity to look not only at a single gene, but at organisms' entire parts list. ^ Here the Sequence Comparison Index (SCI) and the Organism Comparison Index (OCI), algorithms and methods to compare proteins and proteomes, are presented. The complete proteomes of 104 sequenced organisms were compared. Over 280 million full Smith-Waterman alignments were performed on sequence pairs which had a reasonable expectation of being related. From these alignments a whole proteome phylogenetic tree was constructed. This method was also used to compare the small subunit (SSU) rRNA from each organism and a tree constructed from these results. The SSU rRNA tree by the SCI/OCI method looks very much like accepted SSU rRNA trees from sources such as the Ribosomal Database Project, thus validating the method. The SCI/OCI proteome tree showed a number of small but significant differences when compared to the SSU rRNA tree and proteome trees constructed by other methods. Horizontal gene transfer does not appear to affect the SCI/OCI trees until the transferred genes make up a large portion of the proteome. ^ As part of this work, the Database of Related Local Alignments (DaRLA) was created and contains over 81 million rows of sequence alignment information. DaRLA, while primarily used to build the whole proteome trees, can also be applied shared gene content analysis, gene order analysis, and creating individual protein trees. ^ Finally, the standard BLAST method for analyzing shared gene content was compared to the SCI method using 4 spirochetes. The SCI system performed flawlessly, finding all proteins from one organism against itself and finding all the ribosomal proteins between organisms. The BLAST system missed some proteins from its respective organism and failed to detect small ribosomal proteins between organisms. ^

The yeast STM1 protein binds G*G multiplex nucleic acids in vitro and associates with ribosomes in vivo

The formation of triple helical, or triplex DNA has been suggested to occur in several cellular processes such as transcription, replication, and recombination. Our laboratory previously found proteins in HeLa nuclear extracts and in S. cerevisiae whole cell extracts that avidly bound a Purine-motif (Pu) triplex probe in gel shift assays, or EMSA. In order to identify a triplex DNA-binding protein, we used conventional and affinity chromatography to purify the major Pu triplex-binding protein in yeast. Peptide microsequencing and data base searches identified this protein as the product of the STM1 gene. Confirmation that Stm1p is a Pu triplex-binding protein was obtained by EMSA using both recombinant Stm1p and whole cell extracts from stm1Δ yeast. Stm1p had previously been identified as G4p2, a G-quartet DNA- and RNA-binding protein. To study the cellular role and identify the nucleic acid ligand of Stm1p in vivo, we introduced an HA epitope at either the N- or C-terminus of Stm1p and performed immunoprecipitations with the HA.11 mAb. Using peptide microsequencing and Northern analysis, we positively identified a subset of both large and small subunit ribosomal proteins and all four rRNAs as associating with Stm1p. DNase I treatment did not affect the association of Stm1p with ribosomal components, but RNase A treatment abolished the association with all ribosomal proteins and RNA, suggesting this association is RNA-dependent. Sucrose gradient fractionation followed by Western and EMSA analysis confirmed that Stm1p associates with intact 80S monosomes, but not polysomes. The presence of additional, unidentified RNA in the Stm1p-immunoprecipitate, and the absence of tRNAs and elongation factors suggests that Stm1p binds RNA and could be involved in the regulation of translation. Immunofluorescence microscopy data showed Stm1p to be located throughout the cytoplasm, with a specific movement to the bud during the G2 phase of the cell cycle. A dramatically flocculent, large cell phenotype is observed when Stm1p has a C-terminal HA tag in a protease-deficient strain background. When STM1 is deleted in this background, the same phenotype is not observed and the deletion yeast grow very slowly compared to the wild-type. These data suggest that STM1 is not essential, but plays a role in cell growth by interacting with an RNP complex that may contain G*G multiplex RNA. ^

The role of specific regions of ribosomal RNAs in translation termination

The ribosome is a molecular machine that produces proteins in a cell. It consists of RNAs (rRNAs) and proteins. The rRNAs have been implicated in various aspects of protein biosynthesis supporting the idea that they function directly in translation. In this study the direct involvement of rRNA in translation termination was hypothesized and both genetic and biochemical strategies were designed to test this hypothesis. As a result, several regions of rRNAs from both ribosomal subunits were implicated in termination. More specifically, the mutation G1093A in an RNA of the large subunit (23S rRNA) and the mutation C1054A in the small subunit RNA (16S rRNA) of the Escherichia coli ribosome, were shown to affect the binding of the proteins that drive termination, RF1 and RF2. These mutations also caused defects in catalysis of peptidyl-tRNA hydrolysis, the last step of termination. Furthermore, the mutations affected the function of RF2 to a greater extent than that of RF1, a striking result considering the similarity of the RFs. The major defect in RF2 function was consistent with in vivo characteristics of the mutants and can be explained by the inability of the mutant rRNA sites to activate the hydrolytic center, that is the catalytic site for peptidyl-tRNA hydrolysis. Consistent with this explanation is the possibility of a direct interaction between the G1093-region (domain II of 23S rRNA) and the hydrolytic center (most likely domains IV–VI of 23S rRNA). To test that interaction hypothesis selections were performed for mutations in domains IV–VI that compensated for the growth defects caused by G1093A. Several compensatory mutations were isolated which not only restored growth in the presence of G1093A but also appeared to compensate for the termination defects caused by the G1093A. Therefore these results provided genetic evidence for an intramolecular interaction that might lead to peptidyl-tRNA hydrolysis. Finally, a new approach to the study of rRNA involvement in termination was designed. By screening a library of rRNA fragments, a fragment of the 23S rRNA (nt 74-136) was identified that caused readthrough of UGA. The antisense RNA fragment produced a similar effect. The data implicated the corresponding segment of intact 23S rRNA in termination. ^