Evolution of the amino acid substitution in the mammalian myoglobin gene

Summary Multivariate statistical analyses were applied to 16 physical and chemical properties of amino acids. Four of these properties; volume, polarity, isoelectric point (charge), and hydrophobicity were found to explain adequately 96% of the total variance of amino acid attributes. Using these four quantitative measures of amino acid properties, a structural discriminate function in the form of a weighted difference sum of squares equation was developed. The discriminate function is weighted by the location of each particular residue within a given tertiary structure and yields a numerical discriminate or difference value for the replacement of these residues by different amino acids. This resulting discriminate value represents an expression of the perturbation in the local positional environment of a protein when an amino acid substitution occurs. With the use of this structural discriminate function, a residue by residue comparison of the known mammalian myoglobin sequences was carried out in an attempt to elucidate the positions of possible deviations from the known tertiary structure of sperm whale myoglobin. Only 11 of the 153 residue positions in myoglobin demonstrated possible structural deviations. From this analysis, indices of difference were calculated for all amino acid exchanges between the various myoglobins. All comparisons yielded indices of difference that were considerably lower than would be expected if mutations had been fixed at random, even if the organization of the genetic code is taken into consideration. On the basis of these results, it is inferred that some form of selection has acted in the evolution of mammalian myoglobins to favor amino acid substitutions that are compatible with the retention of the original conformation of the protein.     

Comparative analysis of the spatial organization of myoglobins. I. Hydrophobicity profiles]

Hydrophobicity profiles of myoglobins in the animal species far remote in the evolutionary series are considerably similar. A complete coincidence as to the arrangement of hydrophobic zones along the polypeptide chain in myoglobins of the compared species (from a man to mollusc) is revealed at the beginning of alpha-helix of B-segment and in the area corresponding to a cluster which embodies a heme- bound water molecule, distal histidine E7 being directed to this cluster. The mollusc myoglobin with two absent (as compared to myoglobins of other species) hydrophobic sites differs in the profile of hydrophobicity most of all. It is supposed that hydrophobic nuclei forming the heme circumference create a globule "skeleton" thus pre-setting general spatial structure of the myoglobin molecule, which is very significant for its functional activity.     

The antibody response to myoglobin is independent of the immunized species. Analysis in terms of replacements in the antigenic sites and in environmental residues of the cross-reactions of fifteen myoglobins with sperm-whale myoglobin antisera raised in different species.

The recent determination of the entire antigenic structure of sperm-whale myoglobin with rabbit and goat antisera has permitted the examination of whether the antigenic structure recognized by antibodies depends on the species in which the antisera are raised. Also, by knowledge of the antigenic structure, the molecular factors that determine and influence antigenicity can be better understood in terms of the effects of amino acid substitutions occurring in the antigenic sites and in the environmental residues of the sites. In the present work, the myoglobins from finback whale, killer whale, horse, chimpanzee, sheep, goat, bovine, echidna, viscacha, rabbit, dog, cape fox, mouse and chicken were examined for their ability to cross-react with antisera to sperm-whale myoglobin. By immunoadsorbent titration studies with radioiodinated antibodies, each of these myoglobins was able to bind antibodies to sperm-whale myoglobin raised in goat, rabbit, chicken, cat, pig and outbred mouse. It was found that the extent of cross-reaction of a given myoglobin was not dependent on the species in which the antisera were raised. This indicated that the antibody response to sperm-whale myoglobin (i.e. its antigenic structure) is independent of the species in which the antisera are raised and is not directed to regions of sequence differences between the injected myoglobin and the myoglobin of the immunized host. Indeed, in each antiserum from a given species examined, that antiserum reacted with the myoglobin of that species. The extent of this auto-reactivity for a given myoglobin was comparable with the general extent of cross-reactivity shown by that myoglobin with antisera raised in other species. The cross-reactivities and auto-reactivities (both of which are of similar extents for a given myoglobin) can be reasonably rationalized in terms of the effects of amino acid substitutions within the antigenic sites and within the residues close to these sites. These findings confirm that the antigenicity of the sites is inherent in their three-dimensional locations.     

Primary structure of otter (Lutra lutra L.) myoglobin. II. Pepsin peptides of trypsin hydrolysate. Reconstruction of the polypeptide chain of the otter myoglobin globin component

