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.     

Solution 1H NMR study of the active site structure for the double mutant H64Q/V68F cyanide complex from mouse neuroglobin

Solution (1)H NMR spectroscopy has been carried out to investigate the molecular and electronic structures of the active site in H64Q/V68F double mutant mouse neuroglobin in the cyanomet form. Two heme orientations resulting from a 180 degrees rotation about the alpha-gamma-meso axis were observed with a population ratio about 1:1, and the clearly distinguished B isomer was used to perform the study. Based on the analysis of the dipolar shifts and paramagnetic relaxation constants, the distal Gln(64)(E7) side chain is obtained to adopt an orientation that may produce hydrogen bond between the N(epsilon)H(1) and the Fe-bound cyanide. The side chain of Phe(68)(E11) is oriented out of the heme pocket just like that in triple mutant of cyanide complex of sperm whale myoglobin. A 15 degrees rotation of the imidazole ring in axial His(96) is observed, which is close to the varphi angle determined from the crystal structure of NgbCO. The quantitative determinations of the orientation and anisotropies of the paramagnetic susceptibility tensor reveal that cyanide is tilted by 8 degrees from the heme normal which allows for contact to the Gln(64)(E7) N(epsilon)H(1). The E7 and E11 residues appear to control the direction and the extent of tilt of the bound ligand. Furthermore, the tilt of the ligand has no obvious influence on the heme heterogeneity of cyanide ligation for isomer A/B of the wild type and mutant protein, indicating that factors other than steric effects, such as polarity of heme pocket, impacts on ligand binding affinity.     

Role of hydrophilic residues in proton transfer during catalysis by human carbonic anhydrase II

Catalysis by the zinc metalloenzyme human carbonic anhydrase II (HCA II) is limited in maximal velocity by proton transfer between His64 and the zinc-bound solvent molecule. Asn62 extends into the active site cavity of HCA II adjacent to His64 and has been shown to be one of several hydrophilic residues participating in a hydrogen-bonded solvent network within the active site. We compared several site-specific mutants of HCA II with replacements at position 62 (Ala, Val, Leu, Thr, and Asp). The efficiency of catalysis in the hydration of CO 2 for the resulting mutants has been characterized by (18)O exchange, and the structures of the mutants have been determined by X-ray crystallography to 1.5-1.7 A resolution. Each of these mutants maintained the ordered water structure observed by X-ray crystallography in the active site cavity of wild-type HCA II; hence, this water structure was not a variable in comparing with wild type the activities of mutants at residue 62. Crystal structures of wild-type and N62T HCA II showed both an inward and outward orientation of the side chain of His64; however, other mutants in this study showed predominantly inward (N62A, N62V, N62L) or predominantly outward (N62D) orientations of His64. A significant role of Asn62 in HCA II is to permit two conformations of the side chain of His64, the inward and outward, that contributes to maximal efficiency of proton transfer between the active site and solution. The site-specific mutant N62D had a mainly outward orientation of His64, yet the difference in p K a between the proton donor His64 and zinc-bound hydroxide was near zero, as in wild-type HCA II. The rate of proton transfer in catalysis by N62D HCA II was 5% that of wild type, showing that His64 mainly in the outward orientation is associated with inefficient proton transfer compared with His64 in wild type which shows both inward and outward orientations. These results emphasize the roles of the residues of the hydrophilic side of the active site cavity in maintaining efficient catalysis by carbonic anhydrase.     

Tautomerism,acid-base equilibria,and H-bonding of the six histidines in subtilisin BPN' by NMR

