Viscosity-dependent relaxation significantly modulates the kinetics of CO recombination in the truncated hemoglobin TrHbN from Mycobacterium tuberculosis

Kinetic traces were generated for the nanosecond and slower rebinding of photodissociated CO to trHbN in solution and in porous sol-gel matrices as a function of viscosity, conformation, and mutation. TrHbN is one of the two truncated hemoglobins from Mycobacterium tuberculosis. The kinetic traces were analyzed in terms of three distinct phases. These three phases are ascribed to rebinding: (i) from the distal heme pocket, (ii) from the adjacent apolar tunnel prior to conformational relaxation, and (iii) from the apolar tunnel subsequent to conformational relaxation. The fractional content of each of these phases was shown to be a function of the viscosity and, in the case of the sol-gel-encapsulated samples, sample preparation history. The observed kinetic patterns support a model consisting of the following elements: (i) the viscosity and conformation-sensitive dynamics of the Tyr(B10) side chain facilitate diffusion of the dissociated ligand from the distal heme pocket into the adjacent tunnel; (ii) the distal heme pocket architecture determines ligand access from the tunnel back to the heme iron; (iii) the distal heme pocket architecture is governed by a ligand-dependent hydrogen bonding network that limits the range of accessible side chain positions; and (iv) the apolar tunnel linking the heme site to the solvent biases the competition between water and ligand for occupancy of the vacated polar distal heme pocket greatly toward the nonpolar ligand. Implications of these finding with respect to biological function are discussed.     

Heme-ligand tunneling in group I truncated hemoglobins

Truncated hemoglobins (trHbs) are small hemoproteins forming a separate cluster within the hemoglobin superfamily; their functional roles in bacteria, plants, and unicellular eukaryotes are marginally understood. Crystallographic investigations have shown that the trHb fold (a two-on-two alpha-helical sandwich related to the globin fold) hosts a protein matrix tunnel system offering a potential path for ligand diffusion to the heme distal site. The tunnel topology is conserved in group I trHbs, although with modulation of its size/structure. Here, we present a crystallographic investigation on trHbs from Mycobacterium tuberculosis, Chlamydomonas eugametos, and Paramecium caudatum, showing that treatment of trHb crystals under xenon pressure leads to binding of xenon atoms at specific (conserved) sites along the protein matrix tunnel. The crystallographic results are in keeping with data from molecular dynamics simulations, where a dioxygen molecule is left free to diffuse within the protein matrix. Modulation of xenon binding over four main sites is related to the structural properties of the tunnel system in the three trHbs and may be connected to their functional roles. In a parallel crystallographic investigation on M. tuberculosis trHbN, we show that butyl isocyanide also binds within the apolar tunnel, in excellent agreement with concepts derived from the xenon binding experiments. These results, together with recent data on atypical CO rebinding kinetics to group I trHbs, underline the potential role of the tunnel system in supporting diffusion, but also accumulation in multiple copies, of low polarity ligands/molecules within group I trHbs.     

The position 68(E11) side chain in myoglobin regulates ligand capture, bond formation with heme iron, and internal movement into the xenon cavities

After photodissociation, ligand rebinding to myoglobin exhibits complex kinetic patterns associated with multiple first-order geminate recombination processes occurring within the protein and a simpler bimolecular phase representing second-order ligand rebinding from the solvent. A smooth transition from cryogenic-like to solution phase properties can be obtained by using a combination of sol-gel encapsulation, addition of glycerol as a bathing medium, and temperature tuning (-15 --> 65 degrees C). This approach was applied to a series of double mutants, myoglobin CO (H64L/V68X, where X = Ala, Val, Leu, Asn, and Phe), which were designed to examine the contributions of the position 68(E11) side chain to the appearance and disappearance of internal rebinding phases in the absence of steric and polar interactions with the distal histidine. Based on the effects of viscosity, temperature, and the stereochemistry of the E11 side chain, the three major phases, B --> A, C --> A, and D --> A, can be assigned, respectively, to ligand rebinding from the following: (i) the distal heme pocket, (ii) the xenon cavities prior to large amplitude side chain conformational relaxation, and (iii) the xenon cavities after significant conformational relaxation of the position 68(E11) side chain. The relative amplitudes of the B --> A and C --> A phases depend markedly on the size and shape of the E11 side chain, which regulates sterically both ligand return to the heme iron atom and ligand migration to the xenon cavities. The internal xenon cavities provide a transient docking site that allows side chain relaxations and the entry of water into the vacated distal pocket, which in turn slows ligand recombination markedly.     

