Location of amino acid substitutions in human hemoglobin. Mass spectrometric rapid analysis of tryptic peptides

Point mutations alpha 58 His----Tyr (Hb M Boston), beta 6 Glu--lys (Hb C) and beta 26 Glu----Lys (Hb E) have been identified in abnormal hemoglobins by means of tryptic hydrolysis of their alpha- and beta-chains followed by mass-spectrometry coupled with direct extraction of ions from solution. The abnormal hemoglobin Hb M Boston alpha 58 (E7) His----Tyr has been for the first time found in the blood of a patient from the USSR. This express-method is generally applicable for the identification of point mutation in proteins. The amount of protein necessary for the analysis is 100-1000 pmole. The stability, proteolytic degradation of the identified abnormal Hb's and Hb Bart's were investigated. The molecular pathogenesis of the hemoglobinopathies are discussed from the point of view of the observed properties.     

Laser kinetic spectroscopy studies of the bimolecular oxygenation of alpha- and beta-subunits within the R-state of human hemoglobin

Bimolecular oxygenation of tri-liganded R-state human hemoglobin (HbA) is described by bi-exponential kinetics with association rate constants k = 27.2 ± 1.3 (M·sec)^-1 and k = 62.9 ± 1.6 (M·sec)^-1. Both the observed processes have been assigned to the bimolecular oxygenation of - and -subunits of the native tetrameric protein by molecular oxygen. The quantum yields of photodissociation within the completely oxygenated R-state HbA are = 0.0120 ± 0.0017 and = 0.044 ± 0.005 for - and -subunits, respectively. The oxygenation reactions of isolated ^PCMB- and ^PCMB-hemoglobin chains are described by mono-exponential kinetics with the association rate constants k = 44 ± 2 (M·sec)^-1 and k = 51 ± 1 (M·sec)^-1, respectively. The quantum yields of photodissociation of isolated ^PCMB- and ^PCMB-chains (0.056 ± 0.006 and 0.065 ± 0.006, respectively) are greater than that observed for appropriate subunits within the R-state of oxygenated HbA.     

Substitution of the heme binding module in hemoglobin alpha- and beta-subunits. Implication for different regulation mechanisms of the heme proximal structure between hemoglobin and myoglobin

In our previous work, we demonstrated that the replacement of the "heme binding module," a segment from F1 to G5 site, in myoglobin with that of hemoglobin alpha-subunit converted the heme proximal structure of myoglobin into the alpha-subunit type (Inaba, K., Ishimori, K. and Morishima, I. (1998) J. Mol. Biol. 283, 311-327). To further examine the structural regulation by the heme binding module in hemoglobin, we synthesized the betaalpha(HBM)-subunit, in which the heme binding module (HBM) of hemoglobin beta-subunit was replaced by that of hemoglobin alpha-subunit. Based on the gel chromatography, the betaalpha(HBM)-subunit was preferentially associated with the alpha-subunit to form a heterotetramer, alpha(2)[betaalpha(HBM)(2)], just as is native beta-subunit. Deoxy-alpha(2)[betaalpha(HBM)(2)] tetramer exhibited the hyperfine-shifted NMR resonance from the proximal histidyl N(delta)H proton and the resonance Raman band from the Fe-His vibrational mode at the same positions as native hemoglobin. Also, NMR spectra of carbonmonoxy and cyanomet alpha(2)[betaalpha(HBM)(2)] tetramer were quite similar to those of native hemoglobin. Consequently, the heme environmental structure of the betaalpha(HBM)-subunit in tetrameric alpha(2)[betaalpha(HBM)(2)] was similar to that of the beta-subunit in native tetrameric Hb A, and the structural conversion by the module substitution was not clear in the hemoglobin subunits. The contrastive structural effects of the module substitution on myoglobin and hemoglobin subunits strongly suggest different regulation mechanisms of the heme proximal structure between these two globins. Whereas the heme proximal structure of monomeric myoglobin is simply determined by the amino acid sequence of the heme binding module, that of tetrameric hemoglobin appears to be closely coupled to the subunit interactions.     

Binding modes of L35 to alpha- and beta-semihemoglobins: structural insights into the inequivalence of alpha- and beta-subunits of hemoglobin