We have determined by (15)N, (1)H, and (13)C NMR, the chemical behavior of the six histidines in subtilisin BPN' and their PMSF and peptide boronic acid complexes in aqueous solution as a function of pH in the range of from 5 to 11, and have assigned every (15)N, (1)H, C(epsilon 1), and C(delta2) resonance of all His side chains in resting enzyme. Four of the six histidine residues (17, 39, 67, and 226) are neutrally charged and do not titrate. One histidine (238), located on the protein surface, titrates with pK(a) = 7.30 +/- 0.03 at 25 degrees C, having rapid proton exchange, but restricted mobility. The active site histidine (64) in mutant N155A titrates with a pK(a) value of 7.9 +/- 0.3 and sluggish proton exchange behavior, as shown by two-site exchange computer lineshape simulation. His 64 in resting enzyme contains an extremely high C(epsilon 1)-H proton chemical shift of 9.30 parts per million (ppm) owing to a conserved C(epsilon 1)-H(.)O=C H-bond from the active site imidazole to a backbone carbonyl group, which is found in all known serine proteases representing all four superfamilies. Only His 226, and His 64 at high pH, exist as the rare N(delta1)-H tautomer, exhibiting (13)C(delta1) chemical shifts approximately 9 ppm higher than those for N(epsilon 2)-H tautomers. His 64 in the PMSF complex, unlike that in the resting enzyme, is highly mobile in its low pH form, as shown by (15)N-(1)H NOE effects, and titrates with rapid proton exchange kinetics linked to a pK(a) value of 7.47 +/- 0.02.     

Influence of the heme pocket conformation on the structure and vibrations of the Fe-CO bond in myoglobin: a QM/MM density functional study.

The influence of the distal pocket conformation on the structure and vibrations of the heme-CO bond in carbonmonoxy myoglobin (MbCO) is investigated by means of hybrid QM/MM calculations based on density functional theory combined with a classical force field. It is shown that the heme-CO structure (QM treated) is quite rigid and not influenced by the distal pocket conformation (MM treated). This excludes any relation between FeCO distortions and the different CO absorptions observed in the infrared spectra of MbCO (A states). In contrast, both the CO stretch frequency and the strength of the CO...His64 interaction are very dependent on the orientation and tautomerization state of His64. Our calculations indicate that the CO...N(epsilon) type of approach does not contribute to the A states, whereas the CO...H-N(epsilon) interaction is the origin of the A(1) and A(3) states, the His64 residue being protonated at N(epsilon). The strength of the CO...His64 interaction is quantified, in comparison with the analogous O(2)...His64 interaction and with the observed changes in the CO stretch frequency. Additional aspects of the CO...His64 interaction and its biological implications are discussed.     

Coordination structure of the ferric heme iron in engineered distal histidine myoglobin mutants.

Recombinant human myoglobin mutants with the distal His residue (E7, His64) replaced by Leu, Val, or Gln residues were prepared by site-directed mutagenesis and expression in Escherichia coli. Electronic and coordination structures of the ferric heme iron in the recombinant myoglobin proteins were examined by optical absorption, EPR, 1H NMR, magnetic circular dichroism, and x-ray spectroscopy. Mutations, His-->Val and His-->Leu, remove the heme-bound water molecule resulting in a five-coordinate heme iron at neutral pH, while the heme-bound water molecule appears to be retained in the engineered myoglobin with His-->Gln substitution as in the wild-type protein. The distal Val and distal Leu ferric myoglobin mutants at neutral pH exhibited EPR spectra with g perpendicular values smaller than 6, which could be interpreted as an admixture of intermediate (S = 3/2) and high (S = 5/2) spin states. At alkaline pH, the distal Gln mutant is in the same so-called "hydroxy low spin" form as the wild-type protein, while the distal Leu and distal Val mutants are in high spin states. The ligand binding properties of these recombinant myoglobin proteins were studied by measurements of azide equilibrium and cyanide binding. The distal Leu and distal Val mutants exhibited diminished azide affinity and extremely slow cyanide binding, while the distal Gln mutant showed azide affinity and cyanide association rate constants similar to those of the wild-type protein.     

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%.     

The effects of amino acid substitution at position E7 (residue 64) on the kinetics of ligand binding to sperm whale myoglobin

