Modeling electron transfer thermodynamics in protein complexes: interaction between two cytochromes c(3)

Redox protein complexes between type I and type II tetraheme cytochromes c(3) from Desulfovibrio vulgaris Hildenborough are here analyzed using theoretical methodologies. Various complexes were generated using rigid-body docking techniques, and the two lowest energy complexes (1 and 2) were relaxed using molecular dynamics simulations with explicit solvent and subjected to further characterization. Complex 1 corresponds to an interaction between hemes I from both cytochromes c(3). Complex 2 corresponds to an interaction between the heme IV from type I and the heme I from type II cytochrome c(3). Binding free energy calculations using molecular mechanics, Poisson-Boltzmann, and surface accessibility methods show that complex 2 is more stable than complex 1. Thermodynamic calculations on complex 2 show that complex formation induces changes in the reduction potential of both cytochromes c(3), but the changes are larger in the type I cytochrome c(3) (the largest one occurring on heme IV, of approximately 80 mV). These changes are sufficient to invert the global titration curves of both cytochromes, generating directionally in electron transfer from type I to type II cytochrome c(3), a phenomenon of obvious thermodynamic origin and consequences, but also with kinetic implications. The existence of processes like this occurring at complex formation may constitute a natural design of efficient redox chains.     

Functional properties of type I and type II cytochromes c3 from Desulfovibrio africanus

Type I cytochrome c(3) is a key protein in the bioenergetic metabolism of Desulfovibrio spp., mediating electron transfer between periplasmic hydrogenase and multihaem cytochromes associated with membrane bound complexes, such as type II cytochrome c(3). This work presents the NMR assignment of the haem substituents in type I cytochrome c(3) isolated from Desulfovibrio africanus and the thermodynamic and kinetic characterisation of type I and type II cytochromes c(3) belonging to the same organism. It is shown that the redox properties of the two proteins allow electrons to be transferred between them in the physiologically relevant direction with the release of energised protons close to the membrane where they can be used by the ATP synthase.     

Studies of the reduction and protonation behavior of tetraheme cytochromes using atomic detail

A comparative study of tetraheme cytochrome c3 molecules from several species was carried out using recently developed theoretical methods based on continuum electrostatics. The binding joint equilibrium of electrons and protons was simulated, revealing the complete thermodynamic aspects of electron-proton coupling in these molecules. The method yields excellent accuracy in terms of midpoint potentials, giving the correct reduction orders in all molecules examined, except for one heme site. The coupling between electrons and protons is shown to be present and significant at physiological pH in all cases. This phenomenon, known as the redox-Bohr effect, though of thermodynamic nature, is shown to have an intrinsic "dynamic" character at the molecular level (in the sense of the empty/occupied fluctuations at the microscopic level), with the binding states of redox and protonatable sites displaying both correlated averages and correlated fluctuations. The protonatable sites more directly involved in the redox-Bohr effect are identified using, among other properties, the statistical correlation between pairs of sites, which automatically reflects indirect effects mediated by other sites. Several sites are identified in this analysis. Propionate D of heme I seems to be the most interesting, generally showing a high correlation not only with its own heme, but also with heme II, corresponding to an indirect stabilization of the reduced forms of both hemes. Other interesting sites are the free histidines of two of the cytochromes and propionate D of heme IV, the latter being potentially associated with redox-induced structural changes. Among the set of cytochromes c3 analyzed in this study, significant differences are observed for several properties of the acidic cytochrome included in the set, from Desulfovibrio africanus, supporting the hypothesis of a different functional role.     

Structural basis for the molecular properties of cytochrome c6

This is a thorough biochemical, spectroscopic, electrochemical, and structural study of a cytochrome c(6) isolated from the filamentous green alga Cladophora glomerata. The protein sequence, elucidated using chemical and mass spectrometric techniques, features 91 amino acids and the characteristic CXXCH heme-binding motif found in c-type cytochromes. The protein is monomeric in both oxidation forms, thereby putting in question a functional role for protein dimerization. Direct electrochemical measurements established, for the first time, the kinetic and thermodynamic data for the redox process in a cytochrome c(6). In particular, the quasi-reversible and diffusion-controlled redox process is accompanied by negative enthalpy and entropy changes, resulting in an E degrees ' value of 0.352 V at 298 K. The pH-dependent properties of the oxidized protein, detected by UV-visible, NMR, and direct cyclic voltammetry, indicate the presence of two acid-base equilibria occurring in the acidic (pK(a) = 4.5) and alkaline regions (pK(a) = 9.0). NMR and electronic spectra allowed the assignment of these equilibria to deprotonation of heme propionate-7 and to replacement of the axial methionine with another ligand, respectively. The 1.3 A resolution X-ray structure of the oxidized protein, revealing a fold typical for class I cytochromes, suggests that the conserved Lys60 replaces the axial methionine at pH >9. The heme solvent accessibility is low, and no water molecules were found in the vicinity of the axial ligands of the heme Fe. A structure-based alignment of cytochromes c(6), and the direct comparison of their structures, indicate a substantial degree of identity between the tertiary structures and suggest patches involved in protein-protein interaction. In particular, the surface electrostatic potential of cytochromes c(6) features a hydrophobic region around the heme cofactor, and a backside surface rich in negative charges.     

