Structures for protein domains have increased rapidly in recent years owing to advances in structural biology and structural genomics projects. New structures are often similar to those solved previously, and such similarities can give insights into function by linking poorly understood families to those that are better characterized. They also allow the possibility of combing information to find still more proteins adopting a similar structure and sometimes a similar function, and to reprioritize families in structural genomics pipelines. We explore this possibility here by preparing merged profiles for pairs of structurally similar, but not necessarily sequence-similar, domains within the SMART and Pfam database by way of the Structural Classification of Proteins (SCOP). We show that such profiles are often able to successfully identify further members of the same superfamily and thus can be used to increase the sensitivity of database searching methods like HMMer and PSI-BLAST. We perform detailed benchmarks using the SMART and Pfam databases with four complete genomes frequently used as annotation benchmarks. We quantify the associated increase in structural information in Swissprot and discuss examples illustrating the applicability of this approach to understand functional and evolutionary relationships between protein families.

Sequence similarity is the most common measure currently used to infer homology between proteins. Typically, homologous protein domains show sequence similarity over their entire lengths. Here we identify Asp box motifs, initially found as repeats in sialidases and neuraminidases, in new structural and sequence contexts. These motifs represent significantly similar sequences, localized to beta hairpins within proteins that are otherwise different in sequence and three-dimensional structure. By performing a combined sequence- and structure-based analysis we detect Asp boxes in more than nine protein families, including bacterial ribonucleases, sulfite oxidases, reelin, netrins, some lipoprotein receptors, and a variety of glycosyl hydrolases. Although the function common to each of these proteins, if any, remains unclear, we discuss possible functions of Asp boxes on the basis of previously determined experimental results and discuss different evolutionary scenarios for the origin of Asp-box containing proteins.

The sterile alpha motif (SAM) is a protein interaction domain of around 70 amino acids present predominantly in the N- and C-termini of more than 60 diverse proteins that participate in signal transduction and transcriptional repression. SAM domains have been shown to homo- and hetero-oligomerize and to mediate specific protein-protein interactions. A highly conserved subclass of SAM domains is present at the intracellular C-terminus of more than 40 Eph receptor tyrosine kinases that are involved in the control of axonal pathfinding upon ephrin-induced oligomerization and activation in the event of cell-cell contacts. These SAM domains appear to participate in downstream signaling events via interactions with cytosolic proteins. We determined the solution structure of the EphB2 receptor SAM domain and studied its association behavior. The structure consists of five helices forming a compact structure without binding pockets or exposed conserved aromatic residues. Concentration-dependent chemical shift changes of NMR signals reveal two distinct well-separated areas on the domains' surface sensitive to the formation of homotypic oligomers in solution. These findings are supported by analytical ultracentrifugation studies. The conserved Tyr932, which was reported to be essential for the interaction with SH2 domains after phosphorylation, is buried in the hydrophobic core of the structure. The weak capability of the isolated EphB2 receptor SAM domain to form oligomers is supposed to be relevant in vivo when the driving force of ligand binding induces receptor oligomerization. A formation of SAM tetramers is thought to provide an appropriate contact area for the binding of a low-molecular-weight phosphotyrosine phosphatase and to initiate further downstream responses.

More than 60 previously undetected SAM domain-containing proteins have been identified using profile searching methods. Among these are over 40 EPH-related receptor tyrosine kinases (RPTK), Drosophila bicaudal-C, a p53 from Loligo forbesi, and diacyglycerol-kinase isoform delta. This extended dataset suggests that SAM is an evolutionary conserved protein binding domain that is involved in the regulation of numerous developmental processes among diverse eukaryotes. A conserved tyrosine in the SAM sequences of the EPH related RPTKs is likely to mediate cell-cell initiated signal transduction via the binding of SH2 containing proteins to phosphotyrosine.

A thorough sequence analysis of the various members of the eukaryotic protein serine/threonine phosphatase 2C (PP2C) family revealed the conservation of 11 motifs. These motifs could be identified in numerous other sequences, including fungal adenylate cyclases that are predicted to contain a functionally active PP2C domain, and a family of prokaryotic serine/threonine phosphatases including SpoIIE. Phylogenetic analysis of all the proteins indicates a widespread sequence family for which a considerable number of isoenzymes can be inferred.

For the first time, a direct approach for the derivation of an atomic solvation parameter from macromolecular structural data alone is presented. The specific free energy of solvation for hydrophobic surface regions of proteins is delineated from the area distribution of hydrophobic surface patches. The resulting value is 18 cal/(mol.A2), with a statistical uncertainty of +/-2 cal/mol.A2) at the 5% significance level. It compares favorably with the parameters for carbon obtained by other authors who use the the crystal geometry of succinic acid or energies of transfer from hydrophobic solvent to water for small organic compounds. Thus, the transferability of atomic solvation parameters for hydrophobic atoms to macromolecules has been directly demonstrated. A careful statistical analysis demonstrates that surface energy parameters derived from thermodynamic data of protein mutation experiments are clearly less confident.

