Time-resolved resonance Raman study on the binding of carbon monoxide to recombinant human myoglobin and its distal histidine mutants

The Dynamics Behind the Affinity: Controlling Heme-Gas Affinity via Geminate Recombination and Heme Propionate Conformation in the NO Carrier Cytochrome c'

Nitric oxide (NO) sensors are heme proteins which may also bind CO and O₂. Control of heme-gas affinity and their discrimination are achieved by the structural properties and reactivity of the heme and its distal and proximal environments, leading to several energy barriers. In the bacterial NO sensor cytochrome c' from Alcaligenes xylosoxidans (AXCP), the single Leu16Ala distal mutation boosts the affinity for gas ligands by a remarkable 10⁶-10⁸-fold, transforming AXCP from one of the lowest affinity gas binding proteins to one of the highest. Here, we report the dynamics of diatomics after photodissociation from wild type and L16A-AXCP over 12 orders of magnitude in time. For the L16A variant, the picosecond geminate rebinding of both CO and NO appears with an unprecedented 100% yield, and no exit of these ligands from protein to solvent could be observed. Molecular dynamic simulations saliently demonstrate that dissociated CO stays within 4 Å from Fe²⁺, in contrast to wild-type AXCP. The L16A mutation confers a heme propionate conformation and docking site which traps the diatomics, maximizing the probability of recombination and directly explaining the ultrahigh affinities for CO, NO, and O₂. Overall, our results point to a novel mechanism for modulating heme-gas affinities in proteins.

The C-helix in CooA Rolls upon CO Binding to Ferrous Heme

CooA is a homodimeric transcriptional activator from Rhodospirillum rubrum containing one heme in each subunit. CO binding to the heme in its sensor domain activates CooA, facilitating the binding to DNA by its DNA-binding domain. The C-helix links the two domains and shapes an interface between the subunits. To probe the nature of CO activation, residues at positions 112-121 on the C-helix were replaced by Asn or Gln and their effects were evaluated by resonance Raman spectroscopy and by the measurements of CO binding affinity. The nu(Fe-CO) stretching Raman line in CO-bound wild-type CooA was up-shifted by 6 cm(-1) in the L116Q, G117N, and L120Q mutants, indicating unequivocally that these residues are close to the bound CO. Residues Leu116 and Leu120 from each subunit form contacts with the corresponding residues in the opposite subunit, enabling hydrophobic interactions in the inactive ferrous form. Thus, in the CO-bound activated form, both C-helices appear to roll to direct these residues toward the heme, forming a hydrophobic pocket for the bound CO. The CO affinity is approximately one order of magnitude higher in the L112Q, I115Q, L116Q, G117N, L120Q, and T121N mutants but reduced in A114N mutant. The variation indicates that these residues are close to the heme in the ferrous and/or CO-bound forms and are responsible for CooA activation. A roll-and-slide mechanism is proposed for CO activation of CooA.

The Escape Process of Carbon Monoxide from Myoglobin to Solution at Physiological Temperature

The carbon monoxide (CO) docking sites involved in the ligand escape process from the iron atom in hem of myoglobin (Mb) to solution at physiological temperature were studied on the basis of the effect of xenon (Xe) on the ligand escape rate by the transient grating (TG) technique. The TG method provides a direct measurement of the changes in molecular volume. The apparent CO escaping rate and the volume contraction increase with increasing Xe pressure. The pressure dependence of the rate is consistent with that of the Xe population at the Xe(1) site. This result clearly shows that CO is trapped at the Xe(1) site before escaping to solvent in a Xe-free solution at room temperature. It is shown that only CO but not the trapped Xe is released by the photoexcitation of the Xe-trapped MbCO. A dissociation scheme is proposed to explain the enhancement of the escaping rate by the presence of Xe(1). There are two branches for the CO escaping pathway. The dominant part of the dissociated CO escapes to the solvent through the Xe(1) trapping site under the Xe-free condition, and there are at least three intermediate states along this pathway. When a Xe atom blocks the Xe(1) site, the CO escapes through another route.

Carbon monoxide complex of cytochrome b(5) at acidic pH

CO complex of cyt b(5) generated at acidic pH is investigated by absorption, resonance Raman (RR), and far UV CD measurements. The Soret maximum wavelength blue-shifted to 420 nm with other absorption bands observed around 540 and 570 nm for reduced cyt b(5) upon interaction with CO at acidic pH (pH 3.1-3.5). Under this condition, the iron-carbon stretching RR band was observed at 529 cm(-1) (520 cm(-1) for C(18)O), which indicated formation of a heme&bond;CO adduct with a histidine as an axial ligand. Heme dissociated from the reduced cyt b(5) protein at pH approximately 3.5, whereas its rate decreased under CO atmosphere compared with N(2) atmosphere, due to formation of a heme&bond;CO adduct with a histidine as an axial ligand.