Myoglobin mutants giving the largest geminate yield in CO rebinding in the nanosecond time domain.

We have measured the rebinding of carbon monoxide (CO) to some distal mutants of myoglobin (Mb) in the time range from 10(-8) to 10(-1) s by flash photolysis, in which the photodissociated CO rebinds to the heme iron without escaping to the solvent water from the protein matrix. We have found that the double mutants [His64-->Val/Val68-->Thr (H64V/V68T) and His64-->Val/Val68-->Ser (H64V/V68S)] have an extremely large geminate yield (70-80%) in water at 5 degreesC, in contrast to the 7% of the geminate yield of wild-type Mb. The CO geminate yields for these two mutants are the largest in those of Mb mutants reported so far, showing that the two mutants have a unique heme environment that favors CO geminate rebinding. Comparing the crystal structures and 1H-NMR and vibrational spectral data of H64V/V68T and H64V/V68S with those of other mutants, we discuss factors that may control the nanosecond geminate CO rebinding and CO migration in the protein matrix.     

Analysis of the kinetic barriers for ligand binding to sperm whale myoglobin using site-directed mutagenesis and laser photolysis techniques

Time courses for NO, O2, CO, methyl and ethyl isocyanide rebinding to native and mutant sperm whale myoglobins were measured at 20 degrees C following 17-ns and 35-ps laser excitation pulses. His64 (E7) was replaced with Gly, Val, Leu, Phe, and Gln, and Val68 (E11) was replaced with Ala, Ile, and Phe. For both NO and O2, the effective picosecond quantum yield of unliganded geminate intermediates was roughly 0.2 and independent of the amino acids at positions 64 and 68. Geminate recombination of NO was very rapid; 90% rebinding occurred within 0.5-1.0 ns for all of the myoglobins examined; and except for the Gly64 and Ile68 mutants, the fitted recombination rate parameters were little influenced by the size and polarity of the amino acid at position 64 and the size of the residue at position 68. The rates of NO recombination and ligand movement away from the iron atom in the Gly64 mutant increased 3-4-fold relative to native myoglobin. For Ile68 myoglobin, the first geminate rate constant for NO rebinding decreased approximately 6-fold, from 2.3 x 10(10) s-1 for native myoglobin to 3.8 x 10(9) s-1 for the mutant. No picosecond rebinding processes were observed for O2, CO, and isocyanide rebinding to native and mutant myoglobins; all of the observed geminate rate constants were less than or equal to 3 x 10(8) s-1. The rebinding time courses for these ligands were analyzed in terms of a two-step consecutive reaction scheme, with an outer kinetic barrier representing ligand movement into and out of the protein and an inner barrier representing binding to the heme iron atom by ligand occupying the distal portion of the heme pocket. Substitution of apolar amino acids for His64 decreased the absolute free energies of the outer and inner kinetic barriers and the well for non-covalently bound O2 and CO by 1 to 1.5 kcal/mol, regardless of size. In contrast, the His64 to Gln mutation caused little change in the barrier heights for all ligands, showing that the polar nature of His64 inhibits both the bimolecular rate of ligand entry into myoglobin and the unimolecular rate of binding to the iron atom from within the protein. Increasing the size of the position 68(E11) residue in the series Ala to Val (native) to Ile caused little change in the rate of O2 migration into myoglobin or the equilibrium constant for noncovalent binding but did decrease the unimolecular rate for iron-O2 bond formation.(ABSTRACT TRUNCATED AT 400 WORDS)     

