Structural dynamics of myoglobin: effect of internal cavities on ligand migration and binding

Using Fourier transform infrared (FTIR) spectroscopy combined with temperature derivative spectroscopy (TDS) at cryogenic temperatures, we have studied CO binding to the heme and CO migration among cavities in the interior of sperm whale carbonmonoxy myoglobin (MbCO) after photodissociation. Photoproduct intermediates, characterized by CO in different locations, were selectively enhanced by laser illumination at specific temperatures. Measurements were performed on the wild-type protein and a series of mutants (L104W, I107W, I28F, and I28W) in which bulky amino acid side chains were introduced to block passageways between cavities or to fill these sites. Binding of xenon was also employed as an alternative means of filling cavities. In all samples, photolyzed CO ligands were observed to initially bind at primary docking site B in the vicinity of the heme iron, from where they migrate to the secondary docking sites, the Xe4 and/or Xe1 cavities. To examine the relevance of these internal docking sites for physiological ligand binding, we have performed room-temperature flash photolysis on the entire set of proteins in the CO- and O(2)-bound form. Together with the cryospectroscopic results, these data provide a clear picture of the role of the internal sites for ligand escape from and binding to myoglobin.     

Mapping the pathways for O2 entry into and exit from myoglobin

The effects of mutagenesis on geminate and bimolecular O2 rebinding to 90 mutants at 27 different positions were used to map pathways for ligand movement into and out of sperm whale myoglobin. By analogy to a baseball glove, the protein "catches" and then "holds" incoming ligand molecules long enough to allow bond formation with the iron atom. Opening of the glove occurs by outward movements of the distal histidine (His(64)), and the ligands are trapped in the interior "webbing" of the distal pocket, in the space surrounded by Ile(28), Leu(29), Leu(32), Val(68), and Ile(107). The size of this pocket is a major determinant of the rate of ligand entry into the protein. Immediately after photo- or thermal dissociation, O2 moves away from the iron into this interior pocket. The majority of the dissociated ligands return to the active site and either rebind to the iron atom or escape through the His(64) gate. A fraction of the ligands migrate further away from the heme group into cavities that have been defined as Xe binding sites 4 and 1; however, most of these ligands also return to the distal pocket, and net escape through the interior of wild-type myoglobin is <20-25%.     

Modulating distal cavities in the α and β subunits of human HbA reveals the primary ligand migration pathway

The free volume in the active site of human HbA plays a crucial role in governing the bimolecular rates of O(2), CO, and NO binding, the fraction of geminate ligand recombination, and the rate of NO dioxygenation by the oxygenated complex. We have decreased the size of the distal pocket by mutating Leu(B10), Val(E11), and Leu(G8) to Phe and Trp and that of other more internal cavities by filling them with Xe at high gas pressures. Increasing the size of the B10 side chain reduces bimolecular rates of ligand binding nearly 5000-fold and inhibits CO geminate recombination due to both reduction of the capture volume in the distal pocket and direct steric hindrance of Fe-ligand bond formation. Phe and Trp(E11) mutations also cause a decrease in distal pocket volume but, at the same time, increase access to the Fe atom because of the loss of the γ2 CH(3) group of the native Val(E11) side chain. The net result of these E11 substitutions is a dramatic increase in the rate of geminate recombination because dissociated CO is sequestered close to the Fe atom and can rapidly rebind without steric resistance. However, the bimolecular rate constants for binding of ligand to the Phe and Trp(E11) mutants are decreased 5-30-fold, because of a smaller capture volume. Geminate and bimolecular kinetic parameters for Phe and Trp(G8) mutants are similar to those for the native HbA subunits because the aromatic rings at this position cause little change in distal pocket volume and because ligands do not move past this position into the globin interior of wild-type HbA subunits. The latter conclusion is verified by the observation that Xe binding to the α and β Hb subunits has little effect on either geminate or bimolecular ligand rebinding. All of these experimental results argue strongly against alternative ligand migration pathways that involve movements through the protein interior in HbA. Instead, ligands appear to enter through the His(E7) gate and are captured directly in the distal cavity.     

A hierarchy of functionally important relaxations within myoglobin based on solvent effects, mutations and kinetic model

