Allosteric linkage-induced distortions of the prosthetic group in haem proteins as derived by the theoretical interpretation of the depolarization ratio in resonance Raman scattering

Probing the versatility of cytochrome c by spectroscopic means: A Laudatio on resonance Raman spectroscopy

Over the last 50 years resonance Raman spectroscopy has become an invaluable tool for the exploration of chromophores in biological macromolecules. Among them, heme proteins and metal complexes have attracted considerable attention. This interest results from the fact that resonance Raman spectroscopy probes the vibrational dynamics of these chromophores without direct interference from the surrounding. However, the indirect influence via through-bond and through-space chromophore-protein interactions can be conveniently probed and analyzed. This review article illustrates this point by focusing on class 1 cytochrome c, a comparatively simple heme protein generally known as electron carrier in mitochondria. The article demonstrates how through selective excitation of resonance Raman active modes information about the ligation, the redox state and the spin state of the heme iron can be obtained from band positions in the Raman spectra. The investigation of intensities and depolarization ratios emerged as tools for the analysis of in-plane and out-of-plane deformations of the heme macrocycle. The article further shows how resonance Raman spectroscopy was used to characterize partially unfolded states of oxidized cytochrome c. Finally, it describes its use for exploring structural changes due to the protein's binding to anionic surfaces like cardiolipin containing membranes.

Heme-Protein Interactions and Functional Relevant Heme Deformations: The Cytochrome c Case

Heme proteins are known to perform a plethora of biologically important functions. This article reviews work that has been conducted on various class I cytochrome c proteins over a period of nearly 50 years. The article focuses on the relevance of symmetry-lowering heme-protein interactions that affect the function of the electron transfer protein cytochrome c. The article provides an overview of various, mostly spectroscopic studies that explored the electronic structure of the heme group in these proteins and how it is affected by symmetry-lowering deformations. In addition to discussing a large variety of spectroscopic studies, the article provides a theoretical framework that should enable a comprehensive understanding of the physical chemistry that underlies the function not only of cytochrome c but of all heme proteins.

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.

Polarized resonance Raman dispersion spectroscopy on metalporphyrins

Resonance Raman spectroscopy is an ideal tool to investigate the structural properties of chromophores embedded in complex (biological) environments. This holds particularly for metalporphyrins which serve as prosthetic group in various proteins. Resonance Raman dispersion spectroscopy involves the measurement of resonance excitation and depolarization ratios of a large number of Raman lines at various excitation energies covering the spectral region of the chromophore's optical absorption bands. Thus, one obtains resonance excitation profiles and the depolarization ratio dispersion of these bands. While the former contains information about the structure of excited electronic states involved in the Raman scattering process, the latter reflects asymmetric perturbations which lower the porphyrin macrocycle symmetry from ideal D4h. The article introduces and compares different quantum mechanical approaches designed to quantitatively analyze both resonance excitation and the relationship between symmetry lowering and depolarization ratio dispersion.

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.

Optical properties of Yeast Cytochrome c monolayer on gold: an in situ spectroscopic ellipsometry investigation

The adsorption of Yeast Cytochrome c (YCC) on well defined, flat gold substrates has been studied by Spectroscopic Ellipsometry (SE) in the 245-1000 nm wavelength range. The investigation has been performed in aqueous ambient at room temperature, focusing on monolayer-thick films. In situ δΨ and δΔ difference spectra have shown reproducibly well-defined features related to molecular optical absorptions typical of the so-called heme group. The data have been reproduced quantitatively by a simple isotropic optical model, accounting for the molecular absorption spectrum and film-substrate interface effects. The simulations allowed a reliable estimate of the film thickness and the determination of the position and the shape of the so-called Soret absorption peak that, within the experimental uncertainty, is the same found for molecules in liquid. These findings suggest that YCC preserves its native structure upon adsorption. The same optical model was able to reproduce also ex situ results on rinsed and dried samples, dominated by the spectral features associated to the polypeptide chain that tend to overwhelm the heme absorption features.

