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.

We have analyzed the base-pairing probability distributions of 16 S and 16 S-like, and 23 S and 23 S-like ribosomal RNAs of Archaea, Bacteria, chloroplasts, mitochondria and Eukarya, as predicted by the partition function approach for RNA folding introduced by McCaskill. A quantitative analysis of the reliability of RNA folding is done by comparing the base-pairing probability distributions with the structures predicted by comparative sequence analysis (comparative structures). We distinguish two factors that show a relationship to the reliability of RNA minimum free energy structure. The first factor is the dominance of one particular base-pair or the absence of base-pairing for a given base within the base-pairing probability distribution (BPPD). We characterize the BPPD per base, including the probability of not base-pairing, by its Shannon entropy (S). The S value indicates the uncertainty about the base-pairing of a base: low S values result from BPPDs that are strongly dominated by a single base-pair or by the absence of base-pairing. We show that bases with low S values have a relatively high probability that their minimum free energy (MFE) structure corresponds to the comparative structure. The BPPDs of prokaryotes that live at high temperatures (thermophilic Archaea and Bacteria) have, calculated at 37 degrees C, lower S values than the BPPDs of prokaryotes that live at lower temperatures (mesophilic and psychrophilic Archaea and Bacteria). This reflects an adaptation of the ribosomal RNAs to the environmental temperature. A second factor that is important to consider with regard to the reliability of MFE structure folding is a variable degree of applicability of the thermodynamic model of RNA folding for different groups of RNAs. Here we show that among the bases that show low S values, the Archaea and Bacteria have similar, high probabilities (0.96 and 0.94 in 16 S and 0.93 and 0.91 in 23 S, respectively) that the MFE structure corresponds to the comparative structure. These probabilities are lower in the chloroplasts (16 S 0.91, 23 S 0.79), mitochondria (16 S-like 0.89, 23 S-like 0.69) and Eukarya (18 S 0.81, 28 S 0.86).

Grid-free protein folding simulations based on sequence and secondary structure knowledge (using mostly experimentally determined secondary structure information but also analysing results from secondary structure predictions) were investigated using the genetic algorithm, a backbone representation, and standard dihedral angular conformations. Optimal structures are selected according to basic protein building principles. Having previously applied this approach to proteins with helical topology, we have now developed additional criteria and weights for beta-strand-containing proteins, validated them on four small beta-strand-rich proteins with different topologies, and tested the general performance of the method on many further examples from known protein structures with mixed secondary structural type and less than 100 amino acid residues. Topology predictions close to the observed experimental structures were obtained in four test cases together with fitness values that correlated with the similarity of the predicted topology to the observed structures. Root-mean-square deviation values of C alpha atoms in the superposed predicted and observed structures, the latter of which had different topologies, were between 4.5 and 5.5 A(2.9 to 5.1 A without loops). Including 15 further protein examples with unique folds, root-mean-square deviation values ranged between 1.8 and 6.9 A with loop regions and averaged 5.3 A and 4.3 A, including and excluding loop regions, respectively.

The geometry of helix-helix packing in globular proteins is comprehensively analysed within the model of the superposition of two helix lattices which result from unrolling the helix cylinders onto a plane containing points representing each residue. The requirements for the helix geometry (the radius R, the twist angle omega and the rise per residue delta) under perfect match of the lattices are studied through a consistent mathematical model that allows consideration of all possible associations of all helix types (alpha-, pi- and 3(10)). The corresponding equations have three well-separated solutions for the interhelical packing angle, omega, as a function of the helix geometric parameters allowing optimal packing. The resulting functional relations also show unexpected behaviour. For a typically observed alpha-helix (omega = 99.1 degrees, delta = 1.45 A), the three optimal packing angles are omega a,b,c = -37.1 degrees, -97.4 degrees and +22.0 degrees with a periodicity of 180 degrees and respective helix radii Ra,b,c = 3.0 A, 3.5 A and 4.3 A. However, the resulting radii are very sensitive to variations in the twist angle omega. At omega triple = 96.9 degrees, all three solutions yield identical radii at delta = 1.45 A where Rtriple = 3.46 A. This radius is close to that of a poly(Ala) helix, indicating a great packing flexibility when alanine is involved in the packing core, and omega triple is close to the mean observed twist angle. In contrast, the variety of possible theoretical solutions is limited for the other two helix types. Besides the perfect matches, novel suboptimal "knobs into holes" hydrophobic packing patterns as a function of the helix radius are described. Alternative "knobs onto knobs" and mixed models can be applied in cases where salt bridges, hydrogen bonds, disulphide bonds and tight hydrophobic head-to-head contacts are involved in helix-helix associations. An analysis of the experimentally observed packings in proteins confirmed the conclusions of the theoretical model. Nonetheless, the observed alpha-helix packings showed deviations from the 180 degrees periodicity expected from the model. An investigation of the actual three-dimensional geometry of helix-helix packing revealed an explanation for the observed discrepancies where a decisive role was assigned to the defined orientation of the C alpha-C beta vectors of the side-chains. As predicted form the model, helices with different radii (differently sized side-chains in the packing core) were observed to utilize different packing cells (packing patterns). In agreement with the coincidence between Rtriple and the radius of a poly(Ala) helix, Ala was observed to show greatest propensity to build the packing core. The application of the helix lattice superposition model suggests that the packing of amino acid residues is best described by a "knobs into holes" scheme rather than "ridges into grooves". The various specific packing modes made salient by the model should be useful in protein engineering and design.

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.