Identification of Conformational Substates in Oxymyoglobin through the pH-Dependence of the Low-Temperature Photoproduct Yield

The photochemistry of transition metal complexes using density functional theory

The use of density functional theory (DFT) and time-dependent DFT (TD-DFT) to study the photochemistry of metal complexes is becoming increasingly important among chemists. Computational methods provide unique information on the electronic nature of excited states and their atomic structure, integrating spectroscopy observations on transient species and excited-state dynamics. In this contribution, we present an overview on photochemically active transition metal complexes investigated by DFT. In particular, we discuss a representative range of systems studied up to now, which include CO- and NO-releasing inorganic and organometallic complexes, haem and haem-like complexes dissociating small diatomic molecules, photoactive anti-cancer Pt and Ru complexes, Ru polypyridyls and diphosphino Pt derivatives.

QM/MM vibrational mode tracking

Vibrational spectroscopy is a powerful tool to investigate the structure and dynamics of biomolecules. When small subsystems of large molecules such as active centers of enzymes are studied, quantum chemical calculations based on quantum mechanics/molecular mechanics (QM/MM) coupling schemes are a valuable means to interpret the spectra. The goal of this work is a methodological pilot study on how to selectively and thus efficiently extract certain vibrational information for extended molecular systems described by QM/MM methods. This is achieved by an extension of the mode tracking algorithm and a comparison with the partial Hessian diagonalization approach. After validating the methodology for the CO stretching vibration of 2-butanone and a delocalized CO stretch in acetylacetone, the stretching and bending modes of the CO ligand in CO myoglobin are tracked. Such systems represent an ideal application for mode tracking, because only a few strongly localized vibrations are sought for, while the large remainder of the molecule is of interest only as far as it affects these local vibrations. This influence is treated exactly by mode tracking.

Structure and dynamics of dioxygen bound to cobalt and iron heme

In this study we use ab initio molecular dynamics simulations to analyze the structure and dynamics of the oxygen ligand in models of the oxymyoglobin active site and its cobalt-substituted analog. Our calculations are performed for iron-porphyrin and cobalt-porphyrin complexes with imidazole and oxygen as axial ligands, and we investigate the effect of the distal histidine in the structure and dynamics of the metal-oxygen unit (MeO(2), Me = Fe, Co). We find that the interaction between the distal histidine and the oxygen ligand is stronger for the cobalt complex than for the iron one, consistent with the superoxide ion character of the bound O(2). The dynamics of the O(2) ligand can be described as oscillations of the O-O axis projection on the porphyrin plane within a porphyrin quadrant combined with frequent jumps from one quadrant to another. However, the ligand motion is significantly faster for CoO(2) compared to FeO(2). As a result, the iron complex shows localized ligand sites, whereas for cobalt several configurations are possible. This gives support to the highly dynamic motion of the oxygen ligand found in several experiments on cobalt oxymyoglobin and model complexes and underlines the higher mobility of the CoO(2) fragment compared to FeO(2).

Quantum Chemical Evaluation of Protein Control over Heme Ligation: CO/O2 Discrimination in Myoglobin

Control of O2 versus CO binding in myoglobin (Mb) is tuned by a distal histidine residue through steric and H-bonding interactions. These interactions have been evaluated via Car-Parrinello DFT calculations, whose efficiency allows full quantum mechanical treatment of the 13 closest residues surrounding the heme. The small (8 degrees ) deviation of the Fe-C-O bond angle from linearity results from the steric influence of a distal valine residue and not the distal histidine. H-bond energies were evaluated by replacing the distal histidine with the non-H-bonding residue isoleucine. Binding energies for CO and O2 decreased by 0.8 and 4.1 kcal/mol for MbCO and MbO2, in good agreement with experimental H-bond estimates. Ligand discrimination is dominated by distal histidine H-bonding, which is also found to stabilize a metastable side-on isomer of MbO2 that may play a key role in MbO2 photodynamics.

Role of Ligand Bending in the Photodissociation of O2 vs CO-heme: A Time-Dependent Density Functional Study

Time-dependent DFT calculations reveal a strong dependence of low-lying excited states on the <Fe-O-O angle of a model O2-heme adduct and ready access to a metastable side-on isomer. The predicted state crossings provide a pathway for internal conversion to the ground state and account for the low quantum yield for O2 dissociation, while the side-on isomer can account for the observation of an unphotolyzable fraction at low temperature. For CO-heme, the angle-sensitive states lie above the photoaccessed Q state and do not participate in photodissociation.

Kinetics of the reactions of nitrogen monoxide and nitrite with ferryl hemoglobin

