Role of the His64 residue on the properties of the Fe–CO and Fe–O2 bonds in myoglobin. A CHARMM/DFT study

Stationary and Time-Dependent Carbon Monoxide Stretching Mode Features in Carboxy Myoglobin: A Theoretical-Computational Reappraisal

The stationary and time-dependent infrared spectrum (IR) of the CO stretching mode (ν_CO) in carboxymyoglobin (MbCO), a longstanding problem of biophysical chemistry, has been modeled through a theoretical-computational method specifically designed for simulating quantum observables in complex atomic-molecular systems and based on a combined application of long time scale molecular dynamics simulations and quantum-chemical calculations. This study is basically focused on two aspects: (i) the origin of the stationary IR substates (termed as A₀, A₁, and A₃) and (ii) the modeling and the interpretation of the ν_CO energy relaxation. The results, strengthened by a more than satisfactory agreement with the experimental data, concisely indicate that (i) the conformational His64-FeCO relevant substates, i.e., characterized by the formation-disruption of the H-bond between the above moieties, are the main responsible of the presence of two distinct and well separated (A₀ and A₁/A₃) spectroscopic regions; (ii) the characteristic bimodal shape of the A₁/A₃ spectral region, according to our model, is the result of the fluctuation of the electric field pattern as provided by the protein-solvent framework perturbing the bound His64-CO-Heme complex; and (iii) the electric field pattern, in conjunction with the relatively high density of MbCO vibrational states, is also the main determinant of the ν_CO energy relaxation, characterizing its kinetic efficiency.

The unique structural features of carbonmonoxy hemoglobin from the sub-Antarctic fish Eleginops maclovinus

Tetrameric hemoglobins (Hbs) are prototypical systems for the investigations of fundamental properties of proteins. Although the structure of these proteins has been known for nearly sixty years, there are many aspects related to their function/structure that are still obscure. Here, we report the crystal structure of a carbonmonoxy form of the Hb isolated from the sub-Antarctic notothenioid fish Eleginops maclovinus characterised by either rare or unique features. In particular, the distal site of the α chain results to be very unusual since the distal His is displaced from its canonical position. This displacement is coupled with a shortening of the highly conserved E helix and the formation of novel interactions at tertiary structure level. Interestingly, the quaternary structure is closer to the T-deoxy state of Hbs than to the R-state despite the full coordination of all chains. Notably, these peculiar structural features provide a rationale for some spectroscopic properties exhibited by the protein in solution. Finally, this unexpected structural plasticity of the heme distal side has been associated with specific sequence signatures of various Hbs.

Probing the electronic and geometric structure of ferric and ferrous myoglobins in physiological solutions by Fe K-edge absorption spectroscopy

We present an iron K-edge X-ray absorption study of carboxymyoglobin (MbCO), nitrosylmyoglobin (MbNO), oxymyoglobin (MbO2), cyanomyoglobin (MbCN), aquomet myoglobin (metMb) and unligated myoglobin (deoxyMb) in physiological media. The analysis of the XANES region is performed using the full-multiple scattering formalism, implemented within the MXAN package. This reveals trends within the heme structure, absent from previous crystallographic and X-ray absorption analysis. In particular, the iron-nitrogen bond lengths in the porphyrin ring converge to a common value of about 2 Å, except for deoxyMb whose bigger value is due to the doming of the heme. The trends of the Fe-Nε (His93) bond length is found to be consistent with the effect of ligand binding to the iron, with the exception of MbNO, which is explained in terms of the repulsive trans effect. We derive a high resolution description of the relative geometry of the ligands with respect to the heme and quantify the magnitude of the heme doming in the deoxyMb form. Finally, time-dependent density functional theory is used to simulate the pre-edge spectra and is found to be in good agreement with the experiment. The XAS spectra typically exhibit one pre-edge feature which arises from transitions into the unoccupied dσ and dπ - πligand* orbitals. 1s → dπ transitions contribute weakly for MbO2, metMb and deoxyMb. However, despite this strong Fe d contribution these transitions are found to be dominated by the dipole (1s → 4p) moment due to the low symmetry of the heme environment.

Reactivity of deoxy- and oxyferrous dehaloperoxidase B from Amphitrite ornata: identification of compound II and its ferrous-hydroperoxide precursor