Geminate CO rebinding in myoglobin is studied for two viscous solvents, trehalose and sol-gel (bathed in 100% glycerol) at several temperatures. Mutations in key distal hemepocket residues are used to eliminate or enhance specific relaxation modes. The time-resolved data are analyzed with a modified Agmon-Hopfield model which is capable of providing excellent fits in cases where a single relaxation mode is dominant. Using this approach, we determine the relaxation rate constants of specific functionally important modes, obtaining also their Arrhenius activation energies. We find a hierarchy of distal pocket modes controlling the rebinding kinetics. The "heme access mode" (HAM) is responsible for the major slow-down in rebinding. It is a solvent-coupled cooperative mode which restricts ligand return from the xenon cavities. Bulky side-chains, like those His64 and Trp29 (in the L29W mutant), operate like overdamped pendulums which move over and block the binding site. They may be either unslaved (His64) or moderately slaved (Trp29) to the solvent. Small side-chain relaxations, most notably of leucines, are revealed in some mutants (V68L, V68A). They are conjectured to facilitate inter-cavity ligand motion. When all relaxations are arrested (H64L in trehalose), we observe pure inhomogeneous kinetics with no temperature dependence, suggesting that proximal relaxation is not a factor on the investigated timescale.     

The position 68(E11) side chain in myoglobin regulates ligand capture, bond formation with heme iron, and internal movement into the xenon cavities

After photodissociation, ligand rebinding to myoglobin exhibits complex kinetic patterns associated with multiple first-order geminate recombination processes occurring within the protein and a simpler bimolecular phase representing second-order ligand rebinding from the solvent. A smooth transition from cryogenic-like to solution phase properties can be obtained by using a combination of sol-gel encapsulation, addition of glycerol as a bathing medium, and temperature tuning (-15 --> 65 degrees C). This approach was applied to a series of double mutants, myoglobin CO (H64L/V68X, where X = Ala, Val, Leu, Asn, and Phe), which were designed to examine the contributions of the position 68(E11) side chain to the appearance and disappearance of internal rebinding phases in the absence of steric and polar interactions with the distal histidine. Based on the effects of viscosity, temperature, and the stereochemistry of the E11 side chain, the three major phases, B --> A, C --> A, and D --> A, can be assigned, respectively, to ligand rebinding from the following: (i) the distal heme pocket, (ii) the xenon cavities prior to large amplitude side chain conformational relaxation, and (iii) the xenon cavities after significant conformational relaxation of the position 68(E11) side chain. The relative amplitudes of the B --> A and C --> A phases depend markedly on the size and shape of the E11 side chain, which regulates sterically both ligand return to the heme iron atom and ligand migration to the xenon cavities. The internal xenon cavities provide a transient docking site that allows side chain relaxations and the entry of water into the vacated distal pocket, which in turn slows ligand recombination markedly.     

Extended molecular dynamics simulation of the carbon monoxide migration in sperm whale myoglobin

We report the results of an extended molecular dynamics simulation on the migration of photodissociated carbon monoxide in wild-type sperm whale myoglobin. Our results allow following one possible ligand migration dynamics from the distal pocket to the Xe1 cavity via a path involving the other xenon binding cavities and momentarily two additional packing defects along the pathway. Comparison with recent time resolved structural data obtained by Laue crystallography with subnanosecond to millisecond resolution shows a more than satisfactory agreement. In fact, according to time resolved crystallography, CO, after photolysis, can occupy the Xe1 and Xe4 cavities. However, no information on the trajectory of the ligand from the distal pocket to the Xe1 is available. Our results clearly show one possible path within the protein. In addition, although our data refer to a single trajectory, the local dynamics of the ligand in each cavity is sufficiently equilibrated to obtain local structural and thermodynamic information not accessible to crystallography. In particular, we show that the CO motion and the protein fluctuations are strictly correlated: free energy calculations of the migration between adjacent cavities show that the migration is not a simple diffusion but is kinetically or thermodynamically driven by the collective motions of the protein; conversely, the protein fluctuations are influenced by the ligand in such a way that the opening/closure of the passage between adjacent cavities is strictly correlated to the presence of CO in its proximity. The compatibility between time resolved crystallographic experiments and molecular dynamics simulations paves the way to a deeper understanding of the role of internal dynamics and packing defects in the control of ligand binding in heme proteins.     

Structural dynamics of myoglobin: ligand migration among protein cavities studied by Fourier transform infrared/temperature derivative spectroscopy.

