Myoglobin structure and regulation of solvent accessibility of heme pocket

The effects of heme removal on the molecular structure of tuna and sperm whale myoglobin have been investigated by comparing the solvent accessibility to the heme pocket of the two proteins with that of the corresponding apoproteins. Although the heme microenvironment of tuna myoglobin is more polar than that of sperm whale myoglobin, the accessibility of solvent to heme is identical in the two proteins as revealed by thermal perturbation of Soret absorption. The removal of heme produces loss of helical folding and increase of solvent accessibility but the effects are rather different for the two proteins. More precisely, the loss of helical structure upon heme removal is 50% for tuna myoglobin and 15% for sperm whale myoglobin; moreover, the solvent accessibility of the heme pocket of tuna apomyoglobin is 2-3-fold greater than that of sperm whale apomyoglobin. These results have been explained in terms of the lack of helical folding in segment D, the structural organization of which may have a relevant effect in regulating the accessibility of ligands to the heme. The effects produced by charged quenchers reveal that the ligand path from the surface of the molecule to the ion atom of the heme involves a positively charged residue which may reasonably be identified as Arg-45 (sperm whale myoglobin) or Lys-41 (tuna myoglobin) on the basis of recent X-ray crystallographic information.     

Ligand binding to heme proteins: II. Transitions in the heme pocket of myoglobin.

Phenomena occurring in the heme pocket after photolysis of carbonmonoxymyoglobin (MbCO) below about 100 K are investigated using temperature-derivative spectroscopy of the infrared absorption bands of CO. MbCO exists in three conformations (A substrates) that are distinguished by the stretch bands of the bound CO. We establish connections among the A substates and the substates of the photoproduct (B substates) using Fourier-transform infrared spectroscopy together with kinetic experiments on MbCO solution samples at different pH and on orthorhombic crystals. There is no one-to-one mapping between the A and B substates; in some cases, more than one B substate corresponds to a particular A substate. Rebinding is not simply a reversal of dissociation; transitions between B substates occur before rebinding. We measure the nonequilibrium populations of the B substates after photolysis below 25 K and determine the kinetics of B substate transitions leading to equilibrium. Transitions between B substates occur even at 4 K, whereas those between A substates have only been observed above about 160 K. The transitions between the B substates are nonexponential in time, providing evidence for a distribution of substates. The temperature dependence of the B substate transitions implies that they occur mainly by quantum-mechanical tunneling below 10 K. Taken together, the observations suggest that the transitions between the B substates within the same A substate reflect motions of the CO in the heme pocket and not conformational changes. Geminate rebinding of CO to Mb, monitored in the Soret band, depends on pH. Observation of geminate rebinding to the A substates in the infrared indicates that the pH dependence results from a population shift among the substates and not from a change of the rebinding to an individual A substate.     

Structural and dynamic properties of the heme pocket in myoglobin probed by optical spectroscopy

We report the optical absorption spectra of sperm whale deoxy-, oxy-, and carbonmonoxymyoglobin in the temperature range 300–20 K and in 65% glycerol or ethylene glycol–water mixtures. By lowering the temperature, all bands exhibit half-width narrowing and peak frequency shift; moreover, the near-ir bands of deoxymyoglobin show a marked increase of the integrated intensities. Opposed to what has already been reported for human hemoglobin, the temperature dependence of the first moment of the investigated bands does not follow the behavior predicted by the harmonic Franck–Condon approximation and is sizably affected by the solvent composition; this solvent effect is larger in liganded than in nonliganded myoglobin. However, for all the observed bands the behavior of the second moment can be quite well rationalized in terms of the harmonic Franck–Condon approximation and is not dependent on solvent composition. On the basis of these data we put forward some suggestions concerning the structural and dynamic properties of the heme pocket in myoglobin and their dependence upon solvent composition. We also discuss the different behaviors of myoglobin and hemoglobin in terms of the different heme pocket structures and deformabilities of the two proteins.     

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.     

Pressure effects on the proximal heme pocket in myoglobin probed by Raman and near-infrared absorption spectroscopy.

The influence of high pressure on the heme protein conformation of myoglobin in different ligation states is studied using Raman spectroscopy over the temperature range from 30 to 295 K. Photostationary experiments monitoring the oxidation state marker bands demonstrate the change of rebinding rate with pressure. While frequency changes of vibrational modes associated with rigid bonds of the porphyrin ring are <1 cm(-1), we investigate a significant shift of the iron-histidine mode to higher frequency with increasing pressure (approximately 3 cm(-1) for deltaP = 190 MPa in Mb). The observed frequency shift is interpreted structurally as a conformational change affecting the tilt angle between the heme plane and the proximal histidine and the out-of-plane iron position. Independent evidence for iron motion comes from measurements of the redshift of band III in the near-infrared with pressure. This suggests that at high pressure the proximal heme pocket and the protein are altered toward the bound state conformation, which contributes to the rate increase for CO binding. Raman spectra of Mb and photodissociated MbCO measured at low temperature and variable pressure further support changes in protein conformation and are consistent with glasslike properties of myoglobin below 160 K.     

