Control of the redox potential in c-type cytochromes: importance of the entropic contribution

Assessing Thermal Response of Redox Conduction for Anti-Arrhenius Kinetics in a Microbial Cytochrome Nanowire

A micrometers-long helical homopolymer of the outer-membrane cytochrome type S (OmcS) from Geobacter sulfurreducens is proposed to transport electrons to extracellular acceptors in an ancient respiratory strategy of biogeochemical and technological significance. OmcS surprisingly exhibits higher conductivity upon cooling (anti-Arrhenius kinetics), an effect previously attributed to H-bond restructuring and heme redox potential shifts. Herein, the temperature sensitivity of redox conductivity is more thoroughly examined with conventional and constant-redox and -pH molecular dynamics and quantum mechanics/molecular mechanics. A 30 K drop in temperature constituted a weak perturbation to electron transfer energetics, changing electronic couplings (⟨H_mn⟩), reaction free energies (ΔG_mn), reorganization energies (λ_mn), and activation energies (E_a) by at most |0.002|, |0.050|, |0.120|, and |0.045| eV, respectively. Changes in ΔG_mn reflected -0.07 ± 0.03 V shifts in redox potentials that were caused in roughly equal measure by altered electrostatic interactions with the solvent and protein. Changes in intraprotein H-bonding reproduced the earlier observations. Single-particle diffusion and multiparticle steady-state flux models, parametrized with Marcus theory rates, showed that biologically relevant incoherent hopping cannot qualitatively or quantitatively describe electrical conductivity measured by atomic force microscopy in filamentous OmcS. The discrepancy is attributed to differences between solution-phase simulations and solid-state measurements and the need to model intra- and intermolecular vibrations explicitly.

The redox potential of a heme cofactor in Nitrosomonas europaea cytochrome c peroxidase: a polarizable QM/MM study

Redox reactions are crucial to biological processes that protect organisms against oxidative stress. Metalloenzymes, such as peroxidases which reduce excess reactive oxygen species into water, play a key role in detoxification mechanisms. Here we present the results of a polarizable QM/MM study of the reduction potential of the electron transfer heme in the cytochrome c peroxidase of Nitrosomonas europaea. We have found that environment polarization does not substantially affect the computed value of the redox potential. Particular attention has been given to analyzing the role of electrostatic interactions within the protein environment and the solvent on tuning the redox potential of the heme co-factor. We have found that the electrostatic interactions predominantly explain the fluctuations of the vertical ionization/attachment energies of the heme for the sampled configurations, and that the long range electrostatic interactions (up to 40 Å) contribute substantially to the absolute values of the vertical energy gaps.

Mechanistic Studies of Proton-Coupled Electron Transfer in a Calorimetry Cell

Mechanistic studies of proton-coupled electron-transfer (PCET) reactions in proteins are complicated by the challenge of following proton transfer (PT) in these large molecules. Herein we describe the use of isothermal titration calorimetry (ITC) to establish proton involvement in protein redox reactions and the identity of PT sites. We validate this approach with three variants of a heme protein cytochrome c (cyt c) and show that the method yields a wealth of thermodynamic information that is important for characterizing PCET reactions, including reduction potentials, redox-dependent p K_a values, and reaction enthalpies for both electron-transfer (ET) and PT steps. We anticipate that this facile and label-free ITC approach will find widespread applications in studies of other redox proteins and enhance our knowledge of PCET reaction mechanisms.

Fine Tuning of Functional Features of the CuA Site by Loop-Directed Mutagenesis

Here we report the spectroscopic and electrochemical characterization of three novel chimeric Cu_A proteins in which either one or the three loops surrounding the metal ions in the Thermus thermophilus protein have been replaced by homologous human and plant sequences while preserving the set of coordinating amino acids. These conservative modifications mimic basic differences between Cu_A sites from different organisms and allow for fine tuning the energy gap between alternative electronic ground states of Cu_A.. This results in a systematic modulation of thermodynamic and kinetic electron transfer (ET) parameters and in the selection of one of two possible redox-active molecular orbitals, which differ in the ET reorganization energy by a factor of 2. Moreover, the ET mechanism is found to be frictionally controlled, and the modifications introduced into the different chimeras do not affect the frictional activation parameter.

Thermodynamic and kinetic properties of the outer membrane cytochrome OmcF, a key protein for extracellular electron transfer in Geobacter sulfurreducens

Gene knock-out studies on Geobacter sulfurreducens have shown that the monoheme c-type cytochrome OmcF is essential for the extracellular electron transfer pathways involved in the reduction of iron and uranium oxy-hydroxides, as well as, on electricity production in microbial fuel cells. A detailed electrochemical characterization of OmcF was performed for the first time, allowing attaining kinetics and thermodynamic data. The heterogeneous electron transfer rate constant was determined at pH 7 (0.16 ± 0.01 cm s^-1) indicating that the protein displays high electron transfer efficiency compared to other monoheme cytochromes. The pH dependence of the redox potential indicates that the protein has an important redox-Bohr effect in the physiological pH range for G. sulfurreducens growth. The analysis of the structures of OmcF allowed us to assign the redox-Bohr centre to the side chain of His⁴⁷ residue and its pK_a values in the reduced and oxidized states were determined (pK_ox = 6.73; pK_red = 7.55). The enthalpy, entropy and Gibbs free energy associated with the redox transaction were calculated, pointing the reduced form of the cytochrome as the most favourable. The data obtained indicate that G. sulfurreducens cells evolved to warrant a down-hill electron transfer from the periplasm to the outer-membrane associated cytochrome OmcF.

