Structure and ligand binding properties of leucine 29(B10) mutants of human myoglobin

Expression levels of myoglobin in muscle and non-muscle tissues of rainbow trout Oncorhynchus mykiss, a hypoxia intolerant species

Myoglobin (Mb) is one of the most intensively studied intracellular respiratory muscle proteins. Since the discovery of the fascinating fact that Mb is not confined only to oxidative muscle tissues but also is co-localized in different non-muscle tissues of cyprinids, hypoxia tolerant cyprinids have been established as the model teleost. Mb both at mRNA and protein levels have been reported in this study for the first time from a number of muscle and non-muscle tissues of rainbow trout Oncorhynchus mykiss, a hypoxia intolerant species. Mb transcript levels were high in the heart and slow skeletal muscle, and were comparatively high in the gonad and gill among the non-muscle tissues. Western-blotting by using anti-rainbow trout Mb peptide rabbit antibody detected Mb protein in the muscles and several non-muscle tissues. By both RNA in situ hybridization and immunofluorescence, Mb was localized in the cardiomyocytes and oxidative muscle fibers. On the other hand, Mb both at mRNA and protein levels was restricted to the lamellar epithelial cells of the gill, epithelial layers of hepato-biliary duct, neurons and endothelial cells of brain, ooplasm of gonad, kidney tubules, endothelial cells, and epithelial layer of intestine. Neuroglobin isoform 1 and 2 mRNAs along with Mb mRNA were localized in the granular layer of cerebellum. Considering the previous data reported for cyprinids, the expression sites of Mb in the muscle and non-muscle tissues of teleost could be universal, where Mb concerted with the other globins might play meaningful physiological roles.

Application of high pressure laser flash photolysis in studies on selected hemoprotein reactions

This article focuses on the application of high pressure laser flash photolysis for studies on selected hemoprotein reactions with the objective to establish details of the underlying reaction mechanisms. In this context, particular attention is given to the reactions of small molecules such as dioxygen, carbon monoxide, and nitric oxide with selected hemoproteins (hemoglobin, myoglobin, neuroglobin and cytochrome P450(cam)), as well as to photo-induced electron transfer reactions occurring in hemoproteins (particularly in various types of cytochromes). Mechanistic conclusions based on the interpretation of the obtained activation volumes are discussed in this account.

Human myoglobin recognition of oxygen: dynamics of the energy landscape

Femtosecond to nanosecond dynamics of O(2) rebinding to human WT myoglobin and its mutants, V68F and I107F, have been studied by using transient absorption. The results are compared with NO rebinding. Even though the immediate environment around the heme binding site is changed by the mutations, the picosecond geminate rebinding of oxygen is at most minimally affected. On the other hand, the V68F (E11) mutation causes drastic differences in rebinding on the nanosecond time scale, whereas the effect of the I107F (G8) mutation remains relatively small within our 10-ns time window. Unlike traditional homogeneous kinetics and molecular dynamics collisional simulations, we propose a "bifurcation model" for populations of directed and undirected dynamics on the ultrafast time scale, reflecting the distribution of initial protein conformations. The major mutation effect occurs on the time scale on which global protein conformational change is possible, consistent with transitions between the conformations of directed and undirected population playing a role in the O(2) binding. We discuss the relevance of these findings to the bimolecular function of the protein.

Ligand migration in human myoglobin: steric effects of isoleucine 107(G8) on O(2) and CO binding

To investigate the ligand pathway in myoglobin, some mutant myoglobins, in which one of the amino acid residues constituting a putative ligand-docking site, Ile107, is replaced by Ala, Val, Leu, or Phe, were prepared and their structural and ligand binding properties were characterized. The kinetic barrier for the ligand entry to protein inside was lowered by decreasing the side-chain volume at position 107, indicating that the bulky side chain interferes with the formation of the activation state for the ligand migration and the free space near position 107 would be filled with the ligand in the activation state. Another prominent effect of the reduced side-chain volume at position 107 is to stabilize the ligand-binding intermediate state. Because the stabilization can be ascribed to decrease of the positive enthalpy, the enlarged free space near position 107 would relieve unfavorable steric interactions between the ligand and nearby amino acid residues. The side-chain volume at position 107, therefore, is crucial for the kinetic barrier for the ligand migration and free energy of the ligand-binding intermediate state, which allows us to propose that some photodissociated O(2) moves toward position 107 to be trapped and then expelled to the solvent.

Unusual pressure effects on ligand rebinding to the human myoglobin Leucine 29 mutants

