Electron transport in sulfate-reducing bacteria. Molecular modeling and NMR studies of the rubredoxin--tetraheme-cytochrome-c3 complex

Transcriptome-wide marker gene expression analysis of stress-responsive sulfate-reducing bacteria

Sulfate-reducing bacteria (SRB) are terminal members of any anaerobic food chain. For example, they critically influence the biogeochemical cycling of carbon, nitrogen, sulfur, and metals (natural environment) as well as the corrosion of civil infrastructure (built environment). The United States alone spends nearly $4 billion to address the biocorrosion challenges of SRB. It is important to analyze the genetic mechanisms of these organisms under environmental stresses. The current study uses complementary methodologies, viz., transcriptome-wide marker gene panel mapping and gene clustering analysis to decipher the stress mechanisms in four SRB. Here, the accessible RNA-sequencing data from the public domains were mined to identify the key transcriptional signatures. Crucial transcriptional candidate genes of Desulfovibrio spp. were accomplished and validated the gene cluster prediction. In addition, the unique transcriptional signatures of Oleidesulfovibrio alaskensis (OA-G20) at graphene and copper interfaces were discussed using in-house RNA-sequencing data. Furthermore, the comparative genomic analysis revealed 12,821 genes with translation, among which 10,178 genes were in homolog families and 2643 genes were in singleton families were observed among the 4 genomes studied. The current study paves a path for developing predictive deep learning tools for interpretable and mechanistic learning analysis of the SRB gene regulation.

Rubredoxin as a paramagnetic relaxation-inducing probe

The paramagnetic effect due to the presence of a metal center with unpaired electrons is no longer considered a hindrance in protein NMR spectroscopy. In the present work, the paramagnetic effect due to the presence of a metal center with unpaired electrons was used to map the interface of an electron transfer complex. Desulfovibrio gigas cytochrome c(3) was chosen as target to study the effect of the paramagnetic probe, Fe-rubredoxin, which produced specific line broadening in the heme IV methyl resonances M2(1) and M18(1). The rubredoxin binding surface in the complex with cytochrome c(3) was identified in a heteronuclear 2D NMR titration. The identified heme methyls on cytochrome c(3) are involved in the binding interface of the complex, a result that is in agreement with the predicted complexes obtained by restrained molecular docking, which shows a cluster of possible solutions near heme IV. The use of a paramagnetic probe in (1)HNMR titration and the mapping of the complex interface, in combination with a molecular simulation algorithm proved to be a valuable strategy to study electron transfer complexes involving non-heme iron proteins and cytochromes.

Redox Interaction of Cytochrome c3 with [NiFe] Hydrogenase from Desulfovibrio vulgaris Miyazaki F

Cytochrome c3 isolated from a sulfate-reducing bacterium, Desulfovibrio vulgaris Miyazaki F, is a tetraheme protein. Its physiological partner, [NiFe] hydrogenase, catalyzes the reversible oxidoreduction of molecular hydrogen. To elucidate the mechanism of electron transfer between cytochrome c3 and [NiFe] hydrogenase, the transient complex formation by these proteins was investigated by means of NMR. All NH signals of uniformly 15N-labeled ferric cytochrome c3 except N-terminus, Pro, and Gly73 were assigned. 1H-15N HSQC spectra were recorded for 15N-labeled ferric and ferrous cytochrome c3, in the absence and presence of hydrogenase. Chemical shift perturbations were observed in the region around heme 4 in both oxidation states. Additionally, the region between hemes 1 and 3 in ferrous cytochrome c3 was affected in the presence of hydrogenase, suggesting that the mode of interaction is different in each redox state. Heme 3 is probably the electron gate for ferrous cytochrome c3. To investigate the transient complex of cytochrome c3 and hydrogenase in detail, modeling of the complex was performed for the oxidized proteins using a docking program, ZDOCK 2.3, and NMR data. Furthermore, the roles of lysine residues of cytochrome c3 in the interaction with hydrogenase were investigated by site-directed mutagenesis. When the lysine residues around heme 4 were replaced by an uncharged residue, methionine, one by one, the Km of the electron-transfer kinetics increased. The results showed that the positive charges of Lys60, Lys72, Lys95, and Lys101 around heme 4 are important for formation of the transient complex with [NiFe] hydrogenase in the initial stage of the cytochrome c3 reduction. This finding is consistent with the most possible structure of the transient complex obtained by modeling.

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.

