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.     

Hydrophobic distal pocket affects NO-heme geminate recombination dynamics in dehaloperoxidase and H64V myoglobin

The recombination dynamics of NO with dehaloperoxidase (DHP) from Amphitrite ornata following photolysis were measured by femtosecond time-resolved absorption spectroscopy. Singular value decomposition (SVD) analysis reveals two important basis spectra. The first SVD basis spectrum reports on the population of photolyzed NO molecules and has the appearance of the equilibrium difference spectrum between the deoxy and NO forms of DHP. The first basis time course has two kinetic components with time constants of tau(11) approximately 9 ps and tau(12) approximately 50 ps that correspond to geminate recombination. The fast geminate process tau(11) arises from a contact pair with the heme iron in a bound state with S = 3/2 spin. The slow geminate process tau(12) corresponds to the recombination from a more remote docking site >3 A from the heme iron with the greater barrier corresponding to a S = 5/2 spin state. The second SVD basis spectrum represents a time-dependent Soret band shift indicative of heme photophysical processes and protein relaxation with time constants of tau(21) approximately 3 ps and tau(22) approximately 17 ps, respectively. A comparison between the more rapid rate constant of the slow geminate phase in DHP-NO and horse heart myoglobin (HHMbNO) or sperm whale myoglobin (SWMbNO) suggests that protein interactions with photolyzed NO are weaker in DHP than in the wild-type MbNOs, consistent with the hydrophobic distal pocket of DHP. The slower protein relaxation rate tau(22) in DHP-NO relative to HHMbNO implies less effective trapping in the docking site of the distal pocket and is consistent with a greater yield for the fast geminate process. The trends observed for DHP-NO also hold for the H64V mutant of SWMb (H64V MbNO), consistent with a more hydrophobic distal pocket for that protein as well. We examine the influence of solution viscosity on NO recombination by varying the glycerol content in the range from 0% to 90% (v/v). The dominant effect of increasing viscosity is the increase of the rate of the slow geminate process, tau(12), coupled with a population decrease of the slow geminate component. Both phenomena are similar to the effect of viscosity on wild-type Mb due to slowing of protein relaxation resulting from an increased solution viscosity and protein surface dehydration.     

Dynamics of nitric oxide rebinding and escape in horseradish peroxidase

Ultrafast kinetic measurements of NO rebinding to horseradish peroxidase (HRP) are reported for the first time. The geminate kinetics are found to be exponential for all HRP samples studied. The ferric forms of HRP have NO geminate recombination time constants in the range of 15-30 ps, while the ferrous form has a time constant of approximately 7 ps. The simple exponential NO geminate kinetics found for HRP demonstrate that heme relaxation is not the underlying source of the nonexponential NO rebinding in myoglobin (Mb). The NO ligand escape rates from HRP are also determined, and they are found to depend dramatically on the presence or absence of the competitive inhibitor benzohydroxamic acid (BHA). The kinetic results indicate that, in contrast to Mb, there is direct solvent access to the distal heme pocket of HRP.     

Time resolved infrared absorption studies of geminate recombination and vibrational relaxation in OClO photochemistry

Ultrafast time-resolved infrared absorption studies of aqueous chlorine dioxide (OClO) photochemistry are reported. Following photoexcitation at 401 nm, the evolution in optical density at frequencies between 1000 to 1100 cm(-1) is monitored to investigate vibrational energy deposition and relaxation along the asymmetric-stretch coordinate following the reformation of ground-state OClO via geminate recombination of the primary photofragments. The measured kinetics are compared to two proposed models for the vibrational-relaxation dynamics along the asymmetric-stretch coordinate. This comparison demonstrates that the perturbation model derived from molecular dynamics studies is capable of qualitatively reproducing the observed kinetics, where the collisional model employed in previous UV-pump, visible probe experiments demonstrates poor agreement with experiment. The ability of the perturbation model to reproduce the optical-density evolution observed in these studies demonstrates that for aqueous OClO, frequency dependence of the solvent-solute coupling is important in defining the level-dependent vibrational relaxation rates along the asymmetric-stretch coordinate. The absence of optical-density evolution corresponding to the population of higher vibrational levels (n>8) along the asymmetric-stretch coordinate suggests that following geminate recombination, energy is initially deposited into a local Cl-O stretch, with the relaxation of vibrational energy from this coordinate providing for delayed vibrational excitation of the asymmetric- and symmetric-stretch coordinates relative to geminate recombination, as previously observed.     