Fourier transform infrared (FTIR) spectroscopy in the CO stretch bands combined with temperature derivative spectroscopy (TDS) was used to characterize intermediate states obtained by photolysis of two sperm whale mutant myoglobins, YQR (L29(B10)Y, H64(E7)Q, T67(E10)R) and YQRF (with an additional I107(G8)F replacement). Both mutants assume two different bound-state conformations, A(0) and A(3), which can be distinguished by their different CO bands near 1965 and 1933 cm(-1). They most likely originate from different conformations of the Gln-64 side chain. Within each A substate, a number of photoproduct states have been characterized on the basis of the temperature dependence of recombination in TDS experiments. Different locations and orientations of the ligand within the protein can be distinguished by the infrared spectra of the photolyzed CO. Recombination from the primary docking site, B, near the heme dominates below 50 K. Above 60 K, ligand rebinding occurs predominantly from a secondary docking site, C', in which the CO is trapped in the Xe4 cavity on the distal side, as shown by crystallography of photolyzed YQR and L29W myoglobin CO. Another kinetic state (C") has been identified from which rebinding occurs around 130 K. Moreover, a population appearing above the solvent glass transition at approximately 180 K (D state) is assigned to rebinding from the Xe1 cavity, as suggested by the photoproduct structure of the L29W sperm whale myoglobin mutant. For both the YQR and YQRF mutants, rebinding from the B sites near the heme differs for the two A substates, supporting the view that the return of the ligand from the C', C", and D states is not governed by the recombination barrier at the heme iron but rather by migration to the active site. Comparison of YQR and YQRF shows that access to the Xe4 site (C') is severely restricted by introduction of the bulky Phe side chain at position 107.     

Ligand migration and binding in the dimeric hemoglobin of Scapharca inaequivalvis

Using Fourier transform infrared (FTIR) spectroscopy combined with temperature derivative spectroscopy (TDS) at cryogenic temperatures, we have studied CO binding to the heme and CO migration among cavities in the interior of the dimeric hemoglobin of Scapharca inaequivalvis (HbI) after photodissociation. By combining these studies with X-ray crystallography, three transient ligand docking sites were identified: a primary docking site B in close vicinity to the heme iron, and two secondary docking sites C and D corresponding to the Xe4 and Xe2 cavities of myoglobin. To assess the relevance of these findings for physiological binding, we also performed flash photolysis experiments on HbICO at room temperature and equilibrium binding studies with dioxygen. Our results show that the Xe4 and Xe2 cavities serve as transient docking sites for unbound ligands in the protein, but not as way stations on the entry/exit pathway. For HbI, the so-called histidine gate mechanism proposed for other globins appears as a plausible entry/exit route as well.     

Protein-ligand dynamics. A 96 picosecond simulation of a myoglobin-xenon complex

A 96 picosecond dynamics trajectory of myoglobin with five xenon-probe ligands in internal cavities is examined to study the effect of protein motions on ligand motion and internal cavity fluctuations. Average structural and energetic properties indicate that the simulation is well behaved. The average protein volume is similar to the volume of the X-ray model and the main-chain atom root-mean-square deviation between the X-ray model and the average dynamical structure is 1.25 A. The protein volume oscillates 3 to 4% around the volume of the X-ray structure. These fluctuations lead to changes in the internal free volume and in the size, shape and location of atom-sized cavity features. Transient cavities produced in the simulation have a crucial role in the movement of two of the ligands. One of the ligands escapes to the protein surface, whilst a second ligand travels through the protein interior. Complex gating processes involving several protein residues are responsible for producing the necessary pores through which the ligand passes between transient cavities or packing defects.     

Controlling ligand binding in myoglobin by mutagenesis.

A quadruple mutant of sperm whale myoglobin was constructed to mimic the structure found in Ascaris suum hemoglobin. The replacements include His(E7)-->Gln, Leu(B10)-->Tyr, Thr(E10)--> Arg, and Ile(G8)-->Phe. Single, double, and triple mutants were characterized to dissect out the effects of the individual substitutions. The crystal structures of the deoxy and oxy forms of the quadruple mutant were determined and compared with that of native Ascaris hemoglobin. Tyr(B10) myoglobin displays low O(2) affinity, high dissociation rate constants, and heterogeneous kinetic behavior, suggesting unfavorable steric interactions between the B10 phenol side chain and His(E7). In contrast, all mutants containing the Tyr(B10)/Gln(E7) pair show high O(2) affinity, low dissociation rate constants, and simple, monophasic kinetic behavior. Replacement of Ile(107) with Phe enhances nanosecond geminate recombination singly and in combination with the Tyr(B10)/Gln(E7)/Arg(E10) mutation by limiting access to the Xe4 site. These kinetic results and comparisons with native Ascaris hemoglobin demonstrate the importance of distal pocket cavities in governing the kinetics of ligand binding. The approximately 150-fold higher O(2) affinity of Ascaris hemoglobin compared with that for Tyr(B10)/Gln(E7)-containing myoglobin mutants appears to be the result of favorable proximal effects in the Ascaris protein, due to a staggered orientation of His(F8), the lack of a hydrogen bonding lattice between the F4, F7, and F8 residues, and the presence of a large polar Trp(G5) residue in the interior portion of the proximal heme pocket.     

