Resilience of the iron environment in heme proteins

Adiabatic Ligand Binding in Heme Proteins: Ultrafast Kinetics of Methionine Rebinding in Ferrous Cytochrome c

The dynamics of methionine geminate recombination following photodissociation in ferrous cytochrome c is investigated over a broad temperature range. The kinetic response, above the solvent glass transition ( T_g), is nearly monoexponential and displays a weak temperature dependence. Below T_g, the rebinding kinetics are nonexponential and can be explained using a quenched distribution of enthalpic rebinding barriers, arising from a relatively narrow distribution of heme out-of-plane displacements. The Arrhenius prefactor of this (Δ S = 2) reaction is ∼10¹¹ s^-1, which is similar to what has been found for the (Δ S = 1) NO binding reaction in heme proteins. This observation, along with other examples of ultrafast CO binding, provides strong evidence that ligand binding to heme is an adiabatic reaction with a spin-independent prefactor. In order to simultaneously account for the adiabatic nature of the reaction as well as the temperature dependence of both ultrafast CO and methionine geminate rebinding, it is proposed that a spin triplet state intersects and strongly couples to the reactant ( S = 2) and product ( S = 0) state surfaces in the transition state region along the reaction coordinate. It is also suggested that the nature of the intersecting triplet state and the reaction path may depend upon the proximity of the photolyzed ligand relative to the iron atom. At temperatures below ∼60 K, the kinetic data suggest that there is either an unexpected retardation of the heme photoproduct relaxation or that heavy atom quantum mechanical tunneling becomes significant.

A Nuclear Resonance Vibrational Spectroscopic Study of Oxy Myoglobins Reconstituted with Chemically Modified Heme Cofactors: Insights into the Fe-O2 Bonding and Internal Dynamics of the Protein

The molecular mechanism of O₂ binding to hemoglobin (Hb) and myoglobin (Mb) is a long-standing issue in the field of bioinorganic and biophysical chemistry. The nature of Fe-O₂ bond in oxy Hb and Mb had been extensively investigated by resonance Raman spectroscopy, which assigned the Fe-O₂ stretching bands at ∼570 cm^-1. However, resonance Raman assignment of the vibrational mode had been elusive due to the spectroscopic selection rule and to the limited information available about the ground-state molecular structure. Thus, nuclear resonance vibrational spectroscopy was applied to oxy Mbs reconstituted with ⁵⁷Fe-labeled native heme cofactor and two chemically modified ones. This advanced spectroscopy in conjunction with DFT analyses gave new insights into the nature of the Fe-O₂ bond of oxy heme by revealing the effect of heme peripheral substitutions on the vibrational dynamics of heme Fe atom, where the main Fe-O₂ stretching band of the native protein was characterized at ∼420 cm^-1.

SciPhon: a data analysis software for nuclear resonant inelastic X-ray scattering with applications to Fe, Kr, Sn, Eu and Dy

The synchrotron radiation technique of nuclear resonant inelastic X-ray scattering (NRIXS), also known as nuclear resonance vibrational spectroscopy or nuclear inelastic scattering, provides a wealth of information on the vibrational properties of solids. It has found applications in studies of lattice dynamics and elasticity, superconductivity, heme biochemistry, seismology, isotope geochemistry and many other fields. It involves probing the vibrational modes of solids by using the nuclear resonance of Mössbauer isotopes such as ⁵⁷Fe, ⁸³Kr, ¹¹⁹Sn, ¹⁵¹Eu and ¹⁶¹Dy. After data reduction, it provides the partial phonon density of states of the Mössbauer isotope that is investigated, as well as many other derived quantities such as the mean force constant of the chemical bonds and the Debye velocity. The data reduction is, however, not straightforward and involves removal of the elastic peak, normalization and Fourier-Log transformation. Furthermore, some of the quantities derived are highly sensitive to details in the baseline correction. A software package and several novel procedures to streamline and hopefully improve the reduction of the NRIXS data generated at sector 3ID of the Advanced Photon Source have been developed. The graphical user interface software is named SciPhon and runs as a Mathematica package. It is easily portable to other platforms and can be easily adapted for reducing data generated at other beamlines. Several tests and comparisons are presented that demonstrate the usefulness of this software, whose results have already been used in several publications. Here, the SciPhon software is used to reduce Kr, Sn, Eu and Dy NRIXS data, and potential implications for interpreting natural isotopic variations in those systems are discussed.

Long-Range Modulations of Electric Fields in Proteins