Nonequilibrium dynamics simulations of nitric oxide release: comparative study of nitrophorin and myoglobin

Nitrophorin 4 (NP4) is a heme protein that reversibly binds nitric oxide (NO), with release rates modulated by pH change. High-resolution structures of NP4 revealed that pH changes and NO binding induce a large conformational rearrangement in two loops that serve to protect the heme-bound NO molecule from solvent. We used extended (110 ns) molecular dynamics simulations of NP4 at pH 5 and pH 7, modeled by selective deprotonation of acidic groups. Conformational and dynamic changes were observed, consistent with those found in the crystal. Further, major solvent movement and NO escape were observed at pH 7, while the ligand remained in the heme binding pocket at pH 5. As a control, we also performed molecular dynamics (MD) simulations of sperm whale myoglobin, where NO migration into the interior cavities of the protein was observed, consistent with previous reports. We constructed a kinetic model of ligand escape to quantitatively relate the microscopic rate constants to the observed rates, and tested the predictions against the experimental data. The results suggest that release rates of diatomic molecules from heme proteins can be varied by several orders of magnitude through modest adjustments in geminate rebinding and gating behavior.     

Photodissociation and geminate dynamics in solution phase: Picosecond transient absorption studies of tetraphenylhydrazines and diphenyl disulfides

Photodissociation of aromatic molecules and the following geminate dynamics of the photochemically produced fragment radical pair in liquid phase are discussed in relation to the results of ultrafast transient absorption measurements. Photodissociation of tetraphenylhydrazine occurs from the unequilibrated excited singlet state competing with the relaxation to the fluorescence state. In addition to the rapid fragmentation, the N?N bond rupture is clearly demonstrated to take place in the fluorescence state of tetraphenyl-and tetratolyl-hydrazine in tens of picoseconds time scale. In the case of di-p-aminophenyl disulfide in nonpolar solvents, geminate recombination of the p-aminophenylthiyl radical and the formation of the radical dimer were observed in tens of picoseconds. The observed slow geminate recombination can be explained by the repulsion between the large dipole of the fragment radical preventing the reformation of the S?S bond and by the restricted conformation required for the dimer formation at the encounter of radicals.     

CHEMILUMINESCENCE and EPR STUDIES ON THE EXCITATION SITE OF FERRIC-HEME-OXO COMPLEXES OF NATURAL and MODEL HEME SYSTEMS

Chemiluminescence was detected both in the reaction system of H₂O₂ plus heme proteins such as methemo- and metmyoglobin and ferric-protoheme complexes used as a model system. The intensity of chemiluminescence was found to be mediated by ligand binding to the sixth coordination site of the ferric-protoheme compounds, e.g. chemiluminescence was not observed with the bisimidazole ferric-protoheme complex. On the other hand the pentacoordinated histidine ferric-protoheme complex exhibited strong light emission. Comparative studies with various ligand-heme compounds elucidated that light emission was inversely correlated with the binding strength of the respective ligand at the sixth coordination site. The basic reaction mechanism causing the establishment of an excited state was studied by monitoring chemiluminescence and EPR signal formation of ligand-modified heme proteins in the presence of different electron donors. External electron donors such as Trolox C, TMPD and ascorbic acid affected a strong reduction in the development of chemiluminescence suggesting the essential involvement of an inner-molecular electron transfer process. Our results allow the conclusion that chemiluminescence is generated from the decay of an excited state of oxo-heme compounds established as a result of a one electron transfer step from a ligand group to heme iron.     

Towards understanding non-equivalence of α and β subunits within human hemoglobin in conformational relaxation and molecular oxygen rebinding

