Sampling field heterogeneity at the heme of c-type cytochromes by spectral hole burning spectroscopy and electrostatic calculations

Cytochrome c: A Multifunctional Protein Combining Conformational Rigidity with Flexibility

Cytochrome has served as a model system for studying redox reactions, protein folding, and more recently peroxidase activity induced by partial unfolding on membranes. This review illuminates some important aspects of the research on this biomolecule. The first part summarizes the results of structural analyses of its active site. Owing to heme-protein interactions the heme group is subject to both in-plane and out-of-plane deformations. The unfolding of the protein as discussed in detail in the second part of this review can be induced by changes of pH and temperature and most prominently by the addition of denaturing agents. Both the kinetic and thermodynamic folding and unfolding involve intermediate states with regard to all unfolding conditions. If allowed to sit at alkaline pH (11.5) for a week, the protein does not return to its folding state when the solvent is switched back to neutral pH. It rather adopts a misfolded state that is prone to aggregation via domain swapping. On the surface of cardiolipin containing liposomes, the protein can adopt a variety of partially unfolded states. Apparently, ferricytochrome c can perform biological functions even if it is only partially folded.

Hole-burning spectroscopy as a probe of nano-environments and processes in biomolecules: a review

Hole-burning spectroscopy, a high-resolution spectroscopic technique, allows details of heterogeneous nano-environments in biological systems to be obtained from broad absorption bands. Recently, this technique has been applied to proteins, nucleic acids, cells, and substructures of water to probe the electrostatic conditions created by macromolecules and the surrounding solvent. Starting with the factors that obscure the homogeneous linewidth of a chromophore within an inhomogeneously broadened absorption or emission band, we describe properties and processes in biological systems that are reflected in the measured hole spectra. The technique also lends itself to the resolution of perturbation experiments, such as temperature cycling to elucidate energy landscape barriers, applied external electric fields (Stark effect) to measure net internal electric fields, and applied hydrostatic pressure to find the volume compressibility of proteins.

Using spectroscopic tools to probe porphyrin deformation and porphyrin-protein interactions

The reactivity and functionality of heme proteins are to a significant extent determined by the conformation of their functional heme groups and by the interaction of axial ligands with their protein environment. This review focuses on experimental methods and theoretical concepts for elucidating symmetry lowering perturbations of the heme induced by the protein environment of the heme pocket. First, we discuss a variety of methods which can be used to probe the electric field at the heme, including spectral hole burning as well as low temperature absorption and room temperature circular dichroism spectroscopy. Second, we show how heme deformations can be described as superposition of deformations along normal coordinates, thereby using the irreducible representations of the D4h point group as a classification tool. Finally, resonance Raman spectroscopy is introduced as a tool to probe the deformations of metalloprophyrins in solution and in protein matrices by measuring and comparing intensities and depolarization properties rather than wavenumber positions.

Internal electric field in cytochrome C explored by visible electronic circular dichroism spectroscopy

