How Does a Heme Carbene Differ from Diatomic Ligated (NO, CO, and CN-) Analogues in the Axial Bond?

Compared to well studied diatomic ligands (NO, CN-, CO), the axial bonds of carbene hemes is much less known although its significance in biological chemistry. The unusually large quadrupole splitting (Δ EQ = +2.2 mm·s-1) and asymmetric parameter (η = 0.9) of the five-coordinate heme carbene [Fe(TTP)(CCl2)], which is the largest among all known low spin ferrohemes, has driven investigations by means of Mössbauer effect Nuclear Resonance Vibrational Spectroscopy (NRVS). Three distinct measurements on one single crystal (two in-plane and one out-of-plane) have demonstrated comprehensive vibrational structures including stretch (429) and bending modes (472 cm-1) of the axial Fe-CCl2, and revealed iron vibrational anisotropy in three orthogonal directions for the first time. Frontier orbital analysis especially comparisons with diatomic analogues (NO, CN-, CO) suggest that CCl2, similar to NO, has led to strong but anisotropic π bonding in a ligand-based "4C"-coordinate which induced the vibrational anisotropies and very large Mössbauer parameters. This is contrasted to CN- and CO complexes which possess a porphyrin-based "4N"-coordinate electronic and vibrational structures due to inherent on-axis linear ligation.

What Can Be Learned from Nuclear Resonance Vibrational Spectroscopy: Vibrational Dynamics and Hemes.

Nuclear resonance vibrational spectroscopy (NRVS; also known as nuclear inelastic scattering, NIS) is a synchrotron-based method that reveals the full spectrum of vibrational dynamics for Mössbauer nuclei. Another major advantage, in addition to its completeness (no arbitrary optical selection rules), is the unique selectivity of NRVS. The basics of this recently developed technique are first introduced with descriptions of the experimental requirements and data analysis including the details of mode assignments. We discuss the use of NRVS to probe 57Fe at the center of heme and heme protein derivatives yielding the vibrational density of states for the iron. The application to derivatives with diatomic ligands (O2, NO, CO, CN-) shows the strong capabilities of identifying mode character. The availability of the complete vibrational spectrum of iron allows the identification of modes not available by other techniques. This permits the correlation of frequency with other physical properties. A significant example is the correlation we find between the Fe-Im stretch in six-coordinate Fe(XO) hemes and the trans Fe-N(Im) bond distance, not possible previously. NRVS also provides uniquely quantitative insight into the dynamics of the iron. For example, it provides a model-independent means of characterizing the strength of iron coordination. Prediction of the temperature-dependent mean-squared displacement from NRVS measurements yields a vibrational "baseline" for Fe dynamics that can be compared with results from techniques that probe longer time scales to yield quantitative insights into additional dynamical processes.

Wide-dynamic-range kinetic investigations of deep proton tunnelling in proteins.

Directional proton transport along 'wires' that feed biochemical reactions in proteins is poorly understood. Amino-acid residues with high pKa are seldom considered as active transport elements in such wires because of their large classical barrier for proton dissociation. Here, we use the light-triggered proton wire of the green fluorescent protein to study its ground-electronic-state proton-transport kinetics, revealing a large temperature-dependent kinetic isotope effect. We show that 'deep' proton tunnelling between hydrogen-bonded oxygen atoms with a typical donor-acceptor distance of 2.7-2.8 Šfully accounts for the rates at all temperatures, including the unexpectedly large value (2.5 × 10(9) s(-1)) found at room temperature. The rate-limiting step in green fluorescent protein is assigned to tunnelling of the ionization-resistant serine hydroxyl proton. This suggests how high-pKa residues within a proton wire can act as a 'tunnel diode' to kinetically trap protons and control the direction of proton flow.

3D Motions of Iron in Six-Coordinate {FeNO}(7) Hemes by Nuclear Resonance Vibration Spectroscopy.

