Water and ligand entry in myoglobin: assessing the speed and extent of heme pocket hydration after CO photodissociation

Efficient conversion of hemoglobin to a non-vasoactive oxygen carrier by site-specific cross-linking with azido acyl methyl phosphates followed by bio-orthogonal CuAAC with a bis-alkyne

While cross-linked hemoglobin tetramers are functional acellular oxygen carriers, their ability to scavenge endogenous nitric oxide (NO) by endothelial pore penetration results in adverse cardiovascular effects. Animal studies established that cross-linked human hemoglobins, chemically joined into a double protein, avoid NO scavenging, presumably due to their larger size preventing penetration into endothelial regions that produce NO. In the present report, we utilize azide-containing acyl phosphate reagents to form cross-linked hemoglobins then bio-orthogonally click-couple them with a bis-alkyne (CuAAC). The production of these larger oxygen-carrying hemoglobin conjugates is obtained in high yields through subunit-specific cross-linking between each βLys82 ε-amino group. The methyl phosphate leaving groups provide electrostatically induced β-subunit site-selectivity, producing azido-cross-linked hemoglobin that undergoes highly efficient CuAAC compared with previous cross-linkers. The acyl phosphates also efficiently cross-link both T-state and R-state hemoglobin. The resulting bis- and tris-tetrameric hemoglobin conjugates exhibit oxygen affinity and cooperativity that are comparable to those of the native protein. The hemoglobin derivatives from the process we describe can function as sources of oxygen in biomedical applications, such as in ex-vivo donor organ perfusion.

Kinetic mechanisms for O2 binding to myoglobins and hemoglobins

Antonini and Brunori's 1971 book "Hemoglobin and Myoglobin in Their Reactions with Ligands" was a truly remarkable publication that summarized almost 100 years of research on O₂ binding to these globins. Over the ensuing 50 years, ultra-fast laser photolysis techniques, high-resolution and time resolved X-ray crystallography, molecular dynamics simulations, and libraries of recombinant myoglobin (Mb) and hemoglobin (Hb) variants have provided structural interpretations of O₂ binding to these proteins. The resultant mechanisms provide quantitative descriptions of the stereochemical factors that govern overall affinity, including proximal and distal steric restrictions that affect iron reactivity and favorable positive electrostatic interactions that preferentially stabilize bound O₂. The pathway for O₂ uptake and release by Mb and subunits of Hb has been mapped by screening libraries of site-directed mutants in laser photolysis experiments. O₂ enters mammalian Mb and the α and β subunits of human HbA through a channel created by upward and outward rotation of the distal His at the E7 helical position, is non-covalently captured in the interior of the distal cavity, and then internally forms a bond with the heme Fe(II) atom. O₂ dissociation is governed by disruption of hydrogen bonding interactions with His (E7), breakage of the Fe(II)-O₂ bond, and then competition between rebinding and escape through the E7-gate. The structural features that govern the rates of both the individual steps and overall reactions have been determined and provide the framework for: (1) defining the physiological functions of specific globins and their evolution; (2) understanding the clinical features of hemoglobinopathies; and (3) designing safer and more efficient acellular hemoglobin-based oxygen carriers (HBOCs) for transfusion therapy, organ preservation, and other commercially relevant O₂ transport and storage processes.

Distinctive structural properties of THB11, a pentacoordinate Chlamydomonas reinhardtii truncated hemoglobin with N- and C-terminal extensions

Hemoglobins (Hbs) utilize heme b as a cofactor and are found in all kingdoms of life. The current knowledge reveals an enormous variability of Hb primary sequences, resulting in topological, biochemical and physiological individuality. As Hbs appear to modulate their reactivities through specific combinations of structural features, predicting the characteristics of a given Hb is still hardly possible. The unicellular green alga Chlamydomonas reinhardtii contains 12 genes encoding diverse Hbs of the truncated lineage, several of which possess extended N- or C-termini of unknown function. Studies on some of the Chlamydomonas Hbs revealed yet unpredictable structural and biochemical variations, which, along with a different expression of their genes, suggest diverse physiological roles. Chlamydomonas thus represents a promising system to analyze the diversification of Hb structure, biochemistry and physiology. Here, we report the crystal structure, resolved to 1.75 Å, of the heme-binding domain of cyanomet THB11 (Cre16.g662750), one of the pentacoordinate algal Hbs, which offer a free Fe-coordination site in the reduced state. The overall fold of THB11 is conserved, but individual features such as a kink in helix E, a tilted heme plane and a clustering of methionine residues at a putative tunnel exit appear to be unique. Both N- and C-termini promote the formation of oligomer mixtures, and the absence of the C terminus results in reduced nitrite reduction rates. This work widens the structural and biochemical knowledge on the 2/2Hb family and suggests that the N- and C-terminal extensions of the Chlamydomonas 2/2Hbs modulate their reactivity by intermolecular interactions.

Do Osmolytes Impact the Structure and Dynamics of Myoglobin?

Osmolytes are small organic compounds that can affect the stability of proteins in living cells. The mechanism of osmolytes' protective effects on protein structure and dynamics has not been fully explained, but in general, two possibilities have been suggested and examined: a direct interaction of osmolytes with proteins (water replacement hypothesis), and an indirect interaction (vitrification hypothesis). Here, to investigate these two possible mechanisms, we studied myoglobin-osmolyte systems using FTIR, UV-vis, CD, and femtosecond IR pump-probe spectroscopy. Interestingly, noticeable changes are observed in both the lifetime of the CO stretch of CO-bound myoglobin and the spectra of UV-vis, CD, and FTIR upon addition of the osmolytes. In addition, the temperature-dependent CD studies reveal that the protein's thermal stability depends on molecular structure, hydrogen-bonding ability, and size of osmolytes. We anticipate that the present experimental results provide important clues about the complicated and intricate mechanism of osmolyte effects on protein structure and dynamics in a crowded cellular environment.

Cold-Adaptation Signatures in the Ligand Rebinding Kinetics to the Truncated Hemoglobin of the Antarctic Bacterium Pseudoalteromonas haloplanktis TAC125