Rebinding and relaxation in the myoglobin pocket

The infrared stretching bands of carboxymyoglobin (MbCO) and the rebinding of CO to Mb after photodissociation have been studied in the temperature range 10-300 K in a variety of solvents. Four stretching bands imply that MbCO can exist in four substates, A0-A3. The temperature dependences of the intensities of the four bands yield the relative binding enthalpies and and entropies. The integrated absorbances and pH dependences of the bands permit identification of the substates with the conformations observed in the X-ray data (Kuriyan et al., J. Mol. Biol. 192 (1986) 133). At low pH, A0 is hydrogen-bonded to His E7. The substates A0-A3 interconvert above about 180 K in a 75% glycerol/water solvent and above 270 K in buffered water. No major interconversion is seen at any temperature if MbCO is embedded in a solid polyvinyl alcohol matrix. The dependence of the transition on solvent characteristics is explained as a slaved glass transition. After photodissociation at low temperature the CO is in the heme pocket B. The resulting CO stretching bands which are identified as B substates are blue-shifted from those of the A substates. At 40 K, rebinding after flash photolysis has been studied in the Soret, the near-infrared, and the integrated A and B substates. All data lie on the same rebinding curve and demonstrate that rebinding is nonexponential in time from at least 100 ns to 100 ks. No evidence for discrete exponentials is found. Flash photolysis with monitoring in the infrared region shows four different pathways within the pocket B to the bound substates Ai. Rebinding in each of the four pathways B----A is nonexponential in time to at least 10 ks and the four pathways have different kinetics below 180 K. From the time and temperature dependence of the rebinding, activation enthalpy distributions g(HBA) and preexponentials ABA are extracted. No pumping from one A substate to another, or one B substate to another, is observed below the transition temperature of about 180 K. If MbCO is exposed to intense white light for 10-10(3) s before being fully photolyzed by a laser flash, the amplitude of the long-lived states increases. The effect is explained in terms of a hierarchy of substates and substate symmetry breaking. The characteristics of the CO stretching bands and of the rebinding processes in the heme pocket depend strongly on the external parameters of solvent, pH and pressure. This sensitivity suggests possible control mechanisms for protein reactions.     

Waterproofing the heme pocket. Role of proximal amino acid side chains in preventing hemin loss from myoglobin

The ability of myoglobin to bind oxygen reversibly depends critically on retention of the heme prosthetic group. Globin side chains at the Leu(89)(F4), His(97)(FG3), Ile(99)(FG5), and Leu(104)(G5) positions on the proximal side of the heme pocket strongly influence heme affinity. The roles of these amino acids in preventing heme loss have been examined by determining high resolution structures of 14 different mutants at these positions using x-ray crystallography. Leu(89) and His(97) are important surface amino acids that interact either sterically or electrostatically with the edges of the porphyrin ring. Ile(99) and Leu(104) are located in the interior region of the proximal pocket beneath ring C of the heme prosthetic group. The apolar amino acids Leu(89), Ile(99), and Leu(104) "waterproof" the heme pocket by forming a barrier to solvent penetration, minimizing the size of the proximal cavity, and maintaining a hydrophobic environment. Substitutions with smaller or polar side chains at these positions result in exposure of the heme to solvent, the appearance of crystallographically defined water molecules in or near the proximal pocket, and large increases in the rate of hemin loss. Thus, the naturally occurring amino acid side chains at these positions serve to prevent hydration of the His(93)-Fe(III) bond and are highly conserved in all known myoglobins and hemoglobins.     

Proton magnetic resonance study of the influence of heme 2,4 substituents on the exchange rates of labile protons in the heme pocket of myoglobin.

Four exchangeable protons with large hyperfine shifts are assigned in the heme pocket of sperm whale met-cyano myoglobin reconstituted with heme possessing acetyl groups, ethyl groups, bromines, and hydrogens at the 2,4 position, using both relaxation and chemical-shift data. The four protons arise from the ring NH's of the proximal (F8), distal (E7), and FG2 histidines, and the peptide NH of His F8. The similarity of all chemical shifts to those of the native protein as well as the invariance of the relaxation rates of the distal histidyl ring NH dictate essentially the same structure for the heme cavity of both native and reconstituted proteins. The exchange rates with bulk water of the four labile proteins in each modified protein were determined by saturation-transfer and line width methods. All four labile protons were found to have the same exchange rate as in the native protein for acetyl and ethyl 2,4 substituents; the two resolved labile protons in the derivative with 2,4 bromine were also unchanged. The reconstituted protein with hydrogens at the 2,4 position exhibited slower exchange rates for three of the four protons, indicating an increased dynamic stability of the heme pocket in the absence of bulky 2,4 substituents.     

