O2 and CO reactions with heme proteins: quantum yields and geminate recombination on picosecond time scales

Kinetic Study of Ligand Binding and Conformational Changes in Inducible Nitric Oxide Synthase

Nitric oxide synthases (NOSs) are heme enzymes that generate highly reactive nitric oxide from l-arginine (l-Arg) in a complex mechanism that is still only partially understood. We have studied carbon monoxide (CO) binding to the oxygenase domain of murine inducible NOS (iNOS) by using flash photolysis. The P420 and P450 forms of the enzyme, assigned to a protonated and unprotonated proximal cysteine, through which the heme is anchored to the protein, show markedly different CO rebinding properties. The data suggest that P420 has a widely open distal pocket that admits water. CO rebinding to the P450 form strongly depends on the presence of the substrate l-Arg, the intermediate N^ω-hydroxy-l-arginine, and the cofactor tetrahydrobiopterin. After adding these small molecules to the enzyme solution, the CO kinetics change slowly over the hours. This process can be described as a relaxation from a fast rebinding, metastable species to a slowly rebinding, thermodynamically stable species, which we associate with the enzymatically active form. Our results allow us to determine kinetic parameters of l-Arg binding to the ferrous deoxy iNOS protein for the first time and also provide clues regarding the nature of structural differences between the two interconverting species.

Role of Ionic Strength and pH in Modulating Thermodynamic Profiles Associated with CO Escape from Rice Nonsymbiotic Hemoglobin 1

Type 1 nonsymbiotic hemoglobins are found in a wide variety of land plants and exhibit very high affinities for exogenous gaseous ligands. These proteins are presumed to have a role in protecting plant cells from oxidative stress under etiolated/hypoxic conditions through NO dioxygenase activity. In this study we have employed photoacoustic calorimetry, time-resolved absorption spectroscopy, and classical molecular dynamics simulations in order to elucidate thermodynamics, kinetics, and ligand migration pathways upon CO photodissociation from WT and a H73L mutant of type 1 nonsymbiotic hemoglobin from Oryza sativa (rice). We observe a temperature dependence of the resolved thermodynamic parameters for CO photodissociation from CO-rHb1 which we attribute to temperature dependent formation of a network of electrostatic interactions in the vicinity of the heme propionate groups. We also observe slower ligand escape from the protein matrix under mildly acidic conditions in both the WT and H73L mutant (τ = 134 ± 19 and 90 ± 15 ns). Visualization of transient hydrophobic channels within our classical molecular dynamics trajectories allows us to attribute this phenomenon to a change in the ligand migration pathway which occurs upon protonation of the distal His73, His117, and His152. Protonation of these residues may be relevant to the functioning of the protein in vivo given that etiolation/hypoxia can cause a decrease in intracellular pH in plant cells.

The photochemistry of transition metal complexes using density functional theory

The use of density functional theory (DFT) and time-dependent DFT (TD-DFT) to study the photochemistry of metal complexes is becoming increasingly important among chemists. Computational methods provide unique information on the electronic nature of excited states and their atomic structure, integrating spectroscopy observations on transient species and excited-state dynamics. In this contribution, we present an overview on photochemically active transition metal complexes investigated by DFT. In particular, we discuss a representative range of systems studied up to now, which include CO- and NO-releasing inorganic and organometallic complexes, haem and haem-like complexes dissociating small diatomic molecules, photoactive anti-cancer Pt and Ru complexes, Ru polypyridyls and diphosphino Pt derivatives.

Differences between protein dynamics of hemoglobin upon dissociation of oxygen and carbon monoxide

