Electric charge divergence in proteins: insights into the evolution of their three-dimensional properties

Previous studies of protein evolution have identified important mutations in various proteins that affect a small number of residues, but dramatically alter protein function. However, the evolutionary process underlying the three-dimensional protein properties, which are determined by a much larger number of residues, remains unclear. Based on a comparative evolutionary analysis of teleost phosphoglucose isomerases (PGIs; EC 5.3.1.9), we previously demonstrated that the relatively weak selection on many amino acid sites has played an important role in the evolution of protein electric charge as a model of three-dimensional protein properties. To ascertain the generality of this finding, we sought further evidence of this type of protein evolution. For this purpose, we analyzed the vertebrate isoforms of fructose-1,6-bisphosphate aldolase (ALD; EC 4.1.2.13), for which electric charges are known to have diverged after gene duplication. The results showed that the divergence in electric charge between the ALD isoforms was also driven by weak selection on many amino acid sites, as in PGI, confirming the generality of earlier findings. To obtain further insights, ALD and PGI were compared to the proteins pancreatic ribonuclease (EC 3.1.27.5) and triose-phosphate isomerase (EC 5.3.1.1), for which electric charges likely evolved through a well-defined mode of molecular evolution; namely, strong selection on specific amino acid sites. Comparison of the number and composition of amino acids on the protein surface suggested that the absolute number of evolutionarily changeable amino acids in a protein affects the strength of selection pressure acting on individual amino acid sites.     

Rates of Conservative and Radical Nonsynonymous Nucleotide Substitutions in Mammalian Nuclear Genes

To understand the process and mechanism of protein evolution, it is important to know what types of amino acid substitutions are more likely to be under selection and what types are mostly neutral. An amino acid substitution can be classified as either conservative or radical, depending on whether it involves a change in a certain physicochemical property of the amino acid. Assuming Kimura's two-parameter model of nucleotide substitution, I present a method for computing the numbers of conservative and radical nonsynonymous (amino acid altering) nucleotide substitutions per site and estimate these rates for 47 nuclear genes from mammals. The results are as follows. (1) The average radical/conservative rate ratio is 0.81 for charge changes, 0.85 for polarity changes, and 0.49 when both polarity and volume changes are considered. (2) The radical/conservative rate ratio is positively correlated with the nonsynonymous/synonymous rate ratio for charge changes or when both polarity and volume changes are considered. (3) Both the conservative/synonymous rate ratio and the radical/synonymous rate ratio are lower in the rodent lineage than in the primate or artiodactyl lineage, suggesting more intense purifying selection in the rodent lineage, for both conservative and radical nonsynonymous substitutions. (4) Neglecting transition/transversion bias would cause an underestimation of both radical and conservative rates and the ratio thereof. (5) Transversions induce more dramatic genetic alternations than transitions in that transversions produce more amino acid altering changes and among which, more radical changes. Received: 6 April 1999 / Accepted: 16 August 1999     

Phylogenetic network and physicochemical properties of nonsynonymous mutations in the protein-coding genes of human mitochondrial DNA

Theories on molecular evolution predict that phylogenetically recent nonsynonymous mutations should contain more non-neutral amino acid replacements than ancient mutations. We analyzed 840 complete coding-region human mitochondrial DNA (mtDNA) sequences for nonsynonymous mutations and evaluated the mutations in terms of the physicochemical properties of the amino acids involved. We identified 465 distinct missense and 6 nonsense mutations. 48% of the amino acid replacements changed polarity, 26% size, 8% charge, 32% aliphaticity, 13% aromaticity, and 44% hydropathy. The reduced-median networks of the amino acid changes revealed relatively few differences between the major continent-specific haplogroups, but a high variation and highly starlike phylogenies within the haplogroups. Some 56% of the mutations were private, and 25% were homoplasic. Nonconservative changes were more common than expected among the private mutations but less common among the homoplasic mutations. The asymptotic maximum of the number of nonsynonymous mutations in European mtDNA was estimated to be 1,081. The results suggested that amino acid replacements in the periphery of phylogenetic networks are more deleterious than those in the central parts, indicating that purifying selection prevents the fixation of some alleles.     

