A Density Functional Theory Study of Conformers in the Ferrous CO Complex of Horseradish Peroxidase with Distinct Fe−C−O Configurations

Vibrational Stark Effects on Carbonyl, Nitrile, and Nitrosyl Compounds Including Heme Ligands, CO, CN, and NO, Studied with Density Functional Theory

Changes in the matrix electric field in a protein, due for example to mutations or structural fluctuations, can be correlated with changes in the vibrational transition frequencies of suitable chromophores measured by IR spectroscopy through the Stark tuning rate. To make this correlation, the Stark tuning rate must be known from experiment or theory. In this paper, density functional theory at the B3LYP/TZV level of theory is used to compute the Stark tuning rate of adducts of heme porphyrin, namely, -CO, -CN, and -NO+ compounds. The results are compared with the corresponding vibrational frequencies-field dependencies from vibrational Stark spectroscopy of heme-proteins. The zero-field computed Stark tuning rate of 1.3 cm-1/MV/cm for heme-CO is in agreement with experiment, where typically the rate is 2.4/f cm-1/MV/cm for myoglobin, where f is a local field correction between 1.1 and 1.4. Several small nitrile, carbonyl, and dinitrile molecules were studied to rationalize the findings for the heme-adducted models. Here, the higher B3LYP/6-311++G(2d,2p) level of theory could be used so the agreement with recent experimental results is even better than for heme-adducted groups.

Time-resolved FTIR study of CO recombination with horseradish peroxidase

Vibrational changes associated with CO recombination to ferrous horseradish peroxidase were investigated by rapid-scan FTIR (Fourier-transform IR) spectroscopy in the 1200-2200 cm(-1) range. At pH 6.0, two conformers of bound CO are present that appear as negative bands at 1905 and 1934 cm(-1) in photolysis spectra. Their recombination rate constants are identical, confirming that they arise from two substates of bound CO that are in rapid thermal equilibrium, rather than from heterogeneous protein sites. A smaller positive band at 2134 cm(-1) also appears on photolysis and decays with the same rate constant, indicative of an intraprotein geminate site involved in recombination or, possibly, a weak-affinity surface CO-binding site. Other signals arising from protein and haem in the 1700-1200 cm(-1) range can also be time-resolved with similar kinetics.

Probing dynamics of complex molecular systems with ultrafast 2D IR vibrational echo spectroscopy

Ultrafast 2D IR vibrational echo spectroscopy is described and a number of experimental examples are given. Details of the experimental method including the pulse sequence, heterodyne detection, and determination of the absorptive component of the 2D spectrum are outlined. As an initial example, the 2D spectrum of the stretching mode of CO bound to the protein myoglobin (MbCO) is presented. The time dependence of the 2D spectrum of MbCO, which is caused by protein structural evolution, is presented and its relationship to the frequency-frequency correlation function is described and used to make protein structural assignments based on comparisons to molecular dynamics simulations. The 2D vibrational echo experiments on the protein horseradish peroxidase are presented. The time dependence of the 2D spectra of the enzyme in the free form and with a substrate bound at the active site are compared and used to examine the influence of substrate binding on the protein's structural dynamics. The application of 2D vibrational echo spectroscopy to the study of chemical exchange under thermal equilibrium conditions is described. 2D vibrational echo chemical exchange spectroscopy is applied to the study of formation and dissociation of organic solute-solvent complexes and to the isomerization around a carbon-carbon single bond of an ethane derivative.

Infrared absorption study of the heme pocket dynamics of carbonmonoxyheme proteins

The temperature dependencies of the infrared absorption CO bands of carboxy complexes of horseradish peroxidase (HRP(CO)) in glycerol/water mixture at pH 6.0 and 9.3 are interpreted using the theory of optical absorption bandshape. The bands' anharmonic behavior is explained assuming that there is a higher-energy set of conformational substates (CSS(h)), which are populated upon heating and correspond to the protein substates with disordered water molecules in the heme pocket. Analysis of the second moments of the CO bands of the carboxy complexes of myoglobin (Mb(CO)) and hemoglobin (Hb(CO)), and of HRP(CO) with benzohydroxamic acid (HRP(CO)+BHA), shows that the low energy CSS(h) exists also in the open conformation of Mb(CO), where the heme pocket is spacious enough to accommodate a water molecule. In the HRP(CO)+BHA and closed conformations of Mb(CO) and Hb(CO), the heme pocket is packed with BHA and different amino acids, the CSS(h) has much higher energy and is hardly populated even at the highest temperatures. Therefore only motions of these amino acids contribute to the band broadening. These motions are linked to the protein surface and frozen in the glassy matrix, whereas in the liquid solvent they are harmonic. Thus the second moment of the CO band is temperature-independent in glass and is proportional to the temperature in liquid. The temperature dependence of the second moment of the CO peak of HRP(CO) in the trehalose glass exhibits linear coupling to an oscillator. This oscillator can be a moving water molecule locked in the heme pocket in the whole interval of temperatures or a trehalose molecule located in the heme pocket.