Thermodynamic characterization of triheme cytochrome PpcA from Geobacter sulfurreducens: evidence for a role played in e-/H+ energy transduction

The facultative aerobic bacterium Geobacter sulfurreducens produces a small periplasmic c-type triheme cytochrome with 71 residues (PpcA) under anaerobic growth conditions, which is involved in the iron respiration. The thermodynamic properties of the PpcA redox centers and of a protonatable center were determined using NMR and visible spectroscopy techniques. The redox centers have negative and different reduction potentials (-162, -143, and -133 mV for heme I, III, and IV, respectively, for the fully reduced and protonated protein), which are modulated by redox interactions among the hemes (covering a range from 10 to 36 mV) and by redox-Bohr interactions (up to -62 mV) between the hemes and a protonatable center located in the proximity of heme IV. All the interactions between the four centers are dominated by electrostatic effects. The microscopic reduction potential of heme III is the one most affected by the oxidation of the other hemes, whereas heme IV is the most affected by the protonation state of the molecule. The thermodynamic properties of PpcA showed that pH strongly modulates the redox behavior of the individual heme groups. A preferred electron transfer pathway at physiologic pH is defined, showing that PpcA has the necessary thermodynamic properties to perform e-/H+ energy transduction, contributing to a H+ electrochemical potential gradient across the periplasmic membrane that drives ATP synthesis. PpcA is 46% identical in sequence to and shares a high degree of structural similarity with a periplasmic triheme cytochrome c7 isolated from Desulfuromonas acetoxidans, a bacterium closely related to the Geobacteracea family. However, the results obtained for PpcA are quite different from those published for D. acetoxidans c7, and the physiological consequences of these differences are discussed.     

The tetraheme cytochrome from Shewanella oneidensis MR-1 shows thermodynamic bias for functional specificity of the hemes

Bacteria of the genus Shewanella contain an abundant small tetraheme cytochrome in their periplasm when growing anaerobically. Data collected for the protein isolated from S. oneidensis MR-1 and S. frigidimarina indicate differences in the order of oxidation of the hemes. A detailed thermodynamic characterization of the cytochrome from S. oneidensis MR-1 in the physiological pH range was performed, with data collected in the pH range 5.5–9.0 from NMR experiments using partially oxidized samples and from redox titrations followed by visible spectroscopy. These data allow the parsing of the redox and redox–protonation interactions that occur during the titration of hemes. The results show that electrostatic effects dominate the heme–heme interactions, in agreement with modest redox-linked structural modifications, and protonation has a considerable influence on the redox properties of the hemes in the physiological pH range. Theoretical calculations using the oxidized and reduced structures of this protein reveal that the bulk redox–Bohr effect arises from the aggregate fractional titration of several of the heme propionates. This detailed characterization of the thermodynamic properties of the cytochrome shows that only a few of the multiple microscopic redox states that the protein can access are significantly populated at physiological pH. On this basis a functional pathway for the redox activity of the small tetraheme cytochrome from S. oneidensis MR-1 is proposed, where reduction and protonation are thermodynamically coupled in the physiological range. The differences between the small tetraheme cytochromes from the two organisms are discussed in the context of their biological role.     

Proton-assisted two-electron transfer in natural variants of tetraheme cytochromes from Desulfomicrobium Sp

