[Prenylation: from bacteria to eukaryotes].

For their protection from host cell immune defense, intracellular eukaryotic parasites developed a variety of mechanisms, including secretion systems III and IV which inject bacterial effectors directly into eukaryotic cells. These effectors may be posttranslational modified by host cell machinery and may function inside the host cell. Recently, to the list of possible posttranslational modifications of bacterial proteins the prenylation was added. In this work we describe current state of the knowledge about the prenylation of eukaryotic and prokaryotic proteins and its inhibitors. The bioinformatics analyses suggest possibility of prenylation for a number of Francisella genus proteins.

Can molecular mechanisms of biological processes be extracted from expression profiles? Case study: endothelial contribution to tumor-induced angiogenesis.

Whereas the genome contains all potential developmental programs, expression profiles permit the determination of genes that are actively transcribed under defined physiological conditions. In this article, the idea of extracting biological mechanisms from expression data is tested. Molecular processes of the endothelial contribution to angiogenesis are derived from recently published expression profiles. The analysis reveals the sensitivity limits of experimental detection of transcriptional changes and how sequence-analytic techniques can help to identify the function of genes in question. We conclude that the transcripts (http://mendel.imp.univie.ac.at/SEQUENCES/TEMS/) found to be up-regulated in angiogenesis are involved in extracellular matrix remodeling, cellular migration, adhesion, cell-cell communication rather than in angiogenesis initiation or integrative control. Comparison with tissue-specific patterns of EST occurrence shows that, indeed, the presumptive tumor-specific endothelial markers are more generally expressed by cell types involved in migration and matrix remodeling processes. This exemplary study demonstrates how bioinformatics approaches can be helpful in deriving mechanistic information from diverse sources of experimental data.

Taxonomy workbench.

UNLABELLED: At advanced stages of working with user-defined protein and gene sequence collections, it is frequently necessary to link these data to the taxonomic tree and to extract subsets in accordance with taxonomic considerations. Since no general automatic tools had been available, this was a tedious manual effort. Our taxonomy workbench allows processing of sequence sets, mapping of these sets onto the taxonomic tree, collection of taxonomic subsets from them and printing of the whole tree or some part of it. As a side effect, the system enables queries to and navigation within the taxonomy database. AVAILABILITY: An implementation of the taxonomy workbench is accessible for public use as a www-service at http://mendel.imp.univie.ac.at/taxonomy/. Software components for the command-line and for the www-version are available on request. CONTACT: Georg.Schneider@nt.imp.univie.ac.at; Frank.Eisenhaber@nt.imp.univie.ac.at SUPPLEMENTARY INFORMATION: Documentation for the taxonomy workbench can be accessed at http://mendel.imp.univie.ac.at/taxonomy/help.html.

Ubiquitylation in plants: a post-genomic look at a post-translational modification.

In this article, we summarize Arabidopsis genes encoding ubiquitin, ubiquitin-activating enzyme (E1), ubiquitin-conjugating enzymes (E2s) and an additional selected set of proteins related to ubiquitylation. We emphasize comparisons to components from Saccharomyces cerevisiae, with occasional reference to animals. Among the E1 and E2s, Arabidopsis usually has two to four probable orthologs to one yeast gene. Also, Arabidopsis has genes with no likely ortholog in yeast, although they often have potential orthologs in animals. The large number of components with known function in ubiquitylation indicates that this process plays a complex role in cellular physiology.

Molecular cloning and characterization of EndoGlyx-1, an EMILIN-like multisubunit glycoprotein of vascular endothelium.

EndoGlyx-1, the antigen identified with the monoclonal antibody H572, is a pan-endothelial human cell surface glycoprotein complex composed of four different disulfide-bonded protein species with an apparent molecular mass of approximately 500 kDa. Here, we report the purification and peptide analysis of two EndoGlyx-1 subunits, p125 and p140, and the identification of a common, full-length cDNA with an open reading frame of 2847 base pairs. The EndoGlyx-1 cDNA encodes a protein of 949 amino acids with a predicted molecular mass of 105 kDa, found as an entry for an unnamed protein with unknown function in public data bases. A short sequence tag matching the cDNA of this gene was independently discovered by serial analysis of gene expression profiling as a pan-endothelial marker, PEM87. Bioinformatic evaluation classifies EndoGlyx-1 as an EMILIN-like protein composed of a signal sequence, an N-terminal EMI domain, and a C-terminal C1q-like domain, separated from each other by a central coiled-coil-rich region. Biochemical and carbohydrate analysis revealed that p125, p140, and the two additional EndoGlyx-1 subunits, p110 and p200, are exposed on the cell surface. The three smaller subunits show a similar pattern of N-linked and O-linked carbohydrates, as shown by enzyme digestion. Because the two globular domains of EndoGlyx-1 p125/p140 show structural features shared by EMILIN-1 and Multimerin, two oligomerizing glycoproteins implicated in cell-matrix adhesion and hemostasis, it will be of interest to explore similar functions for EndoGlyx-1 in human vascular endothelium.

The Brix domain protein family -- a key to the ribosomal biogenesis pathway?

Six (one archaean and five eukaryotic) protein families have similar domain architecture that includes a central globular Brix domain, and optional N- and obligatory C-terminal segments, both with charged low-complexity regions. Biological data for some proteins in this superfamily suggest a role in ribosome biogenesis and rRNA binding.

Post-translational GPI lipid anchor modification of proteins in kingdoms of life: analysis of protein sequence data from complete genomes.

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.

Molecular cloning and characterization of endosialin, a C-type lectin-like cell surface receptor of tumor endothelium.

Endosialin, the antigen identified with monoclonal antibody FB5, is a highly restricted 165-kDa cell surface glycoprotein expressed by tumor blood vessel endothelium in a broad range of human cancers but not detected in blood vessels or other cell types in many normal tissues. Functional analysis of endosialin has been hampered by a lack of information about its molecular structure. In this study, we describe the purification and partial amino acid sequencing of endosialin, leading to the cloning of a full-length cDNA with an open reading frame of 2274 base pairs. The endosialin cDNA encodes a type I membrane protein of 757 amino acids with a predicted molecular mass of 80.9 kDa. The sequence matches with an expressed sequence tag of unknown function in public data bases, named TEM1, which was independently linked to tumor endothelium by serial analysis of gene expression profiling. Bioinformatic evaluation classifies endosialin as a C-type lectin-like protein, composed of a signal leader peptide, five globular extracellular domains (including a C-type lectin domain, one domain with similarity to the Sushi/ccp/scr pattern, and three EGF repeats), followed by a mucin-like region, a transmembrane segment, and a short cytoplasmic tail. Carbohydrate analysis shows that the endosialin core protein carries abundantly sialylated, O-linked oligosaccharides and is sensitive to O-sialoglycoprotein endopeptidase, placing it in the group of sialomucin-like molecules. The N-terminal 360 amino acids of endosialin show homology to thrombomodulin, a receptor involved in regulating blood coagulation, and to complement receptor C1qRp. This structural kinship may indicate a function for endosialin as a tumor endothelial receptor for as yet unknown ligands, a notion now amenable to molecular investigation.

Regulation of chromatin structure by site-specific histone H3 methyltransferases.

The organization of chromatin into higher-order structures influences chromosome function and epigenetic gene regulation. Higher-order chromatin has been proposed to be nucleated by the covalent modification of histone tails and the subsequent establishment of chromosomal subdomains by non-histone modifier factors. Here we show that human SUV39H1 and murine Suv39h1--mammalian homologues of Drosophila Su(var)3-9 and of Schizosaccharomyces pombe clr4--encode histone H3-specific methyltransferases that selectively methylate lysine 9 of the amino terminus of histone H3 in vitro. We mapped the catalytic motif to the evolutionarily conserved SET domain, which requires adjacent cysteine-rich regions to confer histone methyltransferase activity. Methylation of lysine 9 interferes with phosphorylation of serine 10, but is also influenced by pre-existing modifications in the amino terminus of H3. In vivo, deregulated SUV39H1 or disrupted Suv39h activity modulate H3 serine 10 phosphorylation in native chromatin and induce aberrant mitotic divisions. Our data reveal a functional interdependence of site-specific H3 tail modifications and suggest a dynamic mechanism for the regulation of higher-order chromatin.

brakeless is required for lamina targeting of R1-R6 axons in the Drosophila visual system.

Photoreceptors in the Drosophila eye project their axons retinotopically to targets in the optic lobe of the brain. The axons of photoreceptor cells R1-R6 terminate in the first optic ganglion, the lamina, while R7 and R8 axons project through the lamina to terminate in distinct layers of the second ganglion, the medulla. Here we report the identification of the gene brakeless (bks) and show that its function is required in the developing eye specifically for the lamina targeting of R1-R6 axons. In mosaic animals lacking bks function in the eye, R1-R6 axons project through the lamina to terminate in the medulla. Other aspects of visual system development appear completely normal: photoreceptor and lamina cell fates are correctly specified, R7 axons correctly target the medulla, and both correctly targeted R7 axons and mistargeted R1-R6 axons maintain their retinotopic order with respect to both anteroposterior and dorsoventral axes. bks encodes two unusually hydrophilic nuclear protein isoforms, one of which contains a putative C(2)H(2) zinc finger domain. Transgenic expression of either Bks isoform is sufficient to restore the lamina targeting of R1-R6 axons in bks mosaics, but not to retarget R7 or R8 axons to the lamina. These data demonstrate the existence of a lamina-specific targeting mechanism for R1-R6 axons in the Drosophila visual system, and provide the first entry point in the molecular characterization of this process.

Pds5 cooperates with cohesin in maintaining sister chromatid cohesion.

BACKGROUND: Sister chromatid cohesion depends on a complex called cohesin, which contains at least four subunits: Smc1, Smc3, Scc1 and Scc3. Cohesion is established during DNA replication, is partially dismantled in many, but not all, organisms during prophase, and is finally destroyed at the metaphase-to-anaphase transition. A quite separate protein called Spo76 is required for sister chromatid cohesion during meiosis in the ascomycete Sordaria. Spo76-like proteins are highly conserved amongst eukaryotes and a homologue in Aspergillus nidulans, called BimD, is required for the completion of mitosis. The isolation of the cohesin subunit Smc3 as a suppressor of BimD mutations suggests that Spo76/BimD might function in the same process as cohesin. RESULTS: We show here that the yeast homologue of Spo76, called Pds5, is essential for establishing sister chromatid cohesion and maintaining it during metaphase. We also show that Pds5 co-localizes with cohesin on chromosomes, that the chromosomal association of Pds5 and cohesin is interdependent, that Scc1 recruits Pds5 to chromosomes in G1 and that its cleavage causes dissociation of Pds5 from chromosomes at the metaphase-to-anaphase transition. CONCLUSIONS: Our data show that Pds5 functions as part of the same process as cohesin. Sequence similarities and secondary structure predictions indicate that Pds5 consists of tandemly repeated HEAT repeats, and might therefore function as a protein-protein interaction scaffold, possibly in the cohesin-DNA complex assembly.

Prediction of potential GPI-modification sites in proprotein sequences.

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.

Evaluation of human-readable annotation in biomolecular sequence databases with biological rule libraries.

MOTIVATION: Computer-based selection of entries from sequence databases with respect to a related functional description, e.g. with respect to a common cellular localization or contributing to the same phenotypic function, is a difficult task. Automatic semantic analysis of annotations is not only hampered by incomplete functional assignments. A major problem is that annotations are written in a rich, non-formalized language and are meant for reading by a human expert. This person can extract from the text considerably more information than is immediately apparent due to his extended biological background knowledge and logical reasoning. APPROACH: A technique of automated annotation evaluation based on a combination of lexical analysis and the usage of biological rule libraries has been developed. The proposed algorithm generates new functional descriptors from the annotation of a given entry using the semantic units of the annotation as prepositions for implications executed in accordance with the rule library. RESULTS: The prototype of a software system, the Meta_A(nnotator) program, is described and the results of its application to sequence attribute assignment and sequence selection problems, such as cellular localization and sequence domain annotation of SWISS-PROT entries, are presented. The current software version assigns useful subcellular localization qualifiers to approximately 88% of all SWISS-PROT entries. As shown by demonstrative examples, the combination of sequence and annotation analysis is a powerful approach for the detection of mutual annotation/sequence inconsistencies. AVAILABILITY: Results for the cellular localization assignment can be viewed at the URL http://www.bork. embl-heidelberg.de/CELL_LOC/CELL_LOC.html.

PSIC: profile extraction from sequence alignments with position-specific counts of independent observations.

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).