Raman dispersion spectroscopy probes heme distortions in deoxyHb-trout IV involved in its T-state Bohr effect

The depolarization ratios of heme protein Raman lines arising from vibrations of the heme group exhibit significant dependence on the excitation wavelength. From the analysis of this depolarization ratio dispersion, one obtains information about symmetry-lowering distortions deltaQ(Gamma) of the heme group that can be classified in terms of the symmetry races Gamma = A(1g), B(1g), B(2g), and A(2g) in D(4h) symmetry. The heme-protein interaction can be changed by the protonation of distinct amino acid side chains (i.e., for instance the Bohr groups in hemoglobin derivates), which gives rise to specific static heme distortions for each protonation state. From the Raman dispersion data, it is possible to obtain parameters by fitting to a theoretical expression of the Raman tensor, which provide information on these static distortions and also about the pK values of the involved titrable side chains. We have applied this method to the nu(4) (1,355 cm(-1)) and nu(10) (1,620 cm(-1)) lines of deoxygenated hemoglobin of the fourth component of trout and have measured their depolarization ratio dispersion as a function of pH between 6 and 9. From the pH dependence of the thus derived parameters, we obtain pK values identical to those of the Bohr groups, which were earlier derived from the corresponding O(2)-binding isotherms. These are pK(alpha1) = pK(alpha2) = 8.5 for the alpha and pK(beta1) = 7.5, pK(beta2) = 7.4 for the beta chains. We also obtain the specific distortion parameters for each protonation state. As shown in earlier studies, the nu(4) mode mainly probes distortions from interactions between the proximal histidine and atoms of the heme core (i.e., the nitrogens and the C(alpha) atoms of the pyrroles). Group theoretical argumentation allows us to relate specific changes of the imidazole geometry as determined by its tilt and azimuthal angle and the iron-out-of-plane displacement to distinct variations of the normal distortions deltaQ(Gamma) derived from the Raman dispersion data. Thus, we found that the pH dependence of the heme distortions deltaQ(A1g) (totally symmetric) and deltaQ(B1g) (asymmetric) is caused by variations of the azimuthal rather than the tilt angle of the Fe-His (F8) bond. In contrast to this, the nu(10) line mainly monitors changes resulting from the interaction between peripheral substituents of the porphyrin macrocycle (vinyl). From the pH dependence of the parameters, it is possible to separately identify distortions deltaQ(Gamma) affecting the hemes in the alpha and beta chains, respectively. From this, we find that in the alpha subunit structural changes induced on protonation of the corresponding Bohr groups are mainly transferred via the Fe-N(epsilon) bond and give rise to changes in the azimuthal angle. In the beta subunit, however, in addition, structural changes of the heme pocket arise, which most probably result from protonation of the imidazole of the COOH-terminal His (HC3 beta). This rearranges the net of H bonds between His HC3 beta, Ser (F9 beta), and Glu (F7 beta).

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.

Asymmetric band profile of the Soret band of deoxymyoglobin is caused by electronic and vibronic perturbations of the heme group rather than by a doming deformation

We measured the Soret band of deoxymyoglobin (deoxyMb), myoglobin cyanide (MbCN), and aquo-metmyoglobin (all from horse heart) with absorption and circular dichroism (CD) spectroscopies. A clear non-coincidence was observed between the absorption and CD profiles of deoxyMb and MbCN, with the CD profiles red- and blueshifted with respect to the absorption band position, respectively. On the contrary, the CD and absorption profiles of aquametMb were nearly identical. The observed noncoincidence indicates a splitting of the excited B state due to heme-protein interactions. CD and absorption profiles of deoxyMb and MbCN were self-consistently analyzed by employing a perturbation approach for weak vibronic coupling as well as the relative intensities and depolarization ratios of seven bands in the respective resonance Raman spectra measured with B-band excitation. The respective B(y) component was found to dominate the observed Cotton effect of both myoglobin derivatives. The different signs of the noncoincidences between CD and absorption bands observed for deoxyMb and MbCN are due to different signs of the respective matrix elements of A(1g) electronic interstate coupling, which reflects an imbalance of Gouterman's 50:50 states. The splitting of the B band reflects contributions from electronic and vibronic perturbations of B(1g) symmetry. The results of our analysis suggest that the broad and asymmetric absorption band of deoxyMb results from this band splitting rather than from its dependence on heme doming. Thus, we are able to explain recent findings that the temperature dependences of CO rebinding to myoglobin and the Soret band profile are uncorrelated[Ormos et al., Proc. Natl. Acad. Sci U.S.A. 95, 6762 (1998)].