Cold-adapted organisms have evolved proteins endowed with higher flexibility and lower stability in comparison to their thermophilic homologues, resulting in enhanced reaction rates at low temperatures. In this context, protein-bound water molecules were suggested to play a major role, and their weaker interactions at protein active sites have been associated with cold adaptation. In this work, we tested this hypothesis on truncated hemoglobins (a family of microbial heme-proteins of yet-unclear function) applying molecular dynamics simulations and ligand-rebinding kinetics on a protein from the Antarctic bacterium Pseudoalteromonas haloplanktis TAC125 in comparison with its thermophilic Thermobifida fusca homologue. The CO rebinding kinetics of the former highlight several geminate phases, with an unusually long-lived geminate intermediate. An articulated tunnel with at least two distinct docking sites was identified by analysis of molecular dynamics simulations and was suggested to be at the origin of the unusual geminate rebinding phase. Water molecules are present in the distal pocket, but their stabilization by TrpG8, TyrB10, and HisCD1 is much weaker than in thermophilic Thermobifida fusca truncated hemoglobin, resulting in a faster geminate rebinding. Our results support the hypothesis that weaker water-molecule interactions at the reaction site are associated with cold adaptation.

Structure and function of haemoglobins

Haemoglobin (Hb) is widely known as the iron-containing protein in blood that is essential for O₂ transport in mammals. Less widely recognised is that erythrocyte Hb belongs to a large family of Hb proteins with members distributed across all three domains of life-bacteria, archaea and eukaryotes. This review, aimed chiefly at researchers new to the field, attempts a broad overview of the diversity, and common features, in Hb structure and function. Topics include structural and functional classification of Hbs; principles of O₂ binding affinity and selectivity between O₂/NO/CO and other small ligands; hexacoordinate (containing bis-imidazole coordinated haem) Hbs; bacterial truncated Hbs; flavohaemoglobins; enzymatic reactions of Hbs with bioactive gases, particularly NO, and protection from nitrosative stress; and, sensor Hbs. A final section sketches the evolution of work on the structural basis for allosteric O₂ binding by mammalian RBC Hb, including the development of newer kinetic models. Where possible, reference to historical works is included, in order to provide context for current advances in Hb research.

Role of Heme Pocket Water in Allosteric Regulation of Ligand Reactivity in Human Hemoglobin

Water molecules can enter the heme pockets of unliganded myoglobins and hemoglobins, hydrogen bond with the distal histidine, and introduce steric barriers to ligand binding. The spectrokinetics of photodissociated CO complexes of human hemoglobin and its isolated α and β chains were analyzed for the effect of heme hydration on ligand rebinding. A strong coupling was observed between heme hydration and quaternary state. This coupling may contribute significantly to the 20-60-fold difference between the R- and T-state bimolecular CO binding rate constants and thus to the modulation of ligand reactivity that is the hallmark of hemoglobin allostery. Heme hydration proceeded over the course of several kinetic phases in the tetramer, including the R to T quaternary transition. An initial 150 ns hydration phase increased the R-state distal pocket water occupancy, nw(R), to a level similar to that of the isolated α (∼60%) and β (∼10%) chains, resulting in a modest barrier to ligand binding. A subsequent phase, concurrent with the first step of the R → T transition, further increased the level of heme hydration, increasing the barrier. The final phase, concurrent with the final step of the allosteric transition, brought the water occupancy of the T-state tetramer, nw(T), even higher and close to full occupancy in both the α and β subunits (∼90%). This hydration level could present an even larger barrier to ligand binding and contribute significantly to the lower iron reactivity of the T state toward CO.

A quantitative model for oxygen uptake and release in a family of hemeproteins

MOTIVATION

Hemeproteins have many diverse functions that largely depend on the rate at which they uptake or release small ligands, like oxygen. These proteins have been extensively studied using either simulations or experiments, albeit only qualitatively and one or two proteins at a time.

RESULTS

We present a physical-chemical model, which uses data obtained exclusively from computer simulations, to describe the uptake and release of oxygen in a family of hemeproteins, called truncated hemoglobins (trHbs). Through a rigorous statistical analysis we demonstrate that our model successfully recaptures all the reported experimental oxygen association and dissociation kinetic rate constants, thus allowing us to establish the key factors that determine the rates at which these hemeproteins uptake and release oxygen. We found that internal tunnels as well as the distal site water molecules control ligand uptake, whereas oxygen stabilization by distal site residues controls ligand release. Because these rates largely determine the functions of these hemeproteins, these approaches will also be important tools in characterizing the trHbs members with unknown functions.

CONTACT

lboechi@ic.fcen.uba.ar

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Evolutionary and Functional Relationships in the Truncated Hemoglobin Family

Predicting function from sequence is an important goal in current biological research, and although, broad functional assignment is possible when a protein is assigned to a family, predicting functional specificity with accuracy is not straightforward. If function is provided by key structural properties and the relevant properties can be computed using the sequence as the starting point, it should in principle be possible to predict function in detail. The truncated hemoglobin family presents an interesting benchmark study due to their ubiquity, sequence diversity in the context of a conserved fold and the number of characterized members. Their functions are tightly related to O2 affinity and reactivity, as determined by the association and dissociation rate constants, both of which can be predicted and analyzed using in-silico based tools. In the present work we have applied a strategy, which combines homology modeling with molecular based energy calculations, to predict and analyze function of all known truncated hemoglobins in an evolutionary context. Our results show that truncated hemoglobins present conserved family features, but that its structure is flexible enough to allow the switch from high to low affinity in a few evolutionary steps. Most proteins display moderate to high oxygen affinities and multiple ligand migration paths, which, besides some minor trends, show heterogeneous distributions throughout the phylogenetic tree, again suggesting fast functional adaptation. Our data not only deepens our comprehension of the structural basis governing ligand affinity, but they also highlight some interesting functional evolutionary trends.

Molecular dynamics simulations indicate that tyrosineB10 limits motions of distal histidine to regulate CO binding in soybean leghemoglobin

Myoglobin (Mb) uses strong electrostatic interaction in its distal heme pocket to regulate ligand binding. The mechanism of regulation of ligand binding in soybean leghemoglobin a (Lba) has been enigmatic and more so due to the absence of gaseous ligand bound atomic resolution three-dimensional structure of the plant globin. While the 20-fold higher oxygen affinity of Lba compared with Mb is required for its dual physiological function, the mechanism by which this high affinity is achieved is only emerging. Extensive mutational analysis combined with kinetic and CO-FT-IR spectroscopic investigation led to the hypothesis that Lba depended on weakened electrostatic interaction between distal HisE7 and bound ligand achieved by invoking B10Tyr, which itself hydrogen bonds with HisE7 thus restricting it in a single conformation detrimental to Mb-like strong electrostatic interaction. Such theory has been re-assessed here using CO-Lba in silico model and molecular dynamics simulation. The investigation supports the presence of at least two major conformations of HisE7 in Lba brought about by imidazole ring flip, one of which makes hydrogen bonds effectively with B10Tyr affecting the former's ability to stabilize bound ligand, while the other does not. However, HisE7 in Lba has limited conformational freedom unlike high frequency of imidazole ring flips observed in Mb and in TyrB10Leu mutant of Lba. Thus, it appears that TyrB10 limits the conformational freedom of distal His in Lba, tuning down ligand dissociation rate constant by reducing the strength of hydrogen bonding to bound ligand, which the freedom of distal His of Mb allows.