Electronic circular dichroism (ECD) is a valuable tool to explore the secondary and tertiary structure of proteins. With respect to heme proteins, the corresponding visible ECD spectra, which probe the chirality of the heme environment, have been used to explore functionally relevant structural changes in the heme vicinity. While the physical basis of the obtained ECD signal has been analyzed by Woody and co-workers in terms of multiple electronic coupling mechanism between the electronic transitions of the heme chromophore and of the protein (Hsu, M.C.; Woody, R.W. J. Am. Chem. Soc. 1971, 93, 3515), a theory for a detailed quantitative analysis of ECD profiles has only recently been developed (Schweitzer-Stenner, R.; Gorden, J. P.; Hagarman, A. J. Chem. Phys. 2007, 127, 135103). In the present study this theory is applied to analyze the visible ECD-spectra of both oxidation states of three cytochromes c from horse, cow and yeast. The results reveal that both B- and Q-bands are subject to band splitting, which is caused by a combination of electronic and vibronic perturbations. The B-band splittings are substantially larger than the corresponding Q-band splittings in both oxidation states. For the B-bands, the electronic contribution to the band splitting can be assigned to the internal electric field in the heme pocket, whereas the corresponding Q-band splitting is likely to reflect its gradient (Manas, E. S.; Vanderkooi, J. M.; Sharp, K. A. J. Phys. Chem. B 1999, 103, 6344). We found that the electronic and vibronic splitting is substantially larger in the oxidized than in the reduced state. Moreover, these states exhibit different signs of electronic splitting. These findings suggest that the oxidation process increases the internal electric field and changes its orientation with respect to the molecular coordinate system associated with the N-Fe-N lines of the heme group. For the reduced state, we used our data to calculate electric field strengths between 27 and 31 MV/cm for the investigated cytochrome c species. The field of the oxidized state is more difficult to estimate, owing to the lack of information about its orientation in the heme plane. Based on band splitting and the wavenumber of the band position we estimated a field-strength of ca. 40 MV/cm for oxidized horse heart cytochrome c. The thus derived difference between the field strengths of the oxidized and reduced state would contribute at least -55 kJ/mol to the enthalpic stabilization of the oxidized state. Our data indicate that the corresponding stabilization energy of yeast cytochrome c is smaller.

Conformational substates of horse heart cytochrome c exhibit different thermal unfolding of the heme cavity

The charge transfer (CT) band at 695 nm in the spectrum of ferri-cytochrome c is highly asymmetric, indicating conformational heterogeneity due to the coexistence of different conformational substates. We have measured the respective band profile of horse heart ferri-cytochrome c as a function of temperature between 283 K (10 degrees C) and 333 K (60 degrees C) and found that the well-known decrease of the absorptivity is wavenumber-dependent and exhibits a biphasic behavior. This indicates that the underlying conformational substates differ in their thermodynamic stability with respect to the structural changes associated with the disappearance of the 695 nm band, which eventually (at high temperatures) involves the replacement of M80 by a nearby lysine residue. Our data further indicates that the thermal unfolding process involves two structurally different intermediate states.

Static Normal Coordinate Deformations of the Heme Group in Mutants of Ferrocytochrome c from Saccharomyces cerevisiae Probed by Resonance Raman Spectroscopy

The function of heme proteins is, to a significant extent, influenced by the ligand field probed by the heme iron, which itself can be affected by deformations of the heme macrocycle. The exploration of this field is difficult because the heme structure obtained from X-ray crystallography is not resolved enough to unambiguously identify structural changes on the scale of 10(-2) A. However, asymmetric deformations in this order of magnitude affect the depolarization ratio of the resonance Raman lines assignable to normal vibrations of the heme group. We have measured the dispersion of the depolarization ratios of four structure sensitive Raman bands (i.e., nu4, nu11, nu21, and nu28) in yeast iso-1-ferrocytochrome c and its mutants N52V, Y67F, and N52VY67F with B- and Q-band excitation. The DPR dispersion of all bands indicates the presence of asymmetric in-plane and out-of-plane deformations. The replacement of the polar tyrosine residue at position 67 by phenylalanine significantly increases the triclinic B2g deformation, which involves a distortion of the pyrrole symmetry. We relate this deformation to changes of the electronic structure of pyrrole A, which modulates the interaction between its propionate substituents and the protein environment. This specific heme deformation is eliminated in the double mutant N52VY67F. The additional substitution of N52 by valine induces a tetragonal B1g deformation which involves asymmetric changes of the Fe-N distances and increases the rhombicity of the ligand field probed by the heme iron. This heme deformation might be caused by the elimination of the water-protein hydrogen-bonding network in the heme cavity. The single mutation N52V does not significantly perturb the heme symmetry, but a small B1g deformation is consistent with our data and the heme structure obtained from a 1 ns molecular dynamics simulation of the protein.

Experiments with proteins at low temperature: What do we learn on properties in their functional state?

