Structural dynamics of myoglobin: ligand migration among protein cavities studied by Fourier transform infrared/temperature derivative spectroscopy.

Fourier transform infrared (FTIR) spectroscopy in the CO stretch bands combined with temperature derivative spectroscopy (TDS) was used to characterize intermediate states obtained by photolysis of two sperm whale mutant myoglobins, YQR (L29(B10)Y, H64(E7)Q, T67(E10)R) and YQRF (with an additional I107(G8)F replacement). Both mutants assume two different bound-state conformations, A(0) and A(3), which can be distinguished by their different CO bands near 1965 and 1933 cm(-1). They most likely originate from different conformations of the Gln-64 side chain. Within each A substate, a number of photoproduct states have been characterized on the basis of the temperature dependence of recombination in TDS experiments. Different locations and orientations of the ligand within the protein can be distinguished by the infrared spectra of the photolyzed CO. Recombination from the primary docking site, B, near the heme dominates below 50 K. Above 60 K, ligand rebinding occurs predominantly from a secondary docking site, C', in which the CO is trapped in the Xe4 cavity on the distal side, as shown by crystallography of photolyzed YQR and L29W myoglobin CO. Another kinetic state (C") has been identified from which rebinding occurs around 130 K. Moreover, a population appearing above the solvent glass transition at approximately 180 K (D state) is assigned to rebinding from the Xe1 cavity, as suggested by the photoproduct structure of the L29W sperm whale myoglobin mutant. For both the YQR and YQRF mutants, rebinding from the B sites near the heme differs for the two A substates, supporting the view that the return of the ligand from the C', C", and D states is not governed by the recombination barrier at the heme iron but rather by migration to the active site. Comparison of YQR and YQRF shows that access to the Xe4 site (C') is severely restricted by introduction of the bulky Phe side chain at position 107.     

Structural dynamics of myoglobin: an infrared kinetic study of ligand migration in mutants YQR and YQRF

Recombination of carbon monoxide to myoglobin mutants YQR and YQRF was studied using transient infrared absorption spectroscopy and Fourier transform infrared-temperature derivative spectroscopy (FTIR-TDS). Photoproduct states B, C', C" and D associated with ligands residing in different protein cavities have been identified. After photolysis, ligands migrate to primary docking site B and subsequently rebind or escape to a secondary site (C) within the Xe4 cavity. For YQR, a global analysis of the isothermal rebinding kinetics below 160 K and the TDS data reveal a correlation between the enthalpy barriers governing the two processes. Above 120 K, a protein conformational change in both YQR and YQRF converts photoproduct C' into C" with markedly slowed kinetics. Above approximately 180 K, ligands migrate to the proximal Xe1 site (D) and also exit into the solvent, from where they rebind in a bimolecular reaction.     

Does CO2 permeate through aquaporin-1?

Aquaporins facilitate water permeation across biological membranes. Additionally, glycerol and other small neutral solutes are permeated by related aquaglyceroporins. The role of aquaporins in gas permeation has been a long-standing and controversially discussed issue. We present an extensive set of atomistic molecular dynamics simulations that address the question of CO(2) permeation through human aquaporin-1. Free energy profiles derived from the simulations display a barrier of approximately 23 kJ/mol in the aromatic/arginine constriction region of the water pore, whereas a barrier of approximately 4 kJ/mol was observed for a palmitoyloleoylphosphatidylethanolamine lipid bilayer membrane. The results indicate that significant aquaporin-1-mediated CO(2) permeation is to be expected only in membranes with a low intrinsic CO(2) permeability.     

Ligand migration and protein fluctuations in myoglobin mutant L29W