Electrostatic interactions are essential for controlling the protein structure and function. Whereas so far experimental and theoretical efforts focused on the effect of local electrostatics, this work aims at elucidating the long-range modulation of electric fields in proteins upon binding to charged surfaces. The study is based on cytochrome c (Cytc) variants carrying nitrile reporters for the vibrational Stark effect that are incorporated into the protein via genetic engineering and chemical modification. The Cytc variants were thoroughly characterized with respect to possible structural perturbations due to labeling. For the proteins in solution, the relative hydrogen bond occupancy and the calculated electric fields, both obtained from molecular dynamics (MD) simulations, and the experimental nitrile stretching frequencies were used to develop a relationship for separating hydrogen-bonding and non-hydrogen-bonding electric field effects. This relationship provides an excellent description for the stable Cytc variants in solution. For the proteins bound to Au electrodes coated with charged self-assembled monolayers (SAMs), the underlying MD simulations can only account for the electric field changes Δ E_ads due to the formation of the electrostatic SAM-Cytc complexes but not for the additional contribution, Δ E_int, representing the consequences of the potential drops over the electrode/SAM/protein interfaces. Both Δ E_ads and Δ E_int, determined at distances between 20 and 30 Å with respect to the SAM surface, are comparable in magnitude to the non-hydrogen-bonding electric field in the unbound protein. This long-range modulation of the internal electric field may be of functional relevance for proteins in complexes with partner proteins (Δ E_ads) and attached to membranes (Δ E_ads + Δ E_int).

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 (Δ E_Q = +2.2 mm·s^-1) and asymmetric parameter (η = 0.9) of the five-coordinate heme carbene [Fe(TTP)(CCl₂)], 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-CCl₂, 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 CCl₂, 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.

Ergodicity breaking of iron displacement in heme proteins