It has been thought for several years that the greatly lowered oxygen affinity, high cooperativity, and heterotropic modulation displayed by tetrameric human hemoglobin (Hb) was an exclusive result of the assembly of high affinity alpha(1)beta(1) dimers into alpha(2)beta(2) tetramers. However, in recent times, it has been shown that alpha- and beta-semihemoglobins, namely alpha(heme)beta(apo) and alpha(apo)beta(heme), which are dimers of Hb characterized by a high affinity for oxygen and lack of cooperativity do respond to effectors such as 2-[4-(3,5-dichlorophenylureido) phenoxy]-2-methylpropionic acid (L35), a bezafibrate (BZF) related compound, by decreasing the ligand affinity to a considerable extent (between 60- and 130-fold). In order to shed some light on the structural basis of this phenomenon, we have developed a binding mode of L35 to semihemoglobins through docking analysis using the program GRID. Molecular modelling studies did identify sites on semihemoglobins where favourable interactions with L35 can occur. We found that the effector binds differently to the two semihemoglobins exhibiting high affinity only for the alpha chain heme pocket. The proposed binding models are consistent with the experimental findings and may be rationalized in terms of different hydrophobic and hydrophilic characteristics between alpha- and beta-heme pockets of Hb.     

The mutational spectrum of single base-pair substitutions causing human genetic disease: patterns and predictions

Summary Reports of single base-pair substitutions that cause human genetic disease and that have been located and characterized in an unbiased fashion were collated; 32% of point mutations were CG TG or CG CA transitions consistent with a chemical model of mutation via methylation-mediated deamination. This represents a 12-fold higher frequency than that predicted from random expectation, confirming that CG dinucleotides are indeed hotspots of mutation causing human genetic disease. However, since CG also appears hypermutable irrespective of methylation-mediated deamination, a second mechanism may also be involved in generating CG mutations. The spectrum of point mutations occurring outwith CG dinucleotides is also non-random, at both the mono- and dinucleotide levels. An intrinsic bias in clinical detection was excluded since frequencies of specific amino acid substitutions did not correlate with the chemical difference between the amino acids exchanged. Instead, a strong correlation was observed with the mutational spectrum predicted from the experimentally measured mispairing frequencies of vertebrate DNA polymerases and in vitro. This correlation appears to be independent of any difference in the efficiency of enzymatic proofreading/mismatch-repair mechanisms but is consistent with a physical model of mutation through nucleotide misincorporation as a result of transient misalignment of bases at the replication fork. This model is further supported by an observed correlation between dinucleotide mutability and stability, possibly because transient misalignment must be stabilized long enough for misincorporation to occur. Since point mutations in human genes causing genetic disease neither arise by random error nor are independent of their local sequence environment, predictive models may be considered. We present a computer model (MUTPRED) based upon empirical data; it is designed to predict the location of point mutations within gene coding regions causing human genetic disease. The mutational spectrum predicted for the human factor IX gene was shown to resemble closely the observed spectrum of point mutations causing haemophilia B. Further, the model was able to predict successfully the rank order of disease prevalence and/or mutation rates associated with various human autosomal dominant and sex-linked recessive conditions. Although still imperfect, this model nevertheless represents an initial attempt to relate the variable prevalence of human genetic disease to the mutability inherent in the nucleotide sequences of the underlying genes.     

Association of alpha- and beta-subunits during the biosynthesis of beta-hexosaminidase in cultured human fibroblasts

Subunit association of beta-hexosaminidase was studied in intact fibroblasts using antisera that discriminate between free and associated alpha-chains. These were anti-beta-hexosaminidase A (anti-alpha beta), which precipitated all alpha-chains, free or associated; anti-beta-hexosaminidase B (anti-beta beta), which precipitated those alpha-chains that were associated with beta; and anti-alpha-chains, which recognized only monomeric alpha-chains. After biosynthetic labeling, beta-hexosaminidase or its free alpha-subunit were immuno-precipitated from extracts of cells and medium with the aid of protein A-bearing Staphylococcus aureus, subjected to sodium dodecyl sulfate-polyacrylamide gel electrophoresis, and visualized by fluorography. Pulse-chase labeling showed that the alpha-chains existed predominantly in the monomeric precursor form during the first 5 h, and then began to accumulate in the mature (lysosomal) associated alpha beta form. Precursor alpha beta complexes were secreted, along with some precursor alpha monomers; the latter were catalytically inert. Both alpha- and beta-chains were phosphorylated (a Golgi modification) prior to association. Thus alpha-beta association probably occurred in the Golgi area before transfer to lysosomes and before secretion. Cycloheximide inhibited the association and subsequent maturation of preformed alpha-chains, perhaps by causing a depletion of a pool of beta-chain precursor upstream from the site of subunit association. In fibroblasts from a patient with Sandhoff disease, that produced no beta-chains, the alpha-chains self-associated but their maturation was markedly decreased. We suggest that association with beta-chains is necessary not only for acquisition of catalytic activity but also for transport of alpha-chains to lysosomes.     

