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

Electron transfer in a crystalline cytochrome with four hemes

Diffusion of electrons over distances on the order of 100 μm has been observed in crystals of a small tetraheme cytochrome (STC) from Shewanella oneidensis [J. Huang et al. J. Am. Chem. Soc. 142, 10459-10467 (2020)]. Electron transfer between hemes in adjacent subunits of the crystal is slower and more strongly dependent on temperature than had been expected based on semiclassical electron-transfer theory. We here explore explanations for these findings by molecular-dynamics simulations of crystalline and monomeric STC. New procedures are developed for including time-dependent quantum mechanical energy differences in the gap between the energies of the reactant and product states and for evaluating fluctuations of the electronic-interaction matrix element that couples the two hemes. Rate constants for electron transfer are calculated from the time- and temperature-dependent energy gaps, coupling factors, and Franck-Condon-weighted densities of states using an expression with no freely adjustable parameters. Back reactions are considered, as are the effects of various protonation states of the carboxyl groups on the heme side chains. Interactions with water are found to dominate the fluctuations of the energy gap between the reactant and product states. The calculated rate constant for electron transfer from heme IV to heme Ib in a neighboring subunit at 300 K agrees well with the measured value. However, the calculated activation energy of the reaction in the crystal is considerably smaller than observed. We suggest two possible explanations for this discrepancy. The calculated rate constant for transfer from heme I to II within the same subunit of the crystal is about one-third that for monomeric STC in solution.

Physiological function of soluble cytochrome c-552 from alkaliphilic Pseudomonas alcaliphila AL15-21(T)

It has been found that the alkaliphilic Gram-negative bacterium Pseudomonas alcaliphila AL15-21(T) produces a larger amount of soluble c-type cytochromes at pH 10.0 under air-limited condition than at pH 7.0 under high aeration. Cytochrome c-552 was confirmed as the major c-type cytochrome among three soluble c-type cytochromes in the strain. To understand the physiological function of cytochrome c-552, a P. alcaliphila AL15-21(T) cytochrome c-552 gene deletion mutant without a marker gene was constructed by electrotransformation adjusted in this study for the strain. The maximum specific growth rate and maximum cell turbidity of cells grown at pHs 7.0 and 10.0 under the high-aeration condition did not differ significantly between the wild-type and cytochrome c-552 deletion mutant strains. In the mutant grown at pH 10.0 under low-aeration condition, marked decreases in the maximum specific growth rate (40%) and maximum cell turbidity (25%) compared with the wild type were observed. On the other hand, the oxygen consumption rates of cell suspensions of the mutant obtained by the growth at pH 10 under low-aeration condition were slightly higher than that of the wild type. Considering the high electron-retaining ability of cytochrome c-552, the above observations could be accounted for by cytochrome c-552 acting as an electron sink in the periplasmic space. This may facilitate terminal oxidation in the respiratory system at high pH under air-limited conditions.

Mechanisms underlying dioxygen reduction in laccases. Structural and modelling studies focusing on proton transfer

BACKGROUND

Laccases are enzymes that couple the oxidation of substrates with the reduction of dioxygen to water. They are the simplest members of the multi-copper oxidases and contain at least two types of copper centres; a mononuclear T1 and a trinuclear that includes two T3 and one T2 copper ions. Substrate oxidation takes place at the mononuclear centre whereas reduction of oxygen to water occurs at the trinuclear centre.

RESULTS

In this study, the CotA laccase from Bacillus subtilis was used as a model to understand the mechanisms taking place at the molecular level, with a focus in the trinuclear centre. The structures of the holo-protein and of the oxidised form of the apo-protein, which has previously been reconstituted in vitro with Cu(I), have been determined. The former has a dioxygen moiety between the T3 coppers, while the latter has a monoatomic oxygen, here interpreted as a hydroxyl ion. The UV/visible spectra of these two forms have been analysed in the crystals and compared with the data obtained in solution. Theoretical calculations on these and other structures of CotA were used to identify groups that may be responsible for channelling the protons that are needed for reduction of dioxygen to water.

CONCLUSIONS

These results present evidence that Glu 498 is the only proton-active group in the vicinity of the trinuclear centre. This strongly suggests that this residue may be responsible for channelling the protons needed for the reduction. These results are compared with other data available for these enzymes, highlighting similarities and differences within laccases and multicopper oxidases.

Redox-linked conformational changes of a multiheme cytochrome from Geobacter sulfurreducens

Multiheme c-type cytochromes from members of the Desulfovibrionacea and Geobactereacea families play crucial roles in the bioenergetics of these microorganisms. Thermodynamic studies using NMR and visible spectroscopic techniques on tetraheme cytochromes c(3) isolated from Desulfovibrio spp. and more recently on a triheme cytochrome from Geobacter sulfurreducens showed that the properties of each redox centre are modulated by the neighbouring redox centres enabling these proteins to perform energy transduction and thus contributing to cellular energy conservation. Electron/proton transfer coupling relies on redox-linked conformational changes that were addressed for some multiheme cytochromes from the comparison of protein structure of fully reduced and fully oxidised forms. In this work, we identify for the first time in a multiheme cytochrome the simultaneous presence of two different conformations in solution. This was achieved by probing the different oxidation stages of a triheme cytochrome isolated from G. sulfurreducens using 2D-NMR techniques. The results presented here will be the foundations to evaluate the modulation of the redox centres properties by conformational changes that occur during the reoxidation of a multiheme protein.

Crystal structure of the 16 heme cytochrome from Desulfovibrio gigas: a glycosylated protein in a sulphate-reducing bacterium

Sulphate-reducing bacteria have a wide variety of periplasmic cytochromes involved in electron transfer from the periplasm to the cytoplasm. HmcA is a high molecular mass cytochrome of 550 amino acid residues that harbours 16 c-type heme groups. We report the crystal structure of HmcA isolated from the periplasm of Desulfovibrio gigas. Crystals were grown using polyethylene glycol 8K and zinc acetate, and diffracted beyond 2.1 A resolution. A multiple-wavelength anomalous dispersion experiment at the iron absorption edge enabled us to obtain good-quality phases for structure solution and model building. DgHmcA has a V-shape architecture, already observed in HmcA isolated from Desulfovibrio vulgaris Hildenborough. The presence of an oligosaccharide molecule covalently bound to an Asn residue was observed in the electron density maps of DgHmcA and confirmed by mass spectrometry. Three modified monosaccharides appear at the highly hydrophobic vertex, possibly acting as an anchor of the protein to the cytoplasmic membrane.

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.

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.

Structural evidence for a proton transfer pathway coupled with haem reduction of cytochrome c″ from Methylophilus methylotrophus

The crystal structures of the oxidized and reduced forms of cytochrome c" from Methylophilus methylotrophus were solved from X-ray synchrotron data to atomic resolution. The overall fold of the molecule in the two redox states is very similar and is comparable to that of the oxygen-binding protein from the purple phototrophic bacterium Rhodobacter sphaeroides. However, significant modifications occur near the haem group, in particular the detachment from axial binding of His95 observed upon reduction as well as the adoption of different conformations of some protonatable residues that form a possible proton path from the haem pocket to the protein surface. These changes are associated with the previously well characterized redox-Bohr behaviour of this protein. Furthermore they provide a model for one of the presently proposed mechanisms of proton translocation in the much more complex protein cytochrome c oxidase.

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.

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.

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.

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

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