Picosecond to millisecond laser time-resolved transient absorption spectroscopy was used to study molecular oxygen (O₂) rebinding and conformational relaxation following O₂ photodissociation in the α and β subunits within human hemoglobin in the quaternary R-like structure. Oxy-cyanomet valency hybrids, α₂(Fe²⁺–O₂)β₂(Fe³⁺–CN) and α₂(Fe³⁺–CN)β₂(Fe²⁺–O₂), were used as models for oxygenated R-state hemoglobin. An extended kinetic model for geminate O₂ rebinding in the ferrous hemoglobin subunits, ligand migration between the primary and secondary docking site(s), and nonexponential tertiary relaxation within the R quaternary structure, was introduced and discussed. Significant functional non-equivalence of the α and β subunits in both the geminate O₂ rebinding and concomitant structural relaxation was revealed. For the β subunits, the rate constant for the geminate O₂ rebinding to the unrelaxed tertiary structure and the tertiary transition rate were found to be greater than the corresponding values for the α subunits. The conformational relaxation following the O₂ photodissociation in the α and β subunits was found to decrease the rate constant for the geminate O₂ rebinding, this effect being more than one order of magnitude greater for the β subunits than for the α subunits. Evidence was provided for the modulation of the O₂ rebinding to the individual α and β subunits within human hemoglobin in the R-state structure by the intrinsic heme reactivity through a change in proximal constraints upon the relaxation of the tertiary structure on a picosecond to microsecond time scale. Our results demonstrate that, for native R-state oxyhemoglobin, O₂ rebinding properties and spectral changes following the O₂ photodissociation can be adequately described as the sum of those for the α and β subunits within the valency hybrids. The isolated β chains (hemoglobin H) show similar behavior to the β subunits within the valency hybrids and can be used as a model for the β subunits within the R-state oxyhemoglobin. At the same time, the isolated α chains behave differently to the α subunits within the valency hybrids.

O₂ rebinding and conformational relaxation following O₂ photodissociation were studied on picosecond to millisecond time scale in the α and β subunits within human hemoglobin in the quaternary R-like structure.     

Slow geminate‐charge‐pair recombination dynamics at polymer: Fullerene heterojunctions in efficient organic solar cells

We explore charge recombination dynamics at electron donor‐acceptor heterojunctions, formed between a semiconductor polymer (PCDTBT) and a fullerene derivative (PC₇₀BM), by means of combined time‐resolved photoluminescence and transient absorption spectroscopies. Following prompt exciton dissociation across the heterojunction, a subset of bound electron‐hole pairs recombines with a temperature‐independent rate distribution spanning submicrosecond timescales to produce luminescent charge‐transfer excitons (CTX). At 14 K, this slow mechanism is the dominant geminate charge recombination pathway, whereas we also observe CTX emission on subnanosecond timescales at 293 K. We thus find that at these temperatures, a fraction of the initial charge‐pair population is trapped deeply such that they only recombine slowly over a broad distribution of timescales by quantum tunneling. We identify geminate polaron pairs (GPP) as a reservoir of long‐lived localized states that repopulate the CTX up to microsecond timescales. The observation of such distributed geminate‐charge recombination highlights the importance of the molecular nature of specific donor–acceptor electronic interactions in defining the relaxation pathways of trapped GPP. © 2012 Wiley Periodicals, Inc. J Polym Sci Part B: Polym Phys, 2012     

Ultrafast heme dynamics in ferrous versus ferric cytochrome c studied by time-resolved resonance Raman and transient absorption spectroscopy

Cytochrome c (Cyt c) is a heme protein involved in electron transfer and also in apoptosis. Its heme iron is bisaxially ligated to histidine and methionine side chains and both ferric and ferrous redox states are physiologically relevant, as well as a ligand exchange between internal residue and external diatomic molecule. The photodissociation of internal axial ligand was observed for several ferrous heme proteins including Cyt c, but no time-resolved studies have been reported on ferric Cyt c. To investigate how the oxidation state of the heme influences the primary photoprocesses, we performed a comprehensive comparative study on horse heart Cyt c by subpicosecond time-resolved resonance Raman and femtosecond transient absorption spectroscopy. We found that in ferric Cyt c, in contrast to ferrous Cyt c, the photodissociation of an internal ligand does not take place, and relaxation dynamics is dominated by vibrational cooling in the ground electronic state of the heme. The intermolecular vibrational energy transfer was found to proceed in a single phase with a temperature decay of approximately 7 ps in both ferric and ferrous Cyt c. For ferrous Cyt c, the instantaneous photodissociation of the methionine side chain from the heme iron is the dominant event, and its rebinding proceeds in two phases, with time constants of approximately 5 and approximately 16 ps. A mechanism of this process is discussed, and the difference in photoinduced coordination behavior between ferric and ferrous Cyt c is explained by an involvement of the excited electronic state coupled with conformational relaxation of the heme.     