We present a model of the dynamical transition of atomic displacements in proteins. Increased mean-square displacement at higher temperatures is caused by the softening of the force constant for atomic/molecular displacements by electrostatic and van der Waals forces from the protein-water thermal bath. Displacement softening passes through a nonergodic dynamical transition when the relaxation time of the force-force correlation function enters, with increasing temperature, the instrumental observation window. Two crossover temperatures are identified. The lower crossover, presently connected to the glass transition, is related to the dynamical unfreezing of rotations of water molecules within nanodomains polarized by charged surface residues of the protein. The higher crossover temperature, usually assigned to the dynamical transition, marks the onset of water translations. All crossovers are ergodicity breaking transitions depending on the corresponding observation windows. Allowing stretched exponential relaxation of the protein-water thermal bath significantly improves the theory-experiment agreement when applied to solid protein samples studied by Mössbauer spectroscopy.

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 ⁵⁷Fe 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 (O₂, 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.

Revised CHARMM force field parameters for iron-containing cofactors of photosystem II

Photosystem II is a complex protein-cofactor machinery that splits water molecules into molecular oxygen, protons, and electrons. All-atom molecular dynamics simulations have the potential to contribute to our general understanding of how photosystem II works. To perform reliable all-atom simulations, we need accurate force field parameters for the cofactor molecules. We present here CHARMM bonded and non-bonded parameters for the iron-containing cofactors of photosystem II that include a six-coordinated heme moiety coordinated by two histidine groups, and a non-heme iron complex coordinated by bicarbonate and four histidines. The force field parameters presented here give water interaction energies and geometries in good agreement with the quantum mechanical target data. © 2017 Wiley Periodicals, Inc.

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 (N₂O). 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 ⁵⁷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.

New insights into nitric oxide reductases (NORs) are obtained. Using nuclear resonance vibrational spectroscopy, we probe both iron atoms in mononitrosylated intermediates of the NOR reaction in a biosynthetic protein model that reveal new insights into the structural and electronic features responsible for the NOR activity and its likely mechanism.

Investigations of the low frequency modes of ferric cytochrome c using vibrational coherence spectroscopy

Femtosecond vibrational coherence spectroscopy is used to investigate the low frequency vibrational dynamics of the electron transfer heme protein, cytochrome c (cyt c). The vibrational coherence spectra of ferric cyt c have been measured as a function of excitation wavelength within the Soret band. Vibrational coherence spectra obtained with excitation between 412 and 421 nm display a strong mode at ∼44 cm^–1 that has been assigned to have a significant contribution from heme ruffling motion in the electronic ground state. This assignment is based partially on the presence of a large heme ruffling distortion in the normal coordinate structural decomposition (NSD) analysis of the X-ray crystal structures. When the excitation wavelength is moved into the ∼421–435 nm region, the transient absorption increases along with the relative intensity of two modes near ∼55 and 30 cm^–1. The intensity of the mode near 44 cm^–1 appears to minimize in this region and then recover (but with an opposite phase compared to the blue excitation) when the laser is tuned to 443 nm. These observations are consistent with the superposition of both ground and excited state coherence in the 421–435 nm region due to the excitation of a weak porphyrin-to-iron charge transfer (CT) state, which has a lifetime long enough to observe vibrational coherence. The mode near 55 cm^–1 is suggested to arise from ruffling in a transient CT state that has a less ruffled heme due to its iron d⁶ configuration.

The dynamics of the non-heme iron in bacterial reaction centers from Rhodobacter sphaeroides

We investigate the dynamical properties of the non-heme iron (NHFe) in His-tagged photosynthetic bacterial reaction centers (RCs) isolated from Rhodobacter (Rb.) sphaeroides. Mössbauer spectroscopy and nuclear inelastic scattering of synchrotron radiation (NIS) were applied to monitor the arrangement and flexibility of the NHFe binding site. In His-tagged RCs, NHFe was stabilized only in a high spin ferrous state. Its hyperfine parameters (IS=1.06±0.01mm/s and QS=2.12±0.01mm/s), and Debye temperature (θ(D0)~167K) are comparable to those detected for the high spin state of NHFe in non-His-tagged RCs. For the first time, pure vibrational modes characteristic of NHFe in a high spin ferrous state are revealed. The vibrational density of states (DOS) shows some maxima between 22 and 33meV, 33 and 42meV, and 53 and 60meV and a very sharp one at 44.5meV. In addition, we observe a large contribution of vibrational modes at low energies. This iron atom is directly connected to the protein matrix via all its ligands, and it is therefore extremely sensitive to the collective motions of the RC protein core. A comparison of the DOS spectra of His-tagged and non-His-tagged RCs from Rb. sphaeroides shows that in the latter case the spectrum was overlapped by the vibrations of the heme iron of residual cytochrome c(2), and a low spin state of NHFe in addition to its high spin one. This enabled us to pin-point vibrations characteristic for the low spin state of NHFe.

Adsorption of sulfite oxidase on self-assembled monolayers from molecular dynamics simulations

Sulfite oxidase (SO) is an enzyme catalyzing the terminal step of the metabolism of sulfur-containing amino acids that is essential for almost all living organisms. The catalytic activity of SO in vertebrates strongly depends on the efficiency of the intramolecular electron transfer (IET) between the catalytic Moco domain and the cytochrome b5 (cyt b5) domain. The IET process is assumed to be mediated by large domain motions of the cyt b5 domains within the enzyme. Thus, the interaction of SO with charged surfaces may affect the mobility of the cyt b5 domain required for IET and consequently hinder SO activation. In this study, we present a molecular dynamics approach to investigating the ionic strength dependence of the initial surface adsorption of SO in two different conformations-the crystallographic structure and the model structure for an activated SO-onto mixed amino- and hydroxyl-terminated SAMs. The results show for both conformations at low ionic strengths a strong adsorption of the cyt b5 units onto the SAM, which inhibits the domain motion event required for IET. Under higher ion concentrations, however, the interaction with the surface is weakened by the negatively charged ions acting as a buffer and competing in adsorption with the cathodic cyt b5 domains. This competition prevents the immobilization of the cytochrome b5 units onto the surface, allowing the intramolecular domain motions favoring IET. Our predictions support the interpretation of recent experimental spectroelectrochemical studies on SO.

Nuclear inelastic scattering and Mössbauer spectroscopy as local probes for ligand binding modes and electronic properties in proteins: vibrational behavior of a ferriheme center inside a β-barrel protein

In this work, we present a study of the influence of the protein matrix on its ability to tune the binding of small ligands such as NO, cyanide (CN(-)), and histamine to the ferric heme iron center in the NO-storage and -transport protein Nitrophorin 2 (NP2) from the salivary glands of the blood-sucking insect Rhodnius prolixus. Conventional Mössbauer spectroscopy shows a diamagnetic ground state of the NP2-NO complex and Type I and II electronic ground states of the NP2-CN(-) and NP2-histamine complex, respectively. The change in the vibrational signature of the protein upon ligand binding has been monitored by Nuclear Inelastic Scattering (NIS), also called Nuclear Resonant Vibrational Spectroscopy (NRVS). The NIS data thus obtained have also been calculated by quantum mechanical (QM) density functional theory (DFT) coupled with molecular mechanics (MM) methods. The calculations presented here show that the heme ruffling in NP2 is a consequence of the interaction with the protein matrix. Structure optimizations of the heme and its ligands with DFT retain the characteristic saddling and ruffling only if the protein matrix is taken into account. Furthermore, simulations of the NIS data by QM/MM calculations suggest that the pH dependence of the binding of NO, but not of CN(-) and histamine, might be a consequence of the protonation state of the heme carboxyls.

Electrostatics of the protein-water interface and the dynamical transition in proteins

Atomic displacements of hydrated proteins are dominated by phonon vibrations at low temperatures and by dissipative large-amplitude motions at high temperatures. A crossover between the two regimes is known as a dynamical transition. Recent experiments indicate a connection between the dynamical transition and the dielectric response of the hydrated protein. We analyze two mechanisms of the coupling between the protein atomic motions and the protein-water interface. The first mechanism considers viscoelastic changes in the global shape of the protein plasticized by its coupling to the hydration shell. The second mechanism involves modulations of the local motions of partial charges inside the protein by electrostatic fluctuations. The model is used to analyze mean-square displacements of iron of metmyoglobin reported by Mössbauer spectroscopy. We show that high displacement of heme iron at physiological temperatures is dominated by electrostatic fluctuations. Two onsets, one arising from the viscoelastic response and the second from electrostatic fluctuations, are seen in the temperature dependence of the mean-square displacements when the corresponding relaxation times enter the instrumental resolution window.

Evidence of protein collective motions on the picosecond timescale

We investigate the presence of structural collective motions on a picosecond timescale for the heme protein, cytochrome c, as a function of oxidation and hydration, using terahertz (THz) time domain spectroscopy and molecular dynamics simulations. The THz response dramatically increases with oxidation, with the largest increase for lowest hydrations, and highest frequencies. For both oxidation states the THz response rapidly increases with hydration saturating above ∼25% (g H(2)O/g protein). Quasiharmonic vibrational modes and dipole-dipole correlation functions were calculated from molecular dynamics trajectories. The collective mode density of states alone reproduces the measured hydration dependence, providing strong evidence of the existence of these motions. The large oxidation dependence is reproduced only by the dipole-dipole correlation function, indicating the contrast arises from diffusive motions consistent with structural changes occurring in the vicinity of buried internal water molecules. This source for the observed oxidation dependence is consistent with the lack of an oxidation dependence in nuclear resonant vibrational spectroscopy measurements.

Dynamical transition in proteins and non-Gaussian behavior of low-frequency modes in self-consistent normal mode analysis

Self-consistent normal mode analysis (SCNMA) is applied to heme c type cytochrome f to study temperature-dependent protein motion. Classical normal mode analysis assumes harmonic behavior and the protein mean-square displacement has a linear dependence on temperature. This is only consistent with low-temperature experimental results. To connect the protein vibrational motions between low and physiological temperatures, we have incorporated a fitted set of anharmonic potentials into SCNMA. In addition, quantum harmonic-oscillator theory has been used to calculate the displacement distribution for individual vibrational modes. We find that the modes involving soft bonds exhibit significant non-Gaussian dynamics at physiological temperature, which suggests that it may be the cause of the non-Gaussian behavior of the protein motions probed by elastic incoherent neutron scattering. The combined theory displays a dynamical transition caused by the softening of few "torsional" modes in the low-frequency regime ( <50 cm(-1) or <6 meV or >0.6 ps). These modes change from Gaussian to a classical distribution upon heating. Our theory provides an alternative way to understand the microscopic origin of the protein dynamical transition.

Vibrational dynamics of iron in cytochrome C

Nuclear resonance vibrational spectroscopy (NRVS) and Raman spectroscopy on (54)Fe- and (57)Fe-enriched cytochrome c (cyt c) identify multiple bands involving vibrations of the heme Fe. Comparison with predictions from Fe isotope shifts reveals that 70% of the NRVS signal in the 300-450 cm(-1) frequency range corresponds to vibrations resolved in Soret-enhanced Raman spectra. This frequency range dominates the "stiffness", an effective force constant determined by the Fe vibrational density of states (VDOS), which measures the strength of nearest-neighbor interactions with Fe. The stiffness of the low-spin Fe environment in both oxidation states of cyt c significantly exceeds that for the high-spin Fe in deoxymyoglobin, where the 200-300 cm(-1) frequency range dominates the VDOS. This situation is reflected in the shorter Fe-ligand bond lengths in the former with respect to the latter. The longer Fe-S(Met80) in oxidized cyt c with respect to reduced cyt c leads to a decrease in the stiffness of the iron environment upon oxidation. Comparison with NRVS measurements allows us to assess assignments for vibrational modes resolved in this region of the heme Raman spectrum. We consider the possibility that the 372 cm(-1) band in reduced cyt c involves the Fe-S(Met80) bond.