A kinetic description of dioxygen motion within alpha- and beta-subunits of human hemoglobin in the R-state: geminate and bimolecular stages of the oxygenation reaction

Laser flash photolysis technique is used to study human hemoglobin (HbA) oxygenation. Monomolecular geminate oxygenation of triliganded R-state HbA molecules is described by a function of three exponentials. Geminate oxygenation of the alpha-subunit within R-state HbA is characterized by two components with time constants of 0.14 and 1 ns, while geminate oxygenation of the beta-subunit within HbA is characterized by two components with time constants of 1 and approximately 30 ns. Bimolecular oxygenation of triliganded R-state HbA molecules is described by a biexponential law. Two observed rate constants are assigned to oxygenation of the alpha- and beta-subunit within HbA. The bimolecular association rate constants for O(2) rebinding with the alpha- and beta-subunit within triliganded R-state HbA are k(alpha) = 18.8 +/- 1.3 (microM x s)(-1) and k(beta) = 52 +/- 4 (microM x s)(-1), respectively. The apparent quantum yields of photodissociation of the beta- and alpha-subunit within completely oxygenated R-state HbA differ from each other by a factor of 3.6 and are equal to 0.041 +/- 0.004 and 0.0114 +/- 0.0012, respectively. The apparent quantum yield of photodissociation of completely oxygenated R-state HbA is equal to 0.026 +/- 0.003.     

Effect of conserved intersubunit amino acid substitutions on Hfq protein structure and stability

Hfq is a thermostable RNA-binding bacterial protein that forms a uniquely shaped homohexamer. Based on sequence and structural similarity, Hfq belongs to the like-Sm (LSm) protein family. In spite of a rather high degree of homology between archaeal and eukaryotic LSm proteins, their quaternary structure is different, usually consisting of five to eight monomers. In this work, the importance of conserved intersubunit hydrogen bonds for the Hfq spatial organization was tested. The structures and stabilities for the Gln8Ala, Asn28Ala, Asp40Ala, and Tyr55Ala Hfq mutants were determined. All these proteins have the same hexamer organization, but their stability is different. Elimination of a single intersubunit hydrogen bond due to Gln8Ala, Asp40Ala, and Tyr55Ala substitutions results in decreased stability of the Hfq hexamer. Tyr55Ala Hfq as well as the earlier studied His57Ala Hfq has reduced protein thermostability, which seems to correspond to an opening of the protein hydrophobic core.     

Glucoamylase-like domains in the alpha- and beta-subunits of phosphorylase kinase

Phosphorylase kinase is a four-subunit enzyme involved in the regulation of glycogen breakdown. The traditional textbook view is that only the gamma subunit has enzymatic activity, whereas the other three subunits have a regulatory role. Evidence from homology searches and sequence alignments, however, shows that the alpha- and beta-subunits possess amino-terminal glucoamylase-like domains and suggests that they might possess a previously overlooked amylase activity. If true, this would have important implications for the understanding, diagnosis, and management of glycogen storage diseases. There is thus a clear need to test this hypothesis through enzymatic assays and structural studies.     

The functional impact of 1,570 individual amino acid substitutions in human OTC

1. Download : Download high-res image (149KB)
2. Download : Download full-size image

     

Down-regulation of the alpha- and beta-subunits of the calcium-activated potassium channel in human myometrium with parturition