Association and dissociation rate constants were measured for O2, CO, and alkyl isocyanide binding to a set of genetically engineered sperm whale myoglobins with site-specific mutations at residue 64 (the E7 helical position). Native His was replaced by Gly, Val, Leu, Met, Phe, Gln, Arg, and Asp using the synthetic gene and expression system developed by Springer and Sligar (Springer, B. A., and Sligar, S. G. (1987) Proc. Natl. Acad. Sci. U.S.A. 84, 8961-8965). The His64----Gly substitution produced a sterically unhindered myoglobin that exhibited ligand binding parameters similar to those of chelated protoheme suspended in soap micelles. The order of the association rate constants for isocyanide binding to the mutant myoglobins was Gly64 (approximately 10(7) M-1 s-1) much greater than Val64 approximately Leu64 (approximately 10(6) M-1 s-1) greater than Met64 greater than Phe64 approximately His64 approximately Gln64 (10(5)-10(3) M-1 s-1) and indicates that the barrier to isocyanide entry into the distal pocket is primarily steric in nature. The bimolecular rates of methyl, ethyl, n-propyl, and n-butyl isocyanide binding to the His64----Arg and His64----Asp mutants were abnormally high (1-5 x 10(6) M-1 s-1), suggesting that Arg64 and Asp64 adopt conformations with the charged side chains pointing out toward the solvent creating a less hindered pathway for ligand binding. In contrast to the isocyanide data, the association rate constants for O2 and CO binding exhibited little dependence on the size of the E7 side chain. The values for all the mutants except His64----Gln approached or were larger than those for chelated model heme (i.e. approximately 1 x 10(8) M-1 s-1 for O2 and approximately 1 x 10(7) M-1 s-1 for CO), whereas the corresponding rate parameters for myoglobin containing either Gln64 or His64 were 5- to 10-fold smaller. This result suggests that a major kinetic barrier for O2 and CO binding to native myoglobin may involve disruption of polar interactions between His64 and water molecules found in the distal pocket of deoxymyoglobin. Finally, the rate and equilibrium parameters for O2 and CO binding to the His64----Gln, His64----Val, and His64----Leu mutants were compared to those reported previously for Asian elephant myoglobin (Gln-E7), Aplysia limacina myoglobin (Val-E7), and monomeric Hb II from Glycera dibranchiata (Leu-E7).     

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.     

Structural dynamics in the active site of murine neuroglobin and its effects on ligand binding

We have examined the effects of active site residues on ligand binding to the heme iron of mouse neuroglobin using steady-state and time-resolved visible spectroscopy. Absorption spectra of the native protein, mutants H64L and K67L and double mutant H64L/K67L were recorded for the ferric and ferrous states over a wide pH range (pH 4-11), which allowed us to identify a number of different species with different ligands at the sixth coordination, to characterize their spectroscopic properties, and to determine the pK values of active site residues. In flash photolysis experiments on CO-ligated samples, reaction intermediates and the competition of ligands for the sixth coordination were studied. These data provide insights into structural changes in the active site and the role of the key residues His64 and Lys67. His64 interferes with exogenous ligand access to the heme iron. Lys67 sequesters the distal pocket from the solvent. The heme iron is very reactive, as inferred from the fast ligand binding kinetics and the ability to bind water or hydroxyl ligands to the ferrous heme. Fast bond formation favors geminate rebinding; yet the large fraction of bimolecular rebinding observed in the kinetics implies that ligand escape from the distal pocket is highly efficient. Even slight pH variations cause pronounced changes in the association rate of exogenous ligands near physiological pH, which may be important in functional processes.     

Kinetic analysis of multiple proton shuttles in the active site of human carbonic anhydrase

We have prepared a site-specific mutant of human carbonic anhydrase (HCA) II with histidine residues at positions 7 and 64 in the active site cavity. Using a different isozyme, we have placed histidine residues in HCA III at positions 64 and 67 and in another mutant at positions 64 and 7. Each of these histidine residues can act as a proton transfer group in catalysis when it is the only nonliganding histidine in the active site cavity, except His(7) in HCA III. Using an (18)O exchange method to measure rate constants for intramolecular proton transfer, we have found that inserting two histidine residues into the active site cavity of either isozyme II or III of carbonic anhydrase results in rates of proton transfer to the zinc-bound hydroxide that are antagonistic or suppressive with respect to the corresponding single mutants. The crystal structure of Y7H HCA II, which contains both His(7) and His(64) within the active site cavity, shows the conformation of the side chain of His(64) moved from its position in the wild type and hydrogen-bonded through an intervening water molecule with the side chain of His(7). This suggests a cause of decreased proton transfer in catalysis.     

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.     