Hemoglobin released in the circulation from ruptured red blood cells can be oxidized by hydrogen peroxide or peroxynitrite to generate the highly oxidizing iron(IV)oxo species HbFe(IV)z=O. Nitrogen monoxide, produced in large amounts by activated inducible nitric oxide synthase, can have indirect cytotoxic effects, mainly through the generation of peroxynitrite from its very fast reaction with superoxide. In the present work we have determined the rate constant for the reaction of HbFe(IV)z=O with NO(*), 2.4 x 10(7) M(-1)s(-1) at pH 7.0 and 20 degrees C. The reaction proceeds via the intermediate HbFe(III)ONO, which then dissociates to metHb and nitrite. As these products are not oxidizing and because of its large rate, the reaction of HbFe(IV)z=O with NO(*) may be important to remove the high valent form of hemoglobin, which has been proposed to be at least in part responsible for oxidative lesions. In addition, we have determined that the rate constant for the reaction of HbFe(IV)z=O with nitrite is significantly lower (7.5 x 10(2) M(-1)s(-1) at pH 7.0 and 20 degrees C), but increases with decreasing pH (1.8 x 10(3) M(-1)s(-1) at pH 6.4 and 20 degrees C). Thus, under acidic conditions as found in ischemic tissues, this reaction may also have a physiological relevance.

Myoglobin-CO conformational substate dynamics: 2D vibrational echoes and MD simulations

Two-dimensional (2D) infrared vibrational echoes were performed on horse heart carbonmonoxymyoglobin (MbCO) in water over a range of temperatures. The A(1) and A(3) conformational substates of MbCO are found to have different dephasing rates with different temperature dependences. A frequency-frequency correlation function derived from molecular dynamics simulations on MbCO at 298 K is used to calculate the vibrational echo decay. The calculated decay shows substantial agreement with the experimentally measured decays. The 2D vibrational echo probes protein dynamics and provides an observable that can be used to test structural assignments for the MbCO conformational substates.

Kinetic evidence for three photolyzable taxonomic conformational substates in oxymyoglobin

The kinetics of oxygen geminate binding with the taxonomic substates of MbO2 are reported. The maximum entropy method was used to analyze the rebinding kinetics of MbCO and MbO2 monitored in the Soret. The resulting rate distributions were found to consist of a small number of overlapping bands. A global parametric fit of a series of rate distributions recorded at several temperatures was performed using a Gaussian basis set to resolve the individual enthalpy distributions P(H). This approach was first validated by showing that the well-documented taxonomic substates of MbCO could be recovered. The method was then applied to MbO2. Three taxonomic substates were identified at pH 4.8, whereas only two of them contribute to oxygen geminate rebinding at pH 7.0. These findings show that, similarly to MbCO, MbO2 also exists as three photolyzable and kinetically different taxonomic substates and suggest reconsidering the issue of the photolysis quantum yield of MbO2.

Crystal structures of myoglobin-ligand complexes at near-atomic resolution

We have used x-ray crystallography to determine the structures of sperm whale myoglobin (Mb) in four different ligation states (unligated, ferric aquomet, oxygenated, and carbonmonoxygenated) to a resolution of better than 1.2 A. Data collection and analysis were performed in as much the same way as possible to reduce model bias in differences between structures. The structural differences among the ligation states are much smaller than previously estimated, with differences of <0.25 A root-mean-square deviation among all atoms. One structural parameter previously thought to vary among the ligation states, the proximal histidine (His-93) azimuthal angle, is nearly identical in all the ferrous complexes, although the tilt of the proximal histidine is different in the unligated form. There are significant differences, however, in the heme geometry, in the position of the heme in the pocket, and in the distal histidine (His-64) conformations. In the CO complex the majority conformation of ligand is at an angle of 18 +/- 3 degrees with respect to the heme plane, with a geometry similar to that seen in encumbered model compounds; this angle is significantly smaller than reported previously by crystallographic studies on monoclinic Mb crystals, but still significantly larger than observed by photoselection. The distal histidine in unligated Mb and in the dioxygenated complex is best described as having two conformations. Two similar conformations are observed in MbCO, in addition to another conformation that has been seen previously in low-pH structures where His-64 is doubly protonated. We suggest that these conformations of the distal histidine correspond to the different conformational substates of MbCO and MbO(2) seen in vibrational spectra. Full-matrix refinement provides uncertainty estimates of important structural parameters. Anisotropic refinement yields information about correlated disorder of atoms; we find that the proximal (F) helix and heme move approximately as rigid bodies, but that the distal (E) helix does not.