Effects of Distal Pocket Mutations on the Geminate Recombination of NO with Leghemoglobin on the Picosecond Time Scale

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.

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.

Role of Heme Iron Coordination and Protein Structure in the Dynamics and Geminate Rebinding of Nitric Oxide to the H93G Myoglobin Mutant

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

Ligand Dynamics in Myoglobin: Calculation of Infrared Spectra for Photodissociated NO

We present molecular dynamics simulations of the photodissociated state of MbNO performed at 300 K using a fluctuating charge model for the nitric oxide (NO) ligand. After dissociation, NO is observed to remain mainly in the centre of the distal haem pocket, although some movement towards the primary docking site and the xenon-4 pocket can be seen. We calculate the NO infrared spectrum for the photodissociated ligand within the haem pocket and find a narrow peak in the range 1915-1922 cm(-1). The resulting blue shift of 1 to 8 cm(-1) compared to gas-phase NO is much smaller than the red shifts calculated and observed for carbon monoxide (CO) in Mb. A small splitting, due to NO in the xenon-4 pocket, is also observed. At lower temperatures, the spectra and conformational space explored by the ligand remain largely unchanged, but the electrostatic interactions with residue His64 become increasingly significant in determining the details of the ligand orientation within the distal haem pocket. The investigation of the effect of the L29F mutation reveals significant differences between the behaviour of NO and that of CO, and suggests a coupling between the ligand and the protein dynamics due to the different ligand dipole moments.

Tyrosine B10 inhibits stabilization of bound carbon monoxide and oxygen in soybean leghemoglobin

Detailed comparisons of the carbon monoxide FTIR spectra and ligand-binding properties of a library of E7, E11, and B10 mutants indicate significant differences in the role of electrostatic interactions in the distal pockets of wild-type sperm whale myoglobin and soybean leghemoglobin. In myoglobin, strong hydrogen bonds from several closely related conformations of the distal histidine (His(E7)) side chain preferentially stabilize bound oxygen. In leghemoglobin, the imidazole side chain of His(E7) is confined to a single conformation, which only weakly hydrogen bonds to bound ligands. The phenol side chain of Tyr(B10) appears to "fix" the position of His(E7), probably by donating a hydrogen bond to the Ndelta atom of the imidazole side chain. The proximal pocket of leghemoglobin is designed to favor strong coordination bonds between the heme iron and axial ligands. Thus, high oxygen affinity in leghemoglobin is established by a favorable staggered geometry of the proximal histidine. The interaction between His(E7) and Tyr(B10) prevents overstabilization of bound oxygen. If hydrogen bonding from His(E7) were as strong as it is in mammalian myoglobin, the resultant ultrahigh affinity of leghemoglobin would prevent oxygen transport in root nodules.