Protein--protein docking of electron transfer complexes: cytochrome c oxidase and cytochrome c

Electron transferring protein complexes form only transiently and the crystal structures of electron transfer protein--protein complexes involving cytochrome c could so far be determined only for the pairs of yeast cytochrome c peroxidase (CcP) with iso-1-cytochrome c (iso-1-cyt c) and with horse heart cytochrome c (cyt c). This article presents models from computational docking for complexes of cytochrome c oxidase (COX) from Paracoccus denitrificans with horse heart cytochrome c, and with its physiological counterpart cytochrome c552 (c552). Initial docking is performed with the FTDOCK program, which permits an exhaustive search of translational and rotational space. A filtering procedure is then applied to reduce the number of complexes to a manageable number. In a final step of structural and energetic refinement, the complexes are optimized by rigid-body energy minimization with the molecular mechanics package CHARMM. This methodology was first tested on the CcP:iso-1-cyt c complex, in which the complex with the lowest CHARMM energy has an RMSD from the crystal structure of only 1.8 A (C(alpha) carbon atoms). Notably, the crystal conformation has an even lower energy. The same procedure was then applied to COX:cyt c and COX:c552. The lowest-energy COX:cyt c complex is very similar to a docking model previously described for the complex of bovine cytochrome c oxidase with horse heart cytochrome c. For the COX:c552 complex, cytochrome c552 is found in two different orientations, depending on whether it is docked against COX from a two-subunit or from a four-subunit crystal structure, respectively. Both conformations are discussed critically in the light of the available experimental data.     

Characterization and redox properties of cytochrome c552 from Thermus thermophilus adsorbed on different self-assembled thiol monolayers, used to model the chemical environment of the redox partner

The structure of cytochrome c552 (Cyt-c552) from Thermus thermophilus shows many differences to other c-type cytochromes. The rich lysine domain close to the heme does not exist in this cytochrome, allowing us to postulate that the interaction with its redox partner must be different to the cytochrome c/cytochrome c oxidase interaction. We report a study of Cyt-c552 adsorbed on self-assembled monolayers (SAMs) of functionalized alkanethiols used to mimic the chemical properties of its redox partner (ba3-oxydase). Hydrophilic (-COOH), polar (-OH), hydrophobic (-CH3), and mixed (-OH/-CH3) SAMs grafted on roughened silver electrodes were characterized by X-ray photoelectron spectroscopy. Surface enhanced resonance Raman spectroscopy (SERRS) was employed to determine the structure and the redox properties (E degrees and number of transferred electron) of the heme of Cyt-c552 adsorbed on roughened silver electrodes coated by the different SAMs. The surface that most closely models the environment of the ba3-oxidase is a mixed SAM formed by 50% polar [Ag-(CH2)5-CH2OH] and 50% hydrophobic [Ag-(CH2)5-CH3] alkanethiols. Only the native form B1(6cLS) of Cyt-c552 is detected by SERRS when the protein is adsorbed on such a surface that promotes a protein orientation favorable for the electron transfer (number of transferred electron = 1). We shall discuss the differences and similarities of the electron-transfer mechanism of Cyt-c552 compared to cyt-c.     

基于细胞色素c的胞外电子传递过程

电活性微生物具有独特的胞外电子传递功能,在地球化学循环和环境污染修复中起着重要作用。细胞色素c在电活性微生物胞外电子传递过程中扮演了重要角色,不仅参与直接电子传递途径,还参与电子媒介介导的间接电子传递。其电子传递功能不仅对地球环境中铁、锰、碳等元素的循环具有重要作用,还应用于能源生产、废水处理、生物修复等众多领域,具有良好的应用潜力。本文以电活性微生物的2个模式菌属(希瓦氏菌属和地杆菌属)为例,综述了电活性微生物将电子由胞内转移至胞外的方式和途径,详细阐述了细胞色素c在该胞外电子传递过程中的重要作用,总结了细胞色素c介导的胞外电子传递过程所涉及的分析方法,并对微生物胞外电子传递未来的研究方向提出了展望。     

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.     

Complex structure of cytochrome c–cytochrome c oxidase reveals a novel protein–protein interaction mode