The tetraheme cytochrome c3 isolated from Desulfomicrobium baculatum (DSM 1743)(Dsmb) was cloned, and the sequence analysis showed that this cytochrome differs in just three amino acid residues from the cytochrome c3 isolated from Desulfomicrobium norvegicum (Dsmn): (DsmnXXDsmb) Thr-37 --> Ser, Val-45 --> Ala, and Phe-88 --> Tyr. X-ray crystallography was used to determine the structure of cytochrome c3 from Dsmb, showing that it is very similar to the published structure of cytochrome c3 from Dsmn. A detailed thermodynamic and kinetic characterization of these two tetraheme cytochromes c3 was performed by using NMR and visible spectroscopy. The results obtained show that the network of cooperativities between the redox and protonic centers is consistent with a synergetic process to stimulate the hydrogen uptake activity of hydrogenase. This is achieved by increasing the affinity of the cytochrome for protons through binding electrons and, reciprocally, by favoring a concerted two-electron transfer assisted by the binding of proton(s). The data were analyzed within the framework of the differences in the primary and tertiary structures of the two proteins, showing that residue 88, close to heme I, is the main cause for the differences in the microscopic thermodynamic parameters obtained for these two cytochromes c3. This comparison reveals how replacement of a single amino acid can tune the functional properties of energy-transducing proteins, so that they can be optimized to suit the bioenergetic constraints of specific habitats.     

Redox-coupled conformational alternations in cytochrome c(3) from D. vulgaris Miyazaki F on the basis of its reduced solution structure

Heteronuclear NMR spectroscopy was performed to determine the solution structure of (15)N-labeled ferrocytochrome c(3) from Desulfovibrio vulgaris Miyazaki F (DvMF). Although the folding of the reduced cytochrome c(3) in solution was similar to that of the oxidized one in the crystal structure, the region involving hemes 1 and 2 was different. The redox-coupled conformational change is consistent with the reported solution structure of D. vulgaris Hildenborough ferrocytochrome c(3), but is different from those of other cytochromes c(3). The former is homologous with DvMF cytochrome c(3) in amino acid sequence. Small displacements of hemes 1 and 2 relative to hemes 3 and 4 were observed. This observation is consistent with the unusual behavior of the 2(1)CH(3) signal of heme 3 reported previously. As shown by the (15)N relaxation parameters of the backbone, a region between hemes 1 and 2 has more flexibility than the other regions. The results of this work strongly suggest that the cooperative reduction of hemes 1 and 2 is based on the conformational changes of the C-13 propionate of heme 1 and the aromatic ring of Tyr43, and the interaction between His34 and His 35 through covalent and coordination bonds.     

Molecular basis for redox-Bohr and cooperative effects in cytochrome c3 from Desulfovibrio desulfuricans ATCC 27774: crystallographic and modeling studies of oxidized and reduced high-resolution structures at pH 7.6

The tetraheme cytochrome c3 is a small metalloprotein with ca. 13,000 Da found in sulfate-reducing bacteria, which is believed to act as a partner of hydrogenase. The three-dimensional structure of the oxidized and reduced forms of cytochrome c3 from Desulfovibrio desulfuricans ATCC 27774 at pH 7.6 were determined using high-resolution X-ray crystallography and were compared with the previously determined oxidized form at pH 4.0. Theoretical calculations were performed with both structures, using continuum electrostatic calculations and Monte Carlo sampling of protonation and redox states, in order to understand the molecular basis of the redox-Bohr and cooperativity effects related to the coupled transfer of electrons and protons. We were able to identify groups that showed redox-linked conformational changes. In particular, Glu61, His76, and propionate D of heme II showed important contributions to the redox-cooperativity, whereas His76, propionate A of heme I, and propionate D of heme IV were the key residues for the redox-Bohr effect. Upon reduction, an important movement of the backbone region surrounding hemes I and II was also identified, that, together with a few redox-linked conformational changes in side-chain residues, results in a significant decrease in the solvent accessibility of hemes I and II.     

Resilience of Rhodobacter sphaeroides cytochrome bc1 to heme c1 ligation changes

Typically, c hemes are bound to the protein through two thioether bonds to cysteines and two axial ligands to the heme iron. In high-potential class I c-type cytochromes, these axial ligands are commonly His-Met. A change in this methionine axial ligand is often correlated with a dramatic drop in the heme redox potential and loss of function. Here we describe a bacterial cytochrome c with an unusual tolerance to the alternations in the heme ligation pattern. Substitution of the heme ligating methionine (M185) in cytochrome c1 of the Rhodobacter sphaeroides cytochrome bc1 complex with Lys and Leu lowers the redox midpoint potential but not enough to prevent physiologically competent electron transfer in these fully functional variants. Only when Met-185 is replaced with His is the drop in the redox potential sufficiently large to cause cytochrome bc1 electron transfer chain failure. Functional mutants preserve the structural integrity of the heme crevice: only the nonfunctional His variant allows carbon monoxide to bind to reduced heme, indicating a significant opening of the heme environment. This range of cytochrome c1 ligand mutants exposes both the relative resilience to sixth axial ligand change and the ultimate thermodynamic limits of operation of the cofactor chains in cytochrome bc1.     