Using high pressure flash photolysis, we revealed that the side chain of Leu(29) controls the reaction volume of the ligand migration process in myoglobin, which is the primary factor for the unusual activation volume of ligand binding in some Leu(29) mutants. As we previously reported (Adachi, S., Sunohara, N., Ishimori, K., and Morishima, I. (1992) J. Biol. Chem. 267, 12614-12621), CO bimolecular rebinding in the L29A mutant was unexpectedly decelerated by pressurization, suggesting that the rate-determining step is switched to ligand migration. However, very slow CO bimolecular rebinding of the mutants implies that bond formation is still the rate-determining step. To gain further insights into effects of the side chain on ligand binding, we prepared some new Leu(29) mutants to measure the CO and O(2) rebinding reaction rates under high hydrostatic pressure. CO bimolecular rebinding in the mutants bearing Gly or Ser at position 29 was also decelerated upon pressurization, resulting in apparent positive activation volumes (DeltaV), as observed for O(2) binding. Based on the three-state model, we concluded that the increased space available to ligands in these mutants enhances the volume difference between the geminate and deoxy states (DeltaV(32)), which shifts the apparent activation volume to the positive side, and that the apparent positive activation volume is not due to contribution of the ligand migration process to the rate-determining step.

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.

Structural and functional roles of heme binding module in globin proteins: identification of the segment regulating the heme binding structure

To investigate structural and functional significance of a newly proposed structural unit in globins, the "heme binding module", we synthesized a "heme binding module"-substituted chimeric globin and characterized its function and structure. In our previous study we proposed that the heme binding module, corresponding to the segment from Leu(F1) to Phe(G5) in hemoglobin alpha-subunit, plays a key role in constructing the heme proximal structure in globins. The replacement of the heme binding module in myoglobin with that of hemoglobin alpha-subunit converted the absorption spectra into that of the alpha-subunit, and, in the resonance Raman spectra, the vibration mode characteristic of myoglobin completely disappeared after the module replacement. The hyperfine-shifted NMR resonances for the cyanide-bound form of the module-substituted myoglobin also revealed that the orientation of the axial histidine is close to that of the alpha-subunit rather than that of myoglobin, while the deviations of the resonance positions of the NMR signals from the amino acid residues located in the distal site were subtle, supporting the preferential structural alterations in the heme proximal site. The present finding for the structural alterations in the module-substituted myoglobin confirms that the heme binding module can be a segment regulating the heme proximal structure in globin proteins.

The effects of heme pocket hydrophobicity on the ligand binding dynamics in myoglobin as studied with leucine 29 mutants

To examine the effects of heme pocket hydrophobicity on the ligand binding in myoglobin, some artificial mutants of human myoglobin have been prepared, in which less hydrophobic amino acid residue (Ala, Gly, Ser) is located at the Leu29 (10th residue of the B helix) position. CO rebinding rates for the mutants were markedly decelerated, while the 1H, and 15N NMR spectra of the mutants show that the structural changes around the heme iron for these mutants are rather small. The kinetic and structural properties of the mutants indicate that the ligand binding rate depends on the hydrophobicity inside the heme cavity for these mutants in addition to the volume of the side chain at the 29-position. On the basis of the IR stretching frequency of liganded CO, invasion of water molecules into the heme pocket in the mutants is suggested, which would be induced by the decrease in the hydrophobicity due to the amino acid substitution. A slight red shift of the position of the Soret peak for the serine mutant L29S also supports the reduced hydrophobicity inside the heme cavity. We can concluded that, together with the kinetic properties of the mutants, the hydrophobicity of the heme pocket is one of the key factors in regulating the ligand binding to the heme iron.

Phe-46(CD4) orients the distal histidine for hydrogen bonding to bound ligands in sperm whale myoglobin

The role of Phe-46(CD4) in modulating the functional properties of sperm whale myoglobin was investigated by replacing this residue with Leu, Ile, Val, Ala, Trp, Tyr, and Glu. This highly conserved amino acid almost makes direct contact with the distal histidine and has been postulated to affect ligand binding. The overall association rate constants for CO, O2, and NO binding were little affected by decreasing the size of residue 46 step-wise from Phe to Leu to Val to Ala. In contrast, the rates of CO, O2, and NO dissociation increased 4-, 10-, and 25-fold, respectively, for the same series of mutants, causing large decreases in the affinity of myoglobin for all three diatomic gases. The rates of autooxidation at 37 degrees C, pH 7.0 increased dramatically from approximately 0.1-0.3 h-1 for wild-type, Tyr-46, and Trp-46 myoglobins to 1.5, 5.2, 4.9, and 5.0 h-1 for the Leu-46, Ile-46, Val-46 and Ala-46 mutants, respectively. Rates of NO and O2 geminate recombination were measured using 35 ps and 9 ns laser excitation pulses. Decreasing the size of residue 46 causes significant decreases in the extent of both picosecond and nanosecond rebinding processes. High resolution structures of Leu-46 and Val-46 metmyoglobins, Val-46 CO-myoglobin, and Val-46 deoxymyoglobin were determined by X-ray crystallography. When Phe-46 is replaced by Val, the loss of internal packing volume is compensated by (1) contraction of the CD corner toward the core of the protein, (2) movement of the E-helix toward the mutation site, (3) greater exposure of the distal pocket to intruding solvent molecules, and (4) large disorder in the position of the side chain of the distal histidine (His-64). In wild-type myoglobin, the van der Waals contact between C zeta of Phe-46 and C beta of His-64 appears to restrict rotation of the imidazole side chain. Insertion of Val at position 46 relieves this steric restriction, allowing the imidazole side chain to rotate about the C alpha - C beta bond toward the surface of the globin and about the C beta - C gamma bond toward the space previously occupied by the native Phe-46 side chain. This movement disrupts hydrogen bonding with bound ligands, causing significant decreases in affinity, and opens the distal pocket to solvent water molecules, causing marked increases in the rate of autooxidation.(ABSTRACT TRUNCATED AT 400 WORDS)

