Searching for drug targets in microbial genomes

Comparative analysis of the complete genome sequences of 10 bacterial pathogens available in the public databases offers the first insights into the drug discovery approaches of the near future. Genes that are conserved in different genomes often turn out to be essential, which makes them attractive targets for new broad-spectrum antibiotics. Subtractive genome analysis reveals the genes that are conserved in all or most of the pathogenic bacteria but not in eukaryotes; these are the most obvious candidates for drug targets. Species-specific genes, on the other hand, may offer the possibility to design drugs against a particular, narrow group of pathogens.     

The 64-kilodalton capsid protein homolog of Beet yellows virus is required for assembly of virion tails

The filamentous virion of the closterovirus Beet yellows virus (BYV) consists of a long body formed by the major capsid protein (CP) and a short tail composed of the minor capsid protein (CPm) and the virus-encoded Hsp70 homolog. By using nano-liquid chromatography-tandem mass spectrometry and biochemical analyses, we show here that the BYV 64-kDa protein (p64) is the fourth integral component of BYV virions. The N-terminal domain of p64 is exposed at the virion surface and is accessible to antibodies and mild trypsin digestion. In contrast, the C-terminal domain is embedded in the virion and is inaccessible to antibodies or trypsin. The C-terminal domain of p64 is shown to be homologous to CP and CPm. Mutation of the signature motifs of capsid proteins of filamentous RNA viruses in p64 results in the formation of tailless virions, which are unable to move from cell to cell. These results reveal the dual function of p64 in tail assembly and BYV motility and support the concept of the virion tail as a specialized device for BYV cell-to-cell movement.     

Remarkable interkingdom conservation of intron positions and massive,lineage-specific intron loss and gain in eukaryotic evolution

Sequencing of eukaryotic genomes allows one to address major evolutionary problems, such as the evolution of gene structure. We compared the intron positions in 684 orthologous gene sets from 8 complete genomes of animals, plants, fungi, and protists and constructed parsimonious scenarios of evolution of the exon-intron structure for the respective genes. Approximately one-third of the introns in the malaria parasite Plasmodium falciparum are shared with at least one crown group eukaryote; this number indicates that these introns have been conserved through >1.5 billion years of evolution that separate Plasmodium from the crown group. Paradoxically, humans share many more introns with the plant Arabidopsis thaliana than with the fly or nematode. The inferred evolutionary scenario holds that the common ancestor of Plasmodium and the crown group and, especially, the common ancestor of animals, plants, and fungi had numerous introns. Most of these ancestral introns, which are retained in the genomes of vertebrates and plants, have been lost in fungi, nematodes, arthropods, and probably Plasmodium. In addition, numerous introns have been inserted into vertebrate and plant genes, whereas, in other lineages, intron gain was much less prominent.     

Emergence of diverse biochemical activities in evolutionarily conserved structural scaffolds of proteins

Comparative analysis of numerous protein structures that have become available in the past few years, combined with genome comparison, has yielded new insights into the evolution of enzymes and their functions. In addition to the well-known diversification of substrate specificities, enzymes with several widespread catalytic folds, particularly the TIM barrel, the RRM-like domain and the double-stranded beta-helix (cupin) domain, have been extensively explored in 'reaction space', resulting in the evolution of numerous, diverse catalytic activities supported by the same structural scaffold. Common protein folds differ widely in the diversity of catalyzed reactions. The biochemical plasticity of a fold seems to hinge on the presence of a generic, symmetrical substrate-binding pocket as opposed to highly specialized binding sites.     

Scores of RINGS but no PHDs in ubiquitin signaling

Recently, it has been reported that PHD fingers of MEKK1 kinase and a family of viral and cellular membrane proteins have E3 ubiquitin ligase activity. Here we describe unique sequence and structural signatures that distinguish PHD fingers from RING fingers, which function primarily as E3 ubiquitin ligases, and demonstrate that the Zn-binding modules of the above proteins are distinct versions of the RING domain rather than PHD fingers. Thus, currently available data reveal extreme versatility of RINGs and their derivatives that function as E3 ubiquitin ligases but provide no evidence of this activity among PHD fingers whose principal function appears to involve specific protein-protein and possibly protein-DNA interactions in chromatin.     

