Tetraheme cytochromes

Differences in Bioenergetic Metabolism of Obligately Alkaliphilic Bacillaceae Under High pH Depend on the Aeration Conditions

Alkaliphilic Bacillaceae appear to produce ATP based on the H⁺-based chemiosmotic theory. However, the bulk-based chemiosmotic theory cannot explain the ATP production in alkaliphilic bacteria because the H⁺ concentration required for driving ATP synthesis through the ATPase does not occur under the alkaline conditions. Alkaliphilic bacteria produce ATP in an H⁺-diluted environment by retaining scarce H⁺ extruded by the respiratory chain on the outer surface of the membrane and increasing the potential of the H⁺ for ATP production on the outer surface of the membrane using specific mechanisms of ATP production. Under high-aeration conditions, the high ΔΨ (ca. -170 mV) of the obligate alkaliphilic Evansella clarkii retains H⁺ at the outer surface of the membrane and increases the intensity of the protonmotive force (Δp) per H⁺ across the membrane. One of the reasons for the production of high ΔΨ is the Donnan potential, which arises owing to the induction of impermeable negative charges in the cytoplasm. The intensity of the potential is further enhanced in the alkaliphiles compared with neutralophiles because of the higher intracellular pH (ca. pH 8.1). However, the high ΔΨ observed under high-aeration conditions decreased (∼ -140 mV) under low-aeration conditions. E. clarkii produced 2.5–6.3-fold higher membrane bound cytochrome c in the content of the cell extract under low-aeration conditions than under high-aeration conditions. The predominant membrane-bound cytochrome c in the outer surface of the membrane possesses an extra Asn-rich segment between the membrane anchor and the main body of protein. This structure may influence the formation of an H⁺-bond network that accumulates H⁺ on the outer surface of the membrane. Following accumulation of the H⁺-bond network producing cytochrome c, E. clarkii constructs an H⁺ capacitor to overcome the energy limitation of low aeration at high pH conditions. E. clarkii produces more ATP than other neutralophilic bacteria by enhancing the efficacy per H⁺ in ATP synthesis. In low H⁺ environments, E. clarkii utilizes H⁺ efficiently by taking advantage of its high ΔΨ under high-aeration conditions, whereas under low-aeration conditions E. clarkii uses cytochrome c bound on its outer surface of the membrane as an H⁺ capacitor.

Exploring Coupled Redox and pH Processes with a Force-Field-Based Approach: Applications to Five Different Systems

Coupled redox and pH-driven processes are at the core of many important biological mechanisms. As the distribution of protonation and redox states in a system is associated with the pH and redox potential of the solution, having efficient computational tools that can simulate under these conditions becomes very important. Such tools have the potential to provide information that complement and drive experiments. In previous publications we have presented the implementation of the constant pH and redox potential molecular dynamics (C(pH,E)MD) method in AMBER and we have shown how multidimensional replica exchange can be used to significantly enhance the convergence efficiency of our simulations. In the current work, after an improvement in our C(pH,E)MD approach that allows a given residue to be simultaneously pH- and redox-active, we have employed our methodologies to study five different systems of interest in the literature. We present results for capped tyrosine dipeptide, two maquette systems containing one pH- and redox-active tyrosine (α₃Y and peptide A), and two proteins that contain multiple heme groups (diheme cytochrome c from Rhodobacter sphaeroides and Desulfovibrio vulgaris Hildenborough cytochrome c₃). We show that our results can provide new insights into previous theoretical and experimental findings by using a fully force-field-based and GPU-accelerated approach, which allows the simulations to be executed with high computational performance.

Formation of Proton Motive Force Under Low-Aeration Alkaline Conditions in Alkaliphilic Bacteria