The vibrational spectrum of a six-coordinate nitrosyl iron porphyrinate, monoclinic [Fe(TpFPP)(1-MeIm)(NO)] (TpFPP=tetra-para-fluorophenylporphyrin; 1-MeIm=1-methylimidazole), has been studied by oriented single-crystal nuclear resonance vibrational spectroscopy (NRVS). The crystal was oriented to give spectra perpendicular to the porphyrin plane and two in-plane spectra perpendicular or parallel to the projection of the FeNO plane. These enable assignment of the FeNO bending and stretching modes. The measurements reveal that the two in-plane spectra have substantial differences that result from the strongly bonded axial NO ligand. The direction of the in-plane iron motion is found to be largely parallel and perpendicular to the projection of the bent FeNO on the porphyrin plane. The out-of-plane Fe-N-O stretching and bending modes are strongly mixed with each other, as well as with porphyrin ligand modes. The stretch is mixed with v50 as was also observed for dioxygen complexes. The frequency of the assigned stretching mode of eight Fe-X-O (X=N, C, and O) complexes is correlated with the Fe-XO bond lengths. The nature of highest frequency band at ≈560 cm(-1) has also been examined in two additional new derivatives. Previously assigned as the Fe-NO stretch (by resonance Raman), it is better described as the bend, as the motion of the central nitrogen atom of the FeNO group is very large. There is significant mixing of this mode. The results emphasize the importance of mode mixing; the extent of mixing must be related to the peripheral phenyl substituents.

Recent advances in biosynthetic modeling of nitric oxide reductases and insights gained from nuclear resonance vibrational and other spectroscopic studies.

This Forum Article focuses on recent advances in structural and spectroscopic studies of biosynthetic models of nitric oxide reductases (NORs). NORs are complex metalloenzymes found in the denitrification pathway of Earth's nitrogen cycle where they catalyze the proton-dependent two-electron reduction of nitric oxide (NO) to nitrous oxide (N2O). While much progress has been made in biochemical and biophysical studies of native NORs and their variants, a clear mechanistic understanding of this important metalloenzyme related to its function is still elusive. We report herein UV-vis and nuclear resonance vibrational spectroscopy (NRVS) studies of mononitrosylated intermediates of the NOR reaction of a biosynthetic model. The ability to selectively substitute metals at either heme or nonheme metal sites allows the introduction of independent (57)Fe probe atoms at either site, as well as allowing the preparation of analogues of stable reaction intermediates by replacing either metal with a redox inactive metal. Together with previous structural and spectroscopic results, we summarize insights gained from studying these biosynthetic models toward understanding structural features responsible for the NOR activity and its mechanism. The outlook on NOR modeling is also discussed, with an emphasis on the design of models capable of catalytic turnovers designed based on close mimics of the secondary coordination sphere of native NORs.

Deep proton tunneling in the electronically adiabatic and non-adiabatic limits: comparison of the quantum and classical treatment of donor-acceptor motion in a protein environment.

Analytical models describing the temperature dependence of the deep tunneling rate, useful for proton, hydrogen, or hydride transfer in proteins, are developed and compared. Electronically adiabatic and non-adiabatic expressions are presented where the donor-acceptor (D-A) motion is treated either as a quantized vibration or as a classical "gating" distribution. We stress the importance of fitting experimental data on an absolute scale in the electronically adiabatic limit, which normally applies to these reactions, and find that vibrationally enhanced deep tunneling takes place on sub-ns timescales at room temperature for typical H-bonding distances. As noted previously, a small room temperature kinetic isotope effect (KIE) does not eliminate deep tunneling as a major transport channel. The quantum approach focuses on the vibrational sub-space composed of the D-A and hydrogen atom motions, where hydrogen bonding and protein restoring forces quantize the D-A vibration. A Duschinsky rotation is mandated between the normal modes of the reactant and product states and the rotation angle depends on the tunneling particle mass. This tunnel-mass dependent rotation contributes substantially to the KIE and its temperature dependence. The effect of the Duschinsky rotation is solved exactly to find the rate in the electronically non-adiabatic limit and compared to the Born-Oppenheimer (B-O) approximation approach. The B-O approximation is employed to find the rate in the electronically adiabatic limit, where we explore both harmonic and quartic double-well potentials for the hydrogen atom bound states. Both the electronically adiabatic and non-adiabatic rates are found to diverge at high temperature unless the proton coupling includes the often neglected quadratic term in the D-A displacement from equilibrium. A new expression is presented for the electronically adiabatic tunnel rate in the classical limit for D-A motion that should be useful to experimentalists working near room temperature. This expression also holds when a broad protein conformational distribution of D-A equilibrium distances dominates the spread of the D-A vibrational wavefunction.