Geminate recombination of nitric oxide to endothelial nitric-oxide synthase and mechanistic implications.

The nitric-oxide synthase (NOS) catalyzes the oxidation of L-arginine to L-citrulline and NO through consumption of oxygen bound to the heme. Because NO is produced close to the heme and may bind to it, its subsequent role in a regulatory mechanism should be scrutinized. We therefore examined the kinetics of NO rebinding after photodissociation in the heme pocket of human endothelial NOS by means of time-resolved absorption spectroscopy. We show that geminate recombination of NO indeed occurs and that this process is strongly modulated by L-Arg. This NO rebinding occurs in a multiphasic fashion and spans over 3 orders of magnitude. In both ferric and ferrous states of the heme, a fast nonexponential picosecond geminate rebinding first takes place followed by a slower nanosecond phase. The rates of both phases decreased, whereas their relative amplitudes are changed by the presence of L-Arg; the overall effect is a slow down of NO rebinding. For the isolated oxygenase domain, the picosecond rate is unchanged, but the relative amplitude of the nanosecond binding decreased. We assigned the nanosecond kinetic component to the rebinding of NO that is still located in the protein core but not in the heme pocket. The implications for a mechanism of regulation involving NO binding are discussed.     

Ligand migration in the apolar tunnel of Cerebratulus lacteus mini-hemoglobin

The large apolar tunnel traversing the mini-hemoglobin from Cerebratulus lacteus (CerHb) has been examined by x-ray crystallography, ligand binding kinetics, and molecular dynamic simulations. The addition of 10 atm of xenon causes loss of diffraction in wild-type (wt) CerHbO(2) crystals, but Leu-86(G12)Ala CerHbO(2), which has an increased tunnel volume, stably accommodates two discrete xenon atoms: one adjacent to Leu-86(G12) and another near Ala-55(E18). Molecular dynamics simulations of ligand migration in wt CerHb show a low energy pathway through the apolar tunnel when Leu or Ala, but not Phe or Trp, is present at the 86(G12) position. The addition of 10-15 atm of xenon to solutions of wt CerHbCO and L86A CerHbCO causes 2-3-fold increases in the fraction of geminate ligand recombination, indicating that the bound xenon blocks CO escape. This idea was confirmed by L86F and L86W mutations, which cause even larger increases in the fraction of geminate CO rebinding, 2-5-fold decreases in the bimolecular rate constants for ligand entry, and large increases in the computed energy barriers for ligand movement through the apolar tunnel. Both the addition of xenon to the L86A mutant and oxidation of wt CerHb heme iron cause the appearance of an out Gln-44(E7) conformer, in which the amide side chain points out toward the solvent and appears to lower the barrier for ligand escape through the E7 gate. However, the observed kinetics suggest little entry and escape (≤ 25%) through the E7 pathway, presumably because the in Gln-44(E7) conformer is thermodynamically favored.     

Investigation of laser-induced long-lived states of photolyzed MbCO

We present evidence from resonance Raman and absorption measurements that the extended exposure of MbCO to CW laser light at low temperatures alters the CO rebinding kinetics and leads to a significantly increased population of very long lived states of photolyzed MbCO. This optical "pumping" process is observed for samples frozen in both aqueous buffer and glycerol/buffer and exhibits power law behavior with a very weak temperature dependence. A comparison of the nonexponential rebinding kinetics of CO molecules from the pumped states with the rebinding observed in flash photolysis experiments suggests that the pumped states are distinct geminate states, not observed in flash photolysis experiments. Thus, a four-state model, with two geminate states, is implicated for MbCO. Pumped states may represent "separated geminate pair" states with the CO molecule still in the heme pocket or possibly trapped within a cavity on its way through the protein matrix, consistent with molecular dynamics simulations. The possibility of significant deoxyheme relaxation from a less domed to a more domed configuration, as a result of the multiple photolysis events associated with the pumping process, is also explored. However, the small changes observed in the Soret band line shape and position subsequent to pumping at T less than 180 K tend to rule out this explanation for the pumping process. Since the yield for creating a pumped state is small (e.g., less than 10(-7) for T greater than 100 K), pumping can be observed only after extended illumination and is absent in flash photolysis measurements, even after multiple flashes. At higher temperatures (T greater than 180 K), the escape of the CO molecule to the solvent is observed. Our data are consistent with a "phase transition" of the protein that is coupled to the surrounding matrix. The protein fluctuations are quenched below approximately 185 K for a solvent composed of 70% glycerol and below approximately 260 K for aqueous buffer. We also present the first large amplitude measurements of CO rebinding from the protein exterior, observed below 200 K after freezing the sample under laser illumination.     