In Mitchell’s chemiosmotic theory, a proton (H⁺) motive force across the membrane (Δp), generated by the respiratory chain, drives F₁F_o-ATPase for ATP production in various organisms. The bulk-base chemiosmotic theory cannot account for ATP production in alkaliphilic bacteria. However, alkaliphiles thrive in environments with a H⁺ concentrations that are one-thousandth (ca. pH 10) the concentration required by neutralophiles. This situation is similar to the production of electricity by hydroelectric turbines under conditions of very limited water. Alkaliphiles manage their metabolism via various strategies involving the cell wall structure, solute transport systems and molecular mechanisms on the outer surface membrane. Our experimental results indicate that efficient ATP production in alkaliphilic Bacillus spp. is attributable to a high membrane electrical potential (ΔΨ) generated for an attractive force for H⁺ on the outer surface membrane. In addition, the enhanced F₁F_o-ATPase driving force per H⁺ is derived from the high ΔΨ. However, it is difficult to explain the reasons for high ΔΨ formation based on the respiratory rate. The Donnan effect (which is observed when charged particles that are unable to pass through a semipermeable membrane create an uneven electrical charge) likely contributes to the formation of the high ΔΨ because the intracellular negative ion capacities of alkaliphiles are much higher than those of neutralophiles. There are several variations in the adaptation to alkaline environments by bacteria. However, it could be difficult to utilize high ΔΨ in the low aeration condition due to the low activity of respiration. To explain the efficient ATP production occurring in H⁺-less and air-limited environments in alkaliphilic bacteria, we propose a cytochrome c-associated “H⁺ capacitor mechanism” as an alkaline adaptation strategy. As an outer surface protein, cytochrome c-550 from Bacillusclarkii possesses an extra Asn-rich segment between the region anchored to the membrane and the main body of the cytochrome c. This structure may contribute to the formation of the proton-binding network to transfer H⁺ at the outer surface membrane in obligate alkaliphiles. The H⁺ capacitor mechanism is further enhanced under low-aeration conditions in both alkaliphilic Bacillus spp. and the Gram-negative alkaliphile Pseudomonas alcaliphila.

Design and fine-tuning redox potentials of metalloproteins involved in electron transfer in bioenergetics

Redox potentials are a major contributor in controlling the electron transfer (ET) rates and thus regulating the ET processes in the bioenergetics. To maximize the efficiency of the ET process, one needs to master the art of tuning the redox potential, especially in metalloproteins, as they represent major classes of ET proteins. In this review, we first describe the importance of tuning the redox potential of ET centers and its role in regulating the ET in bioenergetic processes including photosynthesis and respiration. The main focus of this review is to summarize recent work in designing the ET centers, namely cupredoxins, cytochromes, and iron-sulfur proteins, and examples in design of protein networks involved these ET centers. We then discuss the factors that affect redox potentials of these ET centers including metal ion, the ligands to metal center and interactions beyond the primary ligand, especially non-covalent secondary coordination sphere interactions. We provide examples of strategies to fine-tune the redox potential using both natural and unnatural amino acids and native and nonnative cofactors. Several case studies are used to illustrate recent successes in this area. Outlooks for future endeavors are also provided. This article is part of a Special Issue entitled Biodesign for Bioenergetics--the design and engineering of electronic transfer cofactors, proteins and protein networks, edited by Ronald L. Koder and J.L. Ross Anderson.

Unveiling the details of electron transfer in multicenter redox proteins

Metalloproteins modulate the intrinsic properties of transition metals to achieve controlled catalysis, electron transfer, or structural stabilization. Those performing electron transport, redox proteins, are a diverse class of proteins with central roles in numerous metabolic and signaling pathways, including respiration and photosynthesis. Many redox proteins have applications in industry, especially biotechnology, making them the focus of intense research. Redox proteins may contain one or multiple redox centers of the same or a different type. The complexity of proteins with multiple redox centers makes it difficult to establish a detailed molecular mechanism for their activity. Thermodynamic and kinetic information can be interpreted using the molecular structure to elucidate the protein's functional mechanism. This Account reviews experimental strategies developed in recent years to determine the detailed thermodynamic properties of multicenter redox proteins and their kinetic properties during interactions with redox partners. These strategies allow the discrimination of thermodynamic and kinetic properties of each individual redox center. The thermodynamic characterization of the redox transitions results from the combined analysis of data from NMR and UV-visible spectroscopy. Meanwhile, the kinetic characterization of intermolecular electron transfer comes from stopped-flow spectrophotometry. We illustrate an application of these strategies to a particular redox protein, the small tetraheme cytochrome from the periplasmic space of Shewanella oneidensis MR-1. This protein is a convenient prototype for developing methods for the detailed analysis of multicenter electron transfer proteins because hemes have strong UV-visible absorption bands and because heme resonances have exquisite discrimination in NMR spectra. Nonetheless, the methods are fully generalizable. Ultimately, this Account highlights the relevance of detailed characterization of the thermodynamic and kinetic properties of redox proteins. These properties are responsible for the directionality and specificity of the electron transfer process in bioenergetic pathways; a more thorough characterization of these properties should allow better-designed proteins for industrial applications.

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.

Hemoproteins in dissimilatory sulfate- and sulfur-reducing prokaryotes

Dissimilatory sulfate and sulfur reduction evolved billions of years ago and while the bacteria and archaea that use this unique metabolism employ a variety of electron donors, H(2) is most commonly used as the energy source. These prokaryotes use multiheme c-type proteins to shuttle electrons from electron donors, and electron transport complexes presumed to contain b-type hemoproteins contribute to proton charging of the membrane. Numerous sulfate and sulfur reducers use an alternate pathway for heme synthesis and, frequently, uniquely specific axial ligands are used to secure c-type heme to the protein. This review presents some of the types and functional activities of hemoproteins involved in these two dissimilatory reduction pathways.

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.