Mitochondrial cytochrome c oxidase (CcO) transfers electrons from cytochrome c (Cyt.c) to O₂ to generate H₂O, a process coupled to proton pumping. To elucidate the mechanism of electron transfer, we determined the structure of the mammalian Cyt.c–CcO complex at 2.0‐Å resolution and identified an electron transfer pathway from Cyt.c to CcO. The specific interaction between Cyt.c and CcO is stabilized by a few electrostatic interactions between side chains within a small contact surface area. Between the two proteins are three water layers with a long inter‐molecular span, one of which lies between the other two layers without significant direct interaction with either protein. Cyt.c undergoes large structural fluctuations, using the interacting regions with CcO as a fulcrum. These features of the protein–protein interaction at the docking interface represent the first known example of a new class of protein–protein interaction, which we term “soft and specific”. This interaction is likely to contribute to the rapid association/dissociation of the Cyt.c–CcO complex, which facilitates the sequential supply of four electrons for the O₂ reduction reaction.     

Effect of surface coverage of gold(111) electrode with cysteine on the chiral discrimination of DOPA

The enantioselectivity imparted to a gold electrode by modifying its surface with a self-assembled monolayer (SAM) of cysteine (Cys) was investigated for the electrochemical redox reaction of 3,4-dihydroxyphenylalanine (DOPA). A cyclic voltammetric study of the redox reaction revealed that the enantioselectivity was determined by the surface coverage of the gold electrode with Cys molecules. The electrode modified with approximately 1.8 x 10(14) Cys molecules cm(-2) exhibited enantioselectivity in the voltammogram for the oxidation and reduction of DOPA, while the voltammograms obtained by the electrodes with either more or less surface coverages did not exhibit significant enantioselectivity. It is suggested that the accessibility of DOPA to that area of the gold surface which is not blocked by Cys molecules at an optimum surface coverage, is required for the enantioselective redox reaction of DOPA to proceed.     

Replacements in a conserved leucine cluster in the hydrophobic heme pocket of cytochrome c.

A cluster of highly conserved leucine side chains from residues 9, 68, 85, 94, and 98 is located in the hydrophobic heme pocket of cytochrome c. The contributions of two of these, Leu 85 and Leu 94, have been studied using a protein structure-function-mutagenesis approach to probe their roles in the maintenance of overall structural integrity and electron transfer activity. Structural studies of the L85C, L85F, L85M, and L94S mutant proteins show that, in each case, the overall fold of cytochrome c is retained, but that localized conformational shifts are required to accommodate the introduced side chains. In particular, the side chains of Cys 85 and Phe 85 form energetically favorable interactions with Phe 82, whereas Met 85 takes on a more remote conformation to prevent an unfavorable interaction with the phenyl ring of Phe 82. In the case of the L94S mutant protein, the new polar group introduced is found to form hydrogen bonds to nearby carbonyl groups. In all proteins with substitutions at Leu 85, the hydrophobic nature of the heme pocket is preserved and no significant decrease in heme reduction potential is observed. Despite earlier predictions that Leu 85 is an important determinant in cytochrome c electron transfer partner complexation, our studies show this is unlikely to be the case because the considerable surface contour perturbations made by substitutions at this residue do not correspondingly translate into significant changes in electron transfer rates. For the L94S mutant protein, the substitution of a polar hydroxyl group directly into the hydrophobic heme pocket has a larger effect on heme reduction potential, but this is mitigated by two factors. First, the side chain of Ser 94 is rotated away from the heme group and, second, the side chain of Leu 98 shifts into a portion of the new space available, partially shielding the heme group. The Leu 94 Ser substitution does not perturb the highly conserved interface formed by the nearly perpendicular packing of the N- and C-terminal helices of cytochrome c, ruling this out as the cause of this mutant protein becoming thermally labile and having a lower functional activity. Our results show these effects are most likely attributable to disruption of the heme pocket region. Much of the ability of cytochrome c to absorb the introduction of mutations at Leu 85 and Leu 94 appears to be a consequence of the conformational flexibility afforded by the leucine cluster in this region as well as the presence of a nearby internal cavity.(ABSTRACT TRUNCATED AT 400 WORDS)     

On the role of cytochrome c8 in photosynthetic electron transfer of the purple non-sulfur bacterium Rhodoferax fermentans