Conservation of engrailed-like homeobox sequences during vertebrate evolution

The Drosophila melanogaster developmental gene engrailed (en) is a member of a distinct subfamily of homeobox genes with a wide phylogenetic distribution. Here we report the use of reduced stringency polymerase chain reaction (PCR) to amplify and clone 8 genes related to en from 5 vertebrate species, including representatives of the most ancient vertebrate lineages. Nucleotide and deduced amino acid sequence comparisons between mouse, toad, zebrafish, lamprey and hagfish genes reveal extensive evolutionary conservation, and suggests that 2 en-like genes have been retained in most vertebrate lineages.     

What Amino Acid Properties Affect Protein Evolution?

We studied 10 protein-coding mitochondrial genes from 19 mammalian species to evaluate the effects of 10 amino acid properties on the evolution of the genetic code, the amino acid composition of proteins, and the pattern of nonsynonymous substitutions. The 10 amino acid properties studied are the chemical composition of the side chain, two polarity measures, hydropathy, isoelectric point, volume, aromaticity, aliphaticity, hydrogenation, and hydroxythiolation. The genetic code appears to have evolved toward minimizing polarity and hydropathy but not the other seven properties. This can be explained by our finding that the presumably primitive amino acids differed much only in polarity and hydropathy, but little in the other properties. Only the chemical composition (C) and isoelectric point (IE) appear to have affected the amino acid composition of the proteins studied, that is, these proteins tend to have more amino acids with typical C and IE values, so that nonsynonymous mutations tend to result in small differences in C and IE. All properties, except for hydroxythiolation, affect the rate of nonsynonymous substitution, with the observed amino acid changes having only small differences in these properties, relative to the spectrum of all possible nonsynonymous mutations. Received: 2 January 1998 / Accepted: 25 April 1998     

From DNA to fitness differences: sequences and structures of adaptive variants of Colias phosphoglucose isomerase (PGI)

Colias eurytheme butterflies display extensive allozyme polymorphism in the enzyme phosphoglucose isomerase (PGI). Earlier studies on biochemical and fitness effects of these genotypes found evidence of strong natural selection maintaining this polymorphism in the wild. Here we analyze the molecular features of this polymorphism by sequencing multiple alleles and modeling their structures. PGI is a dimer with rotational symmetry. Each monomer provides a critical residue to the other monomer's catalytic center. Sequenced alleles differ at multiple amino acid positions, including cryptic charge-neutral variation, but most consistent differences among the electromorph alleles are at the charge-changing amino acid sites. Principal candidate sites of selection, identified by structural and functional analyses and by their variants' population frequencies, occur in interpenetrating loops across the interface between monomers, where they may alter subunit interactions and catalytic center geometry. Comparison to a second (and basal) species, Colias meadii, also polymorphic for PGI under natural selection, reveals one fixed amino acid difference between their PGIs, which is located in the interpenetrating loop and accompanies functional differences among their variants. We also study nucleotide variability among the PGI alleles, comparing these data to similar data from another glycolytic enzyme gene, glyceraldehyde-3-phosphate dehydrogenase. Despite extensive nonsynonymous and synonymous polymorphism at PGI in each species, the only base changes fixed between species are the two causing the amino acid replacement; this absence of synonymous fixation yields a significant McDonald-Kreitman test. Analyses of these data suggest historical population expansion. Positive peaks of Tajima's D statistic, representing regions of neutral "hitchhiking," are found around the principal candidate sites of selection. This study provides novel views of molecular-structural mechanisms, and beginnings of historical evidence, for a long-persistent balanced enzyme polymorphism at PGI in these and perhaps other species.     

Accelerated evolution of small serum proteins (SSPs)-The PSP94 family proteins in a Japanese viper

