Raman dispersion spectroscopy on the highly saddled nickel(II)-octaethyltetraphenylporphyrin reveals the symmetry of nonplanar distortions and the vibronic coupling strength of normal modes

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.

Molecular simulations of porphyrins and heme proteins

An overview of the use of classical mechanical molecular simulations of porphyrins, hydroporphyrins and heme proteins is given. The topics cover molecular mechanics calculations of structures and conformer energies of porphyrins, energies of barriers for interconversion between stable conformers, molecular dynamics of porphyrins and heme proteins, and normal-coordinate structural analysis of experimental and calculated porphyrin structures. Molecular mechanics and dynamics are currently a fertile area of research on porphyrins. In the future, other computational methods such as Monte Carlo simulations, which have yet to be applied to porphyrins, will come into use and open new avenues of research into molecular simulations of porphyrins.

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.

Normal-coordinate structural decomposition and the vibronic spectra of porphyrins

The connection is made between macrocycle deformations obtained from normal-coordinate structural decomposition (NSD) and the vibronic molecular states and spectra of porphyrins. NSD is a procedure that provides a description of the time-averaged distortion of a porphyrin in terms of normal-coordinate displacements from a D4 h symmetric reference structure. Expressions for the optical absorption spectra with vibrational structure are developed with these NSD-determined deformations as parameters. Using these expressions, the effects of macrocycle distortion on the UV-vis absorption spectra of porphyrins are determined.

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.

Charge-Transfer Interactions in Organic Functional Materials

Our goal in this review is three-fold. First, we provide an overview of a number of quantum-chemical methods that can abstract charge-transfer (CT) information on the excited-state species of organic conjugated materials, which can then be exploited for the understanding and design of organic photodiodes and solar cells at the molecular level. We stress that the Composite-Molecule (CM) model is useful for evaluating the electronic excited states and excitonic couplings of the organic molecules in the solid state. We start from a simple polyene dimer as an example to illustrate how interchain separation and chain size affect the intercahin interaction and the role of the charge transfer interaction in the excited state of the polyene dimers. With the basic knowledge from analysis of the polyene system, we then study more practical organic materials such as oligophenylenevinylenes (OPV_n), oligothiophenes (OT_n), and oligophenylenes (OP_n). Finally, we apply this method to address the delocalization pathway (through-bond and/or through-space) in the lowest excited state for cyclophanes by combining the charge-transfer contributions calculated on the cyclophanes and the corresponding hypothetical molecules with tethers removed. This review represents a step forward in the understanding of the nature of the charge-transfer interactions in the excited state of organic functional materials.

Interchain interactions in organic conjugated dimers: a composite-molecule approach

A theoretical composite-molecule (CM) model is adopted for evaluating the electronic excited states and excitonic couplings of cofacial conjugated dimers where the contributions of charge-transfer (CT) exciton, unavailable by the commonly used supermolecular approach due to the inadequate basis set construction, can be unambiguously identified within this methodology. This method builds up with the basis set of individual molecules and then constructs combined electronic states for the dimer by considering intermolecular interactions including charge-transfer interactions. The dependences of the matrix elements on intermolecular distance and conjugation length are examined. At the short distance region between two of the polyene molecules in the dimer, the CT transitions are apparently mixing to both first and second excited states. Also, some of the matrix elements for the mixing of CT transitions with local transitions which related to the second excited state are found to be considerably larger than the exciton-type elements. An interesting finding is that with increasing the chain size the CT contribution to the second excited state reveals a minimum and indicates HOMO to LUMO charge transfer is not the major CT contribution to the second excited state in the face-to-face polyene dimer with larger chain size and interchain separation in the region of 3.6-4.0 A. A detail analysis reveals that HOMO-1 to LUMO and HOMO to LUMO+1 charge transfers are major CT contributions to the second excited state in the condition under study.

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.

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.

Resonance Raman study of free-base tetraphenylporphine and its dication

Resonance Raman spectra of free-base tetraphenylporphine and its dication obtained with 441.6, 476.5, 488.0 and 514.5 nm excitation lines in the frequency region 100-1625 cm(-1) are reported. Some bands due to in-plane and out-of-plane vibrational modes, which are symmetry forbidden under ideal D(2h), are also seen in the Raman spectra of these molecules. These bands arise due to dynamic and/or static coupling of out-of-plane modes with the allowed in-plane modes. Dynamic coupling may be facilitated by the proton tunneling, while static coupling is due to out-of-plane distortion in the geometrical structure of the molecule. Shift in the positions for certain bands in the Raman spectra of dication are interpreted on the basis of electronic changes due to sharing of electrons of the B(1u) orbital by the two added protons.

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.

Red shifting due to nonplanarity in alkylporphyrins: solid-state polarized UV-vis spectra and ZINDO calculations of two nickel(II)octaethylporphyrins

We present the first demonstration of red shifting upon nonplanarity in alkylporphyrins using two pure conformations having known structures with identical substituents. The traditional view about the relationship of spectral red shifting to nonplanar deformation in porphyrins has been that the deformation from planar to nonplanar forms is in itself the cause of the shifting, but recently this view has been challenged. Among the new arguments is that the substituents required to effect conformational change also bring about nuclear rearrangements in the porphyrin complex which is the actual cause of the red shifting. Octaethylporphyrinatonickel(II), however, exists in both planar and ruffled forms which are determined only by the crystal structure, thus making the issue of different substituents moot. Using a polarized specular reflectance UV-vis microspectrophotometer, we have obtained polarized spectra of pure, solid samples of both forms of NiOEP. We find Soret band red shifting in the solid state that is much larger than previous reports of solution spectra and also report Q-band red shifting. We performed ZINDO calculations on monomers and dimers of both forms of NiOEP, based upon reported structures, and have reproduced the reported solution transition energies and our solid-state spectra as well as the red shifts that we and others have found experimentally. We conclude that, at least in this system, red shifting does indeed result primarily from conformation changes in the porphyrin.

The calculated in vitro and in vivo chlorophyll a absorption bandshape

The room temperature absorption bandshape for the Q transition region of chlorophyll a is calculated using the vibrational frequency modes and Franck-Condon (FC) factors obtained by line-narrowing spectroscopies of chlorophyll a in a glassy (Rebane and Avarmaa, Chem. Phys. 1982; 68:191-200) and in a native environment (Gillie et al., J. Phys. Chem. 1989; 93:1620-1627) at low temperatures. The calculated bandshapes are compared with the absorption spectra of chlorophyll a measured in two different solvents and with that obtained in vivo by a mutational analysis of a chlorophyll-protein complex. It is demonstrated that the measured distributions of FC factors can account for the absorption bandshape of chlorophyll a in a hexacoordinated state, whereas, when pentacoordinated, reduced FC coupling for vibrational frequencies in the range 540-850 cm(-1) occurs. The FC factor distribution for pentacoordinated chlorophyll also describes the native chlorophyll a spectrum but, in this case, either a low-frequency mode (nu < 200 cm(-1)) must be added or else the 262-cm(-1) mode must increase in coupling by about one order of magnitude to describe the skewness of the main absorption bandshape.

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.