We report on the isolation, purification and functional characterization of a soluble c-type cytochrome from light-grown cells of the purple phototroph Rhodoferax fermentans. This cytochrome is basic (pI = 8), has a molecular mass of 12 kDa, and is characterized by a midpoint reduction potential of +285 mV. Partial analysis of the N-terminus amino-acid sequence shows a high similarity with cytochromes of c₈ type (formerly called Pseudomonas cytochrome c-551 type). Time-resolved spectrophotometric studies show that this cytochrome c₈ reduces the tetraheme subunit of the photosynthetic reaction center, in a fast (sub-ms) and a slow (ms) phase. Competition experiments in the presence of both cytochrome c₈ and high potential iron-sulfur protein (HiPIP), isolated from the same microorganism, show that cytochrome c₈ oxidation is decreased upon addition of HiPIP. These observations suggest that cytochrome c₈ and HiPIP might play alternative roles in the photosynthetic electron flow of Rhodoferax fermentans.     

X-ray structure determination of the cytochrome c2: reaction center electron transfer complex from Rhodobacter sphaeroides

In the photosynthetic bacterium Rhodobacter sphaeroides, a water soluble cytochrome c2 (cyt c2) is the electron donor to the reaction center (RC), the membrane-bound pigment-protein complex that is the site of the primary light-induced electron transfer. To determine the interactions important for docking and electron transfer within the transiently bound complex of the two proteins, RC and cyt c2 were co-crystallized in two monoclinic crystal forms. Cyt c2 reduces the photo-oxidized RC donor (D+), a bacteriochlorophyll dimer, in the co-crystals in approximately 0.9 micros, which is the same time as measured in solution. This provides strong evidence that the structure of the complex in the region of electron transfer is the same in the crystal and in solution. X-ray diffraction data were collected from co-crystals to a maximum resolution of 2.40 A and refined to an R-factor of 22% (R(free)=26%). The structure shows the cyt c2 to be positioned at the center of the periplasmic surface of the RC, with the heme edge located above the bacteriochlorophyll dimer. The distance between the closest atoms of the two cofactors is 8.4 A. The side-chain of Tyr L162 makes van der Waals contacts with both cofactors along the shortest intermolecular electron transfer pathway. The binding interface can be divided into two domains: (i) A short-range interaction domain that includes Tyr L162, and groups exhibiting non-polar interactions, hydrogen bonding, and a cation-pi interaction. This domain contributes to the strength and specificity of cyt c2 binding. (ii) A long-range, electrostatic interaction domain that contains solvated complementary charges on the RC and cyt c2. This domain, in addition to contributing to the binding, may help steer the unbound proteins toward the right conformation.     

Interpretation of ionic strength dependencies of the electron transfer in the Cytochrome c-Cytochrome b 5 system for the wild-type and lysine-mutated yeast Cytochrome c

Ionic strength dependencies of electron transfer between Cytochrome b ₅ and variants of yeast Cytochrome c were analyzed by curve fitting to the simple model of the electrostatic interaction between the two proteins assuming the process to be non-diffusion-controlled. Mutagenesis of Lys79, but not Lys72, leads to an increase of effective radius of the interacting charged species, suggesting that the mutation effects of the two residues on the electrostatic field distribution near the contact site are different, even within the crude electrostatic model used. Extrapolation of the ionic strength dependencies to infinite ionic strength resulted in similar values, around (2–3)×10^-6 for all Cyt c variants considered thus showing the lysine residue mutations to primarily affect protein association rather than the electron transfer directly. Based on the ionic strength dependencies of binding constants of the two proteins into an electrostatically stabilized complex, the monomolecular electron transfer rate constant was estimated to be 1.1×10⁴–1.6×10⁵ s^-1. The electrostatic part of the binding energy of the complex at I=0.19 was estimated to be-2.4 kcal/mol, strongly at variance with the values-13.0 and-6.4 kcal/mol reported for the two types of complexes identified using Brownian dynamics techniques. Possible reasons for this discrepancy are discussed.     

Kinetic evidence for the re-definition of electron transfer pathways from cytochrome c to O2 within cytochrome oxidase

