Protein Matrix and Dielectric Effect in Cytochromec

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.

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.

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.

Local stability of Rhodobacter capsulatus cytochrome c2 probed by solution phase hydrogen/deuterium exchange and mass spectrometry

The hydrogen/deuterium exchange kinetics of Rhodobacter capsulatus cytochrome c2 have been determined using mass spectrometry. As expected, the relative domain stability was generally similar to that of the cytochrome c2 structural homolog, horse heart cytochrome c, but we were able to find evidence to support the presence of a second, small beta-sheet not found in the horse cytochrome, which stabilizes a structural region dominated by Omega loops. Importantly, we find that the so-called hinge region, comprised of 15 amino acids, which include the methionine sixth heme ligand (M96), is destabilized on oxidation, and this destabilization is propagated to a portion of the second Omega loop, most likely through perturbation of two hydrogen bonds that couple these two domains in the three dimensional structure. The mutation of a lysine at position 93 to proline amplifies the destabilization observed on oxidation of the wild-type cytochrome c2 and results in further destabilization observed in regions 52-60, 75-82, and 83-97. This suggests that hydrogen bond interactions involving two bound waters, the T94 hydroxyl, the front heme propionate and the Y75 hydroxyl, are significantly compromised upon mutation. In summary, these observations are consistent with the approximately 20-fold increase in the movement of the hinge away from the heme face in the oxidized cytochrome c2 as determined by ligand binding kinetics. Thus, H/D exchange kinetics can be used to identify relatively subtle structural features and at least in some cases facilitate the understanding of the structural basis of the dynamic properties of proteins.

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.

Structure and redox properties of the haem centre in the C357M mutant of cytochrome P450cam

The effects of site-specific mutation of the axial cysteine (C357M) to a methionine residue in cytochrome P450cam on the enzyme's coordination geometry and redox potential have been investigated. The absorption spectra of the haem centre in the C357M mutant of the enzyme showed close similarity to those of cytochrome c both in the oxidised and reduced forms. A well-defined absorption peak at 695 nm, similar to that seen in the case of cytochrome c and characteristic of methionine ligation to the ferric haem, was observed. The results indicated that the haem of C357M cytochrome P450cam is possibly axially coordinated to a methionine and a histidine, analogously to cytochrome c. The circular dichroism spectra in the visible and the far-UV regions suggested that the tertiary structure of the haem cavity in the C357M mutant cytochrome P450cam was distinctly different from that in the wild-type enzyme or in cytochrome c, although the secondary structure of the mutant remained identical to that of the wild-type cytochrome P450cam. Comparison of the natures of the CD spectra in the 400 nm and 695 nm regions of the C357M mutant of cytochrome P450cam with those of horse cytochrome c suggested (R) chirality at the sulfur atom of the iron-bound methionine residue in the mutant. The redox potential of the haem centre, estimated by redox titration of the C357M mutant, was found to be +260 mV, which is much higher than that in the wild-type enzyme and similar to the redox potential of cytochrome c. This supported the concept that axial ligation of the haem plays the major role in tuning the redox potential of the haem centre in haem proteins.

Control of the Redox Potential of Pseudomonas aeruginosa Cytochrome c551 through the Fe−Met Coordination Bond Strength and pKa of a Buried Heme Propionic Acid Side Chain

Pseudomonas aeruginosa cytochrome c(551) and a series of its mutants exhibiting various thermostabilities have been studied by paramagnetic (1)H NMR and cyclic voltammetry in an effort to elucidate the molecular mechanisms responsible for control of the redox potentials (E degrees ') of the proteins. The study revealed that the E degrees ' value of the protein is regulated by two molecular mechanisms operating independently of each other. One is based on the Fe-Met coordination bond strength in the protein, which is determined by the amino acid side chain packing in the protein, and the other on the pK(a) of the heme 17-propionic acid side chain, which is affected by the electrostatic environment. The former mechanism alters the magnitude of the E degrees ' value throughout the entire pH range, and the latter regulates the pK values reflected by the pH profile of the E degrees ' value. These findings provide novel insights into functional regulation of the protein, which could be utilized for tuning the E degrees ' value of the protein by means of protein engineering.

Distance dependence of interactions between charged centres in proteins with common structural features

Data collected for interactions among redox centres, and interactions between redox centres and acid-base residues in a family of small multihaem cytochromes are analysed. The distance dependent attenuation of the interactions between non-surface charges, for separations that range from 8 to 23 angstroms, can be described by a simple function derived from the Debye-Huckel formalism, fit to 9.5 and 7.6 as values for the relative dielectric constant and Debye length, respectively. However, there is considerable scatter in the data despite the structural similarities among the proteins, which is discussed in the framework of using such simple models in predicting properties of novel proteins.

Inferring functional constraints and divergence in protein families using 3D mapping of phylogenetic information

Comparative sequence analysis has been used to study specific questions about the structure and function of proteins for many years. Here we propose a knowledge-based framework in which the maximum likelihood rate of evolution is used to quantify the level of constraint on the identity of a site. We demonstrate that site-rate mapping on 3D structures using datasets of rhodopsin-like G-protein receptors and alpha- and beta-tubulins provides an excellent tool for pinpointing the functional features shared between orthologous and paralogous proteins. In addition, functional divergence within protein families can be inferred by examining the differences in the site rates, the differences in the chemical properties of the side chains or amino acid usage between aligned sites. Two novel analytical methods are introduced to characterize rate- independent functional divergence. These are tested using a dataset of two classes of HMG-CoA reductases for which only one class can perform both the forward and reverse reaction. We show that functionally divergent sites occur in a cluster of sites interacting with the catalytic residues and that this information should facilitate the design of experimental strategies to directly test functional properties of residues.

Probing electrostatic interactions in cytochrome c using site-directed chemical modification

This communication reports the generation of an electrostatic probe using chemical modification of methionine side chains. The alkylation of methionine by iodoacetamide was achieved in a set of Saccharomyces cerevisiae iso-1-cytochrome c mutants, introducing the nontitratable, nondelocalized positive charge of a carboxyamidomethylmethionine sulfonium (CAMMS) ion at five surface and one buried site in the protein. Changes in redox potential and its variation with temperature were used to calculate microscopic effective dielectric constants operating between the probe and the heme iron. Dielectric constants (epsilon) derived from deltadeltaG values were not useful due to entropic effects, but epsilon(deltadeltaH) gave results that supported the theory. The effect on biological activity of surface derivatization was interpreted in terms of protein-protein interactions. The introduction of an electrostatic probe in cytochrome c often resulted in marked effects on activity with one of two physiological partners: cytochrome c reductase, especially if introduced at position 65, and cytochrome c oxidase, if at position 28.

Protein dynamics and cytochrome c: correlations between ligand vibrations and redox activity

A recently developed method to directly observe specific protein vibrations, based on deuteration, has been employed to examine the redox-dependent structural and fluctional properties of cytochrome c. The dynamics of the protein-based methionine heme ligand were examined by selectively deuterating the ligand's methyl group. The frequency and line width of the C-D bonds were easily observable and shown to be sensitive to mutation-induced changes in the protein redox potential. However, of seven mutants examined, the C-D line widths were independent of the redox-state of the protein. Therefore, although the ligand dynamics depend on the protein's redox state, there are no detected differences in protein dynamics of the oxidized and reduced proteins.