O2 and Water Migration Pathways between the Solvent and Heme Pockets of Hemoglobin with Open and Closed Conformations of the Distal HisE7

Hemoglobin transports O2 by binding the gas at its four hemes. Hydrogen bonding between the distal histidine (HisE7) and heme-bound O2 significantly increases the affinity of human hemoglobin (HbA) for this ligand. HisE7 is also proposed to regulate the release of O2 to the solvent via a transient E7 channel. To reveal the O2 escape routes controlled by HisE7 and to evaluate its role in gating heme access, we compare simulations of O2 diffusion from the distal heme pockets of the T and R states of HbA performed with HisE7 in its open (protonated) and closed (neutral) conformations. Irrespective of HisE7's conformation, we observe the same four or five escape routes leading directly from the α- or β-distal heme pockets to the solvent. Only 21-53% of O2 escapes occur via these routes, with the remainder escaping through routes that encompass multiple internal cavities in HbA. The conformation of the distal HisE7 controls the escape of O2 from the heme by altering the distal pocket architecture in a pH-dependent manner, not by gating the E7 channel. Removal of the HisE7 side chain in the GlyE7 variant exposes the distal pockets to the solvent, and the percentage of O2 escapes to the solvent directly from the α- or β-distal pockets of the mutant increases to 70-88%. In contrast to O2, the dominant water route from the bulk solvent is gated by HisE7 because protonation and opening of this residue dramatically increase the rate of influx of water into the empty distal heme pockets. The occupancy of the distal heme site by a water molecule, which functions as an additional nonprotein barrier to binding of the ligand to the heme, is also controlled by HisE7. Overall, analysis of gas and water diffusion routes in the subunits of HbA and its GlyE7 variant sheds light on the contribution of distal HisE7 in controlling polar and nonpolar ligand movement between the solvent and the hemes.

Identification of Mutational Hot Spots for Substrate Diffusion: Application to Myoglobin

The pathways by which small molecules (substrates or inhibitors) access active sites are a key aspect of the function of enzymes and other proteins. A key problem in designing or altering such proteins is to identify sites for mutation that will have the desired effect on the substrate transport properties. While specific access channels have been invoked in the past, molecular simulations suggest that multiple routes are possible, complicating the analysis. This complexity, however, can be captured by a Markov State Model (MSM) of the ligand diffusion process. We have developed a sensitivity analysis of the resulting rate matrix, which identifies the locations where mutations should have the largest effect on the diffusive on rate. We apply this method to myoglobin, which is the best characterized example both from experiment and simulation. We validate the approach by translating the sensitivity parameter obtained from this method into the CO binding rates in myoglobin upon mutation, resulting in a semi-quantitative correlation with experiments. The model is further validated against an explicit simulation for one of the experimental mutants.

Full kinetics of CO entry, internal diffusion, and exit in myoglobin from transition-path theory simulations

We use Markovian milestoning molecular dynamics (MD) simulations on a tessellation of the collective variable space for CO localization in myoglobin to estimate the kinetics of entry, exit, and internal site-hopping. The tessellation is determined by analysis of the free-energy surface in that space using transition-path theory (TPT), which provides criteria for defining optimal milestones, allowing short, independent, cell-constrained MD simulations to provide properly weighted kinetic data. We coarse grain the resulting kinetic model at two levels: first, using crystallographically relevant internal cavities and their predicted interconnections and solvent portals; and second, as a three-state side-path scheme inspired by similar models developed from geminate recombination experiments. We show semiquantitative agreement with experiment on entry and exit rates and in the identification of the so-called "histidine gate" at position 64 through which ≈90% of flux between solvent and the distal pocket passes. We also show with six-dimensional calculations that the minimum free-energy pathway of escape through the histidine gate is a "knock-on" mechanism in which motion of the ligand and the gate are sequential and interdependent. In total, these results suggest that such TPT simulations are indeed a promising approach to overcome the practical time-scale limitations of MD to allow reliable estimation of transition mechanisms and rates among metastable states.

Solvation and cavity occupation in biomolecules

BACKGROUND

Solvation density locations are important for protein dynamics and structure. Knowledge of the preferred hydration sites at biomolecular interfaces and those in the interior of cavities can enhance understanding of structure and function. While advanced X-ray diffraction methods can provide accurate atomic structures for proteins, that technique is challenged when it comes to providing accurate hydration structures, especially for interfacial and cavity bound solvent molecules.

METHODS

Advances in integral equation theories which include more accurate methods for calculating the long-ranged Coulomb interaction contributions to the three-dimensional distribution functions make it possible to calculate angle dependent average solvent structure, accurately, around and inside irregular molecular conformations. The proximal radial distribution method provides another approximate method to determine average solvent structures for biomolecular systems based on a proximal or near neighbor solvent distribution that can be constructed from previously collected solvent distributions. These two approximate methods, along with all-atom molecular dynamics simulations are used to determine the solvent density inside the myoglobin heme cavity.

DISCUSSION AND RESULTS

Myoglobin is a good test system for these methods because the cavities are many and one is large, tens of Å(3), but is shown to have only four hydration sites. These sites are not near neighbors which implies that the large cavity must have more than one way in and out.

CONCLUSIONS

Our results show that main solvation sites are well reproduced by all three methods. The techniques also produce a clearly identifiable solvent pathway into the interior of the protein.

GENERAL SIGNIFICANCE

The agreement between molecular dynamics and less computationally demanding approximate methods is encouraging. This article is part of a Special Issue entitled Recent developments of molecular dynamics.

Ligand uptake modulation by internal water molecules and hydrophobic cavities in hemoglobins