Transient complexes of redox proteins: structural and dynamic details from NMR studies

Redox proteins participate in many metabolic routes, in particular those related to energy conversion. Protein-protein complexes of redox proteins are characterized by a weak affinity and a short lifetime. Two-dimensional NMR spectroscopy has been applied to many redox protein complexes, providing a wealth of information about the process of complex formation, the nature of the interface and the dynamic properties of the complex. These studies have shown that some complexes are non-specific and exist as a dynamic ensemble of orientations while in other complexes the proteins assume a single orientation. The binding interface in these complexes consists of a small hydrophobic patch for specificity, surrounded by polar, uncharged residues that may enhance dissociation, and, in most complexes, a ring or patch of charged residues that enhances the association by electrostatic interactions. The entry and exit port of the electrons is located within the hydrophobic interaction site, ensuring rapid electron transfer from one redox centre to the next.

The cytochrome c3-[Fe]-hydrogenase electron-transfer complex: structural model by NMR restrained docking

Cytochrome c(3) (M(r) 13000) is a low redox potential cytochrome specific of the anaerobic metabolism in sulfate-reducing bacteria. This tetrahemic cytochrome is an intermediate between the [Fe]-hydrogenase and the cytochrome Hmc in Desulfovibrio vulgaris Hildenborough strain. The present work describes the structural model of the cytochrome c(3)-[Fe]-hydrogenase complex obtained by nuclear magnetic resonance restrained docking. This model connects the distal cluster of the [Fe]-hydrogenase to heme 4 of the cytochrome, the same heme found in the interaction with cytochrome Hmc. This result gives evidence that cytochrome c(3) is an electron shuttle between the periplasmic hydrogenase and the Hmc membrane-bound complex.

Three-dimensional structure of the nonaheme cytochrome c from Desulfovibrio desulfuricans Essex in the Fe(III) state at 1.89 A resolution

A nine heme group containing cytochrome c isolated from the soluble and membrane fractions of Desulfovibrio desulfuricans Essex, termed nonaheme cytochrome c, was crystallized, and the structure was solved using the multiple wavelength anomalous dispersion (MAD) phasing method. Refinement was carried out to a resolution of 1.89 A, and anisotropic temperature factors were addressed to the iron and sulfur atoms in the model. The structure revealed two cytochrome c(3) like domains with the typical arrangement of four heme centers. Both domains flanked an extra heme buried under the protein surface. This heme is held in position by loop extensions in each of the two domains. Although both the N- and C-terminal tetraheme domains exhibit a fold and heme arrangement very similar to that of cytochrome c(3), they differ considerably in their loop extensions and electrostatic surface. Analysis of the structure provides evidence for a different function of both domains, namely, anchoring the protein in a transmembranous complex with the N-terminal domain and formation of an electron-transfer complex with hydrogenase by the C-terminal domain.

Solution structure of a zinc substituted eukaryotic rubredoxin from the cryptomonad alga Guillardia theta

The rubredoxin from the cryptomonad Guillardia theta is one of the first examples of a rubredoxin encoded in a eukaryotic organism. The structure of a soluble zinc-substituted 70-residue G. theta rubredoxin lacking the membrane anchor and the thylakoid targeting sequence was determined by multidimensional heteronuclear NMR, representing the first three-dimensional (3D) structure of a eukaryotic rubredoxin. For the structure calculation a strategy was applied in which information about hydrogen bonds was directly inferred from a long-range HNCO experiment, and the dynamics of the protein was deduced from heteronuclear nuclear Overhauser effect data and exchange rates of the amide protons. The structure is well defined, exhibiting average root-mean-square deviations of 0.21 A for the backbone heavy atoms and 0.67 A for all heavy atoms of residues 7-56, and an increased flexibility toward the termini. The structure of this core fold is almost identical to that of prokaryotic rubredoxins. There are, however, significant differences with respect to the charge distribution at the protein surface, suggesting that G. theta rubredoxin exerts a different physiological function compared to the structurally characterized prokaryotic rubredoxins. The amino-terminal residues containing the putative signal peptidase recognition/cleavage site show an increased flexibility compared to the core fold, but still adopt a defined 3D orientation, which is mainly stabilized by nonlocal interactions to residues of the carboxy-terminal region. This orientation might reflect the structural elements and charge pattern necessary for correct signal peptidase recognition of the G. theta rubredoxin precursor.