S-nitrosylation-induced conformational change in blackfin tuna myoglobin

S-nitrosylation is a post-translational protein modification that can alter the function of a variety of proteins. Despite the growing wealth of information that this modification may have important functional consequences, little is known about the structure of the moiety or its effect on protein tertiary structure. Here we report high-resolution x-ray crystal structures of S-nitrosylated and unmodified blackfin tuna myoglobin, which demonstrate that in vitro S-nitrosylation of this protein at the surface-exposed Cys-10 directly causes a reversible conformational change by "wedging" apart a helix and loop. Furthermore, we have demonstrated in solution and in a single crystal that reduction of the S-nitrosylated myoglobin with dithionite results in NO cleavage from the sulfur of Cys-10 and rebinding to the reduced heme iron, showing the reversibility of both the modification and the conformational changes. Finally, we report the 0.95-A structure of ferrous nitrosyl myoglobin, which provides an accurate structural view of the NO coordination geometry in the context of a globin heme pocket.     

Solvent accessibility of the thrombin-thrombomodulin interface

The kinetics of solvent accessibility at the protein-protein interface between thrombin and a fragment of thrombomodulin, TMEGF45, have been monitored by amide hydrogen/deuterium (H/2H) exchange detected by MALDI-TOF mass spectrometry. The interaction is rapid and reversible, requiring development of theory and experimental methods to distinguish H/2H exchange due to solvent accessibility at the interface from H/2H exchange due to complex dissociation. Association and dissociation rate constants were measured by surface plasmon resonance and amide H/2H exchange rates were measured at different pH values and concentrations of TMEGF45. When essentially 100% of the thrombin was bound to TMEGF45, two segments of thrombin became completely solvent-inaccessible, as evidenced by the pH insensitivity of the amide H/2H exchange rates. These segments form part of anion-binding exosite I and contain the residues for which alanine substitution abolishes TM binding. Several other regions of thrombin showed slowing of amide exchange upon TMEGF45 binding, but the exchange remained pH-dependent, suggesting that these regions of thrombin were rendered only partially solvent-inaccessible by TMEGF45 binding. These partially inaccessible regions of thrombin form both surface and buried contacts into the active site of thrombin and contain residues implicated in allosteric changes in thrombin upon TM binding.     

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

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

Solvent accessibility to flavin in oxynitrilase

Oxynitrilase containing 2-thioFAD [C(2) = S] in place of FAD exhibits catalytic activity similar to that of native enzyme. Reaction of methyl methanethiolsulfonate with 2-thioFAD bound to oxynitrilase results in the formation of the corresponding flavin disulfide [C(2)-SSCH3]. Normal flavin [C(2) = O] is formed by reacting 2-thioFAD oxynitrilase with m-chloroperoxybenzoate or H2O2. Both reactions proceed via a spectrally detectable flavin 2-S-oxide intermediate [C(2) = S+-O-], but sizable amounts of this intermediate accumulate only in the m-chloroperoxybenzoate reaction (about 40%). While similar reactions have been reported with free 2-thioflavin, kinetic and other data indicate that the oxynitrilase reactions occur with intact enzyme. This shows that the 2-position of the pyrimidine ring in the bound coenzyme is accessible to solvent. The data are consistent with previous studies on the reaction of peroxides with oxynitrilase-bound 5-deazaFAD which show that the pyrimidine ring is accessible at position 4. Analogous studies indicate that the pyrimidine ring is buried in the case of flavin bound to lactate oxidase, since the data indicate that both positions 2 and 4 are inaccessible to solvent.     

Solvent accessibility studies on polysaccharides

The solvent accessible surface area (ASA) of the polysaccharides, namely (i) carrageenan (1CAR); (ii) agarose (1AGA); (iii) guaran (GUR); (iv) capsular polysaccharide (1CAP); and (v) hyaluronan (1HUA), have been computed using the solvent accessibility technique of Lee and Richards. The results show that the average variation of ASA for the various atoms in the molecules lie in the range 1-30 A(2). Irrespective of position of sulfation, either at two or four in the sugar residues in 1CAR, the charged groups interact almost equally with the solvent. The ASA values for the chains A and B in 1AGA and 1CAR indicate that there are not much interchain interactions and the chains in both the molecules interact equally with the solvent. Residue-wise analysis indicates that the ASAs of residues vary alternately, high-low-high value pattern that is similar to that of the hydrophobic behaviour of beta-strands in proteins. The results also suggest that in these polysaccharides D-configuration residues have higher ASA than L-configuration residues.     