Structural dynamics of myoglobin: an infrared kinetic study of ligand migration in mutants YQR and YQRF

Recombination of carbon monoxide to myoglobin mutants YQR and YQRF was studied using transient infrared absorption spectroscopy and Fourier transform infrared-temperature derivative spectroscopy (FTIR-TDS). Photoproduct states B, C', C" and D associated with ligands residing in different protein cavities have been identified. After photolysis, ligands migrate to primary docking site B and subsequently rebind or escape to a secondary site (C) within the Xe4 cavity. For YQR, a global analysis of the isothermal rebinding kinetics below 160 K and the TDS data reveal a correlation between the enthalpy barriers governing the two processes. Above 120 K, a protein conformational change in both YQR and YQRF converts photoproduct C' into C" with markedly slowed kinetics. Above approximately 180 K, ligands migrate to the proximal Xe1 site (D) and also exit into the solvent, from where they rebind in a bimolecular reaction.     

Ligand migration and protein fluctuations in myoglobin mutant L29W

We have determined eight X-ray structures of myoglobin mutant L29W at various experimental conditions. In addition, infrared spectroscopic experiments are presented, which are discussed in the light of the X-ray structures. Two distinct conformations of the CO-ligated protein were identified, giving rise to two stretching bands of heme-bound CO. If L29W MbCO crystals are illuminated around 180 K, a deoxy species is formed. The CO molecules migrate to the proximal side of the heme and remain trapped in the so-called Xe1 cavity upon temperature decrease to 105 K. The structure of this photoproduct is almost identical to the equilibrium high-temperature deoxy Mb structure. If the temperature is cycled to increasingly higher values, CO recombination is observed. Three intermediate structures have been determined during the rebinding process. Efficient recombination occurs only above 180 K, the characteristic temperature for the onset of protein dynamics. Rebinding is remarkably slow because bulky residues His64 and Trp29 block important migration pathways of the CO molecule.     

Inversion of axial coordination in myoglobin to create a "proximal" ligand binding pocket

A ligand binding pocket has been created on the proximal side of the heme in porcine myoglobin by site-directed mutagenesis. Our starting point was the H64V/V68H double mutant which has been shown to have bis-histidine (His68 and His93) heme coordination [Dou, Y., Admiraal, S. J., Ikeda-Saito, M., Krzywda, S., Wilkinson, A. J., Li, T., Olson, J. S., Prince, R. C., Pickering, I. J., George, G. N. (1995) J. Biol. Chem. 270, 15993-16001]. The replacement of the proximal His93 ligand by noncoordinating Ala (H64V/V68H/H93A) or Gly (H64V/V68H/H93G) residues resulted unexpectedly in a six-coordinate low-spin species in both ferric and ferrous states. To test the hypothesis that the sixth coordinating ligand in the triple mutants was the imidazole of His97, this residue was mutated to Phe, in the quadruple mutants, H64V/V68H/H93A/H97F and H64V/V68H/H93G/H97F. The ferric quadruple mutants show a clear water/hydroxide alkaline transition and high cyanide and CO affinities, characteristics similar to those of wild-type myoglobin. The nu(Fe-CO) and nu(C-O) stretching frequencies in the ferrous-CO state of the quadruple mutants indicate that the "proximal" ligand binding heme pocket is less polar than the distal pocket in the wild-type protein. Thus, we conclude that the proximal heme pocket in the quadruple mutants has a similar affinity for exogenous ligands to the distal pocket of wild-type myoglobin but that the two pockets have different polarities. The quadruple mutants open up new approaches for developing heme chemistry on the myoglobin scaffold.     

Myoglobin cavities provide interior ligand pathway