The structure of PccH from Geobacter sulfurreducens - a novel low reduction potential monoheme cytochrome essential for accepting electrons from an electrode

The structure of cytochrome c (GSU3274) designated as PccH from Geobacter sulfurreducens was determined at a resolution of 2.0 Å. PccH is a small (15 kDa) cytochrome containing one c-type heme, found to be essential for the growth of G. sulfurreducens with respect to accepting electrons from graphite electrodes poised at -300 mV versus standard hydrogen electrode. with fumarate as the terminal electron acceptor. The structure of PccH is unique among the monoheme cytochromes described to date. The structural fold of PccH can be described as forming two lobes with the heme sandwiched in a cleft between the two lobes. In addition, PccH has a low reduction potential of -24 mV at pH 7, which is unusual for monoheme cytochromes. Based on difference in structure, together with sequence phylogenetic analysis, we propose that PccH can be regarded as a first characterized example of a new subclass of class I monoheme cytochromes. The low reduction potential of PccH may enable the protein to be redox active at the typically negative potential ranges encountered by G. sulfurreducens. Because PccH is predicted to be located in the periplasm of this bacterium, it could not be involved in the first step of accepting electrons from the electrode but is very likely involved in the downstream electron transport events in the periplasm.

Domain-swapped dimer of Pseudomonas aeruginosa cytochrome c551: structural insights into domain swapping of cytochrome c family proteins

Cytochrome c (cyt c) family proteins, such as horse cyt c, Pseudomonas aeruginosa cytochrome c₅₅₁ (PA cyt c₅₅₁), and Hydrogenobacter thermophilus cytochrome c₅₅₂ (HT cyt c₅₅₂), have been used as model proteins to study the relationship between the protein structure and folding process. We have shown in the past that horse cyt c forms oligomers by domain swapping its C-terminal helix, perturbing the Met–heme coordination significantly compared to the monomer. HT cyt c₅₅₂ forms dimers by domain swapping the region containing the N-terminal α-helix and heme, where the heme axial His and Met ligands belong to different protomers. Herein, we show that PA cyt c₅₅₁ also forms domain-swapped dimers by swapping the region containing the N-terminal α-helix and heme. The secondary structures of the M61A mutant of PA cyt c₅₅₁ were perturbed slightly and its oligomer formation ability decreased compared to that of the wild-type protein, showing that the stability of the protein secondary structures is important for domain swapping. The hinge loop of domain swapping for cyt c family proteins corresponded to the unstable region specified by hydrogen exchange NMR measurements for the monomer, although the swapping region differed among proteins. These results show that the unstable loop region has a tendency to become a hinge loop in domain-swapped proteins.

Importance of polarization effect in the study of metalloproteins: application of polarized protein specific charge scheme in predicting the reduction potential of azurin

Molecular dynamics (MD) simulation is commonly used in the study of protein dynamics, and in recent years, the extension of MD simulation to the study of metalloproteins is gaining much interest. Choice of force field is crucial in MD studies, and the inclusion of metal centers complicates the process of accurately describing the electrostatic environment that surrounds the redox centre. Herein, we would like to explore the importance of including electrostatic contribution from both protein and solvent in the study of metalloproteins. MD simulations with the implementation of thermodynamic integration will be conducted to model the reduction process of azurin from Pseudomonas aeruginosa. Three charge schemes will be used to derive the partial charges of azurin. These charge schemes differ in terms of the amount of immediate environment, respective to copper, considered during charge fitting, which ranges from the inclusion of copper and residues in the first coordination sphere during density functional theory charge fitting to the comprehensive inclusion of protein and solvent effect surrounding the metal centre using polarized protein-specific charge scheme. From the simulations conducted, the relative reduction potential of the mutated azurins respective to that of wild-type azurin (ΔEcal) were calculated and compared with experimental values. The ΔEcal approached experimental value with increasing consideration of environmental effect hence substantiating the importance of polarization effect in the study of metalloproteins. This study also attests the practicality of polarized protein-specific charge as a computational tool capable of incorporating both protein environment and solvent effect into MD simulations.

The role of key residues in structure, function, and stability of cytochrome-c