Predicting function: from genes to genomes and back.

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-based fold predictions for Mycoplasma genitalium proteins.

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.

Are knowledge-based potentials derived from protein structure sets discriminative with respect to amino acid types?

The parametric description of residue environments through solvent accessibility, backbone conformation, or pairwise residue-residue distances is the key to the comparison between amino acid types at protein sequence positions and residue locations in structural templates (condition of protein sequence-structure match). For the first time, the research results presented in this study clarify and allow to quantify, on a rigorous statistical basis, to what extent the amino acid type-specific distributions of commonly used environment parameters are discriminative with respect to the 20 amino acid types. Relying on the Bahadur theory, we estimate the probability of error in a single-sequence-structure alignment based on weak or absent discriminative power in a learning database of protein structure. We present the results for many residue environment variables and demonstrate that each fold description parameter is sensitive with respect to only a few amino acid types while indifferent to most of the other amino acid types. Even complex structural characteristics combining solvent-accessible surface area, backbone conformation, and pairwise distances distinguish only some amino acid types, whereas the others remain nondiscriminated. We find that the knowledge-based potentials currently in use treat especially Ala, Asp, Gln, His, Ser, Thr, and Tyr as essentially "average" amino acids. Thus, highly discriminative amino acid types define the alignment register in gapless sequence-structure alignments. The introduction of gaps leads to alignment ambiguities at sequence positions occupied by nondiscriminated amino acid types. Therefore, local sequence-structure alignments produced by techniques with gaps cannot be reliable. Conceptionally new and more sensitive environment parameters must be invented.