ARM and HEAT motifs are tandemly repeated sequences of approximately 50 amino acid residues that occur in a wide variety of eukaryotic proteins. An exhaustive search of sequence databases detected new family members and revealed that at least 1 in 500 eukaryotic protein sequences contain such repeats. It also rendered the similarity between ARM and HEAT repeats, believed to be evolutionarily related, readily apparent. All the proteins identified in the database searches could be clustered by sequence similarity into four groups: canonical ARM-repeat proteins and three groups of the more divergent HEAT-repeat proteins. This allowed us to build improved sequence profiles for the automatic detection of repeat motifs. Inspection of these profiles indicated that the individual repeat motifs of all four classes share a common set of seven highly conserved hydrophobic residues, which in proteins of known three-dimensional structure are buried within or between repeats. However, the motifs differ at several specific residue positions, suggesting important structural or functional differences among the classes. Our results illustrate that ARM and HEAT-repeat proteins, while having a common phylogenetic origin, have since diverged significantly. We discuss evolutionary scenarios that could account for the great diversity of repeats observed.

Short protein repeats, frequently with a length between 20 and 40 residues, represent a significant fraction of known proteins. Many repeats appear to possess high amino acid substitution rates and thus recognition of repeat homologues is highly problematic. Even if the presence of a certain repeat family is known, the exact locations and the number of repetitive units often cannot be determined using current methods. We have devised an iterative algorithm based on optimal and sub-optimal score distributions from profile analysis that estimates the significance of all repeats that are detected in a single sequence. This procedure allows the identification of homologues at alignment scores lower than the highest optimal alignment score for non-homologous sequences. The method has been used to investigate the occurrence of eleven families of repeats in Saccharomyces cerevisiae, Caenorhabditis elegans and Homo sapiens accounting for 1055, 2205 and 2320 repeats, respectively. For these examples, the method is both more sensitive and more selective than conventional homology search procedures. The method allowed the detection in the SwissProt database of more than 2000 previously unrecognised repeats belonging to the 11 families. In addition, the method was used to merge several repeat families that previously were supposed to be distinct, indicating common phylogenetic origins for these families.

We provide statistically reliable sequence evidence indicating that at least 12 of 23 SCOP (betaalpha)(8) (TIM) barrel superfamilies share a common origin. This includes all but one of the known and predicted TIM barrels found in central metabolism. The statistical evidence is complemented by an examination of the details of protein structure, with certain structural locations favouring catalytic residues even though the nature of their molecular function may change. The combined analysis of sequence, structure and function also enables us to propose a phylogeny of TIM barrels. Based on these data, we are able to examine differing theories of pathway and enzyme evolution, by mapping known TIM barrel folds to the pathways of central metabolism. The results favour widespread recruitment of enzymes between pathways, rather than a "backwards evolution" model, and support the idea that modern proteins may have arisen from common ancestors that bound key metabolites.

Glycosylphosphatidylinositol (GPI) lipid anchoring is a common posttranslational modification known mainly from extracellular eukaryotic proteins. Attachment of the GPI moiety to the carboxyl terminus (omega-site) of the polypeptide follows after proteolytic cleavage of a C-terminal propeptide. For the first time, a new prediction technique locating potential GPI-modification sites in precursor sequences has been applied for large-scale protein sequence database searches. The composite prediction function (with separate parametrisation for metazoan and protozoan proteins) consists of terms evaluating both amino acid type preferences at sequence positions near a supposed omega-site as well as the concordance with general physical properties encoded in multi-residue correlation within the motif sequence. The latter terms are especially successful in rejecting non-appropriate sequences from consideration. The algorithm has been validated with a self-consistency and two jack-knife tests for the learning set of fully annotated sequences from the SWISS-PROT database as well as with a newly created database "big-Pi" (more than 300 GPI-motif mutations extracted from original literature sources). The accuracy of predicting the effect of mutations in the GPI sequence motif was above 83 %. Lists of potential precursor proteins which are non-annotated in SWISS-PROT and SPTrEMBL are presented on the WWW-page http://www.embl-heidelberg.de/beisenha/gpi/gpi_p rediction. html The algorithm has been implemented in the prototype software "big-Pi predictor" which may find application as a genome annotation and target selection tool.

The multidomain Mac-2 binding protein (M2BP) is present in serum and in the extracellular matrix in the form of linear and ring-shaped oligomers, which interact with galectin-3, fibronectin, collagens, integrins and other large glycoproteins. Domain 1 of M2BP (M2BP-1) shows homology with the cysteine-rich SRCR domain of scavanger receptor. Domains 2 and 3 are related to the dimerization domains BTB/POZ and IVR of the Drosophila kelch protein. Recombinant M2BP, its N-terminal domain M2BP-1 and a fragment consisting of putative domains 2, 3 and 4 (M2BP-2,3,4) were investigated by scanning transmission electron microscopy, transmission electron microscopy, analytical ultracentrifugation and binding assays. The ring oligomers formed by the intact protein are comprised of approximately 14 nm long segments composed of two 92 kDa M2BP monomers. Although the rings vary in size, decamers predominate. The various linear oligomers also observed are probably ring precursors, dimers predominate. M2BP-1 exhibits a native fold, does not oligomerize and is inactive in cell attachment. M2BP-2,3,4 aggregates to heterogeneous, protein filled ring-like structures as shown by metal shadowed preparations. These aggregates retain the cell-adhesive potential indicating native folding. It is hypothesized that the rings provide an interaction pattern for multivalent interactions of M2BP with target molecules or complexes of ligands.

