Analysis of the kinetic barriers for ligand binding to sperm whale myoglobin using site-directed mutagenesis and laser photolysis techniques

Time courses for NO, O2, CO, methyl and ethyl isocyanide rebinding to native and mutant sperm whale myoglobins were measured at 20 degrees C following 17-ns and 35-ps laser excitation pulses. His64 (E7) was replaced with Gly, Val, Leu, Phe, and Gln, and Val68 (E11) was replaced with Ala, Ile, and Phe. For both NO and O2, the effective picosecond quantum yield of unliganded geminate intermediates was roughly 0.2 and independent of the amino acids at positions 64 and 68. Geminate recombination of NO was very rapid; 90% rebinding occurred within 0.5-1.0 ns for all of the myoglobins examined; and except for the Gly64 and Ile68 mutants, the fitted recombination rate parameters were little influenced by the size and polarity of the amino acid at position 64 and the size of the residue at position 68. The rates of NO recombination and ligand movement away from the iron atom in the Gly64 mutant increased 3-4-fold relative to native myoglobin. For Ile68 myoglobin, the first geminate rate constant for NO rebinding decreased approximately 6-fold, from 2.3 x 10(10) s-1 for native myoglobin to 3.8 x 10(9) s-1 for the mutant. No picosecond rebinding processes were observed for O2, CO, and isocyanide rebinding to native and mutant myoglobins; all of the observed geminate rate constants were less than or equal to 3 x 10(8) s-1. The rebinding time courses for these ligands were analyzed in terms of a two-step consecutive reaction scheme, with an outer kinetic barrier representing ligand movement into and out of the protein and an inner barrier representing binding to the heme iron atom by ligand occupying the distal portion of the heme pocket. Substitution of apolar amino acids for His64 decreased the absolute free energies of the outer and inner kinetic barriers and the well for non-covalently bound O2 and CO by 1 to 1.5 kcal/mol, regardless of size. In contrast, the His64 to Gln mutation caused little change in the barrier heights for all ligands, showing that the polar nature of His64 inhibits both the bimolecular rate of ligand entry into myoglobin and the unimolecular rate of binding to the iron atom from within the protein. Increasing the size of the position 68(E11) residue in the series Ala to Val (native) to Ile caused little change in the rate of O2 migration into myoglobin or the equilibrium constant for noncovalent binding but did decrease the unimolecular rate for iron-O2 bond formation.(ABSTRACT TRUNCATED AT 400 WORDS)     

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.     

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.     

Ligand binding to a hemoprotein lacking the distal histidine. The myoglobin from aplysia limacina (Val(E7))