Dehaloperoxidase (DHP) from the terebellid polychaete Amphitrite ornata is a bifunctional enzyme that possesses both hemoglobin and peroxidase activities. The bifunctional nature of DHP as a globin peroxidase appears to be at odds with the traditional starting oxidation state for each individual activity. Namely, reversible oxygen binding is only mediated via a ferrous heme in globins, and peroxidase activity is initiated from ferric centers and to the exclusion of the oxyferrous oxidation state from the peroxidase cycle. Thus, to address what appears to be a paradox, herein we report the details of our investigations into the DHP catalytic cycle when initiated from the deoxy- and oxyferrous states using biochemical assays, stopped-flow UV-visible, and rapid-freeze-quench electron paramagnetic resonance spectroscopies, and anaerobic methods. We demonstrate the formation of Compound II directly from deoxyferrous DHP B upon its reaction with hydrogen peroxide and show that this occurs both in the presence and in the absence of trihalophenol. Prior to the formation of Compound II, we have identified a new species that we have preliminarily attributed to a ferrous-hydroperoxide precursor that undergoes heterolysis to generate the aforementioned ferryl intermediate. Taken together, the results demonstrate that the oxyferrous state in DHP is a peroxidase competent starting species, and an updated catalytic cycle for DHP is proposed in which the ferric oxidation state is not an obligatory starting point for the peroxidase catalytic cycle of dehaloperoxidase. The data presented herein provide a link between the peroxidase and oxygen transport activities, which furthers our understanding of how this bifunctional enzyme is able to unite its two inherent functions in one system.

Nature of the Fe-O2 bonding in oxy-myoglobin: effect of the protein

The nature of the Fe-O2 bonding in oxy-myoglobin was probed by theoretical calculations: (a) QM/MM (hybrid quantum mechanical/molecular mechanical) calculations using DFT/MM and CASSCF/MM methods and (b) gas-phase calculations using DFT (density functional theory) and CASSCF (complete active space self-consistent field) methods. Within the protein, the O2 is hydrogen bonded by His64 and the complex feels the bulk polarity of the protein. Removal of the protein causes major changes in the complex. Thus, while CASSCF/MM and DFT/MM are similar in terms of state constitution, degree of O2 charge, and nature of the lowest triplet state, the gas-phase CASSCF(g) species is very different. Valence bond (VB) analysis of the CASSCF/MM wave function unequivocally supports the Weiss bonding mechanism. This bonding arises by electron transfer from heme-Fe(II) to O2 and the so formed species coupled then to a singlet state Fe(III)-O2(-) that possesses a dative sigma(Fe-O) bond and a weakly coupled pi(Fe-O2) bond pair. The bonding mechanism in the gas phase is similar, but now the sigma(Fe-O) bond involves higher back-donation from O2(-) to Fe(III), while the constituents of pi(Fe-O2) bond pair have greater delocalization tails. The protein thus strengthens the Fe(III)-O2(-) character of the complex and thereby affects its bonding features and the oxygen binding affinity of Mb. The VB model is generalized, showing how the protein or the axial ligand of the oxyheme complex can determine the nature of its bonding in terms of the blend of the three bonding models: Weiss, Pauling, and McClure-Goddard.

A DFT study on the relative affinity for oxygen of the alpha and beta subunits of hemoglobin

DFT calculations are carried out on computational models of the active center of the alpha and beta subunits of hemoglobin in both its oxygenated (R) and deoxygenated (T) states. The computational models are defined by the full heme group, including all porphyrin substituents, and the four amino acids closer to it. The role of the protein environment is introduced by freezing the position of the alpha carbon atom of each of the four amino acids to the positions they have in the available PDB structures. Oxygen affinity is then evaluated by computing the energy difference between the optimized structures of the oxygenated and deoxygenated forms of each model. The results indicate a higher affinity of the alpha subunits over the beta ones. Analysis of the computed structures points out to the strength of the hydrogen bond between the distal histidine and the oxygen molecule as a key factor in discriminating the different systems.

Ultrafast dynamics of myoglobin without the distal histidine: stimulated vibrational echo experiments and molecular dynamics simulations