A chemometric analysis of the resonance Raman spectra of mutant carbonmonoxy-myoglobins reveals the effects of polarity

Resonance Raman spectra of 10 carbonmonoxy-myoglobins have been obtained, including sperm whale native, pig wild-type, and the mutants H64L, H64A, V68T, V68N, H64V/V68T, F43W, F46V, and L29F. This series was chosen in order to study the effect of ligand binding pocket polarity on the positions of the v(Fe-CO) and delta (Fe-C-O) bands. Spectra of both 12CO and 13CO isotopic forms have been obtained and a detailed analysis has facilitated the identification of both the ligand-specific bands and six underlying porphyrin bands which are insensitive to this isotopic substitution. Along with a band-fitting analysis of infrared spectra, these resonance Raman data provide a comprehensive evaluation of the vibrations of the FeCO unit. The band positions of the ligand-specific modes are found to depend on the structure of the ligand binding pocket, arising from the strength of back-bonding within the FeCO unit, and clear correlations exist between the v(Fe-CO), delta (Fe-C-O), and v(C-O) band positions which characterize this synergic bonding. The results are consistent with the proposal that the vibration band positions are determined primarily by the electrostatic potential at the ligand. Five discrete band sets are observed for this set of mutants, suggesting that 5 discrete conformations occur.

Theoretical study of the electrostatic and steric effects on the spectroscopic characteristics of the metal-ligand unit of heme proteins. 2. C-O vibrational frequencies, 17O isotropic chemical shifts, and nuclear quadrupole coupling constants

The quantum chemical calculations, vibronic theory of activation, and London-Pople approach are used to study the dependence of the C-O vibrational frequency, 17O isotropic chemical shift, and nuclear quadrupole coupling constant on the distortion of the porphyrin ring and geometry of the CO coordination, changes in the iron-carbon and iron-imidazole distances, magnitude of the iron displacement out of the porphyrin plane, and presence of the charged groups in the heme environment. It is shown that only the electrostatic interactions can cause the variation of all these parameters experimentally observed in different heme proteins, and the heme distortions could modulate this variation. The correlations between the theoretically calculated parameters are shown to be close to the experimentally observed ones. The study of the effect of the electric field of the distal histidine shows that the presence of the four C-O vibrational bands in the infrared absorption spectra of the carbon monoxide complexes of different myoglobins and hemoglobins can be caused by the different orientations of the different tautomeric forms of the distal histidine. The dependence of the 17O isotropic chemical shift and nuclear quadrupole coupling constant on pH and the distal histidine substitution can be also explained from the same point of view.

Recombinant [Phe(beta)63]hemoglobin shows rapid oxidation of the beta chains and low-affinity, non-cooperative oxygen binding to the alpha subunits

We have engineered alpha2beta2 [Phe63]hemoglobin by changing the highly conserved distal histidine of the beta chains to a phenylalanine. The mutant tetramer binds four high-affinity ligands, such as CO or NO, to the ferrous form, or CN to the oxidized iron; however, it binds only two low-affinity ligands, oxygen and azide. The absorption spectrum of the ferrous deoxy or ferric forms are not normal, displaying an enhanced absorption of the visible band near 560 nm. Half of the autooxidation process, attributed to the mutated beta subunits, is over 1000-fold faster than for Hb A. The mutant Hb exhibits non-cooperative binding of two oxygens with an affinity about fivefold lower than those of HbA valency hybrids (alpha met beta)2. Functional properties of this mutant Hb resemble those of Hb Saskatoon ([Tyr63]Hb) [Suzuki, T., Hayashi, A., Shimizu, A. & Yamamura, Y. (1966) Biochim. Biophys. Acta 127, 280-282]. Flash-photolysis experiments also indicate non-cooperative behaviour: the CO-recombination kinetics were independent of the fraction dissociated. Furthermore, the amplitude of the CO bimolecular phase was the same for the (alpha(CO)metbeta)2 valency hybrid or the (alphaCO betaCO)2 form, suggesting mainly geminate CO-recombination kinetics to the beta chains. EPR and Resonance Raman spectra did not show evidence for a hemichrome, normally considered as a six-coordinated iron with low-spin character. The EPR and resonance Raman spectra for the mutated beta subunits demonstrate the presence of a high-spin compound in the ferric and deoxy ferrous forms. In particular, the ferrous mutated beta subunits are penta-coordinated. The abnormal absorption spectra are possibly due to an interaction between the porphyrin and the phenyl ring in the distal position rather than to direct binding to the iron.

Nanosecond time-resolved spectroscopy of biomolecular processes

