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.     

Geminate rebinding in R-state hemoglobin: kinetic and computational evidence for multiple hydrophobic pockets

Biphasic geminate rebinding of CO to myoglobin upon flash photolysis has been associated to ligand distribution in hydrophobic cavities, structurally detected by time-resolved crystallography, xenon occupancy, and molecular simulations. We show that the time course of CO rebinding to human hemoglobin also exhibits a biphasic geminate rebinding when the protein is entrapped in wet nanoporous silica gel. A simple branched kinetic scheme, involving the bound state A, the primary docking site C, and a secondary binding site B was used to calculate the microscopic rates and the time-dependent population of the intermediate species. The activation enthalpies of the associated transitions were determined in the absence and presence of 80% glycerol. Potential hydrophobic docking cavities within the alpha and beta chains of hemoglobin were identified by computational modeling using xenon as a probe. A hydrophobic pocket on the distal side of the heme, corresponding to Xe4 in Mb, and a nearby site that does not have a correspondence in Mb were detected. Neither potential xenon sites on the proximal side nor a migration channel from the distal to proximal site was located. The small enthalpic barriers between states B and C are in very good agreement with the location of the xenon sites on the distal side. Furthermore, the connection between the two xenon sites is relatively open, explaining why the decreased mobility of the protein with viscosity only slightly perturbs the energetics of ligand migration between the two sites.     

Free-energy barriers in MbCO rebinding

The rebinding of CO to myoglobin (Mb) from locations around the active site is studied using a combination of molecular dynamics and stochastic simulations for native and L29F mutant Mb. The interaction between the dissociated ligand and the protein environment is described by the recently developed fluctuating three-point charge model for the CO molecule. Umbrella sampling along trajectories, previously found to sample the binding site (B) and the Xe4 pocket, is used to construct free-energy profiles for the ligand escape. On the basis of the Smoluchowski equation, the relaxation of different initial population distributions is followed in space and time. For native Mb at room temperature, the calculated rebinding times are in good agreement with experimental values and give an inner barrier of 4.3 kcal/mol between the docking site B (Mb...CO) and the A state (bound MbCO), compared to an effective barrier, Heff, of 4.5 kcal/mol and barriers into the majority conformation A1 and the minority conformation A3 of 2.4 and 4.3 kcal/mol, respectively. In the case of the L29F mutant, the free-energy surface is flatter and the dynamics is much more rapid. As was found in experiment, escape to the Xe4 pocket is facile for L29F whereas, for native Mb, the barriers to this site are larger. At lower temperatures, the rebinding dynamics is delayed by orders of magnitude also due to increased barriers between the docking sites.     

Diffusive dynamics on multidimensional rough free energy surfaces

The dynamics of processes relevant to chemistry and biophysics on rough free energy landscapes is investigated using a recently developed algorithm to solve the Smoluchowski equation. Two different processes are considered: ligand rebinding in MbCO and protein folding. For the rebinding dynamics of carbon monoxide (CO) to native myoglobin (Mb) from locations around the active site, the two-dimensional free energy surface (FES) is constructed using extensive molecular dynamics simulations. The surface describes the minima in the A state (bound MbCO), CO in the distal pocket and in the Xe4 pocket, and the transitions between these states and allows to study the diffusion of CO in detail. For the folding dynamics of protein G, a previously determined two-dimensional FES was available. To follow the diffusive dynamics on these rough free energy surfaces, the Smoluchowski equation is solved using the recently developed hierarchical discrete approximation method. From the relaxation of the initial nonequilibrium distribution, experimentally accessible quantities such as the rebinding time for CO or the folding time for protein G can be calculated. It is found that the free energy barrier for CO in the Xe4 pocket and in the distal pocket (B state) closer to the heme iron is approximately 6 kcal/mol which is considerably larger than the inner barrier which separates the bound state and the B state. For the folding of protein G, a barrier of approximately 10 kcal/mol between the unfolded and the folded state is consistent with folding times of the order of milliseconds.     

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.     