The authors compared the spectral response of Zn-substituted horseradish peroxidase in a glycerol/water solvent to hydrostatic pressure at 2 K and ambient temperature. The low temperature experiments clearly demonstrate the presence of at least three different conformations with drastically different elastic properties. However, the main conformation, which determines the fluorescence spectrum at ambient temperature, did not show any significant difference between low and high temperature and pressure. The authors conclude that the local compressibility of the heme pocket of the protein depends only very weakly on temperature.

The importance of vibronic perturbations in ferrocytochrome c spectra: a reevaluation of spectral properties based on low-temperature optical absorption, resonance Raman, and molecular-dynamics simulations

We have measured and analyzed the low-temperature (T=10 K) absorption spectrum of reduced horse heart and yeast cytochrome c. Both spectra show split and asymmetric Q(0) and Q(upsilon) bands. The spectra were first decomposed into the individual split vibronic sidebands assignable to B(1g) (nu15) and A(2g) (nu19, nu21, and nu22) Herzberg-Teller active modes due to their strong intensity in resonance Raman spectra acquired with Q(0) and Q(upsilon) excitations. The measured band splittings and asymmetries cannot be rationalized solely in terms of electronic perturbations of the heme macrocycle. On the contrary, they clearly point to the importance of considering not only electronic perturbations but vibronic perturbations as well. The former are most likely due to the heterogeneity of the electric field produced by charged side chains in the protein environment, whereas the latter reflect a perturbation potential due to multiple heme-protein interactions, which deform the heme structure in the ground and excited states. Additional information about vibronic perturbations and the associated ground-state deformations are inferred from the depolarization ratios of resonance Raman bands. The results of our analysis indicate that the heme group in yeast cytochrome c is more nonplanar and more distorted along a B(2g) coordinate than in horse heart cytochrome c. This conclusion is supported by normal structural decomposition calculations performed on the heme extracted from molecular-dynamic simulations of the two investigated proteins. Interestingly, the latter are somewhat different from the respective deformations obtained from the x-ray structures.

Optical band splitting and electronic perturbations of the heme chromophore in cytochrome C at room temperature probed by visible electronic circular dichroism spectroscopy

We have measured the electronic circular dichroism (ECD) of the ferri- and ferro-states of several natural cytochrome c derivatives (horse heart, chicken, bovine, and yeast) and the Y67F mutant of yeast in the region between 300 and 750 nm. Thus, we recorded the ECD of the B- and Q-band region as well as the charge-transfer band at approximately 695 nm. The B-band region of the ferri-state displays a nearly symmetric couplet at the B0-position that overlaps with a couplet 790 cm-1 higher in energy, which we assigned to a vibronic side-band transition. For the ferro-state, the couplet is greatly reduced, but still detectable. The B-band region is dominated by a positive Cotton effect at energies lower than B0 that is attributed to a magnetically allowed iron-->heme charge-transfer transition as earlier observed for nitrosyl myoglobin and hemoglobin. The Q-band region of the ferri-state is poorly resolved, but displays a pronounced positive signal at higher wavenumbers. This must result from a magnetically allowed transition, possibly from the methionine ligand to the dxy-hole of Fe3+. For the ferro-state, the spectra resolve the vibronic structure of the Qv-band. A more detailed spectral analysis reveals that the positively biased spectrum can be understood as a superposition of asymmetric couplets of split Q0 and Qv-states. Substantial qualitative and quantitative differences between the respective B-state and Q-state ECD spectra of yeast and horse heart cytochrome c can clearly be attributed to the reduced band splitting in the former, which results from a less heterogeneous internal electric field. Finally, we investigated the charge-transfer band at 695 nm in the ferri-state spectrum and found that it is composed of at least three bands, which are assignable to different taxonomic substates. The respective subbands differ somewhat with respect to their Kuhn anisotropy ratio and their intensity ratios are different for horse and yeast cytochrome c. Our data therefore suggests different substate populations for these proteins, which is most likely assignable to a structural heterogeneity of the distal Fe-M80 coordination of the heme chromophore.

Functionally Relevant Electric-Field Induced Perturbations of the Prosthetic Group of Yeast Ferrocytochrome c Mutants Obtained from a Vibronic Analysis of Low-Temperature Absorption Spectra

