GenTHREADER: an efficient and reliable protein fold recognition method for genomic sequences

MODBASE, a database of annotated comparative protein structure models

MODBASE is a queryable database of annotated comparative protein structure models. The models are derived by MODPIPE, an automated modeling pipeline relying on the programs PSI-BLAST and MODELLER. The database currently contains 3D models for substantial portions of approximately 17 000 proteins from 10 complete genomes, including those of Caenorhabditis elegans, Saccharomyces cerevisiae and Escherichia coli, as well as all the available sequences from Arabidopsis thaliana and Homo sapiens. The database also includes fold assignments and alignments on which the models were based. In addition, special care is taken to assess the quality of the models. ModBase is accessible through a web interface at http://guitar.rockefeller.edu/modbase/

Beta propellers: structural rigidity and functional diversity

Recently solved structures and proposed models have helped to reveal the structural characteristics of the beta-propeller fold, as well as the features that contribute to its high rigidity and stability. Possible strategies for identifying beta-propeller proteins in newly characterised sequences are helping to overcome the problems of predicting the beta-propeller fold from amino acid sequences.

Mutations of acidic residues in RAG1 define the active site of the V(D)J recombinase

The RAG1 and RAG2 proteins collaborate to initiate V(D)J recombination by binding recombination signal sequences (RSSs) and making a double-strand break between the RSS and adjacent coding DNA. Like the reactions of their biochemical cousins, the bacterial transposases and retroviral integrases, cleavage by the RAG proteins requires a divalent metal ion but does not involve a covalent protein/DNA intermediate. In the transposase/integrase family, a triplet of acidic residues, commonly called a DDE motif, is often found to coordinate the metal ion used for catalysis. We show here that mutations in each of three acidic residues in RAG1 result in mutant derivatives that can bind the RSS but whose ability to catalyze either of the two chemical steps of V(D)J cleavage (nicking and hairpin formation) is severely impaired. Because both chemical steps are affected by the same mutations, a single active site appears responsible for both reactions. Two independent lines of evidence demonstrate that at least two of these acidic residues are directly involved in coordinating a divalent metal ion: The substitution of Cys for Asp allows rescue of some catalytic function, whereas an alanine substitution is no longer subject to iron-induced hydroxyl radical cleavage. Our results support a model in which the RAG1 protein contains the active site of the V(D)J recombinase and are interpreted in light of predictions about the structure of RAG1.

Benchmarking PSI-BLAST in genome annotation

The recognition of remote protein homologies is a major aspect of the structural and functional annotation of newly determined genomes. Here we benchmark the coverage and error rate of genome annotation using the widely used homology-searching program PSI-BLAST (position-specific iterated basic local alignment search tool). This study evaluates the one-to-many success rate for recognition, as often there are several homologues in the database and only one needs to be identified for annotating the sequence. In contrast, previous benchmarks considered one-to-one recognition in which a single query was required to find a particular target. The benchmark constructs a model genome from the full sequences of the structural classification of protein (SCOP) database and searches against a target library of remote homologous domains (<20 % identity). The structural benchmark provides a reliable list of correct and false homology assignments. PSI-BLAST successfully annotated 40 % of the domains in the model genome that had at least one homologue in the target library. This coverage is more than three times that if one-to-one recognition is evaluated (11 % coverage of domains). Although a structural benchmark was used, the results equally apply to just sequence homology searches. Accordingly, structural and sequence assignments were made to the sequences of Mycoplasma genitalium and Mycobacterium tuberculosis (see http://www.bmm.icnet. uk). The extent of missed assignments and of new superfamilies can be estimated for these genomes for both structural and functional annotations.

The extracellular regions of PSMA and the transferrin receptor contain an aminopeptidase domain: implications for drug design

The Aeromonas proteolytica aminopeptidase (AMP), Pseudomonas sp. (RS-16) carboxypeptidase G2 (CPG2), and Streptomyces griseus aminopeptidase (SGAP) are zinc dependent proteolytic enzymes with cocatalytic zinc ion centers and a conserved aminopeptidase fold. A BLAST search with the sequence of the solved AMP structure indicated that a similar domain could be found in prostate-specific membrane antigen (PSMA) and the transferrin receptor (TfR). When the PSMA or TfR sequence was input into the THREADER program, the top structural matches were SGAP and AMP confirming that these are structurally conserved domains. Optimal sequence alignment of PSMA and TfR using the known three-dimensional structures of AMP, CPG2, and SGAP shows that the critical amino acids involved in forming the catalytic pocket are conserved in PSMA but absent in the TfR. The specificity pocket in AMP is formed from four aromatic side chains and the equivalent region in CPG2/PSMA has a changed sequence pattern. Since CPG2 and PSMA are folate hydrolases, the changed specificity pocket leaves space to accommodate the large pteroate moiety of folic acid. In contrast, no enzyme function has been ascribed to the TfR.

Advances in structural genomics

New computational techniques have allowed protein folds to be assigned to all or parts of between a quarter (Caenorhabditis elegans) and a half (Mycoplasma genitalium) of the individual protein sequences in different genomes. These assignments give a new perspective on domain structures, gene duplications, protein families and protein folds in genome sequences.