GeneSeeker: extraction and integration of human disease-related information from web-based genetic databases.

The identification of genes underlying human genetic disorders requires the combination of data related to cytogenetic localization, phenotypes and expression patterns, to generate a list of candidate genes. In the field of human genetics, it is normal to perform this combination analysis by hand. We report on GeneSeeker (http://www.cmbi.ru.nl/GeneSeeker/), a web server that gathers and combines data from a series of databases. All database searches are performed via the web interfaces provided with the original databases, guaranteeing that the most recent data are queried, and obviating data warehousing. GeneSeeker makes the same selection of candidate genes as the human geneticists would have performed, and thus reducing the time-consuming process to a few minutes. GeneSeeker is particularly well suited for syndromes in which the disease gene displays altered expression patterns in the affected tissue(s).

Cloning and sequence analysis of a hypothalamic cDNA encoding a D1c dopamine receptor in tilapia.

Physiological and pharmacological studies have indicated that during acid stress a D1-like dopamine receptor becomes functional on intermediate pituitary melanocyte-stimulating hormone cells of tilapia (Oreochromis mossambicus). As a first step towards physiological expression studies we isolated a D1-like dopamine receptor from a tilapia hypothalamus cDNA library. Construction of a phylogenetic tree of most of the D1-like receptors known in human, rat, Xenopus, goldfish and Drosophila revealed that the here presented clone is most likely the tilapia equivalent of the Xenopus D1c dopamine receptor.

Average orientation of aromatic residues in proteins determined from linear dichroism spectroscopy. A comparison of results on bovine gamma-crystallins with X-ray data.

The structures of the two very closely related proteins, bovine gamma II- and gamma IVa-crystallin have been studied by means of near-ultra-violet linear dichroism spectroscopy on squeezed polyacrylamide gel systems. The crystallin spectra are discussed in terms of the spectra of the aromatic chromophores present in these proteins and provide detailed information on the average orientation of these residues in the proteins. A comparison of our results with information based on crystallographic X-ray experiments shows excellent agreement, reflecting even some of the minor differences between the two proteins studied. Since linear dichroism measurements as performed here take a few days only, and can be done on most aqueous protein solutions, linear dichroism spectroscopy may give a valuable contribution to structural studies on proteins.

Evidence for limited within-person evolution of the V3 domain of the HIV-1 envelope in the amsterdam population.

OBJECTIVE: To study the development of the V3 region of the HIV-1 envelope over time, both within subjects and population-wide. METHODS: Direct V3 sequences were obtained from viral RNA from seroconversion samples of 138 individuals [32 intravenous drug users (IVDU), 106 homosexual men], as well as from 5-year follow-up samples of 45 of these individuals (11 IVDU, 34 homosexual men). RESULTS: The population-wide variation of the V3 region in both the seroconversion samples and the 5-year samples steadily increased over consecutive years and were of similar magnitude in each calendar year. The variation in the IVDU group was slightly lower (presumably lagging behind) than in the homosexual group, but also increased over time. The consensus sequence, representing the centre of the swarm of variants, remained almost stationary in 10 years of evolution. The V3 sequences from virions in serum collected 5 years after seroconversion still resembled those from the seroconversion sample, either in overall similarity or in specific (signature) amino acids. Seroconversion and late sequences from a donor-recipient pair were also very similar. CONCLUSIONS: The variation in V3 sequences from seroconversion samples is as large as that in 5-year follow-up samples from the same calendar year, suggesting that there is no strong selection for a particular V3 genotype at transmission. The HIV-1 subtype B quasispecies in a naive population appears to evolve through unbiased expansion around a stationary consensus sequence. Despite its large variability, the V3 region retains many of its individual characteristics after 5 years of infection. This indicates that the sampling moment (relative to the seroconversion data) will not greatly influence the results of phylogenetic analyses.

Subtilases: the superfamily of subtilisin-like serine proteases.

Subtilases are members of the clan (or superfamily) of subtilisin-like serine proteases. Over 200 subtilases are presently known, more than 170 of which with their complete amino acid sequence. In this update of our previous overview (Siezen RJ, de Vos WM, Leunissen JAM, Dijkstra BW, 1991, Protein Eng 4:719-731), details of more than 100 new subtilases discovered in the past five years are summarized, and amino acid sequences of their catalytic domains are compared in a multiple sequence alignment. Based on sequence homology, a subdivision into six families is proposed. Highly conserved residues of the catalytic domain are identified, as are large or unusual deletions and insertions. Predictions have been updated for Ca(2+)-binding sites, disulfide bonds, and substrate specificity, based on both sequence alignment and three-dimensional homology modeling.

ProSearch: fast searching of protein sequences with regular expression patterns related to protein structure and function.

This report describes the implementation of ProSearch, a computer program that can efficiently search for motifs in protein sequences. ProSearch currently uses motifs that are contained in the PROSITE database, but user-developed patterns can easily be added to any search. ProSearch can generate a report identifying the patterns present in a given protein sequence, their locations and, if desired, a short description of the identified patterns. The program is written in AWK (a small interpreted computer language), which can run on all computer platforms commonly found in laboratories. ProSearch can search a 348-amino acid protein for 690 patterns in less than 5 s on a typical workstation.

Cloning and expression of two proopiomelanocortin mRNAs in the common carp (Cyprinus carpio L.).

Proopiomelanocortin (POMC) is the precursor for a number of biologically active peptides such as adrenocorticotropic hormone (ACTH), alpha-melanocyte-stimulating hormone (alpha-MSH) and beta-endorphin. It is well known that these peptides are involved in the stress response in fish as well as in mammals. We have cloned two different carp POMC cDNAs called, POMC-I and POMC-II. The nucleotide sequences of 955 bp for POMC-I and 959 bp for POMC-II share 93.5% identity in their cDNAs, and the deduced amino acid sequences (both 222 amino acids) are 91.4% identical. In the ACTH and beta-MSH domain, two amino acid substitutions are found, whereas alpha-MSH and beta-endorphin are identical. For beta-MSH, the serine replacement (in POMC-I) by a glycine (in POMC-II) results in a putative amidation site Pro-X-Gly for POMC-II. We used RT-PCR to show that both POMC mRNAs are expressed in the hypophysis, hypothalamus and other parts of the brain of a single fish. Furthermore, in a phylogenetic tree based on POMC sequences the divergence of carp POMC-I and -II from tetraploid animals (salmon, trout and xenopus) is demonstrated.

Oreochromis mossambicus (tilapia) corticotropin-releasing hormone: cDNA sequence and bioactivity.

Although hypothalamic corticotropin-releasing hormone (CRH) is involved in the stress response in all vertebrate groups, only a limited number of studies on this neuroendocrine peptide deals with non-mammalian neuroendocrine systems. We determined the cDNA sequence of the CRH precursor of the teleost Oreochromis mossambicus (tilapia) and studied the biological potency of the CRH peptide in a homologous teleost bioassay. Polymerase chain reaction (PCR) with degenerate and specific primers yielded fragments of tilapia CRH cDNA. Full-length CRH cDNA (988 nucleotides) was obtained by screening a tilapia hypothalamus cDNA library with the tilapia CRH PCR products. The precursor sequence (167 amino acids) contains a signal peptide, the CRH peptide and a motif conserved among all vertebrate CRH precursors. Tilapia CRH (41 aa) displays between 63% and 80% amino acid sequence identity to CRH from other vertebrates, whereas the degree of identity to members of the urotensin I/urocortin lineage is considerably lower. In a phylogenetic tree, based on alignment of all full CRH peptide precursors presently known, the three teleost CRH precursors (tilapia; sockeye salmon, Oncorhynchus nerka; white sucker, Catostomus commersoni) form a monophyletic group distinct from amphibian and mammalian precursors. Despite the differences between the primary structures of tilapia and rat CRH, maximally effective concentrations of tilapia and rat CRH were equally potent in stimulating adrenocorticotropic hormone (ACTH) and alpha-MSH release by tilapia pituitaries in vitro. The tilapia and salmon CRH sequences show that more variation exists between orthologous vertebrate CRH structures, and teleost CRHs in particular than previously recognized. Whether the structural differences reflect different mechanisms of action of this peptide in the stress response remains to be investigated.

Genealogy of the alpha-crystallin--small heat-shock protein superfamily.

Sequences of 40 very diverse representatives of the alpha-crystallin-small heat-shock protein (alpha-Hsp) superfamily are compared. Their characteristic C-terminal 'alpha-crystallin domain' of 80-100 residues contains short consensus sequences that are highly conserved from prokaryotes to eukaryotes. There are, in addition, some positions that clearly distinguish animal from non-animal alpha-Hsps. The alpha-crystallin domain is predicted to consist of two hydrophobic beta-sheet motifs, separated by a hydrophilic region which is variable in length. Combination of a conserved alpha-crystallin domain with a variable N-terminal domain and C-terminal extension probably modulates the properties of the various alpha-Hsps as stress-protective and structural oligomeric proteins. Phylogeny reconstruction indicates that multiple alpha-Hsps were already present in the last common ancestor of pro- and eukaryotes. It is suggested that during eukaryote evolution, animal and non-animal alpha-Hsps originated from different ancestral gene copies. Repeated gene duplications gave rise to the multiple alpha-Hsps present in most organisms.

Divergence in codon usage of Lactobacillus species.

We have analyzed codon usage patterns of 70 sequenced genes from different Lactobacillus species. Codon usage in lactobacilli is highly biased. Both inter-species and intra-species heterogeneity of codon usage bias was observed. Codon usage in L. acidophilus is similar to that in L. helveticus, but dissimilar to that in L. bulgaricus, L. casei, L. pentosus and L. plantarum. Codon usage in the latter three organisms is not significantly different, but is different from that in L. bulgaricus. Inter-species differences in codon usage can, at least in part, be explained by differences in mutational drift. L. bulgaricus shows GC drift, whereas all other species show AT drift. L. acidophilus and L. helveticus rarely use NNG in family-box (a set of synonymous) codons, in contrast to all other species. This result may be explained by assuming that L. acidophilus and L. helveticus, but not other species examined, use a single tRNA species for translation of family-box codons. Differences in expression level of genes are positively correlated with codon usage bias. Highly expressed genes show highly biased codon usage, whereas weakly expressed genes show much less biased codon usage. Codon usage patterns at the 5'-end of Lactobacillus genes is not significantly different from that of entire genes. The GC content of codons 2-6 is significantly reduced compared with that of the remainder of the gene. The possible implications of a reduced GC content for the control of translation efficiency are discussed.

Phylogenetic trees constructed from hydrophobicity values of protein sequences.

Information about conformational properties of a protein is contained in the hydrophobicity values of the amino acids in its primary sequence. We have investigated the possibility of extracting meaningful evolutionary information from the comparison of the hydrophobicity values of the corresponding amino acids in the sequences of homologous proteins. Distance matrices for six families of homologous proteins were made on the basis of the differences in hydrophobicity values of the amino acids. The phylogenetic trees constructed from such matrices were at least as good (as judged from their faithful reflection of evolutionary relationships), as trees constructed from the usual minimum mutation distance matrix.

The expanding small heat-shock protein family, and structure predictions of the conserved "alpha-crystallin domain".

The ever-increasing number of proteins identified as belonging to the family of small heat-shock proteins (shsps) and alpha-crystallins enables us to reassess the phylogeny of this ubiquitous protein family. While the prokaryotic and fungal representatives are not properly resolved, most of the plant and animal shsps and related proteins are clearly grouped in distinct clades, reflecting a history of repeated gene duplications. The members of the shsp family are characterized by the presence of a conserved homologous "alpha-crystallin domain," which sometimes is present in duplicate. Predictions are made of secondary structure and solvent accessibility of this domain, which together with hydropathy profiles and intron positions support the presence of two similar hydrophobic beta-sheet-rich motifs, connected by a hydrophilic alpha-helical region. Together with an overview of the newly characterized members of the shsp family, these data help to define this family as being involved as stable structural proteins and as molecular chaperones during normal development and induced under pathological and stressful conditions.

Homology modelling and protein engineering strategy of subtilases, the family of subtilisin-like serine proteinases.

Subtilases are members of the family of subtilisin-like serine proteases. Presently, greater than 50 subtilases are known, greater than 40 of which with their complete amino acid sequences. We have compared these sequences and the available three-dimensional structures (subtilisin BPN', subtilisin Carlsberg, thermitase and proteinase K). The mature enzymes contain up to 1775 residues, with N-terminal catalytic domains ranging from 268 to 511 residues, and signal and/or activation-peptides ranging from 27 to 280 residues. Several members contain C-terminal extensions, relative to the subtilisins, which display additional properties such as sequence repeats, processing sites and membrane anchor segments. Multiple sequence alignment of the N-terminal catalytic domains allows the definition of two main classes of subtilases. A structurally conserved framework of 191 core residues has been defined from a comparison of the four known three-dimensional structures. Eighteen of these core residues are highly conserved, nine of which are glycines. While the alpha-helix and beta-sheet secondary structure elements show considerable sequence homology, this is less so for peptide loops that connect the core secondary structure elements. These loops can vary in length by greater than 150 residues. While the core three-dimensional structure is conserved, insertions and deletions are preferentially confined to surface loops. From the known three-dimensional structures various predictions are made for the other subtilases concerning essential conserved residues, allowable amino acid substitutions, disulphide bonds, Ca(2+)-binding sites, substrate-binding site residues, ionic and aromatic interactions, proteolytically susceptible surface loops, etc. These predictions form a basis for protein engineering of members of the subtilase family, for which no three-dimensional structure is known.

Dogfish alpha-crystallin sequences. Comparison with small heat shock proteins and Schistosoma egg antigen.

The amino acid sequences of the alpha-crystallin A and B chains of the dogfish, Squalus acanthias, have been determined. Comparison with alpha-crystallins from other species reveals that charged amino acid replacements have been strongly avoided in the evolution of this lens protein. The homology of alpha-crystallins with the small heat shock proteins is pronounced throughout the major part of the proteins, starting from the position of the first intron in the alpha-crystallin genes, but is also detectable in the amino-terminal sequences of human, Xenopus, and Drosophila small heat shock proteins. In addition, a remarkable short sequence similarity is present only in the amino termini of dogfish alpha B and Drosophila HSP22. The Schistosoma egg antigen p40 turns out to have a tandemly repeated region of homology with the common sequence domain of alpha-crystallins and small heat shock proteins. Comparison of hydropathy profiles indicates the conservation of conformation of the common domains in these three families of proteins. Construction of phylogenetic trees suggests that the alpha A and alpha B genes apparently originated from a single ancestral small heat shock protein gene and indicates that introns have been lost during the evolution of the heat shock protein genes.

Evolution of the alpha-crystallin/small heat-shock protein family.

The common characteristic of the alpha-crystallin/small heat-shock protein family is the presence of a conserved homologous sequence of 90-100 residues. Apart from the vertebrate lens proteins--alpha A- and alpha B-crystallin--and the ubiquitous group of 15-30-kDa heat-shock proteins, this family also includes two mycobacterial surface antigens and a major egg antigen of Schistosoma mansoni. Multiple small heat-shock proteins are especially present in higher plants, where they can be distinguished in at least two classes of cytoplasmic proteins and a chloroplast-located class. The alpha-crystallins have recently been found in many tissues outside the lens, and alpha B-crystallin, in particular, behaves in many respects like a small heat-shock protein. The homologous sequences constitute the C-terminal halves of the proteins and probably represent a structural domain with a more variable C-terminal extension. These domains must be responsible for the common structural and functional properties of this protein family. Analysis of the phylogenetic tree and comparison of the biological properties of the various proteins in this family suggest the following scenario for its evolution: The primordial role of the small heat-shock protein family must have been to cope with the destabilizing effects of stressful conditions on cellular integrity. The alpha-crystallin-like domain appears to be very stable, which makes it suitable both as a surface antigen in parasitic organisms and as a long-living lens protein in vertebrates. It has recently been demonstrated that, like the other heat-shock proteins, the alpha-crystallins and small heat-shock proteins function as molecular chaperones, preventing undesired protein-protein interactions and assisting in refolding of denatured proteins. Many of the small heat-shock proteins are differentially expressed during normal development, and there is good evidence that they are involved in cytomorphological reorganizations and in degenerative diseases. In conjunction with the stabilizing, thermoprotective role of alpha-crystallins and small heat-shock proteins, they may also be involved in signal transduction. The reversible phosphorylation of these proteins appears to be important in this respect.

Molecular evolution of interleukin-3.

Chimpanzee, tamarin, and marmoset interleukin-3 (IL-3) genes were cloned, sequenced, and expressed. Western blot analysis demonstrated that functional genes were isolated. IL-3 sequences were compared with those of mouse, rat, rhesus monkey, gibbon, and man. Multiple alignment of the IL-3 coding regions showed that only a few regions had been conserved during mammalian evolution, which are likely associated with functional domains of the IL-3 protein. Substitution rates for the various lineages were calculated and the numbers of synonymous and nonsynonymous substitutions were estimated separately. Distance matrices of the IL-3 coding regions were used to construct phylogenetic trees which revealed large differences in IL-3 evolution rate as well as a more rapid substitution rate for rodents and a rate slowdown during hominoid evolution. Extremes were rhesus monkey IL-3, which accumulated few synonymous substitutions, and gibbon IL-3, which had almost exclusively synonymous substitutions. In rhesus monkey IL-3, nonsynonymous substitutions outnumbered synonymous substitutions, which could not be readily explained by a random process of substitutions. We assume that during evolution of IL-3, the majority of the amino acid replacements and the impaired interspecies functional cross-reactivity originate from selection mechanisms with the most likely selective force being the structure of the heterodimeric IL.3 cell-surface receptor. Insight into IL-3 architecture and structural analysis of the IL-3 receptor are needed to analyze the unusually fast evolution of IL-3 in more detail.

Ratites as oldest offshoot of avian stem--evidence from alpha-crystallin A sequences.

One of the most disputed issues in avian phylogeny is the origin of the ratites, the large flightless birds of the Southern Hemisphere (reviewed in refs 1-3). It is still not generally agreed whether the ostriches, rheas, emus and cassowaries, and probably kiwis, form a natural, monophyletic group, although much recent evidence supports this view. Also, their phylogenetic relationship with the other avian orders remains unresolved, comparative protein sequence studies might shed new light on this problem. Therefore, we determined the amino acid sequence of the eye lens protein alpha-crystallin A in ostrich, rhea and emu, and in representatives of 13 other avian orders. Comparison of these sequences with known alpha A sequences of mammals, reptiles, frog and dogfish provides strong evidence that the ratites, as a monophyletic assemblage, represent the first offshoot of the avain line.

Distinct frequency-distributions of homopolymeric DNA tracts in different genomes.

The unusual base composition of the genome of the human malaria parasite Plasmodium falciparum prompted us to systematically investigate the occurrence of homopolymeric DNA tracts in the P. falciparum genome and, for comparison, in the genomes of Homo sapiens , Saccharomyces cerevisiae , Caenorhabditis elegans , Arabidopsis thaliana , Escherichia coli and Mycobacterium tuberculosis. Comparison of theobserved frequencies with the frequencies as expected for random DNA revealed that homopolymeric (dA:dT) tracts occur well above chance in the eukaryotic genome. In the majority of these genomes, (dA:dT) tract overrepresentation proved to be an exponential function of the tract length. (dG:dC) tract overrepresentation was absent or less pronounced in both prokaryotic and eukaryotic genomes. On the basis of our results, we propose that homopolymeric (dA:dT) tracts are expanded via replication slippage. This slippage-mediated expansion does not operate on tracts with lengths below a critical threshold of 7-10 bp.

Extreme differences in charge changes during protein evolution.

The maintenance of a proper distribution of charged amino acid residues might be expected to be an important factor in protein evolution. We therefore compared the inferred changes in charge during the evolution of 43 protein families with the changes expected on the basis of random base substitutions. It was found that certain proteins, like the eye lens crystallins and most histones, display an extreme avoidance of changes in charge. Other proteins, like phospholipase A2 and ferredoxin, apparently have sustained more charged replacements than expected, suggesting a positive selection for changes in charge. Depending on function and structure of a protein, charged residues apparently can be important targets for selective forces in protein evolution. It appears that actual biased codon usage tends to decrease the proportion of charged amino acid replacements. The influence of nonrandomness of mutations is more equivocal. Genes that use the mitochondrial instead of the universal code lower the probability that charge changes will occur in the encoded proteins.