Internal water molecules play an active role in ligand uptake regulation, since displacement of retained water molecules from protein surfaces or cavities by incoming ligands can promote favorable or disfavorable effects over the global binding process. Detection of these water molecules by X-ray crystallography is difficult given their positional disorder and low occupancy. In this work, we employ a combination of molecular dynamics simulations and ligand rebinding over a broad time range to shed light into the role of water molecules in ligand migration and binding. Computational studies on the unliganded structure of the thermostable truncated hemoglobin from Thermobifida fusca (Tf-trHbO) show that a water molecule is in the vicinity of the iron heme, stabilized by WG8 with the assistance of YCD1, exerting a steric hindrance for binding of an exogenous ligand. Mutation of WG8 to F results in a significantly lower stabilization of this water molecule and in subtle dynamical structural changes that favor ligand binding, as observed experimentally. Water is absent from the fully hydrophobic distal cavity of the triple mutant YB10F-YCD1F-WG8F (3F), due to the lack of residues capable of stabilizing it nearby the heme. In agreement with these effects on the barriers for ligand rebinding, over 97% of the photodissociated ligands are rebound within a few nanoseconds in the 3F mutant case. Our results demonstrate the specific involvement of water molecules in shaping the energetic barriers for ligand migration and binding.

Long-range electron transfer with myoglobin immobilized at Au/mixed-SAM junctions: mechanistic impact of the strong protein confinement

Horse muscle myoglobin (Mb) was tightly immobilized at Au-deposited ~15-Å-thick mixed-type (1:1) alkanethiol SAMs, HS-(CH₂)₁₁-COOH/HS-(CH₂)₁₁-OH, and placed in contact with buffered H₂O or D₂O solutions. Fast-scan cyclic voltammetry (CV) and a Marcus-equation-based analysis were applied to determine unimolecular standard rate constants and reorganization free energies for electron transfer (ET), under variable-temperature (15-55 °C) and -pressure (0.01-150 MPa) conditions. The CV signal was surprisingly stable and reproducible even after multiple temperature and pressure cycles. The data analysis revealed the following values: standard rate constant, 33 s⁻¹ (25 °C, 0.01 MPa, H₂O); reorganization free energy, 0.5 ± 0.1 eV (throughout); activation enthalpy, 12 ± 3 kJ mol⁻¹; activation volume, -3.1 ± 0.2 cm³ mol⁻¹; and pH-dependent solvent kinetic isotope effect (k(H)⁰/k(D)⁰), 0.7-1.4. Furthermore, the values for the rate constant and reorganization free energy are very similar to those previously found for cytochrome c electrostatically immobilized at the monocomponent Au/HS-(CH₂)₁₁-COOH junction. In vivo, Mb apparently forms a natural electrostatic complex with cytochrome b₅ (cyt-b₅) through the "dynamic" (loose) docking pattern, allowing for a slow ET that is intrinsically coupled to the water's removal from the "defective" heme iron (altogether shaping the biological repair mechanism for Mb's "met" form). In contrary, our experiments rather mimic the case of a "simple" (tight) docking of the redesigned (mutant) Mb with cyt-b₅ (Nocek et al. J. Am. Chem. Soc. 2010, 132, 6165-6175). According to our analysis, in this configuration, Mb's distal pocket (linked to the "ligand channel") seems to be arrested within the restricted configuration, allowing the rate-determining reversible ET process to be coupled only to the inner-sphere reorganization (minimal elongation/shortening of an Fe-OH₂ bond) rather than the pronounced detachment (rebinding) of water and, hence, to be much faster.

Internal water and microsecond dynamics in myoglobin

Myoglobin (Mb) binds diatomic ligands, like O₂, CO, and NO, in a cavity that is only transiently accessible. Crystallography and molecular simulations show that the ligands can migrate through an extensive network of transiently connected cavities but disagree on the locations and occupancy of internal hydration sites. Here, we use water ²H and ¹⁷O magnetic relaxation dispersion (MRD) to characterize the internal water molecules in Mb under physiological conditions. We find that equine carbonmonoxy Mb contains 4.5 ± 1.0 ordered internal water molecules with a mean survival time of 5.6 ± 0.5 μs at 25 °C. The likely locations of these water molecules are the four polar hydration sites, including one of the xenon-binding cavities, that are fully occupied in all high-resolution crystal structures of equine Mb. The finding that water escapes from these sites, located 17–31 Å apart in the protein, on the same μs time scale suggests a global exchange mechanism. We propose that this mechanism involves transient penetration of the protein by H-bonded water chains. Such a mechanism could play a functional role by eliminating trapped ligands. In addition, the MRD results indicate that 2 or 3 of the 11 histidine residues of equine Mb undergo intramolecular hydrogen exchange on a μs time scale.

Functional consequences of the open distal pocket of dehaloperoxidase-hemoglobin observed by time-resolved X-ray crystallography

Using time-resolved X-ray crystallography, we contrast a bifunctional dehaloperoxidase-hemoglobin (DHP) with previously studied examples of myoglobin and hemoglobin to understand the functional role of the distal pocket of globins. One key functional difference between DHP and other globins is the requirement that H2O2 enter the distal pocket of oxyferrous DHP to displace O2 from the heme Fe atom and thereby activate the heme for the peroxidase function. The open architecture of DHP permits more than one molecule to simultaneously enter the distal pocket of the protein above the heme to facilitate the unique peroxidase cycle starting from the oxyferrous state. The time-resolved X-ray data show that the distal pocket of DHP lacks a protein valve found in the two other globins that have been studied previously. The photolyzed CO ligand trajectory in DHP does not have a docking site; rather, the CO moves immediately to the Xe-binding site. From there, CO can escape but can also recombine an order of magnitude more rapidly than in other globins. The contrast with DHP dynamics and function more precisely defines the functional role of the multiple conformational states of myoglobin. Taken together with the high reduction potential of DHP, the open distal site helps to explain how a globin can also function as a peroxidase.

Reactivity and dynamics of H2S, NO, and O2 interacting with hemoglobins from Lucina pectinata