Geminate carbon monoxide rebinding to a c-type haem

A chemically modified form of cytochrome c(cyt. c), termed carboxymethyl cytochrome c(cm cyt. c), possesses a vacant sixth coordination site to the haem iron that is available to bind external ligands. We present data on the rapid flash photolysis of CO from the ferrous haem iron of cm cyt. c and describe the kinetics and spectral transitions that accompany the recombination. This was achieved using 30-femtosecond laser pulses and a white light continuum to monitor spectral transitions. Whereas the photo-dissociation quantum yield is close to 1, the yield of CO escape from the protein (the apparent quantum yield, varphi) relative to myoglobin (varphi=1) is small due to rapid geminate recombination of CO. On ligand photo-dissociation the haem undergoes a spin-state transition from low-spin ferrous CO bound to penta-coordinate high-spin. Subsequently the system reverts to the CO bound form. The data were fitted with a minimum number of exponentials using global analysis. Recombination of CO with the haem iron of cm cyt. c is multiphasic (tau=16 ps, 120 ps and 1 ns), involving three spectrally distinct components. The fraction of haem (0.11) not recombining with CO within 4 ns is similar to the value of varphi(0.12) measured on the same preparation by the "pulse method" (M. Brunori, G. Giacometti, E. Antonini and J. Wyman, Proc. Natl. Acad. Sci. USA, 1973, 70, 3141-3144, ). This implies that no further geminate recombination occurs at t>4 ns. This unusually efficient CO-haem geminate recombination indicates the sterically hindered ("caged") nature of the distal haem pocket in cm cyt. c from which it is difficult for CO to escape. The large geminate phase may be contrasted with the behaviour of myoglobin in which geminate recombination is small. This is in general agreement with the well-documented extensive structural dynamics in myoglobin that allow ligand passage, and a higher structural rigidity in cyt. c imposed by the restraints of minimising reorganisation energy for electron transfer (M. Brunori, D. Bourgeois and D. Vallone, J. Struct. Biol., 2004, 147, 223-234, ). The high pH ferrous form of cm cyt. c is a low-spin species having a lysine bound to the central iron atom of the haem (M. Brunori, M. Wilson and E. Antonini, J. Biol. Chem., 1972, 247, 6076-6081; G. Silkstone, G. Stanway, P. Brzezinski and M. Wilson, Biophys. Chem., 2002, 98, 65-77, ). This high pH (pH approximately 8) form of deoxy cm cyt. c undergoes photo-dissociation of lysine (although the proximal histidine is possible) after photo-excitation. Recombination occurs with a time constant (tau) of approximately 7 ps. This is similar to that observed for the geminate rebinding of the Met80 residue in native ferrous cyt. c(tau approximately 6 ps) following its photo-dissociation (S. Cianetti, M. Negrerie, M. Vos, J.-L. Martin and S. Kruglik, J. Am. Chem. Soc., 2004, 126, 13 932-13 933; W. Wang, X. Ye, A. Demidov, F. Rosca, T. Sjodin, W. Cao, M. Sheeran and P. Champion, J. Phys. Chem., 2000, 104, 10 789-10 801, ).     

Determination of microscopic rate constants for CO binding and migration in myoglobin encapsulated in silica gels

CO rebinding kinetics after nanosecond photolysis of myoglobin encapsulated in wet silica gels exhibits an enhanced geminate phase that allows the determination of the microscopic rate constants and the activation barriers for distinct ligand docking sites inside the protein matrix. Using a maximum entropy method, we demonstrate that the geminate phase can be well-described by a biphasic lifetime distribution, reflecting rebinding from the distal and proximal sites. Microscopic rates and activation barriers were estimated using a four-state model.     

Unusual ligand discrimination by a myoglobin reconstituted with a hydrophobic domain-linked heme

