Conformational component in the coupled transfer of multiple electrons and protons in a monomeric tetraheme cytochrome

Structural and functional insights of GSU0105, a unique multiheme cytochrome from G. sulfurreducens

Geobacter sulfurreducens possesses over 100 cytochromes that assure an effective electron transfer to the cell exterior. The most abundant group of cytochromes in this microorganism is the PpcA family, composed of five periplasmic triheme cytochromes with high structural homology and identical heme coordination (His-His). GSU0105 is a periplasmic triheme cytochrome synthetized by G. sulfurreducens in Fe(III)-reducing conditions but is not present in cultures grown on fumarate. This cytochrome has a low sequence identity with the PpcA family cytochromes and a different heme coordination, based on the analysis of its amino acid sequence. In this work, amino acid sequence analysis, site-directed mutagenesis, and complementary biophysical techniques, including ultraviolet-visible, circular dichroism, electron paramagnetic resonance, and nuclear magnetic resonance spectroscopies, were used to characterize GSU0105. The cytochrome has a low percentage of secondary structural elements, with features of α-helices and β-sheets. Nuclear magnetic resonance shows that the protein contains three low-spin hemes (Fe(II), S = 0) in the reduced state. Electron paramagnetic resonance shows that, in the oxidized state, one of the hemes becomes high-spin (Fe(III), S = 5/2), whereas the two others remain low-spin (Fe(III), S = 1/2). The data obtained also indicate that the heme groups have distinct axial coordination. The apparent midpoint reduction potential of GSU0105 (-154 mV) is pH independent in the physiological range. However, the pH modulates the reduction potential of the heme that undergoes the low- to high-spin interconversion. The reduction potential values of cytochrome GSU0105 are more distinct compared to those of the PpcA family members, providing the protein with a larger functional working redox potential range. Overall, the results obtained, together with an amino acid sequence analysis of different multiheme cytochrome families, indicate that GSU0105 is a member of a new group of triheme cytochromes.

A Heme Propionate Staples the Structure of Cytochrome c for Methionine Ligation to the Heme Iron

Ligand-switch reactions at the heme iron are common in biological systems, but their mechanisms and the features of the polypeptide fold that support dual ligation are not well understood. In cytochrome c (cyt c), two low-stability loops (Ω-loop C and Ω-loop D) are connected by the heme propionate HP6. At alkaline pH, the native Met80 ligand from Ω-loop D switches to a Lys residue from the same loop. Deprotonation of an as yet unknown group triggers the alkaline transition. We have created the two cyt c variants T49V/K79G and T78V/K79G with altered connections of these two loops to HP6. Electronic absorption, NMR, and EPR studies demonstrate that at pH 7.4 ferric forms of these variants are Lys-ligated, whereas ferrous forms maintain the native Met80 ligation. Measurements of protein stability, cyclic voltammetry, pH-jump and gated electron-transfer kinetics have revealed that these Thr to Val substitutions greatly affect the alkaline transition in both ferric and ferrous proteins. The substitutions modify the stability of the Met-ligated species and reduction potentials of the heme iron. The kinetics of ligand-switch processes are also altered, and analyses of these effects implicate redox-dependent differences in metal-ligand interactions and the role of the protein dynamics, including cross-talk between the two Ω-loops. With the two destabilized variants, it is possible to map energy levels for the Met- and Lys-ligated species in both ferric and ferrous proteins and assess the role of the protein scaffold in redox-dependent preferences for these two ligands. The estimated shift in the heme iron reduction potential upon deprotonation of the "trigger" group is consistent with those associated with deprotonation of an HP, suggesting that HP6, on its own or as a part of a hydrogen-bonded cluster, is a likely "trigger" for the Met to Lys ligand switch.

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.

The role of intramolecular interactions in the functional control of multiheme cytochromes c