Over the past two decades, nanosecond absorption and vibrational spectroscopies have developed into powerful tools for monitoring the secondary, tertiary, and quaternary structural relaxations of biological macromolecules under near-physiological conditions of solvent and temperature. Observed through such methods, the dynamic response of a biomolecule to photoinitiated excursions from equilibrium can reveal valuable information about the structure-function relationship, information beyond that obtained from the static structures provided by X-ray crystallography, nuclear magnetic resonance spectroscopy, and other steady-state methods. Most recently, the development of ultra-sensitive polarization techniques for absorption spectroscopy has greatly enhanced the amount of time-resolved structural information that can be obtained from the broadened electronic spectra of biomolecules. This review examines nanosecond absorption, vibrational, and polarized absorption methods, and their applications to protein function and folding, emphasizing the complementary nature of information obtained from electronic and vibrational spectra measured on the nanosecond time scale.

Picosecond geminate recombination of CO to the complexes calmodulin*heme-CO and calmodulin*heme-CO*melittin

Picosecond CO recombination kinetics have been measured after photodissociation of the artificial complexes calmodulin*heme-CO and calmodulin*heme-CO*melittin. These systems show an enhancement of the geminate fraction of kinetics relative to unbound heme-CO, due in part to fast geminate kinetics (tau=50ps for the initial phase), as well as a decrease in the rate of migration of CO away from the binding site. This indicates that calmodulin provides a complete pocket around the heme group. Rather than competing with the hemes for binding to calmodulin, the melittin seems to act as a cap to further enclose the hemes; melittin increases the affinity of calmodulin for heme-CO, but only weakly affects the CO recombination kinetics.

Theoretical study of the distal-side steric and electrostatic effects on the vibrational characteristics of the FeCO unit of the carbonylheme proteins and their models

The vibronic theory of activation and quantum chemical intermediate neglect of differential overlap (INDO) calculations are used to study the activation of carbon monoxide (change of the C-O bond index and force field constant) by the imidazole complex with heme in dependence on the distortion of the porphyrin ring, geometry of the CO coordination, iron-carbon and iron-imidazole distances, iron displacement out of the porphyrin plane, and presence of the charged groups in the heme environment. It is shown that the main contribution to the CO activation stems from the change in the sigma donation from the 5 sigma CO orbital to iron, and back-bonding from the iron to the 2 pi orbital of CO. It follows from the results that none of the studied distortions can explain, by itself, the wide variation of the C-O vibrational frequency in the experimentally studied model compounds and heme proteins. To study the dependence of the properties of the FeCO unit on the presence of charged groups in the heme environment, the latter are simulated by the homogeneous electric field and point charges of different magnitude and location. The results show that charged groups can strongly affect the strength of the C-O bond and its vibrational frequency. It is found that the charges located on the distal side of the heme plane can affect the Fe-C and C-O bond indexes (and, consequently, the Fe-C and C-O vibrational frequencies), both in the same and in opposite directions, depending on their position. The theoretical results allow us to understand the peculiarities of the effect of charged groups on the properties of the FeCO unit both in heme proteins and in their model compounds.

Binding of nitric oxide and carbon monoxide to soluble guanylate cyclase as observed with Resonance raman spectroscopy

Resonance Raman spectra have been recorded for the ferrous heme of soluble guanylate cyclase (sGC), the only receptor known thus far for .NO. On the basis of the frequencies of porphyrin core sensitive vibrations in the high frequency region of the Raman spectrum, we conclude that the ferrous heme is five-coordinate, high spin, when no exogenous ligands are present. We assign a prominent vibration that occurs at 204 cm-1 in the reduced enzyme to the heme Fe(2+)-proximal histidine stretching vibration. In the .NO bound form of the enzyme, the heme Fe2+ retains a five-coordinate geometry. Assuming that .NO binds to the distal side of the heme, this observation indicates that the weak Fe-His bond breaks when .NO binds. Two isotope-sensitive vibrations are observed in the .NO bound enzyme, one at 1677 cm-1, attributed to the N-O stretching vibration, and one at 525 cm-1, attributed to the Fe-NO stretching vibration. When CO is bound to the ferrous heme, the heme ligation is six-coordinate. From this, we conclude that the Fe-His bond is still intact and that, if cleavage of the Fe-proximal ligand bond is necessary for complete activation of sGC, then CO should only weakly activate the enzyme, which has been shown to be the case. In the carbonmonoxy enzyme, the Fe-CO stretching vibration is observed at 472 cm-1 and the Fe-C-O bending vibration is detected at 562 cm-1. These frequencies are the lowest yet observed for the Fe-CO stretching and Fe-C-O bending modes in heme proteins or model systems with imidazole as the proximal ligand and suggest that there is significant negative polarity in the distal pocket. The negative polarity and the low frequency of the Fe-His stretching vibration may account for the very low O2 affinity of sGC.

Origin of the pH-dependent spectroscopic properties of pentacoordinate metmyoglobin variants