New, reconstituted horse heart myoglobins possessing a hydrophobic domain at the terminal of the two heme propionate side chains were constructed. The O2 and CO bindings for the reconstituted deoxymyoglobins were examined in detail by laser flash photolysis and stopped-flow rapid mixing techniques. The artificially created domain worked as a barrier against exogenous ligand penetration into the heme pocket, whereas the bound O2 was stabilized in the reconstituted myoglobin as well as in the native one. In contrast, the CO dissociation rate for the reconstituted myoglobin increased by 20-fold compared to the native protein, suggesting that the incorporation of the hydrophobic domain onto the heme pocket perturbs the distal-site structure of the reconstituted myoglobin. As a result, the substantial ligand selectivity for the reconstituted myoglobin significantly increases in favor of O2 over CO with the M' value (= KCO/KO2) of 0.88, whereas, to the best of our knowledge, there is no myoglobin mutant in which the O2 affinity exceeds the CO one. The present work concludes that the O2 selectivity of myoglobin over CO is markedly improved by chemically modifying the heme propionates without any mutation of the amino acid residues in the distal site.     

Orientational distribution of CO before and after photolysis of MbCO and HbCO: a determination using time-resolved polarized Mid-IR spectroscopy

The technique of time-resolved polarized mid-IR spectroscopy was used to probe the orientational distribution of carbon monoxide (CO) bound to and docked within horse myoglobin, sperm whale myoglobin, and human hemoglobin A in neutral pH solution at 283 K. An accurate determination of the orientation required that the experimentally measured polarization anisotropy be corrected for the effects of fractional photolysis in an optically thick sample. The experimental method measures the direction of the transition dipole, which is parallel to the CO bond axis when docked and nearly parallel when bound to the heme. The polarization anisotropy of bound CO is virtually the same for all protein systems investigated and is unchanging across its inhomogeneously broadened mid-IR absorption spectrum. From these results, it was concluded that the transition dipole moment of bound CO is oriented 相似文献   

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.     

Non-exponential non-Arrhenius relaxation in the course of CO rebinding to heme proteins

The dispersive transport model for relaxation of photolyzed heme proteins has been improved to take into account the coupling of the ligand-heme geminate recombination and the non-Gaussian diffusive dynamics of conformational changes in heme proteins. Contrary to the earlier deterministic version of the model, the present more rigorous formulation is based on the stochastic approach to the problem. This implies that the time evolution of protein conformations should be described in terms of the transient distribution which satisfies the Smoluchowski-type differential equation with a time-dependent diffusion coefficient. The obtained analytical solution of this equation enables us to relate main kinetic parameters of the geminate recombination and quantities characterizing the ligand-heme interaction. The derived expressions demonstrate that the reaction barrier shifts with time towards higher values following the near-stretched exponential behavior in agreement with experiment. Such a behavior is governed by the non-exponential non-Arrhenius conformational relaxation. The latter process can be identified by the characteristics “footprint” left on the experimental rebinding curve and is shown to be responsible for some kinetically different phases of the ligand-heme geminate recombination observed within distinct temperature ranges.     

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.     

O2 and CO binding properties of artificial hemoproteins formed by complexing iron protoporphyrin IX with human serum albumin mutants