Cytochrome-c (cyt-c), a multi-functional protein, plays a significant role in the electron transport chain, and thus is indispensable in the energy-production process. Besides being an important component in apoptosis, it detoxifies reactive oxygen species. Two hundred and eighty-five complete amino acid sequences of cyt-c from different species are known. Sequence analysis suggests that the number of amino acid residues in most mitochondrial cyts-c is in the range 104 ± 10, and amino acid residues at only few positions are highly conserved throughout evolution. These highly conserved residues are Cys14, Cys17, His18, Gly29, Pro30, Gly41, Asn52, Trp59, Tyr67, Leu68, Pro71, Pro76, Thr78, Met80, and Phe82. These are also known as "key residues", which contribute significantly to the structure, function, folding, and stability of cyt-c. The three-dimensional structure of cyt-c from ten eukaryotic species have been determined using X-ray diffraction studies. Structure analysis suggests that the tertiary structure of cyt-c is almost preserved along the evolutionary scale. Furthermore, residues of N/C-terminal helices Gly6, Phe10, Leu94, and Tyr97 interact with each other in a specific manner, forming an evolutionary conserved interface. To understand the role of evolutionary conserved residues on structure, stability, and function, numerous studies have been performed in which these residues were substituted with different amino acids. In these studies, structure deals with the effect of mutation on secondary and tertiary structure measured by spectroscopic techniques; stability deals with the effect of mutation on T m (midpoint of heat denaturation), ∆G D (Gibbs free energy change on denaturation) and folding; and function deals with the effect of mutation on electron transport, apoptosis, cell growth, and protein expression. In this review, we have compiled all these studies at one place. This compilation will be useful to biochemists and biophysicists interested in understanding the importance of conservation of certain residues throughout the evolution in preserving the structure, function, and stability in proteins.

Photosynthetic cytochrome c550

Cytochrome c550 (cyt c550) is a membrane component of the PSII complex in cyanobacteria and some eukaryotic algae, such as red and brown algae. Cyt c550 presents a bis-histidine heme coordination which is very unusual for monoheme c-type cytochromes. In PSII, the cyt c550 with the other extrinsic proteins stabilizes the binding of Cl(-) and Ca(2+) ions to the oxygen evolving complex and protects the Mn(4)Ca cluster from attack by bulk reductants. The role (if there is one) of the heme of the cyt c550 is unknown. The low midpoint redox potential (E(m)) of the purified soluble form (from -250 to -314mV) is incompatible with a redox function in PSII. However, more positive values for the Em have been obtained for the cyt c550 bound to the PSII. A very recent work has shown an E(m) value of +200mV. These data open the possibility of a redox function for this protein in electron transfer in PSII. Despite the long distance (22Å) between cyt c550 and the nearest redox cofactor (Mn(4)Ca cluster), an electron transfer reaction between these components is possible. Some kind of protective cycle involving a soluble redox component in the lumen has also been proposed. The aim of this article is to review previous studies done on cyt c550 and to consider its function in the light of the new results obtained in recent years. The emphasis is on the physical properties of the heme and its redox properties. This article is part of a Special Issue entitled: Photosynthesis Research for Sustainability: from Natural to Artificial.

Modulating heme redox potential through protein-induced porphyrin distortion

Hemoproteins are ubiquitous in biology and are commonly involved in critical processes such as electron transfer, oxidative phosphorylation, and signal transduction. Both the protein environment and the heme cofactor contribute to generate the range of chemical properties needed for these diverse functions. Using the heme nitric oxide/oxygen binding (H-NOX) protein from the thermophilic bacterium Thermoanaerobacter tengcongensis, we have shown that heme electronic properties can be modulated by porphyrin distortion within the same protein scaffold without changing the heme ligation state or heme environment. The degree of heme distortion was found to be directly correlated to the electron density at the heme iron, as evidenced by dramatic changes in the heme redox potential and pK(a) of the distal ligand ((-)OH vs H(2)O). Protein-induced porphyrin distortion represents a new strategy to rationally tune the electronic properties of protein-bound porphyrins and could be used to engineer proteins with desired reactivity or functionality.

Electrochemistry of nanozeolite-immobilized cytochrome c in aqueous and nonaqueous solutions

The electrochemical properties of cytochrome c (cyt c) immobilized on multilayer nanozeolite-modified electrodes have been examined in aqueous and nonaqueous solutions. Layers of Linde type-L zeolites were assembled on indium tin oxide (ITO) glass electrodes followed by the adsorption of cyt c, primarily via electrostatic interactions, onto modified ITO electrodes. The heme protein displayed a quasi-reversible response in aqueous solution with a redox potential of +324 mV (vs NHE), and the surface coverage (Gamma*) increased linearly for the first four layers and then gave a nearly constant value of 200 pmol cm(-2). On immersion of the modified electrodes in 95% (v/v) nonaqueous solutions, the redox potential decreased significantly, a decrease that originated from changes in both the enthalpy and entropy of reduction. On reimmersion of the modified electrode in buffer, the faradic response immediately returned to its original value. These results demonstrate that nanozeolites are potential stable supports for redox proteins and enzymes.

Control of reduction thermodynamics in [2Fe-2S] ferredoxins Entropy-enthalpy compensation and the influence of surface mutations