Protein dynamics of human adult hemoglobin (HbA) upon ligand photolysis of oxygen (O(2)) and carbon monoxide (CO) was investigated using time-resolved resonance Raman (TR(3)) spectroscopy. The TR(3) spectra of the both photoproducts at 1-ns delay differed from that of the equilibrium deligated form (deoxy form) in the frequencies of the iron-histidine stretching [ν(Fe-His)] and methine wagging (γ(7)) modes, and the band intensity of pyrrole stretching and substituent bending (ν(8)) modes. Spectral changes of the O(2) photoproduct in the submicrosecond region were faster than those of the CO photoproduct, indicating that the structural dynamics following the photodissociation is ligand dependent for HbA. In contrast, no ligand dependence of the dynamics was observed for myoglobin, which has a structure similar to that of the subunit of HbA. The structural dynamics and relevance to the functionality of HbA also are discussed.

Ligand migration in human indoleamine-2,3 dioxygenase

Human indoleamine 2,3-dioxygenase (hIDO), a monomeric heme enzyme, catalyzes the oxidative degradation of L-tryptophan (L-Trp) and other indoleamine derivatives. Its activity follows typical Michaelis-Menten behavior only for L-Trp concentrations up to 50 μM; a further increase in the concentration of L-Trp causes a decrease in the activity. This substrate inhibition of hIDO is a result of the binding of a second L-Trp molecule in an inhibitory substrate binding site of the enzyme. The molecular details of the reaction and the inhibition are not yet known. In the following, we summarize the present knowledge about this heme enzyme.

Investigations of vibrational coherence in the low-frequency region of ferric heme proteins

Femtosecond coherence spectroscopy is applied to a series of ferric heme protein samples. The low-frequency vibrational spectra that are revealed show dominant oscillations near 40 cm(-1). MbCN is taken as a typical example of a histidine-ligated, six-coordinate, ferric heme and a comprehensive spectroscopic analysis is carried out. The results of this analysis reveal a new heme photoproduct species, absorbing near 418 nm, which is consistent with the photolysis of the His(93) axial ligand. The photoproduct undergoes subsequent rebinding/recovery with a time constant of approximately 4 ps. The photoproduct lineshapes are consistent with a photolysis quantum yield of 75-100%, although the observation of a relatively strong six-coordinate heme coherence near 252 cm(-1) (assigned to nu(9) in the MbCN Raman spectrum) suggests that the 75% lower limit is much more likely. The phase and amplitude excitation profiles of the low-frequency mode at 40 cm(-1) suggest that this mode is strongly coupled to the MbCN photoproduct species and it is assigned to the doming mode of the transient penta-coordinated material. The absolute phase of the 40 cm(-1) mode is found to be pi/2 on the red side of 418 nm and it jumps to 3pi/2 as excitation is tuned to the blue side of 418 nm. The absolute phase of the 40 cm(-1) signal is not explained by the standard theory for resonant impulsive stimulated Raman scattering. New mechanisms that give a dominant momentum impulse to the resonant wavepacket, rather than a coordinate displacement, are discussed. The possibilities of heme iron atom recoil after photolysis, as well as ultrafast nonradiative decay, are explored as potential ways to generate the strong momentum impulse needed to understand the phase properties of the 40 cm(-1) mode.

Quantum Chemical Evaluation of Protein Control over Heme Ligation: CO/O2 Discrimination in Myoglobin

Control of O2 versus CO binding in myoglobin (Mb) is tuned by a distal histidine residue through steric and H-bonding interactions. These interactions have been evaluated via Car-Parrinello DFT calculations, whose efficiency allows full quantum mechanical treatment of the 13 closest residues surrounding the heme. The small (8 degrees ) deviation of the Fe-C-O bond angle from linearity results from the steric influence of a distal valine residue and not the distal histidine. H-bond energies were evaluated by replacing the distal histidine with the non-H-bonding residue isoleucine. Binding energies for CO and O2 decreased by 0.8 and 4.1 kcal/mol for MbCO and MbO2, in good agreement with experimental H-bond estimates. Ligand discrimination is dominated by distal histidine H-bonding, which is also found to stabilize a metastable side-on isomer of MbO2 that may play a key role in MbO2 photodynamics.