Ultrafast protein dynamics of the CO adduct of a myoglobin mutant with the polar distal histidine replaced by a nonpolar valine (H64V) have been investigated by spectrally resolved infrared stimulated vibrational echo experiments and molecular dynamics (MD) simulations. In aqueous solution at room temperature, the vibrational dephasing rate of CO in the mutant is reduced by approximately 50% relative to the native protein. This finding confirms that the dephasing of the CO vibration in the native protein is sensitive to the interaction between the ligand and the distal histidine. The stimulated vibrational echo observable is calculated from MD simulations of H64V within a model in which vibrational dephasing is driven by electrostatic forces. In agreement with experiment, calculated vibrational echoes show slower dephasing for the mutant than for the native protein. However, vibrational echoes calculated for H64V do not show the quantitative agreement with measurements demonstrated previously for the native protein.

Cytochrome c552 Mutants: Structure and Dynamics at the Active Site Probed by Multidimensional NMR and Vibration Echo Spectroscopy

Spectrally resolved infrared stimulated vibrational echo experiments are used to measure the vibrational dephasing of a CO ligand bound to the heme cofactor in two mutated forms of the cytochrome c552 from Hydrogenobacter thermophilus. The first mutant (Ht-M61A) is characterized by a single mutation of Met61 to an Ala (Ht-M61A), while the second variant is doubly modified to have Gln64 replaced by an Asn in addition to the M61A mutation (Ht-M61A/Q64N). Multidimensional NMR experiments determined that the geometry of residue 64 in the two mutants is consistent with a non-hydrogen-bonding and hydrogen-bonding interaction with the CO ligand for Ht-M61A and Ht-M61A/Q64N, respectively. The vibrational echo experiments reveal that the shortest time scale vibrational dephasing of the CO is faster in the Ht-M61A/Q64N mutant than that in Ht-M61A. Longer time scale dynamics, measured as spectral diffusion, are unchanged by the Q64N modification. Frequency-frequency correlation functions (FFCFs) of the CO are extracted from the vibrational echo data to confirm that the dynamical difference induced by the Q64N mutation is primarily an increase in the fast (hundreds of femtoseconds) frequency fluctuations, while the slower (tens of picoseconds) dynamics are nearly unaffected. We conclude that the faster dynamics in Ht-M61A/Q64N are due to the location of Asn64, which is a hydrogen bond donor, above the heme-bound CO. A similar difference in CO ligand dynamics has been observed in the comparison of the CO derivative of myoglobin (MbCO) and its H64V variant, which is caused by the difference in axial residue interactions with the CO ligand. The results suggest a general trend for rapid ligand vibrational dynamics in the presence of a hydrogen bond donor.

Dynamics of proteins encapsulated in silica sol-gel glasses studied with IR vibrational echo spectroscopy

Spectrally resolved infrared stimulated vibrational echo spectroscopy is used to measure the fast dynamics of heme-bound CO in carbonmonoxy-myoglobin (MbCO) and -hemoglobin (HbCO) embedded in silica sol-gel glasses. On the time scale of approximately 100 fs to several picoseconds, the vibrational dephasing of the heme-bound CO is measurably slower for both MbCO and HbCO relative to that of aqueous protein solutions. The fast structural dynamics of MbCO, as sensed by the heme-bound CO, are influenced more by the sol-gel environment than those of HbCO. Longer time scale structural dynamics (tens of picoseconds), as measured by the extent of spectral diffusion, are the same for both proteins encapsulated in sol-gel glasses compared to that in aqueous solutions. A comparison of the sol-gel experimental results to viscosity-dependent vibrational echo data taken on various mixtures of water and fructose shows that the sol-gel-encapsulated MbCO exhibits dynamics that are the equivalent of the protein in a solution that is nearly 20 times more viscous than bulk water. In contrast, the HbCO dephasing in the sol-gel reflects only a 2-fold increase in viscosity. Attempts to alter the encapsulating pore size by varying the molar ratio of silane precursor to water (R value) used to prepare the sol-gel glasses were found to have no effect on the fast or steady-state spectroscopic results. The vibrational echo data are discussed in the context of solvent confinement and protein-pore wall interactions to provide insights into the influence of a confined environment on the fast structural dynamics experienced by a biomolecule.