Five small serum proteins (SSPs) with molecular masses of 6.5-10 kDa were detected in Habu (Trimeresurus flavoviridis) serum; this included two novel proteins SSP-4 and SSP-5. The amino acid sequences of these proteins and of SSP-1, SSP-2, and SSP-3, which were reported previously, were determined on the basis of the nucleotide sequences of their cDNAs. Although these proteins exhibited only limited sequence identity to mammalian prostatic secretory protein of 94 amino acids (PSP94), the topological pattern of disulfide bonds in SSPs was identical to that of the mammalian proteins. SSP-3 and SSP-4 lacked approximately 30 residues at the C-terminal. Each of the full-length cDNAs encoded a mature protein of 62-90 residues and a highly conserved signal peptide. The evolutionary distances between SSPs estimated on the basis of the amino acid changes were significantly greater than those of the synonymous nucleotide substitutions; these finding, together with results from analyses of nonsynonymous to synonymous rates of change (dN/dS) suggest that snake SSPs have endured substantial accelerated adaptive protein evolution. Such accelerated positive selection in SSPs parallels other findings of similar molecular evolution in snake venom proteins and suggests that diversifying selection on both systems may be linked, and that snake SSP genes may have evolved by gene duplication and rapid diversification to facilitate the acquisition of various functions to block venom activity within venomous snakes.     

Rapid evolution of immunoglobulin superfamily C2 domains expressed in immune system cells

To test the hypothesis that proteins expressed in cells of the vertebrate immune system evolve unusually rapidly, 107 orthologous immunoglobulin C2 domains were compared between human and murine rodent. The analysis showed that the rate of nonsynonymous (amino-acid- altering) nucleotide substitution in these domains was correlated with factors associated with protein structure and with breadth of tissue expression, as well as with the rate of synonymous substitution. However, when such factors were controlled for statistically, there remained a strong positive association between expression in the immune system and nonsynonymous rate, with the highest rates being seen in genes expressed in the immune system only. Certain immune system genes are known to be subject to positive selection favoring diversity at the amino acid level; most of these genes encode receptors that interact directly with foreign antigens. The observed acceleration of the rate of nonsynonymous evolution in C2 domains of immune system proteins may be explained by either (1) reduced constraint at the amino acid level on molecules interacting with immune system receptors that are themselves evolving rapidly due to positive diversifying selection or (2) positive selection favoring amino acid changes correlated with changes in the immune system receptors.      

Bifunctional phosphoglucose/phosphomannose isomerases from the Archaea Aeropyrum pernix and Thermoplasma acidophilum constitute a novel enzyme family within the phosphoglucose isomerase superfamily

The hyperthermophilic crenarchaeon Aeropyrum pernix contains phosphoglucose isomerase (PGI) activity. However, obvious homologs with significant identity to known PGIs could not be identified in the sequenced genome of this organism. The PGI activity from A. pernix was purified and characterized. Kinetic analysis revealed that, unlike all known PGIs, the enzyme catalyzed reversible isomerization not only of glucose 6-phosphate but also of epimeric mannose 6-phosphate at similar catalytic efficiency, thus defining the protein as bifunctional phosphoglucose/phosphomannose isomerase (PGI/PMI). The gene pgi/pmi encoding PGI/PMI (open reading frame APE0768) was identified by matrix-assisted laser desorption ionization time-of-flight analyses; the gene was overexpressed in Escherichia coli as functional PGI/PMI. Putative PGI/PMI homologs were identified in several (hyper)thermophilic archaea and two bacteria. The homolog from Thermoplasma acidophilum (Ta1419) was overexpressed in E. coli, and the recombinant enzyme was characterized as bifunctional PGI/PMI. PGI/PMIs showed low sequence identity to the PGI superfamily and formed a distinct phylogenetic cluster. However, secondary structure predictions and the presence of several conserved amino acids potentially involved in catalysis indicate some structural and functional similarity to the PGI superfamily. Thus, we propose that bifunctional PGI/PMI constitutes a novel protein family within the PGI superfamily.     

A new method for estimating nonsynonymous substitutions and its applications to detecting positive selection