Role of His505 in the soluble fumarate reductase from Shewanella frigidimarina

The X-ray structure of the soluble fumarate reductase from Shewanella frigidimarina [Taylor, P., Pealing, S. L., Reid, G. A., Chapman, S. K., and Walkinshaw, M. D. (1999) Nat. Struct. Biol. 6, 1108-1112] clearly shows the presence of an internally bound sodium ion. This sodium ion is coordinated by one solvent water molecule (Wat912) and five backbone carbonyl oxygens from Thr506, Met507, Gly508, Glu534, and Thr536 in what is best described as octahedral geometry (despite the rather long distance from the sodium ion to the backbone oxygen of Met507 (3.1 A)). The water ligand (Wat912) is, in turn, hydrogen bonded to the imidazole ring of His505. This histidine residue is adjacent to His504, a key active-site residue thought to be responsible for the observed pK(a) of the enzyme. Thus, it is possible that His505 may be important in both maintaining the sodium site and in influencing the active site. Here we describe the crystallographic and kinetic characterization of the H505A and H505Y mutant forms of the Shewanella fumarate reductase. The crystal structures of both mutant forms of the enzyme have been solved to 1.8 and 2.0 A resolution, respectively. Both show the presence of the sodium ion in the equivalent position to that found in the wild-type enzyme. The structure of the H505A mutant shows the presence of two water molecules in place of the His505 side-chain which form part of a hydrogen-bonding network with Wat48, a ligand to the sodium ion. The structure of the H505Y mutant shows the hydroxyl group of the tyrosine side-chain hydrogen-bonding to a water molecule which is also a ligand to the sodium ion. Apart from these features, there are no significant structural alterations as a result of either substitution. Both the mutant enzymes are catalytically active but show markedly different pH profiles compared to the wild-type enzyme. At high pH (above 8.5), the wild type and mutant enzymes have very similar activities. However, at low pH (6.0), the H505A mutant enzyme is some 20-fold less active than wild-type. The combined crystallographic and kinetic results suggest that His505 is not essential for sodium binding but does affect catalytic activity perhaps by influencing the pK(a) of the adjacent His504.     

Roles of the heme proximal side residues tryptophan409 and tryptophan421 of neuronal nitric oxide synthase in the electron transfer reaction

Nitric oxide synthase (NOS) has an oxygenase domain with a thiol-coordinated heme active side similar to cytochrome P450. In contrast to cytochrome P450, however, conserved aromatic amino acids are situated in the heme proximal side of NOS. For example, in endothelial NOS (eNOS), the indole-ring nitrogen of Trp180 hydrogen-binds to the thiol of Cys186, the internal axial ligand to the heme. And, the aromatic side chain of Trp192 forms a bridge between this residue and the protein. Trp180 and Trp192 of eNOS correspond to Trp409 and Trp421 of neuronal NOS (nNOS), respectively. In order to understand the roles of the aromatic amino acids in catalysis, we generated Trp409His, Trp409Leu, Trp421His and Trp421Leu mutants of nNOS and determined their catalytic parameters. The Trp409Leu mutant was very poorly expressed in E. coli and was easily denatured during purification procedures. The NO formation activities of the Trp409His and Trp421Leu mutants were 11 and 25 micromol/min per micromol heme, respectively, and are lower than that (44 micromol/min per micromol heme) of the wild type. The activity (46 micromol/min per micromol heme) of the Trp421His mutant was comparable to that of the wild-type enzyme. However, NADPH oxidation rates of Trp421His (230 micromol/min per micromol heme) and Trp421Leu (104 micromol/min per microol heme) in the presence of L-Arg were much larger than those observed for the wild type (65 micromol/min per micromol heme) and the Trp409His mutant (43 micromol/min per micromol heme). The cytochrome c reduction rate of the Trp421His mutant was 6-fold larger than that of the wild type. The heme reduction rate with NADPH for the Trp421His mutant (0.09 min(-1)) was much lower than that (1.0 min(-1)) of the wild type. Taken together, it appears that Trp421 may be involved in inter-domain/inter-subunit electron transfer reactions.     