The time course of ligand recombination to the myoglobin from Aplysia limacina, which has Val(E7), was measured following photolysis by flashes of 35 ps to 300 ns with a time resolution of 10 ps or 1 ns. CO shows only biomolecular recombination. O2 has a small geminate reaction with a half-time of tens of picoseconds, but no nanosecond geminate reaction. NO has two picosecond relaxations with half-times of 70 ps (15%) and 1 ns (80%) and one nanosecond relaxation with a half-time of 4.6 ns. The biomolecular rates for O2 and NO are the same: 2 x 10(7) M-1 s-1. Methyl and ethyl isonitriles have a geminate reaction with a half-time of 35 ps. Ethyl isonitrile has, in addition, a nanosecond relaxation (25%) with a half-time of 100 ns. t-Butyl isonitrile has four geminate relaxations (10 ps, 35 ps, 1 ns, and 1 microseconds). Analysis of the results suggests much easier movement of ligand between the heme pocket and the exterior than in sperm whale myoglobin (His(E7]. The reactivity of the heme is little different, placing the effect of the differences from sperm whale myoglobin on the distal side of the heme.     

Ligand recombination to the alpha and beta subunits of human hemoglobin

The rebinding of CO, O2, NO, methyl, ethyl, n-propyl, and n-butyl isocyanide to isolated alpha and beta chains and intact hemoglobin at pH 7, 20 degrees C was examined both during and after a 30-ns dye laser pulse. The resultant absorbance changes were analyzed in terms of a linear three-step reaction scheme: Hb + X in equilibrium with C in equilibrium with B in equilibrium with A or HbX, where A is the final bound state, and C and B are geminate states. Rate constants were assigned for each of the transitions in this mechanism using fitting procedures described previously for analyzing ligand rebinding to sperm whale myoglobin at room temperature (Gibson, Q. H., Olson, J. S., McKinnie, R. E., and Rohlfs, R. J. (1986) J. Biol. Chem. 261, 10228-10239). Five major conclusions were obtained. First, initial geminate recombination phases for the NO and O2 complexes of hemoglobin and its isolated subunits exhibit half-times equal to approximately 12 and approximately 440 ps, respectively. These values are in excellent agreement with more direct, picosecond measurements of the geminate recombination of HbNO (Cornelius, P. A., Hochstrasser, R. M., and Steele, A. W. (1983) J. Mol. Biol. 163, 119-128) and HbO2 (Friedman, J. M., Scott, T. W., Fisanick, G. J., Simon, S. R., Findsen, E. W., Ondrias, M. R., and MacDonald, V. W. (1985) Science 229, 187-229) following extremely short laser pulses. Second, the correspondence between our nanosecond measurements and the published picosecond data suggests strongly that the intrinsic photochemical yield of all ferrous, hexacoordinate heme complexes approaches one. Third, the major differences between the isolated alpha and beta chains involve the rate of ligand migration to the solvent, kC----X and the extent of recombination from the second geminate state, C, as measured by the ratio kC----B/kC----X. Fourth, for both isolated chains and intact hemoglobin, the rate and equilibrium constants for the formation of the initial O2 geminate state starting from ligand in the solvent (i.e. kX----B and KX----B) are 5-10 times greater than the corresponding parameters for the formation of the first CO geminate state. Fifth, the rate-limiting step for NO, O2, and isonitrile binding to hemoglobin and its isolated subunits is ligand migration up to the initial geminate state (i.e. kX----B). In the case of CO binding, both migration to state B and iron-ligand bond formation (kB----A) affect the overall, bimolecular association rate constant.     

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.     

Contributions of residue 45(CD3) and heme-6-propionate to the biomolecular and geminate recombination reactions of myoglobin

Overall association and dissociation rate constants were measured at 20 degrees C for O2, CO, and alkyl isocyanide binding to position 45 (CD3) mutants of pig and sperm whale myoglobins and to sperm whale myoglobin reconstituted with protoheme IX dimethyl ester. In pig myoglobin, Lys45(CD3) was replaced with Arg, His, Ser, and Glu; in sperm whale myoglobin, Arg45(CD3) was replaced with Ser and Gly. Intramolecular rebinding of NO, O2, and methyl isocyanide to Arg45, Ser45, Glu45, and Lys45(native) pig myoglobins was measured following 35-ps and 17-ns excitation pulses. The shorter, picosecond laser flash was used to examine ligand recombination from photochemically produced contact pairs, and the longer, nanosecond flash was used to measure the rebinding of ligands farther removed from the iron atom. Mutations at position 45 or esterification of the heme did not change significantly (less than or equal to 2-fold) the overall association rate constants for NO, CO, and O2 binding at room temperature. These data demonstrate unequivocally that Lys(Arg)45 makes little contribution to the outer kinetic barrier for the entry of diatomic gases into the distal pocket of myoglobin, a result that contradicts a variety of previous structural and theoretical interpretations. However, the rates of geminate recombination of NO and O2 and the affinity of myoglobin for O2 were dependent upon the basicity of residue 45. The series of substitutions Arg45, Lys45, Ser45, and Glu45 in pig myoglobin led to a 3-fold decrease in the initial rate for the intramolecular, picosecond rebinding of NO and 4-fold decrease in the geminate rate constant for the nanosecond rebinding of O2. (ABSTRACT TRUNCATED AT 250 WORDS)     

A comparison of the geminate recombination kinetics of several monomeric heme proteins

The geminate rate constants for CO, O2, NO, methyl, ethyl, n-propyl, and n-butyl isocyanide rebinding to soybean leghemoglobin and monomeric component II of Glycera dibranchiata hemoglobin were measured at pH 7, 20 degrees C using a dye laser with a 30-ns square-wave pulse. The results were compared to the corresponding parameters for sperm whale myoglobin and the isolated alpha and beta subunits of human hemoglobin (Olson, J.S., Rohlfs, R.J., and Gibson, Q.H. (1987) J. Biol. Chem., 262, 12930-12938). The rate-limiting step for O2, NO, and isonitrile binding to all five proteins is ligand migration up to the initial geminate state, and the rate of this process determines the overall bimolecular association rate constant for these ligands. In contrast, iron-ligand bond formation limits the overall bimolecular rate for CO binding. The distal pockets in leghemoglobin and in Glycera HbII are approximately 10 times more accessible kinetically to diatomic ligands than that in sperm whale myoglobin. This difference accounts for the much larger association rate constants (1-2 x 10(8) M-1 s-1) that are observed for O2 and NO binding to leghemoglobin and Glycera HbII. The rates of isonitrile migration through leghemoglobin are also very large and indicate a very fluid or open distal structure near the sixth coordination position. In contrast, there is a marked decrease in the rate of migration up to and away from the sixth coordination position in Glycera HbII with increasing ligand size. These results were also used to interpret previously published rate constants and quantum yields for the high (R) and low (T) affinity states of human hemoglobin. In contrast to the differences between the monomeric proteins, the differences between the CO-, O2-, and NO-binding parameters for R and T state hemoglobin appear to be due to a decrease in the geminate reactivity of the heme iron atom, with little or no change in the accessibility of the distal pocket.     

Proximal and distal influences on ligand binding kinetics in microperoxidase and heme model compounds

We use laser flash photolysis and time-resolved Raman spectroscopy of CO-bound heme complexes to study proximal and distal influences on ligand rebinding kinetics. We report kinetics of CO rebinding to microperoxidase (MP) and 2-methylimidazole ligated Fe protoporphyrin IX in the 10 ns to 10 ms time window. We also report CO rebinding kinetics of MP in the 150 fs to 140 ps time window. For dilute, micelle-encapsulated (monodisperse) samples of MP, we do not observe the large amplitude geminate decay at approximately 100 ps previously reported in time-resolved IR measurements on highly concentrated samples [Lim, M., Jackson, T. A., and Anfinrud, P. A. (1997) J. Biol. Inorg. Chem. 2, 531-536]. However, for high concentration aggregated samples, we do observe the large amplitude picosecond CO geminate rebinding and find that it is correlated with the absence of the iron-histidine vibrational mode in the time-resolved Raman spectrum. On the basis of these results, the energetic significance of a putative distal pocket CO docking site proposed by Lim et al. may need to be reconsidered. Finally, when high concentration samples of native myoglobin (Mb) were studied as a control, an analogous increase in the geminate rebinding kinetics was not observed. This verifies that studies of Mb under dilute conditions are applicable to the more concentrated regime found in the cellular milieu.     

The amino acid sequence of human chorionic gonadotropin. The alpha subunit and beta subunit.

The amino acid sequences of both the alpha and beta subunits of human chorionic gonadotropin have been determined. The amino acid sequence of the alpha subunit is: Ala - Asp - Val - Gln - Asp - Cys - Pro - Glu - Cys-10 - Thr - Leu - Gln - Asp - Pro - Phe - Ser - Gln-20 - Pro - Gly - Ala - Pro - Ile - Leu - Gln - Cys - Met - Gly-30 - Cys - Cys - Phe - Ser - Arg - Ala - Tyr - Pro - Thr - Pro-40 - Leu - Arg - Ser - Lys - Lys - Thr - Met - Leu - Val - Gln-50 - Lys - Asn - Val - Thr - Ser - Glu - Ser - Thr - Cys - Cys-60 - Val - Ala - Lys - Ser - Thr - Asn - Arg - Val - Thr - Val-70 - Met - Gly - Gly - Phe - Lys - Val - Glu - Asn - His - Thr-80 - Ala - Cys - His - Cys - Ser - Thr - Cys - Tyr - Tyr - His-90 - Lys - Ser. Oligosaccharide side chains are attached at residues 52 and 78. In the preparations studied approximately 10 and 30% of the chains lack the initial 2 and 3 NH2-terminal residues, respectively. This sequence is almost identical with that of human luteinizing hormone (Sairam, M. R., Papkoff, H., and Li, C. H. (1972) Biochem. Biophys. Res. Commun. 48, 530-537). The amino acid sequence of the beta subunit is: Ser - Lys - Glu - Pro - Leu - Arg - Pro - Arg - Cys - Arg-10 - Pro - Ile - Asn - Ala - Thr - Leu - Ala - Val - Glu - Lys-20 - Glu - Gly - Cys - Pro - Val - Cys - Ile - Thr - Val - Asn-30 - Thr - Thr - Ile - Cys - Ala - Gly - Tyr - Cys - Pro - Thr-40 - Met - Thr - Arg - Val - Leu - Gln - Gly - Val - Leu - Pro-50 - Ala - Leu - Pro - Gin - Val - Val - Cys - Asn - Tyr - Arg-60 - Asp - Val - Arg - Phe - Glu - Ser - Ile - Arg - Leu - Pro-70 - Gly - Cys - Pro - Arg - Gly - Val - Asn - Pro - Val - Val-80 - Ser - Tyr - Ala - Val - Ala - Leu - Ser - Cys - Gln - Cys-90 - Ala - Leu - Cys - Arg - Arg - Ser - Thr - Thr - Asp - Cys-100 - Gly - Gly - Pro - Lys - Asp - His - Pro - Leu - Thr - Cys-110 - Asp - Asp - Pro - Arg - Phe - Gln - Asp - Ser - Ser - Ser - Ser - Lys - Ala - Pro - Pro - Pro - Ser - Leu - Pro - Ser-130 - Pro - Ser - Arg - Leu - Pro - Gly - Pro - Ser - Asp - Thr-140 - Pro - Ile - Leu - Pro - Gln. Oligosaccharide side chains are found at residues 13, 30, 121, 127, 132, and 138. The proteolytic enzyme, thrombin, which appears to cleave a limited number of arginyl bonds, proved helpful in the determination of the beta sequence.     

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.     

Characterization of the signal for rapid internalization of the bovine mannose 6-phosphate/insulin-like growth factor-II receptor.

The signal for rapid internalization of the mannose 6-phosphate/insulin-like growth factor II receptor has been localized to the sequence Tyr-Lys-Tyr-Ser-Lys-Val in positions 24-29 of its 163-residue cytoplasmic tail. Most of the activity of this signal is mediated by the carboxyl 4 amino acids, especially Tyr26 and Val29 (Canfield, W. M., Johnson, K. F., Ye, R. D., Gregory, W. and Kornfeld, S. (1991) J. Biol. Chem. 266, 5682-5688). In this study, we have tested the effect of a series of mutations on the internalization rate of a mutant receptor that contains a 29-amino acid cytoplasmic tail terminating with the 4-amino acid internalization sequence Tyr-Ser-Lys-Val. Replacement of Tyr26 with Phe or Trp gave rise to mutant receptors that were internalized at 10% the wild-type rate, while receptors with Ala, Leu, Ile, Val, or Asn at this position were totally inactive. Val29 could be replaced by other large hydrophobic residues (Phe, Leu, Ile, or Met) with no loss of activity, but the presence of Ala, Gly, Arg, Gln, or Tyr in this position inactivated the signal. Ser27 could be effectively replaced by many different amino acids, but not by Pro or Gly. However, Gly27 could be tolerated if the residues at positions 28 and 29 were also changed. A change in the 2-residue spacing between Tyr26 and Val29 destroyed the signal. These data show that the essential elements of this signal are an aromatic residue, especially a Tyr in the first position, separated from a large hydrophobic residue in the last position by 2 amino acids. The residues in positions 2 and 3 of the signal may have a modulating effect on its activity. The Tyr-Ser-Lys-Val signal could be moved to a more proximal region of the cytoplasmic tail with only a modest loss of activity. In addition, the signal could be effectively replaced by the putative 4-residue signals of seven other receptors and membrane proteins known to undergo rapid endocytosis, including the Tyr-Thr-Arg-Phe sequence of the transferrin receptor, a Type II membrane protein. These results are compatible with the 4-residue signals of this type being interchangeable, even among Type I and Type II membrane proteins.     

Distal Val346Ile mutation in inducible NO synthase promotes substrate-dependent NO confinement

The function of inducible NO synthase (WT iNOS) depends on the release of NO from the ferric heme before the enzyme is reduced. Key parameters controlling ligand dynamics include the distal and proximal heme pocket amino acids, as well as the inner solvent molecules. In this work, we tested how a point mutation in the distal heme side of WT iNOS affected the geminate rebinding of NO by ultrafast kinetics and molecular dynamics simulations. The mutation sequestered much of the photodissociated NO close to the heme compared to WT iNOS, with a main picosecond phase accounting for 78% of the rebinding to the arginine-bound Val346Ile protein. Consequently, the probability of NO release from Val346Ile decreased as compared to that from WT iNOS, provided the substrate binding site is filled. These data are rationalized by a steric effect of the Ile methyl group inducing events mediated by the substrate, transmitted via the propionates to the NO and the protein. This model is consistent with the role of the H-bonding network involving the heme, the substrate, and the BH4 cofactor in controlling NO release, with a key role of the heme propionates [Gautier et al. (2006) Nitric Oxide 15, 312]. These data support the effect of Val346Ile mutation in decreasing NO release and slowing down NO synthesis compared to WT iNOS determined by single turnover catalysis [Wang et al. (2004) J. Biol. Chem. 279, 19018].     

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.     

Ligand binding and protein relaxation in heme proteins: a room temperature analysis of NO geminate recombination

Ultrafast absorption spectroscopy is used to study heme-NO recombination at room temperature in aqueous buffer on time scales where the ligand cannot leave its cage environment. While a single barrier is observed for the cage recombination of NO with heme in the absence of globin, recombination in hemoglobin and myoglobin is nonexponential. Examination of hemoglobin with and without inositol hexaphosphate points to proximal constraints as important determinants of the geminate rebinding kinetics. Molecular dynamics simulations of myoglobin and heme-imidazole subsequent to ligand dissociation were used to investigate the transient behavior of the Fe-proximal histidine coordinate and its possible involvement in geminate recombination. The calculations, in the context of the absorption measurements, are used to formulate a distinction between nonexponential rebinding that results from multiple protein conformations (substates) present at equilibrium or from nonequilibrium relaxation of the protein triggered by a perturbation such as ligand dissociation. The importance of these two processes is expected to depend on the time scale of rebinding relative to equilibrium fluctuations and nonequilibrium relaxation. Since NO rebinding occurs on the picosecond time scale of the calculated myoglobin relaxation, a time-dependent barrier is likely to be an important factor in the observed nonexponential kinetics. The general implications of the present results for ligand binding in heme proteins and its time and temperature dependence are discussed. It appears likely that, at low temperatures, inhomogeneous protein populations play an important role and that as the temperature is raised, relaxation effects become significant as well.     

Picosecond geminate recombination of nitrosylmyoglobins

The kinetics of NO geminate recombination to sperm whale and elephant myoglobins has been studied on the picosecond time scale using an amplified colliding-pulse mode-locked ring dye laser. The dynamics of ligand rebinding are shown to be affected by the distal structure of the protein surrounding the heme pocket.     

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 amino acid sequence of hemoglobin II from the symbiont-harboring clamLucina pectinata

The cytoplasmic hemoglobin II from the gill of the clamLucina pectinata consists of 150 amino acid residues, has a calculatedM _m of 17,476, including heme and an acetylated N-terminal residue. It retains the invariant residues Phe 44 at position CD1 and His 65 at the proximal position F8, as well as the highly conserved Trp 15 at position A12 and Pro 38 at position C2. The most likely candidate for the distal residue at position E7, based on the alignment with other globins, is Gln 65. However, optical and EPR spectroscopic studies of the ferri Hb II (Kraus, D. W., Wittenberg, J. B., Lu, J. F., and Peisach, J.,J. Biol. Chem. 265, 16054–16059, 1990) have implicated a tyrosinate oxygen as the distal ligand. Modeling of theLucina Hb II sequence, using the crystal structure of sperm whale aquometmyoglobin, showed that Tyr 30 substituting for the Leu located at position B10 can place its oxygen within 2.8 Å of the water molecule occupying the distal ligand position. This structural alteration is facilitated by the coordinate mutation of the residue at position CD4, from Phe 46 in the sperm whale myoglobin sequence to Leu 47 inLucina Hb II.     

Structural dynamics of ligand diffusion in the protein matrix: A study on a new myoglobin mutant Y(B10) Q(E7) R(E10)

A triple mutant of sperm whale myoglobin (Mb) [Leu(B10) --> Tyr, His(E7) --> Gln, and Thr(E10) --> Arg, called Mb-YQR], investigated by stopped-flow, laser photolysis, crystallography, and molecular dynamics (MD) simulations, proved to be quite unusual. Rebinding of photodissociated NO, O2, and CO from within the protein (in a "geminate" mode) allows us to reach general conclusions about dynamics and cavities in proteins. The 3D structure of oxy Mb-YQR shows that bound O2 makes two H-bonds with Tyr(B10)29 and Gln(E7)64; on deoxygenation, these two residues move toward the space occupied by O2. The bimolecular rate constant for NO binding is the same as for wild-type, but those for CO and O2 binding are reduced 10-fold. While there is no geminate recombination with O2 and CO, geminate rebinding of NO displays an unusually large and very slow component, which is pretty much abolished in the presence of xenon. These results and MD simulations suggest that the ligand migrates in the protein matrix to a major "secondary site," located beneath Tyr(B10)29 and accessible via the motion of Ile(G8)107; this site is different from the "primary site" identified by others who investigated the photolyzed state of wild-type Mb by crystallography. Our hypothesis may rationalize the O2 binding properties of Mb-YQR, and more generally to propose a mechanism of control of ligand binding and dissociation in hemeproteins based on the dynamics of side chains that may (or may not) allow access to and direct temporary sequestration of the dissociated ligand in a docking site within the protein. This interpretation suggests that very fast (picosecond) fluctuations of amino acid side chains may play a crucial role in controlling O2 delivery to tissue at a rate compatible with physiology.     

Role of heme iron coordination and protein structure in the dynamics and geminate rebinding of nitric oxide to the H93G myoglobin mutant: implications for nitric oxide sensors

The influence of the heme iron coordination on nitric oxide binding dynamics was investigated for the myoglobin mutant H93G (H93G-Mb) by picosecond absorption and resonance Raman time-resolved spectroscopies. In the H93G-Mb, the glycine replacing the proximal histidine does not interact with the heme iron so that exogenous substituents like imidazole may coordinate to the iron at the proximal position. Nitrosylation of H93G-Mb leads to either 6- or 5-coordinate species depending on the imidazole concentration. At high concentrations, (imidazole)-(NO)-6-coordinate heme is formed, and the photoinduced rebinding kinetics reveal two exponential picosecond phases ( approximately 10 and approximately 100 ps) similar to those of wild type myoglobin. At low concentrations, imidazole is displaced by the trans effect leading to a (NO)-5-coordinate heme, becoming 4-coordinate immediately after photolysis as revealed from the transient Raman spectrum. In this case, NO rebinding kinetics remain bi-exponential with no change in time constant of the fast component whose amplitude increases with respect to the 6-coordinate species. Bi-exponential NO geminate rebinding in 5-coordinate H93G-Mb is in contrast with the single-exponential process reported for nitrosylated soluble guanylate cyclase (Negrerie, M., Bouzhir, L., Martin, J. L., and Liebl, U. (2001) J. Biol. Chem. 276, 46815-46821). Thus, our data show that the iron coordination state or the heme iron out-of-plane motion are not at the origin of the bi-exponential kinetics, which depends upon the protein structure, and that the 4-coordinate state favors the fast phase of NO geminate rebinding. Consequently, the heme coordination state together with the energy barriers provided by the protein structure control the dynamics and affinity for NO-binding enzymes.