The standard methods for computing the number of nonsynonymous substitutions (Ka) lump all amino acid changes into one single class, even though their rates of substitution vary by at least 10-fold (Tang et al., 2004). Classifying these changes by their physicochemical properties has not been suitably effective in isolating the fastest evolving classes of changes. We now propose to use the Universal index U of Tang et al. (2004) to classify the 75 elementary amino acid changes (codons differing by 1 bp) by their evolutionary exchangeability. Let Ki denote the Ka value of each class (i = 1, ..., 75 from the most to the least exchangeable). The cumulative Ki for the top 10 classes, denoted Kh (for high-exchangeability types), has two important properties: (1) Kh usually accounts for 25%-30% of total amino acid changes and (2) when the observed number of amino acid substitutions is large, Kh is predictably twice the value of Ka. This shall be referred to as the twofold approximation. The new method for estimating Kh is applied to the comparisons between human and macaque and between mouse and rat. The twofold approximation holds well in these data sets, and the signature of positive selection can be more easily discerned using the Kh statistic than using Ka. Many genes with Ka/Ks > 0.5 can now be shown to have Kh/Ks > 1 and to have evolved adaptively, at least for the high-exchangeability group of amino acid changes.     

Positive selection signatures in the TLR7 family

While adaptive immunity genes evolve rapidly under the influence of positive selection, innate immune system genes are known to evolve slowly due to strong purifying selection. Among the sensors of the innate immune system, Toll-like receptors (TLRs) are particularly important due to their ability to recognize and respond to pathogen-associated molecular patterns (PAMP), such as lipopolysaccharides, peptidoglycans, and nucleic acids from bacteria or viruses. In the present study, we examine the evolutionary process that has operated on the TLR7 family genes TLR7, TLR8, and TLR9. The results demonstrate that the average Ka/Ks (the ratio between nonsynonymous and synonymous substitution rates) of each TLR family gene is far lower than one regardless of estimating methods, supporting previous observations of strong purifying selection in this gene family. Interestingly, however, analysis of Ka/Ks ratios along the coding regions of TLR7 family genes by sliding-window analysis reveals a few narrow high peaks (Ka/Ks > 1). The most prominent peak corresponds to a specific region in the ectodomain, which exists only in the TLR7 family, suggesting that this unique structure of the TLR7 family might have been a target of positive selection in a variety of lineages. Furthermore, maximum likelihood model tests suggest that positive selection is the best explanation for a certain fraction of the amino acid substitutions in the TLR9.     

A novel phosphoglucose isomerase (PGI)/phosphomannose isomerase from the crenarchaeon Pyrobaculum aerophilum is a member of the PGI superfamily: structural evidence at 1.16-A resolution

The crystal structure of a dual specificity phosphoglucose isomerase (PGI)/phosphomannose isomerase from Pyrobaculum aerophilum (PaPGI/PMI) has been determined in native form at 1.16-A resolution and in complex with the enzyme inhibitor 5-phosphoarabinonate at 1.45-A resolution. The similarity of its fold, with the inner core structure of PGIs from eubacterial and eukaryotic sources, confirms this enzyme as a member of the PGI superfamily. The almost total conservation of amino acids in the active site, including the glutamate base catalyst, shows that PaPGI/PMI uses the same catalytic mechanisms for both ring opening and isomerization for the interconversion of glucose 6-phosphate (Glc-6-P) to fructose 6-phosphate (Fru-6-P). The lack of structural differences between native and inhibitor-bound enzymes suggests this activity occurs without any of the conformational changes that are the hallmark of the well characterized PGI family. The lack of a suitable second base in the active site of PaPGI/PMI argues against a PMI mechanism involving a trans-enediol intermediate. Instead, PMI activity may be the result of additional space in the active site imparted by a threonine, in place of a glutamine in other PGI enzymes, which could permit rotation of the C-2-C-3 bond of mannose 6-phosphate.     

Within-host evolution of CD8+-TL epitopes encoded by overlapping and non-overlapping reading frames of simian immunodeficiency virus