The binding properties of O2 and CO to recombinant human serum albumin (rHSA) mutants with a prosthetic heme group have been physicochemically and kinetically characterized. Iron(III) protoporphyrin IX (hemin) is bound in subdomain IB of wild-type rHSA [rHSA(wt)] with weak axial coordination by Tyr-161. The reduced ferrous rHSA(wt)-heme under an Ar atmosphere exists in an unusual mixture of four- and five-coordinate complexes and is immediately autoxidized by O2. To confer O2 binding capability on this naturally occurring hemoprotein, a proximal histidine was introduced into position Ile-142 or Leu-185 by site-directed mutagenesis. A single mutant (I142H) and three double mutants (I142H/Y161L, I142H/Y161F, and Y161L/L185H) were prepared. Both rHSA(I142H/Y161L)-heme and rHSA(I142H/Y161F)-heme formed ferrous five-N-coordinate high-spin complexes with axial ligation of His-142 under an Ar atmosphere. These artificial hemoproteins bind O2 at room temperature. Mutation at the other side of the porphyrin, Y161L/L185H, also allowed O2 binding to the heme. In contrast, the single mutant rHSA(I142H)-heme could not bind O2, suggesting that removal of Y161 is necessary to confer reversible O2 binding. Laser flash photolysis experiments showed that the kinetics of CO recombination with the rHSA(mutant)-heme were biphasic, whereas O2 rebinding exhibited monophasic kinetics. This could be due to the two different geometries of the axial imidazole coordination arising from the two orientations of the porphyrin plane in the heme pocket. The O2 binding affinities of the rHSA(mutant)-heme were significantly lower than those of hemoglobin and myoglobin, principally due to the high O2 dissociation rates. Changing Leu-161 to Phe-161 at the distal side increased the association rates of both O2 and CO, which resulted in enhanced binding affinity.     

The reactivity with CO of AHb1 and AHb2 from Arabidopsis thaliana is controlled by the distal HisE7 and internal hydrophobic cavities

The nonsymbiotic hemoglobins, AHb1 and AHb2, have recently been isolated from Arabidopsis thaliana. Using steady-state and time-resolved spectroscopic methods, we show that Fe2+ AHb1 contains a mixture of penta- and hexacoordinated heme, while Fe2+ AHb2 is fully hexacoordinated. In the CO complexes, polar interactions and H-bonds with the ligand are stronger for AHb1 than for AHb2. The ligand binding kinetics are substantially different, reflecting the distribution between the penta- and hexacoordinated species, and indicate that protein dynamics and ligand migration pathways are very specific for each of the two proteins. In particular, a very small, non-exponential geminate rebinding observed in AHb1 suggests that the distal heme cavity is connected with the exterior by a relatively open channel. The large, temperature-dependent geminate rebinding observed for AHb2 implies a major role of protein dynamics in the ligand migration from the distal cavity to the solvent. The structures of AHb1 and AHb2, modeled on the basis of the homologous rice hemoglobin, exhibit a different cavity system that is fully compatible with the observed ligand binding kinetics. Overall, these kinetic and structural data are consistent with the putative NO-dioxygenase activity previously attributed to AHb1, whereas the role of AHb2 remains elusive.     

Exciplex Formation in Lipid-bound Escherichia coli Flavohemoglobin

Flavohemoglobins have the particular capability of binding unsaturated and cyclopropanated fatty acids as free acids or phospholipids. Fatty acid binding to the ferric heme results in a weak but direct bonding interaction. Ferrous and ferric protein, in presence or absence of a bound lipid molecule, have been characterized by transient absorption spectroscopy. Measurements have been also carried out both on the ferrous deoxygenated and on the CO bound protein to investigate possible long-range interaction between the lipid acyl chain moiety and the ferrous heme. After excitation of the deoxygenated derivatives the relaxation process reveals a slow dynamics (350 ps) in lipid-bound protein but is not observed in the lipid-free protein. The latter feature and the presence of an extra contribution in the absorption spectrum, indicates that the interaction of iron heme with the acyl chain moiety occurs only in the excited electronic state and not in the ground electronic state. Data analysis highlights the formation of a charge-transfer complex in which the iron ion of the lipid-bound protein in the expanded electronic excited state, possibly represented by a high spin Fe III intermediate, is able to bind to the sixth coordination ligand placed at a distance of at 3.5 Å from the iron. A very small nanosecond geminate rebinding is observed for CO adduct in lipid-free but not in the lipid-bound protein. The presence of the lipid thus appears to inhibit the mobility of CO in the heme pocket.     