1H-NMR spectroscopic manifestations of ligand binding to the kringle 4 domain of human plasminogen

Structural aspects of the binding of the linear ligands N alpha-acetyl-L-lysine (AcLys) and epsilon-aminocaproic acid (epsilon ACA) and of the cyclic analogs trans-(aminomethyl)-cyclohexanecarboxylic acid (AMCHA) and p-benzylaminesulfonic acid (BASA) to the intact plasminogen kringle 4 domain have been investigated by 1H-NMR spectroscopy at 300 and 600 MHz. Ligand binding results in consistent shifts of the His-II (His31), Trp-I (Trp25?), Trp-II (Trp62?), Trp-III (Trp72), Tyr-II (Tyr50), and Phe64 ring signals. BASA tends to induce larger shifts than elicited by the aliphatic ligands, most noticeably on Trp-II and on Trp72, suggesting that the ligand aromatic ring interacts with the two indole groups. Trp-II and, to lesser extent, Trp-I interact with an acidic side chain group, in a manner that is blocked by BASA. BASA binding also perturbs Tyr-II (Tyr50), Tyr-III (Tyr41), and Tyr-IV (Tyr74) over a wide pH range and lowers the pKa* of His31 from approximately 4.8 to approximately 4.6. His-III (His33) responds to BASA and AMCHA but is relatively insensitive to the linear ligands. His33 carries a sterically shielded side chain which, in conjunction with Leu46, Trp-I, Tyr50, and Tyr74, participates in structuring the kringle hydrophobic core, contiguous to the binding site. Pronounced shifts are observed for aliphatic resonances stemming from the kringle-bound molecules of AMCHA, AcLys, and epsilon ACA. It is proposed that the lysine-binding site is mostly supported by the loop that extends from Cys51 through Cys71 and that aromatic residues, which include Trp-II, Trp72, and Phe64, play a major role in interacting with the nonpolar segment of the ligand molecule. The binding site also encompasses Tyr50, Tyr74, His31, and His33 although it is not clear the extent to which these residues interact directly with the ligand.     

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.     

Structural Assignment of Spectra by Characterization of Conformational Substates in Bound MbCO

Residue motions of the distal heme pocket and bound CO ligand of carbonmonoxy Myoglobin are studied using a combination of molecular dynamics simulations and quantum chemical methods. Using mixed quantum mechanics/molecular mechanics calculations together with sampling from molecular dynamics simulations (QM/MM(MD)), the experimentally observed spectroscopic A₀ and A₁ substates of the bound CO ligand are assigned to the open and closed conformation of His⁶⁴ and the His_ɛ⁶⁴ tautomer, respectively. Several previously proposed origins of the A₃ substate, including rotamers of the doubly protonated His⁶⁴H⁺ side chain, His⁶⁴H⁺ inside the distal pocket, and cooperative motions with Arg⁴⁵, are investigated with QM/MM(MD). However, the signatures of the calculated infrared spectra do not agree with the experimentally observed ones. For additional insight on this, extensive molecular dynamics simulations are used together with improved electrostatics for the bound ligand. A CO fluctuating charge model is developed to describe the ab initio dipole and quadrupole moments of the bound ligand. CO absorption spectra are then obtained directly from the dynamics simulations. Finally, the electrostatics of the heme pocket is examined in detail in an attempt to determine the structural origins of the observed spectroscopic A-states from MD simulations. However, contrary to related simulations for unbound CO in myoglobin, the shifts and splittings for carbonmonoxy Myoglobin are generally small and difficult to relate to structural change. This suggests that coupling of the CO motion to other degrees of freedom, such as the Fe-CO stretching and bending, is important to correctly describe the dynamics of bound CO in myoglobin.     