The reaction with O2 of equimolar mixtures of cytochrome c and cytochrome c oxidase in high and low ionic strength buffers has been examined by flow-flash spectrophotometry at room temperature. In low ionic strength media where cytochrome c and the oxidase are bound in an electrostatic, 1:1 complex some of the cytochrome c is oxidised at a faster rate than a metal centre of the oxidase. In contrast, when cytochrome c and cytochrome c oxidase are predominantly dissociated at high ionic strength cytochrome c oxidation occurs only slowly (t1/2 = 5 s) following the complete oxidation of the oxidase. These results demonstrate that maximal rates of electron transfer from cytochrome c to O2 occur when both substrates are present on the enzyme. The heterogeneous oxidation of cytochrome c observed in the complex implies more than one route for electron transfer within the enzyme. Possibilities for new electron transfer pathways from cytochrome c to O2 are proposed.     

High resolution X-ray crystallographic structure of bovine heart cytochrome c and its application to the design of an electron transfer biosensor

Cytochrome c is one of the most studied proteins probably due to its electron-transfer properties in aerobic and anaerobic respiration. Particularly, cytochrome c from bovine heart is a small protein, M(r) 12,230 Da, globular (hydrodynamic diameter of 3.4 nm), soluble in different buffer solutions, and commercially available. Despite being a quite well-studied protein and relatively easy to manipulate from the biochemical and electrochemical viewpoint, its 3D structure has never been published. In this work, the purification, crystallization and 3D structure of one of the cytochrome c isoforms is presented to 1.5 A resolution. It is also shown how the presence of isoforms made both the purification and crystallization steps difficult. Finally, a new approach for protein electrocrystallization and design of biosensors is presented.     

Structural and functional studies on the tetraheme cytochrome subunit and its electron donor proteins: the possible docking mechanisms during the electron transfer reaction

The photosynthetic reaction centers (RCs) classified as the group II possess a peripheral cytochrome (Cyt) subunit, which serves as the electron mediator to the special-pair. In the cycle of the photosynthetic electron transfer reactions, the Cyt subunit accepts electrons from soluble electron carrier proteins, and re-reduces the photo-oxidized special-pair of the bacteriochlorophyll. Physiologically, high-potential cytochromes such as the cytochrome c₂ and the high-potential iron–sulfur protein (HiPIP) function as the electron donors to the Cyt subunit. Most of the Cyt subunits possess four heme c groups, and it was unclear which heme group first accepts the electron from the electron donor. The most distal heme to the special-pair, the heme-1, has a lower redox potential than the electron donors, which makes it difficult to understand the electron transfer mechanism mediated by the Cyt subunit. Extensive mutagenesis combined with kinetic studies has made a great contribution to our understanding of the molecular interaction mechanisms, and has demonstrated the importance of the region close to the heme-1 in the electron transfer. Moreover, crystallographic studies have elucidated two high-resolution three-dimensional structures for the RCs containing the Cyt subunit, the Blastochloris viridis and Thermochromatium tepidum RCs, as well as the structures of their electron donors. An examination of the structural data also suggested that the binding sites for both the cytochrome c₂ and the HiPIP are located adjacent to the solvent-accessible edge of the heme-1. In addition, it is also indicated by the structural and biochemical data that the cytochrome c₂ and the HiPIP dock with the Cyt subunit by different mechanisms although the two electron donors utilize the same region for the interactions; cytochrome c₂ is recognized through electrostatic interactions while hydrophobic interactions are important in the HiPIP docking.     

Crystal structure of a ferredoxin reductase for the CYP199A2 system from Rhodopseudomonas palustris