We present a method based on hierarchical self-organizing maps (SOMs) for recognizing patterns in protein sequences. The method is fully automatic, does not require prealigned sequences, is insensitive to redundancy in the training set, and works surprisingly well even with small learning sets. Because it uses unsupervised neural networks, it is able to extract patterns that are not present in all of the unaligned sequences of the learning set. The identification of these patterns in sequence databases is sensitive and efficient. The procedure comprises three main training stages. In the first stage, one SOM is trained to extract common features from the set of unaligned learning sequences. A feature is a number of ungapped sequence segments (usually 4-16 residues long) that are similar to segments in most of the sequences of the learning set according to an initial similarity matrix. In the second training stage, the recognition of each individual feature is refined by selecting an optimal weighting matrix out of a variety of existing amino acid similarity matrices. In a third stage of the SOM procedure, the position of the features in the individual sequences is learned. This allows for variants with feature repeats and feature shuffling. The procedure has been successfully applied to a number of notoriously difficult cases with distinct recognition problems: helix-turn-helix motifs in DNA-binding proteins, the CUB domain of developmentally regulated proteins, and the superfamily of ribokinases. A comparison with the established database search procedure PROFILE (and with several others) led to the conclusion that the new automatic method performs satisfactorily.

Study of the most conserved region in many beta/alpha-barrels, the phosphate-binding site, revealed a sequence motif in a few beta/alpha-barrels with known tertiary structure, namely glycolate oxidase (GOX), cytochrome b2 (Cyb2), tryptophan synthase alpha subunit (TrpA), and the indoleglycerolphosphate synthase (TrpC). Database searches identified this motif in numerous other enzyme families: (1) IMP dehydrogenase (IMPDH) and GMP reductase (GuaC); (2) phosphoribosylformimino-5-aminoimidazol carboxamide ribotide isomerase (HisA) and the cyclase-producing D-erythro-imidazole-glycerolphosphate (HisF) of the histidine biosynthetic pathway; (3) dihydroorotate dehydrogenase (PyrD); (4) glutamate synthase (GltB); (5) ThiE and ThiG involved in the biosynthesis of thiamine as well as related proteins; (6) an uncharacterized open reading frame from Erwinia herbicola; and (7) a glycerol uptake operon antiterminator regulatory protein (GlpP). Secondary structure predictions of the different families mentioned above revealed an alternating order of beta-strands and alpha-helices in agreement with a beta/alpha-barrel-like topology. The putative phosphate-binding site is always found near the C-terminus of the enzymes, which are all at least about 200 amino acids long. This is compatible with its assumed location between strand 7 and helix 8. The identification of a significant motif in functionally diverse enzymes suggests a divergent evolution of at least a considerable fraction of beta/alpha-barrels. In addition to the known accumulation of beta/alpha-barrels in the tryptophan biosynthetic pathway, we observe clusters of these enzymes in histidine biosynthesis, purine metabolism, and apparently also in thiamine biosynthesis. The substrates are mostly heterocyclic compounds.(ABSTRACT TRUNCATED AT 250 WORDS)

Using a variety of homology search methods and multiple alignments, a new extracellular module was identified in (1) agrin, (2) enterokinase, (3) a 63-kDa sea urchin sperm protein, (4) perlecan, (5) the breast cancer marker MUCI (episialin), (6) the cell surface antigen 114/A10, and (7/8) two functionally uncharacterized, probably extracellular, Caenorhabditis elegans proteins. Despite the functional diversity of these adhesive proteins, a common denominator seems to be their existence in heavily glycosylated environments. In addition, the better characterized proteins mentioned above contain all O-glycosidic-linked carbohydrates such as heparan sulfate that contribute considerably to their molecular masses. The common module might regulate or assist binding to neighboring carbohydrate moieties.

The analysis of the 269 open reading frames of yeast chromosome VIII by computational methods has yielded 24 new significant sequence similarities to proteins of known function. The resulting predicted functions include three particularly interesting cases of translation-associated proteins: peptidyl-tRNA hydrolase, a ribosome recycling factor homologue, and a protein similar to cytochrome b translational activator CBS2. The methodological limits of the meaningful transfer of functional information between distant homologues are discussed.