Structural bases for heme binding and diatomic ligand recognition in truncated hemoglobins

Truncated hemoglobins (trHbs) are low-molecular-weight oxygen-binding heme-proteins distributed in eubacteria, cyanobacteria, unicellular eukaryotes, and in higher plants, constituting a distinct group within the hemoglobin (Hb) superfamily. TrHbs display amino acid sequences 20-40 residues shorter than classical (non)vertebrate Hbs and myoglobins, to which they are scarcely related by sequence similarity. The trHb tertiary structure is based on a 2-on-2 alpha-helical sandwich, which represents a striking editing of the highly conserved 3-on-3 alpha-helical globin fold, achieved through deletion/truncation of alpha-helices and specific residue substitutions. Despite their 'minimal' polypeptide chain span, trHbs display an inner tunnel/cavity system held to support ligand diffusion to/from the heme distal pocket, accumulation of heme ligands within the protein matrix, and/or multiligand reactions. Moreover, trHbs bind and effectively stabilize the heme and recognize diatomic ligands (i.e., O2, CO, NO, and cyanide), albeit with varying thermodynamic and kinetic parameters. Here, structural bases for heme binding and diatomic ligand recognition by trHbs are reviewed.     

Picosecond geminate recombination of nitrosylmyoglobins

The kinetics of NO geminate recombination to sperm whale and elephant myoglobins has been studied on the picosecond time scale using an amplified colliding-pulse mode-locked ring dye laser. The dynamics of ligand rebinding are shown to be affected by the distal structure of the protein surrounding the heme pocket.     

Proximal and distal influences on ligand binding kinetics in microperoxidase and heme model compounds

We use laser flash photolysis and time-resolved Raman spectroscopy of CO-bound heme complexes to study proximal and distal influences on ligand rebinding kinetics. We report kinetics of CO rebinding to microperoxidase (MP) and 2-methylimidazole ligated Fe protoporphyrin IX in the 10 ns to 10 ms time window. We also report CO rebinding kinetics of MP in the 150 fs to 140 ps time window. For dilute, micelle-encapsulated (monodisperse) samples of MP, we do not observe the large amplitude geminate decay at approximately 100 ps previously reported in time-resolved IR measurements on highly concentrated samples [Lim, M., Jackson, T. A., and Anfinrud, P. A. (1997) J. Biol. Inorg. Chem. 2, 531-536]. However, for high concentration aggregated samples, we do observe the large amplitude picosecond CO geminate rebinding and find that it is correlated with the absence of the iron-histidine vibrational mode in the time-resolved Raman spectrum. On the basis of these results, the energetic significance of a putative distal pocket CO docking site proposed by Lim et al. may need to be reconsidered. Finally, when high concentration samples of native myoglobin (Mb) were studied as a control, an analogous increase in the geminate rebinding kinetics was not observed. This verifies that studies of Mb under dilute conditions are applicable to the more concentrated regime found in the cellular milieu.     

Structural dynamics in the active site of murine neuroglobin and its effects on ligand binding