Surface-induced intramolecular electron transfer in multi-centre redox metalloproteins: the di-haem protein cytochrome c(4) in homogeneous solution and at electrochemical surfaces

Intramolecular electron transfer (ET) between transition metal centres is a core feature of biological ET and redox enzyme function. The number of microscopic redox potentials and ET rate constants is, however, mostly prohibitive for experimental mapping, but two-centre proteins offer simple enough communication networks for complete mapping to be within reach. At the same time, multi-centre redox proteins operate in a membrane environment where conformational dynamics and ET patterns are quite different from the conditions in a homogeneous solution. The bacterial respiratory di-haem protein Pseudomonas stutzeri cytochrome c(4) offers a prototype target for environmental gating of intra-haem ET. ET between P. stutzeri cyt c(4) and small molecular reaction partners in solution appears completely dominated by intermolecular ET of each haem group/protein domain, with no competing intra-haem ET, for which accompanying propionate-mediated proton transfer is a further barrier. The protein can, however, be immobilized on single-crystal, modified Au(111) electrode surfaces with either the low-potential N terminal or the high-potential C terminal domain facing the surface, clearly with fast intramolecular ET as a key feature in the electrochemical two-ET process. This dual behaviour suggests a pattern for multi-centre redox metalloprotein function. In a homogeneous solution, which is not the natural environment of cyt c(4), the two haem group domains operate largely independently with conformations prohibitive for intramolecular ET. Binding to a membrane or electrochemical surface, however, triggers conformational opening of intramolecular ET channels. The haem group orientation in P. stutzeri cyt c(4) is finally noted to offer a case for orientation dependent electronic rectification between a substrate and a tip in electrochemical in situ scanning tunnelling microscopy or nanoscale electrode configurations.

Thermodynamic and kinetic characterisation of individual haems in multicentre cytochromes c3

The characterisation of individual centres in multihaem proteins is difficult due to the similarities in the redox and spectroscopic properties of the centres. NMR has been used successfully to distinguish redox centres and allow the determination of the microscopic thermodynamic parameters in several multihaem cytochromes c(3) isolated from different sulphate-reducing bacteria. In this article we show that it is also possible to discriminate the kinetic properties of individual centres in multihaem proteins, if the complete microscopic thermodynamic characterisation is available and the system displays fast intramolecular equilibration in the time scale of the kinetic experiment. The deconvolution of the kinetic traces using a model of thermodynamic control provides a reference rate constant for each haem that does not depend on driving force and can be related to structural factors. The thermodynamic characterisation of three tetrahaem cytochromes and their kinetics of reduction by sodium dithionite are reported in this paper. Thermodynamic and kinetic data were fitted simultaneously to a model to obtain microscopic reduction potentials, haem-haem and haem-proton interacting potentials, and reference rate constants for the haems. The kinetic information obtained for these cytochromes and recently published data for other multihaem cytochromes is discussed with respect to the structural factors that determine the reference rates. The accessibility for the reducing agent seems to play an important role in controlling the kinetic rates, although is clearly not the only factor.

Desulfovibrio desulfuricans G20 Tetraheme Cytochrome Structure at 1.5Å and Cytochrome Interaction with Metal Complexes

The structure of the type I tetraheme cytochrome c(3) from Desulfovibrio desulfuricans G20 was determined to 1.5 Angstrom by X-ray crystallography. In addition to the oxidized form, the structure of the molybdate-bound form of the protein was determined from oxidized crystals soaked in sodium molybdate. Only small structural shifts were obtained with metal binding, consistent with the remarkable structural stability of this protein. In vitro experiments with pure cytochrome showed that molybdate could oxidize the reduced cytochrome, although not as rapidly as U(VI) present as uranyl acetate. Alterations in the overall conformation and thermostability of the metal-oxidized protein were investigated by circular dichroism studies. Again, only small changes in protein structure were documented. The location of the molybdate ion near heme IV in the crystal structure suggested heme IV as the site of electron exit from the reduced cytochrome and implicated Lys14 and Lys56 in binding. Analysis of structurally conserved water molecules in type I cytochrome c(3) crystal structures identified interactions predicted to be important for protein stability and possibly for intramolecular electron transfer among heme molecules.

Heterologous expression of dodecaheme “nanowire” cytochromes c from Geobacter sulfurreducens