13 peptic peptides have been isolated from the insoluble (at pH 5.0) fraction of the tryptic hydrolysate of main chromatographic component of otter myoglobin and their amino acid composition and N-terminal amino acid sequences have been determined. The isolated peptides contain in total 40 amino acid residues. The results obtained, along with those on tryptic peptides and the comparison with homologous portions of myoglobins of the known primary structure, allowed reconstructing the complete amino acid sequence of otter myoglobin.     

Amino acid sequence of Ondatra myoglobin (Ondatra zibethica). Characteristic features of myoglobins from semiaquatic animals

Based on the amino acid composition of globin, amino acid analysis and N-terminal sequencing of peptides as well as a comparative analysis of the primary structure of beaver, coypu rat and otter myoglobins with the use of the fingerprinting technique, the amino acid sequence of the major component of ondatra myoglobin including 153 amino acid residues was reconstructed. The results of a comparative analysis of the primary structure of myoglobin and the peculiarities of the functional morphology of myoglobins from semi-aquatic animals and sperm whale and the role of amino acid substitutions in the spatial structure of ondatra myoglobin are discussed.     

Primary structure of myoglobins from 31 species of birds

Primary structure of myoglobins (Mbs) from 31 avian species of 15 orders were reported, although portions of the structures in the 2 species could not be determined. At least 68 of the total 153 amino acid sites were invariant all through the avian, reptilian and human Mbs, and 20 of these sites were "internal", forming the internal hydrophobic cavities in which the heme group remains wrapped. Furthermore, at 27 sites, if replaced, the replacements were mostly conservative, and 13 of the conservative sites were "internal". Thus the all 33 "internal" sites, important for structural and functional stability of the protein, have been well preserved, either invariant or conserved, during evolution from reptiles to birds and mammals. The residue 71 (E14) in 4 penguin species was not deleted as previously reported in emperor penguin Mb but occupied by Gln. The residue 121 (GH3) was deleted in all 3 species studied of Falconiformes. Out of 9 anseriforms, 5 species of different genera showed the identical structure. Secondary structures as viewed by hydropathy profiles were highly similar throughout the reptilian, avian and mammalian Mbs.     

Amino acid sequence of myoglobin from the molluscBursatella leachii

The complete amino acid sequence of myoglobin from the triturative stomach of gastropodic molluscBursatella leachii has been determined. It is composed of 146 amino acid residues, is acetylated at the N-terminus, and contains a single histidine residue at position 95 which corresponds to the heme-binding proximal histidine. The E7 distal histidine, which is conserved widely in myoglobins and hemoglobins, is replaced by valine inBursatella myoglobin. The amino acid sequence ofBursatella myoglobin shows strong homology (73–84%) with those ofAplysia andDolabella myoglobins.     

Myoglobin of the shark Heterodontus portusjacksoni: isolation and amino acid sequence

Myoglobin isolated from red muscle of the shark H. portusjacksoni was purified by ion-exchange chromatography on sulfopropyl-Sephadex and gel-filtration. Amino acid analysis and sequence determination showed 148 amino acid residues. The amino terminal residue is acetylated as shown by mass spectrographic analysis of N-terminal peptides. There is a deletion of four residues at the amino terminal end as well as one residue in the CD interhelical area relative to other myoglobins. The complete amino acid sequence has been determined following digestion with trypsin, chymotrypsin, pepsin and staphylococcal protease. Sequences of the purified peptides were determined by the dansyl-Edman procedure. The amino acid sequence showed approximately 85 differences from mammalian, monotreme and bird myoglobins. The date of divergence of the shark H. portusjacksoni from these other orders was estimated at 450 +/- 16 million years, based on the number of amino acid differences between species and allowing for multiple mutations during the evolutionary period. This estimate agrees well with similar estimates made using alpha- and beta-globin sequences, in contrast to widely differing estimates of dates of divergence for monotremes using the same three globin chains. Compared with myoglobins from species previously studied, there are many more differences in amino acid sequences, and in many positions residues are found that are more characteristic of alpha- and beta-globins, suggesting a conservation of residues over a long period of evolutionary time. There are fewer stabilizing hydrogen bonds and salt-linkages than in other myoglobins.     

Structural and functional aspects of the heart ventricle myoglobin of bluefin tuna