Cytochrome P450‐199A2 from Rhodopseudomonas palustris oxidizes para‐substituted benzoic acids and may play a role in lignin and aromatic acid degradation pathways in the bacterium. CYP199A2 has an associated [2Fe‐2S] ferredoxin, palustrisredoxin (Pux) but not a ferredoxin reductase. A genome search identified the palustrisredoxin reductase (PuR) gene. PuR was produced in Escherichia coli and shown to be a flavin‐dependent protein that supports efficient electron transfer from NADH to Pux, thus reconstituting CYP199A2 monooxygenase activity (k_cat = 37.9 s^–1 with 4‐methoxybenzoic acid). The reduction of Pux by PuR shows K_m = 4.2 μM and k_cat = 262 s^–1 in 50 mM Tris, pH 7.4. K_m is increased to 154 μM in the presence of 200 mM KCl, indicating the importance of ionic interactions in PuR/Pux binding. The crystal structure of PuR has been determined at 2.2 Å resolution and found to be closely related to that of other oxygenase‐coupled NADH‐dependent ferredoxin reductases. Residues on the surface that had been proposed to be involved in ferredoxin reductase‐ferredoxin binding are conserved in PuR. However, Lys328 in PuR lies over the FAD isoalloxazine ring and, together with His11 and Gln41, render the electrostatic potential of the surface more positive and may account for the greater involvement of electrostatic interactions in ferredoxin binding by PuR. Consistent with these observations the K328G mutation weakened Pux binding and virtually eliminated the dependence of PuR/Pux binding on salt concentration, thus confirming that the FAD si side surface in the vicinity of Lys328 is the ferredoxin binding site. Proteins 2009. © 2009 Wiley‐Liss, Inc.     

Intermonomer electron transfer in the bc1 complex dimer is controlled by the energized state and by impaired electron transfer between low and high potential hemes

The cytochrome bc(1) complex (commonly called Complex III) is the central enzyme of respiratory and photosynthetic electron transfer chains. X-ray structures have revealed the bc(1) complex to be a dimer, and show that the distance between low potential (b(L)) and high potential (b(H)) hemes, is similar to the distance between low potential hemes in different monomers. This suggests that electron transfer between monomers should occur at the level of the b(L) hemes. Here, we show that although the rate constant for b(L)-->b(L) electron transfer is substantial, it is slow compared to the forward rate from b(L) to b(H), and the intermonomer transfer only occurs after equilibration within the first monomer. The effective rate of intermonomer transfer is about 2-orders of magnitude slower than the direct intermonomer electron transfer.     

Elucidating the mechanisms of coupled electron transfer and catalytic reactions by protein film voltammetry

Protein film voltammetry, the direct electrochemistry of redox enzymes and proteins, provides precise and comprehensive information on complicated reaction mechanisms. By controlling the driving force for a reaction (using the applied potential) and monitoring the reaction in real time (using the current), it allows thermodynamic and kinetic information to be determined simultaneously. Two challenges are inherent to protein film voltammetry: (i) to adsorb the protein or enzyme in a native and active configuration on the electrode surface, and (ii) to understand and interpret voltammetric results on both a qualitative and quantitative level, allowing mechanistic models to be proposed and rigorous experiments to test these models to be devised. This review focuses on the second of these two challenges. It describes how to use protein film voltammetry to derive mechanistic and biochemically relevant information about redox proteins and enzymes, and how to evaluate and interpret voltammetric results. Selected key studies are described in detail, to illustrate their underlying principles, strategies and physical interpretations.     

Independent tyrosyl contributions to the CD of Ff gene 5 protein and the distinctive effects of Y41H and Y41F mutants on protein-protein cooperative interactions

The gene 5 protein (g5p) of the Ff virus contains five Tyr, individual mutants of which have now all been characterized by CD spectroscopy. The protein has a dominant tyrosyl 229-nm L(a) CD band that is shown to be approximately the sum of the five individual Tyr contributions. Tyr41 is particularly important in contributing to the high cooperativity with which the g5p binds to ssDNA, and Y41F and Y41H mutants are known to differ in dimer-dimer packing interactions in crystal structures. We compared the solution structures and binding properties of the Y41F and Y41H mutants using CD spectroscopy. Secondary structures of the mutants were similar by CD analyses and close to those derived from the crystal structures. However, there were significant differences in the binding properties of the two mutant proteins. The Y41H protein had an especially low binding affinity and perturbed the spectrum of poly[d(A)] in 2 mM Na(+) much less than did Y41F and the wild-type gene 5 proteins. Moreover, a change in the Tyr 229 nm band, assigned to the perturbation of Tyr34 at the dimer-dimer interface, was absent in titrations with the Y41H mutant under low salt conditions. In contrast, titrations with the Y41H mutant in 50 mM Na(+) exhibited typical CD changes of both the nucleic acid and the Tyr 229-nm band. Thus, protein-protein and g5p-ssDNA interactions appeared to be mutually influenced by ionic strength, indicative of correlated changes in the ssDNA binding and cooperativity loops of the protein or of indirect structural constraints.