Solvent accessibility study on tRNAPhe

In order to assess the solvent–solute association in the tRNA^Phe molecule, solvent accessibility calculations were carried out for its crystalline and completely extended states following the method of Lee and Richards. To do this, results from the calculations on model trinucleotide systems pApXpA with different bases at position X were used. In the folded form of the molecule, it was found that the oxygen atoms O(I) and O(II) of almost all the phosphate groups and the O(2′) atoms of the sugar rings situated throughout the backbone were highly exposed to the solvent. The amount of reduction found in the solvent accessibilities of the various atoms in going from the extended state to the folded state of the molecule indicates the kind of compactness of the tertiary structure in tRNA^Phe. The results give quantitative support to many characteristics of the tRNA molecule, such as loop sections, buried/exposed residues, hydrophobic interactions, etc., which were thought to be due to other factors.     

Mutational study of the bacterial hemoglobin distal heme pocket

Ligand binding experiments on three mutants in the distal heme pocket of Vitreoscilla hemoglobin (GlnE7His, ProE8Ala, and GlnE7His,ProE8Ala) were used to probe the role of GlnE7 and ProE8 in the pocket's unusual structure. The oxygen dissociation constants for the wild type, E8Ala mutant, and E7His mutant proteins were 4.5, 4.7, and 1.7microM, respectively; the K(d) for the double mutant was not determinable by our technique. Visible-Soret spectra of the carbonyl and cyanyl forms and FT-IR of the carbonyl form of the E8 mutant were similar to those of the wild type; the opposite was true for the GlnE7His and GlnE7His,ProE8Ala mutants, which also differed from wild type in the visible-Soret spectra of their oxidized forms. Models of the effects of the mutations on distal pocket structure were consistent with the experimental findings, particularly the larger effects of the GlnE7His change.     

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.     

Spectroscopic characterization of heme A reconstituted myoglobin.

The focus of this study was to examine the functional role of the unusual peripheral substitution of heme A. The effects of heme A stereochemistry on the reconstitution of the porphyrin have been examined in the heme A-apo-myoglobin complex using optical absorption and resonance Raman and electron paramagnetic resonance spectroscopies. The addition of one equivalent of heme A to apo-Mb produces a complex which displays spectroscopic signals consistent with a distribution of high- and low-spin heme chromophores. These results indicate that the incorporation of heme A into apo-Mb significantly perturbs the protein refolding.     

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.     

1H NMR structure of the heme pocket of HNO-myoglobin

The unique (1)H NMR signal of nitrosyl hydride at 14.8 ppm is used to obtain a solution structure of the distal pocket of Mb-HNO, a rare nitroxyl adduct with a half-life of several months at room temperature. (1)H NMR, NOESY and TOCSY data were obtained under identical experimental conditions on solutions of the diamagnetic HNO and CO complexes of equine Mb, allowing direct comparison of NMR data to a crystallographically characterized structure. Twenty NOEs between the nitrosyl hydride and protein and heme-based signals were observed. The HNO orientation obtained by modeling the experimental (1)H NMR NOESY data yielded an orientation of ca. -104 degrees referenced to the N-Fe-N vector between alpha and beta mesoprotons. An essentially identical orientation was obtained by simple energy minimization of the HNO adduct using ESFF potentials, suggesting steric control of the orientation. Differences in chemical shifts are seen for protons on residues Phe43(CD1) and Val68(E11), but both exhibit virtually identical NOESY contacts to other residues, and thus are attributed to small movements of ca. 0.1 A within the strong ring current. The most significant differences are seen in the NOESY peak intensities and chemical shifts for the ring non-labile protons of the distal His64(E7). The orientation of the His64(E7) in Mb-HNO was analyzed on the basis of the NOESY cross-peak changes and chemical shift changes, predicting a ca. 20 degrees rotation about the beta-gamma bond. The deduced HNO and His64(E7) orientations result in geometry where the His64(E7) ring can serve as the donor for a significant H-bond to the oxygen atom of the bound HNO.     

Structural and functional aspects of the heart ventricle myoglobin of bluefin tuna

The heart ventricle myoglobin of bluefin tuna has been purified to an apparent homogeneity. The amino acid analysis has revealed only a limited number of substitutions between the myoglobins of yellowfin and bluefin tuna. The alpha-helix content of tuna myoglobin has been found considerably lower than that of mammalian myoglobin. No correlation has been discovered between the conformational stability and alpha-helix content. Denaturation experiments have shown that the whole structure of tuna myoglobin results from the interaction of two structural units which represent the product of independent folding processes. The structure of tuna myoglobin has been found more open and disorganized than that of sperm whale. This result has been related to the low content of electrostatic interactions and explained in terms of evolutive adaptations.