Nonplanar heme deformations and excited state displacements in nickel porphyrins detected by Raman spectroscopy at soret excitation

We have correlated the Raman intensities of out-of-plane modes of nickel porphyrins with the nonplanar deformations of specific symmetries, i.e., static normal coordinate deformations (SNCDs) expressed in terms of irreducible representations of the unperturbed D(4h) point group. The model porphyrins Ni(II) octaethyltetraphenylporphyrin (NiOETPP), Ni(II) tetra(isopropyl)porphyrin (NiT((i)Pr)P), Ni(II) tetra(tert-butyl)porphyrin (NiT((t)Bu)P), and Ni(II) meso-tetraphenylporphyrin (NiTPP) were chosen because they exhibit significant out-of-plane deformations of different symmetries. At B-band excitation, the Raman scattering of out-of-plane modes becomes activated mostly via the Franck-Condon mechanism. Some characteristic bands from out-of-plane modes in the spectra were identified as reliable predictors of the type and magnitude of out-of-plane deformation. The gamma(10)-gamma(13) bands are indicators of ruffling (B(1u)) deformations for porphyrins, as confirmed by data for NiTPP, NiT((i)Pr)P, and NiT((t)Bu)P, where the Raman intensity increases with the magnitude of the ruffling deformation. The gamma(15)-gamma(17) bands are indicators of saddling (B(2u)) deformations, as shown by data for NiOETPP, which is highly saddled. By comparing the relative intensities of these prominent Raman bands we estimated the vibronic coupling parameters using a self-consistent analysis, and showed that they reproduce the respective B-band absorption profiles. We also identified the deformations along the lowest wavenumber normal coordinates as the predominant reason for the Raman activity of out-of-plane modes. Our results suggest that some of the normal coordinates (gamma(10) and gamma(13)) may be used as tools to quantitatively probe the nonplanar deformations of metalloporphyrins with alkyl substituents at the meso-positions. Out-of-plane deformations also increase the vibronic coupling strength of some low frequency in-plane Raman modes, namely, nu(7) and nu(8). Generally, the Raman data suggest that the excited B-state is substantially more nonplanar than the ground state. The overall larger vibronic coupling of ruffled porphyrins yields substantially larger dipole strengths for the vibronic sidebands associated with the B-state transition, so that the relative absorptivity of the B(v) band can be used as a convenient tool to probe the nonplanarity of the porphyrin macrocycle.

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.

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.

Symmetry Properties of Vibrational Modes in Mesoporphyrin IX Dimethyl Ester Investigated by Polarization-Sensitive Resonance Raman and CARS Spectroscopy

The symmetry properties of selected vibrational modes of mesoporphyrin IX dimethyl ester (MP-IX-DME) in solution are investigated under different electronic resonance conditions. The Raman band parameters of the macrocycle modes nu(2), nu(10), nu(11), and nu(19) are determined from a quantitative analysis of polarized spontaneous resonance Raman (RR) and polarization-sensitive (PS) multiplex coherent anti-Stokes Raman scattering (CARS) spectra obtained with pre-resonant B band and resonant Qx band excitation, respectively. Additionally, the molecular geometry and the vibrational modes of MP-IX-DME are calculated by employing density functional theory (DFT) on the B3LYP/6-31G(d) level. Both the DFT-derived structure and the Raman spectroscopic parameters of MP-IX-DME indicate minor deviations from an ideal D2h macrocycle symmetry. To assess the influence of the beta substitution pattern on the in-plane symmetry, calculated normal-mode vectors and several experimentally detected parameters, such as peak positions, depolarization ratios, and coherent phases, are analyzed. The effects of the macrocycle substitution pattern are different for the selected vibrational modes: nu(2) in particular is very sensitive to subtle perturbations of the in-plane symmetry. The considerable activity of totally symmetric vibrations observed in the PS CARS spectra of MP-IX-DME and the correlation of mode symmetries with coherent phases confirm earlier PS CARS results on octaethylporphine (OEP) acquired under the same electronic resonance conditions.

