Theoretical study of the discrimination between O(2) and CO by myoglobin

Combined quantum chemical and molecular mechanics geometry optimisations have been performed on myoglobin without or with O(2) or CO bound to the haem group. The results show that the distal histidine residue is protonated on the N(epsilon 2) atom and forms a hydrogen bond to the haem ligand both in the O(2) and the CO complexes. We have also re-refined the crystal structure of CO[bond]myoglobin by a combined quantum chemical and crystallographic refinement. Thereby, we probably obtain the most accurate available structure of the active site of this complex, showing a Fe[bond]C[bond]O angle of 171 degrees, and Fe[bond]C and C[bond]O bond lengths of 170-171 and 116-117 pm. The resulting structures have been used to calculate the strength of the hydrogen bond between the distal histidine residue and O(2) or CO in the protein. This amounts to 31-33 kJ/mol for O(2) and 2-3 kJ/mol for CO. The difference in hydrogen-bond strength is 21-22 kJ/mol when corrected for entropy effects. This is slightly larger than the observed discrimination between O(2) or CO by myoglobin, 17 kJ/mol. We have also estimated the strain of the active site inside the protein. It is 2-4 kJ/mol larger for the O(2) complex than for the CO complex, independent of which crystal structure the calculations are based on. Together, these results clearly show that myoglobin discriminates between O(2) and CO mainly by electrostatic interactions, rather than by steric strain.     

Will the real FeCO please stand up?

 The paradigm that nature protects us from CO poisoning by forcing the bound CO to bend over in heme proteins, thereby reducing its binding affinity, is now in textbooks, but is nevertheless problematic. Results from vibrational spectroscopy give no evidence for bent CO, although X-ray crystallography continues to indicate appreciable distortions in myoglobin. However, the energetic significance of the discrepancy is doubtful, since new Density Functional Theory calculations indicate that much less energy is required to distort the CO than had been thought, perhaps 2 kcal/mol or less. Binding studies on site-directed mutants of myoglobin suggest that steric hindrance by the distal histidine is worth ca. 1 kcal/mol. While the distal histidine does account for the discrimination by Mb against CO and in favor of O₂, most of the effect is due to its H-bond with bound O₂. Received, accepted: 23 May 1997     

Myoglobin discriminates between O2, NO, and CO by electrostatic interactions with the bound ligand

 Most biological substrates have distinctive sizes, shapes, and charge distributions which can be recognized specifically by proteins. In contrast, myoglobin must discriminate between the diatomic gases O₂, CO, and NO which are apolar and virtually the same size. Selectivity occurs at the level of the covalent Fe-ligand complexes, which exhibit markedly different bond strengths and electrostatic properties. By pulling a water molecule into the distal pocket, His64(E7)¹ inhibits the binding of all three ligands by a factor of ∼10 compared to that observed for protoheme-imidazole complexes in organic solvents. In the case of O₂ binding, this unfavorable effect is overcome by the formation of a strong hydrogen bond between His64(E7) and the highly polar FeO₂ complex. This favorable electrostatic interaction stabilizes the bound O₂ by a factor of ∼1000, and the net result is a 100-fold increase in overall affinity compared to model hemes or mutants with an apolar residue at position 64. Electrostatic interaction between FeCO and His64 is very weak, resulting in only a two- to three-fold stabilization of the bound state. In this case, the inhibitory effect of distal pocket water dominates, and a net fivefold reduction in K _CO is observed for the wild-type protein compared to mutants with an apolar residue at position 64. Bound NO is stabilized ∼tenfold by hydrogen bonding to His64. This favorable interaction with FeNO exactly compensates for the tenfold inhibition due to the presence of distal pocket water, and the net result is little change in K _NO when the distal histidine is replaced with apolar residues. Thus, it is the polarity of His64 which allows discrimination between the diatomic gases. Direct steric hindrance by this residue plays a minor role as judged by: (1) the independence of K _O2, K _CO, and K _NO on the size of apolar residues inserted at position 64, and (2) the observation of small decreases, not increases, in CO affinity when the mobility of the His64 side chain is increased. Val68(E11) does appear to hinder selectively the binding of CO. However, the extent is no more than a factor of 2–5, and much smaller than electrostatic stabilization of bound O₂ by the distal histidine. Received, accepted: 23 May 1997     