Multiheme cytochromes c are difficult to produce in heterologous systems. The genome of delta-proteobacterium Geobacter sulfurreducens contains more than a hundred genes coding for c-type cytochromes. Among those are two dodecaheme cytochromes c representing a new class of multiheme cytochromes, whose putative structure is a one-dimensional array of small highly homologous domains that contain three hemes and are covalently bound by short linkers. They are likely to form "nanowires" that are part of the electron transfer chain. We cloned the genes coding for the two cytochromes into a vector we developed for ligation-independent cloning of proteins targeted to the Escherichia coli periplasmic space. We expressed the proteins in E. coli co-transformed with a plasmid harboring the cytochrome c maturation genes. Expression levels were optimized by varying IPTG concentrations used for induction. Although both proteins appeared insoluble or strongly associated with cell membranes, they were solubilized using 0.5 M sodium chloride which was more selective than conventional solubilizing agents, such as HEGA-10 or beta-octylglucoside. The solubilized proteins were dialyzed and purified by cation exchange chromatography followed by gel filtration. Mass-spectrometry analysis confirmed that both purified proteins contained the complete set of covalently attached hemes, 12 per molecule. Their visible spectra were typical of c-type cytochromes. Both proteins were successfully crystallized.

Solution structures of tetrahaem ferricytochrome c3 from Desulfovibrio vulgaris (Hildenborough) and its K45Q mutant: The molecular basis of cooperativity

The NMR structure of the oxidised wild-type cytochrome c3 from Desulfovibrio vulgaris Hildenborough was determined in solution. Using a newly developed methodology, NMR data from the K45Q mutant was then grafted onto data from the wild-type protein to determine the structure in the region of the mutation. The structural origins of the redox-Bohr effect and haem-haem cooperativities are discussed with respect to the redox-related conformational changes observed in solution.

Cytochrome c and bioenergetic hypothetical model for alkaliphilic Bacillus spp

Although a bioenergetic parameter is unfavorable for production of ATP (DeltapH<0), the growth rate and yield of alkaliphilic Bacillus strains are higher than those of neutralophilic Bacillus subtilis. This finding suggests that alkaliphiles possess a unique energy-producing machinery taking advantage of the alkaline environment. Expected bioenergetic parameters for the production of ATP (DeltapH and DeltaPsi) do not reflect the actual parameters for energy production. Certain strains of alkaliphilic Bacillus spp. possess large amounts of cytochrome c when grown at a high pH. The growth rate and yield are higher at pH 10 than at pH 7 in facultative alkaliphiles. These findings suggest that a large amount of cytochrome c at high pHs (e.g., pH 10) may be advantageous for sustaining growth. To date, isolated cytochromes c of alkaliphiles have a very low midpoint redox potential (less than +100 mV) compared with those of neutralophiles (approximately +220 mV). On the other hand, the redox potential of the electron acceptor from cytochrome c, that is, cytochrome c oxidase, seems to be normal (redox potential of cytochrome a=+250 mV). This large difference in midpoint redox potential between cytochrome c and cytochrome a concomitant with the configuration (e.g., a larger negative ion capacity at the inner surface membrane than at the outer surface for the attraction of H+ to the intracellular membrane and a large amount of cyrochrome c) supporting H+-coupled electron transfer of cytochrome c may have an important meaning in the adaptation of alkaliphiles at high pHs. This respiratory system includes a more rapid and efficient H+ and e- flow across the membrane in alkaliphiles than in neutralophiles.

Resonance Raman fingerprinting of multiheme cytochromes from the cytochrome c 3 family

Resonance Raman (RR) spectroscopy was used to investigate conformational characteristics of the hemes of several ferricytochromes of the cytochrome c3 family, electron transfer proteins isolated from the periplasm and membranes of sulfate-reducing bacteria. Our analysis concentrated on the low-frequency region of the RR spectra, a fingerprint region that includes vibrations for heme-protein C-S bonds [nu(C(a)S)]. It has been proposed that these bonds are directly involved in the electron transfer process. The three groups of tetraheme cytochrome c3 analyzed, namely Type I cytochrome c (3) (TpIc (3)s), Type II cytochrome c (3) (TpIIc (3)s) and Desulfomicrobium cytochromes c3, display different frequency separations for the two nu(C(a)S) lines that are similar among members of each group. These spectral differences correlate with differences in protein structure observed among the three groups of cytochromes c3. Two larger cytochromes of the cytochrome c3 family display RR spectral characteristics for the nu(C(a)S) lines that are closer to TpIIc3 than to TpIc3. Two other multiheme cytochromes from Desulfovibrio that do not belong to the cytochrome c3 family display nu(C(a)S) lines with reverse relative areas in comparison with the latter family. This RR study shows that the small differences in protein structure observed among these cytochrome c3 correlate to differences on the heme-protein bonds, which are likely to have an impact upon the protein function, making RR spectroscopy a sensitive and useful tool for characterizing these cytochromes.

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.

Reorganization and conformational changes in the reduction of tetraheme cytochromes