Proximal effects in the modulation of nitric oxide synthase reactivity: a QM-MM study

Nitric oxide synthases (NOS) are heme proteins that have a cysteine residue as axial ligand, which generates nitric oxide (NO). The proximal environment, specifically H-bonding between tryptophan (Trp) 178 and thiolate, has been proposed to play a fundamental role in the modulation of NOS activity. We analyzed the molecular basis of this modulation by performing electronic structure calculations on isolated model systems and hybrid quantum-classical computations of the active sites in the protein environment for wild-type and mutant (Trp 178 x Gly) proteins. Our results show that in the ferrous proteins NO exhibits a considerable trans effect. We also showed that in the ferrous (Fe(+2)) mutant NOS the absence of Trp, experimentally associated to a protonated cysteine, weakens the Fe-S bond and yields five coordinate complexes. In the ferric (Fe(+3)) state, the NO dissociation energy is shown to be slightly smaller in the mutant NOS, implying that the Fe(+3)-NO complex has a shorter half-life. We found computational evidence suggesting that ferrous NOS is favored in wild-type NOS when compared to the Trp mutant, consistently with the fact that Trp mutants have been shown to accumulate less Fe(+2)-NO dead end species. We also found that the heme macrocycle showed a significant distortion in the wild-type protein, due to the presence of the nearby Trp 178. This may also play a role in the subtle tuning of the electronic structure of the heme moiety.

Heme structural perturbation of PEG-modified horseradish peroxidase C in aromatic organic solvents probed by optical absorption and resonance Raman dispersion spectroscopy

The heme structure perturbation of poly(ethylene glycol)-modified horseradish peroxidase (HRP-PEG) dissolved in benzene and toluene has been probed by resonance Raman dispersion spectroscopy. Analysis of the depolarization ratio dispersion of several Raman bands revealed an increase of rhombic B(1g) distortion with respect to native HRP in water. This finding strongly supports the notion that a solvent molecule has moved into the heme pocket where it stays in close proximity to one of the heme's pyrrole rings. The interactions between the solvent molecule, the heme, and the heme cavity slightly stabilize the hexacoordinate high spin state without eliminating the pentacoordinate quantum mixed spin state that is dominant in the resting enzyme. On the contrary, the model substrate benzohydroxamic acid strongly favors the hexacoordinate quantum mixed spin state and induces a B(2g)-type distortion owing to its position close to one of the heme methine bridges. These results strongly suggest that substrate binding must have an influence on the heme geometry of HRP and that the heme structure of the enzyme-substrate complex (as opposed to the resting state) must be the key to understanding the chemical reactivity of HRP.

The endogenous calcium ions of horseradish peroxidase C are required to maintain the functional nonplanarity of the heme