Detailed thermodynamic and structural data measured in soluble monomeric multiheme cytochromes c provided the basis to investigate the functional significance of interactions between redox co-factors. The steep decay of intramolecular interactions with distance means that close proximity of the redox centers is necessary to modulate the intrinsic reduction potentials in a significant way. This ensures selection of specific populations during redox activity in addition to maintaining fast intramolecular electron transfer. Therefore, intramolecular interactions between redox co-factors play an important role in establishing the biological function of the protein by controlling how electrons flow through and are distributed among the co-factors.

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.

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.

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.

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.

On the Use of Different Dielectric Constants for Computing Individual and Pairwise Terms in Poisson−Boltzmann Studies of Protein Ionization Equilibrium

Poisson-Boltzmann (PB) models are a fast and common tool for studying electrostatic processes in proteins, particularly their ionization equilibrium (protonation and/or reduction), often yielding quite good results when compared with more detailed models. Yet, they are conceptually very simple and necessarily approximate, their empirical character being most evident when it comes to the choice of the dielectric constant assigned to the protein region. The present study analyzes several factors affecting the ability of PB-based methods to model protein ionization equilibrium. We give particular attention to a suggestion made by Warshel and co-workers (e.g., Sham et al. J. Phys. Chem. B 1997, 101, 4458) of using different protein dielectric constants for computing the individual (site) and the pairwise (site-site) terms of the ionization free energies. Our prediction of pK(a) values for several proteins indicates that no advantage is obtained by such a procedure, even for sites that are buried and/or display large pK(a) shifts relative to the solution values. In particular, the present methodology gives the best predictions using a dielectric constant around 20, for shifted/buried and nonshifted/exposed sites alike. The similarities and differences between the PB model and Warshel's PDLD/S model are discussed, as well as the reasons behind their apparently discrepant results. The present PB model is shown to predict also good reduction potentials in redox proteins.

The Protonation State of a Heme Propionate Controls Electron Transfer in CytochromecOxidase†

In cytochrome c oxidase (CcO), exergonic electron transfer reactions from cytochrome c to oxygen drive proton pumping across the membrane. Elucidation of the proton pumping mechanism requires identification of the molecular components involved in the proton transfer reactions and investigation of the coupling between internal electron and proton transfer reactions in CcO. While the proton-input trajectory in CcO is relatively well characterized, the components of the output pathway have not been identified in detail. In this study, we have investigated the pH dependence of electron transfer reactions that are linked to proton translocation in a structural variant of CcO in which Arg481, which interacts with the heme D-ring propionates in a proposed proton output pathway, was replaced with Lys (RK481 CcO). The results show that in RK481 CcO the midpoint potentials of hemes a and a(3) were lowered by approximately 40 and approximately 15 mV, respectively, which stabilizes the reduced state of Cu(A) during reaction of the reduced CcO with O(2). In addition, while the pH dependence of the F --> O rate in wild-type CcO is determined by the protonation state of two protonatable groups with pK(a) values of 6.3 and 9.4, only the high-pK(a) group influences this rate in RK481 CcO. The results indicate that the protonation state of the Arg481 heme a(3) D-ring propionate cluster having a pK(a) of approximately 6.3 modulates the rate of internal electron transfer and may act as an acceptor of pumped protons.

Binding of ligands originates small perturbations on the microscopic thermodynamic properties of a multicentre redox protein

NMR and visible spectroscopy coupled to redox measurements were used to determine the equilibrium thermodynamic properties of the four haems in cytochrome c3 under conditions in which the protein was bound to ligands, the small anion phosphate and the protein rubredoxin with the iron in the active site replaced by zinc. Comparison of these results with data for the isolated cytochrome shows that binding of ligands causes only small changes in the reduction potentials of the haems and their pairwise interactions, and also that the redox-sensitive acid-base centre responsible for the redox-Bohr effect is essentially unaffected. Although neither of the ligands tested is a physiological partner of cytochrome c3, the small changes observed for the thermodynamic properties of cytochrome c3 bound to these ligands vs. the unbound state, indicate that the thermodynamic properties measured for the isolated protein are relevant for a physiological interpretation of the role of this cytochrome in the bioenergetic metabolism of Desulfovibrio.