Cross-linking between cytochrome c3 and flavodoxin from Desulfovibrio gigas

Tetraheme cytochrome c3 (13 kDa) and flavodoxin (16 kDa), are small electron transfer proteins that have been used to mimic, in vitro, part of the electron-transfer chain that operates between substract electron donors and respiratory electron acceptors partners in Desulfovibrio species (Palma, N., Moura, I., LeGall, J., Van Beeumen, J., Wampler, J., Moura, J. J. G. (1994) Biochemistry 33, 6394-6407). The electron transfer between these two proteins is believed to occur through the formation of a specific complex where electrostatic interaction is the main driving force (Stewart, D., LeGall, J., Moura, I., Moura, J.J.G., Peck, H.D., Xavier, A.V., Weiner, P.K. and Wampler, J.E. (1988) Biochemistry 27, 2444-2450, Stewart, D., LeGall, J., Moura, I., Moura, J.J.G., Peck, H.D., Xavier, A.V., Weiner, P., Wampler, J. (1989) Eur. J. Biochem. 185, 695-700). In order to obtain structural information of the pre-complex, a covalent complex between the two proteins was prepared. A water-soluble carbodiimide [EDC (1-ethyl-3(3 dimethylaminopropyl) carbodiimide hydrochloride] was used for the cross linking reaction. The reaction was optimized varying a wide number of experimental parameters such as ionic strength, protein and cross linker concentration, and utilization of different cross linkers and reaction time between the crosslinker and proteins.

The primary and three-dimensional structures of a nine-haem cytochrome c from Desulfovibrio desulfuricans ATCC 27774 reveal a new member of the Hmc family

BACKGROUND

Haem-containing proteins are directly involved in electron transfer as well as in enzymatic functions. The nine-haem cytochrome c (9Hcc), previously described as having 12 haem groups, was isolated from cells of Desulfovibrio desulfuricans ATCC 27774, grown under both nitrate- and sulphate-respiring conditions.

RESULTS

Models for the primary and three-dimensional structures of this cytochrome, containing 292 amino acid residues and nine haem groups, were derived using the multiple wavelength anomalous dispersion phasing method and refined using 1.8 A diffraction data to an R value of 17.0%. The nine haem groups are arranged into two tetrahaem clusters, with Fe-Fe distances and local protein fold similar to tetrahaem cytochromes c3, while the extra haem is located asymmetrically between the two clusters.

CONCLUSIONS

This is the first known three-dimensional structure in which multiple copies of a tetrahaem cytochrome c3-like fold are present in the same polypeptide chain. Sequence homology was found between this cytochrome and the C-terminal region (residues 229-514) of the high molecular weight cytochrome c from Desulfovibrio vulgaris Hildenborough (DvH Hmc). A new haem arrangement in domains III and IV of DvH Hmc is proposed. Kinetic experiments showed that 9Hcc can be reduced by the [NiFe] hydrogenase from D. desulfuricans ATCC 27774, but that this reduction is faster in the presence of tetrahaem cytochrome c3. As Hmc has never been found in D. desulfuricans ATCC 27774, we propose that 9Hcc replaces it in this organism and is therefore probably involved in electron transfer across the membrane.

Structural and kinetic studies of the Y73E mutant of octaheme cytochrome c3 (Mr = 26 000) from Desulfovibrio desulfuricans Norway

A combination of structural, kinetic, and interaction experiments has been used to study the role of a highly conserved aromatic residue, Tyr73, parallel to the sixth heme axial ligand of heme 4 in multiheme cytochrome c3 (Mr = 26 000), also called cytochrome cc3 or octaheme cytochrome, from Desulfovibrio desulfuricans Norway. This residue is expected to be involved in intermolecular electron transfer and protein-protein interaction, since heme 4 is described to be the interaction site between physiological partners. The kinetic experiments show that the Y73E replacement provokes no significant change in the electron-transfer reaction with the physiological partner, the [NiFeSe] hydrogenase, but that the protein-protein interaction between cytochrome c3 (Mr = 26 000) and hydrogenase is strongly affected by the mutation. The aromatic residue does not play a role in maintaining the axial heme ligand in a particular orientation, since the mutation did not affect the orientation of histidine 77, the sixth axial ligand of heme 4. The structural analysis by X-ray crystallography clearly shows that a rearrangement of the charged residues in the vicinity of the mutation site is responsible for the change in protein-protein interaction, which is of an electrostatic nature. Lys22 and Arg66, residues which are located at the interacting surface, are twisted toward the mutated position Glu73 in order to compensate for the negative charge and therefore are no longer accessible for the docking with a physiological partner. Tyr73 has instead a structural function and probably a role in maintaining the hydrophobic environment of the heme 4 cavity rather than a function in the intermolecular electron transfer with the physiological partners.

Further characterization of the two tetraheme cytochromes c3 from Desulfovibiro africanus: nucleotide sequences, EPR spectroscopy and biological activity

The genes encoding the basic and acidic tetraheme cytochromes c3 from Desulfovibrio africanus have been sequenced. The corresponding amino acid sequences of the basic and acidic cytochromes c3 indicate that the mature proteins consist of a single polypeptide chain of 117 and 103 residues, respectively. Their molecular masses, 15102 and 13742 Da, respectively, determined by mass spectrometry, are in perfect agreement with those calculated from their amino acid sequences. Both D. africanus cytochromes c3 are synthesized as precursor proteins with signal peptides of 23 and 24 residues for the basic and acidic cytochromes, respectively. These cytochromes c3 exhibit the main structural features of the cytochrome c3 family and contain the 16 strictly conserved cysteine + histidine residues directly involved in the heme binding sites. The D. africanus acidic cytochrome c3 differs from all the other homologous cytochromes by its low content of basic residues and its distribution of charged residues in the amino acid sequence. The presence of four hemes per molecule was confirmed by EPR spectroscopy in both cytochromes c3. The g-value analysis suggests that in both cytochromes, the angle between imidazole planes of the axial histidine ligands is close to 90 degrees for one heme and much lower for the three others. Moreover, an unusually high exchange interaction (approximately 10[-2] cm[-1]) was evidenced between the highest potential heme (-90 mV) and one of the low potential hemes in the basic cytochrome c3. The reactivity of D. africanus cytochromes c3 with heterologous [NiFe] and [Fe] hydrogenases was investigated. Only the basic one interacts with the two types of hydrogenase to achieve efficient electron transfer, whereas the acidic cytochrome c3 exchanges electrons specifically with the basic cytochrome c3. The difference in the specificity of the two D. africanus cytochromes c3 has been correlated with their highly different content of basic and acidic residues.

Crystal structure of a dimeric octaheme cytochrome c3 (M(r) 26,000) from Desulfovibrio desulfuricans Norway

BACKGROUND

The octaheme cytochrome C3 (M(r) 26,000; cc3) from Desulfovibrio desulfuricans Norway is a dimeric cytochrome made up of two identical subunits, each containing four heme groups. It is involved in the redox transfer chain of sulfate-reducing bacteria, which links the periplasmic oxidation of hydrogen to the cytoplasmic reduction of sulfate. The amino-acid sequence of cc3 shows similarities to that of the tetraheme cytochrome c3 (M(r) 13,000; c3) from the same bacteria. Structural analysis of cc3 forms a basis for understanding the precise roles of the multiheme-containing redox proteins and the reason for the presence of several different multiheme cytochromes in one bacterial strain.

RESULTS

The crystal structure of cytochrome cc3 has been determined at 2.16 A resolution. The subunits display the c3 structural fold with significant amino-acid substitutions, relative to the tetraheme cytochromes c3, in the regions of the dimer interface. The identical subunits are related by a crystallographic twofold axis, with one heme of each subunit in close contact. The overall structure and the environments of the different heme groups are compared with those of the tetraheme cytochromes c3.

CONCLUSIONS

A common scheme for interactions between these types of cytochrome and their redox partners involves the interaction of a heme crevice, surrounded by positively charged lysine residues, with acidic residues surrounding the redox partner's functional group. Despite the relatively acidic character of cytochrome cc3, the crevice of one heme is surrounded by a high number of positively charged residues, in the same manner as has been reported for cytochromes c3. The environment of this heme is formed by four flexible surface loops which are variable in length and orientation in the different c3-type cytochromes although the overall structural folds are very similar. It has been proposed that this region, adapted in topology and charge, is the interaction site for physiological partners and is also most likely to be the interaction site in the dimeric cytochrome cc3.

Kinetic studies on the electron-transfer reaction between cytochrome c3 and flavodoxin from Desulfovibrio vulgaris strain Hildenborough

The kinetic properties of the electron-transfer process between reduced Desulfovibrio vulgaris cytochrome c3 and D. vulgaris flavodoxin have been studied by anaerobic stopped-flow techniques. Anaerobic titrations of reduced cytochrome c3 with oxidized flavodoxin show a stoichiometry of 4 mol of flavodoxin required to oxidize the tetraheme cytochrome. Flavodoxin neutral semiquinone and oxidized cytochrome c3 are the only observable products of the reaction. At pH 7.5, the four-electron-transfer reaction is biphasic. Both the rapid and the slow phases exhibit limiting rates as the flavodoxin concentration is increased with respective rates of 73.4 and 18.5 s-1 and respective Kd values of 65.9 +/- 9.4 microM and 54.5 +/- 13 microM. A biphasic electron-transfer rate is observed when the ionic strength is increased to 100 mM KCl; however, the observed rate is no longer saturable, and relative second-order rate constants of 5.3 x 10(5) and 8.5 x 10(4) M-1 s-1 are calculated. The magnitude of the rapid phase of electron transfer diminishes with the level of heme reduction when varying reduced levels of the cytochrome are mixed with oxidized flavodoxin. No rapid phase is observed when 0.66e(-)-reduced cytochrome c3 reacts with an approximately 25-fold molar excess of flavodoxin. At pH 6.0, the electron-transfer reaction is monophasic with a limiting rate of 42 +/- 1.4 s-1 and a Kd value of approximately 8 microM. Increasing the ionic strength of the pH 6.0 solution to 100 microM KCl results in a biphasic reaction with relative second-order rate constants of 5.3 x 10(5) and 1.1 x 10(4) M-1 s-1. Azotobacter vinelandii flavodoxin reacts with reduced D. vulgaris cytochrome c3 in a slow, monophasic manner with limiting rate of electron transfer of 1.2 +/- 0.06 s-1 and a Kd value of 80.9 +/- 10.7 microM. These results are discussed in terms of two equilibrium conformational states for the cytochrome which are dependent on the pH of the medium and the level of heme reduction [Catarino et al. (1991) Eur. J. Biochem. 207, 1107-1113].

Involvement of electrostatic interactions in cytochrome c complex formations

Structural studies on various electron transfer complexes involving the tetrahemic cytochrome c3 provided evidence that one of the hemes (heme 4) is the interacting site on the molecule. The reactivity of this particular heme is allocated to the positive charges found around the heme group which are strongly involved in the electrostatic interaction processes. Electrostatic and hydrophobic effects in complex formation are considered on the basis of two electron transfer complex examples: the soluble cytochrome c-cytochrome c peroxidase and the membrane bound photosynthetic reaction center.

Aldehyde oxidoreductase activity in Desulfovibrio gigas: in vitro reconstitution of an electron-transfer chain from aldehydes to the production of molecular hydrogen

The molybdenum [iron-sulfur] protein, first isolated from Desulfovibrio gigas by Moura et al. [Moura, J. J. G., Xavier, A. V., Bruschi, M., Le Gall, J., Hall, D. O., & Cammack, R. (1976) Biochem. Biophys. Res. Commun. 72, 782-789], was later shown to mediate the electronic flow from salicylaldehyde to a suitable electron acceptor, 2,6-dichlorophenolindophenol (DCPIP) [Turner, N., Barata, B., Bray, R. C., Deistung, J., LeGall, J., & Moura, J. J. G. (1987) Biochem. J. 243, 755-761]. The DCPIP-dependent aldehyde oxidoreductase activity was studied in detail using a wide range of aldehydes and analogues. Steady-state kinetic analysis (KM and Vmax) was performed for acetaldehyde, propionaldehyde, benzaldehyde, and salicylaldehyde in excess DCPIP concentration, and a simple Michaelis-Menten model was shown to be applicable as a first kinetic approach. Xanthine, purine, allopurinol, and N1-methylnicotinamide (NMN) could not be utilized as enzyme substrates. DCPIP and ferricyanide were shown to be capable of cycling the electronic flow, whereas other cation and anion dyes [O2 and NAD(P)+] were not active in this process. The enzyme showed an optimal pH activity profile around 7.8. This molybdenum hydroxylase was shown to be part of an electron-transfer chain comprising four different soluble proteins from D. gigas, with a total of 11 discrete redox centers, which is capable of linking the oxidation of aldehydes to the reduction of protons.

Voltammetric studies of the catalytic electron-transfer process between the Desulfovibrio gigas hydrogenase and small proteins isolated from the same genus

The kinetics of electron transfer between the Desulfovibrio gigas hydrogenase and several electron-transfer proteins from Desulfovibrio species were investigated by cyclic voltammetry, square-wave voltammetry and chronoamperometry. The cytochrome c3 from Desulfovibrio vulgaris (Hildenborough), Desulfovibrio desulfuricans (Norway 4), Desulfovibrio desulfuricans (American Type Culture Collection 27774) and D. gigas (NCIB 9332) were used as redox carriers. They differ in their redox potentials and isoelectric point. Depending on the pH, all the reduced forms of these cytochromes were effective in electron exchange with hydrogenase. Other small electron-transfer proteins such as ferredoxin I, ferredoxin II and rubredoxin from D. gigas were tentatively used as redox carriers. Only ferredoxin II was effective in mediating electron exchange between hydrogenase and the working electrode. The second-order rate constants k for the reaction between reduced proteins and hydrogenase were calculated based on the theory of the simplest electrocatalytic mechanism [Moreno, C., Costa, C., Moura, I., Le Gall, J., Liu, M. Y., Payne, W. J., van Dijk, C. & Moura, J. J. G. (1993) Eur. J. Biochem. 212, 79-86] and the results obtained by cyclic voltammetry were compared with those obtained by chronoamperometry. Values for k of 10(5)-10(6) M-1 s-1 (cytochrome c3 as electron carrier) and 10(4) M-1 s-1 (ferredoxin II as the electron carrier) were determined. The rate-constant values are discussed in terms of the existence of an electrostatic interaction between the electrode surface and the redox carrier and between the redox carrier and a positively charged part of the enzyme.

Kinetic study of the reduction mechanism for Desulfovibrio gigas cytochrome c3

The kinetic aspects of the reduction process in cytochrome c3 from Desulfovibrio gigas have been investigated over a wide range of pH values ranging between pH 5.8 and pH 9.8. The data have been analyzed in the framework of an I2H4 interaction network coupled to a proton-linked equilibrium between two tertiary structures (Cornish-Bowden, A. & Koshland, D.E. Jr (1970) J. Biol. Chem. 245, 6241-6250). The kinetic rate constants for the reduction of the four hemes for the two tertiary conformations have been characterized in the framework of the thermodynamic network obtained from the equilibrium analysis (Coletta, M., Catarino, T., LeGall, J.J. & Xavier, A.V. (1991) Eur. J. Biochem. 202, 1101-1106). The intrinsic reduction rate constants determined by reaction with sodium dithionite for two hemes (namely heme 4 and heme 1) are significantly faster than those for the other two heme residues. In view of the equilibrium redox properties, heme 4 (with the fastest reduction rate) may then work as the kinetic electron-capturing site for the electrons from sodium dithionite. The transfer to hemes 2 and 3 then occurs by virtue of their free-energy levels at equilibrium. At our experimental conditions, there is also transfer of electrons to hemes 2 and 3 from heme 1, which is reduced at a slower rate than heme 4, thus contributing to the biphasic kinetics observed for the overall process. The kinetic parameters obtained are discussed in terms of the mechanism proposed for the coupling between the electron and proton transfer, as induced by the heme/heme cooperativity network.

A thermodynamic model for the cooperative functional properties of the tetraheme cytochrome c3 from Desulfovibrio gigas

A thermodynamic model is presented to describe the redox behaviour of the tetraheme cytochrome c3 from Desulfovibrio gigas. This molecule displays different intrinsic redox potentials for the four hemes and during the redox titration process, interactions among different hemes occur, thus altering the values of redox potentials according to which of the hemes are oxidized [Santos, H., Moura, J.J.G., Moura, I., LeGall, J. & Xavier, A.V. (1984) Eur. J. Biochem. 141, 283-296]. This complex cooperative behaviour [Xavier, A.V. (1986) J. Inorg. Biochem. 28, 239-243] has been analyzed here using an I2H4-interaction network [Cornish-Bowden, A. & Koshland, D.E. Jr (1970) J. Biol. Chem. 245, 6241-6250] coupled to a proton-linked equilibrium between two tertiary structures. Such a formalism, which requires a reduced number of parameters, is able to fully account quantitatively for the pH dependence of the NMR redox-titration curves. The 'redox-Bohr' effect is discussed in terms of the available structure and thermodynamic data and a functional mechanism is proposed.

Molecular dynamics simulations of the cytochrome c3-rubredoxin complex from Desulfovibrio vulgaris

Molecular dynamics simulations have been carried out on the complex formed between the tetraheme cytochrome c3 and the iron protein rubredoxin from the sulfate-reducing bacterium Desulfovibrio vulgaris. These simulations were performed both with explicit solvent water molecules included, and without solvent molecules using a distance-dependent dielectric constant to approximate the screening effects of solvent. The results of both simulations are strikingly different, indicating that the representation of environmental effects is important in such simulations. For example, a striking adaptation of the two proteins seen in the nonsolvated simulation is not seen when explicit solvent water is included; in fact, the complex appears to become weaker in the solvated simulation. Nonetheless, the iron-iron distance decreases more significantly in the solvated simulation than in the nonsolvated simulation. It was found that in both cases molecular dynamics optimized the structures further than energy minimization alone.

Identification of the site of interaction between cytochrome c3 and ferredoxin using peptide mapping of the cross-linked complex

Structural studies carried out on a cross-linked complex between cytochrome c3 and ferredoxin I, both isolated from Desulfovibrio desulfuricans Norway, allowed the identification of the site of interaction between the two redox proteins. Staphylococcus aureus proteinase and chymotrypsin digestions led to characterization of peptides containing both cytochrome c3 and ferredoxin sequences. The cytochrome c3 sequences involved in the three isolated cross-linked peptides contained several lysine residues localized around the heme 4 crevice. This analysis stressed the peculiar role of lysines 100, 101, 103, 104 and 113, which could be considered as major cross-link sites, as opposed to the lysines 75, 79 and 82, which could be considered as minor cross-link sites. One cross-linked peptide, containing two ferredoxin sequences joined to one cytochrome c3 sequence, had been isolated, suggesting the possibility of more than one cross-link per covalent complex. All these results led to the identification of heme 4 of cytochrome c3 as the site of interaction for the ferredoxin I. This study confirms the proposal that could be deduced from the hypothetical structure of the complex built by computer graphics modelling (Cambillau, C., Frey, M., Mosse, J., Guerlesquin, F. and Bruschi, M. (1988) Proteins: struct., funct. genet. 4, 63-70).

Structural and functional approach toward a classification of the complex cytochrome c system found in sulfate-reducing bacteria

Following the discovery of the tetraheme cytochrome c3 in the strict anaerobic sulfate-reducing bacteria (Postgate, J.R. (1954) Biochem. J. 59, xi; Ishimoto et al. (1954) Bull. Chem. Soc. Japan 27, 564-565), a variety of c-type cytochromes (and others) have been reported, indicating that the array of heme proteins in these bacteria is complex. We are proposing here a tentative classification of sulfate- (and sulfur-) reducing bacteria cytochromes c based on: number of hemes per monomer, heme axial ligation, heme spin state and primary structures (whole or fragmentary). Different and complementary spectroscopic tools have been used to reveal the structural features of the heme sites.

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.