Hemoglobin HbI from the clam Lucina pectinata is involved in H2S transport, whereas homologous heme protein HbII/III is involved in O2 metabolism. Despite similar tertiary structures, HbI and HbII/III exhibit very different reactivity toward heme ligands H2S, O2, and NO. To investigate this reactivity at the heme level, we measured the dynamics of ligand interaction by time-resolved absorption spectroscopy in the picosecond to nanosecond time range. We demonstrated that H2S can be photodissociated from both ferric and ferrous HbI. H2S geminately rebinds to ferric and ferrous out-of-plane iron with time constants (τgem) of 12 and 165 ps, respectively, with very different proportions of photodissociated H2S exiting the protein (24% in ferric and 80% in ferrous HbI). The Gln(E7)His mutation considerably changes H2S dynamics in ferric HbI, indicating the role of Gln(E7) in controling H2S reactivity. In ferric HbI, the rate of diffusion of H2S from the solvent into the heme pocket (kentry) is 0.30 μM(-1) s(-1). For the HbII/III-O2 complex, we observed mainly a six-coordinate vibrationally excited heme-O2 complex with O2 still bound to the iron. This explains the low yield of O2 photodissociation and low koff from HbII/III, compared with those of HbI and Mb. Both isoforms behave very differently with regard to NO and O2 dynamics. Whereas the amplitude of geminate rebinding of O2 to HbI (38.5%) is similar to that of myoglobin (34.5%) in spite of different distal heme sites, it appears to be much larger for HbII/III (77%). The distal Tyr(B10) side chain present in HbII/III increases the energy barrier for ligand escape and participates in the stabilization of bound O2 and NO.

Ligand migration through hemeprotein cavities: insights from laser flash photolysis and molecular dynamics simulations

The presence of cavities and tunnels in the interior of proteins, in conjunction with the structural plasticity arising from the coupling to the thermal fluctuations of the protein scaffold, has profound consequences on the pathways followed by ligands moving through the protein matrix. In this perspective we discuss how quantitative analysis of experimental rebinding kinetics from laser flash photolysis, trapping of unstable conformational states by embedding proteins within the nanopores of silica gels, and molecular simulations can synergistically converge to gain insight into the migration mechanism of ligands. We show how the evaluation of the free energy landscape for ligand diffusion based on the outcome of computational techniques can assist the definition of sound reaction schemes, leading to a comprehensive understanding of the broad range of chemical events and time scales that encompass the transport of small ligands in hemeproteins.

Investigations of ferric heme cyanide photodissociation in myoglobin and horseradish peroxidase

The photodissociation of cyanide from ferric myoglobin (MbCN) and horseradish peroxidase (HRPCN) has definitively been observed. This has implications for the interpretation of ultrafast IR (Helbing et al. Biophys. J. 2004, 87, 1881-1891) and optical (Gruia et al. Biophys. J. 2008, 94, 2252-2268) studies that had previously suggested the Fe-CN bond was photostable in MbCN. The photolysis of ferric MbCN takes place with a quantum yield of ~75%, and the resonance Raman spectrum of the photoproduct observed in steady-state experiments as a function of laser power and sample spinning rate is identical to that of ferric Mb (metMb). The data are quantitatively analyzed using a simple model where cyanide is photodissociated and, although geminate rebinding with a rate of kBA ≈ (3.6 ps)(-1) is the dominant process, some CN(-) exits from the distal heme pocket and is replaced by water. Using independently determined values for the CN(-) association rate, we find that the CN(-) escape rate from the ferric myoglobin pocket to the solution at 293 K is kout ≈ (1-2) × 10(7) s(-1). This value is very similar to, but slightly larger than, the histidine gated escape rate of CO from Mb (1.1 × 10(7) s(-1)) under the same conditions. The analysis leads to an escape probability kout/(kout + kBA) ~ 10(-4), which is unobservable in most time domain kinetic measurements. However, the photolysis is surprisingly easy to detect in Mb using cw resonance Raman measurements. This is due to the anomalously slow CN(-) bimolecular association rate (170 M(-1) s(-1)), which arises from the need for water to exchange at the ferric heme binding site of Mb. In contrast, ferric HRP does not have a heme bound water molecule and its CN(-) bimolecular association rate is larger by ~10(3), making the CN(-) photolysis more difficult to observe.

Small ligand-globin interactions: reviewing lessons derived from computer simulation

In this work we review the application of classical and quantum-mechanical atomistic computer simulation tools to the investigation of small ligand interaction with globins. In the first part, studies of ligand migration, with its connection to kinetic association rate constants (kon), are presented. In the second part, we review studies for a variety of ligands such as O2, NO, CO, HS(-), F(-), and NO2(-) showing how the heme structure, proximal effects, and the interactions with the distal amino acids can modulate protein ligand binding. The review presents mainly results derived from our previous works on the subject, in the context of other theoretical and experimental studies performed by others. The variety and extent of the presented data yield a clear example of how computer simulation tools have, in the last decade, contributed to our deeper understanding of small ligand interactions with globins. This article is part of a Special Issue entitled: Oxygen Binding and Sensing Proteins.

A computational study of water and CO migration sites and channels inside myoglobin

Pathways are computed for transport of H₂O and CO in myoglobin (Mb), using the single sweep and zero-temperature string methods in a fully atomistic, explicitly solvated model system. Our predictions of sites and barriers in the pathways for CO transport agree with previous studies. For H₂O, we predict a binding site in the distal pocket (DP), in agreement with crystallographic observations, and another one close to Leu 29 which explains the importance of this residue in controlling the pocket's hydrophobicity, as well as disordered minima in the largely apolar xenon cavities. In particular, H₂O can occupy and transition among the xenon cavities, Xe4, Xe2, and Xe3. Our results support the hypothesis that the thermodynamically most favorable entry/exit portal for H₂O is the so-called histidine gate (HG), the same as for CO. This result, along with the observation of water occupation of both DP and apolar Xe cavities, suggest that water and small gas molecules like CO compete for access to the protein interior, and therefore models of gas molecule transport within proteins should also explicitly consider water transport.

Following ligand migration pathways from picoseconds to milliseconds in type II truncated hemoglobin from Thermobifida fusca

CO recombination kinetics has been investigated in the type II truncated hemoglobin from Thermobifida fusca (Tf-trHb) over more than 10 time decades (from 1 ps to ∼100 ms) by combining femtosecond transient absorption, nanosecond laser flash photolysis and optoacoustic spectroscopy. Photolysis is followed by a rapid geminate recombination with a time constant of ∼2 ns representing almost 60% of the overall reaction. An additional, small amplitude geminate recombination was identified at ∼100 ns. Finally, CO pressure dependent measurements brought out the presence of two transient species in the second order rebinding phase, with time constants ranging from ∼3 to ∼100 ms. The available experimental evidence suggests that the two transients are due to the presence of two conformations which do not interconvert within the time frame of the experiment. Computational studies revealed that the plasticity of protein structure is able to define a branched pathway connecting the ligand binding site and the solvent. This allowed to build a kinetic model capable of describing the complete time course of the CO rebinding kinetics to Tf-trHb.

Molecular probes of the mechanism of cytochrome P450. Oxygen traps a substrate radical intermediate