We have examined the effects of active site residues on ligand binding to the heme iron of mouse neuroglobin using steady-state and time-resolved visible spectroscopy. Absorption spectra of the native protein, mutants H64L and K67L and double mutant H64L/K67L were recorded for the ferric and ferrous states over a wide pH range (pH 4-11), which allowed us to identify a number of different species with different ligands at the sixth coordination, to characterize their spectroscopic properties, and to determine the pK values of active site residues. In flash photolysis experiments on CO-ligated samples, reaction intermediates and the competition of ligands for the sixth coordination were studied. These data provide insights into structural changes in the active site and the role of the key residues His64 and Lys67. His64 interferes with exogenous ligand access to the heme iron. Lys67 sequesters the distal pocket from the solvent. The heme iron is very reactive, as inferred from the fast ligand binding kinetics and the ability to bind water or hydroxyl ligands to the ferrous heme. Fast bond formation favors geminate rebinding; yet the large fraction of bimolecular rebinding observed in the kinetics implies that ligand escape from the distal pocket is highly efficient. Even slight pH variations cause pronounced changes in the association rate of exogenous ligands near physiological pH, which may be important in functional processes.     

Ligand binding to heme proteins. VI. Interconversion of taxonomic substates in carbonmonoxymyoglobin.

The kinetic properties of the three taxonomic A substates of sperm whale carbonmonoxy myoglobin in 75% glycerol/buffer are studied by flash photolysis with monitoring in the infrared stretch bands of bound CO at nu(A0) approximately 1967 cm-1, nu(A1) approximately 1947 cm-1, and nu(A3) approximately 1929 cm-1 between 60 and 300 K. Below 160 K the photodissociated CO rebinds from the heme pocket, no interconversion among the A substates is observed, and rebinding in each A substate is nonexponential in time and described by a different temperature-independent distribution of enthalpy barriers with a different preexponential. Measurements in the electronic bands, e.g., the Soret, contain contributions of all three A substates and can, therefore, be only approximately modeled with a single enthalpy distribution and a single preexponential. The bond formation step at the heme is fastest for the A0 substate, intermediate for the A1 substate, and slowest for A3. Rebinding between 200 and 300 K displays several processes, including geminate rebinding, rebinding after ligand escape to the solvent, and interconversion among the A substates. Different kinetics are measured in each of the A bands for times shorter than the characteristic time of fluctuations among the A substates. At longer times, fluctuational averaging yields the same kinetics in all three A substates. The interconversion rates between A1 and A3 are determined from the time when the scaled kinetic traces of the two substates merge. Fluctuations between A1 and A3 are much faster than those between A0 and either A1 or A3, so A1 and A3 appear as one kinetic species in the exchange with A0. The maximum-entropy method is used to extract the distribution of rate coefficients for the interconversion process A0 <--> A1 + A3 from the flash photolysis data. The temperature dependencies of the A substate interconversion processes are fitted with a non-Arrhenius expression similar to that used to describe relaxation processes in glasses. At 300 K the interconversion time for A0 <--> A1 + A3 is 10 microseconds, and extrapolation yields approximately 1 ns for A1 <--> A3. The pronounced kinetic differences imply different structural rearrangements. Crystallographic data support this conclusion: They show that formation of the A0 substate involves a major change of the protein structure; the distal histidine rotates about the C(alpha)-C(beta) bond, and its imidazole sidechain swings out of the heme pocket into the solvent, whereas it remains in the heme pocket in the A1 <--> A3 interconversion. The fast A1 <--> A3 exchange is inconsistent with structural models that involve differences in the protonation between A1 and A3.     

Ligand binding to truncated hemoglobin N from Mycobacterium tuberculosis is strongly modulated by the interplay between the distal heme pocket residues and internal water