The heart ventricle myoglobin of bluefin tuna has been purified to an apparent homogeneity. The amino acid analysis has revealed only a limited number of substitutions between the myoglobins of yellowfin and bluefin tuna. The alpha-helix content of tuna myoglobin has been found considerably lower than that of mammalian myoglobin. No correlation has been discovered between the conformational stability and alpha-helix content. Denaturation experiments have shown that the whole structure of tuna myoglobin results from the interaction of two structural units which represent the product of independent folding processes. The structure of tuna myoglobin has been found more open and disorganized than that of sperm whale. This result has been related to the low content of electrostatic interactions and explained in terms of evolutive adaptations.     

The conformation of T4 bacteriophage dihydrofolate reductase from circular dichroism

The secondary and tertiary structure of T4 bacteriophage dihydrofolate reductase is investigated by vacuum ultraviolet circular dichroism (CD) spectroscopy and probability analysis of the primary amino acid sequence. The far ultraviolet CD spectrum of the enzyme in the range of 260-178 nm is analyzed by the generalized inverse and variable selection methods developed by our laboratory. Variable selection yields an average content of 26% alpha-helix, 21% antiparallel beta-sheet, 10% parallel beta-sheet, 20% beta-turns, and 32% "other" structures within the T4 protein. The characteristic peaks of the CD spectrum indicate that the enzyme has a lot of antiparallel beta-sheet, which is typical of the alpha + beta tertiary class of globular proteins. The secondary structure of the protein is also analyzed by using four statistical methods on the amino acid sequence. Although the secondary structures predicted by each individual statistical method vary to a considerable extent, the fractions of each structure jointly predicted by a majority of the methods are in excellent agreement with our CD analysis. The alternating arrangement for some segments of alpha-helix and beta-sheet predicted from primary structure to be within the enzyme is characteristic of proteins containing parallel beta-sheet. This supports our conclusion that the protein contains both parallel and antiparallel beta-sheet structures, but finding both types of beta-sheet also means that the protein may have the variation on alpha/beta tertiary structure recently found in EcoRI endonuclease and thymidylate synthase. These observations, in conjunction with other physical properties of the T4 reductase, suggest that the enzyme perhaps shares an evolution in common with the dihydrofolate reductases derived from type I R-plasmids rather than with the host-cell protein.     

A method for detecting centres of natural selection in protein structures: potential for predicting the location of functional areas

A simple method has been developed to detect protein microenvironments which are likely to be the focus of natural selection, and thereby important for function. It relies on distinguishing between the probability of an amino acid type arising by genetic mutation and the probability that it will be chosen by natural selection. When applied to proteins of known tertiary structure, the method revealed that major differences exist in the balance between neutral and selected change, and also that functional sites can be highlighted.     

Nearest-neighbour analysis of myoglobin antigenic sites. Nearest-neighbour residues whose replacement can alter the environment of binding-site residue(s) and thus change their characteristics and binding capability.

By using the known antigenic structure of sperm-whale myoglobin previously determined in this laboratory and the X-ray co-ordinates for the myoglobin molecule, we have calculated the nearest-atom distances between each of the residues of the antigenic sites and all the other amino acids of the myoglobin molecule. These calculations have enabled us to identify the nearest-neighbour residues to each of the residues in the five antigenic sites, and which thus describe the immediate molecular environment of the sites. The influences of chemical changes or replacements in these environmental residues on the binding capacity of an antigenic site, when considered together with replacements directly in the antigenic sites, are expected to account for the major effects and will be extremely useful in explaining the cross-reactions of myoglobins from various species. However, it is stressed that the analysis has limitations due to the qualitative estimates of the effects, the influences of substitutions of once-removed or even at more distant locations (especially when they are cumulative) and finally the influences of any conformational re-adjustments when these occur as a result of the replacement(s).     

3-D Model of the bee venom acid phosphatase: Insights into allergenicity

The acid phosphatase Api m 3 is the major allergen of the honeybee venom. Except for the amino acid sequence, no other structural information for the enzyme is available. We applied homology modeling to assign the three-dimensional structure of Api m 3. The structure of the homodimeric human prostatic acid phosphatase was used to model the Api m 3 tertiary structure. IgE epitopes and antigenic sites were predicted using programs based on the structure of known epitopes and analysis of the 3-D model. The model of Api m 3 revealed an active site similar to those of the histidine—type acid phosphatases with conservation of the catalytically important residues. The observed substitutions in the phosphate ion binding site suggest differences in the substrate specificity in comparison to other acid phosphatases. The analysis of the Api m 3 three-dimensional model revealed a very likely mechanism of enzyme action.     