Phyletic distributions of eukaryotic signalling domains were studied using recently developed sensitive methods for protein sequence analysis, with an emphasis on the detection and accurate enumeration of homologues in bacteria and archaea. A major difference was found between the distributions of enzyme families that are typically found in all three divisions of cellular life and non-enzymatic domain families that are usually eukaryote-specific. Previously undetected bacterial homologues were identified for# plant pathogenesis-related proteins, Pad1, von Willebrand factor type A, src homology 3 and YWTD repeat-containing domains. Comparisons of the domain distributions in eukaryotes and prokaryotes enabled distinctions to be made between the domains originating prior to the last common ancestor of all known life forms and those apparently originating as consequences of horizontal gene transfer events. A number of transfers of signalling domains from eukaryotes to bacteria were confidently identified, in contrast to only a single case of apparent transfer from eukaryotes to archaea.

Using a new method for construction and database searches of sequence consensus strings, we have identified a new superfamily of protein modules comprising laminin G, thrombospondin N and the pentraxin families. The conserved patterns correspond mainly to hydrophobic core residues located in central beta strands of the known three-dimensional structures of two pentraxins, the human C-reactive protein and the serum amyloid P-component. Thus, we predict a similar jellyroll fold for all members of this superfamily. In addition, the conservation of two exposed aspartate residues in the majority of superfamily members suggests hitherto unrecognised functional sites.

Predicting function from sequence using computational tools is a highly complicated procedure that is generally done for each gene individually. This review focuses on the added value that is provided by completely sequenced genomes in function prediction. Various levels of sequence annotation and function prediction are discussed, ranging from genomic sequence to that of complex cellular processes. Protein function is currently best described in the context of molecular interactions. In the near future it will be possible to predict protein function in the context of higher order processes such as the regulation of gene expression, metabolic pathways and signalling cascades. The analysis of such higher levels of function description uses, besides the information from completely sequenced genomes, also the additional information from proteomics and expression data. The final goal will be to elucidate the mapping between genotype and phenotype.

Homology search techniques based on the iterative PSI-BLAST method in combination with various filters for low sequence complexity are applied to assign folds to all Mycoplasma genitalium proteins. The resulting procedure (implemented as a web server) is able to predict at least one domain in 37% of these proteins automatically, with an estimated accuracy higher than 98%. Taking structural features such as coiled coil or transmembrane regions aside, folds can be assigned to more than half of the globular proteins in a bacterium just by iterative sequence comparison.

A systematic screening of sequence databases with a motif hitherto found only in animal and poxvirus proteins has revealed a trail leading back to prokaryotes. Fortuitously, an X-ray structure is available for one of the identified sequences and shows the fundamental fold to be a set of circularly arranged beta sheets. This structure may be very widely distributed throughout the biological world in sialidases and some other enzymes. In bacteria, a mobile noncatalytic domain is often associated with these same enzymes.

Since the first crystal structure of an immunoglobulin revealed a modular architecture, the characteristic beta-sheet fold of the immunoglobulin domain has been found in many other proteins of diverse biological function. Here, a systematic comparison of 23 Ig domain structures with less than 25% pairwise residue identity was performed using automatic structural alignment and analysis of beta-sheet and loop topology. Sequence consensus patterns were identified for nine distinct families with at most marginal similarity to each other. The analysis reveals a common structural core of only four beta-strands (b, c, e and f), embedded in an antiparallel curled beta-sheet sandwich with a total of three to five additional strands (a, c', c'', d, g) and a characteristic intersheet angle. The variation in the position of the edge strands (a, c', c'', d and g) relative to the common core defines four different topological subtypes that correlate with the length of the intervening sequence between strands c and e, the most variable region in sequence. The switch of strand c' from one sheet to the other in seven-stranded domains appears to result from short c-e segments, rather than being a major structural discriminator. The high degree of structural flexibility outside the common core and the extreme variability of side-chain packing inside the core do not support a protein folding pathway common to all members of the structural class. Mutation rates of immunoglobulin-like domains in different proteins vary considerably. Disulfide bridges, thought to contribute to structural stability, are not necessarily invariant in number and location within a subclass.

Sequence analysis has revealed the presence of 31 copies of an extracellular domain, here called CUB, in 16 functionally diverse proteins such as the dorso-ventral patterning protein tolloid, bone morphogenetic protein 1, a family of spermadhesins, complement subcomponents Cls/Clr and the neuronal recognition molecule A5. Most of them are known to be involved in developmental processes. Our analysis of this new family includes the identification of seven previously undescribed members, the characterization of conserved features and a topology prediction of this approximately 110 residue spanning domain, which suggests an antiparallel beta-barrel similar to those in immunoglobulins.

Cloning and sequence analysis of a new open reading frame from Bacillus cereus reveals the relationship to a recently identified family of putative eukaryotic transcription activators similar to the yeast SNF2 gene product. As a result of comparative analysis of sequence features conserved in all members of this family, a gene from a chilo iridescent virus, as well as a putative helicase from Escherichia coli (hepA), can also be grouped into this family. The unexpected presence of prokaryotic and viral sequences in the previously purely eukaryotic SNF2 family suggests a defined subgroup of DNA helicases present in all species, with specific function in transcription activation.