The pH dependence of the electronic and EPR spectra of two variants of horse heart myoglobin (Mb) in which the distal His64 ligand has been replaced by either Thr or Ile has been studied. Both of these variants exhibit spectroscopic changes with pH that are indicative of a transition between two ferric high-spin forms that occurs with a pKa of 9.49 for the His64Thr variant and 9.26 for the His64Ile variant and that is distinctly different from the pH-dependent spectroscopic changes related to titration of the distal aquo ligand of wild-type Mb. The electronic and EPR spectra of both variants at all values of pH studied are consistent with the presence of a pentacoordinate heme iron center. For the His64Thr variant, a high-resolution (1.9 A) structure determination establishes the lack of the distal aquo ligand and demonstrates an out-of-plane movement of the ferric iron toward the proximal histidine together with a decrease of the Fe-His bond length. Investigation of this pH-linked equilibrium by EPR spectroscopy reveals rhombically split high-spin signals at both pH 7 and 11 with a greater degree of rhombicity exhibited by the alkaline species. We propose that the pH-linked spectroscopic transition exhibited by these distal histidine variants results from the deprotonation of the proximal His93 residue to produce imidazolate ligation at alkaline pH.

Alteration of axial coordination by protein engineering in myoglobin. Bisimidazole ligation in the His64-->Val/Val68-->His double mutant

Pig and human myoglobin have been engineered to reverse the positions of the distal histidine and valine (i.e. His64(E7)-->Val and Val68(E11)-->His). Spectroscopic and ligand binding properties have been measured for human and pig H64V/V68H myoglobin, and the structure of the pig H64V/V68H double mutant has been determined to 2.07-A resolution by x-ray crystallography. The crystal structure shows that the N epsilon of His68 is located 2.3 A away from the heme iron, resulting in the formation of a hexacoordinate species. The imidazole plane of His68 is tilted relative to the heme normal; moreover it is not parallel to that of His93, in agreement with our previous proposal (Qin, J., La Mar, G. N., Dou, Y., Admiraal, S. J., and Ikeda-Saito, M. (1994) J. Biol. Chem. 269, 1083-1090). At cryogenic temperatures, the heme iron is in a low spin state, which exhibits a highly anisotropic EPR spectrum (g1 = 3.34, g2 = 2.0, and g3 < 1), quite different from that of the imidazole complex of metmyoglobin. The mean iron-nitrogen distance is 2.01 A for the low spin ferric state as determined by x-ray spectroscopy. The ferrous form of H64V/V68H myoglobin shows an optical spectrum that is similar to that of b-type cytochromes and consistent with the hexacoordinate bisimidazole hemin structure determined by the x-ray crystallography. The double mutation lowers the ferric/ferrous couple midpoint potential from +54 mV of the wild-type protein to -128 mV. Ferrous H64V/V68H myoglobin binds CO and NO to form stable complexes, but its reaction with O2 results in a rapid autooxidation to the ferric species. All of these results demonstrate that the three-dimensional positions of His64 and Val68 in the wild-type myoglobin are as important as the chemical nature of the side chains in facilitating reversible O2 binding and inhibiting autooxidation.

Distal residue-CO interaction in carbonmonoxy myoglobins: a molecular dynamics study of three distal mutants

Six 90-ps molecular dynamics trajectories, two for each of three distal mutants of sperm whale carbonmonoxy myoglobin, are reported; solvent waters within 16 A of the active site have been included. In both His64GIn trajectories, the distal side chain remains part of the heme pocket, forming a "closed" conformation similar to that of the wild type 64N delta H tautomer. Despite a connectivity more closely resembling the N epsilon H histidine tautomer, close interactions with the carbonyl ligand similar to those observed for the wild type 64N epsilon H tautomer are prevented in this mutant by repulsive interactions between the carbonyl O and the 64O epsilon. The aliphatic distal side chain of the His64Leu mutant shows little interaction with the carbonyl ligand in either His64Leu trajectory. Solvent water molecules move into and out of the active site in the His64Gly mutant trajectories; during all the other carbonmonoxy myoglobin trajectories, including the wild type distal tautomers considered in an earlier work, solvent molecules rarely encroach closer than 6 A of the active site. These results are consistent with a recent structural interpretation of the wild type infrared spectrum, and the current reinterpretation that the distal-ligand interaction in carbonmonoxy myoglobin is largely electrostatic, not steric, in nature.

Identification of the iron-carbonyl stretch in distal histidine mutants of carbonmonoxymyoglobin

Soret-excitation resonance Raman (RR) spectra are reported for six distal histidine mutants of carbonmonoxymyoglobin including H64A, H64V, H64L, H64I, H64W, and H64W/L29F. Based on 13CO isotope shifts, the iron-carbonyl stretching vibrations are unambiguously identified. The correct assignment of these modes eliminates the differences in the conformational substate occupations predicted by the RR versus IR data.

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.