The myoglobin protein binds oxygen and catalyzes NO oxidation. As a key model protein, its dynamics have been well studied by spectroscopy and by crystallography as well as by simulation. Nonetheless, visualization of the mechanism of movement of ligands within myoglobin has been difficult. Coordinates of the A1 and A3 taxonomic spectral states of myoglobin from the 1 A crystal structure (1a6g) are generated as consistent sets of correlated clusters of residues with A or B crystal alternates. Analysis of cavities in these A1 and A3 conformations clarifies the pathway of ligand motion from distal entry through interior movement to the proximal side of the heme. Cavities opened up by buried alternate conformations link the distal to the proximal side of the heme. Structural conservation highlights the relevance of this pathway to human neuroglobin. Cavity migration via myoglobin crystal alternates provides a specific link of protein structure to protein dynamics and protein function and demonstrates the relevance of substates (discrete disorder) to function for all proteins.     

Molecular dynamics simulation of carboxy and deoxy human cytoglobin in solution

Cytoglobin (Cgb), the fourth member of the vertebrate heme globin family, is widely expressed in mammalian tissues, and reversibly binds to CO, O₂ and other small ligands. The diverse functions of Cgb may include ligand transport, redox reactions and enzymatic catalysis. Recent studies indicate that Cgb is a potential gene medicine for fibrosis and cancer therapy. In the present work, molecular dynamics (MD) simulations were performed to investigate the functionally related structural properties and dynamic characteristics in carboxy and deoxy human Cgb. The simulation results showed that the loop regions and internal cavities were significantly affected through the binding of an exogenous ligand. The AB, GH and EF loops were found to undergo significant rearrangement and this led to distinct cavity adjustments in Xe2, Xe4 and the distal pocket. In addition, solvent accessibility and torsion angle analyses revealed an interactive distal network comprised of His⁸¹(E7), Leu⁴⁶(B10) and Arg⁸⁴(E10). The MD study of carboxy and deoxy human Cgb revealed that CO-ligated Cgb modulates the protein conformation primarily by loop and cavity rearrangements rather than the heme sliding mechanism found in neuroglobin (Ngb). The significant differences between Cgb and Ngb in the loop and cavity properties are presumably linked to their various biological functions.     

Structural dynamics of myoglobin: spectroscopic and structural characterization of ligand docking sites in myoglobin mutant L29W

We have studied CO binding to the heme and CO migration among protein internal cavities after photodissociation in sperm whale carbonmonoxy myoglobin (MbCO) mutant L29W using Fourier transform infrared (FTIR) spectroscopy combined with temperature derivative spectroscopy (TDS) and kinetic experiments at cryogenic temperatures. Photoproduct intermediates, characterized by CO at particular locations in the protein, were selectively enhanced by applying special laser illumination protocols. These studies were performed on the L29W mutant protein and a series of double mutants constructed so that bulky amino acid side chains block passageways between cavities or fill these sites. Binding of xenon was also employed as an alternative means of occluding cavities. All mutants exhibit two conformations, A(I) and A(II), with distinctly different photoproduct states and ligand binding properties. These differences arise mainly from different positions of the W29 and H64 side chains in the distal heme pocket [Ostermann, A., et al. (2000) Nature 404, 205-208]. The detailed knowledge of the interplay between protein structure, protein dynamics, and ligand migration at cryogenic temperatures allowed us to develop a dynamic model that explains the slow CO and O(2) bimolecular association observed after flash photolysis at ambient temperature.     

The leghemoglobin proximal heme pocket directs oxygen dissociation and stabilizes bound heme

Sperm whale myoglobin (Mb) and soybean leghemoglobin (Lba) are two small, monomeric hemoglobins that share a common globin fold but differ widely in many other aspects. Lba has a much higher affinity for most ligands, and the two proteins use different distal and proximal heme pocket regulatory mechanisms to control ligand binding. Removal of the constraint provided by covalent attachment of the proximal histidine to the F-helices of these proteins decreases oxygen affinity in Lba and increases oxygen affinity in Mb, mainly because of changes in oxygen dissociation rate constants. Hence, Mb and Lba use covalent constraints in opposite ways to regulate ligand binding. Swapping the F-helices of the two proteins brings about similar effects, highlighting the importance of this helix in proximal heme pocket regulation of ligand binding. The F7 residue in Mb is capable of weaving a hydrogen-bonding network that holds the proximal histidine in a fixed orientation. On the contrary, the F7 residue in Lba lacks this property and allows the proximal histidine to assume a conformation favorable for higher ligand binding affinity. Geminate recombination studies indicate that heme iron reactivity on picosecond timescales is not the dominant cause for the effects observed in each mutation. Results also indicate that in Lba the proximal and distal pocket mutations probably influence ligand binding independently. These results are discussed in the context of current hypotheses for proximal heme pocket structure and function.     