Abalone myoglobins evolved from indoleamine dioxygenase: The cDNA-derived amino acid sequence of myoglobin fromNordotis madaka

The cDNA for the unusual 41 kD myoglobin of the abaloneNordotis madaka was amplified by polymerase chain reaction (PCR), and the cDNA-derived amino acid sequence of 378 residues was determined. As with the myoglobin of the related abaloneSulculus diversicolor (Suzuki and Takagi,J. Mol. Biol. 228, 698–700, 1992), the sequence ofNordotis myoglobin showed no significant homology with any other globins, but showed high homology (35% identity) with vertebrate indoleamine 2,3-dioxygenase, a tryptophan degrading enzyme containing heme. The amino acid sequence homology betweenNordotis andSulculus myoglobins was 87%. These results support our previous idea that the abalone myoglobins evolved from a gene for indoleamine dioxygenase, but not from a globin gene, and therefore all of the hemoglobins and myoglobins are not homologous. Thus, abalone myoglobins appear to be a typical case of convergent evolution.     

Missense mutation in the lacI gene of Escherichia coli. Inferences on the structure of the repressor protein

The lac repressor has been studied extensively but a precise three-dimensional structure remains unknown. Studies using mutational data can complement other information and provide insight into protein structure. We have been using the lacI gene-repressor protein system to study the mutational specificity of spontaneous and induced mutation. The sequencing of over 6000 lacI- mutations has revealed 193 missense mutations generating 189 amino acid replacements at 102 different sites within the lac repressor. Replacement sites are not distributed evenly throughout the protein, but are clustered in defined regions. Almost 40% of all sites and over one-half of all substitutions found occur within the amino-terminal 59 amino acid residues, which constitute the DNA-binding domain. The core domain (residues 60 to 360) is less sensitive to amino acid replacement. Here, substitution is found in regions involved in subunit aggregation and at sites surrounding residues that are implicated in sugar-binding. The distribution and nature of missense mutational sites directs attention to particular amino acid residues and residue stretches.     

The stepwise mutation model: an experimental evaluation utilizing hemoglobin variants

The stepwise mutation model of Ohta and Kimura (1973) was proposed to explain patterns of genetic variability revealed by means of electrophoresis. The assumption that electrophoretic mobility was principally determined by unit changes in net molecular charge has been criticized by Johnson (1974, 1977). This assumption has been tested directly using hemoglobin. Twenty-seven human hemoglobin variants with known amino acid substitutions, and 26 nonhuman hemoglobins with known sequences were studied by starch gel electrophoresis. Of these hemoglobins, 60 to 70% had electrophoretic mobilities that could be predicted solely on the basis of net charge calculated from the amino acid composition alone, ignoring tertiary structure. Only four hemoglobins showed a mobility that was clearly different from an expected mobility calculated using only the net charge of the molecule. For the remaining 30% of hemoglobins studied, mobility was determined by a combination of net charge and other unidentified components, probably reflecting changes in ionization of some amino acid residues as a result of small alterations in tertiary structure due to the amino acid substitution in the variant. For the nonhuman hemoglobins, the deviation of a sample from its expected mobility increased with increasing amino acid divergence from human hemoglobin A.-It is concluded that the net electrostatic charge of a molecule is the principal determinant of electrophoretic mobility under the conditions studied. However, because of the significant deviation from strict stepwise mobility detected for 30 to 40% of the variants studied, it is further concluded that the infinite-allele model of Kimura and Crow (1964) or a "mixed model" such as that proposed by Li (1976) may be more appropriate than the stepwise mutation model for the analysis of much of the available electrophoretic data from natural populations.     

Comparison of dissimilarity patterns of E coli, yeast and mammalian tRNAs.