The diagnostic substrate tetramethylcyclopropane (TMCP) has been reexamined as a substrate with three drug- and xenobiotic-metabolizing cytochrome P450 enzymes, human CYP2E1, CYP3A4 and rat CYP2B1. The major hydroxylation product in all cases was the unrearranged primary alcohol along with smaller amounts of a rearranged tertiary alcohol. Significantly, another ring-opened product, diacetone alcohol, was also observed. With CYP2E1 this product accounted for 20% of the total turnover. Diacetone alcohol also was detected as a product from TMCP with a biomimetic model catalyst, FeTMPyP, but not with a ruthenium porphyrin catalyst. Lifetimes of the intermediate radicals were determined from the ratios of rearranged and unrearranged products to be 120, 13 and 1ps for CYP2E1, CYP3A4 and CYP2B1, respectively, corresponding to rebound rates of 0.9×10(10)s(-1), 7.2×10(10)s(-1) and 1.0×10(12)s(-1). For the model iron porphyrin, FeTMPyP, a radical lifetime of 81ps and a rebound rate of 1.2×10(10)s(-1) were determined. These apparent radical lifetimes are consistent with earlier reports with a variety of CYP enzymes and radical clock substrates, however, the large amounts of diacetone alcohol with CYP2E1 and the iron porphyrin suggest that for these systems a considerable amount of the intermediate carbon radical is trapped by molecular oxygen. These results add to the view that cage escape of the intermediate carbon radical in [Fe(IV)-OH ()R] can compete with cage collapse to form a C-O bond. The results could be significant with regard to our understanding of iron-catalyzed C-H hydroxylation, the observation of P450-dependent peroxidation and the development of oxidative stress, especially for CYP2E1.

Time resolved thermodynamics associated with ligand photorelease in heme peroxidases and globins: Open access channels versus gated ligand release

Heme proteins represent a diverse class of biomolecules responsible for an extremely diverse array of physiological functions including electron transport, monooxygenation, ligand transport and storage, cellular signaling, respiration, etc. An intriguing aspect of these proteins is that such functional diversity is accomplished using a single type of heme macrocycle based upon iron protoporphyrin IX. The functional diversity originates from a delicate balance of inter-molecular interactions within the protein matrix together with well choreographed dynamics that modulate the heme electronic structure as well as ligand entry/exit pathways from the bulk solvent to the active site. Of particular interest are the dynamics and energetics associated with the entry/exit of ligands as this process plays a significant role in regulating the rates of heme protein activity. Time-resolved photoacoustic calorimetry (PAC) has emerged as a powerful tool through which to probe the underlying energetics associated with small molecule dissociation and release to the bulk solvent in heme proteins on time scales from tens of nanoseconds to several microseconds. In this review, the results of PAC studies on various classes of heme proteins are summarized highlighting how different protein structures affect the thermodynamics of ligand migration. This article is part of a Special Issue entitled: Protein Dynamics: Experimental and Computational Approaches.

Blocking the gate to ligand entry in human hemoglobin

His(E7) to Trp replacements in HbA lead to markedly biphasic bimolecular CO rebinding after laser photolysis. For isolated mutant subunits, the fraction of fast phase increases with increasing [CO], suggesting a competition between binding to an open conformation with an empty E7 channel and relaxation to blocked or closed, slowly reacting states. The rate of conformational relaxation of the open state is ∼18,000 s(-1) in α subunits and ∼10-fold faster in β subunits, ∼175,000 s(-1). Crystal structures were determined for tetrameric α(WT)β(Trp-63) HbCO, α(Trp-58)β(WT) deoxyHb, and Trp-64 deoxy- and CO-Mb as controls. In Trp-63(E7) βCO, the indole side chain is located in the solvent interface, blocking entry into the E7 channel. Similar blocked Trp-64(E7) conformations are observed in the mutant Mb crystal structures. In Trp-58(E7) deoxy-α subunits, the indole side chain fills both the channel and the distal pocket, forming a completely closed state. The bimolecular rate constant for CO binding, k'(CO), to the open conformations of both mutant Hb subunits is ∼80-90 μm(-1) s(-1), whereas k'(CO) for the completely closed states is 1000-fold slower, ∼0.08 μm(-1) s(-1). A transient intermediate with k'(CO) ≈ 0.7 μm(-1) s(-1) is observed after photolysis of Trp-63(E7) βCO subunits and indicates that the indole ring blocks the entrance to the E7 channel, as observed in the crystal structures of Trp(E7) deoxyMb and βCO subunits. Thus, either blocking or completely filling the E7 channel dramatically slows bimolecular binding, providing strong evidence that the E7 channel is the major pathway (≥90%) for ligand entry in human hemoglobin.

Mechanisms of peroxynitrite interactions with heme proteins

Oxygenated heme proteins are known to react rapidly with nitric oxide (NO) to produce peroxynitrite (PN) at the heme site. This process could lead either to attenuation of the effects of NO or to nitrosative protein damage. PN is a powerful nitrating and oxidizing agent that has been implicated in a variety of cell injuries. Accordingly, it is important to delineate the nature and variety of reaction mechanisms of PN interactions with heme proteins. In this Forum, we survey the range of reactions of PN with heme proteins, with particular attention to myoglobin and cytochrome c. While these two proteins are textbook paradigms for oxygen binding and electron transfer, respectively, both have recently been shown to have other important functions that involve NO and PN. We have recently described direct evidence that ferrylmyolgobin (ferrylMb) and nitrogen dioxide (NO(2)) are both produced during the reaction of PN and metmyolgobin (metMb) (Su, J.; Groves, J. T. J. Am. Chem. Soc. 2009, 131, 12979-12988). Kinetic evidence indicates that these products evolve from the initial formation of a caged radical intermediate [Fe(IV) horizontal lineO.NO(2)]. This caged pair reacts mainly via internal return with a rate constant k(r) to form metMb and nitrate in an oxygen-rebound scenario. Detectable amounts of ferrylMb are observed by stopped-flow spectrophotometry, appearing at a rate consistent with the rate, k(obs), of heme-mediated PN decomposition. Freely diffusing NO(2), which is liberated concomitantly from the radical pair (k(e)), preferentially nitrates myoglobin Tyr103 and added fluorescein. For cytochrome c, Raman spectroscopy has revealed that a substantial fraction of cytochrome c converts to a beta-sheet structure, at the expense of turns and helices at low pH (Balakrishnan, G.; Hu, Y.; Oyerinde, O. F.; Su, J.; Groves, J. T.; Spiro, T. G. J. Am. Chem. Soc., 2007, 129, 504-505). It is proposed that a short beta-sheet segment, comprising residues 37-39 and 58-61, extends itself into the large 37-61 loop when the latter is destabilized by protonation of H26, which forms an anchoring hydrogen bond to loop residue P44. This conformation change ruptures the Met80-Fe bond, as revealed by changes in ligation-sensitive Raman bands. It also induces peroxidase activity with the same temperature profile. This process is suggested to model the apoptotic peroxidation of cardiolipin by cytochrome c.