Molecular dynamics simulation (MD) constitutes an alternative to time-consuming experiments for studying conformational changes. We apply MD on a redox system where experimental information exists for the fully oxidized and fully reduced states: tetraheme cytochrome c3. Instead of doing one simulation for each state, we apply 10 4-ns replicas for both states, which provides robust statistics to characterize the redox changes. Besides these long simulations, we perform 120 short ones (50 ps), where an equilibrated oxidized state is perturbed to a reduced state. This allows the application of a nonequilibrium method, the subtraction technique, which makes it possible to characterize the different timescales of conformational changes. Reduction induces conformational changes in the N-terminus and on the loops spanning residues 36-42 and 88-93, which correlate very well with experiments, demonstrating the applicability of this methodology. We also analyze the effect of reduction on hydrogen bonds, solvent accessible surface and bound water, the changes being found to involve the hemes and propionate groups. Redox-induced protonation is also investigated, by protonating the propionates D from hemes I and IV. Although this change in the former does not have major conformational consequences, it induces in the latter conformational changes beyond the ones obtained with reduction.

Interaction and electron transfer between the high molecular weight cytochrome and cytochrome c3 from Desulfovibrio vulgaris Hildenborough: Kinetic, microcalorimetric, EPR and electrochemical studies

The complex formation between the tetraheme cytochrome c3 and hexadecaheme high molecular weight cytochrome c (Hmc), the structure of which has recently been resolved, has been characterized by cross-linking experiments, EPR, electrochemistry and kinetic analysis, and some key parameters of the interaction were determined. The analysis of electron transfer between [Fe] hydrogenase, cytochrome c3 and Hmc demonstrates a redox-shuttling role of cytochrome c3 in the pathway from hydrogenase to Hmc, and shows an effect of redox state on the interaction between the two cytochromes. The role of polyheme cytochromes in electron transfer from periplasmic hydrogenase to membrane redox proteins is assessed. A model with cytochrome c3 as an intermediate between hydrogenase and various polyheme cytochromes is proposed and its physiological consequences are discussed.

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.

Structure of a novel c7-type three-heme cytochrome domain from a multidomain cytochrome c polymer

The structure of a novel c(7)-type cytochrome domain that has two bishistidine coordinated hemes and one heme with histidine, methionine coordination (where the sixth ligand is a methionine residue) was determined at 1.7 A resolution. This domain is a representative of domains that form three polymers encoded by the Geobacter sulfurreducens genome. Two of these polymers consist of four and one protein of nine c(7)-type domains with a total of 12 and 27 hemes, respectively. Four individual domains (termed A, B, C, and D) from one such multiheme cytochrome c (ORF03300) were cloned and expressed in Escherichia coli. The domain C produced diffraction quality crystals from 2.4 M sodium malonate (pH 7). The structure was solved by MAD method and refined to an R-factor of 19.5% and R-free of 21.8%. Unlike the two c(7) molecules with known structures, one from G. sulfurreducens (PpcA) and one from Desulfuromonas acetoxidans where all three hemes are bishistidine coordinated, this domain contains a heme which is coordinated by a methionine and a histidine residue. As a result, the corresponding heme could have a higher potential than the other two hemes. The apparent midpoint reduction potential, E(app), of domain C is -105 mV, 50 mV higher than that of PpcA.

Role of the aromatic ring of Tyr43 in tetraheme cytochrome c(3) from Desulfovibrio vulgaris Miyazaki F

Tyrosine 43 is positioned parallel to the fifth heme axial ligand, His34, of heme 1 in the tetraheme cytochrome c(3). The replacement of tyrosine with leucine increased the redox potential of heme 1 by 44 and 35 mV at the first and last reduction steps, respectively; its effects on the other hemes are small. In contrast, the Y43F mutation hardly changed the potentials. It shows that the aromatic ring at this position contributes to lowering the redox potential of heme 1 locally, although this cannot be the major contribution to the extremely low redox potentials of cytochrome c(3). Furthermore, temperature-dependent line-width broadening in partially reduced samples established that the aromatic ring at position 43 participates in the control of the kinetics of intramolecular electron transfer. The rate of reduction of Y43L cytochrome c(3) by 5-deazariboflavin semiquinone under partially reduced conditions was significantly different from that of the wild type in the last stage of the reduction, supporting the involvement of Tyr43 in regulation of reduction kinetics. The mutation of Y43L, however, did not induce a significant change in the crystal structure.

A novel iron centre in the split-Soret cytochrome c from Desulfovibrio desulfuricans ATCC 27774