Using computer methods for database search and multiple alignment, statistically significant sequence similarities were identified between several nitrilases with distinct substrate specificity, cyanide hydratases, aliphatic amidases, beta-alanine synthase, and a few other proteins with unknown molecular function. All these proteins appear to be involved in the reduction of organic nitrogen compounds and ammonia production. Sequence conservation over the entire length, as well as the similarity in the reactions catalyzed by the known enzymes in this family, points to a common catalytic mechanism. The new family of enzymes is characterized by several conserved motifs, one of which contains an invariant cysteine that is part of the catalytic site in nitrilases. Another highly conserved motif includes an invariant glutamic acid that might also be involved in catalysis.

Using computer methods for multiple alignment, sequence motif search, and tertiary structure modeling, we show that eukaryotic translation elongation factor 1 gamma (EF1 gamma) contains an N-terminal domain related to class theta glutathione S-transferases (GST). GST-like proteins related to class theta comprise a large group including, in addition to typical GSTs and EF1 gamma, stress-induced proteins from bacteria and plants, bacterial reductive dehalogenases and beta-etherases, and several uncharacterized proteins. These proteins share 2 conserved sequence motifs with GSTs of other classes (alpha, mu, and pi). Tertiary structure modeling showed that in spite of the relatively low sequence similarity, the GST-related domain of EF1 gamma is likely to form a fold very similar to that in the known structures of class alpha, mu, and pi GSTs. One of the conserved motifs is implicated in glutathione binding, whereas the other motif probably is involved in maintaining the proper conformation of the GST domain. We predict that the GST-like domain in EF1 gamma is enzymatically active and that to exhibit GST activity, EF1 gamma has to form homodimers. The GST activity may be involved in the regulation of the assembly of multisubunit complexes containing EF1 and aminoacyl-tRNA synthetases by shifting the balance between glutathione, disulfide glutathione, thiol groups of cysteines, and protein disulfide bonds. The GST domain is a widespread, conserved enzymatic module that may be covalently or noncovalently complexed with other proteins. Regulation of protein assembly and folding may be 1 of the functions of GST.

Kinases that catalyze phosphorylation of sugars, called here sugar kinases, can be divided into at least three distinct nonhomologous families. The first is the hexokinase family, which contains many prokaryotic and eukaryotic sugar kinases with diverse specificities, including a new member, rhamnokinase from Salmonella typhimurium. The three-dimensional structure of hexokinase is known and can be used to build models of functionally important regions of other kinases in this family. The second is the ribokinase family, of unknown three-dimensional structure, and comprises pro- and eukaryotic ribokinases, bacterial fructokinases, the minor 6-phosphofructokinase 2 from Escherichia coli, 6-phosphotagatokinase, 1-phosphofructokinase, and, possibly, inosine-guanosine kinase. The third family, also of unknown three-dimensional structure, contains several bacterial and yeast galactokinases and eukaryotic mevalonate and phosphomevalonate kinases and may have a substrate binding region in common with homoserine kinases. Each of the three families of sugar kinases appears to have a distinct three-dimensional fold, since conserved sequence patterns are strikingly different for the three families. Yet each catalyzes chemically equivalent reactions on similar or identical substrates. The enzymatic function of sugar phosphorylation appears to have evolved independently on the three distinct structural frameworks, by convergent evolution. In addition, evolutionary trees reveal that (1) fructokinase specificity has evolved independently in both the hexokinase and ribokinase families and (2) glucose specificity has evolved independently in different branches of the hexokinase family. These are examples of independent Darwinian adaptation of a structure to the same substrate at different evolutionary times. The flexible combination of active sites and three-dimensional folds observed in nature can be exploited by protein engineers in designing and optimizing enzymatic function.

With the completion of the first phase of the European yeast genome sequencing project, the complete DNA sequence of chromosome III of Saccharomyces cerevisiae has become available (Oliver, S. G., et al., 1992, Nature 357, 38-46). We have tested the predictive power of computer sequence analysis of the 176 probable protein products of this chromosome, after exclusion of six problem cases. When the results of database similarity searches are pooled with prior knowledge, a likely function can be assigned to 42% of the proteins, and a predicted three-dimensional structure to a third of these (14% of the total). The function of the remaining 58% remains to be determined. Of these, about one-third have one or more probable transmembrane segments. Among the most interesting proteins with predicted functions are a new member of the type X polymerase family, a transcription factor with an N-terminal DNA-binding domain related to GAL4, a "fork head" DNA-binding domain previously known only in Drosophila and in mammals, and a putative methyltransferase. Our analysis increased the number of known significant sequence similarities on chromosome III by 13, to now 67. Although the near 40% success rate of identifying unknown protein function by sequence analysis is surprisingly high, the information gap between known protein sequences and unknown function is expected to widen and become a major bottleneck of genome projects in the near future. Based on the experience gained in this test study, we suggest that the development of an automated computer workbench for protein sequence analysis must be an important item in genome projects.