Spectroscopic study of Ser92 mutants of human myoglobin: hydrogen bonding effect of Ser92 to proximal His93 on structure and property of myoglobin

Neutron diffraction studies have demonstrated that the hydroxyl group oxygen of Ser92(F7) is hydrogen bonded to the proximal His93(48) N epsilon H proton in myoglobin (Mb) [Cheng, X., & Shoenborn, B. P. (1991) J. Mol. Biol. 220, 381-399]. In order to examine the importance of this hydrogen bond, Ser92 was replaced with Ala and Asp in human Mb. By comparing the optical, 1H-NMR, resonance Raman, and IR spectra of Mb(S92A) in several spin and oxidation states with those of wild-type Mb, it was found that the mutation causes a structural change on the heme proximal side but not on the distal side. Comparison of the NMR spectra of the cyanomet form of Mb(S92A) and Mb(WT) suggests that the imidazole plane of His93 rotates somewhat around the Fe-N delta (His93) bond upon loss of the hydrogen bond between His93 and Ser92. The 2D 1H-NMR measurements of the CO complexes show that mutation of Ser92 to Ala changes the relative position of the His97 imidazole group to the heme plane, but the change is not so drastic as reported in the crystal data of Ser92 mutant of pig Mb [Smerdon et al. (1993) Biochemistry 32, 5132-5138]. On the other hand, ligand (CO, O2) binding is only slightly affected by this mutation. From these results, we conclude that the Ser92-His93 hydrogen bond maintains the protein structure of the proximal heme pocket, but it does not strongly affect the electronic structure of the heme as well as of the His93 imidazole ring.(ABSTRACT TRUNCATED AT 250 WORDS)

FTIR analysis of the interaction of azide with horse heart myoglobin variants

The interaction of azide with variants of horse heart myoglobin (Mb) has been characterized by Fourier transform infrared (FTIR), electron paramagnetic resonance (EPR), and UV-VIS absorption spectroscopy and by molecular modeling calculations. Distal histidine variants (His64Thr, His64Ile, His64Lys) and charged surface variants (Val67Arg, Lys45Glu, Lys45Glu/Lys63Glu) were included in this study. All variants, with the exception of Val67Arg, have a lower azide affinity than the wild-type protein. Analysis of the temperature dependence of the FTIR spectra (277-313 K) revealed that the wild-type protein and all variants exhibit a high-spin/low-spin equilibrium. Introduction of positively charged amino acid residues shifts nu max for the low-spin form to higher energy while negatively charged residues shifted this maximum to lower energy. The low azide binding affinity exhibited by the His64Thr and His64Ile variants is accompanied by a shift of the nu max for the low-spin infrared band to lower energy and by a significant increase in the corresponding half-bandwidths. This observation indicates greater mobility of the bound azide ligand in these variants. The His64Lys variant exhibits two infrared bands attributable to low-spin forms that are assigned to two different conformations of the lysyl residue. In one conformation, the lysine is proposed to form a hydrogen bond with the bound azide similar to that proposed to occur between the distal histidine and bound azide, and in the other conformation no interaction occurs.

Application of molecular dynamics and free energy perturbation methods to metalloporphyrin-ligand systems II: CO and dioxygen binding to myoglobin