We have determined eight X-ray structures of myoglobin mutant L29W at various experimental conditions. In addition, infrared spectroscopic experiments are presented, which are discussed in the light of the X-ray structures. Two distinct conformations of the CO-ligated protein were identified, giving rise to two stretching bands of heme-bound CO. If L29W MbCO crystals are illuminated around 180 K, a deoxy species is formed. The CO molecules migrate to the proximal side of the heme and remain trapped in the so-called Xe1 cavity upon temperature decrease to 105 K. The structure of this photoproduct is almost identical to the equilibrium high-temperature deoxy Mb structure. If the temperature is cycled to increasingly higher values, CO recombination is observed. Three intermediate structures have been determined during the rebinding process. Efficient recombination occurs only above 180 K, the characteristic temperature for the onset of protein dynamics. Rebinding is remarkably slow because bulky residues His64 and Trp29 block important migration pathways of the CO molecule.     

Structural dynamics of myoglobin: FTIR-TDS study of NO migration and binding

By using Fourier transform infrared photolysis difference spectroscopy combined with temperature derivative spectroscopy at cryogenic temperatures, we have measured infrared spectra of the stretching absorption on nitric oxide (NO) in the heme-bound and photodissociated states of ferrous and ferric nitrosyl myoglobin (MbNO) and a few site-specific Mb mutants. The NO absorption was utilized as a sensitive local probe of ligand interactions with active-site residues and movements within the protein. By comparison with results obtained in previous spectroscopic and structural studies of carbonmonoxy myoglobin (MbCO), the MbNO data were interpreted in structural terms. In the NO-bound state, conformational heterogeneity was inferred from the appearance of multiple bands, arising from different electrostatic interactions with active site residues, most importantly, His-64. In ferrous MbNO, a primary photoproduct site similar to site B of MbCO was found, as indicated by a characteristic NO stretching spectrum. In ferric MbNO, the His-64 side chain appears to interfere with trapping of NO in this site; only a very weak photoproduct spectrum was observed in Mb variants in which His-64 was present. Upon extended illumination, the photoproduct spectrum changed in a characteristic way, indicating that NO readily migrates to a secondary docking site C, the Xe4 cavity, in which the ligand performs librational motions on the picosecond time scale. This docking site may play a role in the physiological NO scavenging reaction. Surprisingly, NO cannot be trapped at all in secondary docking site D, the Xe1 cavity.     

Theoretical investigation of infrared spectra and pocket dynamics of photodissociated carbonmonoxy myoglobin

Molecular dynamics simulations of the photodissociated state of carbonmonoxy myoglobin (MbCO) are presented using a fluctuating charge model for CO. A new three-point charge model is fitted to high-level ab initio calculations of the dipole and quadrupole moment functions taken from the literature. The infrared spectrum of the CO molecule in the heme pocket is calculated using the dipole moment time autocorrelation function and shows good agreement with experiment. In particular, the new model reproduces the experimentally observed splitting of the CO absorption spectrum. The splitting of 3-7 cm(-1) (compared to the experimental value of 10 cm(-1)) can be directly attributed to the two possible orientations of CO within the docking site at the edge of the distal heme pocket (the B states), as previously suggested on the basis of experimental femtosecond time-resolved infrared studies. Further information on the time evolution of the position and orientation of the CO molecule is obtained and analyzed. The calculated difference in the free energy between the two possible orientations (Fe...CO and Fe...OC) is 0.3 kcal mol(-1) and agrees well with the experimentally estimated value of 0.29 kcal mol(-1). A comparison of the new fluctuating charge model with an established fixed charge model reveals some differences that may be critical for the correct prediction of the infrared spectrum and energy barriers. The photodissociation of CO from the myoglobin mutant L29F using the new model shows rapid escape of CO from the distal heme pocket, in good agreement with recent experimental data. The effect of the protein environment on the multipole moments of the CO ligand is investigated and taken into account in a refined model. Molecular dynamics simulations with this refined model are in agreement with the calculations based on the gas-phase model. However, it is demonstrated that even small changes in the electrostatics of CO alter the details of the dynamics.     

A Spectroscopic Study of Structural Heterogeneity and Carbon Monoxide Binding in Neuroglobin

