To investigate the occurrence of glycosylphosphatidylinositol (GPI) lipid anchor modification in various taxonomic ranges, potential substrate proteins have been searched for in completely sequenced genomes. We applied the big-pi predictor for the recognition of propeptide cleavage and anchor attachment sites with a new, generalized analytical form of the extreme-value distribution for evaluating false-positive prediction rates. (i) We find that GPI modification is present among lower and higher Eukaryota (approximately 0.5% of all proteins) but it seems absent in all eubacterial and three archaeobacterial species studied. Four other archaean genomes appear to encode such a fraction of substrate proteins (in the range of eukaryots) that they cannot be explained as false-positive predictions. This result supports the possible existence of GPI anchor modification in an archaean subgroup. (ii) The frequency of GPI-modified proteins on various chromosomes of a given eukaryotic species is different. (iii) Lists of potentially GPI-modified proteins in complete genomes with their predicted cleavage sites are available at http://mendel.imp.univie.ac.at/gpi/gpi_genomes.html. (iv) Orthologues of known transamidase subunits have been found only for EUKARYA: Inconsistencies in domain structure among homologues some of which may indicate sequencing errors are described. We present a refined model of the transamidase complex.

Sequence weighting techniques are aimed at balancing redundant observed information from subsets of similar sequences in multiple alignments. Traditional approaches apply the same weight to all positions of a given sequence, hence equal efficiency of phylogenetic changes is assumed along the whole sequence. This restrictive assumption is not required for the new method PSIC (position-specific independent counts) described in this paper. The number of independent observations (counts) of an amino acid type at a given alignment position is calculated from the overall similarity of the sequences that share the amino acid type at this position with the help of statistical concepts. This approach allows the fast computation of position-specific sequence weights even for alignments containing hundreds of sequences. The PSIC approach has been applied to profile extraction and to the fold family assignment of protein sequences with known structures. Our method was shown to be very productive in finding distantly related sequences and more powerful than Hidden Markov Models or the profile methods in WiseTools and PSI-BLAST in many cases. The profile extraction routine is available on the WWW (http://www.bork.embl-heidelberg. de/PSIC or http://www.imb.ac.ru/PSIC).

Glycosylphosphatidylinositol (GPI) anchoring is a common post-translational modification of extracellular eukaryotic proteins. Attachment of the GPI moiety to the carboxyl terminus (omega-site) of the polypeptide occurs after proteolytic cleavage of a C-terminal propeptide. In this work, the sequence pattern for GPI-modification was analyzed in terms of physical amino acid properties based on a database analysis of annotated proprotein sequences. In addition to a refinement of previously described sequence signals, we report conserved sequence properties in the regions omega - 11...omega - 1 and omega + 4...omega + 5. We present statistical evidence for volume-compensating residue exchanges with respect to the positions omega - 1...omega + 2. Differences between protozoan and metazoan GPI-modification motifs consist mainly in variations of preferences to amino acid types at the positions near the omega-site and in the overall motif length. The variations of polypeptide substrates are exploited to suggest a model of the polypeptide binding site of the putative transamidase, the enzyme catalyzing the GPI-modification. The volume of the active site cleft accommodating the four residues omega - 1...omega + 2 appears to be approximately 540 A3.

Specific residue interactions as revealed from a few and readily available experiments can be quite important in shaping a protein's tertiary topology by complementing basic and general folding principles. This experimental information is employed in structure prediction (mainchain topology) based on sequence knowledge and the genetic algorithm with its ability to optimize simultaneously many parameters. Examples investigated include the distribution of cysteinyl S-S bonds, protein side-chain ligands to iron-sulfur cages, cofactor-ligands, crosslinks amongst side-chains, and conserved hydrophobic and catalytic residues. Such interactions yield an improvement in the predicted topology (0.4-6.6 A root mean square deviation in the positions of the backbone C alpha-atoms relative to those observed) compared with those resulting from simulations relying only on basic protein folding principles. For several examples the resultant topology depended critically on knowledge of the few and specific interactions such that the relationship between predicted and observed C alpha-positions was near random without their use. The combined methodology (experimental data and the genetic algorithm) should prove helpful in settings where experiment and theory can cooperate in successive steps to elucidate an unknown structure.

The hydrophobic part of the solvent-accessible surface of a typical monomeric globular protein consists of a single, large interconnected region formed from faces of apolar atoms and constituting approximately 60% of the solvent-accessible surface area. Therefore, the direct delineation of the hydrophobic surface patches on an atom-wise basis is impossible. Experimental data indicate that, in a two-state hydration model, a protein can be considered to be unified with its first hydration shell in its interaction with bulk water. We show that, if the surface area occupied by water molecules bound at polar protein atoms as generated by AUTOSOL is removed, only about two-thirds of the hydrophobic part of the protein surface remains accessible to bulk solvent. Moreover, the organization of the hydrophobic part of the solvent-accessible surface experiences a drastic change, such that the single interconnected hydrophobic region disintegrates into many smaller patches, i.e. the physical definition of a hydrophobic surface region as unoccupied by first hydration shell water molecules can distinguish between hydrophobic surface clusters and small interconnecting channels. It is these remaining hydrophobic surface pieces that probably play an important role in intra- and intermolecular recognition processes such as ligand binding, protein folding and protein-protein association in solution conditions. These observations have led to the development of an accurate and quick analytical technique for the automatic determination of hydrophobic surface patches of proteins. This technique is not aggravated by the limiting assumptions of the methods for generating explicit water hydration positions. Formation of the hydrophobic surface regions owing to the structure of the first hydration shell can be computationally simulated by a small radial increment in solvent-accessible polar atoms, followed by calculation of the remaining exposed hydrophobic patches. We demonstrate that a radial increase of 0.35-0.50 A resembles the effect of tightly bound water on the organization of the hydrophobic part of the solvent-accessible surface.