Horseradish peroxidase C (HRPC) binds 2 mol calcium per mol of enzyme with binding sites located distal and proximal to the heme group. The effect of calcium depletion on the conformation of the heme was investigated by combining polarized resonance Raman dispersion spectroscopy with normal coordinate structural decomposition analysis of the hemes extracted from models of Ca(2+)-bound and Ca(2+)-depleted HRPC generated and equilibrated using molecular dynamics simulations. Results show that calcium removal causes reorientation of heme pocket residues. We propose that these rearrangements significantly affect both the in-plane and out-of-plane deformations of the heme. Analysis of the experimental depolarization ratios are clearly consistent with increased B(1g)- and B(2g)-type distortions in the Ca(2+)-depleted species while the normal coordinate structural decomposition results are indicative of increased planarity for the heme of Ca(2+)-depleted HRPC and of significant changes in the relative contributions of three of the six lowest frequency deformations. Most noteworthy is the decrease of the strong saddling deformation that is typical of all peroxidases, and an increase in ruffling. Our results confirm previous work proposing that calcium is required to maintain the structural integrity of the heme in that we show that the preferred geometry for catalysis is lost upon calcium depletion.

Structure of poly(ethylene glycol)-modified horseradish peroxidase in organic solvents: infrared amide I spectral changes upon protein dehydration are largely caused by protein structural changes and not by water removal per se

Fourier transform infrared (FTIR) spectroscopy has emerged as a powerful tool to guide the development of stable lyophilized protein formulations by providing information on the structure of proteins in amorphous solids. The underlying assumption is that IR spectral changes in the amide I and III region upon protein dehydration are caused by protein structural changes. However, it has been claimed that amide I IR spectral changes could be the result of water removal per se. Here, we investigated whether such claims hold true. The structure of horseradish peroxidase (HRP) and poly(ethylene glycol)-modified HRP (HRP-PEG) has been investigated under various conditions (in aqueous solution, the amorphous dehydrated state, and dissolved/suspended in toluene and benzene) by UV-visible (UV-Vis), FTIR, and resonance Raman spectroscopy. The resonance Raman and UV-Vis spectra of dehydrated HRP-PEG dissolved in neat toluene or benzene were very similar to that of HRP in aqueous buffer, and thus the heme environment (heme iron spin, coordination, and redox state) was essentially the same under both conditions. Therefore, the three-dimensional structure of HRP-PEG dissolved in benzene and toluene was similar to that in aqueous solution. The amide I IR spectra of HRP-PEG in aqueous buffer and of dehydrated HRP-PEG dissolved in neat benzene and toluene were also very similar, and the secondary structure compositions (percentages of alpha-helices and beta-sheets) were within the standard error the same. These results are irreconcilable with recent claims that water removal per se could cause substantial amide I IR spectral changes (M. van de Weert, P.I. Haris, W.E. Hennink, and D.J. Crommelin. 2001. Anal. Biochem. 297:160-169). On the contrary, amide I IR spectral changes upon protein dehydration are caused by perturbations in the secondary structure.

Protein dynamics in an intermediate state of myoglobin: optical absorption, resonance Raman spectroscopy, and x-ray structure analysis