Kinetic spectroscopy of heme hydration and ligand binding in myoglobin and isolated hemoglobin chains: an optical window into heme pocket water dynamics

The entry of a water molecule into the distal heme pocket of pentacoordinate heme proteins such as myoglobin and the alpha,beta chains of hemoglobin can be detected by time-resolved spectroscopy in the heme visible bands after photolysis of the CO complex. Reviewing the evidence from spectrokinetic studies of Mb variants, we find that this optical method measures the occupancy of non(heme)coordinated water in the distal pocket, n(w), with high fidelity. This evidence further suggests that perturbation of the kinetic barrier presented by distal pocket water is often the dominant mechanism by which active site mutations affect the bimolecular rate constant for CO binding. Water entry into the heme pockets of isolated hemoglobin subunits was detected by optical methods. Internal hydration is higher in the native alpha chains than in the beta chains, in agreement with previous crystallographic results for the subunits within Hb tetramers. The kinetic parameters obtained from modeling of the water entry and ligand rebinding in Mb mutants and native Hb chains are consistent with an inverse dependence of the bimolecular association rate constant on the water occupancy factor. This correlation suggests that water and ligand mutually exclude one another from the distal pockets of both types of hemoglobin chains and myoglobin.

Ultrafast dynamics of diatomic ligand binding to nitrophorin 4

Nitrophorin 4 (NP4) is a heme protein that stores and delivers nitric oxide (NO) through pH-sensitive conformational change. This protein uses the ferric state of a highly ruffled heme to bind NO tightly at low pH and release it at high pH. In this work, the rebinding kinetics of NO and CO to NP4 are investigated as a function of iron oxidation state and the acidity of the environment. The geminate recombination process of NO to ferrous NP4 at both pH 5 and pH 7 is dominated by a single approximately 7 ps kinetic phase that we attribute to the rebinding of NO directly from the distal pocket. The lack of pH dependence explains in part why NP4 cannot use the ferrous state to fulfill its function. The kinetic response of ferric NP4NO shows two distinct phases. The relative geminate amplitude of the slower phase increases dramatically as the pH is raised from 5 to 8. We assign the fast phase of NO rebinding to a conformation of the ferric protein with a closed hydrophobic pocket. The slow phase is assigned to the protein in an open conformation with a more hydrophilic heme pocket environment. Analysis of the ultrafast kinetics finds the equilibrium off-rate of NO to be proportional to the open state population as well as the pH-dependent amplitude of escape from the open pocket. When both factors are considered, the off-rate increases by more than an order of magnitude as the pH changes from 5 to 8. The recombination of CO to ferrous NP4 is observed to have a large nonexponential geminate amplitude with rebinding time scales of approximately 10(-11)-10(-9) s at pH 5 and approximately 10(-10)-10(-8) s at pH 7. The nonexponential CO rebinding kinetics at both pH 5 and pH 7 are accounted for using a simple model that has proven effective for understanding CO binding in a variety of other heme systems (Ye, X.; et al. Proc. Natl. Acad. Sci. U.S.A. 2007, 104, 14682).

Direct detection of the oxygen rebound intermediates, ferryl Mb and NO2, in the reaction of metmyoglobin with peroxynitrite

Oxygenated hemoproteins are known to react rapidly with nitric oxide (NO) to produce peroxynitrite (PN) at the heme site. This process could lead either to attenuation of the effects of NO or to nitrosative protein damage. Peroxynitrite is a powerful nitrating and oxidizing agent that has been implicated in a variety of cell injuries. Accordingly, it is important to delineate the nature and variety of reaction mechanisms of PN reactions with heme proteins. Here, we present direct evidence that ferrylMb and NO(2) are both produced during the reaction of PN and metmyoglobin (metMb). Kinetic evidence indicates that these products evolve from initial formation of a caged radical intermediate [Fe(IV)=O *NO(2)]. This caged pair reacts mainly via internal return with a rate constant k(r) to form metMb and nitrate in an oxygen rebound scenario. Detectable amounts of ferrylMb are observed by stopped-flow spectrophotometry, appearing at a rate consistent with the rate, k(obs), of heme-mediated PN decomposition. Freely diffusing NO(2), which is liberated concomitantly from the radical pair (k(e)), preferentially nitrates Tyr103 in horse heart myoglobin. The ratio of the rates of in-cage rebound and cage escape, k(r)/k(e), was found to be approximately 10 by examining the nitration yields of fluorescein, an external NO(2) trap. This rebound/escape model for the metMb/PN interaction is analogous to the behavior of alkyl hyponitrites and the well-studied geminate recombination processes of deoxymyoglobin with O(2), CO, and NO. The scenario is also similar to the stepwise events of substrate hydroxylation by cytochrome P450 and other oxygenases. It is likely, therefore, that the reaction of metMb with ONOO(-) and that of oxyMb with NO proceed through the same [Fe(IV)=O *NO(2)] caged radical intermediate and lead to similar outcomes. The results indicate that while oxyMb may reduce the concentration of intracellular NO, it would not eliminate the formation of NO(2) as a decomposition product of peroxynitrite.

Optical detection of disordered water within a protein cavity

Internal water molecules are important to protein structure and function, but positional disorder and low occupancies can obscure their detection by X-ray crystallography. Here, we show that water can be detected within the distal cavities of myoglobin mutants by subtle changes in the absorbance spectrum of pentacoordinate heme, even when the presence of solvent is not readily observed in the corresponding crystal structures. A well-defined, noncoordinated water molecule hydrogen bonded to the distal histidine (His64) is seen within the distal heme pocket in the crystal structure of wild type (wt) deoxymyoglobin. Displacement of this water decreases the rate of ligand entry into wt Mb, and we have shown previously that the entry of this water is readily detected optically after laser photolysis of MbCO complexes. However, for L29F and V68L Mb no discrete positions for solvent molecules are seen in the electron density maps of the crystal structures even though His64 is still present and slow rates of ligand binding indicative of internal water are observed. In contrast, time-resolved perturbations of the visible absorption bands of L29F and V68L deoxyMb generated after laser photolysis detect the entry and significant occupancy of water within the distal pockets of these variants. Thus, the spectral perturbation of pentacoordinate heme offers a potentially robust system for measuring nonspecific hydration of the active sites of heme proteins.