Light-induced relaxation of photolyzed carbonmonoxy myoglobin: a temperature-dependent x-ray absorption near-edge structure (XANES) study

X-ray absorption near-edge structure (XANES) spectra at the Fe K-edge have been measured and compared on solution samples of horse carbonmonoxy-myoglobin and its photoproducts, prepared by two different photolysis protocols: 1), extended illumination at low temperature (15 K) by white light; and 2), slow-cool from 140 to 10 K at a rate of 0.5 K/min while illuminating the sample with a 532-nm continuous-wave laser source. CO recombination has been followed while increasing the temperature at a rate of 1.2 K/min. After extended illumination at 15 K, a single process is observed, corresponding to CO recombination from a completely photolyzed species with CO bound to the primary docking site (formally B-state, in agreement with previous x-ray diffraction studies). The temperature peak for this single process is approximately 50 K. Using slow-cool illumination, data show a two-state recombination curve, the two temperature peaks being roughly assigned to 50 K and 110 K. These results are in good agreement with previous FTIR studies using temperature-derivative spectroscopy. The XANES spectroscopic markers probe structural differences between the photoproduct induced by extended illumination at 15 K and the photoproduct induced by slow-cool illumination. These differences in the XANES data have been interpreted as due to light-induced Fe-heme relaxation that does not involve CO migration from the B-state. A quantitative description of the unrelaxed and relaxed B-states, including the measurements of the Fe-N(p), Fe-N(His), and Fe-CO distances, and the out-of-plane Fe displacement, has been obtained via a procedure (MXAN) recently developed by us. This work shows that XANES, being able to extract both kinetic and structural parameters in a single experiment, is a powerful tool for structural dynamic studies of proteins.

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

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

Role of Ligand Bending in the Photodissociation of O2 vs CO-heme: A Time-Dependent Density Functional Study

Time-dependent DFT calculations reveal a strong dependence of low-lying excited states on the <Fe-O-O angle of a model O2-heme adduct and ready access to a metastable side-on isomer. The predicted state crossings provide a pathway for internal conversion to the ground state and account for the low quantum yield for O2 dissociation, while the side-on isomer can account for the observation of an unphotolyzable fraction at low temperature. For CO-heme, the angle-sensitive states lie above the photoaccessed Q state and do not participate in photodissociation.

Measurements of the photodissociation quantum yields of MbNO and MbO(2) and the vibrational relaxation of the six-coordinate heme species

The (t approximately 0) photodissociation quantum yields (Y(0)) of MbNO and MbO(2) are measured to be 50 +/- 5 and 28 +/- 6%, respectively, using MbCO (Y(0) = 100%) as a reference. When photolysis does not take place, we find that a significant portion of the photon energy contributes to heating of the residual six-coordinate heme (MbNO and MbO(2)). The time constant for vibrational relaxation of the six-coordinate ligand-bound heme is found to be close to 1 ps for both samples. The MbO(2) sample also shows a approximately 4-ps optical response that is assigned to a rapid phase (25-30% amplitude) of O(2) geminate rebinding. We observe no additional geminate recombination in the MbO(2) sample out to 120 ps. In contrast, the MbNO sample displays significant geminate recombination over the first 120 ps, which can be adequately fit with two exponentials whose amplitudes and time constants appear to depend weakly on the pump wavelength. This more complex kinetic behavior conceivably arises due to heating of the photodissociated heme and its effect on the geminate recombination as the system cools. Overall, the data are consistent with a hypothesis that distortions along the iron-ligand bending coordinate play a key role in the photodissociation process. The transient formation of an unphotolyzable FeO(2) side-on binding geometry is suggested to be responsible for the lowered quantum yield of MbO(2) relative to MbNO.

Kinetic evidence for three photolyzable taxonomic conformational substates in oxymyoglobin