Purifying and directional selection in overlapping prokaryotic genes

In overlapping genes, the same DNA sequence codes for two proteins using different reading frames. Analysis of overlapping genes can help in understanding the mode of evolution of a coding region from noncoding DNA. We identified 71 pairs of convergent genes, with overlapping 3' ends longer than 15 nucleotides, that are conserved in at least two prokaryotic genomes. Among the overlap regions, we observed a statistically significant bias towards the 123:132 phase (i.e. the second codon base in one gene facing the degenerate third position in the second gene). This phase ensures the least mutual constraint on nonconservative amino acid replacements in both overlapping coding sequences. The excess of this phase is compatible with directional (positive) selection acting on the overlapping coding regions. This could be a general evolutionary mode for genes emerging from noncoding sequences, in which the protein sequence has not been subject to selection.     

Microevolutionary genomics of bacteria

The availability of multiple complete genome sequences from the same species can facilitate attempts to systematically address basic questions in genome evolution. We refer to such efforts as "microevolutionary genomics". We report the results of comparative analyses of complete intraspecific genome (and proteome) sequences from four bacterial species--Chlamydophila pneumoniae, Escherichia coli, Helicobacter pylori and Neisseria meningitidis. Comparisons of average synonymous (K(s)) and nonsynonymous (K(a)) substitution rates were used to assess the influence of various biological factors on the rate of protein evolution. For example, E. coli experiences the most intense purifying selection of the species analyzed, and this may be due to the relatively larger population size of this species. In addition, essential genes were shown to be more evolutionarily conserved than nonessential genes in E. coli and duplicated genes have higher rates of evolution than unique genes for all species studied except C. pneumoniae. Different functional categories of genes were shown to evolve at significantly different rates emphasizing the role of category-specific functional constraints in determining evolutionary rates. Finally, functionally characterized genes tend to be conserved between strains, while uncharacterized genes are over-represented among the unique, strain-specific genes. This suggests the possibility that nonessential genes are responsible for driving the evolutionary diversification between strains.     

Identification of the catalytic sites of a papain-like cysteine proteinase of murine coronavirus.

The murine coronavirus mouse hepatitis virus gene 1 is expressed as a polyprotein, which is cleaved into multiple proteins posttranslationally. One of the proteins is p28, which represents the amino-terminal portion of the polyprotein and is presumably generated by the activity of an autoproteinase domain of the polyprotein (S. C. Baker, C. K. Shieh, L. H. Soe, M.-F. Chang, D. M. Vannier, and M. M. C. Lai, J. Virol. 63:3693-3699, 1989). In this study, the boundaries and the critical amino acid residues of this putative proteinase domain were characterized by deletion analysis and site-directed mutagenesis. Proteinase activity was monitored by examining the generation of p28 during in vitro translation in rabbit reticulocyte lysates. Deletion analysis defined the proteinase domain to be within the sequences encoded from the 3.6- to 4.4-kb region from the 5' end of the genome. A 0.7-kb region between the substrate (p28) and proteinase domain could be deleted without affecting the proteolytic cleavage. However, a larger deletion (1.6 kb) resulted in the loss of proteinase activity, suggesting the importance of spacing sequences between proteinase and substrate. Computer-assisted analysis of the amino acid sequence of the proteinase domain identified potential catalytic cysteine and histidine residues in a stretch of sequence distantly related to papain-like cysteine proteinases. The role of these putative catalytic residues in the proteinase activity was studied by site-specific mutagenesis. Mutations of Cys-1137 or His-1288 led to a complete loss of proteinase activity, implicating these residues as essential for the catalytic activity. In contrast, most mutations of His-1317 or Cys-1172 had no or only minor effects on proteinase activity. This study establishes that mouse hepatitis virus gene 1 encodes a proteinase domain, in the region from 3.6 to 4.4 kb from the 5' end of the genome, which resembles members of the papain family of cysteine proteinases and that this proteinase domain is responsible for the cleavage of the N-terminal peptide.