Neuroglobin (Ngb) is a small globular protein that binds diatomic ligands like oxygen, carbon monoxide (CO) and nitric oxide at a heme prosthetic group. We have performed FTIR spectroscopy in the infrared stretching bands of CO and flash photolysis with monitoring in the electronic heme absorption bands to investigate structural heterogeneity at the active site of Ngb and its effects on CO binding and migration at cryogenic temperatures. Four CO stretching bands were identified; they correspond to discrete conformations that differ in structural details and CO binding properties. Based on a comparison of bound-state and photoproduct IR spectra of the wild-type protein, Ngb distal pocket mutants and myoglobin, we have provided structural interpretations of the conformations associated with the different CO bands. We have also studied ligand migration to the primary docking site, B. Rebinding from this site is governed by very low enthalpy barriers (∼1 kJ/mol), indicating an extremely reactive heme iron. Moreover, we have observed ligand migration to a secondary docking site, C, from which CO rebinding involves higher enthalpy barriers.     

Direct probing of ion pair formation using a symmetric triangulenium dye

The 2,6,10-tris(dialkylamino)trioxatriangulenium dyes (ATOTA(+)) are highly stabilised cationic chromophores with D(3h) symmetry. The symmetry gives rise to a degeneracy of the main electronic transition. In low polarity solvents significant splitting of this degenerate transition is observed and assigned to ion pair formation. Ion pairing of the 2,6,10-tris(dioctylamino)trioxatriangulenium ion with Cl(-), BF(4)(-), PF(6)(-) and TRISPHAT anions was studied using absorption spectroscopy. A clear correlation is found between the size of the anion and the splitting of the ATOTA(+) transitions. In benzene the Cl(-) salt displays a splitting of 1955 cm(-1), while the salt of the much larger TRISPHAT ion has a splitting of 1543 cm(-1). TD-DFT calculations confirm the splitting of the states and provide a detailed insight into the electronic structure of the ion pairs. The different degree of splitting in different ion pairs is found to correlate with the magnitude of the electric field generated in each ion pair, thus leading to the conclusion that the effect seen is an internal Stark effect. By insertion of an amphiphilic derivative of the ATOTA(+) chromophore in an oriented lamellar liquid crystal, it was possible to resolve the two bands of the double peak spectrum and show their perpendicular orientation in the molecular framework, as predicted by the calculations.     

Noncontact dipole effects on channel permeation. V. Computed potentials for fluorinated gramicidin

Experimental and theoretical calculations indicate that the dipole moment of the four Trp side chains in gramicidin A (gA) channels modify channel conductance through long-range electrostatic interactions. Electrostatic ion/side-chain interaction energies along the channel were computed with CHARMM using ab initio atom charges for native and 4-, 5-, or 6-fluorinated Trp side chains. The bulk water reaction to the polar side chains was included using the method of images as implemented by, and channel waters in idealized structures were included. Ion/Trp interaction energies were approximately -0.6 kcal/mol throughout the channel for all four of the native Trp pairs. Channel waters produced a modest reduction in the magnitude of interactions, essentially offsetting images representing the bulk water outside the channel. The effects of side-chain fluorination depended on ring position and, to a lesser extent, residue number. Compared with native Trp, 5-fluorination reduces the translocation barrier with minor effects on the exit barrier. In contrast, 6-fluorination primarily reduces exit barrier. 4-Fluorination produces a more complex double-well energy profile. Effects of measured side-chain movements resulting from fluorination or change in lipid bilayer were negligible whereas thermal side chain librations cause large effects, especially in the region of the ion-binding sites.     

Resonance transduction of low level periodic signals by an enzyme: an oscillatory activation barrier model.