Large-conductance, calcium-dependent potassium (BKCa) channels are implicated in maintaining uterine quiescence during pregnancy. The mechanisms whereby calcium sensitivity of the BKCa channel is dramatically removed at parturition remain unknown. The aim of the present study was to investigate whether this loss of calcium sensitivity of the BKCa channel with the onset of labor is associated with changes in the protein expression of the alpha- and/or beta-subunit or arises from a physical dissociation of the alpha-subunit from the beta-subunit. The beta-subunit is a key determinant of BKCa-channel Ca2+ sensitivity. Western blot analysis, using alpha- and beta-subunit-specific antibodies, detected bands of 110-125 and 36 kDa, respectively. Protein expression levels of the alpha-subunit in term labor myometrium were significantly reduced compared with term pregnancy without labor. Furthermore, alpha-subunit levels at term pregnancy were significantly increased relative to the nonpregnant state, whereas levels at preterm gestations were unchanged. Densitometric analysis demonstrated significantly decreased beta-subunit levels in term and preterm labor samples compared with term nonlabor samples. Immunoprecipitation studies revealed the presence of both the alpha- and beta-subunits in samples taken before or after the onset of labor. We conclude that during labor, the alpha-subunit is not physically uncoupled from the beta-subunit, but a decline occurs in the level of beta-subunit protein, which may underlie the loss of calcium and voltage sensitivity of the BKCa channel with labor. Furthermore, reduced beta-subunit protein in preterm labor myometrium implies that ion channels may also contribute to pathophysiological labor.     

Residual structure in large fragments of staphylococcal nuclease: effects of amino acid substitutions

In an attempt to develop a model of the denatured state of staphylococcal nuclease that can be analyzed experimentally under physiological conditions, a series of four large fragments of this small protein which extend from residues 1 to 103, 1 to 112, 1 to 128, and 1 to 136 have been generated through the overexpression of nuclease genes containing stop codons at defined positions. Large amounts of protein fragments were accumulated in induced cells and were purified by carrying out all fractionation steps in the presence of 6 M urea. The far-ultraviolet circular dichroism spectra of all four fragments suggested the presence of small to moderate amounts of residual structure. When the CD spectra were monitored as a function of concentrations of the tight-binding ligands Ca2+ and thymidine 3',5'-bisphosphate and the known affinity constants for wild-type nuclease (1-149) were used, apparent equilibrium constants of 160 and 2000 for the reversible denaturation reaction for fragments 1-136 and 1-128, respectively, were estimated. Four single and two double mutations, all of which exhibit unusual behavior in the full-length protein on solvent denaturation [Shortle, D., & Meeker, A. K. (1986) Proteins: Struct., Funct., Genet. 1, 81-89] and thermal denaturation [Shortle, D., Meeker, A. K., & Freire, E. (1988) Biochemistry 27, 4761-4768], were recombined into the 1-136 and 1-128 fragment expression vectors, and purified mutant fragments were characterized.(ABSTRACT TRUNCATED AT 250 WORDS)     

Estimation of the number of amino acid substitutions per site when the substitution rate varies among sites

A general model for estimating the number of amino acid substitutions per site (d) from the fraction of identical residues between two sequences (q) is proposed. The well-known Poisson-correction formula q = e ^–d corresponds to a site-independent and amino-acid-independent substitution rate. Equation q = (1 – e ^–2d )/2d, derived for the case of substitution rates that are site-independent, but vary among amino acids, approximates closely the empirical method, suggested by Dayhoff et al. (1978). Equation q = 1/(1 + d) describes the case of substitution rates that are amino acid-independent but vary among sites. Lastly, equation q = [ln(1 + 2d)]/2d accounts for the general case where substitution rates can differ for both amino acids and sites.     

Reciprocal amino acid substitutions in the evolution of homologous peptides

Amino acid sequences of peptides are often inferred from their amino acid compositions by comparison with homologous peptides of known sequence. The probabilities are considered that by such an approach errors are made due to the occurrence of balanced double changes, i.e. reciprocal substitutions, between two homologous peptides of identical compositions. Formulae are derived for the calculation of these probabilities, depending on peptide length and evolutionary distance. However, such calculations requiring too much computer time, the probabilities for reciprocal substitutions are estimated by simulation of evolutionary changes in peptides. It can be concluded from the resulting data that for many purposes the possible errors in amino acid sequences partially inferred from amino acid compositions are acceptably small.     

Estimation of the number of amino acid substitutions per site when the substitution rate varies among sites

Correspondence to: M. Nei 0871     

Effects of single-point amino acid substitutions on the structure and function neuraminidase proteins in influenza A virus