A number of experimental approaches have been developed for identification of recognition (identity) sites in tRNAs. Along with them a theoretical methodology has been proposed by McClain et al that is based on concomitant analysis of all tRNA sequences from a given species. This approach allows an evaluation of nucleotide combinations present in isoacceptor tRNAs specific for the given amino acid, and not present in equivalent positions in cloverleaf structure in other tRNAs of the same organism. These elements predicted from computer analysis of the databank could be tested experimentally for their participation in forming recognition sites. The correlation between theoretical predictions and experimental data appeared promising. The aim of the present work consisted of introducing further improvements into McClain's procedure by: i), introducing into analysis a variable region in tRNAs which had not been previously considered; to accomplish this, 'normalization' of variable nucleotides was suggested, based on primary and tertiary structures of tRNAs; ii), developing a new procedure for comparison of patterns for synonymous and non-synonymous tRNAs from different organisms; iii), analysis of 3- and 4-positional contacts between tRNAs and enzymes in addition to a formerly used 2-positional model. A systematic application of McClain's procedure to mammalian, yeast and E coli tRNAs led to the following results: i), imitancy patterns for non-synonymous tRNAs of any amino acid specificity and from any organisms analysed so far overlap by no more than 30%, providing a structural basis for discrimination with high fidelity between cognate and non-cognate tRNAs; ii), the predicted identity sites are non-randomly distributed within tRNA molecules; the dominant role is ascribed to only two regions--anticodon and amino acid stem which are located far apart from one another at extremes of all tRNA molecules; iii), the imitancy patterns for synonymous tRNAs in lower (yeast) and higher (mammalian) eukaryotes are similar but not identical; iv), distribution of predicted identity sites in the cloverleaf structure in prokaryotes and eukaryotes is essentially different: in eubacterial tRNAs the major role in recognition plays anticodon and/or amino acid acceptor stem, whereas in eukaryotic (both unicellular and multicellular) tRNAs the remaining part of the molecules is also involved in recognition; v), the imitancy patterns of synonymous tRNAs from prokaryotes and eukaryotes are dissimilar, this observation leads to the prediction that the tRNA identity sites for the same amino acid in prokaryotes and eukaryotes may differ.     

NMR study of Galeorhinus japonicus myoglobin. 1H-NMR study of molecular structure of the heme cavity

The molecular structure of the active site of myoglobin from the shark, Galeorhinus japonicus, has been studied by 1H-NMR. Some hyperfine-shifted amino acid proton resonances in the met-cyano form of G. japonicus myoglobin have been unambiguously assigned by the combined use of various two-dimensional NMR techniques; they were compared with the corresponding resonances in Physter catodon myoglobin. The orientations of ThrE10 and IleFG5 residues relative to the heme in G. japonicus met-cyano myoglobin were semiquantitatively estimated from the analysis of their shifts using the magnetic susceptibility tensor determined by a method called MATDUHM (magnetic anisotropy tensor determination utilizing heme methyls) [Yamamoto, Y., Nanai, N. & Ch?j?, R. (1990) J. Chem. Soc., Chem. Commun., 1556-1557] and the results were compared with the crystal structure of P. catodon carbonmonoxy myoglobin [Hanson, J. C. & Schoenborn, B. P. (1981) J. Mol. Biol. 153, 117-124]. In spite of a substantial difference in shift between the corresponding amino acid proton resonances for the two proteins, the orientations of these amino acid residues relative to the heme in the active site of both myoglobins were found to be highly alike.     

Amino acid sequence of myoglobin from the molluscBursatella leachii

The complete amino acid sequence of myoglobin from the triturative stomach of gastropodic molluscBursatella leachii has been determined. It is composed of 146 amino acid residues, is acetylated at the N-terminus, and contains a single histidine residue at position 95 which corresponds to the heme-binding proximal histidine. The E7 distal histidine, which is conserved widely in myoglobins and hemoglobins, is replaced by valine inBursatella myoglobin. The amino acid sequence ofBursatella myoglobin shows strong homology (73–84%) with those ofAplysia andDolabella myoglobins.     

棉铃虫多核型多角体病毒38k基因的克隆、序列分析及其蛋白高级结构的预测

采用半补齐方法建立棉铃虫多核衣壳型多角体病毒基因组文库,通过对插入片段进行克隆鉴定和序列分析,获得了38k基因.该基因上游具有晚期调控保守序列TTAAG,是一个晚期表达基因,基因阅读框为903 bp,共编码300个氨基酸,氨基酸序列同源性分析结果表明其与α类杆状病毒的同源性较高,有较近的亲缘关系.氨基酸高级结构的分析表明其与与磷酸酶结构相似性达到95%,与病毒核衣壳的组装有关.