Electric-field-induced redox potential shifts of tetraheme cytochromes c3 immobilized on self-assembled monolayers: surface-enhanced resonance Raman spectroscopy and simulation studies

The tetraheme protein cytochrome c(3) (Cyt-c(3)) from Desulfovibrio gigas, immobilized on a self-assembled monolayer (SAM) of 11-mercaptoundecanoic acid, is studied by theoretical and spectroscopic methods. Molecular dynamics simulations indicate that the protein docks to the negatively charged SAM via its lysine-rich domain around the exposed heme IV. Complex formation is associated with only little protein structural perturbations. This finding is in line with the resonance Raman and surface-enhanced resonance Raman (SERR) spectroscopic results that indicate essentially the same heme pocket structures for the protein in solution and adsorbed on SAM-coated Ag electrodes. Electron- and proton-binding equilibrium calculations reveal substantial negative shifts of the redox potentials compared to the protein in solution. The magnitude of these shifts decreases in the order heme IV (-161 mV) > heme III (-73 mV) > heme II (-57 mV) > heme I (-26 mV), resulting in a change of the order of reduction. These shifts originate from the distance-dependent electrostatic interactions between the SAM headgroups and the individual hemes, leading to a stabilization of the oxidized forms. The results of the potential-dependent SERR spectroscopic analyses are consistent with the theoretical predictions and afford redox potential shifts of -160 mV (heme IV), -90 mV (heme III), -70 mV (heme II), and +20 mV (heme I) relative to the experimental redox potentials for Cyt-c(3) in solution. SERR spectroscopic experiments reveal electric-field-induced changes of the redox potentials also for the structurally very similar Cyt-c(3) from Desulfovibrio vulgaris, although the shifts are somewhat smaller compared to Cyt-c(3) from D. gigas. This study suggests that electric-field-induced redox potential shifts may also occur upon binding to biomembranes or partner proteins and thus may affect biological electron transfer processes.

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.

Distance dependence of interactions between charged centres in proteins with common structural features

Data collected for interactions among redox centres, and interactions between redox centres and acid-base residues in a family of small multihaem cytochromes are analysed. The distance dependent attenuation of the interactions between non-surface charges, for separations that range from 8 to 23 angstroms, can be described by a simple function derived from the Debye-Huckel formalism, fit to 9.5 and 7.6 as values for the relative dielectric constant and Debye length, respectively. However, there is considerable scatter in the data despite the structural similarities among the proteins, which is discussed in the framework of using such simple models in predicting properties of novel proteins.

Thermodynamic and choreographic constraints for energy transduction by cytochrome c oxidase

Cooperative effects are fundamental for electroprotonic energy transduction processes, crucial to sustain much of life chemistry. However, the primary cooperative mechanism by which transmembrane proteins couple the downhill transfer of electrons to the uphill activation (acidification) of protic groups is still a matter of great controversy. To understand cooperative processes fully, it is necessary to obtain the microscopic thermodynamic parameters of the functional centres and relate them to the relevant structural features, a task difficult to achieve for large proteins. The approach discussed here explores how this may be done by extrapolation from mechanisms used by simpler proteins operative in similar processes. The detailed study of small, soluble cytochromes performing electroprotonic activation has shown how they use anti-electrostatic effects to control the synchronous movement of charges. These include negative e(-)/H(+) (redox-Bohr effect) cooperativities. This capacity is the basis to discuss an unorthodox mechanism consistent with the available experimental data on the process of electroprotonic energy transduction performed by cytochrome c oxidase (CcO).

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.

Investigation of protonatable residues in Rhodothermus marinus caa3 haem-copper oxygen reductase: comparison with Paracoccus denitrificans aa3 haem-copper oxygen reductase