Modulating carbon monoxide binding affinity and kinetics in myoglobin: the roles of the distal histidine and the heme pocket docking site

 Myoglobin has long served as a model system for understanding the relations between protein structure, dynamics, and function. Its ability to discriminate between toxic CO and vital O₂, two small ligands that are almost equivalent in size and dipole moment, has attracted much attention. To understand discrimination and reversible ligand-binding in Mb, both the bound state and the "docked" state that leads to binding need to be studied. We have reported previously the nearly linear Fe–C–O geometry of bound CO and the nearly orthogonal geometry of docked CO [Lim et al. (1995), Science 269 : 962]. With the exception of X-ray structures, a preponderance of evidence points to a nearly linear Fe–C–O geometry and calls into question the proposal that the highly conserved distal histidine forces CO to bind in a nonoptimal geometry. The differences between the bound CO structures determined using IR and X-ray methods might arise from a water molecule hydrogen bonded to the distal histidine in some of the unit cells. Discrimination by Mb is manifested not only thermodynamically but also kinetically. Time-resolved CO rebinding studies that compare Mb with microperoxidase suggest that the heme pocket docking site in Mb exerts steric control of the ligand rebinding rate, slowing the rate of CO binding by a factor of more than 10⁴. Received, accepted: 23 May 1997     

Effects of local protein environment on the binding of diatomic molecules to heme in myoglobins. DFT and dispersion-corrected DFT studies

The heme-AB binding energies (AB = CO, O₂) in a wild-type myoglobin (Mb) and two mutants (H64L, V68N) of Mb have been investigated in detail with both DFT and dispersion-corrected DFT methods, where H64L and V68N represent two different, opposite situations. Several dispersion correction approaches were tested in the calculations. The effects of the local protein environment were accounted for by including the five nearest surrounding residues in the calculated systems. The specific role of histidine-64 in the distal pocket was examined in more detail in this study than in other studies in the literature. Although the present calculated results do not change the previous conclusion that the hydrogen bonding by the distal histidine-64 residue plays a major role in the O₂/CO discrimination by Mb, more details about the interaction between the protein environment and the bound ligand have been revealed in this study by comparing the binding energies of AB to a porphyrin and the various myoglobins. The changes in the experimental binding energies from one system to another are well reproduced by the calculations. Without constraints on the residues in geometry optimization, the dispersion correction is necessary, since it improves the calculated structures and energetic results significantly.     

Myoglobin models and steric origins of the discrimination between O2 and CO

 Synthetic models of the myoglobin active site have provided much insight into factors that affect CO and O₂ binding in the proteins. "Capped" and "pocket" metal porphyrin systems have been developed to probe how steric factors affect ligand binding and ultimately to elucidate important aspects of the mechanism of CO discrimination in the proteins. These model porphyrins are among the most thoroughly characterized systems to date. From the twenty-one known crystal structures, analysis of the types of distortion that occur upon ligand binding under the cap, including porphyrin doming and ruffling, lateral and horizontal movement of the cap, and bending and tilting of the Fe–C–O bond, provides an indication of how steric interactions will affect structure in Hb and Mb. The model porphyrin systems discussed range from those that discriminate against O₂ binding compared to biological systems to those with similar CO and O₂ binding strength to myoglobin, and also to those that bind both O₂ and CO very weakly or not at all. The primary type of distortion observed upon CO binding is vertical or lateral movement of the cap and some ruffling of the porphyrin plane. Minimal bending or tilting of the M–C–O bond is observed, suggesting that the Fe–C–O bending that has been found from crystal structures of the hemoproteins is unlikely. Received, accepted: 23 May 1997     

The dimer-monomer conversion of Cerithidea myoglobin coupled with the heme iron oxidation

With regard to the distal (E7) residue, gastropod sea mollusc contains both types of myoglobin, one with and the other lacking the distal histidine. We have isolated a myoglobin from the radular muscle of Cerithidea rhizophorarum, a small whelk found on the Japanese coast. Unlike Aplysia myoglobin having a single histidine residue at position 95, Cerithidea myoglobin contains three histidines at positions 48, 66 and 98. Moreover, Cerithidea MbO₂ exists as homodimers and is oxidized, not to the usual form of metMb but to the hemichrome monomers. It was also found that the hemichrome monomers thus produced can easily be converted back to the dimerized oxy-form, if the ferric protein was reduced carefully with a slight excess of sodium hydrosulfite. This dimer-monomer conversion coupled with the heme-iron oxidation in Cerithidea myoglobin is very unique, and the distal histidine at position 66 is probably responsible for its reversible formation of hemichrome from the ferric met-form that occurred transiently in due course of the oxidation reaction of Cerithidea MbO₂     

The distal residue-CO interaction in carbonmonoxy myoglobins: a molecular dynamics study of two distal histidine tautomers.