The reaction thermodynamics for the one-electron reduction of the [2Fe-2S] cluster of both human ferredoxin and various surface point mutants, in which each of the negatively charged residues Asp72, Glu73, Asp76, and Asp79 were converted to Ala, have been determined by variable temperature spectroelectrochemical measurements. The above are conserved residues that have been implicated in interactions between the vertebrate-type ferredoxins and their redox partners. In all cases, and similar to other 2Fe-ferredoxins, the reduction potentials are negative as a result of both an enthalpic and entropic stabilization of the oxidized state. Although all Hs Fd mutants, with the exception of Asp72Ala, show slightly higher E degrees ' values than that of wild type Hs Fd, according to expectations for a purely electrostatic model, they exhibit changes in the H degrees '(rc) values that are electrostatically counter-intuitive. The observation of enthalpy-entropy compensation within the protein series indicates that the mutation-induced changes in H degrees '(rc) and S degrees '(rc) are dominated by reduction-induced solvent reorganization effects. Protein-based entropic effects are likely to be responsible for the low E degrees ' value of D72A.

Ultrahigh-resolution study on Pyrococcus abyssi rubredoxin: II. Introduction of an O–H···Sγ–Fe hydrogen bond increased the reduction potential by 65 mV

The effect of D-H...S(gamma)-Fe hydrogen bonding on the reduction potential of rubredoxin was investigated by the introduction of an O-H...S(gamma)-Fe hydrogen bond on the surface of Pyrococcus abyssi rubredoxin. The formation of a weak hydrogen bond between Ser44-O(gamma) and Cys42-S(gamma) in mutant W4L/R5S/A44S increased the reduction potential by 56 mV. When side effects of the mutation were taken into account, the contribution of the additional cluster hydrogen bond to the reduction potential was estimated to be +65 mV. The structural analysis was based on ultrahigh-resolution structures of oxidized P. abyssi rubredoxin W4L/R5S and W4L/R5S/A44S refined to 0.69 and 0.86 A, respectively.

Redox reactivity of the heme Fe3+/Fe2+ couple in native myoglobins and mutants with peroxidase-like activity