The Factors Influencing Nonlinear Characteristics of the Short‐Circuit Current in Dye‐Sensitized Solar Cells Investigated by a Numerical Model

A numerical model for interpretation of the light‐intensity‐dependent nonlinear characteristics of the short‐circuit current in dye‐sensitized solar cells is suggested. The model is based on the continuity equation and includes the influences of the nongeminate recombination between electrons and electron acceptors in the electrolyte and the geminate recombination between electrons and oxidized dye molecules. The influences of the order and rate constant of the nongeminate recombination reaction, the light‐absorption coefficient of the dye, the film thickness, the rate constant of geminate recombination, and the regeneration rate constant on the nonlinear characteristics of the short‐circuit current are simulated and analyzed. It is proposed that superlinear and sublinear characteristics of the short‐circuit current should be attributed to low electron‐collection efficiency and low dye‐regeneration efficiency, respectively. These results allow a deep understanding of the origin of the nonlinear characteristics of the short‐circuit current in solar cells.     

Ultrafast Heme Dynamics of Ferric Cytochrome c in Different Environments: Electronic,Vibrational, and Conformational Relaxation

The excited‐state dynamics of ferric cytochrome c (Cyt c), an important electron‐transfer heme protein, in acidic to alkaline medium and in its unfolded form are investigated by using femtosecond pump–probe spectroscopy, exciting the heme and Tryptophan (Trp) to understand the electronic, vibrational, and conformational relaxation of the heme. At 390 nm excitation, the electronic relaxation of heme is found to be ≈150 fs at different pH values, increasing to 480 fs in the unfolded form. Multistep vibrational relaxation dynamics of the heme, including fast and slow processes, are observed at pH 7. However, in the unfolded form and at pH 2 and 11, fast phases of vibrational relaxation dominate, revealing the energy dissipation occurring through the covalent bond interaction between the heme and the nearest amino acids. A significant shortening of the excited‐state lifetime of Trp is observed at various pH values at 280 nm excitation due to resonance energy transfer to the heme. The longer time constant (25 ps) observed in the unfolded form is attributed to a complete global conformational relaxation of Cyt c.     

Intermolecular hydrogen bonding in chlorine dioxide photochemistry: A time-resolved resonance Raman study

The geminate-recombination and vibrational-relaxation dynamics of chlorine dioxide (OClO) dissolved in ethanol and 2,2,2-trifluoroethanol (TFE) are investigated using time-resolved resonance Raman spectroscopy. Stokes spectra are measured as a function of time following photoexcitation using degenerate pump and probe wavelengths of 398 nm. For OClO dissolved in ethanol, subpicosecond geminate recombination occurs resulting in the reformation of ground-state OClO with a quantum yield of 0.5±0.1. Following recombination, intermolecular-vibrational relaxation of OClO occurs with a time constant of 31±10 ps. For OClO dissolved in TFE, recombination occurs with a time constant of 1.8±0.8 ps and a quantum yield of only 0.3±0.1. The intermolecular-vibrational-relaxation time constant of OClO in TFE is 79±27 ps. The reduced geminate-recombination quantum yield, delayed recombination, and slower vibrational relaxation for OClO in TFE is interpreted in terms of greater self-association of the solvent. Degenerate pump–probe experiments are also presented that demonstrate decay of the Cl-solvent charge-transfer complex on the ∼1-ns time scale in ethanol and TFE. This time is significantly longer than the abstraction times observed for other systems demonstrating that Cl hydrogen abstraction from alcohols occurs in the presence of a significant energy barrier.     

Spin relaxation effects in photochemically induced dynamic nuclear polarization spectroscopy of nuclei with strongly anisotropic hyperfine couplings