In order to clarify the effect of amino acid substitutions on the structure and function of the neuraminidase (NA) protein of influenza A virus, we introduced single-point amino acid substitutions into the NA protein of the A/Tokyo/3/67 (H2N2) strain using PCR-based random mutation. The rate of tolerant random one amino acid substitutions in the NA protein was 47%. Rates of tolerant substitutions for the stalk and for the surface and inner portion of the head region of the NA protein were 79, 54, and 19%, respectively. Deleterious changes, such as those causing the NA protein to stop at the Golgi/endoplasmic reticulum, were scattered throughout the protein. On the other hand, the ratio of mutations with which the NA protein lost neuraminidase activity, but was transported to the cell surface, decreased in proportion to the distance from the structural center of enzyme active site. In order to investigate the effect of accumulated amino acid substitutions on the structural character of the N2NA protein during evolution, the same amino acid substitutions were introduced by site-directed mutagenesis at 23 homologous positions on N2 proteins of A/Tokyo/3/67, A/Bangkok/15/85 (H3N2), and A/Mie/1/2004 (H3N2). The results showed a shift, or discordance, in tolerance at some of the positions. An increase in discordance was correlated with the interval in years between virus strains, and the discordance rate was estimated to be 0.6-0.7% per year.     

Stepping-stone spatial structure causes slow decay of linkage disequilibrium and shifts the site frequency spectrum

The symmetric island model with D demes and equal migration rates is often chosen for the investigation of the consequences of population subdivision. Here we show that a stepping-stone model has a more pronounced effect on the genealogy of a sample. For samples from a small geographical region commonly used in genetic studies of humans and Drosophila, there is a shift of the frequency spectrum that decreases the number of low-frequency-derived alleles and skews the distribution of statistics of Tajima, Fu and Li, and Fay and Wu. Stepping-stone spatial structure also changes the two-locus sampling distribution and increases both linkage disequilibrium and the probability that two sites are perfectly correlated. This may cause a false prediction of cold spots of recombination and may confuse haplotype tests that compute probabilities on the basis of a homogeneously mixing population.     

Haem degradation in human haemoglobin in vitro. Separation of the contribution of the alpha- and beta-subunits.

As part of an investigation into whether alpha 1-foetoprotein (alpha 1-FP) plays the same transport role in foetal serum as albumin does in the adult, the binding properties of both proteins were compared with respect to the binding of a series of compounds known to be bound by albumin's specific drug-binding sites. The binding of warfarin, phenylbutazone, azapropazone, diazepam, digitoxin and cholic acid by rat alpha 1-FP and serum albumin was studied by equilibrium dialysis at 4 degrees C. Rat alpha 1-FP was shown to have neither albumin's high-affinity site II (diazepam as marker) nor its site III (digitoxin and cholic acid as markers). High-affinity binding by alpha 1-FP was found for the specific markers (warfarin, phenylbutazone, azapropazone) of albumin's drug-binding site I. However, instead of albumin's one high-affinity site/molecule, a mean value of 0.5 site/molecule was obtained with rat alpha 1-FP. Charcoal treatment at neutral pH of rat serum albumin did not affect its measured binding properties, but treatment of the alpha 1-FP led to an increased affinity for warfarin, phenylbutazone and azapropazone without a change in the measured number of sites, indicating competition for binding at this site by (an) endogenous ligand(s). These results are discussed in terms of the structures of the two proteins and with respect to the physiological implications of the differences found.     

Code dependent conservation of the physico-chemical properties in amino acid substitutions

The frequency of amino acid replacements in families of typical proteins has been elegantly analyzed by Argyle (1980) showing that the most frequent replacements involve a conservation of the amino acid chemical properties. The cyclic arrangement of the twenty amino acids resulting from the most frequent replacements has been described as an amino acid chemical ring. In this work, a novel amino acid replacement frequency ring is proposed, for which a conservation of over 90% of the most general physico-chemical properties can be deduced. The amino acid chemical similarity ring is also analyzed in terms of the genetic code base probability changes, showing that the discrepancy that exists between the standard deviation value of the amino acid replacement frequency matrix and its respective ideal value is almost equal to that deduced from the corresponding base codon replacement probability matrices. These differences are finally evaluated and discussed in terms of the restrictions imposed by the structure of the genetic code and the physico-chemical dissimilarities between some codons of amino acids which are chemically similar.