The kinetics of oxygen geminate binding with the taxonomic substates of MbO2 are reported. The maximum entropy method was used to analyze the rebinding kinetics of MbCO and MbO2 monitored in the Soret. The resulting rate distributions were found to consist of a small number of overlapping bands. A global parametric fit of a series of rate distributions recorded at several temperatures was performed using a Gaussian basis set to resolve the individual enthalpy distributions P(H). This approach was first validated by showing that the well-documented taxonomic substates of MbCO could be recovered. The method was then applied to MbO2. Three taxonomic substates were identified at pH 4.8, whereas only two of them contribute to oxygen geminate rebinding at pH 7.0. These findings show that, similarly to MbCO, MbO2 also exists as three photolyzable and kinetically different taxonomic substates and suggest reconsidering the issue of the photolysis quantum yield of MbO2.

Heme photolysis occurs by ultrafast excited state metal-to-ring charge transfer

Ultrafast time-resolved resonance Raman spectra of carbonmonoxy hemoglobin (Hb), nitroxy Hb, and deoxy Hb are compared to determine excited state decay mechanisms for both ligated and unligated hemes. Transient absorption and Raman data provide evidence for a sequential photophysical relaxation pathway common to both ligated and unligated forms of Hb* (photolyzed heme), in which the excited state 1Q decays sequentially: 1Q-->Hb*I-->Hb*II-->Hb ground state. Consistent with the observed kinetics, the lifetimes of these states are <50 fs, approximately 300 fs, and approximately 3 ps for 1Q, Hb*I, and Hb*II, respectively. The transient absorption data support the hypothesis that the Hb*I state results from an ultrafast iron-to-porphyrin ring charge transfer process. The Hb*II state arises from porphyrin ring-to-iron back charge transfer to produce a porphyrin ground state configuration a nonequilibrium iron d-orbital population. Equatorial d-pi* back-bonding of the heme iron to the porphyrin during the lifetime of the Hb*II state accounts for the time-resolved resonance Raman shifts on the approximately 3 ps time scale. The proposed photophysical pathway suggests that iron-to-ring charge transfer is the key event in the mechanism of photolysis of diatomic ligands following a porphyrin ring pi-pi* transition.

Crystal structures of myoglobin-ligand complexes at near-atomic resolution

We have used x-ray crystallography to determine the structures of sperm whale myoglobin (Mb) in four different ligation states (unligated, ferric aquomet, oxygenated, and carbonmonoxygenated) to a resolution of better than 1.2 A. Data collection and analysis were performed in as much the same way as possible to reduce model bias in differences between structures. The structural differences among the ligation states are much smaller than previously estimated, with differences of <0.25 A root-mean-square deviation among all atoms. One structural parameter previously thought to vary among the ligation states, the proximal histidine (His-93) azimuthal angle, is nearly identical in all the ferrous complexes, although the tilt of the proximal histidine is different in the unligated form. There are significant differences, however, in the heme geometry, in the position of the heme in the pocket, and in the distal histidine (His-64) conformations. In the CO complex the majority conformation of ligand is at an angle of 18 +/- 3 degrees with respect to the heme plane, with a geometry similar to that seen in encumbered model compounds; this angle is significantly smaller than reported previously by crystallographic studies on monoclinic Mb crystals, but still significantly larger than observed by photoselection. The distal histidine in unligated Mb and in the dioxygenated complex is best described as having two conformations. Two similar conformations are observed in MbCO, in addition to another conformation that has been seen previously in low-pH structures where His-64 is doubly protonated. We suggest that these conformations of the distal histidine correspond to the different conformational substates of MbCO and MbO(2) seen in vibrational spectra. Full-matrix refinement provides uncertainty estimates of important structural parameters. Anisotropic refinement yields information about correlated disorder of atoms; we find that the proximal (F) helix and heme move approximately as rigid bodies, but that the distal (E) helix does not.

Altered ligand rebinding kinetics due to distal-side effects in hemoglobin chico (Lysbeta66(E10) --> thr)