A 1H NMR comparison of the met-cyano complexes of elephant and sperm whale myoglobin. Assignment of labile proton resonances in the heme cavity and determination of the distal glutamine orientation from relaxation data

The met-cyano complex of elephant myoglobin has been investigated by high field 1H NMR spectroscopy, with special emphasis on the use of exchangeable proton resonances in the heme cavity to obtain structural information on the distal glutamine. Analysis of the distance dependence of relaxation rates and the exchange behavior of the four hyperfine shifted labile proton resonances has led to the assignment of the proximal His-F8 ring and peptide NHs and the His-FG3 ring NH and the distal Gln-E7 amide NH. The similar hyperfine shift patterns for both the apparent heme resonances as well as the labile proton peaks of conserved resonances in elephant and sperm whale met-cyano myoglobins support very similar electronic/molecular structures for their heme cavities. The essentially identical dipolar shifts and dipolar relaxation times for the distal Gln-E7 side chain NH and the distal His-E7 ring NH in sperm whale myoglobin indicate that those labile protons occupy the same geometrical position relative to the iron and heme plane. This geometry is consistent with the distal residue hydrogen bonding to the coordinated ligand. The similar rates and identical mechanisms of exchange with bulk water of the labile protons for the three conserved residues in the elephant and sperm whale heme cavity indicate that the dynamic stability of the proximal side of the heme pocket is unaltered upon the substitution (His----Gln). The much slower exchange rate (by greater than 10(4] of the distal NH in elephant relative to sperm whale myoglobin supports the assignment of the resonance to the intrinsically less labile amide side chain.     

Crystal structure of cytoglobin: the fourth globin type discovered in man displays heme hexa-coordination

Cytoglobin is a recently discovered hemeprotein belonging to the globin superfamily together with hemoglobin, myoglobin and neuroglobin. Although distributed in almost all human tissues, cytoglobin has not been ascribed a specific function. Human cytoglobin is composed of 190 amino acid residues. Sequence alignments show that a protein core region (about 150 residues) is structurally related to hemoglobin and myoglobin, being complemented by about 20 extra residues both on the N and C termini. In the absence of exogenous ligands (e.g. O2), the cytoglobin distal HisE7 residue is coordinated to the heme Fe atom, thus decreasing the ligand affinity. The crystal structure of human cytoglobin (2.1 A resolution, 21.3% R-factor) highlights a three-over-three alpha-helical globin fold, covering residues 18-171; the 1-17 N-terminal, and the 172-190 C-terminal residue segments are disordered in both molecules of the crystal asymmetric unit. Heme hexa-coordination is evident in one of the two cytoglobin chains, whereas alternate conformation for the heme distal region, achieving partial heme penta-coordination, is observed in the other. Human cytoglobin displays a large apolar protein matrix cavity, next to the heme, not related to the myoglobin cavities recognized as temporary ligand docking stations. The cavity, which may provide a heme ligand diffusion pathway, is connected to the external space through a narrow tunnel nestled between the globin G and H helices.     

H93G myoglobin cavity mutant as versatile template for modeling heme proteins: magnetic circular dichroism studies of thiolate- and imidazole-ligated complexes

Recent ligand binding and spectroscopic investigations of the myoglobin H93G cavity mutant are reviewed, revealing it to be a versatile template for the preparation of model heme complexes of defined structure. The H93G myoglobin cavity mutant is shown to be capable of forming mixed ligand adducts because of the difference in accessibility of the two sides of the ferric heme iron. With imidazole bound in the proximal cavity, H93G myoglobin also forms reasonably stable oxyferrous and oxoferryl derivatives, thereby providing a potential system to use for the study of such complexes with proximal ligands other than imidazole. In addition, thiolate-ligated ferric H93G derivatives are described that serve as spectroscopic models for the high-spin ferric state of cytochrome P450. All of the complexes described are characterized with magnetic circular dichroism spectroscopy, and they are compared to the appropriate derivatives of native myoglobin and P450.