Cytochrome c-554 from Methylosinus trichosporium OB3b; a protein that belongs to the cytochrome c2 family and exhibits a HALS-Type EPR signal

A small soluble cytochrome c-554 purified from Methylosinus trichosporium OB3b has been purified and analyzed by amino acid sequencing, mass spectrometry, visible, CD and EPR spectroscopies. It is found to be a mono heme protein with a characteristic cytochrome c fold, thus fitting into the class of cytochrome c(2), which is the bacterial homologue of mitochondrial cytochrome c. The heme iron has a Histidine/Methionine axial ligation and exhibits a highly anisotropic/axial low spin (HALS) EPR signal, with a g(max) at 3.40, and ligand field parameters V/ξ = 0.99, Δ/ξ = 4.57. This gives the rhombicity V/Δ = 0.22. The structural basis for this HALS EPR signal in Histidine/Methionine ligated hemes is not resolved. The ligand field parameters observed for cytochrome c-554 fits the observed pattern for other cytochromes with similar ligation and EPR behaviour.     

Thermodynamic and kinetic characterization of trihaem cytochrome c3 from Desulfuromonas acetoxidans.

Trihaem cytochrome c3 (also known as cytochrome c551.5 and cytochrome c7) is isolated from the periplasmic space of Desulfuromonas acetoxidans, a sulfur-reducing bacterium. Thermodynamic and kinetic data for the trihaem cytochrome c3 are presented and discussed in the context of the possible physiological implications of its functional properties with respect to the natural habitat of D. acetoxidans, namely as a symbiont with green sulfur bacteria working as a mini-sulfuretum. The thermodynamic properties were determined through the fit of redox titration data, followed by NMR and visible spectroscopy, to a model of four functional centres that describes the network of cooperativities between the three haems and one protolytic centre. The kinetics of trihaem cytochrome c3 reduction by sodium dithionite were studied using the stopped-flow technique and the data were fitted to a kinetic model that makes use of the thermodynamic properties to obtain the rate constants of the individual haems. This analysis indicates that the electrons enter the cytochrome mainly via haem I. The reduction potentials of the haems in this cytochrome show little variation with pH within the physiological range, and the kinetic studies show that the rates of reduction are also independent of pH in the range studied. Thus, although the trihaem cytochrome c3 is readily reduced by hydrogenases from Desulfovibrio sp. and its haem core is similar to that of the homologous tetrahaem cytochromes c3, its physico-chemical properties are quite different, which suggests that these multihaem cytochromes with similar structures perform different functions.     

Ferredoxin electron transfer site on cytochrome c3. Structural hypothesis of an intramolecular electron transfer pathway within a tetra-heme cytochrome.

To specify electron exchanges involving Desulfovibrio desulfuricans Norway tetra-heme cytochrome c3, the chemical modification of arginine 73 residue, was performed. Biochemical and biophysical studies have shown that the modified cytochrome retains its ability to both interact and act as an electron carrier with its redox partners, ferredoxin and hydrogenase. Moreover, the chemical modification effects on the cytochrome c3 1H NMR spectrum were similar to that induced by the presence of ferredoxin. This suggests that arginine 73 is localized on the cytochrome c3 ferredoxin interacting site. The identification of heme 4, the closest heme to arginine 73, as the ferredoxin interacting heme helps us to hypothesize about the role of the three other hemes in the molecule. A structural hypothesis for an intramolecular electron transfer pathway, involving hemes 4, 3 and 1, is proposed on the basis of the crystal structures of D. vulgaris Miyazaki and D. desulfuricans Norway cytochromes c3. The unique role of some structural features (alpha helix, aromatic residues) intervening between the heme groups, is proposed.     

Crystal structure of the oxidised and reduced acidic cytochrome c3from Desulfovibrio africanus.