Comprehensive Fe-ligand vibration identification in {FeNO}6 hemes.

Oriented single-crystal nuclear resonance vibrational spectroscopy (NRVS) has been used to obtain all iron vibrations in two {FeNO}(6) porphyrinate complexes, five-coordinate [Fe(OEP)(NO)]ClO4 and six-coordinate [Fe(OEP)(2-MeHIm)(NO)]ClO4. A new crystal structure was required for measurements of [Fe(OEP)(2-MeHIm)(NO)]ClO4, and the new structure is reported herein. Single crystals of both complexes were oriented to be either parallel or perpendicular to the porphyrin plane and/or axial imidazole ligand plane. Thus, the FeNO bending and stretching modes can now be unambiguously assigned; the pattern of shifts in frequency as a function of coordination number can also be determined. The pattern is quite distinct from those found for CO or {FeNO}(7) heme species. This is the result of unchanging Fe-N(NO) bonding interactions in the {FeNO}(6) species, in distinct contrast to the other diatomic ligand species. DFT calculations were also used to obtain detailed predictions of vibrational modes. Predictions were consistent with the intensity and character found in the experimental spectra. The NRVS data allow the assignment and observation of the challenging to obtain Fe-Im stretch in six-coordinate heme derivatives. NRVS data for this and related six-coordinate hemes with the diatomic ligands CO, NO, and O2 reveal a strong correlation between the Fe-Im stretch and Fe-N(Im) bond distance that is detailed for the first time.

The diagnostic vibrational signature of pentacoordination in heme carbonyls.

Heme-carbonyl complexes are widely exploited for the insight they provide into the structural basis of function in heme-based proteins, by revealing the nature of their bonded and nonbonded interactions with the protein. This report presents two novel results which clearly establish a FeCO vibrational signature for crystallographically verified pentacoordination. First, anisotropy in the NRVS density of states for ν(Fe-C) and δ(FeCO) in oriented single crystals of [Fe(OEP)(CO)] clearly reveals that the Fe-C stretch occurs at higher frequency than the FeCO bend and considerably higher than any previously reported heme carbonyl. Second, DFT calculations on a series of heme carbonyls reveal that the frequency crossover occurs near the weak trans O atom donor, furan. As ν(Fe-C) occurs at lower frequencies than δ(FeCO) in all heme protein carbonyls reported to date, the results reported herein suggest that they are all hexacoordinate.

Anisotropic iron motion in nitrosyl iron porphyrinates: natural and synthetic hemes.

The vibrational spectra of two five-coordinate nitrosyl iron porphyrinates, [Fe(OEP)(NO)] (OEP = dianion of 2,3,7,8,12,13,17,18-octaethylporphyrin) and [Fe(DPIX)(NO)] (DPIX = deuteroporphyrin IX), have been studied by oriented single-crystal nuclear resonance vibrational spectroscopy. Single crystals (both are in the triclinic crystal system) were oriented to give vibrational spectra perpendicular to the porphyrin plane. Additionally, two orthogonal in-plane measurements that were also either perpendicular or parallel to the projection of the FeNO plane onto the porphyrin plane yield the complete set of vibrations with iron motion. In addition to cleanly enabling the assignment of the FeNO bending and stretching modes, the measurements reveal that the two in-plane spectra from the parallel and perpendicular in-plane directions for both compounds have substantial differences. The assignment of these in-plane vibrations were aided by density functional theory predictions. The differences in the two in-plane directions result from the strongly bonded axial NO ligand. The direction of the in-plane iron motion is thus found to be largely parallel and perpendicular to the projection of the FeNO plane on the porphyrin plane. These axial ligand effects on the in-plane iron motion are related to the strength of the axial ligand-to-iron bond.