Haemoglobin polymorphisms affect the oxygen-binding properties in Atlantic cod populations

A major challenge in evolutionary biology is to identify the genes underlying adaptation. The oxygen-transporting haemoglobins directly link external conditions with metabolic needs and therefore represent a unique system for studying environmental effects on molecular evolution. We have discovered two haemoglobin polymorphisms in Atlantic cod populations inhabiting varying temperature and oxygen regimes in the North Atlantic. Three-dimensional modelling of the tetrameric haemoglobin structure demonstrated that the two amino acid replacements Met55beta1Val and Lys62beta1Ala are located at crucial positions of the alpha1beta1 subunit interface and haem pocket, respectively. The replacements are proposed to affect the oxygen-binding properties by modifying the haemoglobin quaternary structure and electrostatic feature. Intriguingly, the same molecular mechanism for facilitating oxygen binding is found in avian species adapted to high altitudes, illustrating convergent evolution in water- and air-breathing vertebrates to reduction in environmental oxygen availability. Cod populations inhabiting the cold Arctic waters and the low-oxygen Baltic Sea seem well adapted to these conditions by possessing the high oxygen affinity Val55-Ala62 haplotype, while the temperature-insensitive Met55-Lys62 haplotype predominates in the southern populations. The distinct distributions of the functionally different haemoglobin variants indicate that the present biogeography of this ecologically and economically important species might be seriously affected by global warming.

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).

The pH dependence of heme pocket hydration and ligand rebinding kinetics in photodissociated carbonmonoxymyoglobin

We monitored the occupancy of a functionally important non-coordinated water molecule in the distal heme pocket of sperm whale myoglobin over the pH range 4.3-9.4. Water occupancy was assessed by using time-resolved spectroscopy to detect the perturbation of the heme visible band absorption spectrum caused by water entry after CO photodissociation ( Goldbeck, R. A., Bhaskaran, S., Ortega, C., Mendoza, J. L., Olson, J. S., Soman, J., Kliger, D. S., and Esquerra, R. M. (2006) Proc. Natl. Acad. Sci. U. S. A. 103, 1254-1259 ). We found that the water occupancy observed during the time interval between ligand photolysis and diffusive recombination decreased by nearly 20% as the pH was lowered below 6. This decrease accounted for most of the concomitant increase in the observed CO bimolecular recombination rate constant, as the lower water occupancy presented a smaller kinetic barrier to CO entry into the pocket at lower pH. These results were consistent with a model in which the distal histidine, which stabilizes the water molecule within the distal pocket by accepting a hydrogen bond, tends to swing out of the pocket upon protonation and destabilize the water occupancy at low pH. Extrapolation of this model to lower pH suggests that the additional increase in ligand association rate constant observed previously in stopped-flow studies at pH 3 may also be due in part to reduced distal water occupancy concomitant with further His64 protonation and coupled protein conformational change.

Temperature-dependent heme kinetics with nonexponential binding and barrier relaxation in the absence of protein conformational substates

We present temperature-dependent kinetic measurements of ultrafast diatomic ligand binding to the "bare" protoheme (L(1)-FePPIX-L(2), where L(1) = H(2)O or 2-methyl imidazole and L(2) = CO or NO). We found that the binding of CO is temperature-dependent and nonexponential over many decades in time, whereas the binding of NO is exponential and temperature-independent. The nonexponential nature of CO binding to protoheme, as well as its relaxation above the solvent glass transition, mimics the kinetics of CO binding to myoglobin (Mb) but on faster time scales. This demonstrates that the nonexponential kinetic response observed for Mb is not necessarily due to the presence of protein conformational substates but rather is an inherent property of the solvated heme. The nonexponential kinetic data were analyzed by using a linear coupling model with a distribution of enthalpic barriers that fluctuate on slower time scales than the heme-CO recombination time. Below the solvent glass transition (T(g) approximately 180 K), the average enthalpic rebinding barrier for H(2)O-PPIX-CO was found to be approximately 1 kJ/mol. Above T(g), the barrier relaxes and is approximately 6 kJ/mol at 290 K. Values for the first two moments of the heme doming coordinate distribution extracted from the kinetic data suggest significant anharmonicity above T(g). In contrast to Mb, the protoheme shows no indication of the presence of "distal" enthalpic barriers. Moreover, the wide range of Arrhenius prefactors (10(9) to 10(11) s(-1)) observed for CO binding to heme under differing conditions suggests that entropic barriers may be an important source of control in this class of biochemical reactions.

The contribution of heme propionate groups to the conformational dynamics associated with CO photodissociation from horse heart myoglobin

Photoacoustic calorimetry and transient absorption spectroscopy were used to study conformational dynamics associated with CO photodissociation from horse heart myoglobin (Mb) reconstituted with either Fe protoporphyrin IX dimethylester (FePPDME), Fe octaethylporphyrin (FeOEP), or with native Fe protoporphyrin IX (FePPIX). The volume and enthalpy changes associated with the Fe-CO bond dissociation and formation of a transient deoxyMb intermediate for the reconstituted Mbs were found to be similar to those determined for native Mb (DeltaV1 = -2.5+/-0.6 ml mol(-1) and DeltaH1 = 8.1+/-3.0 kcal mol(-1)). The replacement of FePPIX by FeOEP significantly alters the conformational dynamics associated with CO release from protein. Ligand escape from FeOEP reconstituted Mb was determined to be roughly a factor of two faster (tau=330 ns) relative to native protein (tau=700 ns) and accompanying reaction volume and enthalpy changes were also found to be smaller (DeltaV2 = 5.4+/-2.5 ml mol(-1) and DeltaH2 = 0.7+/-2.2 kcal mol(-1)) than those for native Mb (DeltaV2 = 14.3+/-0.8 ml mol(-1) and DeltaH2 = 7.8+/-3.5 kcal mol(-1)). On the other hand, volume and enthalpy changes for CO release from FePPIX or FePPDME reconstituted Mb were nearly identical to those of the native protein. These results suggest that the hydrogen bonding network between heme propionate groups and nearby amino acid residues likely play an important role in regulating ligand diffusion through protein matrix. Disruption of this network leads to a partially open conformation of protein with less restricted ligand access to the heme binding pocket.