An engineered heme–copper center in myoglobin: CO migration and binding

We have investigated CO migration and binding in Cu_BMb, a copper-binding myoglobin double mutant (L29H–F43H), by using Fourier transform infrared spectroscopy and flash photolysis over a wide temperature range. This mutant was originally engineered with the aim to mimic the catalytic site of heme–copper oxidases. Comparison of the wild-type protein Mb and Cu_BMb shows that the copper ion in the distal pocket gives rise to significant effects on ligand binding to the heme iron. In Mb and copper-free Cu_BMb, primary and secondary ligand docking sites are accessible upon photodissociation. In copper-bound Cu_BMb, ligands do not migrate to secondary docking sites but rather coordinate to the copper ion. Ligands entering the heme pocket from the outside normally would not be captured efficiently by the tight distal pocket housing the two additional large imidazole rings. Binding at the Cu ion, however, ensures efficient trapping in Cu_BMb. The Cu ion also restricts the motions of the His64 side chain, which is the entry/exit door for ligand movement into the active site, and this restriction results in enhanced geminate and slow bimolecular CO rebinding. These results support current mechanistic views of ligand binding in hemoglobins and the role of the Cu_B in the active of heme–copper oxidases. This article is part of a Special Issue entitled: Oxygen Binding and Sensing Proteins.     

A novel site-directed mutant of myoglobin with an unusually high O2 affinity and low autooxidation rate.

Mutants of sperm whale myoglobin were constructed at position 29 (B10 in helix notation) to examine the effects of distal pocket size on the rates of ligand binding and autooxidation. Leu29 was replaced with Ala, Val, and Phe using the synthetic gene and Escherichia coli expression system of Springer and Sligar (Springer, B. A., and Sligar, S. G. (1987) Proc. Natl. Acad. Sci. U. S. A. 84, 8961-8965). Structures of the ferric forms of Val29 and Phe29, and the oxy form of Phe29 myoglobin were determined to 1.7 A by x-ray crystallography. The ferric mutant proteins are remarkably isomorphous with the wild type protein except in the immediate vicinity of residue 29. Thus, the protein structure in the distal pocket of myoglobin can accommodate either a large "hole" (i.e. Ala or Val) or a large side chain (i.e. Phe) at position 29 without perturbation of tertiary structure. Phe29 oxymyoglobin is also identical to the native oxy protein in terms of overall structure and interactions between the bound O2 and His64, Val68, Phe43, and Ile107. The distance between the nearest side chain atom of residue 29 and the second atom of the bound oxygen molecule is 3.2 A in the Phe29 protein and 4.9 A in native myoglobin. The equilibrium constants for O2 binding to Ala29, Val29, and Leu29 (native) myoglobin are the same, approximately 1.0 x 10(6) M-1 at 20 degrees C, whereas that for the Phe29 protein is markedly greater, 15 x 10(6) M-1. This increase in affinity is due primarily to a 10-fold decrease in the O2 dissociation rate constant for the Phe29 mutant and appears to be the result of stabilizing interactions between the negative portion of the bound O2 dipole and the partially positive edge of the phenyl ring. Increasing the size of residue 29 causes large decreases in the rate of autooxidation of myoglobin: k(ox) = 0.24, 0.23, 0.055, and 0.005 h-1 for Ala29, Val29, Leu29 (native), and Phe29 myoglobin, respectively, in air at 37 degrees C. Thus, the Leu29----Phe mutation produces a reduced protein that is remarkably stable and is expressed in E. coli as 100% MbO2. The selective pressure to conserve Leu29 at the B10 position probably represents a compromise between reducing the rate of autooxidation and maintaining a large enough O2 dissociation rate constant to allow rapid oxygen release during respiration.