Spectroscopic and computational study of a nonheme iron nitrosyl center in a biosynthetic model of nitric oxide reductase.

A major barrier to understanding the mechanism of nitric oxide reductases (NORs) is the lack of a selective probe of NO binding to the nonheme FeB center. By replacing the heme in a biosynthetic model of NORs, which structurally and functionally mimics NORs, with isostructural ZnPP, the electronic structure and functional properties of the FeB nitrosyl complex was probed. This approach allowed observation of the first S=3/2 nonheme {FeNO}(7) complex in a protein-based model system of NOR. Detailed spectroscopic and computational studies show that the electronic state of the {FeNO}(7) complex is best described as a high spin ferrous iron (S=2) antiferromagnetically coupled to an NO radical (S=1/2) [Fe(2+)-NO(.)]. The radical nature of the FeB -bound NO would facilitate N-N bond formation by radical coupling with the heme-bound NO. This finding, therefore, supports the proposed trans mechanism of NO reduction by NORs.

Analytical Harmonic Vibrational Frequencies for the Green Fluorescent Protein Computed with ONIOM: Chromophore Mode Character and Its Response to Environment.

A systematic comparison of different environmental effects on the vibrational modes of the 4-hydroxybenzylidene-2,3-dimethylimidazolinone (HBDI) chromophore using the ONIOM method allows us to model how the molecule's spectroscopic transitions are modified in the Green Fluorescent Protein (GFP). ONIOM(QM:MM) reduces the expense of normal mode calculations when computing the majority of second derivatives only at the MM level. New developments described here for the efficient solution of the CPHF equations, including contributions from electrostatic interactions with environment charges, mean that QM model systems of ∼100 atoms can be embedded within a much larger MM environment of ∼5000 atoms. The resulting vibrational normal modes, their associated frequencies, and dipole derivative vectors have been used to interpret experimental difference spectra (GFPI2-GFPA), chromophore vibrational Stark shifts, and changes in the difference between electronic and vibrational transition dipoles (mode angles) in the protein environment.

Probing heme vibrational anisotropy: an imidazole orientation effect?

The complete iron vibrational spectrum of the five-coordinate high-spin complex [Fe(OEP)(2-MeHIm)], where OEP = octaethylporphyrinato and 2-MeHIm = 2-methylimidazole, has been obtained by oriented single-crystal nuclear resonance vibrational spectroscopy (NRVS) data. Measurements have been made in three orthogonal directions, which provides quantitative information for all iron motion. These experimental data, buttressed by density functional theory (DFT) calculations, have been used to define the effects of the axial ligand orientation. Although the axial imidazole removes the degeneracy in the in-plane vibrations, the imidazole orientation does not appear to control the direction of the in-plane iron motion. This is in contrast to the effect of the imidazolate ligand, as defined by DFT calculations, which does have substantial effects on the direction of the in-plane iron motion. The axial NO ligand has been found to have the strongest orientational effect (Angew. Chem., Int. Ed., 2010, 49, 4400). Thus the strength of the directional properties are in the order NO > imidazolate > imidazole, consistent with the varying strength of the Fe-ligand bond.

Effects of imidazole deprotonation on vibrational spectra of high-spin iron(II) porphyrinates.

The effects of the deprotonation of coordinated imidazole on the vibrational dynamics of five-coordinate high-spin iron(II) porphyrinates have been investigated using nuclear resonance vibrational spectroscopy. Two complexes have been studied in detail with both powder and oriented single-crystal measurements. Changes in the vibrational spectra are clearly related to structural differences in the molecular structures that occur when imidazole is deprotonated. Most modes involving the simultaneous motion of iron and imidazolate are unresolved, but the one mode that is resolved is found at higher frequency in the imidazolates. These out-of-plane results are in accord with earlier resonance Raman studies of heme proteins. We also show the imidazole vs imidazolate differences in the in-plane vibrations that are not accessible to resonance Raman studies. The in-plane vibrations are at lower frequency in the imidazolate derivatives; the doming mode shifts are inconclusive. The stiffness, an experimentally determined force constant that averages the vibrational details to quantify the nearest-neighbor interactions, confirms that deprotonation inverts the relative strengths of axial and equatorial coordination.

