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

A brief survey of the "cytochromome"

Multihaem cytochromes c are widespread in nature where they perform numerous roles in diverse anaerobic metabolic pathways. This is achieved in two ways: multihaem cytochromes c display a remarkable diversity of ways to organize multiple hemes within the protein frame; and the hemes possess an intrinsic reactive versatility derived from diverse spin, redox and coordination states. Here we provide a brief survey of multihaem cytochromes c that have been characterized in the context of their metabolic role. The contribution of multihaem cytochromes c to dissimilatory pathways handling metallic minerals, nitrogen compounds, sulfur compounds, organic compounds and phototrophism are described. This aims to set the stage for the further exploration of the vast unknown "cytochromome" that can be anticipated from genomic databases.

Charge parametrization of the DvH-c3 heme group: validation using constant-(pH,E) molecular dynamics simulations

We studied the effect of using different heme group charge parametrization methods and schemes (Merz-Kollman, CHELPG, and single- and multiconformational RESP) on the quality of the results produced by the constant-(pH,E) MD method, applied to the redox titration of Desulfovibrio vulgaris Hildenborough cytochrome c(3). These new and more accurate charge sets enabled us to overcome the previously reported dependence of the method's performance on the dielectric constant, ε, assigned to the protein region. In particular, we found a systematic, clear shift of the E(mod) toward more negative values than those previously reported, in agreement with an electrostatics based reasoning. The simulations showed strong coupling between protonating/redox sites. We were also able to capture significant direct and, especially, indirect interactions between hemes, such as those mediated by histidine 67. Our results highlight the importance of having a good quantum description of the system before deriving atomic partial charges for classic force fields.

Composite films of lecithin and heme proteins with electrochemical and electrocatalytic activities