The Rhodothermus marinus caa(3 )haem-copper oxygen reductase contains all the residues of the so-called D- and K-proton channels, with the notable exception of the helix VI glutamate residue (Glu278(I) in Paracoccus denitrificans aa(3)), being nevertheless a true oxygen reductase reducing O(2) to water, and an efficient proton pump. Instead, in the same helix, but one turn below, it has a tyrosine residue (Tyr256(I), R. marinus caa(3) numbering), whose hydroxyl group occupies the same spatial position as the carboxylate group of Glu278(I), as deduced by comparative modelling techniques. Therefore, we proposed previously that this tyrosine residue could play an important role in the proton pathway. In this article we further study this hypothesis, by investigating the equilibrium thermodynamics of protonation in R. marinus caa(3), using theoretical methodologies based on the structural model previously obtained. Control calculations are also performed for the P. denitrificans aa(3) oxygen reductase. In both oxygen reductases we find several residues that are proton active (i.e., that display partial protonation) at physiological pH, some of them being redox sensitive (i.e., sensitive to the protein redox state). However, the caa(3 )Tyr256(I) is not proton active at physiological pH, in contrast to the aa(3) Glu278(I) which is both proton active at physiological pH and shows a high redox sensitivity. In R. marinus caa(3) we do not find any other residues in the same protein zone that can have this property. Therefore, there are no putative D-channel residues that are proton active in this oxidase. The protonatable residues of the K-channel are much more functionally conserved in both oxygen reductases than the same type of residues in the D-channel. Two (Tyr262(I) and Lys336(I), caa(3) numbering) out of three protonatable K-channel residues are proton active and redox sensitive in both proteins.

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

The tetraheme cytochrome c3 is a small metalloprotein with ca. 13,000 Da found in sulfate-reducing bacteria, which is believed to act as a partner of hydrogenase. The three-dimensional structure of the oxidized and reduced forms of cytochrome c3 from Desulfovibrio desulfuricans ATCC 27774 at pH 7.6 were determined using high-resolution X-ray crystallography and were compared with the previously determined oxidized form at pH 4.0. Theoretical calculations were performed with both structures, using continuum electrostatic calculations and Monte Carlo sampling of protonation and redox states, in order to understand the molecular basis of the redox-Bohr and cooperativity effects related to the coupled transfer of electrons and protons. We were able to identify groups that showed redox-linked conformational changes. In particular, Glu61, His76, and propionate D of heme II showed important contributions to the redox-cooperativity, whereas His76, propionate A of heme I, and propionate D of heme IV were the key residues for the redox-Bohr effect. Upon reduction, an important movement of the backbone region surrounding hemes I and II was also identified, that, together with a few redox-linked conformational changes in side-chain residues, results in a significant decrease in the solvent accessibility of hemes I and II.

Redox-Bohr and other cooperativity effects in the nine-heme cytochrome C from Desulfovibrio desulfuricans ATCC 27774: crystallographic and modeling studies

The nine-heme cytochrome c is a monomeric multiheme cytochrome found in Desulfovibrio desulfuricans ATCC 27774. The polypeptide chain comprises 296 residues and wraps around nine hemes of type c. It is believed to take part in the periplasmic assembly of proteins involved in the mechanism of hydrogen cycling, receiving electrons from the tetraheme cytochrome c3. With the purpose of understanding the molecular basis of electron transfer processes in this cytochrome, we have determined the crystal structures of its oxidized and reduced forms at pH 7.5 and performed theoretical calculations of the binding equilibrium of protons and electrons in these structures. This integrated study allowed us to observe that the reduction process induced relevant conformational changes in several residues, as well as protonation changes in some protonatable residues. In particular, the surroundings of hemes I and IV constitute two areas of special interest. In addition, we were able to ascertain the groups involved in the redox-Bohr effect present in this cytochrome and the conformational changes that may underlie the redox-cooperativity effects on different hemes. Furthermore, the thermodynamic simulations provide evidence that the N- and C-terminal domains function in an independent manner, with the hemes belonging to the N-terminal domain showing, in general, a lower redox potential than those found in the C-terminal domain. In this way, electrons captured by the N-terminal domain could easily flow to the C-terminal domain, allowing the former to capture more electrons. A notable exception is heme IX, which has low redox potential and could serve as the exit path for electrons toward other proteins in the electron transfer pathway.