Ground-state proton transfer in the photoswitching reactions of the fluorescent protein Dronpa.

The reversible photoswitching between the 'on' and 'off' states of the fluorescent protein Dronpa involves photoisomerization as well as protein side-chain rearrangements, but the process of interconversion remains poorly characterized. Here we use time-resolved infrared measurements to monitor the sequence of these structural changes, but also of proton transfer events, which are crucial to the development of spectroscopic contrast. Light-induced deprotonation of the chromophore phenolic oxygen in the off state is a thermal ground-state process, which follows ultrafast (9 ps) trans-cis photoisomerization, and so does not involve excited-state proton transfer. Steady-state infrared difference measurements exclude protonation of the imidazolinone nitrogen in both the on and off states. Pump-probe infrared measurements of the on state reveal a weakening of the hydrogen bonding between Arg66 and the chromophore C=O, which could be central to initiating structural rearrangement of Arg66 and His193 coinciding with the low quantum yield cis-trans photoisomerization.

Vibrational probes and determinants of the S = 0 ⇌ S = 2 spin crossover in five-coordinate [Fe(TPP)(CN)]-.

The low-frequency vibrational characterization of the spin-crossover complex, five-coordinate cyano(tetraphenylporphyrinato)iron(II), [Fe(TPP)(CN)](-), is reported. Nuclear resonance vibrational spectroscopy has been used to measure all low-frequency vibrations involving iron at several temperatures; this yields vibrational spectra of both the low- (S = 0) and high-spin (S = 2) states. Multitemperature oriented single-crystal measurements facilitate assignments of the vibrational character of all modes and are consistent with the DFT-predicted spectra. The availability of the entire iron vibrational spectrum allows for the complete correlation of the modes between the two spin states. These data demonstrate that not only do the frequencies of the vibrations shift to lower values for the high-spin species as would be expected owing to the weaker bonds in the high-spin state, but also the mixing of iron modes with ligand modes changes substantially. Diagrams illustrating the changing character of the modes and their correlation are given. The reduced iron-ligand frequencies are the primary factor in the entropic stabilization of the high-spin state responsible for the spin crossover.

Vibrational dynamics of oxygenated heme proteins.

Advanced spectroscopic techniques coupled with DFT calculations reveal the vibrational dynamics of the iron in stable dioxygen complexes with myoglobin and with a mutant engineered to model the catalytic site of heme-copper oxidases. The unprecedented level of detail will constrain computational modelling of reactions with oxygen.

Spectroscopic characterization of (57)Fe-enriched cytochrome c.

Investigation of the heme iron dynamics in cytochrome c with Mössbauer spectroscopy and especially nuclear resonance vibrational spectroscopy requires the replacement of the natural abundant heme iron with the (57)Fe isotope. For demetallization, we use a safer and milder ferrous sulfate-hydrochloric acid method in addition to the harsher commonly used hydrofluoric acid-based procedure. The structural integrity of the (57)Fe-reconstituted protein in both oxidation states is confirmed from absorption spectra and a detailed analysis of the rich resonance Raman spectra. These results reinforce the application of metal-substituted heme c proteins as reliable models for the native proteins.

Structural changes that occur upon photolysis of the Fe(II)(a3)-CO complex in the cytochrome ba(3)-oxidase of Thermus thermophilus: a combined X-ray crystallographic and infrared spectral study demonstrates CO binding to Cu(B).