Resonance Raman spectroscopy of oxoiron(IV) porphyrin π-cation radical and oxoiron(IV) hemes in peroxidase intermediates

The catalytic cycle intermediates of heme peroxidases, known as compounds I and II, have been of long standing interest as models for intermediates of heme proteins, such as the terminal oxidases and cytochrome P450 enzymes, and for non-heme iron enzymes as well. Reports of resonance Raman signals for compound I intermediates of the oxo-iron(IV) porphyrin pi-cation radical type have been sometimes contradictory due to complications arising from photolability, causing compound I signals to appear similar to those of compound II or other forms. However, studies of synthetic systems indicated that protein based compound I intermediates of the oxoiron(IV) porphyrin pi-cation radical type should exhibit vibrational signatures that are different from the non-radical forms. The compound I intermediates of horseradish peroxidase (HRP), and chloroperoxidase (CPO) from Caldariomyces fumago do in fact exhibit unique and characteristic vibrational spectra. The nature of the putative oxoiron(IV) bond in peroxidase intermediates has been under discussion in the recent literature, with suggestions that the Fe(IV)O unit might be better described as Fe(IV)-OH. The generally low Fe(IV)O stretching frequencies observed for proteins have been difficult to mimic in synthetic ferryl porphyrins via electron donation from trans axial ligands alone. Resonance Raman studies of iron-oxygen vibrations within protein species that are sensitive to pH, deuteration, and solvent oxygen exchange, indicate that hydrogen bonding to the oxoiron(IV) group within the protein environment contributes to substantial lowering of Fe(IV)O frequencies relative to those of synthetic model compounds.

A study of the horseradish peroxidase catalytic site by FTIR spectroscopy

Vibrational changes in the catalytic site of horseradish peroxidase were investigated by FTIR (Fourier-transform infrared) spectroscopy in the 1000-2500 cm(-1) range. Difference spectra were generated by photolysis of the haemII-CO compound at different pH/pD values. The spectra report on the fine structure around the catalytic site and show vibrational changes of protein backbone, amino acid residues and cofactors. Assignments of the FTIR vibrations can be made based upon known crystal structures, comparisons with absorption frequencies and extinction coefficients of model amino acids and cofactors, effects of H2O/2H2O exchange and changes of pH/pD. Concomitant with the photolysis of the CO ligand are changes due to haem and protein vibrations, predominant among which are arginine and histidine residue vibrations.

Temperature excursion infrared (TEIR) spectroscopy used to study hydrogen bonding between water and biomolecules

Water is a highly polar molecule that is capable of making four H-bonding linkages. Stability and specificity of folding of water-soluble protein macromolecules are determined by the interplay between water and functional groups of the protein. Yet, under some conditions, water can be replaced with sugar or other polar protic molecules with retention of protein structure. Infrared (IR) spectroscopy allows one to probe groups on the protein that interact with solvent, whether the solvent is water, sugar or glycerol. The basis of the measurement is that IR spectral lines of functional groups involved in H-bonding show characteristic spectral shifts with temperature excursion, reflecting the dipolar nature of the group and its ability to H-bond. For groups involved in H-bonding to water, the stretching mode absorption bands shift to lower frequency, whereas bending mode absorption bands shift to higher frequency as temperature decreases. The results indicate increasing H-bonding and decreasing entropy occurring as a function of temperature, even at cryogenic temperatures. The frequencies of the amide group modes are temperature dependent, showing that as temperature decreases, the amide group H-bonds to water strengthen. These results are relevant to protein stability as a function of temperature. The influence of solvent relaxation is demonstrated for tryptophan fluorescence over the same temperature range where the solvent was examined by infrared spectroscopy.