Four independent 90 ps molecular dynamics simulations of sperm-whale wild-type carbonmonoxy myoglobin (MbCO) have been calculated using a new AMBER force field for the haem prosthetic group. Two trajectories have the distal 64N delta nitrogen protonated, and two have the 64N epsilon nitrogen protonated; all water molecules within 16 A of the carbonyl O are included. In three trajectories, the distal residue remains part of the haem pocket, with the protonated distal nitrogen pointing into the active site. This is in contrast with the neutron diffraction crystal structure, but is consistent with the solution phase CO stretching frequencies (upsilon CO) of MbCO and various of its mutants. There are significant differences in the "closed" pocket structures found for each tautomer: the 64N epsilon H trajectories both show stable distal-CO interactions, whereas the 64N delta H tautomer) has a weaker interaction resulting in a more mobile distal side chain. One trajectory (a 64N delta H tautomer) has the distal histidine moving out into the "solvent", leaving the pocket in an "open" structure, with a large unhindered entrance to the active site. These trajectories suggest that the three upsilon CO frequencies observed for wild-type MbCO in solution, rather than representing significantly different Fe-C-O geometries as such, arise from three different haem pocket structures, each with different electric fields at the ligand. Each pocket structure corresponds to a different distal histidine conformer: the A3 band to the 64N epsilon H tautomer, the A1,2 band to the 64N delta H tautomer, and the A0 band to the absence of any significant interaction with the distal side chain.     

Discrimination between oxygen and carbon monoxide and inhibition of autooxidation by myoglobin. Site-directed mutagenesis of the distal histidine

Sperm whale myoglobin mutants were constructed using site-directed mutagenesis to replace the highly conserved distal histidine residue (His(E7)-64). His-64 was substituted with Gly, Val, Phe, Cys, Met, Lys, Arg, Asp, Thr, and Tyr, and all 10 mutant proteins expressed to approximately 10% of the total soluble cell protein in Escherichia coli as heme containing myoglobin. With the exception of His-64----Tyr, which did not form a stable oxygen (O2) complex, all mutant proteins could be reduced and bound O2 and carbon monoxide (CO) reversibly. However, removal of the distal histidine increased the rate of autooxidation 40-350-fold. The His-64----Gly, Val, Phe, Met, and Arg mutants all showed markedly increased O2 dissociation rate constants which were approximately 50-1500-fold higher than those for wild-type myoglobin and increased O2 association rate constants which were approximately 5-15-fold higher than those for the native protein. All mutants studied (except His-64----Tyr) showed approximately 10-fold increased CO association rates and relatively unchanged CO dissociation rates. These altered O2 and CO association and dissociation rate constants resulted in 3-14-fold increased CO affinities, 10-200-fold decreased O2 affinities, and 50-380-fold greater M (KCO/KO2) values for the mutants compared to the wild-type protein. Thus, the distal histidine of myoglobin discriminates between CO and O2 binding by both sterically hindering bound CO and stabilizing bound O2 through hydrogen bonding. The increased autooxidation rates observed for the mutants appear to be due to a decrease in oxygen affinity and an increase in solvent anion accessibility to the distal pocket.     

EPR characterization of the stereochemistry of the distal heme pocket of the engineered human myoglobin mutants.

Recombinant human myoglobin mutants with the distal histidine residue replaced by Leu, Val, or Gln residues have been prepared by site-directed mutagenesis and expression in Escherichia coli. The recombinant apomyoglobin proteins have been successfully reconstituted with cobaltous protoporphyrin IX to obtain cobalt myoglobin mutant proteins, and the role of the distal histidine residue on the interaction between the bound ligand and the myoglobin molecule has been studied by EPR spectroscopy. We found that the distal histidine residue is significant in the orientation of the bound oxygen molecule. Low temperature photolysis experiments on both oxy cobalt proteins and ferric nitric oxide complexes indicated that the nature of the photolyzed form depends on the steric crowding of the distal heme pocket. To our surprise, the distal Leu mutant has a less restricted, less sterically crowded distal heme pocket than that of the distal Val mutant myoglobin, despite the fact that Leu has a larger side chain volume than Val. Our results demonstrate that the distal heme pocket steric crowding is not necessarily related to the side chain volume of the E7 residue.     

Modification of the heme distal side in myoglobin by cyanogen bromide. Heme environmental structures and ligand binding properties of the modified myoglobin