The purpose of the work was to provide a crystallographic demonstration of the venerable idea that CO photolyzed from ferrous heme-a(3) moves to the nearby cuprous ion in the cytochrome c oxidases. Crystal structures of CO-bound cytochrome ba(3)-oxidase from Thermus thermophilus, determined at ~2.8-3.2Å resolution, reveal a Fe-C distance of ~2.0Å, a Cu-O distance of 2.4Å and a Fe-C-O angle of ~126°. Upon photodissociation at 100K, X-ray structures indicate loss of Fe(a3)-CO and appearance of Cu(B)-CO having a Cu-C distance of ~1.9Å and an O-Fe distance of ~2.3Å. Absolute FTIR spectra recorded from single crystals of reduced ba(3)-CO that had not been exposed to X-ray radiation, showed several peaks around 1975cm(-1); after photolysis at 100K, the absolute FTIR spectra also showed a significant peak at 2050cm(-1). Analysis of the 'light' minus 'dark' difference spectra showed four very sharp CO stretching bands at 1970cm(-1), 1977cm(-1), 1981cm(-1), and 1985cm(-1), previously assigned to the Fe(a3)-CO complex, and a significantly broader CO stretching band centered at ~2050cm(-1), previously assigned to the CO stretching frequency of Cu(B) bound CO. As expected for light propagating along the tetragonal axis of the P4(3)2(1)2 space group, the single crystal spectra exhibit negligible dichroism. Absolute FTIR spectrometry of a CO-laden ba(3) crystal, exposed to an amount of X-ray radiation required to obtain structural data sets before FTIR characterization, showed a significant signal due to photogenerated CO(2) at 2337cm(-1) and one from traces of CO at 2133cm(-1); while bands associated with CO bound to either Fe(a3) or to Cu(B) in "light" minus "dark" FTIR difference spectra shifted and broadened in response to X-ray exposure. In spite of considerable radiation damage to the crystals, both X-ray analysis at 2.8 and 3.2Å and FTIR spectra support the long-held position that photolysis of Fe(a3)-CO in cytochrome c oxidases leads to significant trapping of the CO on the Cu(B) atom; Fe(a3) and Cu(B) ligation, at the resolutions reported here, are otherwise unaltered.

Nuclear resonance vibrational spectra of five-coordinate imidazole-ligated iron(II) porphyrinates.

Nuclear resonance vibrational spectra have been obtained for six five-coordinate imidazole-ligated iron(II) porphyrinates, [Fe(Por)(L)] (Por = tetraphenylporphyrinate, octaethylporphyrinate, tetratolylporphyrinate, or protoporphyrinate IX and L = 2-methylimidazole or 1,2-dimethylimidazole). Measurements have been made on both powder and oriented crystal samples. The spectra are dominated by strong signals around 200-300 cm(-1). Although the in-plane and out-of-plane vibrations are seriously overlapped, oriented crystal spectra allow their deconvolution. Thus, oriented crystal experimental data, along with density functional theory (DFT) calculations, enable the assignment of key vibrations in the spectra. Molecular dynamics are also discussed. The nature of the Fe-N(Im) vibrations has been elaborated further than was possible from resonance Raman studies. Our study suggests that the Fe motions are coupled with the porphyrin core and peripheral groups motions. Both peripheral groups and their conformations have significant influence on the vibrational spectra (position and shape).

New perspectives on iron-ligand vibrations of oxyheme complexes.

We report our studies of the vibrational dynamics of iron for three imidazole-ligated oxyheme derivatives that mimic the active sites of histidine-ligated heme proteins complexed with dioxygen. The experimental vibrational data are obtained from nuclear resonance vibrational spectroscopy (NRVS) measurements conducted on both powder samples and oriented single crystals, and which includes several in-plane (ip) and out-of-plane (oop) measurements. Vibrational spectral assignments have been made through a combination of the oriented sample spectra and predictions based on density functional theory (DFT) calculations. The two Fe-O(2) modes that have been previously observed by resonance Raman spectroscopy in heme proteins are clearly shown to be very strongly mixed and are not simply either a bending or stretching mode. In addition, a third Fe-O(2) mode, not previously reported, has been identified. The long-sought Fe-Im stretch, not observed in resonance Raman spectra, has been identified and compared with the frequencies observed for the analogous CO and NO species. The studies also suggest that the in-plane iron motion is anisotropic and is controlled by the orientation of the Fe-O(2) group and not sensitive to the in-plane Fe-N(p) bonds and/or imidazole orientations.