Understanding the mechanism of proton movement linked to oxygen reduction in cytochrome c oxidase: lessons from other proteins

Cytochrome c oxidase is a large intrinsic membrane protein designed to use the energy of electron transfer and oxygen reduction to pump protons across a membrane. The molecular mechanism of the energy conversion process is not understood. Other proteins with simpler, better resolved structures have been more completely defined and offer insight into possible mechanisms of proton transfer in cytochrome c oxidase. Important concepts that are illustrated by these model systems include the ideas of conformational change both close to and at a distance from the triggering event, and the formation of a transitory water-linked proton pathway during a catalytic cycle. Evidence for the applicability of these concepts to cytochrome c oxidase is discussed.

A mechano-chemical model for energy transduction in cytochrome c oxidase: the work of a Maxwell's god

Cytochrome c3 has a central role in the energetics of Desulfovibrio sp., where it performs an electroprotonic energy transduction step. This process uses a network of cooperativities, largely based on anti-Coulomb components, resulting from a mechano-chemical energy coupling mechanism. This mechanism provides a model coherent with the data available for the redox chemistry of haem a of cytochrome c oxidase and its link to the activation of protons. A crucial feature of the model is an anti-Coulomb effect that sets the stage for a molecular ratchet, ensuring vectoriality for the redox-driven localised movement of protons across the membrane, against an electrochemical gradient.

Sulfate respiration in Desulfovibrio vulgaris Hildenborough. Structure of the 16-heme cytochrome c HmcA AT 2.5-A resolution and a view of its role in transmembrane electron transfer

The crystal structure of the high molecular mass cytochrome c HmcA from Desulfovibrio vulgaris Hildenborough is described. HmcA contains the unprecedented number of sixteen hemes c attached to a single polypeptide chain, is associated with a membrane-bound redox complex, and is involved in electron transfer from the periplasmic oxidation of hydrogen to the cytoplasmic reduction of sulfate. The structure of HmcA is organized into four tetraheme cytochrome c(3)-like domains, of which the first is incomplete and contains only three hemes, and the final two show great similarity to the nine-heme cytochrome c from Desulfovibrio desulfuricans. An isoleucine residue fills the vacant coordination space above the iron atom in the five-coordinated high-spin Heme 15. The characteristics of each of the tetraheme domains of HmcA, as well as its surface charge distribution, indicate this cytochrome has several similarities with the nine-heme cytochrome c and the Type II cytochrome c(3) molecules, in agreement with their similar genetic organization and mode of reactivity and further support an analogous physiological function for the three cytochromes. Based on the present structure, the possible electron transfer sites between HmcA and its redox partners (namely Type I cytochrome c(3) and other proteins of the Hmc complex), as well as its physiological role, are discussed.

Molecular dynamics study of Desulfovibrio africanus cytochrome c3 in oxidized and reduced forms

A 5-ns molecular dynamics study of a tetraheme cytochrome in fully oxidized and reduced forms was performed using the CHARMM molecular modeling program, with explicit water molecules, Langevin dynamics thermalization, Particle Mesh Ewald long-range electrostatics, and quantum mechanical determination of heme partial charges. The simulations used, as starting points, crystallographic structures of the oxidized and reduced forms of the acidic cytochrome c(3) from Desulfovibrio africanus obtained at pH 5.6. In this paper we also report structures for the two forms obtained at pH 8. In contrast to previous cytochrome c(3) dynamics simulations, our model is stable. The simulation structures agree reasonably well with the crystallographic ones, but our models show higher flexibility and the water molecules are more labile. We have compared in detail the differences between the simulated and experimental structures of the two redox states and observe that the hydration structure is highly dependent on the redox state. We have also analyzed the interaction energy terms between the hemes, the protein residues, and water. The direct electrostatic interaction between hemes is weak and nearly insensitive to the redox state, but the remaining terms are large and contribute in a complex way to the overall potential energy differences that we see between the redox states.

Thermodynamic and kinetic characterization of trihaem cytochrome c3 from Desulfuromonas acetoxidans