Hb Chico is an unusual human hemoglobin variant that has lowered oxygen affinity, but unaltered cooperativity and anion sensitivity. Previous studies showed these features to be associated with distal-side heme pocket alterations that confer increased structural rigidity on the molecule and that increase water content in the beta-chain heme pocket. We report here that the extent of nanosecond geminate rebinding of oxygen to the variant and its isolated beta-chains is appreciably decreased. Structural alterations in this variant decrease its oxygen recombination rates without significantly altering rates of migration out of the heme pocket. Data analysis indicates that one or more barriers that impede rebinding of oxygen from docking sites in the heme pocket are increased, with less consequence for CO rebinding. Resonance Raman spectra show no significant alterations in spectral regions sensitive to interactions between the heme iron and the proximal histidine residue, confirming that the functional differences in the variant are due to distal-side heme pocket alterations. These effects are discussed in the context of a schematic representation of heme pocket wells and barriers that could aid the design of novel hemoglobins with altered ligand affinity without loss of the normal allosteric responses that facilitate unloading of oxygen to respiring tissues.

On the origin of heme absorption band shifts and associated protein structural relaxation in myoglobin following flash photolysis

The role of the protein structural change monitored by absorption band shifts following flash photolysis of CO from myoglobin is discussed in terms of structure-function relationships. Evidence is presented that the Soret band shift does not depend primarily on the covalent linkage of the heme iron to the protein by using the mutation H93G(L) in which the proximal histidine 93 is replaced by glycine and an exogenous ligand L, which coordinates the heme iron but is not covalently bound to the globin. While CO rebinding kinetics depend strongly on the nature of the exogenous ligand L in H93G(L), the magnitude and time evolution of the Soret band shift in a viscous buffer on the nanosecond time scale are hardly perturbed in all cases studied. Comparison of the Soret band and charge transfer Band III shifts demonstrates that both have a similar time dependence on the nanosecond to microsecond time scale following flash photolysis in viscous solvents. We conclude that the nonexponential kinetics of protein relaxation probed by transient absorption band position shifts involves primarily distal coordinates prior to ligand escape. This result agrees with earlier measurements of Soret band shifts in distal pocket mutants of myoglobin (1). We suggest that the band shifts are primarily a response to changes in the electrostatic field around the heme (a transient Stark shift) associated with changes in protein structure that occur following ligand photodissociation.

Global mapping of structural solutions provided by the extended X-ray absorption fine structure ab initio code FEFF 6.01: structure of the cryogenic photoproduct of the myoglobin-carbon monoxide complex

X-ray methods based on synchrotron technology have the promise of providing time-resolved structural data based on the high flux and brightness of the X-ray beams. One of the most closely examined problems in this area of time-resolved structure determination has been the examination of intermediates in ligand binding to myoglobin. Recent crystallographic experiments using synchrotron radiation have identified the protein tertiary and heme structural changes that occur upon photolysis of the myoglobin--carbon monoxide complex at cryogenic temperatures [Schlichting, I., Berendzen, J., Phillips, G., & Sweet, R. (1994) Nature 371, 808--812]. However, the precision of protein crystallographic data (approximately 0.2 A) is insufficient to provide precise metrical details of the iron--ligand bond lengths. Since bond length changes on this scale can trigger reactivity changes of several orders of magnitude, such detail is critical to a full understanding of metalloprotein structure--function relationships. Extended X-ray absorption fine structure (EXAFS) spectroscopy has the potential for analyzing bond distances to a precision of 0.02 A but is hampered by its relative insensitivity to the geometry of the backscattering atoms. Thus, it is often unable to provide a unique solution to the structure without ancillary structural information. We have developed a suite of computer programs that incorporate this ancillary structural information and compute the expected experimental spectra for a wide ranging series of Cartesian coordinate sets (global mapping). The programs systematically increment the distance of the metal to various coordinating ligands (along with their associated higher shells). Then, utilizing the ab initio EXAFS code FEFF 6.01, simulated spectra are generated and compared to the actual experimental spectra, and the differences are computed. Finally, the results for hundreds of simulations can be displayed (and compared) in a single plot. The power of this approach is demonstrated in the examination of high signal to noise EXAFS data from a photolyzed solution sample of the myoglobin--carbon monoxide complex at 10 K. Evaluation of these data using our global mapping procedures placed the iron to pyrrole nitrogen average distances close to the value for deoxymyoglobin (2.05 +/- 0.01 A), while the distance from iron to the proximal histidine nitrogen is seen to be 2.20 +/- 0.04 A. It is also shown that one cannot uniquely position the CO ligand on the basis of the EXAFS data alone, as a number of reasonable minima (from the perspective of the EXAFS) are observed. This provides a reasonable explanation for the multiplicity of solutions that have been previously reported. The results presented here are seen to be in complete agreement with the crystallographic results of Schlichting et al. (1994) within the respective errors of the two techniques; however, the extended X-ray absorption fine structure data allow the iron--ligand bond lengths to be precisely defined. An examination of the available spectroscopic data, including EXAFS, shows that the crystallographic results of Schlichting et al. (1994) are highly relevant to the physiological solution state and must be taken into account in any attempt to understand the incomplete relaxation process of the heme iron for the Mb*CO photoproduct at low temperature.