We describe experimental results and theoretical models for nuclear and electron spin relaxation processes occurring during the evolution of 19F-labeled geminate radical pairs on a nanosecond time scale. In magnetic fields of over 10 T, electron-nucleus dipolar cross-relaxation and longitudinal DeltaHFC-Deltag (hyperfine coupling anisotropy--g-tensor anisotropy) cross-correlation are shown to be negligibly slow. The dominant relaxation process is transverse DeltaHFC-Deltag cross-correlation, which is shown to lead to an inversion in the geminate 19F chemically induced dynamic nuclear polarization (CIDNP) phase for sufficiently large rotational correlation times. This inversion has recently been observed experimentally and used as a probe of local mobility in partially denatured proteins (Khan, F.; et al. J. Am. Chem. Soc. 2006, 128, 10729-10737). The essential feature of the spin dynamics model employed here is the use of the complete spin state space and the complete relaxation superoperator. On the basis of the results reported, we recommend this approach for reliable treatment of magnetokinetic systems in which relaxation effects are important.     

Structural relaxation in complex liquids: non-Markovian dynamics in a bistable potential

The time correlation function C(t) identical with of the distance fluctuations of a particle moving in a bistable potential under the action of fractional Gaussian noise (fGn) is calculated from a Smoluchowski-type equation derived from a generalized Langevin equation (GLE). The time derivative of this function, dC(t)dt, is compared with data from optical Kerr effect measurements of liquid crystal dynamics in the vicinity of the isotropic-to-nematic transition, which are related to the time derivative of an orientational correlation function. A number of characteristic features of the experimental decay curves, including short and intermediate time power law behavior and long time exponential relaxation, are qualitatively reproduced by the analytical calculations, even though the latter do not explicitly treat orientational degrees of freedom. The GLE formalism with fGn was, in fact, originally proposed as a model of protein conformational fluctuations, so the present results suggest that it may also serve more generally as a model of structural relaxation in complex condensed phase media.     

Strong ligand-protein interactions revealed by ultrafast infrared spectroscopy of CO in the heme pocket of the oxygen sensor FixL

In heme-based sensor proteins, ligand binding to heme in a sensor domain induces conformational changes that eventually lead to changes in enzymatic activity of an associated catalytic domain. The bacterial oxygen sensor FixL is the best-studied example of these proteins and displays marked differences in dynamic behavior with respect to model globin proteins. We report a mid-IR study of the configuration and ultrafast dynamics of CO in the distal heme pocket site of the sensor PAS domain FixLH, employing a recently developed method that provides a unique combination of high spectral resolution and range and high sensitivity. Anisotropy measurements indicate that CO rotates toward the heme plane upon dissociation, as is the case in globins. Remarkably, CO bound to the heme iron is tilted by ~30° with respect to the heme normal, which contrasts to the situation in myoglobin and in present FixLH-CO X-ray crystal structure models. This implies protein-environment-induced strain on the ligand, which is possibly at the origin of a very rapid docking-site population in a single conformation. Our observations likely explain the unusually low affinity of FixL for CO that is at the origin of the weak ligand discrimination between CO and O(2). Moreover, we observe orders of magnitude faster vibrational relaxation of dissociated CO in FixL than in globins, implying strong interactions of the ligand with the distal heme pocket environment. Finally, in the R220H FixLH mutant protein, where CO is H-bonded to a distal histidine, we demonstrate that the H-bond is maintained during photolysis. Comparison with extensively studied globin proteins unveils a surprisingly rich variety in both structural and dynamic properties of the interaction of a diatomic ligand with the ubiquitous b-type heme-proximal histidine system in different distal pockets.     

Model for proton transport coupled to protein conformational change: application to proton pumping in the bacteriorhodopsin photocycle

A modeling method is presented for protein systems in which proton transport is coupled to conformational change, as in proton pumps and in motors driven by the proton-motive force. Previously developed methods for calculating pKa values in proteins using a macroscopic dielectric model are extended beyond the equilibrium case to a master-equation model for the time evolution of the system through states defined by ionization microstate and a discrete set of conformers. The macroscopic dielectric model supplies free energy changes for changes of protonation microstate, while the method for obtaining the energetics of conformational change and the relaxation rates, the other ingredients needed for the master equation, are system dependent. The method is applied to the photoactivated proton pump, bacteriorhodopsin, using conformational free energy differences from experiment and treating relaxation rates through three adjustable parameters. The model is found to pump protons with an efficiency relatively insensitive to parameter choice over a wide range of parameter values, and most of the main features of the known photocycle from very early M to the return to the resting state are reproduced. The boundaries of these parameter ranges are such that short-range proton transfers are faster than longer-range ones, which in turn are faster than conformational changes. No relaxation rates depend on conformation. The results suggest that an "accessibility switch", while not ruled out, is not required and that vectorial proton transport can be achieved through the coupling of the energetics of ionization and conformational states.     