A metastable state of myoglobin is produced by reduction of metmyoglobin at low temperatures. This is done either by irradiation with x-rays at 80 K or by electron transfer from photoexcited tris(2, 2'-bipyridine)-ruthenium(II) at 20 K. At temperatures above 150 K, the conformational transition toward the equilibrium deoxymyoglobin is observed. X-ray crystallography, Raman spectroscopy, and temperature-dependent optical absorption spectroscopy show that the metastable state has a six-ligated iron low-spin center. The x-ray structure at 115K proves the similarity of the metastable state with metmyoglobin. The Raman spectra yield the high-frequency vibronic modes and give additional information about the distortion of the heme. Analysis of the temperature dependence of the line shape of the Soret band reveals that a relaxation within the metastable state starts at approximately 120 K. Parameters representative of static properties of the intermediate state are close to those of CO-ligated myoglobin, while parameters representative of dynamics are close to deoxymyoglobin. Thus within the metastable state the relaxation to the equilibrium is initiated by changes in the dynamic properties of the active site.

Heme symmetry, vibronic structure, and dynamics in heme proteins: ferrous nicotinate horse myoglobin and soybean leghemoglobin

We report the visible and Soret absorption bands, down to cryogenic temperatures, of the ferrous nicotinate adducts of native and deuteroheme reconstituted horse heart myoglobin in comparison with soybean leghemoglobin-a. The band profile in the visible region is analyzed in terms of vibronic coupling of the heme normal modes to the electronic transition in the framework of the Herzberg-Teller approximation. This theoretical approach makes use of the crude Born-Oppenheimer states and therefore neglects the mixing between electronic and vibrational coordinates; however, it takes into account the vibronic nature of the visible absorption bands and allows an estimate of the vibronic side bands for both Condon and non-Condon vibrational modes. In this framework, an x-y splitting of the Q transition for native and deuteroheme reconstituted horse myoglobin is clearly assessed and attributed to electronic perturbations that, in turn, are caused by a reduction of the typical D(4h) symmetry of the system due to heme distortions of B(1g)-type symmetry and/or to an x-y asymmetric position of the nicotinate ring; in deuteroheme reconstituted horse myoglobin the asymmetric heme peripheral substituents add to the above effect(s). On the contrary, in leghemoglobin-a no spectral splitting upon nicotinate binding is observed, pointing to a planar heme configuration in which only distortions of A(1g)-type symmetry are effective and to which the nicotinate ring is bound in an x - y symmetric position. The local dynamic properties of the heme pocket of the three proteins are investigated through the temperature dependence of spectral line broadening. Leghemoglobin-a behaves as a softer matrix with respect to horse myoglobin, thus validating the hypothesis of a looser heme pocket conformation in the former protein, which allows a nondistorted heme configuration and a symmetric binding of the bulky nicotinate ligand.

Thermal fluctuations between conformational substates of the Fe(2+)-HisF8 linkage in deoxymyoglobin probed by the Raman active Fe-N epsilon (HisF8) stretching vibration

We have measured the VFe-His Raman band of horse heart deoxymyoglobin dissolved in an aqueous solution as a function of temperature between 10 and 300 K. The minimal model to which these data can be fitted in a statistically significant and physically meaningful way comprises four different Lorentzian bands with frequencies at 197, 209, 218, and 226 cm-1, and a Gaussian band at 240 cm-1, which exhibit halfwidths between 10 and 12.5 cm-1. All these parameters were assumed to be independent of temperature. The temperature dependence of the apparent total band shape's frequency is attributed to an intensity redistribution of the subbands at omega 1 = 209 cm-1, omega 2 = 218 cm-1, and omega 3 = 226 cm-1, which are assigned to Fe-N epsilon (HisF8) stretching modes in different conformational substrates of the Fe-HisF8 linkage. They comprise different out-of-plane displacements of the heme iron. The two remaining bands at 197 and 240 cm-1 result from porphyrin modes. Their intensity ratio is nearly temperature independent. The intensity ratio I3/I2 of the vFe-His subbands exhibits a van't Hoff behavior between 150 and 300 K, bending over in a region between 150 and 80 K, and remains constant between 80 and 10 K, whereas I2/I1 shows a maximum at 170 K and approaches a constant value at 80 K. These data can be fitted by a modified van't Hoff expression, which accounts for the freezing into a non-equilibrium distribution of substates below a distinct temperature Tf and also for the linear temperature dependence of the specific heat of proteins. The latter leads to a temperature dependence of the entropic and enthalpic differences between conformational substates. The fits to the intensity ratios of the vFe-His subbands yield a freezing temperature of Tf = 117 K and a transition region of delta T = 55 K. In comparison we have utilized the above thermodynamic model to reanalyze earlier data on the temperature dependence of the ratio Ao/A1 of two subbands underlying the infrared absorption band of the CO stretching vibration in CO-ligated myoglobin (A. Ansari, J. Berendzen, D. Braunstein, B. R. Cowen, H. Frauenfelder, M. K. Kong, I. E. T. Iben, J. Johnson, P. Ormos, T. B. Sauke, R. Scholl, A. Schulte, P. J. Steinbach, R. D. Vittitow, and R. D. Young, 1987, Biophys. Chem. 26:237-335). This yields thermodynamic parameters, in particular the freezing temperature (Tf = 231 K) and the width of the transition region (AT =8 K), which are significantly different from the corresponding parameters obtained from the above vFe-His data, but very close to values describing the transition of protein bound water from a liquid into an amorphous state. These findings and earlier reported data on the temperature dependence exhibited by the Soret absorption bands of various deoxy and carbonmonoxymyoglobins led us to the conclusion that the fluctuations between conformational substates of the heme environment in carbonmonoxymyoglobin are strongly coupled to motions within the hydration shell, whereas the thermal motions between the substates of the Fe-HisF8 linkage in deoxymyoglobin proceed on an energy landscape that is mainly determined by the intrinsic properties of the protein. The latter differ from protein fluctuations monitored by Mossbauer experiments ondeoxymyoglobin crystals which exhibit a strong coupling to the protein bound water and most probably reflect a higher tier in the hierarchical arrangement of substates and equilibrium fluctuations.

pH-induced conformational changes of the Fe(2+)-N epsilon (His F8) linkage in deoxyhemoglobin trout IV detected by the Raman active Fe(2+)-N epsilon (His F8) stretching mode

To investigate heme-protein coupling via the Fe(2+)-N epsilon (His F8) linkage we have measured the profile of the Raman band due to the Fe(2+)-N epsilon (His F8) stretching mode (nu Fe-His) of deoxyHb-trout IV and deoxyHbA at various pH between 6.0 and 9.0. Our data establish that the band of this mode is composed of five different sublines. In deoxyHb-trout IV, three of these sublines were assigned to distinct conformations of the alpha-subunit (omega alpha 1 = 202 cm-1, omega alpha 2 = 211 cm-1, omega alpha 3 = 217 cm-1) and the other two to distinct conformations of the beta-subunit (omega beta 1 = 223 cm-1 and omega beta 2 = 228 cm-1). Human deoxyHbA exhibits two alpha-chain sublines at omega alpha 1 = 203 cm-1, omega alpha 2 = 212 cm-1 and two beta-chain sublines at omega beta 1 = 217 cm-1 and omega beta 2 = 225 cm-1. These results reveal that each subunit exists in different conformations. The intensities of the nu Fe-His sublines in deoxyHb-trout IV exhibit a significant pH dependence, whereas the intensities of the corresponding sublines in the deoxyHbA spectrum are independent on pH. This finding suggests that the structural basis of the Bohr effect is different in deoxyHbA and deoxyHb-trout IV. To analyse the pH dependence of the deoxyHb-trout IV sublines we have applied a titration model describing the intensity of each nu Fe-His subline as an incoherent superposition of the intensities from sub-sublines with the same frequency but differing intrinsic intensities due to the different protonation states of the respective subunit. The molar fractions of these protonation states are determined by the corresponding Bohr groups (i.e., pK alpha 1 = pK alpha 2 = 8.5, pK beta 1 = 7.5, pK beta 2 = 7.4) and pH. Hence, the intensities of these sublines reflect the pH dependence of the molar fractions of the involved protonation states. Fitting this model to the pH-dependent line intensities yields a good reproduction of the experimental data. To elucidate the structural basis of the observed results we have employed models proposed by Bangchoroenpaurpong, O., K. T. Schomaker, and P. M. Champion. (1984. J. Am. Chem. Soc. 106:5688-5698) and Friedman, J. M., B. F. Campbell, and R. W. Noble. (1990. Biophys. Chem. 37:43-59) which describe the coupling between the sigma *orbitals of the Fe2+-NJ(His F8) bond and the phi * orbitals of the pyrrole nitrogens in terms of the tilt angle theta between its Fe2+-N,(HisF8)-bond and the heme normal and the azimuthal angle phi between the Fe2+-N.(His F8) projection on the heme and the N1-N3 axis.Our results indicate that each subconformation reflected by different frequencies of the VFe His-subline is related to different tilt angles theta, whereas the pH-induced intensity variations of each VFe His subline of the deoxy Hb trout IV spectrum are caused by changes of the azimuthal angle phi.