Met, deoxy, and CO forms of myoglobin (Mb) react with a stoichiometric amount of cyanogen bromide (BrCN) to cause substantial changes in the 1H NMR, optical absorption, and infrared spectra. These spectral changes were interpreted as arising from the substantial alterations in the heme environments, most probably due to the modification of the histidine residue at the heme distal side. It is also revealed that the modified Mb does not combine with some exogenous ligands such as CN-, CH3NH2, and O2, although it does with N-3 or CO. These unique ligand binding properties are also discussed with relevance to a role of the distal histidine in stabilizing the coordinated ligand through a hydrogen bond and to a steric constraint.     

Spin-state equilibria and axial ligand bonding in FixL hydroxide: a resonance raman study

 Vibrational assignments for the Fe-OH unit of ferric alkaline forms of two deletion derivatives of Rhizobium meliloti FixL, FixL*, a functional O₂-sensing heme kinase, and FixLN, which contains only the heme domain, are made. Appearance of ²H- and ¹⁸O-sensitive Raman bands indicates that the heme group of FixL binds hydroxide as a distal ligand to form a six-coordinate complex. The alkaline FixLs are distributed between high- and low-spin states. The high- and low-spin bands corresponding to the ν (Fe-OH) modes occur at 479 and 539 cm^–1, respectively. Low temperature favors formation of the low-spin complex, indicative of a thermal spin-state equilibrium. The ν (Fe-OH) frequencies of FixLN and FixL* are 11 to 18 cm^–1 lower than those observed for the respective vibrations in alkaline myoglobin and hemoglobin. The weaker Fe-OH bond in the FixLs is attributed to a lack of hydrogen bonding on the distal side of the heme pocket. Received: 20 November 1997 / Accepted: 2 March 1998     

The role of entropy in the discrimination between CO and O2 in myoglobin

Using stopped-flow rapid mixing and flash photolysis techniques, the dissociation rate coefficients of horse carbonmonoxy myoglobin (hMbCO) and oxygenated myoglobin (hMbO₂) in aqueous solution have been determined as a function of temperature between 274 and 342 K. From the Arrhenius plot, an activation enthalpy for dissociation of 74 kJ/mol was obtained for both ligands. The pronounced kinetic differences arise from markedly different pre-exponentials. We compare the Arrhenius parameters with those of the association reaction, as measured at cryogenic temperatures. In our analysis we conclude that the entropy loss upon binding of O₂ is twice as large as that for CO. Taking reasonable estimates for the frequency factor, the transition state entropy in hMbO₂ is located roughly half way in between the entropies of the bound and unbound states. By contrast, the entropy of the transition state in hMbCO appears to be identical to that of the bound state. Possible structural reasons for the different behavior are discussed. Received: 13 January 1997 / Accepted: 24 April 1997     

Resonance Raman studies of CO and O2 binding to elephant myoglobin (distal His(E7)----Gln)

Carbon monoxide and dioxygen were employed as resonance Raman-visible ligands for probing the nature of the heme-binding site in elephant myoglobin, which has glutamine in the distal position (E7) instead of the usual histidine. The distal histidine (E7) residue has been thought to be responsible for weakening carbon monoxide binding to hemoproteins. It is of interest to see how the His(E7)----Gln replacement affects such parameters as nu(Fe-N epsilon), nu(Fe-CO), delta(Fe-C-O), nu(C-O), delta(Fe-O-O), and nu(O-O) vibrational frequencies and relative intensities. Elephant myoglobin has a CO affinity approximately 6 times higher than that for human/sperm whale myoglobin (Mb). If this enhanced affinity were solely due to the removal of some of the steric hindrance that normally tilts the CO off the heme axis, one would expect the nu(Fe-CO) frequency to decrease and the nu(C-O) frequency to increase relative to the corresponding values in sperm whale Mb. However, the opposite was found. In addition, strong enhancement of the Fe-C-O bending mode was observed. These results suggest that the Fe-C-O linkage remains distorted. In elephant Mb, new interactions resulting from the conformational change accompanying ligand binding may be responsible for the increased CO binding. Similar spectra were obtained for elephant and sperm whale oxymyoglobin. This suggests that the interactions of bound O2 are not markedly affected by the glutamine replacement.     

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.     

Structural features promoting dioxygen production by Dechloromonas aromatica chlorite dismutase