The protein contribution to the relative binding affinity of the ligands CO and O2 toward myoglobin (Mb) has been simulated using free energy perturbation calculations. The tautomers of the His E7 residue are different for the oxymyoglobin (MbO2) and carboxymyoglobin (MbCO) systems. This was modeled by performing two-step calculations that mutate the ligand and mutate the His E7 tautomers in separate steps. Differences in hydrogen bonding to the O2 and CO ligands were incorporated into the model. The O2 complex was calculated to be 2-3 kcal/mol more stable than the corresponding CO complex when compared to the same difference in an isolated heme control. This value agrees well with the experimental value of 2.0 kcal/mol. In qualitative agreement with experiments, the Fe-C-O bond is found to be bent (theta = 159.8 degrees) with a small tilt (theta = 6.2 degrees). The contributions made by each of the 29 residues--within the 9.0-A radius of the iron atom--to the free energy difference are separated into van der Waals and electrostatic contributions; the latter contributions are dominant. Aside from the proximal histidine and the heme group, the residues having the largest difference in free energy in mutating MbO2-->MbCO are His E7, Phe CD1, Phe CD4, Val E11, and Thr E10.

Interpretation of hyperfine shift patterns in ferricytochromes b5 in terms of angular position of the heme: a sensitive probe for peripheral heme protein interactions

The 1H-NMR hyperfine shift pattern of the heme in a variety of low-spin ferricytochromes b5 has been analyzed in terms of the angular position of the prosthetic group within a structurally and magnetically-conserved protein matrix. A simple model is presented in which the changes in the spread of the predominantly contact shifted methyl and predominantly dipolar shifted meso-H signals of the heme, as well as shift trends for individual signals, provide sensitive indicators of the orientation of the heme relative to the orbital hole (singly-occupied d orbital), which in turn is related to the rhombic magnetic axes. The invariance of the axial His and non-coordinated residue hyperfine shifts show that it is the heme within a relatively rigid protein matrix, rather than the magnetic coordinate system, which is displaced angularly about the heme normal in order to accommodate variations in the polypeptide, orientation of the heme about the alpha,gamma-meso axis, and the length of heme carboxylate chains. Native heme shows increased counterclockwise rotation about the heme normal in the order rat-->beef-->chicken ferricytochrome b5, which is attributed largely to increased bulk of a variable sequence hydrophobic cluster consisting of residues 23, 25 and 32. The two alternate heme orientations about the alpha,gamma-meso axis are shown to also differ by rotation about the heme normal. A semiquantitative estimate of the degree of angular accommodation based on the spread of the meso-H rhombic dipolar shifts indicate rotations of 2-10 degrees. Possible functional consequences of such angular accommodation in relation to the role of these proteins in electron transfer are discussed.

Investigations of the thermostability of rubredoxin models using molecular dynamics simulations

The affects of differences in amino acid sequence on the temperature stability of the three-dimensional structure of the small beta-sheet protein, rubredoxin (Rd), was revealed when a set of homology models was subjected to molecular dynamics simulations at relatively high temperatures. Models of Rd from the hyperthermophile, Pyrococcus furiosus (Pf), an organism that grows optimally at 100 degrees C, were compared to three mesophilic Rds of known X-ray crystal structure. Simulations covering the limits of known Rd thermostabilities were carried out at temperatures of 300 K, 343 K, 373 K, and 413 K. They suggest that Rd stability is correlated with structural dynamics. Because the dynamic behavior of three Pf Rd models was consistently different from the dynamic behavior of the three mesophilic Rd structures, detailed analysis of the temperature-dependent dynamic behavior was carried out. The major differences between the models of the protein from the hyperthermophile and the others were: (1) an obvious temperature-dependent transition in the mesophilic structures not seen with the Pf Rd models, (2) consistent AMBER energy for the Pf Rd due to differences in nonbonded interaction terms, (3) less variation in the average conformations for the Pf Rd models with temperature, and (4) the presence of more extensive secondary structure for the Pf Rd models. These unsolvated dynamics simulations support a simple, general hypothesis to explain the hyperthermostability of Pf Rd. Its structure simplifies the conformational space to give a single minimum accessible over an extreme range of temperatures, whereas the mesophilic proteins sample a more complex conformational space with two or more minima over the same temperature range.

Transition state of an unfolding step in human cyanomet myoglobin

We studied an unfolding step of cyanomet myoglobin (Mb) unfolding, for demonstrating dynamical structural changes in the transition state of the process. Three leucine-->alanine mutant Mbs (L29A, L72A and L104A) were prepared for this study. The urea-induced largely monophasic process was monitored by absorption spectroscopy. Linear relations between [urea] and the activation energy (delta G not equal to) of the relaxation for all the Mbs showed that the slope m not equal to urea (= delta(delta G not equal to)/delta[urea])) was altered by either reduction of pH or the L-->A mutations. Thermodynamic interpretations of the changes in m not equal to urea led to a conclusion that the exposed surface area of Mb in the transition state was determined by both protein-core stability and pH conditions. We also performed urea- and acid-denaturation experiments, and gave some inspections on differences between mutational effects on the structure of the transition state and the denatured state.