Unique among sulphate-reducing bacteria, Desulfovibrio africanus has two periplasmic tetraheme cytochromes c3, one with an acidic isoelectric point which exhibits an unusually low reactivity towards hydrogenase, and another with a basic isoelectric point which shows the usual cytochrome c3reactivity. The crystal structure of the oxidised acidic cytochrome c3of Desulfovibrio africanus (Dva.a) was solved by the multiple anomalous diffraction (MAD) method and refined to 1.6 A resolution. Its structure clearly belongs to the same family as the other known cytochromes c3, but with weak parentage with those of the Desulfovibrio genus and slightly closer to the cytochromes c3of Desulfomicrobium norvegicum. In Dva.a, one edge of heme I is completely exposed to the solvent and surrounded by a negatively charged protein surface. Heme I thus seems to play an important role in electron exchange, in addition to heme III or heme IV which are the electron exchange ports in the other cytochromes c3. The function of Dva.a and the nature of its redox partners in the cell are thus very likely different.By alignment of the seven known 3D structures including Dva.a, it is shown that the structure which is most conserved in all cytochromes c3is the four-heme cluster itself. There is no conserved continuous protein structure which could explain the remarkable invariance of the four-heme cluster. On the contrary, the proximity of the heme edges is such that they interact directly by hydrophobic and van der Waals contacts. This direct interaction, which always involves a pyrrole CA-CB side-chain and its bound protein cysteine Sgammaatom, is probably the main origin of the four-heme cluster stability. The same kind of interaction is found in the chaining of the hemes in other multihemic redox proteins.The crystal structure of reduced Dva. a was solved at 1.9 A resolution. The comparison of the oxidised and reduced structures reveals changes in the positions of water molecules and polar residues which probably result from changes in the protonation state of amino acids and heme propionates. Water molecules are found closer to the hemes and to the iron atoms in the reduced than in the oxidised state. A global movement of a chain fragment in the vicinity of hemes III and IV is observed which result very likely from the electrostatic reorganization of the polypeptide chain induced by reduction.     

Cytochrome c3 from Desulfovibrio gigas: crystal structure at 1.8 A resolution and evidence for a specific calcium-binding site.

Crystals of the tetraheme cytochrome c3 from sulfate-reducing bacteria Desulfovibrio gigas (Dg) (MW 13 kDa, 111 residues, four heme groups) were obtained and X-ray diffraction data collected to 1.8 A resolution. The structure was solved by the method of molecular replacement and the resulting model refined to a conventional R-factor of 14.9%. The three-dimensional structure shows many similarities to other known crystal structures of tetraheme c3 cytochromes, but it also shows some remarkable differences. In particular, the location of the aromatic residues around the heme groups, which may play a fundamental role in the electron transfer processes of the molecule, are well conserved in the cases of hemes I, III, and IV. However, heme II has an aromatic environment that is completely different to that found in other related cytochromes c3. Another unusual feature is the presence of a Ca2+ ion coordinated by oxygen atoms supplied by the protein within a loop near the N-terminus. It is speculated that this loop may be stabilized by the presence of this Ca2+ ion, may contribute to heme-redox perturbation, and might even be involved in the specificity of recognition with its redox partner.     

Equilibrium unfolding of a small low-potential cytochrome, cytochrome c553 from Desulfovibrio vulgaris.

To understand general aspects of stability and folding of c-type cytochromes, we have studied the folding characteristics of cytochrome c553 from Desulfovibrio vulgaris (Hildenborough). This cytochrome is structurally similar but lacks sequence homology to other heme proteins; moreover, it has an abnormally low reduction potential. Unfolding of oxidized and reduced cytochrome c553 by guanidine hydrochloride (GuHCl) was monitored by circular dichroism (CD) and Soret absorption; the same unfolding curves were obtained with both methods supporting that cytochrome c553 unfolds by an apparent two-state process. Reduced cytochrome c553 is 7(3) kJ/mol more stable than the oxidized form; accordingly, the reduction potential of unfolded cytochrome c553 is 100(20) mV more negative than that of the folded protein. In contrast to many other unfolded cytochrome c proteins, upon unfolding at pH 7.0 both oxidized and reduced heme in cytochrome c553 become high-spin. The lack of heme misligation in unfolded cytochrome c553 implies that its unfolded structure is less constrained than those of cytochromes c with low-spin, misligated hemes.     

EPR determination of interaction redox potentials in a multiheme cytochrome: cytochrome c3 from Desulfovibrio desulfuricans Norway