We have measured the low temperature (T = 20 K) absorption spectra of the N52A, N52V, N52I, Y67F, and N52AY67F mutants of ferrous Saccharomyces cerevisiae (baker's yeast) cytochrome c. All the bands in the Q0- and Q(v)-band region are split, and the intensity distributions among the split bands are highly asymmetric. The spectra were analyzed by a decomposition into Voigtian profiles. The spectral parameters thus obtained were further analyzed in terms of the vibronic coupling model of Schweitzer-Stenner and Bigman (Schweitzer-Stenner, R.; Bigman, D. J. Phys. Chem. B 2001, 7064-7073) to identify parameters related to electronic and vibronic perturbations of the heme macrocycle. We report that the electronic perturbation is of B(1g) symmetry and reflects the heterogeneity of the electric field at the heme, that is, the difference between the gradients along the perpendicular N-Fe-N axis of the heme core. We found that all the investigated mutations substantially increase this electronic perturbation, so that the spectral properties become similar to those of horse heart cytochrome c. Moreover, the electronic perturbation was found to correlate nonlinearly with the enthalpy changes associated with the reduction of the heme iron. Group theoretical arguments are invoked to propose a simple model which explains how a perturbation of the obtained symmetry can stabilize the reduced state of the heme iron. Finally, vibronic coupling parameters obtained from the analysis of the Q(v)-band region suggest that the investigated mutations decrease the nonplanar deformations of the heme group. This finding was reproduced by a normal mode structural decomposition (NSD) analysis of the N52V and N52VY67F heme conformations obtained from a 1 ns molecular dynamics simulation. We argue that the reduced nonplanarity contributes to the stabilization of the reduced state.

Optical and IR absorption as probe of dynamics of heme proteins

The spectroscopy of horseradish peroxidase with and without the substrate analogue benzohydroxamic acid (BHA) was monitored in different solvents as a function of the temperature in the interval from 10 to 300 K. Thermal broadening of the Q(0,0) optical absorption band arises mainly from interaction of the electronic pi --> pi transition with the heme vibrations. In contrast, the width of the IR absorption band of CO bound to heme is controlled by the coupling of the CO transition moment to the electric field of the protein matrix. The IR bandwidth of the substrate free enzyme in the glycerol/H2O solvent hardly changes in the glassy matrix and strongly increases upon heating above the glass transition. Heating of the same enzyme in the trehalose/H2O glass considerably broadens the band. The binding of the substrate strongly diminishes the temperature broadening of the CO band. This result is consistent with the view that the BHA strongly reduces the amplitude of vibrations of the heme pocket environment. Unusually strong thermal broadening of the CO band above the glass transition is interpreted to be caused by thermal population of a very flexible excited conformational substate. The thermal broadening of the same band in the trehalose glass is caused by an increase of the protein vibrational amplitude in each of the conformational substates, their population being independent of the temperature in the glassy matrix.

Molten-globule and other conformational forms of zinc cytochrome C. Effect of partial and complete unfolding of the protein on its electron-transfer reactivity

To test the effect of protein conformation on reactivity, we use laser flash photolysis to compare the electron-transfer properties of the triplet state of zinc-substituted cytochrome c, designated (3)Zncyt, in the folded forms at low (F(low)) and high (F(high)) ionic strength, molten-globule (MG) form, and the forms unfolded by acid (U(acid)) and urea (U(urea)) toward the following four oxidative quenchers: Fe(CN)(6)(3-), Co(acac)(3), Co(phen)(3)(3+), and iron(III) cytochrome c. We characterize the conformational forms of Zncyt on the basis of the far-UV circular dichroism, Soret absorption, and rate constant for natural decay of the triplet state. This rate constant in the absence of quencher increases in the order F(high) < F(low) < MG < U(acid) < U(urea) because the exposure of porphyrin to solvent increases as Zncyt unfolds. Bimolecular rate constants for the reaction of (3)Zncyt with the four quenchers show significant effects on reactivity of electrostatic interactions and porphyrin exposure to solvent. This rate constant at the ionic strength of 20 mM increases upon unfolding by urea and acid, respectively, as follows: 1340-fold and 466-fold when the quencher is Co(phen)(3)(3+) and 168-fold and 36-fold when the quencher is cyt(III). To compare reactivity of (3)Zncyt in the F(low), F(high), MG, U(acid), and U(urea) forms without complicating effects of electrostatic interactions, we used the electroneutral quencher Co(acac)(3). Indeed, reactivity of folded (3)Zncyt with Co(acac)(3) was independent of ionic strength. Reactivity of (3)Zncyt with Co(acac)(3) upon partial and complete unfolding increases 10-fold, 54-fold, and 64-fold in the molten-globule, urea-unfolded, and acid-unfolded forms.

Influence of static and dynamic disorder on the visible and infrared absorption spectra of carbonmonoxy horseradish peroxidase

Spectroscopy of horseradish peroxidase with and without the substrate analog, benzohydroxamic acid, was monitored in a glycerol/water solvent as a function of temperature. It was determined from the water infrared (IR) absorption that the solvent has a glass transition at 170-180 K. In the absence of substrate, both the heme optical Q(0,0) absorption band and the IR absorption band of CO bound to heme broaden markedly upon heating from 10-300 K. The Q(0,0) band broadens smoothly in the whole temperature interval, whereas the IR bandwidth is constant in the glassy matrix and increases from 7 to 16 cm(-1) upon heating above the glass transition. Binding of substrate strongly diminishes temperature broadening of both the bands. The results are consistent with the view that the substrate strongly reduces the amplitude of motions of amino acids forming the heme pocket. The main contribution to the Q(0,0) bandwidth arises from the heme vibrations that are not affected by the phase transition. The CO band thermal broadening stems from the anharmonic coupling with motions of the heme environment, which, in the glassy state, are frozen in. Unusually strong temperature broadening of the CO band is interpreted to be caused by thermal population of a very flexible excited conformational substrate. Analysis of literature data on the thermal broadening of the A(0) band of Mb(CO) (Ansari et al., 1987. Biophys. Chem. 26:337-355) shows that such a state presents itself also in myoglobin.

Aromatic substrate specificity of horseradish peroxidase C studied by a combined fluorescence line narrowing/energy minimization approach: the effect of localized side-chain reorganization

Horseradish peroxidase C binds a wide variety of small H-donor compounds such as benzohydroxamic acid (BHA) and 2-naphthohydroxamic acid (NHA). In this work, we use the Mg(II)-mesoporphyrin prosthetic group derivative as a spectroscopic probe of the active site and of the interaction with the substrates. We report on high-resolution fluorescence line-narrowed spectra which show that the effects of substrate binding on the electronic transitions are similar for both substrates and present data on the normal vibrational modes that are active in the vibronic spectra. Analysis of the vibrational frequencies shows that the Mg(II) ion is 5-coordinate in all cases, thus ruling out a solvent water as sixth ligand. The frequency shifts observed as a result of substrate binding are also indicative of a more rigid prosthetic group upon substrate binding. We present models for MgMP-HRP and its complexes with both substrates and compare the resulting structures on the basis of a modeling approach combining energy minimization to finite difference Poisson--Boltzmann calculations which partitions the various relative protein contributions to substrate binding. We show that the electrostatic potential of the prosthetic group is modified by the binding event. Analysis of the unbound and bound energy-minimized structures shows that the enzyme modulates substrate binding by subtle charge reorganization in the vicinity of the catalytic site and that this rearrangement is not attributable to significant secondary structure conformational changes but to side-chain reorganization.

Macromolecular electrostatics: continuum models and their growing pains

Theoretical understanding of macromolecular electrostatics has advanced substantially over the past year. Continuum models have given promising results for calculating protein-ligand binding free energy differences, as well as pK(a)s and redox properties, particularly with explicit treatment of multiple conformers. Generalized Born and other techniques have led to the first molecular dynamics simulations of proteins and RNA with continuum solvent. Continuum and microscopic descriptions of dielectric relaxation have been critically compared.