Simulation of loss mechanisms in organic solar cells: A description of the mesoscopic Monte Carlo technique and an evaluation of the first reaction method

In this letter we evaluate the accuracy of the first reaction method (FRM) as commonly used to reduce the computational complexity of mesoscale Monte Carlo simulations of geminate recombination and the performance of organic photovoltaic devices. A wide range of carrier mobilities, degrees of energetic disorder, and applied electric field are considered. For the ranges of energetic disorder relevant for most polyfluorene, polythiophene, and alkoxy poly(phenylene vinylene) materials used in organic photovoltaics, the geminate separation efficiency predicted by the FRM agrees with the exact model to better than 2%. We additionally comment on the effects of equilibration on low-field geminate separation efficiency, and in doing so emphasize the importance of the energy at which geminate carriers are created upon their subsequent behavior.     

Spin-forbidden ligand binding to the ferrous-heme group: ab initio and DFT studies

The potential energy surfaces (PESs) and associated energy barriers that characterize the spin-forbidden recombination reactions of the gas-phase ferrous deoxy-heme group with CO, NO, and H2O ligands have been calculated using density functional theory (DFT). The bond energy for binding of O2 has also been calculated. Extensive large basis set CCSD(T) calculations on two small models of the heme group have been used to calibrate the accuracy of different DFT functionals for treating these systems. Pure functionals are shown to overestimate the stability of the low-spin forms of the deoxy-heme model, and to overestimate the binding energy of H2O and CO, whereas hybrid functionals such as B3PW91 and B3LYP yield accurate results. Accordingly, the latter functionals have been used to explore the PESs for binding. CO binding is found to involve a significant barrier of ca. 3 kcal mol-1 due to the need to change from the deoxy-heme quintet ground state to the bound singlet state. Binding of water does not involve a barrier, but the resulting bond is weak and may be further weakened in the protein environment, which should explain why water binding is not usually observed in heme proteins such as myoglobin. NO binding involves a low barrier, which is consistent with observed rapid geminate recombination. The calculated bond energies are in good agreement with previous reported values and in fair agreement with experiment for CO and O2. The value for NO is significantly lower than the experimentally derived bond energy, suggesting that B3LYP is less accurate in this case.     

Temperature-dependent studies of NO recombination to heme and heme proteins

The rebinding kinetics of NO to the heme iron of myoglobin (Mb) is investigated as a function of temperature. Below 200 K, the transition-state enthalpy barrier associated with the fastest (approximately 10 ps) recombination phase is found to be zero and a slower geminate phase (approximately 200 ps) reveals a small enthalpic barrier (approximately 3 +/- 1 kJ/mol). Both of the kinetic rates slow slightly in the myoglobin (Mb) samples above 200 K, suggesting that a small amount of protein relaxation takes place above the solvent glass transition. When the temperature dependence of the NO recombination in Mb is studied under conditions where the distal pocket is mutated (e.g., V68W), the rebinding kinetics lack the slow phase. This is consistent with a mechanism where the slower (approximately 200 ps) kinetic phase involves transitions of the NO ligand into the distal heme pocket from a more distant site (e.g., in or near the Xe4 cavity). Comparison of the temperature-dependent NO rebinding kinetics of native Mb with that of the bare heme (PPIX) in glycerol reveals that the fast (enthalpically barrierless) NO rebinding process observed below 200 K is independent of the presence or absence of the proximal histidine ligand. In contrast, the slowing of the kinetic rates above 200 K in MbNO disappears in the absence of the protein. Generally, the data indicate that, in contrast to CO, the NO ligand binds to the heme iron through a "harpoon" mechanism where the heme iron out-of-plane conformation presents a negligible enthalpic barrier to NO rebinding. These observations strongly support a previous analysis (Srajer et al. J. Am. Chem. Soc. 1988, 110, 6656-6670) that primarily attributes the low-temperature stretched exponential rebinding of MbCO to a quenched distribution of heme geometries. A simple model, consistent with this prior analysis, is presented that explains a variety of MbNO rebinding experiments, including the dependence of the kinetic amplitudes on the pump photon energy.