The facultative sulfate/nitrate-reducing bacterium Desulfovibrio desulfuricans ATCC 27774 harbours a split-Soret cytochrome c. This cytochrome is a homodimeric protein, having two bis-histidinyl c-type haems per monomer. It has an unique architecture at the haem domain: each haem has one of the coordinating histidines provided by the other monomer, and in each monomer the haems are parallel to each other, almost in van der Waals contact. This work reports the cloning and sequencing of the gene encoding for this cytochrome and shows, by transcriptional analysis, that it is more expressed in nitrate-grown cells than in sulfate-grown ones. In addition, the gene-deduced amino acid sequence revealed two new cysteine residues that could be involved in the binding of a non-haem iron centre. Indeed, the presence of a novel type of an iron-sulfur centre (possibly of the [2Fe-2S] type) was demonstrated by EPR spectroscopy, and putative models for its localization and structure in the cytochrome molecule are proposed on the basis of the so-far-known 3D crystallographic structure of the aerobically purified split-Soret cytochrome, which lacks this centre.

The crystal structure of the hexadeca-heme cytochrome Hmc and a structural model of its complex with cytochrome c(3)

Sulfate-reducing bacteria contain a variety of multi-heme c-type cytochromes. The cytochrome of highest molecular weight (Hmc) contains 16 heme groups and is part of a transmembrane complex involved in the sulfate respiration pathway. We present the 2.42 A resolution crystal structure of the Desulfovibrio vulgaris Hildenborough cytochrome Hmc and a structural model of the complex with its physiological electron transfer partner, cytochrome c(3), obtained by NMR restrained soft-docking calculations. The Hmc is composed of three domains, which exist independently in different sulfate-reducing species, namely cytochrome c(3), cytochrome c(7), and Hcc. The complex involves the last heme at the C-terminal region of the V-shaped Hmc and heme 4 of cytochrome c(3), and represents an example for specific cytochrome-cytochrome interaction.

Production and preliminary characterization of a recombinant triheme cytochrome c(7) from Geobacter sulfurreducens in Escherichia coli

Multiheme cytochromes c have been found in a number of sulfate- and metal ion-reducing bacteria. Geobacter sulfurreducens is one of a family of microorganisms that oxidize organic compounds, with Fe(III) oxide as the terminal electron acceptor. A triheme 9.6 kDa cytochrome c(7) from G. sulfurreducens is a part of the metal ion reduction pathway. We cloned the gene for cytochrome c(7) and expressed it in Escherichia coli together with the cytochrome c maturation gene cluster, ccmABCDEFGH, on a separate plasmid. We designed two constructs, with and without an N-terminal His-tag. The untagged version provided a good yield (up to 6 mg/l of aerobic culture) of the fully matured protein, with all three hemes attached, while the N-terminal His-tag appeared to be detrimental for proper heme incorporation. The recombinant protein (untagged) is properly folded, it has the same molecular weight and displays the same absorption spectra, both in reduced and in oxidized forms, as the protein isolated from G. sulfurreducens and it is capable of reducing metal ions in vitro. The shape parameters for the recombinant cytochrome c(7) determined by small angle X-ray scattering are in good agreement with the ones calculated from a homologous cytochrome c(7) of known structure.

Cloning, sequencing and expression of the tetraheme cytochrome c(3) from Desulfovibrio gigas

The gene encoding the tetraheme cytochrome c(3) from Desulfovibrio gigas was cloned and sequenced from a 2.7-kb EcoRI-PstI insert of D. gigas DNA. The derived amino acid sequence showed that the D. gigas cytochrome c(3) is synthesized as a precursor protein with an N-terminal signal peptide sequence of 25 residues and allowed the correction of the previous reported amino acid sequence (Matias et al. Protein Science 5 (1996) 1342-1354). Expression in D. vulgaris (Hildenborough) was possible by conjugal transfer of a recombinant broad-host-range vector pSUP104 containing a SmaI fragment of the D. gigas cytochrome c(3) gene. Biochemical, immunological and spectroscopic analysis of the purified protein showed that the recombinant cytochrome is identical to that isolated from D. gigas.

Structural basis for the network of functional cooperativities in cytochrome c(3) from Desulfovibrio gigas: solution structures of the oxidised and reduced states