The survival of Mycobacterium tuberculosis requires detoxification of host *NO. Oxygenated Mycobacterium tuberculosis truncated hemoglobin N catalyzes the rapid oxidation of nitric oxide to innocuous nitrate with a second-order rate constant (k'(NOD) approximately 745 x 10(6) m(-1) x s(-1)), which is approximately 15-fold faster than the reaction of horse heart myoglobin. We ask what aspects of structure and/or dynamics give rise to this enhanced reactivity. A first step is to expose what controls ligand/substrate binding to the heme. We present evidence that the main barrier to ligand binding to deoxy-truncated hemoglobin N (deoxy-trHbN) is the displacement of a distal cavity water molecule, which is mainly stabilized by residue Tyr(B10) but not coordinated to the heme iron. As observed in the Tyr(B10)/Gln(E11) apolar mutants, once this kinetic barrier is lowered, CO and O(2) binding is very rapid with rates approaching 1-2 x 10(9) m(-1) x s(-1). These large values almost certainly represent the upper limit for ligand binding to a heme protein and also indicate that the iron atom in trHbN is highly reactive. Kinetic measurements on the photoproduct of the *NO derivative of met-trHbN, where both the *NO and water can be directly followed, revealed that water rebinding is quite fast (approximately 1.49 x 10(8) s(-1)) and is responsible for the low geminate yield in trHbN. Molecular dynamics simulations, performed with trHbN and its distal mutants, indicated that in the absence of a distal water molecule, ligand access to the heme iron is not hindered. They also showed that a water molecule is stabilized next to the heme iron through hydrogen-bonding with Tyr(B10) and Gln(E11).     

Ligand binding to heme proteins: II. Transitions in the heme pocket of myoglobin.

Phenomena occurring in the heme pocket after photolysis of carbonmonoxymyoglobin (MbCO) below about 100 K are investigated using temperature-derivative spectroscopy of the infrared absorption bands of CO. MbCO exists in three conformations (A substrates) that are distinguished by the stretch bands of the bound CO. We establish connections among the A substates and the substates of the photoproduct (B substates) using Fourier-transform infrared spectroscopy together with kinetic experiments on MbCO solution samples at different pH and on orthorhombic crystals. There is no one-to-one mapping between the A and B substates; in some cases, more than one B substate corresponds to a particular A substate. Rebinding is not simply a reversal of dissociation; transitions between B substates occur before rebinding. We measure the nonequilibrium populations of the B substates after photolysis below 25 K and determine the kinetics of B substate transitions leading to equilibrium. Transitions between B substates occur even at 4 K, whereas those between A substates have only been observed above about 160 K. The transitions between the B substates are nonexponential in time, providing evidence for a distribution of substates. The temperature dependence of the B substate transitions implies that they occur mainly by quantum-mechanical tunneling below 10 K. Taken together, the observations suggest that the transitions between the B substates within the same A substate reflect motions of the CO in the heme pocket and not conformational changes. Geminate rebinding of CO to Mb, monitored in the Soret band, depends on pH. Observation of geminate rebinding to the A substates in the infrared indicates that the pH dependence results from a population shift among the substates and not from a change of the rebinding to an individual A substate.     

Investigations of photolysis and rebinding kinetics in myoglobin using proximal ligand replacements

We use laser flash photolysis and time-resolved Raman spectroscopy of CO-bound H93G myoglobin (Mb) mutants to study the influence of the proximal ligand on the CO rebinding kinetics. In H93G mutants, where the proximal linkage with the protein is eliminated and the heme can bind exogenous ligands (e.g., imidazole, 4-bromoimidazole, pyridine, or dibromopyridine), we observe significant effects on the CO rebinding kinetics in the 10 ns to 10 ms time window. Resonance Raman spectra of the various H93G Mb complexes are also presented to aid in the interpretation of the kinetic results. For CO-bound H93G(dibromopyridine), we observe a rapid large-amplitude geminate phase with a fundamental CO rebinding rate that is approximately 45 times faster than for wild-type MbCO at 293 K. The absence of an iron proximal ligand vibrational mode in the 10 ns photoproduct Raman spectrum of CO-bound H93G(dibromopyridine) supports the hypothesis that proximal ligation has a significant influence on the kinetics of diatomic ligand binding to the heme.     

Viscosity-dependent energy barriers and equilibrium conformational fluctuations in oxygen recombination with hemerythrin

The recombination kinetics of photo-dissociated oxyhemerythrin (Sipunculus nudus) have been investigated between 298 K and 90 K. Fast geminate recombinations compete with oxygen escape into the solvent, from which a subsequent slower bimolecular rebinding takes place. In phosphate buffer (pH 7.7) at 278 K, the fast and slow processes are exponential and have comparable amplitudes. Whereas the oxygen escape rate rapidly decreases upon increasing the viscosity, the inward rate from the solvent is found to be independent of viscosity, up to about 50 cP (50 mPa.s). The data suggest that a Brownian-motion-driven displacement of one or several side-chain residues is implied in oxygen escape from within the protein and also that hemerythrin undergoes a conformational change in the deoxy state. At higher viscosities and lower temperature only the geminate phase is observed and the kinetics progressively depart from an exponential. Below about 130 K, the kinetics resemble those reported in the literature for heme proteins. They are consistent with a temperature-independent non-equilibrium frozen distribution of conformational substates. However, between 190 K and 130 K, the profile of the kinetics is invariant on a log/log plot and the results simply differ by a translation along the log t axis. It is shown that this property is expected only for a temperature-dependent distribution of substates in a Boltzmann equilibrium. From room temperature, where rebinding is exponential, down to the 'freezing' temperature, the geminate recombinations display a variety of kinetic laws. It can be shown, however, that for a broad class of substate distributions, the initial slope of the kinetic plot follows an Arrhenius relationship. The activation energy is equal to that of the exponential rate constant measured at high temperature. This result establishes the conditions under which protein data obtained from low-temperature kinetics can be extrapolated to physiological temperature.     

Ligand migration and hexacoordination in type 1 non-symbiotic rice hemoglobin

Type 1 non-symbiotic rice hemoglobin (rHb1) shows bis-histidyl heme hexacoordination and is capable of binding diatomic ligands reversibly. The biological function is as yet unclear, but the high oxygen affinity makes it unlikely to be involved in oxygen transport. In order to gain insight into possible physiological roles, we have studied CO rebinding kinetics after laser flash photolysis of rHb1 in solution and encapsulated in silica gel. CO rebinding to wt rHb1 in solution occurs through a fast geminate phase with no sign of rebinding from internal docking sites. Encapsulation in silica gel enhances migration to internal cavities. Site-directed mutagenesis of FB10, a residue known to have a key role in the regulation of hexacoordination and ligand affinity, resulted in substantial effects on the rebinding kinetics, partly inhibiting ligand exit to the solvent, enhancing geminate rebinding and enabling ligand migration within the internal cavities. The mutation of HE7, one of the histidyl residues involved in the hexacoordination, prevents hexacoordination, as expected, but also exposes ligand migration through a complex system of cavities. This article is part of a Special Issue entitled: Protein Dynamics: Experimental and Computational Approaches.     

Unraveling the molecular basis for ligand binding in truncated hemoglobins: The trHbO Bacillus subtilis case

Truncated hemoglobins (trHbs) are heme proteins present in bacteria, unicellular eukaryotes, and higher plants. Their tertiary structure consists in a 2‐over‐2 helical sandwich, which display typically an inner tunnel/cavity system for ligand migration and/or storage. The microorganism Bacillus subtilis contains a peculiar trHb, which does not show an evident tunnel/cavity system connecting the protein active site with the solvent, and exhibits anyway a very high oxygen association rate. Moreover, resonant Raman results of CO bound protein, showed that a complex hydrogen bond network exists in the distal cavity, making it difficult to assign unambiguously the residues involved in the stabilization of the bound ligand. To understand these experimental results with atomistic detail, we performed classical molecular dynamics simulations of the oxy, carboxy, and deoxy proteins. The free energy profiles for ligand migration suggest that there is a key residue, GlnE11, that presents an alternate conformation, in which a wide ligand migration tunnel is formed, consistently with the kinetic data. This tunnel is topologically related to the one found in group I trHbs. On the other hand, the results for the CO and O₂ bound protein show that GlnE11 is directly involved in the stabilization of the cordinated ligand, playing a similar role as TyrB10 and TrpG8 in other trHbs. Our results not only reconcile the structural data with the kinetic information, but also provide additional insight into the general behaviour of trHbs. Proteins 2010. © 2009 Wiley‐Liss, Inc.