Functional composite films made from lecithin micelles and the two heme proteins of met-myoglobin (Mb) and met-hemoglobin (Hb) are reported in this paper. Proteins in functional composite films have much higher rates of electron transfer than proteins in solutions on carbon paste (CP) electrodes. Cyclic voltammograms (CVs) all give a pair of well-defined and quasi-reversible peaks, corresponding to the heme FeIII/FeII redox couple of proteins. Differential pulse voltammograms (DPVs) also show the same formal potential (E0') values of proteins under identical conditions. Electronic and vibrational spectra indicate that proteins in these films are not denatured, but their conformational differences from native states may exist. The E0' value for Mb in the lecithin film is found to be pH dependent. The Mb lecithin film can catalytically reduce O2 and H2O2, and its analytical application to H2O2 determination is established.

Photoinduced hydrogen production by direct electron transfer from photosystem I cross-linked with cytochrome c3 to [NiFe]-hydrogenase

The photosynthetic reaction center is an efficient molecular device for the conversion of light energy to chemical energy. In a previous study, we synthesized the hydrogenase/photosystem I (PSI) complex, in which Ralstonia hydrogenase was linked to the cytoplasmic side of Synechocystis PSI, to modify PSI so that it photoproduced molecular hydrogen (H2). In that study, hydrogenase was fused with a PSI subunit, PsaE, and the resulting hydrogenase-PsaE fusion protein was self-assembled with PsaE-free PSI to give the hydrogenase/PSI complex. Although the hydrogenase/PSI complex served as a direct light-to-H2 conversion system in vitro, the activity was totally suppressed by adding physiological PSI partners, ferredoxin (Fd) and ferredoxin-NADP+-reductase (FNR). In the present study, to establish an H2 photoproduction system in which the activity is not interrupted by Fd and FNR, position 40 of PsaE from Synechocystis sp. PCC6803, corresponding to the Fd-binding site on PSI, was selected and targeted for the cross-linking with cytochrome c3 (cytc3) from Desulfovibrio vulgaris. The covalent adduct of cytc3 and PsaE was stoichiometrically assembled with PsaE-free PSI to form the cytc3/PSI complex. The NADPH production by the cytc3/PSI complex coupled with Fd and FNR decreased to approximately 20% of the original activity, whereas the H2 production by the cytc3/PSI complex coupled with hydrogenase from Desulfovibrio vulgaris was enhanced 7-fold. Consequently, in the simultaneous presence of hydrogenase, Fd, and FNR, the light-driven H2 production by the hydrogenase/cytc3/PSI complex was observed (0.30 pmol Hz/mg chlorophyll/h). These results suggest that the cytc3/PSI complex may produce H2 in vivo.

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.

The Tmc complex from Desulfovibrio vulgaris hildenborough is involved in transmembrane electron transfer from periplasmic hydrogen oxidation

Three membrane-bound redox complexes have been reported in Desulfovibrio spp., whose genes are not found in the genomes of other sulfate reducers such as Desulfotalea psycrophila and Archaeoglobus fulgidus. These complexes contain a periplasmic cytochrome c subunit of the cytochrome c(3) family, and their presence in these organisms probably correlates with the presence of a pool of periplasmic cytochromes c(3), also absent in the two other sulfate reducers. In this work we report the isolation and characterization of the first of such complexes, Tmc from D. vulgaris Hildenborough, which is associated with the tetraheme type II cytochrome c(3). The isolated Tmc complex contains four subunits, including the TpIIc(3) (TmcA), an integral membrane cytochrome b (TmcC), and two cytoplasmically predicted proteins, an iron-sulfur protein (TmcB) and a tryptophan-rich protein (TmcD). Spectroscopic studies indicate the presence of eight hemes c and two hemes b in the complex pointing to an alpha(2)betagammadelta composition (TmcA(2)BCD). EPR analysis reveals the presence of a [4Fe4S](3+) center and up to three other iron-sulfur centers in the cytoplasmic subunit. Nearly full reduction of the redox centers in the Tmc complex could be obtained upon incubation with hydrogenase/TpIc(3), supporting the role of this complex in transmembrane transfer of electrons resulting from periplasmic oxidation of hydrogen.

Proton thrusters: overview of the structural and functional features of soluble tetrahaem cytochromes c 3

Tetrahaem cytochromes c (3) from sulfate-reducing bacteria have revealed exquisite complexity in their ligand binding properties and they couple the cooperative binding of two electrons with the binding of protons. In this review, the molecular mechanisms for these cooperative effects are described, and the functional consequences of these cooperativities are discussed in the context of the general mechanisms of biological energy transduction and the specific physiological metabolism of Desulfovibrio.

Roles of charged residues in pH-dependent redox properties of cytochrome c3 from Desulfovibrio vulgaris Miyazaki F

Complicated pH-properties of the tetraheme cytochrome c3 (cyt c3) from Desulfovibrio vulgaris Miyazaki F (DvMF) were examined by the pH titrations of 1H-15N HSQC spectra in the ferric and ferrous states. The redox-linked pKa shift for the propionate group at C13 of heme 1 was observed as the changes of the NH signals around it. This pKa shift is consistent with the redox-linked conformational alteration responsible for the cooperative reduction between hemes 1 and 2. On the other hand, large chemical shift changes caused by the protonation/deprotonation of Glu41 and/or Asp42, and His67 were redox-independent. Nevertheless, these charged residues affect the redox properties of the four hemes. Furthermore, one of interesting charged residues, Glu41, was studied by site-directed mutagenesis. E41K mutation increased the microscopic redox potentials of heme 1 by 46 and 34 mV, and heme 2 by 35 and 30 mV at the first and last reduction steps, respectively. Although global folding in the crystal structure of E41K cyt c3 is similar to that of wild type, local change was observed in 1H NMR spectrum. Glu41 is important to keep the stable conformation in the region between hemes 1 and 2, controlling the redox properties of DvMF cyt c3. In contrast, the kinetic parameters for electron transfer from DvMF [NiFe] hydrogenase were not influenced by E41K mutation. This suggests that the region between hemes 1 and 2 is not involved in the interaction with [NiFe] hydrogenase, and it supports the idea that heme 4 is the exclusive entrance gate to accept the electron in the initial reduction stage.

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.

Sulphate respiration from hydrogen in Desulfovibrio bacteria: a structural biology overview

Sulphate-reducing organisms are widespread in anaerobic enviroments, including the gastrointestinal tract of man and other animals. The study of these bacteria has attracted much attention over the years, due also to the fact that they can have important implications in industry (in biocorrosion and souring of oil and gas deposits), health (in inflamatory bowel diseases) and the environment (bioremediation). The characterization of the various components of the electron transport chain associated with the hydrogen metabolism in Desulfovibrio has generated a large and comprehensive list of studies. This review summarizes the more relevant aspects of the current information available on the structural data of various molecules associated with hydrogen metabolism, namely hydrogenases and cytochromes. The transmembrane redox complexes known to date are also described and discussed. Redox-Bohr and cooperativity effects, observed in a few cytochromes, and believed to be important for their functional role, are discussed. Kinetic studies performed with these redox proteins, showing clues to their functional inter-relationship, are also addressed. These provide the groundwork for the application of a variety of molecular modelling approaches to understanding electron transfer and protein interactions among redox partners, leading to the characterization of several transient periplasmic complexes. In contrast to the detailed understanding of the periplasmic hydrogen oxidation process, very little is known about the cytoplasmic side of the respiratory electron transfer chain, in terms of molecular components (with exception of the terminal reductases), their structure and the protein-protein interactions involved in sulphate reduction. Therefore, a thorough understanding of the sulphate respiratory chain in Desulfovibrio remains a challenging task.

The type I/type II cytochrome c3 complex: an electron transfer link in the hydrogen-sulfate reduction pathway

In Desulfovibrio metabolism, periplasmic hydrogen oxidation is coupled to cytoplasmic sulfate reduction via transmembrane electron transfer complexes. Type II tetraheme cytochrome c3 (TpII-c3), nine-heme cytochrome c (9HcA) and 16-heme cytochrome c (HmcA) are periplasmic proteins associated to these membrane-bound redox complexes and exhibit analogous physiological function. Type I tetraheme cytochrome c3 (TpI-c3) is thought to act as a mediator for electron transfer from hydrogenase to these multihemic cytochromes. In the present work we have investigated Desulfovibrio africanus (Da) and Desulfovibrio vulgaris Hildenborough (DvH) TpI-c3/TpII-c3 complexes. Comparative kinetic experiments of Da TpI-c3 and TpII-c3 using electrochemistry confirm that TpI-c3 is much more efficient than TpII-c3 as an electron acceptor from hydrogenase (second order rate constant k = 9 x 10(8) M(-1) s(-1), K(m) = 0.5 microM as compared to k = 1.7 x 10(7) M(-1) s(-1), K(m) = 40 microM, for TpI-c3 and TpII-c3, respectively). The Da TpI-c3/TpII-c3 complex was characterized at low ionic strength by gel filtration, analytical ultracentrifugation and cross-linking experiments. The thermodynamic parameters were determined by isothermal calorimetry titrations. The formation of the complex is mainly driven by a positive entropy change (deltaS = 137(+/-7) J mol(-1) K(-1) and deltaH = 5.1(+/-1.3) kJ mol(-1)) and the value for the association constant is found to be (2.2(+/-0.5)) x 10(6) M(-1) at pH 5.5. Our thermodynamic results reveal that the net increase in enthalpy and entropy is dominantly produced by proton release in combination with water molecule exclusion. Electrostatic forces play an important role in stabilizing the complex between the two proteins, since no complex formation is detected at high ionic strength. The crystal structure of Da TpI-c3 has been solved at 1.5 angstroms resolution and structural models of the complex have been obtained by NMR and docking experiments. Similar experiments have been carried out on the DvH TpI-c3/TpII-c3 complex. In both complexes, heme IV of TpI-c3 faces heme I of TpII-c3 involving basic residues of TpI-c3 and acidic residues of TpII-c3. A secondary interacting site has been observed in the two complexes, involving heme II of Da TpII-c3 and heme III of DvH TpI-c3 giving rise to a TpI-c3/TpII-c3 molar ratio of 2:1 and 1:2 for Da and DvH complexes, respectively. The physiological significance of these alternative sites in multiheme cytochromes c is discussed.

Redox and spectroscopic properties of the Escherichia coli nitric oxide-detoxifying system involving flavorubredoxin and its NADH-oxidizing redox partner

Under anaerobic conditions, the flavodiiron NO-reductase from Escherichia coli (flavorubredoxin, FlRd) constitutes one of the major protective enzymes against nitric oxide. The redox and spectroscopic properties of the rubredoxin (Rd), non-heme diiron, and FMN sites of flavorubredoxin were determined, which was complemented by the study of truncated versions of FlRd: one consisting only of its rubredoxin module, and another consisting of its flavodiiron structural core (lacking the Rd domain). The studies here reported were performed by a combination of potentiometry with visible and EPR spectroscopies. Moreover, we present the first direct EPR evidence for the presence of the non-heme diiron site in the flavodiiron proteins family. Also, the redox properties of the FlRd physiological partner, the NADH:flavorubredoxin oxidoreductase (FlRd-Red), were determined. It is further shown that the redox properties of this complex electron transfer system are fine-tuned upon interaction of the two enzymes.

Direct electrochemistry and Raman spectroscopy of sol-gel-encapsulated myoglobin

The direct electrochemistry of myoglobin (Mb) has been observed at a glassy carbon (GC) electrode coated with silica sol-gel-encapsulated Mb film. A well-behaved cyclic voltammogram is observed with a midpoint potential (E(1/2)) of -0.25 V vs Ag/AgCl in a pH 7.0 phosphate buffer. This potential, which is pH-dependent, is 70-90 mV more negative than the formal potential values obtained by using the spectroeletrochemical titration method at the same pH. Square wave voltametry (SWV) also shows a peak potential of -0.25 V for the reduction of Mb under the same experimental conditions. Both cathodic and anodic peak currents have a linear relationship with the scan rate. The midpoint potential decreases with pH, having a slope of -30 mV/pH. UV-vis and resonance Raman spectroscopic studies reveal that the sol-gel provides a bio-compatible environment where Mb retains a structure similar to its solution form, a 6-coordinated aquomet myoglobin. These results suggest that the silica sol-gel is a useful matrix for studying direct electrochemistry of other heme proteins.

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.

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.

The putidaredoxin reductase-putidaredoxin electron transfer complex: theoretical and experimental studies

Interaction and electron transfer between putidaredoxin reductase (Pdr) and putidaredoxin (Pdx) from Pseudomonas putida was studied by molecular modeling, mutagenesis, and stopped flow techniques. Based on the crystal structures of Pdr and Pdx, a complex between the proteins was generated using computer graphics methods. In the model, Pdx is docked above the isoalloxazine ring of FAD of Pdr with the distance between the flavin and [2Fe-2S] of 14.6 A. This mode of interaction allows Pdx to easily adjust and optimize orientation of its cofactor relative to Pdr. The key residues of Pdx located at the center, Asp(38) and Trp(106), and at the edge of the protein-protein interface, Tyr(33) and Arg(66), were mutated to test the Pdr-Pdx computer model. The Y33F, Y33A, D38N, D38A, R66A, R66E, W106F, W106A, and Delta106 mutations did not affect assembly of the [2Fe-2S] cluster and resulted in a marginal change in the redox potential of Pdx. The electron-accepting ability of Delta106 Pdx was similar to that of the wild-type protein, whereas electron transfer rates from Pdr to other mutants were diminished to various degrees with the smallest and largest effects on the kinetic parameters of the Pdr-to-Pdx electron transfer reaction caused by the Trp(106) and Tyr(33)/Arg(66) substitutions, respectively. Compared with wild-type Pdx, the binding affinity of all studied mutants to Pdr was significantly higher. Experimental results were in agreement with theoretical predictions and suggest that: (i) Pdr-Pdx complex formation is mainly driven by steric complementarity, (ii) bulky side chains of Tyr(33), Arg(66), and Trp(106) prevent tight binding of oxidized Pdx and facilitate dissociation of the reduced iron-sulfur protein from Pdr, and (iii) transfer of an electron from FAD to [2Fe-2S] can occur with various orientations between the cofactors through multiple electron transfer pathways that do not involve Trp(106) but are likely to include Asp(38) and Cys(39).