Chlorite dismutase (Cld) is a heme enzyme capable of rapidly and selectively decomposing chlorite (ClO₂ ⁻) to Cl⁻ and O₂. The ability of Cld to promote O₂ formation from ClO₂ ⁻ is unusual. Heme enzymes generally utilize ClO₂ ⁻ as an oxidant for reactions such as oxygen atom transfer to, or halogenation of, a second substrate. The X-ray crystal structure of Dechloromonas aromatica Cld co-crystallized with the substrate analogue nitrite (NO₂ ⁻) was determined to investigate features responsible for this novel reactivity. The enzyme active site contains a single b-type heme coordinated by a proximal histidine residue. Structural analysis identified a glutamate residue hydrogen-bonded to the heme proximal histidine that may stabilize reactive heme species. A solvent-exposed arginine residue likely gates substrate entry to a tightly confined distal pocket. On the basis of the proposed mechanism of Cld, initial reaction of ClO₂ ⁻ within the distal pocket generates hypochlorite (ClO⁻) and a compound I intermediate. The sterically restrictive distal pocket probably facilitates the rapid rebound of ClO⁻ with compound I forming the Cl⁻ and O₂ products. Common to other heme enzymes, Cld is inactivated after a finite number of turnovers, potentially via the observed formation of an off-pathway tryptophanyl radical species through electron migration to compound I. Three tryptophan residues of Cld have been identified as candidates for this off-pathway radical. Finally, a juxtaposition of hydrophobic residues between the distal pocket and the enzyme surface suggests O₂ may have a preferential direction for exiting the active site.     

Nitric oxide myoglobin: Crystal structure and analysis of ligand geometry

The structure of the ferrous nitric oxide form of native sperm whale myoglobin has been determined by X-ray crystallography to 1.7 Å resolution. The nitric oxide ligand is bent with respect to the heme plane: the Fe-N-O angle is 112°. This angle is smaller than those observed in model compounds and in lupin leghemoglobin. The exact angle appears to be influenced by the strength of the proximal bond and hydrogen bonding interactions between the distal histidine and the bound ligand. Specifically, the N_ϵ atom of histidine⁶⁴ is located 2.8 Å away from the nitrogen atom of the bound ligand, implying electrostatic stabilization of the FeNO complex. This interpretation is supported by mutagenesis studies. When histidine⁶⁴ is replaced with apolar amino acids, the rate of nitric oxide dissociation from myoglobin increases tenfold. Proteins 30:352–356, 1998. © 1998 Wiley-Liss, Inc.     

Hemoprotein models based on a covalent helix-heme-helix sandwich. 3. Coordination properties, reactivity and catalytic application of Fe(III)- and Fe(II)-mimochrome I

 The coordination state of Fe(III)- and Fe(II)-mimochrome I, a covalent peptide-deuteroheme sandwich involving two nonapeptides bearing a histidine residue in a central position, was studied by UV-visible, EPR, and resonance Raman spectroscopy. The ferric and ferrous states of this new iron species mainly exist, at pH 7, in a low-spin hexacoordinate form with two axial histidine ligands coming from the peptide chains. A minor amount of high-spin form for the ferric state is also present at pH 7. However, it is mainly high-spin at pH 2 or in DMSO. Fe(II)-mimochrome I binds CO with an affinity comparable to that of myoglobin and hemoglobin. Fe(III)-mimochrome I reacts with alkylhydroxylamine and arylhydrazines, leading to the corresponding Fe(II)-nitrosoalkyl and Fe(III)-σ-aryl complexes, respectively. These reactions were greatly dependent on the solvent used and on the pH, and were much slower than the corresponding reactions performed by deuterohemin in the presence of excess imidazole. All these results indicate that the reactivity of iron-mimochrome I is controlled by the binding of the peptide chains to the iron. The reactivity shown by this complex at neutral pH is intermediate between that observed for iron porphyrins in the presence of excess imidazole and that of hemoproteins characterized by a strong bis-histidine axial coordination, such as cytochrome b ₅. Fe(III)-mimochrome I is able to catalyze styrene epoxidation by using a [Fe(III)-mimochrome I]/[H₂O₂]/[stryrene] ratio of 1 : 10 : 2000 in phosphate buffer solution (pH 7.2) containing 2% CTAB both under strictly anaerobic conditions and in the presence of oxygen, at 0  °C. Received: 26 May 1998 / Accepted: 20 August 1998     

Spectral properties unique to the myoglobins lacking the usual distal histidine residue

Myoglobins can be divided into two groups. One group contains the usual myoglobins that have histidine at the distal (E7) position, and the other contains a few, but interesting myoglobins that lack the usual distal histidine residue. Spectroscopic examinations have shown that there is a remarkable difference in the Soret band between the two types of myoglobin, and an absorbance ratio of the Soret peak of the acidic met-form to that of the oxy-form seems to be very useful as a simple criterion for predicting whether or not a myoglobin has the usual distal histidine residue.     

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.