Cytochrome c(3) is a 14 kDa tetrahaem protein that plays a central role in the bioenergetic metabolism of Desulfovibrio spp. This involves an energy transduction mechanism made possible by a complex network of functional cooperativities between redox and redox/protolytic centres (the redox-Bohr effect), which enables cytochrome c(3) to work as a proton activator. The three-dimensional structures of the oxidised and reduced Desulfovibrio gigas cytochrome c(3) in solution were solved using 2D (1)H-NMR data. The reduced protein structures were calculated using INDYANA, an extended version of DYANA that allows automatic calibration of NOE data. The oxidised protein structure, which includes four paramagnetic centres, was solved using the program PARADYANA, which also includes the structural paramagnetic parameters. In this case, initial structures were used to correct the upper and lower volume restraints for paramagnetic leakage, and angle restraints derived from (13)C Fermi contact shifts of haem moiety substituents were used for the axial histidine ligands. Despite the reduction of the NOE intensities by paramagnetic relaxation, the final family of structures is of similar precision and accuracy to that obtained for the reduced form. Comparison of the two structures shows that, although the global folds of the two families of structures are similar, significant localised differences occur upon change of redox state, some of which could not be detected by comparison with the X-ray structure of the oxidised state: (1) there is a redox-linked concerted rearrangement of Lys80 and Lys90 that results in the stabilisation of haem moieties II and III when both molecules are oxidised or both are reduced, in agreement with the previously measured positive redox cooperativity between these two haem moieties. This cooperativity regulates electron transfer, enabling a two-electron step adapted to the function of cytochromes c(3) as the coupling partner of hydrogenase; and (2) the movement of haem I propionate 13 towards the interior of the protein upon reduction explains the positive redox-Bohr effect, establishing the structural basis for the redox-linked proton activation mechanism necessary for energy conservation, driving ATP synthesis.

Expression of a Desulfovibrio tetraheme cytochrome c in Escherichia coli

A tetraheme cytochrome c was successfully overexpressed for the first time in Escherichia coli. Desulfovibrio desulfuricans ATCC 27774 tetraheme cytochrome c(3) was expressed in aerobically grown Escherichia coli cotransformed with Escherichia coli ccm gene cluster (Arslan et al. (1998) Bioch. Biophys. Res. Commun. 251, 744-747). The analysis of the produced cytochrome showed that the signal peptide was correctly cleaved, the four heme groups were inserted and the electronic structure around the heme irons was conserved, i.e., the recombinant tetraheme cytochrome was identical to that isolated from the native source. Contradicting previous results which indicated that Escherichia coli was only capable of producing apocytochrome c(3) (Pollock et al. (1989) J. Gen. Microbiol. 135, 2319-2328), the present work proves unequivocally that the holoform can also be obtained.

Cytochrome c(3) mutants of Desulfovibrio desulfuricans

To explore the physiological role of tetraheme cytochrome c(3) in the sulfate-reducing bacterium Desulfovibrio desulfuricans G20, the gene encoding the preapoprotein was cloned, sequenced, and mutated by plasmid insertion. The physical analysis of the DNA from the strain carrying the integrated plasmid showed that the insertion was successful. The growth rate of the mutant on lactate with sulfate was comparable to that of the wild type; however, mutant cultures did not achieve the same cell densities. Pyruvate, the oxidation product of lactate, served as a poor electron source for the mutant. Unexpectedly, the mutant was able to grow on hydrogen-sulfate medium. These data support a role for tetraheme cytochrome c(3) in the electron transport pathway from pyruvate to sulfate or sulfite in D. desulfuricans G20.

A sequential electron transfer from hydrogenases to cytochromes in sulfate-reducing bacteria

A central step in the energy metabolism of sulfate-reducing bacteria is the oxidation of molecular hydrogen, catalyzed by a periplasmic hydrogenase. The resulting electrons are then transferred to various electron transport chains and used for cytoplasmic sulfate reduction. The complex formation between [NiFeSe] hydrogenase and the soluble periplasmic polyheme cytochromes from Desulfomicrobium norvegicum was characterized by cross-linking experiments, BIAcore and kinetics analysis. Analysis of electron transfer between [NiFeSe] hydrogenase and octaheme cytochrome c(3) (M(r) 26¿ omitted¿000) pointed out that this cytochrome is reduced faster in the presence of catalytic amounts of tetraheme cytochrome c(3) (M(r) 13¿ omitted¿000) isolated from the same organism. The activation of the hydrogenase-dependent reduction of polyheme cytochromes by cytochrome c(3) (M(r) 13¿ omitted¿000), which is now described in both Desulfovibrio and Desulfomicrobium, is proposed as a general mechanism. During this process, cytochrome c(3) (M(r) 13¿ omitted¿000) would act as an electron shuttle in between hydrogenase and the polyheme cytochromes and its conductivity appears to be an important factor.

A proton-NMR investigation of the fully reduced cytochrome c7 from Desulfuromonas acetoxidans. Comparison between the reduced and the oxidized forms