In cytochromes c3 which contain four hemes per molecule, the redox properties of each heme may depend upon the redox state of the others. This effect can be described in terms of interaction redox potentials between the hemes and must be taken into account in the characterization of the redox properties of the molecule. We present here a method of measurement of these interactions based on the EPR study of the redox equilibria of the protein. The microscopic and macroscopic midpoint potentials and the interaction potentials are deduced from the analysis of the redox titration curves of the intensity and the amplitude of the EPR spectrum. This analysis includes a precise simulation of the spectrum of the protein in the oxidized state in order to determine the relative contribution of each heme to the spectral amplitude. Using our method on cytochrome c3 from D. desulfuricans Norway, we found evidence for the existence of weak interaction potentials between the hemes. The three interaction potentials which have been measured are characterized by absolute values lower than 20 mV in contrast with the values larger than 40-50 mV which have been reported for cytochrome c3 from D. gigas. Simulations of the spectra of samples poised at different potentials indicate a structural modification of the heme with the most negative potential during the first step of reduction. The correspondence between the redox sites as characterized by the EPR potentiometric titration and the hemes in the tridimensional structure is discussed.     

Two types of cytochrome cd1 in the aerobic photosynthetic bacterium, Erythrobacter sp. OCh 114

Components I and II of cytochrome cd1 which had different spectral features were purified from the aerobic photosynthetic bacterium, Erythrobacter sp. strain OCh 114. Component I showed an absorption maxima at 700 and 406 nm in the oxidized form, and at 621, 552.5, 548 and 416 nm in the reduced form. Component II showed an absorption maxima at 635 and 410 nm in the oxidized form and at 628, 552.5, 548 and 417 nm in the reduced form. The relative molecular mass, Mr, of both cytochromes was determined to be 135,000 with two identical subunits. Components I and II showed pI values of 7.6 and 6.8, respectively. The redox potential of hemes ranged from +234 mV to +242 mV, except for the heme d1 of component I (Em7 = +134 mV). Components I and II showed both cytochrome c oxidase and nitrite reductase activities. Cytochrome c oxidase activity was strongly inhibited by a low concentration of nitrite and cyanide. Erythrobacter cytochromes c-551 and c-552 were utilized as electron donors for the cytochrome c oxidase reaction. The high affinity of cytochrome c-552 to component II (Km = 1.27 microM) suggested a physiological significance for this cytochrome. Erythrobacter cytochromes cd1 are unique in their presence in cells grown under aerobic conditions as compared to other bacterial cytochromes cd1 which are formed only under denitrifying conditions.     

Binding of ligands originates small perturbations on the microscopic thermodynamic properties of a multicentre redox protein

NMR and visible spectroscopy coupled to redox measurements were used to determine the equilibrium thermodynamic properties of the four haems in cytochrome c3 under conditions in which the protein was bound to ligands, the small anion phosphate and the protein rubredoxin with the iron in the active site replaced by zinc. Comparison of these results with data for the isolated cytochrome shows that binding of ligands causes only small changes in the reduction potentials of the haems and their pairwise interactions, and also that the redox-sensitive acid-base centre responsible for the redox-Bohr effect is essentially unaffected. Although neither of the ligands tested is a physiological partner of cytochrome c3, the small changes observed for the thermodynamic properties of cytochrome c3 bound to these ligands vs. the unbound state, indicate that the thermodynamic properties measured for the isolated protein are relevant for a physiological interpretation of the role of this cytochrome in the bioenergetic metabolism of Desulfovibrio.     

Ionic Strength-Dependent Physicochemical Factors in Cytochrome c3 Regulating the Electron Transfer Rate

The effect of ionic strength on the macroscopic and microscopic redox potentials and the heme environment of cytochrome c₃ from Desulfovibrio vulgaris Miyazaki F have been investigated by NMR and electrochemical methods. The redox potentials of this tetraheme protein are found to be ionic strength-dependent. Especially, the microscopic redox potentials of hemes 2 and 3 at the fourth reduction step increase significantly with increasing ionic strength, which is in contradiction to the theoretical expectation. The coordinated imidazole proton signals are unaffected by ionic strength. However, the methyl and propionate proton signals of hemes 1 and 4 showed significant ionic strength dependencies that are distinct from those for hemes 2 and 3. This heme classification is the same as that found in the ionic strength dependencies of the microscopic redox potentials at the fourth reduction step. Furthermore, the effect of ionic strength on the electrostatic potentials at the heme irons has been examined on the theoretical basis. The electrostatic potential at heme 4 changes up to 1 M ionic strength, which was not expected from the observations reported on cytochromes so far. These results are discussed in connection with the reported anomalous ionic strength dependency of the reduction rate of cytochrome c₃.