Evidence for two geminate rebinding states following laser photolysis of R state hemoglobin encapsulated in wet silica gels

In this letter we report the first experimental evidence for CO rebinding to human hemoglobin from multiple geminate states. The analysis of the rebinding kinetics using a maximum entropy method allowed the identification of two distinct rebinding states within the protein matrix, which become populated under conditions of increased viscosity in a silica gel at high glycerol concentration. Our findings suggest the presence of at least two distinct docking sites for the photolyzed ligand. Assuming a minimal four-state model, we estimate the microscopic rates and the activation energies for the elementary processes.     

Iron hemiporphycene as a functional prosthetic group for myoglobin

The iron complex of hemiporphycene, a molecular hybrid of porphyrin with porphycene, was incorporated into the apomyoglobin pocket to examine ligand binding ability of the iron atom in the novel porphyrinoid. Apomyoglobin was successfully coupled with a stoichiometric amount of ferric hemiporphycene to afford the reconstituted myoglobin equipped with the iron coordination structure of native protein. Cyanide, imidazole, and fluoride coordinated to the ferric protein with affinities comparable with those for native myoglobin. The ferrous myoglobin was functionally active to bind O(2) and CO reversibly at pH 7.4 and 20 degrees C. The O(2) affinity is 12-fold higher than that of native myoglobin while the CO affinity is slightly lower, suggesting decreased discrimination between O(2) and CO in the heme pocket. The functional anomaly was interpreted to reflect increased sigma-bonding character in the Fe(II)-O(2) bond. In contrast with 6-coordinate native NO protein, the NO myoglobin containing ferrous hemiporphycene is in a mixed 5- and 6-coordinate state. This observation suggests that the in-plane configuration of the iron atom in hemiporphycene is destabilized by NO. Influence of the core deformation was also detected with both the infrared absorption for the ferrous CO derivative and electron paramagnetic resonance for ferric imidazole complex. Anomalies in the ferric and ferrous derivatives were ascribed to the modified iron-N(pyrrole) interactions in the asymmetric metallo core of hemiporphycene.     

An electrostatic model for the frequency shifts in the carbonmonoxy stretching band of myoglobin: correlation of hydrogen bonding and the stark tuning rate

The effect of internal and applied external electric fields on the vibrational stretching frequency for bound CO (nu(CO)) in myoglobin mutants was studied using density functional theory. Geometry optimization and frequency calculations were carried out for an imidazole-iron-porphine-carbonmonoxy adduct with various small molecule hydrogen-bonding groups. Over 70 vibrational frequency calculations of different model geometries and hydrogen-bonding groups were compared to derive overall trends in the C-O stretching frequency (nu(CO)) in terms of the C-O bond length and Mulliken charge. Simple linear functions were derived to predict the Stark tuning rate using an approach analogous to the vibronic theory of activation.(1) Potential energy calculations show that the strongest interaction occurs for C-H or N-H hydrogen bonding nearly perpendicular to the Fe-C-O bond axis. The calculated frequencies are compared to the structural data available from 18 myoglobin crystal structures, supporting the hypothesis that the vast majority of hydrogen-bonding interactions with CO occur from the side, rather than the end, of the bound CO ligand. The nu(CO) frequency shifts agree well with experimental frequency shifts for multiple bands, known as A states, and site-directed mutations in the distal pocket of myoglobin. The model calculations quantitatively explain electrostatic effects in terms of specific hydrogen-bonding interactions with bound CO in heme proteins.     

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.     

CO cage recombination in hemoglobin: Picosecond photolysis and nanosecond observation

Carboxy hemoglobin in aqueous solution was photodissociated by laser pulses of 30 ps at 532 nm. Kinetic studies show that only upon complete photodissociation can the pure CO geminate binding process be revealed. The protein region of the iron cage, where the geminate ligand diffuses before reaching the heme, is smaller than a single subunit and larger than the heme pocket.