The overall rate of an enzyme catalyzed reaction is determined by the activation barrier of a rate-limiting step. If the barrier is oscillatory due to the intrinsic properties of a fluctuating enzyme, this enzymatic reaction will be influenced by a low level periodic electric field through the resonance transduction between the applied field and the oscillatory activation barrier. The ATP hydrolysis activity of a highly purified, detergent solubilized Ecto-ATPase from chicken oviduct was used to test the above concept. At 37 degrees C, this activity (1,800 mumols mg-1 min-1) was stimulated up to 47% (to 2,650 mumols mg-1 min-1) by an alternating electric field (AC), with a frequency window at 10 kHz. The maximal stimulation occurred at 5.0 V (peak-to-peak) cm-1. The potential drop across the dimension of the enzyme was approximately 10 microV (micelle diameter 20 nm). The activation barrier, or the Arrhenius activation energy, of the ATP splitting was measured to be 30 kT and the maximal barrier oscillation was calculated to be approximately 2.5 kT according to the oscillatory activation barrier (OAB) model. With the optimal AC field, full impact of the electric stimulation could be effected in much less than a second. The OAB model is many orders of magnitude more sensitive for deciphering low level periodic signals than the electroconformational coupling (ECC) model, although the latter has the ability to actively transduce energy while the former does not.(ABSTRACT TRUNCATED AT 250 WORDS)     

Carbinolamine formation and dehydration in a DNA repair enzyme active site

