Evolution of Structural Shape in Bacterial Globin-Related Proteins

A model for protein sequence evolution based on selective pressure for protein stability: application to hemoglobins

Negative selection against protein instability is a central influence on evolution of proteins. Protein stability is maintained over evolution despite changes in underlying sequences. An empirical all-site stability-based model of evolution was developed to focus on the selection of residues arising from their contributions to protein stability. In this model, site rates could vary. A structure-based method was used to predict stationary frequencies of hemoglobin residues based on their propensity to promote protein stability at a site. Sites with destabilizing residues were shown to change more rapidly in hemoglobins than sites with stabilizing residues. For diverse proteins the results were consistent with stability-based selection. Maximum likelihood studies with hemoglobins supported the stability-based model over simple Poisson-based methods. These observations are consistent with suggestions that purifying selection to maintain protein structural stability plays a dominant role in protein evolution.

Spatial autocorrelation of amino Acid replacement rates in the vasopressin receptor family

Evolutionary rates of sites can be independent of one another or correlated in some fashion. Significant spatial autocorrelation was observed for site amino acid replacement rates in vasopressin receptor family proteins (VPRs). Spatial autocorrelation of rates is the propensity of residues to lie near other residues of similar rate in the folded protein structure. Optimal correlation occurred at a distance suggesting that residues in contact had correlated rates. As another way to study the same phenomenon, VPR was partitioned into >40 x 10 A(3) contiguous spatial clusters for amino acid replacement rate estimation. Partitioning was done without preconception of functional regions of the protein and with a random partition control. Cluster rates exhibited an overdispersed distribution suggesting that rates were not randomly distributed in the spatial partitions. In tests, cluster partitioning improved maximum likelihood and Bayesian likelihood models for VPR evolution. Spatial clusters with outlier rates, or lineage-specific clusters differing in rate, proved to contain VPR features likely to be under selection. Thus the spatial autocorrelation observed is probably not just a statistical finding, but likely has an evolutionary basis in protein function.

A model of globin evolution

Putative globins have been identified in 426 bacterial, 32 Archaeal and 67 eukaryote genomes. Among these sequences are the hitherto unsuspected presence of single domain sensor globins within Bacteria, Fungi, and a Euryarchaeote. Bayesian phylogenetic trees suggest that their occurrence in the latter two groups could be the result of lateral gene transfer from Bacteria. Iterated psiblast searches based on groups of globin sequences indicate that bacterial flavohemoglobins are closer to metazoan globins than to the other two lineages, the 2-over-2 globins and the globin-coupled sensors. Since Bacteria is the only kingdom to have all the subgroups of the three globin lineages, we propose a working model of globin evolution based on the assumption that all three lineages originated and evolved only in Bacteria. Although the 2-over-2 globins and the globin-coupled sensors recognize flavohemoglobins, there is little recognition between them. Thus, in the first stage of globin evolution, we favor a flavohemoglobin-like single domain protein as the ancestral globin. The next stage comprised the splitting off to single domain 2-over-2 and sensor-like globins, followed by the covalent addition of C-terminal domains resulting in the chimeric flavohemoglobins and globin-coupled sensors. The last stage encompassed the lateral gene transfers of some members of the three globin lineages to specific groups of Archaea and Eukaryotes.