Trihaem cytochrome c3 (also known as cytochrome c551.5 and cytochrome c7) is isolated from the periplasmic space of Desulfuromonas acetoxidans, a sulfur-reducing bacterium. Thermodynamic and kinetic data for the trihaem cytochrome c3 are presented and discussed in the context of the possible physiological implications of its functional properties with respect to the natural habitat of D. acetoxidans, namely as a symbiont with green sulfur bacteria working as a mini-sulfuretum. The thermodynamic properties were determined through the fit of redox titration data, followed by NMR and visible spectroscopy, to a model of four functional centres that describes the network of cooperativities between the three haems and one protolytic centre. The kinetics of trihaem cytochrome c3 reduction by sodium dithionite were studied using the stopped-flow technique and the data were fitted to a kinetic model that makes use of the thermodynamic properties to obtain the rate constants of the individual haems. This analysis indicates that the electrons enter the cytochrome mainly via haem I. The reduction potentials of the haems in this cytochrome show little variation with pH within the physiological range, and the kinetic studies show that the rates of reduction are also independent of pH in the range studied. Thus, although the trihaem cytochrome c3 is readily reduced by hydrogenases from Desulfovibrio sp. and its haem core is similar to that of the homologous tetrahaem cytochromes c3, its physico-chemical properties are quite different, which suggests that these multihaem cytochromes with similar structures perform different functions.

Influence of structure, pH and membrane potential on proton movement in cytochrome oxidase

Cytochrome c oxidase (CcO) reconstituted into phospholipid vesicles and subject to a membrane potential, exhibits different characteristics than the free enzyme, with respect to effects of mutations, pH, inhibitors, and native structural differences between CcO from different species. The results indicate that the membrane potential influences the conformation of CcO and the direction of proton movement in the exit path. The importance of the protein structure above the hemes in proton exit, back leak and respiratory control is discussed.

Plasticity of proton pathways in haem-copper oxygen reductases

Oxygen reductases are the final enzymes in the aerobic respiratory chains catalysing the reduction of dioxygen to water, with the concomitant translocation of protons through the bacterial cytoplasmatic or mitochondrial membranes. Most of these enzymes belong to the family of haem-copper oxygen reductases. Intraprotein proton-conducting pathways are needed for the chemical reaction and for the translocated protons. Based on sequence and structural analyses, and site-directed mutagenesis, two proton channels were established for the mitochondrial-like oxygen reductases. However, the amino acid residues forming these channels are not conserved among the family members. Most importantly, many oxygen reductases do not contain ionisable amino acid residues in the putative proton pathways nor in alternative positions. The diversity of channels in haem-copper oxygen reductases exemplifies the plasticity of proton pathways that occurred throughout evolution, and strongly suggests a substantial role for water as the main proton carrier.

Redox-coupled conformational alternations in cytochrome c(3) from D. vulgaris Miyazaki F on the basis of its reduced solution structure

Heteronuclear NMR spectroscopy was performed to determine the solution structure of (15)N-labeled ferrocytochrome c(3) from Desulfovibrio vulgaris Miyazaki F (DvMF). Although the folding of the reduced cytochrome c(3) in solution was similar to that of the oxidized one in the crystal structure, the region involving hemes 1 and 2 was different. The redox-coupled conformational change is consistent with the reported solution structure of D. vulgaris Hildenborough ferrocytochrome c(3), but is different from those of other cytochromes c(3). The former is homologous with DvMF cytochrome c(3) in amino acid sequence. Small displacements of hemes 1 and 2 relative to hemes 3 and 4 were observed. This observation is consistent with the unusual behavior of the 2(1)CH(3) signal of heme 3 reported previously. As shown by the (15)N relaxation parameters of the backbone, a region between hemes 1 and 2 has more flexibility than the other regions. The results of this work strongly suggest that the cooperative reduction of hemes 1 and 2 is based on the conformational changes of the C-13 propionate of heme 1 and the aromatic ring of Tyr43, and the interaction between His34 and His 35 through covalent and coordination bonds.