Search

Haptoglobin (Hp) is a hemoglobin-binding plasma protein consisting of two types of chains, called alpha and beta, which originate from a common polypeptide. In humans, but not in other mammals, Hp has been shown to occur in two allelic forms, Hp1 and Hp2, which differ in the length of the alpha-chain. The longer alpha-chain (in Hp2) seems to have arisen by an internal duplication of a gene segment coding for almost the entire alpha-chain of Hp1. In this article we show that Hp of cow (Bos taurus) contains an alpha-chain, the structure of which is similar to that of the human Hp2 alpha-chain. Furthermore, comparison of the structure of bovine Hp and human Hp2 suggests that the bovine gene arose by a duplication of the gene segment homologous to that duplicated in human Hp2. However, a phylogenetic analysis indicates that the two genes were formed independently. The evolutionary pressure that has led to the fixation of the Hps with a longer alpha-chain is not known.

Functionally linked genes generally evolve at similar rates and the knowledge of this particular feature of genomic evolution has been used as the basis for the phylogenetic profiling method. We illustrate here an exception to this rule in the evolution of the cytochrome bd complex. This is a two-component oxidase complex, with the subunits I and II known to be widely present in bacteria. The subunits within the cytochrome bd complex are under the same evolutionary pressure and most likely behave in the same evolutionary manner. However, the sequence similarity of genes encoding subunit II varies considerably across species. Genes encoding subunit II evolve 1.2 times faster on most of the branches of their phylogeny than subunit I genes. Furthermore, the genes encoding subunit II in Oceanobacillus iheyensis, Bacillus halodurans, and Staphylococcus species do not have detectable homologues within E. coli due to their large divergence. Together, the two subunits of cytochrome bd reveal an interesting example of an asymmetric pattern of evolutionary change.

Fasciola spp., commonly known as liver fluke, are significant trematode parasites of livestock and humans. They secrete several cathepsin L-like cysteine proteases, some of which differ in enzymatic properties and timing of expression in the parasite's life cycle. A detailed sequence and evolutionary analysis is presented, based on 18 cathepsin L-like enzymes isolated from Fasciola spp. (including a novel clone identified in this study). The enzymes form a monophyletic group which has experienced several gene duplication events over the last approximately 135 million years, giving rise to the present-day enzymatic repertoire of the parasite. This timing of these duplications appears to correlate with important points in the evolution of the mammalian hosts. Furthermore, the dates suggest that Fasciola hepatica and Fasciola gigantica diverged around 19 million years ago. A novel analysis, based on the pattern of amino acid diversity, was used to identify sites in the enzyme that are predicted to be subject to positive adaptive evolution. Many of these sites occur within the active site cleft of the enzymes, and hence would be expected to lead to differences in substrate specificity. Using homology modeling, with reference to previously obtained biochemical data, we are able to predict S2 subsite specificity for these enzymes: specifically those that can accommodate bulky hydrophobic residues in the P2 position and those that cannot. A number of other positions subject to evolutionary pressure and potentially significant for enzyme function are also identified, including sites anticipated to diminish cystatin binding affinity.

The common assumption of operons as composed of genes that cooperate in a biological process is confirmed here by showing that Escherichia coli operons tend to be composed of genes that belong to the same general class of cellular function. Furthermore, the comparison between the genomic organization of E. coli and that of Bacillus subtilis shows that the genes that are homologous to genes that belong to experimentally characterized E. coli operons tend to cluster in neighboring regions of the genome. This tendency is greater for the subset of E. coli operons whose genes belong to a single functional class. These observations indicate strong evolutionary pressure that, translated into functional constraints, leads to the inclusion of many essential functions in conserved operons and clusters in these two distant species.

Recent evidence suggests that gamete recognition proteins may be subjected to directed evolutionary pressure that enhances sequence variability. We evaluated whether diversity enhancing selection is operating on a marine invertebrate fertilization protein by examining the intraspecific DNA sequence variation of a 273-base pair region located at the 5' end of the sperm bindin locus in 134 adult red sea urchins (Strongylocentrotus franciscanus). Bindin is a sperm recognition protein that mediates species-specific gamete interactions in sea urchins. The region of the bindin locus examined was found to be polymorphic with 14 alleles. Mean pairwise comparison of the 14 alleles indicates moderate sequence diversity (p-distance = 1.06). No evidence of diversity enhancing selection was found. It was not possible to reject the null hypothesis that the sequence variation observed in S. franciscanus bindin is a result of neutral evolution. Statistical evaluation of expected proportions of replacement and silent nucleotide substitutions, observed versus expected proportions of radical replacement substitutions, and conformance to the McDonald and Kreitman test of neutral evolution all indicate that random mutation followed by genetic drift created the polymorphisms observed in bindin. Observed frequencies were also highly similar to results expected for a neutrally evolving locus, suggesting that the polymorphism observed in the 5' region of S. franciscanus bindin is a result of neutral evolution.

Using a maximum-likelihood formalism, we have developed a method with which to reconstruct the sequences of ancestral proteins. Our approach allows the calculation of not only the most probable ancestral sequence but also of the probability of any amino acid at any given node in the evolutionary tree. Because we consider evolution on the amino acid level, we are better able to include effects of evolutionary pressure and take advantage of structural information about the protein through the use of mutation matrices that depend on secondary structure and surface accessibility. The computational complexity of this method scales linearly with the number of homologous proteins used to reconstruct the ancestral sequence.

Structural transitions in oligomeric proteins due to ligand binding are important in biomolecular regulatory processes. The transitions may occur on the secondary, tertiary or quarternary structure levels. Detailed consideration of the time sequence of ligand binding to the oligomer shows that there is an intrinsic dynamic asymmetry in all oligomer transitions, even if the initial and the final state are completely symmetric. This asymmetry has important bearing on the evolution and the divergence of the primary structure (amino acid sequence) of oligomeric proteins. It may explain (at least in part) the occurrence of oligomeric proteins with similar but not identical protomers. Certain specific groups of oligomers are shown to be under greater evolutionary pressure for protomer structure divergence. The dynamic asymmetry of oligomer transitions also results in higher complexity in reaction kinetics. Some implications on ribosome structural evolution are discussed.