The reaction enthalpy and entropy for the one-electron reduction of the ferric heme in horse heart and sperm whale aquometmyoglobins (Mb) have been determined exploiting a spectroelectrochemical approach. Also investigated were the T67R, T67K, T67R/S92D and T67R/S92D Mb-H variants (the latter containing a protoheme-L: -histidine methyl ester) of sperm whale Mb, which feature peroxidase-like activity. The reduction potential (E degrees ') in all species consists of an enthalpic term which disfavors Fe(3+) reduction and a larger entropic contribution which instead selectively stabilizes the reduced form. This behavior differs from that of the heme redox enzymes and electron transport proteins investigated so far. The reduction thermodynamics in the series of sperm whale Mb variants show an almost perfect enthalpy-entropy compensation, indicating that the mutation-induced changes in DeltaH(o')(rc) and DeltaS(o')(rc) are dominated by reduction-induced solvent reorganization effects. The modest changes in E degrees ' originate from the enthalpic effects of the electrostatic interactions of the heme with the engineered charged residues. The small influence that the mutations exert on the reduction potential of myoglobin suggests that the increased peroxidase activity of the variants is not related to changes in the redox reactivity of the heme iron, but are likely related to a more favored substrate orientation within the distal heme cavity.

The redox couple of the cytochrome c cyanide complex: the contribution of heme iron ligation to the structural stability, chemical reactivity, and physiological behavior of horse cytochrome c

Contrary to most heme proteins, ferrous cytochrome c does not bind ligands such as cyanide and CO. In order to quantify this observation, the redox potential of the ferric/ferrous cytochrome c-cyanide redox couple was determined for the first time by cyclic voltammetry. Its E0' was -240 mV versus SHE, equivalent to -23.2 kJ/mol. The entropy of reaction for the reduction of the cyanide complex was also determined. From a thermodynamic cycle that included this new value for the cyt c cyanide complex E0', the binding constant of cyanide to the reduced protein was estimated to be 4.7 x 10(-3) L M(-1) or 13.4 kJ/mol (3.2 kcal/mol), which is 48.1 kJ/mol (11.5 kcal/mol) less favorable than the binding of cyanide to ferricytochrome c. For coordination of cyanide to ferrocytochrome c, the entropy change was earlier experimentally evaluated as 92.4 J mol(-1) K(-1) (22.1 e.u.) at 25 K, and the enthalpy change for the same net reaction was calculated to be 41.0 kJ/mol (9.8 kcal/mol). By taking these results into account, it was discovered that the major obstacle to cyanide coordination to ferrocytochrome c is enthalpic, due to the greater compactness of the reduced molecule or, alternatively, to a lower rate of conformational fluctuation caused by solvation, electrostatic, and structural factors. The biophysical consequences of the large difference in the stabilities of the closed crevice structures are discussed.

Axial ligation and polypeptide matrix effects on the reduction potential of heme proteins probed on their cyanide adducts

The enthalpic and entropic changes accompanying the reduction reaction of the six-coordinate cyanide adducts of cytochrome c, microperoxidase-11 and a few plant peroxidases were measured electrochemically. Once the compensating changes in reduction enthalpy and entropy due to solvent reorganization effects are factorized out, it is found that cyanide binding stabilizes enthalpically the ferriheme following the order: cyochrome c > peroxidase > microperoxidase-11. The effect is inversely correlated to the solvent accessibility of the heme. Comparison of the reduction thermodynamics for the cyanide adducts of cytochrome c and plant peroxidases with those for microperoxidase-11 and myoglobin, respectively, yielded an estimate of the consequences of protein encapsulation and of the anionic character of the proximal histidine on the reduction potential of the heme-cyanide group. Insertion of the heme-CN group into the folded peptide chain of cyt c induces an enthalpy-based decrease in E degrees ' of approximately 100 mV, consistent with the lower net charge of the oxidized as compared to the reduced iron center, whereas a full imidazolate character of the proximal histidine stabilizes enthalpically the ferriheme by approximately 400 mV. The latter value should be best considered as an upper limit since it also includes some solvation effects arising from the nature of the protein systems being compared.

Changes in hydrogen-bond strengths explain reduction potentials in 10 rubredoxin variants

The rubredoxin from Clostridium pasteurianum (CpRd) provides an excellent system for investigating how the protein sequence modulates the reduction potential of the active site in an iron-sulfur protein. (15)N NMR spectroscopy has allowed us to determine with unprecedented accuracy the strengths of all six key hydrogen bonds between protein backbone amides and the sulfur atoms of the four cysteine residues that ligate the iron in the oxidized (Fe(III)) and reduced (Fe(II)) forms of wild-type CpRd and nine mutants (V44G, V44A, V44I, V44L, V8G, V8A, V8I, V8L, and V8G/V44G). The length (or strength) of each hydrogen bond was inferred from the magnitude of electron spin delocalized across the hydrogen bond from the iron atom onto the nitrogen. The aggregate lengths of these six hydrogen bonds are shorter in both oxidation states in variants with higher reduction potential than in those with lower reduction potential. Differences in aggregate hydrogen bonding upon reduction correlate linearly with the published reduction potentials for the 10 CpRd variants, which span 126 mV. Sequence effects on the reduction potential can be explained fully by their influence on hydrogen-bond strengths.

Direct electrochemistry of engineered cytochrome b562 molecules with a ligand binding pocket

The rapid and reversible electron transfer reaction of cytochrome b562 was observed at an In2O3 electrode. The estimated heterogeneous electron transfer rate constant (k0') was k0' > or = 5.0 x 10(-3) cm s(-1) at pH 6.5. When the methionine-7 (Met-7) residue, which coordinates to the heme iron as an axial ligand, of the wild-type cytochrome b562 was replaced by an Ala or Gly residue, a water molecule bound to the heme iron and the electron transfer rate constants decreased to 1.3 x 10(-3) and 1.8 x 10(-3) cm s(-1), respectively. This decrease in the electron transfer rate would be due to the larger reorganization energy for the structural change at the redox site. The midpoint potential of cytochrome b562 was shifted negatively by approximately 135 mV by replacing Met-7 with Ala or Gly. Similar dissociation kinetics of cyanide for the mutated molecules as compared to native myoglobin was obtained.

Electric-field-induced redox potential shifts of tetraheme cytochromes c3 immobilized on self-assembled monolayers: surface-enhanced resonance Raman spectroscopy and simulation studies

The tetraheme protein cytochrome c(3) (Cyt-c(3)) from Desulfovibrio gigas, immobilized on a self-assembled monolayer (SAM) of 11-mercaptoundecanoic acid, is studied by theoretical and spectroscopic methods. Molecular dynamics simulations indicate that the protein docks to the negatively charged SAM via its lysine-rich domain around the exposed heme IV. Complex formation is associated with only little protein structural perturbations. This finding is in line with the resonance Raman and surface-enhanced resonance Raman (SERR) spectroscopic results that indicate essentially the same heme pocket structures for the protein in solution and adsorbed on SAM-coated Ag electrodes. Electron- and proton-binding equilibrium calculations reveal substantial negative shifts of the redox potentials compared to the protein in solution. The magnitude of these shifts decreases in the order heme IV (-161 mV) > heme III (-73 mV) > heme II (-57 mV) > heme I (-26 mV), resulting in a change of the order of reduction. These shifts originate from the distance-dependent electrostatic interactions between the SAM headgroups and the individual hemes, leading to a stabilization of the oxidized forms. The results of the potential-dependent SERR spectroscopic analyses are consistent with the theoretical predictions and afford redox potential shifts of -160 mV (heme IV), -90 mV (heme III), -70 mV (heme II), and +20 mV (heme I) relative to the experimental redox potentials for Cyt-c(3) in solution. SERR spectroscopic experiments reveal electric-field-induced changes of the redox potentials also for the structurally very similar Cyt-c(3) from Desulfovibrio vulgaris, although the shifts are somewhat smaller compared to Cyt-c(3) from D. gigas. This study suggests that electric-field-induced redox potential shifts may also occur upon binding to biomembranes or partner proteins and thus may affect biological electron transfer processes.

Redox and redox-coupled processes of heme proteins and enzymes at electrochemical interfaces

Modern bioelectrochemical methods rely upon the immobilisation of redox proteins and enzymes on electrodes coated with biocompatible materials to prevent denaturation. However, even when protein denaturation is effectively avoided, heterogeneous protein electron transfer is often coupled to non-Faradaic processes like reorientation, conformational transitions or acid-base equilibria. Disentangling these processes requires methods capable of probing simultaneously the structure and reaction dynamics of the adsorbed species. Here we provide an overview of the recent developments in Raman and infrared surface-enhanced spectroelectrochemical techniques applied to the study of soluble and membrane bound redox heme proteins and enzymes. Possible biological implications of the findings are critically discussed.

Electron transfer in ferredoxin: are tunneling pathways evolutionarily conserved?

A theoretical study of electron transfer (ET) pathways in a recently crystallized Clostridium acidurici ferredoxin is reported. The electronic structure of the protein complex is treated at the semiempirical extended Hückel level, and the tunneling pathways are calculated with the rigorous quantum mechanical method of tunneling currents. The model predicts two pathways between the two [4Fe-4S] cubanes: a strong one running directly from Cys(14) to Cys(43) and a weaker one from Cys(14) via Ile(23) to Cys(18), whereas other amino acids do not play a significant role in the electron tunneling. The cysteine ligands conduct almost all of the current when Ile(23) is mutated to valine in silico, so that there is no appreciable change in the ET rate. The calculated value of the transfer matrix element is consistent with the experimentally determined rate of transfer. Results of the sequence analysis performed on this ferredoxin reveal that Ile(23) is a highly variable amino acid compared with the cubane-ligating cysteine amino acids, even though Ile(23) lies directly between the donor and acceptor complexes. We further argue that the homologous proteins with a [3Fe-4S] cofactor, which does not have one of the four cysteine ligands, use the same tunneling pathways as those in this ferredoxin, on the basis of the high homology as well as the absolute conservation of Cys(14) and Cys(43) which serve as the main tunneling conduit. Our results explain why mutation of amino acids around and between the donor and acceptor cubane clusters, including that of Ile(23), does not appreciably affect the rate of transfer and add support to the proposal that there exist evolutionarily conserved electron tunneling pathways in biological ET reactions.

Redox properties of cytochrome c

The redox properties of cytochromes (cyt) c, a ubiquitous class of heme-containing electron transport proteins, have been extensively investigated over the last two decades. The reduction potential (E degrees') is central to the chemistry of cyt c for two main reasons. First, E degrees' influences both the thermodynamic and kinetic aspects of the electron exchange reaction with redox partners. Second, this thermodynamic parameter is remarkably sensitive to changes in the properties of the heme and the protein matrix, and hence can be profitably used for the investigation of the solution chemistry of cyt c. This research area owes much to the exploitation of voltammetric techniques for the determination of E degrees' for metalloproteins, which dates back to the late 1970s. Since then, much effort has been devoted to the comprehension of the molecular factors that control E degrees' in cyt c, which include first coordination sphere effects on the heme iron, the interactions of the heme group with the surrounding polypeptide chain and the solvent, and also include medium effects related to the nature and ionic composition of the solvent, pH, the presence of potential protein ligands, and the temperature. This article provides an overview of the most significant advances made in this field recently.

Adrenodoxin: structure, stability, and electron transfer properties

Adrenodoxin is an iron-sulfur protein that belongs to the broad family of the [2Fe-2S]-type ferredoxins found in plants, animals and bacteria. Its primary function as a soluble electron carrier between the NADPH-dependent adrenodoxin reductase and several cytochromes P450 makes it an irreplaceable component of the steroid hormones biosynthesis in the adrenal mitochondria of vertebrates. This review intends to summarize current knowledge about structure, function, and biochemical behavior of this electron transferring protein. We discuss the recently solved first crystal structure of the vertebrate-type ferredoxin, the truncated adrenodoxin Adx(4-108), that offers the unique opportunity for better understanding of the structure-function relationships and stabilization of this protein, as well as of the molecular architecture of [2Fe-2S] ferredoxins in general. The aim of this review is also to discuss molecular requirements for the formation of the electron transfer complex. Essential comparison between bacterial putidaredoxin and mammalian adrenodoxin will be provided. These proteins have similar tertiary structure, but show remarkable specificity for interactions only with their own cognate cytochrome P450. The discussion will be largely centered on the protein-protein recognition and kinetics of adrenodoxin dependent reactions.

Temperature dependence of the formal reduction potential of putidaredoxin

Putidaredoxin (Pdx), a [2Fe-2S] redox protein of size M(r) 11,600, transfers two electrons in two separate steps from the flavin containing putidaredoxin reductase to the heme protein, cytochrome CYP101 in the P450cam catalytic cycle. It has recently come to light, through NMR measurements, that there can be appreciable differences in the Pdx conformational dynamics between its reduced and oxidized states. The redox reaction entropy, deltaS(0')rc = (S(0')Pdx(r)-S(0')Pdx(0)), as determined from measurements of the variation in formal potential with temperature, E0'(T), provides a measure of the strength of this influence on Pdx function. We designed a spectroelectrochemical cell using optically transparent tin oxide electrodes, without fixed or diffusible mediators, to measure E0'(T) over the temperature range 0-40 degrees C. The results indicate that the redox reaction entropy for Pdx is biphasic, decreasing from -213 +/- 27 J mol(-1) K(-1) over 0-27 degrees C, to -582 +/- 150 J mol(-1) K (-1) over 27-40 degrees C. These redox reaction entropy changes are significantly more negative than the changes reported for most cytochromes, although our measurement over the temperature interval 0-27 degrees C is in the range reported for other iron-sulfur proteins. This suggests that Pdx (and other ferredoxins) is a less rigid system than monohemes, and that redox-linked changes in conformation, and/or conformational dynamics, impart to these proteins the ability to interact with a number of redox partners.

Simulation of electron-proton coupling with a Monte Carlo method: application to cytochrome c3 using continuum electrostatics

A new method is presented for simulating the simultaneous binding equilibrium of electrons and protons on protein molecules, which makes it possible to study the full equilibrium thermodynamics of redox and protonation processes, including electron-proton coupling. The simulations using this method reflect directly the pH and electrostatic potential of the environment, thus providing a much closer and realistic connection with experimental parameters than do usual methods. By ignoring the full binding equilibrium, calculations usually overlook the twofold effect that binding fluctuations have on the behavior of redox proteins: first, they affect the energy of the system by creating partially occupied sites; second, they affect its entropy by introducing an additional empty/occupied site disorder (here named occupational entropy). The proposed method is applied to cytochrome c3 of Desulfovibrio vulgaris Hildenborough to study its redox properties and electron-proton coupling (redox-Bohr effect), using a continuum electrostatic method based on the linear Poisson-Boltzmann equation. Unlike previous studies using other methods, the full reduction order of the four hemes at physiological pH is successfully predicted. The sites more strongly involved in the redox-Bohr effect are identified by analysis of their titration curves/surfaces and the shifts of their midpoint redox potentials and pKa values. Site-site couplings are analyzed using statistical correlations, a method much more realistic than the usual analysis based on direct interactions. The site found to be more strongly involved in the redox-Bohr effect is propionate D of heme I, in agreement with previous studies; other likely candidates are His67, the N-terminus, and propionate D of heme IV. Even though the present study is limited to equilibrium conditions, the possible role of binding fluctuations in the concerted transfer of protons and electrons under nonequilibrium conditions is also discussed. The occupational entropy contributions to midpoint redox potentials and pKa values are computed and shown to be significant.

15N-labelling and preliminary heteronuclear NMR study of Desulfovibrio vulgaris Hildenborough cytochrome c553

When using heteronuclear NMR, 15N-labelling is necessary for structural analysis, dynamic studies and determination of complex formation. The problems that arise with isotopic labelling of metalloproteins are due to their complex maturation process, which involves a large number of factors. Cytochromes c are poorly expressed in Escherichia coli and the overexpression that is necessary for 15N-labelling, requires an investigation of the expression host and special attention to growth conditions. We have succeeded in the heterologous expression and the complete and uniform isotopic 15N-labelling of the cytochrome c553 from Desulfovibrio vulgaris Hildenborough, in a sulphate-reducing bacterium, D. desulfuricans G200, by using a growth medium combining 15N-ammonium chloride and 15N-Celtone. These conditions allowed us to obtain approximately 0.8 mg x L-1 of pure labelled cytochrome c553. 1H and 15N-assignments for both the oxidized and the reduced states of cytochrome c553 were obtained from two-dimensional heteronuclear experiments. Pseudocontact effects due to the haem Fe3+ have been analysed for the first time through 15N and 1H chemical shifts in a c-type cytochrome.

Electrochemical behaviour of cytochrome c at electrically heated microelectrodes

The structural changes in cytochrome c with temperature have been been followed using a recently developed electrically-heated microelectrode sensor. Differential pulse voltammetry was used to perform electrochemical measurements of cytochrome c oxidation at different temperatures at heated bare gold electrodes contained in phosphate-buffered cytochrome c solution at room temperature. The voltammetric response shows the onset of unfolding and a marked dependence of the signal on electrode temperature. This augurs well for applications of heated electrodes as local probes in the study of the temperature dependence of electron transfer processes of other redox proteins, avoiding problems of bulk deterioration.

Using entropies of reaction to predict changes in protein stability: tyrosine-67-phenylalanine variants of rat cytochrome c and yeast Iso-1 cytochromes c

Using the voltammetric method of square-wave voltammetry, a direct electrochemical examination was made of the wild type and Tyr67Phe mutant of both rat cytochrome c and yeast iso-1-cytochrome c. In addition to determining the equilibrium reduction potential (E0') for each cytochrome, the entropy of reaction, deltaS0'(Rxn)(deltaS0'(Rxn) = S0'(Red) - S0'(Ox)), for the reduction process was determined via the non-isothermal method. Having determined deltaS0'(Rxn) and E0', deltaH0' was calculated. For rat cytochrome c, it was found that deltaS0'(Rxn) = -43 J mol(-1) K(-1) for the wild type and -53 J mol(-1) K(-1) for the Tyr67Phe variant, with the deltaH0' for both the wild type and variant nearly identical, indicating that the changes in reduction potential and probably stability are due to changes in deltaS0'(Rxn). In contrast the measured deltaS0'(Rxn) for yeast iso-1-cytochrome c demonstrated significant changes in both entropic and enthalpic contributions in going from wild type to mutant cytochrome c. The entropy of reaction provides information regarding the relative degree of solvation, and very likely the degree of compactness, of the oxidized state versus the reduced state of the redox protein. A thermodynamic scheme and stability derivation are presented that show how the entropies of reaction of wild type versus variant cytochromes contribute to and predict changes in stability in going from oxidized to reduced protein. For yeast iso-1-cytochrome c, the thermodynamically predicted change in stability was very close to the experimentally observed value, based on previous differential scanning calorimetric stability measurements. While such data is not available for rat cytochrome c, consideration of the enormously increased local stability of the rat oxidized cytochrome c variant predicts that the reduced rat variant will be even more stable than the already stabilized oxidized variant.

Structural organization in peptide fragments of cytochrome c by heme binding

Prosthetic groups are often important structural organizers of proteins as well as essential functional components. Insertion of prosthetic groups is usually spontaneous, and implies an apoprotein that is partially preorganized to provide a recognition surface for specific binding. Cytochrome c is distinguished by having its heme attached by a dedicated heme lyase through thioether links to cysteine side-chains, and the apoprotein shows no evidence of preorganization under physiological conditions. Nevertheless, addition of heme to two short fragments of cytochrome c enhances helical structure substantially (from approximately 8% to approximately 22%), an effect that depends on iron ligation but not thioether linkage. The helical segments in the corresponding parts of the native holoprotein have little contact surface with heme, implying that the increased helical structure in the fragment complex may depend on tertiary interactions. The absence of the intervening polypeptide chain suggests that the complex represents a relatively independent folded subdomain.

Molecular dynamics simulation of cytochrome c3: studying the reduction processes using free energy calculations

The tetraheme cytochrome c3 from Desulfovibrio vulgaris Hildenborough is studied using molecular dynamics simulation studies in explicit solvent. The high heme content of the protein, which has its core almost entirely made up of c-type heme, presents specific problems in the simulation. Instability in the structure is observed in long simulations above 1 ns, something that does not occur in a monoheme cytochrome, suggesting problems in heme parametrization. Given these stability problems, a partially restrained model, which avoids destruction of the structure, was created with the objective of performing free energy calculations of heme reduction, studies that require long simulations. With this model, the free energy of reduction of each individual heme was calculated. A correction in the long-range electrostatic interactions of charge groups belonging to the redox centers had to be made in order to make the system physically meaningful. Correlation is obtained between the calculated free energies and the experimental data for three of four hemes. However, the relative scale of the calculated energies is different from the scale of the experimental free energies. Reasons for this are discussed. In addition to the free energy calculations, this model allows the study of conformational changes upon reduction. Even if the precise details of the structural changes that take place in this system upon individual heme reduction are probably out of the reach of this study, it appears that these structural changes are small, similarly to what is observed for other redox proteins. This does not mean that their effect is minor, and one example is the conformational change observed in propionate D from heme I when heme II becomes reduced. A motion of this kind could be the basis of the experimentally observed cooperativity effects between heme reduction, namely positive cooperativity.

Entropies of redox reactions between proteins and mediators: the temperature dependence of reversible electrode potentials in aqueous buffers

The temperature dependencies of the reversible electrode potentials for a number of charge transfer reactions of redox mediators were used to evaluate the corresponding charge transfer entropies in Tris-HCl (pH 8) buffer. The redox mediator thermodynamic data, along with reaction enthalpy data for mediator redox protein electron transfer, were used to evaluate the charge transfer entropy for the cytochrome c redox couple [(cytc)ox/(cytc)red] in Tris-HCl (pH 8) buffer and were found to be equal to -16 cal/degrees K mol. Reversible electrode potentials at 298 degrees K for the redox mediator half-reactions were observed to vary from -528 to +657 mV (vs NHE). Charge transfer entropies were observed to depend upon the structure of the redox mediators and to vary from -13.8 to -29.7 cal/degrees K mol for a closely related series of organic dications (viologens) and a value of -43.6 cal/degrees K mol was observed for the [Fe(CN)6]3-/[Fe(CN)6]4- couple under the same conditions. A procedure for determining charge transfer entropies of protein redox couples which cannot be studied by direct electrochemical methods is outlined. The factors contributing to the magnitude of the charge transfer entropies are discussed.