Picosecond infrared study of the photodynamics of carbonmonoxy-cytochrome c oxidase

Time-resolved infrared (TRIR) techniques have been employed to study the reactions of carbon monoxide with the cytochrome alpha 3-Cu(B) site of cytochrome c oxidase (CcO). The ligation dynamics immediately following photodissociation have been investigated using picosecond TRIR spectroscopy and linear dichroism. The rate of photoinitiated transfer of CO from cytochrome alpha 3 to CuB was measured directly by monitoring the development of the transient CuBCO absorption. In less than 1 ps, a stationary CuBCO spectrum develops, which together with the CO infrared linear dichroism is constant until the CO dissociates from CuB on a microsecond time scale. These observations indicate that the CO is transferred between metals and reaches its equilibrium conformation in less than 1 ps. This unprecedented ligand transfer rate has profound implications with regard to the structure and dynamics of the cytochrome alpha 3-CuB site, the functional architecture of the protein, and coordination dynamics in general.

Time-resolved absorption and magnetic circular dichroism spectroscopy of cytochrome c3 from Desulfovibrio

The UV-visible absorption and magnetic circular dichroism (MCD) spectra of the ferric, ferrous, CO-ligated forms and kinetic photolysis intermediates of the tetraheme electron-transfer protein cytochrome c3 (Cc3) are reported. Consistent with bis-histidinyl axial coordination of the hemes in this Class III c-type cytochrome, the Soret and visible region MCD spectra of ferric and ferrous Cc3 are very similar to those of other bis-histidine axially coordinated hemeproteins such as cytochrome b5. The MCD spectra indicate low spin state for both the ferric (S = 1/2) and ferrous (S = 0) oxidation states. CO replaces histidine as the axial sixth ligand at each heme site, forming a low-spin complex with an MCD spectrum similar to that of myoglobin-CO. Photodissociation of Cc3-CO (observed photolysis yield = 30%) produces a transient five-coordinate, high-spin (S = 2) species with an MCD spectrum similar to deoxymyoglobin. The recombination kinetics of CO with heme Fe are complex and appear to involve at least five first-order or pseudo first-order rate processes, corresponding to time constants of 5.7 microseconds, 62 microseconds, 425 microseconds, 2.9 ms, and a time constant greater than 1 s. The observed rate constants were insensitive to variation of the actinic photon flux, suggesting noncooperative heme-CO rebinding. The growing in of an MCD signal characteristic of bis-histidine axial ligation within tens of microseconds after photodissociation shows that, although heme-CO binding is thermodynamically favored at 1 atm CO, binding of histidine to the sixth axial site competes kinetically with CO rebinding.