In order to suggest detailed mechanistic hypotheses for the formation and dehydration of a key carbinolamine intermediate in the T4 pyrimidine dimer glycosylase (T4PDG) reaction, we have investigated these reactions using steered molecular dynamics with a coupled quantum mechanics-molecular mechanics potential (QM/MM). We carried out simulations of DNA abasic site carbinolamine formation with and without a water molecule restrained to remain within the active site quantum region. We recovered potentials of mean force (PMF) from thirty replicate reaction trajectories using Jarzynski averaging. We demonstrated feasible pathways involving water, as well as those independent of water participation. The water-independent enzyme-catalyzed reaction had a bias-corrected Jarzynski-average barrier height of approximately (6.5 kcal mol(-1) (27.2 kJ mol(-1)) for the carbinolamine formation reaction and 44.5 kcal mol(-1) (186 kJ mol(-1)) for the reverse reaction at this level of representation. When the proton transfer was facilitated with an intrinsic quantum water, the barrier height was approximately 15 kcal mol(-1) (62.8 kJ mol(-1)) in the forward (formation) reaction and 19 kcal mol(-1) (79.5 kJ mol(-1)) for the reverse. In addition, two modes of unsteered (free dynamics) carbinolamine dehydration were observed: in one, the quantum water participated as an intermediate proton transfer species, and in the other, the active site protonated glutamate hydrogen was directly transferred to the carbinolamine oxygen. Water-independent unforced proton transfer from the protonated active site glutamate carboxyl to the unprotonated N-terminal amine was also observed. In summary, complex proton transfer events, some involving water intermediates, were studied in QM/MM simulations of T4PDG bound to a DNA abasic site. Imine carbinolamine formation was characterized using steered QM/MM molecular dynamics. Dehydration of the carbinolamine intermediate to form the final imine product was observed in free, unsteered, QM/MM dynamics simulations, as was unforced acid-base transfer between the active site carboxylate and the N-terminal amine.     

The effect of ligand dynamics on heme electronic transition band III in myoglobin.

Band III is a near-infrared electronic transition at ~13,000 cm(-1) in heme proteins that has been studied extensively as a marker of protein conformational relaxation after photodissociation of the heme-bound ligand. To examine the influence of the heme pocket structure and ligand dynamics on band III, we have studied carbon monoxide recombination in a variety of myoglobin mutants after photolysis at 3 K using Fourier transform infrared temperature-derivative spectroscopy with monitoring in three spectral ranges, (1) band III, the mid-infrared region of (2) the heme-bound CO, and (3) the photodissociated CO. Here we present data on mutant myoglobins V68F and L29W, which both exhibit pronounced ligand movements at low temperature. From spectral and kinetic analyses in the mid-infrared, a small number of photoproduct populations can be distinguished, differing in their distal heme pocket conformations and/or CO locations. We have decomposed band III into its individual photoproduct contributions. Each photoproduct state exhibits a different "kinetic hole-burning" (KHB) effect, a coupling of the activation enthalpy for rebinding to the position of band III. The analysis reveals that the heme pocket structure and the photodissociated CO markedly affect the band III transition. A strong kinetic hole-burning effect results only when the CO ligand resides in the docking site on top of the heme group. Migration of CO away from the heme group leads to an overall blue shift of band III. Consequently, band III can be used as a sensitive tool to study ligand dynamics after photodissociation in heme proteins.     

High pressure fourier transform infrared (FT-IR) spectroscopic studies on inducible nitric oxide (NO) synthase active site: a comparison to cytochrome p450CAM.

High pressure Fourier transform infrared (FT-IR) spectroscopy is performed for the first time to analyse the active site of inducible nitric oxide synthase (iNOSox) using the carbon monoxide (CO) heme iron ligand stretch mode (nuCO) as spectroscopic probe. A membrane-driven sapphire anvil high-pressure cell is used. Three major conformational substates exist in substrate-free iNOSox which are characterized by nuCO at approximately 1936, 1945 and 1952 cm(-1). High pressure favors the 1936 cm(-1) substate with a volume difference to the 1945 substate of approximately -21 cm3/mol. The pressure induced cytochrome P420 formation with a reaction volume of approximately -80 cm3/mol is observed. Arginine binding produces a very low nuCO at approximately 1905 cm(-1) caused by the H-bond from the substrate to CO. nuCO for the substates in the substrate-free and arginine-bound proteins shift linearly with pressure which is qualitatively similar to the observation on cytochrome P450cam. The slightly smaller positive slope of the shift in substrate-free iNOSox compared to substrate-free P450cam is interpreted as a slightly lesser compressible heme pocket. In contrast, the significant slower negative slope for arginine-bound iNOSox compared to camphor-bound P450cam results from the different kind of interactions to the CO ligand (electrostatic interaction in P450cam, H-bond in iNOSox).     

Vibrational Stark effects calibrate the sensitivity of vibrational probes for electric fields in proteins

Infrared spectroscopy is widely used to probe local environments and dynamics in proteins. The introduction of a unique vibration at a specific site of a protein or more complex assembly offers many advantages over observing the spectra of an unmodified protein. We have previously shown that infrared frequency shifts in proteins can arise from differences in the local electric field at the probe vibration. Thus, vibrational frequencies can be used to map electric fields in proteins at many sites or to measure the change in electric field due to a perturbation. The Stark tuning rate gives the sensitivity of a vibrational frequency to an electric field, and for it to be useful, the Stark tuning rate should be as large as possible. Vibrational Stark effect spectroscopy provides a direct measurement of the Stark tuning rate and allows a quantitative interpretation of frequency shifts. We present vibrational Stark spectra of several bond types, extending our work on nitriles and carbonyls and characterizing four additional bond types (carbon-fluorine, carbon-deuterium, azide, and nitro bonds) that are potential probes for electric fields in proteins. The measured Stark tuning rates, peak positions, and extinction coefficients provide the primary information needed to design amino acid analogues or labels to act as probes of local environments in proteins.     

Investigations of optical line shapes and kinetic hole burning in myoglobin.

We present the results of an extensive investigation of the optical line shapes of deoxymyoglobin (Mb), the ligand-bound form (MbCO), and the low-temperature photoproduct (Mb*). The thermal properties and the pH dependence of the Soret band and the near infrared band III (approximately 760 nm) are analyzed, taking into account the underlying vibrational properties of the absorption bands. The strong temperature dependence associated with the Soret band of MbCO and band III of Mb indicates significant coupling to low-frequency modes that may not be directly observed in the resonance Raman spectra. On the basis of analogous line-shape studies in a variety of heme systems, we assign the low-frequency coupling in MbCO to torsional motions of the CO molecule. The low-frequency mode coupled to band III (approximately 70 cm-1) is found to lie quite close to the value for the heme-doming motion (approximately 50 cm-1) calculated by using the kinetically determined value of the force constant (17 N/m). Significant inhomogeneous broadening in the Soret region of Mb and Mb* is found to be due to a "nonkinetic" coordinate that we associate with the orientation of the proximal histidine. A "kinetic" coordinate, associated with the equilibrium displacement of the iron atom from the porphyrin plane (a) is found to contribute to the inhomogeneous broadening of both the Soret band and band III. The relaxation of the heme as the system evolves from from Mb* to Mb is followed optically as a function of temperature, and a sharp transition temperature is found at 185 K. The blue shifts of the Soret band and band III as Mb* evolves to Mb are found to be nearly identical (delta v*ABS approximately 140 cm-1) and attributed to changes in the mean value of a between Mb* (a*0) and Mb (a0 = 0.45 A). A simple quadratic model for the coordinate coupling that simultaneously accounts for the observed shift, delta v*ABS, the low-temperature kinetics and the kinetic hole burning predicts a*0 = 0.2 +/- 0.05 A and EA = 16 +/- 2 kJ/mol for the room temperature Arrhenius barrier height at the heme. A simple quantitative method for the analysis of kinetic hole-burning experiments is also developed and applied to recent studies involving quaternary and subunit-specific hemoglobin structures.     

Structural dynamics controls nitric oxide affinity in nitrophorin 4

Nitrophorin 4 (NP4) is one of seven nitric oxide (NO) transporting proteins in the blood-sucking insect Rhodnius prolixus. In its physiological function, NO binds to a ferric iron centered in a highly ruffled heme plane. Carbon monoxide (CO) also binds after reduction of the heme iron. Here we have used Fourier transform infrared spectroscopy at cryogenic temperatures to study CO and NO binding and migration in NP4, complemented by x-ray cryo-crystallography on xenon-containing NP4 crystals to identify cavities that may serve as ligand docking sites. Multiple infrared stretching bands of the heme-bound ligands indicate different active site conformations with varying degrees of hydrophobicity. Narrow infrared stretching bands are observed for photodissociated CO and NO; temperature-derivative spectroscopy shows that these bands are associated with ligand docking sites close to the extremely reactive heme iron. No rebinding from distinct secondary sites was detected, although two xenon binding cavities were observed in the x-ray structure. Photolysis studies at approximately 200 K show efficient NO photoproduct formation in the more hydrophilic, open NP4 conformation. This result suggests that ligand escape is facilitated in this conformation, and blockage of the active site by water hinders immediate reassociation of NO to the ferric iron. In the closed, low-pH conformation, ligand escape from the active site of NP4 is prevented by an extremely reactive heme iron and the absence of secondary ligand docking sites.     

Infrared Study of Carbon Monoxide Migration among Internal Cavities of Myoglobin Mutant L29W

Myoglobin, a small globular heme protein that binds gaseous ligands such asO₂, CO and NO reversibly at the heme iron, provides an excellent modelsystem for studying structural and dynamic aspects of protein reactions. Flashphotolysis experiments, performed over wide ranges in time and temperature, reveal a complex ligand binding reaction with multiple kinetic intermediates, resulting from protein relaxation and movements of the ligand within the protein. Our recent studies of carbonmonoxy-myoglobin (MbCO) mutant L29W, using time-resolved infrared spectroscopy in combination with x-ray crystallography, have correlated kinetic intermediates with photoproduct structures that are characterized by the CO residing in different internal protein cavities, so-called xenon holes. Here we have used Fourier transform infrared temperature derivative spectroscopy (FTIR-TDS) to further examine the role of internal cavities in the dynamics. Different cavities can be accessed by the CO ligands at different temperatures, and characteristic infrared absorption spectra have been obtained for the different locations of the CO ligand within the protein, enabling us to monitor ligand migration through the protein as well as conformational changes of the protein.     

Estimation of cytoplasmic nitrate and its electrochemical potential in barley roots using 13NO3 and compartmental analysis

13NO3 was used to determine the intracellular compartmentation of NO3 in barley roots (Hordeum vulgare cv. Klondike), followed by a thermodynamic analysis of nitrate transport.Plants were grown in one-tenth Johnson's medium with 1 mol m3 NO3 (NO3-grown plants) or 1 mol m3 NH4NO3 (NH4NO3-grown plants).The cytoplasmic concentrations of NO3 in roots were only approx. 3-6 mol m3 (half-time for exchange approx. 21 s) in both NO3 and NH4NO3 plants. These pool sizes are consistent with published nitrate microelectrode data, but not with previous compartmental analyses.The electrochemical potential gradient for nitrate across the plasmalemma was +26 +/-1 kJ mol1 in both NO3- and NH4NO3-grown plants, indicating active uptake of nitrate. At an external pH of 6, the plasmalemma electrochemical potential for protons would be approx. -29 +/- 4 kJ mol1. If the cytoplasmic pH was 7.3 +/- 0.2, then 2H+/1NO3 cotransport, or a primary ATP-driven pump (2NO3/1ATP), are both thermodynamically possible. NO3 is also actively transported across the tonoplast (approx. +6 to 7 kJ mol1).     

Theoretical study of the discrimination between O(2) and CO by myoglobin

Combined quantum chemical and molecular mechanics geometry optimisations have been performed on myoglobin without or with O(2) or CO bound to the haem group. The results show that the distal histidine residue is protonated on the N(epsilon 2) atom and forms a hydrogen bond to the haem ligand both in the O(2) and the CO complexes. We have also re-refined the crystal structure of CO[bond]myoglobin by a combined quantum chemical and crystallographic refinement. Thereby, we probably obtain the most accurate available structure of the active site of this complex, showing a Fe[bond]C[bond]O angle of 171 degrees, and Fe[bond]C and C[bond]O bond lengths of 170-171 and 116-117 pm. The resulting structures have been used to calculate the strength of the hydrogen bond between the distal histidine residue and O(2) or CO in the protein. This amounts to 31-33 kJ/mol for O(2) and 2-3 kJ/mol for CO. The difference in hydrogen-bond strength is 21-22 kJ/mol when corrected for entropy effects. This is slightly larger than the observed discrimination between O(2) or CO by myoglobin, 17 kJ/mol. We have also estimated the strain of the active site inside the protein. It is 2-4 kJ/mol larger for the O(2) complex than for the CO complex, independent of which crystal structure the calculations are based on. Together, these results clearly show that myoglobin discriminates between O(2) and CO mainly by electrostatic interactions, rather than by steric strain.     

Ligand migration and binding in the dimeric hemoglobin of Scapharca inaequivalvis

Using Fourier transform infrared (FTIR) spectroscopy combined with temperature derivative spectroscopy (TDS) at cryogenic temperatures, we have studied CO binding to the heme and CO migration among cavities in the interior of the dimeric hemoglobin of Scapharca inaequivalvis (HbI) after photodissociation. By combining these studies with X-ray crystallography, three transient ligand docking sites were identified: a primary docking site B in close vicinity to the heme iron, and two secondary docking sites C and D corresponding to the Xe4 and Xe2 cavities of myoglobin. To assess the relevance of these findings for physiological binding, we also performed flash photolysis experiments on HbICO at room temperature and equilibrium binding studies with dioxygen. Our results show that the Xe4 and Xe2 cavities serve as transient docking sites for unbound ligands in the protein, but not as way stations on the entry/exit pathway. For HbI, the so-called histidine gate mechanism proposed for other globins appears as a plausible entry/exit route as well.