In order to understand the impact of overlapping reading frames on natural selection by host CD8+ T lymphocytes (CD8(+)-TL), we analyzed the pattern of nucleotide substitution in simian immunodeficiency virus (SIV) genomes sampled from populations at time of death in 35 rhesus monkeys. Both the mean number of nonsynonymous nucleotide substitutions per nonsynonymous site (d(N)) and the mean number of synonymous nucleotide substitutions per synonymous site (d(S)) were elevated in overlap regions in comparison to non-overlap regions. Mean d(N) exceeded mean d(S) in CD8(+)-TL epitopes restricted by the host's class I major histocompatibility complex molecules. This pattern, which is indicative of positive Darwinian selection favoring amino acid changes in these epitopes, was seen in both overlap and non-overlap regions; but mean d(N) was particularly elevated in restricted CD8(+)-TL epitopes encoded in overlap regions. Amino acid changes from the inoculum were defined as parallel if the same amino acid change occurred at the same site independently in two or more monkeys, and a surprisingly high proportion (71.9%) of observed amino acid changes throughout the SIV genome occurred in parallel in different monkeys. The proportion of parallel changes in restricted epitopes encoded by overlapping reading frames was still higher (80%), supporting the hypothesis that the interaction of positive selection and overlapping reading frames enhances the probability of convergent or parallel amino acid change.     

Evolution of prothrombin: Isolation and characterization of the cDNAs encoding chicken and hagfish prothrombin

The cDNA sequences of chicken and hagfish prothrombin have been determined. The sequences predict that prothrombin from both species is synthesized as a prepro-protein consisting of a putative Gla domain, two kringle domains, and a two-chain protease domain. Chicken and hagfish prothrombin share 51.6% amino acid sequence identity (313/627 residues). Both chicken and hagfish prothrombin are structurally very similar to human, bovine, rat, and mouse prothrombin and all six species share 41% amino acid sequence identity. Amino acid sequence alignments of human, bovine, rat, mouse, chicken, and hagfish prothrombin suggest that the thrombin B-chain and the propeptide-Gla domain are the regions most constrained for the common function(s) of vertebrate prothrombins.The nucleotide sequences reported in this paper have been submitted to the EMBL/Genbank database under the following secession numbers: M 81391 for Gallus gallus, M 81393 for Eptatretus stouti.Correspondence to: R.T.A. MacGillivray     

Adaptations to in situ feeding: novel nutrient acquisition pathways in an ancient vertebrate

During feeding, hagfish may immerse themselves in the body cavities of decaying carcasses, encountering high levels of dissolved organic nutrients. We hypothesized that this feeding environment might promote nutrient acquisition by the branchial and epidermal epithelia. The potential for Pacific hagfish, Eptatretus stoutii, to absorb amino acids from the environment across the skin and gill was thus investigated. l-alanine and glycine were absorbed via specific transport pathways across both gill and skin surfaces, the first such documentation of direct organic nutrient acquisition in a vertebrate animal. Uptake occurred via distinct mechanisms with respect to concentration dependence, sodium dependence and effects of putative transport inhibitors across each epithelium. Significant differences in the absorbed amino acid distribution between the skin of juveniles and adults were noted. The ability to absorb dissolved organic matter across the skin and gill may be an adaptation to a scavenging lifestyle, allowing hagfish to maximize sporadic opportunities for organic nutrient acquisition. From an evolutionary perspective, hagfish represent a transitory state between the generalized nutrient absorption pathways of aquatic invertebrates and the more specialized digestive systems of aquatic vertebrates.     

The complete nucleotide sequences of the 5 genetically distinct plastid genomes of Oenothera, subsection Oenothera: II. A microevolutionary view using bioinformatics and formal genetic data