The solution structure via 1H NMR of the fully reduced form of cytochrome c7 has been obtained. The protein sample was kept reduced by addition of catalytic amounts of Desulfovibrio gigas iron hydrogenase in H2 atmosphere after it had been checked that the presence of the hydrogenase did not affect the NMR spectrum. A final family of 35 conformers with rmsd values with respect to the mean structure of 8.7 +/- 1.5 nm and 12.4 +/- 1.3 nm for the backbone and heavy atoms, respectively, was obtained. A highly disordered loop involving residues 54-61 is present. If this loop is ignored, the rmsd values are 6.2 +/- 1.1 nm and 10.2 +/- 1.0 nm for the backbone and heavy atoms, respectively, which represent a reasonable resolution. The structure was analyzed and compared with the already available structure of the fully oxidized protein. Within the indetermination of the two solution structures, the result for the two redox forms is quite similar, confirming the special structural features of the three-heme cluster. A useful comparison can be made with the available crystal structures of cytochromes c3, which appear to be highly homologous except for the presence of a further heme. Finally, an analysis of the factors affecting the reduction potentials of the heme irons was performed, revealing the importance of net charges in differentiating the reduction potential when the other parameters are kept constant.

Sequencing the gene encoding desulfovibrio desulfuricans ATCC 27774 nine-heme cytochrome c

Contradicting early suggestions, the sequencing of the gene encoding the Desulfovibrio desulfuricans (ATCC 27774) nine-heme cytochrome c proves that this cytochrome is not the product of the degradation of the 16-heme containing cytochrome c [Coelho et al. (1996) Acta Cryst. D52, 1202-1208]. However, preliminary data indicate that the cytochrome gene is part of an operon similar to the DvH hmc operon, which contains the gene coding for the 16-heme cytochrome c [Rossi et al. (1993) J. Bacteriol. 175, 4699-4711]. Also, the amino acid sequence deduced from the DNA sequence shows four residues in the C-terminal not predicted in the amino acid sequence obtained by X-ray methods [Matias et al. (1999) Structure 7, 119-130].

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

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

Multiheme cytochromes from the sulfur-reducing bacterium Desulfuromonas acetoxidans

Two new multiheme cytochromes were isolated from the anaerobic sulfur reducing bacterium Desulfuromonas acetoxidans. They have monomeric molecular masses of 50 and 65 kDa and contain six and eight hemes, respectively. Visible and EPR spectroscopies, in the as-isolated (oxidised) cytochromes, show the presence of only low-spin hemes in the 50-kDa cytochrome, and of high-spin and low-spin hemes in the 65-kDa cytochrome. The EPR spectra of the native 65-kDa cytochrome indicate multiple heme-heme interactions, including integer-spin systems as judged by parallel-mode EPR. The 50-kDa cytochrome has a complex redox pattern, as shown by EPR redox titrations, and contains one heme with unusual characteristics. Both cytochromes cover an extremely wide range of reduction potentials, which go from +100 mV to -375 mV for the 50-kDa cytochrome, and +185 mV to -235 mV for the 65-kDa cytochrome. The two cytochromes were tested for hydroxylamine oxidoreductase activity and polysulfide reductase activity, but neither displayed any activity. In contrast, it was found for the first time that the previously characterised cytochrome c551.5, from the same bacterium is very active in the reduction of polysulfide, which suggests that it acts as a terminal reductase in D. acetoxidans.

NMR characterization and solution structure determination of the oxidized cytochrome c7 from Desulfuromonas acetoxidans

The solution structure of the three-heme electron transfer protein cytochrome c7 from Desulfuromonas acetoxidans is reported. The determination of the structure is obtained through NMR spectroscopy on the fully oxidized, paramagnetic form. The richness of structural motifs and the presence of three prosthetic groups in a protein of 68 residues is discussed in comparison with the four-heme cytochromes c3 already characterized through x-ray crystallography. In particular, the orientation of the three hemes present in cytochrome c7 is similar to that of three out of four hemes of cytochromes c3. The reduction potentials of the individual hemes, which have been obtained through the sequence-specific assignment of the heme resonances, are discussed with respect to the properties of the protein matrix. This information is relevant for any attempt to understand the electron transfer pathway.

Structural and functional characterization of cytochrome c3 from D. desulfuricans ATCC 27774 by 1H-NMR

Cooperativity between redox and protonation centres is known to be crucial for the function of complex proteins, but it is often difficult to describe in terms of thermodynamic parameters. Cytochrome c3 is a good model for these studies since, while retaining the overall complexity of larger systems, it is suitable for detailed crystallographic and spectroscopic studies. Assignment of the haem substituent NMR resonances, together with NMR redox titrations of cytochrome c3 from D. desulfuricans ATCC 27774, was used to correlate relative redox potentials to specific haems in the structure: haem II approximately equal to haem I < haem IV < haem III. This order is different from that determined for the homologous proteins studied and in disagreement with that previously reported for this cytochrome (Morais, J., Palma, N., Frazäo, C., Caldeira, J., LeGall, J., Moura, I., Moura, J.J.G. and Carrondo, M.A. (1995) Biochemistry 34, 12830-12841).