A unique combination of genetic features and a rich stock of information make the flowering plant genus Oenothera an appealing model to explore the molecular basis of speciation processes including nucleus-organelle coevolution. From representative species, we have recently reported complete nucleotide sequences of the 5 basic and genetically distinguishable plastid chromosomes of subsection Oenothera (I-V). In nature, Oenothera plastid genomes are associated with 6 distinct, either homozygous or heterozygous, diploid nuclear genotypes of the 3 basic genomes A, B, or C. Artificially produced plastome-genome combinations that do not occur naturally often display interspecific plastome-genome incompatibility (PGI). In this study, we compare formal genetic data available from all 30 plastome-genome combinations with sequence differences between the plastomes to uncover potential determinants for interspecific PGI. Consistent with an active role in speciation, a remarkable number of genes have high Ka/Ks ratios. Different from the Solanacean cybrid model Atropa/tobacco, RNA editing seems not to be relevant for PGIs in Oenothera. However, predominantly sequence polymorphisms in intergenic segments are proposed as possible sources for PGI. A single locus, the bidirectional promoter region between psbB and clpP, is suggested to contribute to compartmental PGI in the interspecific AB hybrid containing plastome I (AB-I), consistent with its perturbed photosystem II activity.     

Degradation, receptor binding affinity, and potency of insulin from the Atlantic hagfish (Myxine glutinosa) determined in isolated rat fat cells

Insulin from the Atlantic hagfish, Myxine glutinosa, a primitive vertebrate, was studied with respect to degradation, receptor binding, and stimulation of glucose transport and metabolism in isolated rat adipocytes. The degradation was studied in a concentrated suspension with about 100mul of cells/ml of suspension. 125I-labeled hagfish insulin and 125I-labeled pig insulin were degraded at the same rate when present in concentrations of 0.3nM. Native hagfish insulin inhibited the rate of degradation of 125I-labeled pig insulin half-maximally at a concentration of 12+/-2 nM (S.D., n=6) as compared to 130+/-32 nM (S.D.,n=6) for pig insulin. Native hagfish insulin in a concentration of 130 nM was biologically inactivated at a rate several times slower than pig insulin in the same concentration. The results indicate that the maximal velocity (Vmax) of degradation of hagfish insulin as well as the concentration causing half-maximal velocity (Km) are about 10 times lower for hagfish insulin than for pig insulin. The receptor binding and the biological effects of hagfish insulin were studied in dilute cell suspensions where the degradation of hormone in the medium was negligible. The receptor binding affinity of hagfish insulin was 23+/-7 per cent (S.D., n=10) of that of pig insulin. Hagfish insulin was able to elicit the same maximal stimulation of both 3-o-methylglucose exchange and lipogenesis from glucose as pig insulin. However, the potency of hagfish insulin with respect to activation of lipogenesis was only 4.6+/-0.6 per cent (S.D., n=15) of that of pig insulin. Hagfish insulin thus constitutes the first described insulin which exhibits a discrepancy between relative binding affinity and relative potency. This discrepancy was not due to the methionine residue (B31) at the COOH-terminal end of the B chain of hagfish insulin, since removal of this residue caused no marked change in the binding affinity or the potency. The results indicate that the receptor occupancy must be 5 times higher with hagfish insulin than with pig insulin to cause a particular degree of activation of lipogenesis. Hagfish insulin might therefore be characterized as a "partial antagonist" on the receptors. However, it was not possible to demonstrate antagonistic properties of hagfish insulin on the cells. The effect of hagfish insulin plus pig insulin in submaximally stimulating concentrations was additive. Furthermore, the decay of activation of adipocytes after incubation with hagfish insulin followed the same time course as the decay of activation after incubation with pig insulin in a concentration of equal potency. These phenomena are in agreement with the concept that adipocytes possess a large excess of receptors which can mediate the effect of insulin on lipogenesis from glucose.     

SRY基因在人猿超科和旧大陆猴中具有不同的进化规律

通过ＰＣＲ扩增、测序,得到了白臀叶猴和红面猴的ＳＲＹ基因全序列。结合现有的灵长类其他物种序列进行分析,验证了ＨＭＧ盒的保守性。通过构建系统发育树,比较旧大陆猴和人猿超科两个类群内和类群间ＨＭＧ盒侧翼序列Ｋａ／Ｋｓ的比率。有趣的是,人猿超科两物种比较呈现较高的Ｋａ／Ｋｓ比值,但在旧大陆猴中及旧大陆猴与狨猴间的Ｋａ／Ｋｓ比值显著低于人猿超科的,呈现很不同的格局。同时,对于ＨＭＧ盒序列,Ｋａ／Ｋｓ比值在