NMR structure of the haem core of a novel tetrahaem cytochrome isolated from Shewanella frigidimarina: identification of the haem-specific axial ligands and order of oxidation

A Biochemical Deconstruction-Based Strategy to Assist the Characterization of Bacterial Electric Conductive Filaments

Periplasmic nanowires and electric conductive filaments made of the polymeric assembly of c-type cytochromes from Geobacter sulfurreducens bacterium are crucial for electron storage and/or extracellular electron transfer. The elucidation of the redox properties of each heme is fundamental to the understanding of the electron transfer mechanisms in these systems, which first requires the specific assignment of the heme NMR signals. The high number of hemes and the molecular weight of the nanowires dramatically decrease the spectral resolution and make this assignment extremely complex or unattainable. The nanowire cytochrome GSU1996 (~42 kDa) is composed of four domains (A to D) each containing three c-type heme groups. In this work, the individual domains (A to D), bi-domains (AB, CD) and full-length nanowire were separately produced at natural abundance. Sufficient protein expression was obtained for domains C (~11 kDa/three hemes) and D (~10 kDa/three hemes), as well as for bi-domain CD (~21 kDa/six hemes). Using 2D-NMR experiments, the assignment of the heme proton NMR signals for domains C and D was obtained and then used to guide the assignment of the corresponding signals in the hexaheme bi-domain CD. This new biochemical deconstruction-based procedure, using nanowire GSU1996 as a model, establishes a new strategy to functionally characterize large multiheme cytochromes.

One simple step in the identification of the cofactors signals, one giant leap for the solution structure determination of multiheme proteins

Multiheme proteins play major roles in various biological systems. Structural information on these systems in solution is crucial to understand their functional mechanisms. However, the presence of numerous proton-containing groups in the heme cofactors and the magnetic properties of the heme iron, in particular in the oxidised state, complicates significantly the assignment of the NMR signals. Consequently, the multiheme proteins superfamily is extremely under-represented in structural databases, which constitutes a severe bottleneck in the elucidation of their structural-functional relationships. In this work, we present a strategy that simplifies the assignment of the NMR signals in multiheme proteins and, concomitantly, their solution structure determination, using the triheme cytochrome PpcA from the bacterium Geobacter sulfurreducens as a model. Cost-effective isotopic labeling was used to double label (13C/15N) the protein in its polypeptide chain, with the correct folding and heme post-translational modifications. The combined analysis of 1H-13C HSQC NMR spectra obtained for labeled and unlabeled samples of PpcA allowed a straight discrimination between the heme cofactors and the polypeptide chain signals and their confident assignment. The results presented here will be the foundations to assist solution structure determination of multiheme proteins, which are still very scarce in the literature.

Structural insights into the modulation of the redox properties of two Geobacter sulfurreducens homologous triheme cytochromes

The redox properties of a periplasmic triheme cytochrome, PpcB from Geobacter sulfurreducens, were studied by NMR and visible spectroscopy. The structure of PpcB was determined by X-ray diffraction. PpcB is homologous to PpcA (77% sequence identity), which mediates cytoplasmic electron transfer to extracellular acceptors and is crucial in the bioenergetic metabolism of Geobacter spp. The heme core structure of PpcB in solution, probed by 2D-NMR, was compared to that of PpcA. The results showed that the heme core structures of PpcB and PpcA in solution are similar, in contrast to their crystal structures where the heme cores of the two proteins differ from each other. NMR redox titrations were carried out for both proteins and the order of oxidation of the heme groups was determined. The microscopic properties of PpcB and PpcA redox centers showed important differences: (i) the order in which hemes become oxidized is III-I-IV for PpcB, as opposed to I-IV-III for PpcA; (ii) the redox-Bohr effect is also different in the two proteins. The different redox features observed between PpcB and PpcA suggest that each protein uniquely modulates the properties of their co-factors to assure effectiveness in their respective metabolic pathways. The origins of the observed differences are discussed.

Functional properties of type I and type II cytochromes c3 from Desulfovibrio africanus

Type I cytochrome c(3) is a key protein in the bioenergetic metabolism of Desulfovibrio spp., mediating electron transfer between periplasmic hydrogenase and multihaem cytochromes associated with membrane bound complexes, such as type II cytochrome c(3). This work presents the NMR assignment of the haem substituents in type I cytochrome c(3) isolated from Desulfovibrio africanus and the thermodynamic and kinetic characterisation of type I and type II cytochromes c(3) belonging to the same organism. It is shown that the redox properties of the two proteins allow electrons to be transferred between them in the physiologically relevant direction with the release of energised protons close to the membrane where they can be used by the ATP synthase.

Redox behaviour of the haem domain of flavocytochromec3fromShewanella frigidimarinaprobed by NMR

Flavocytochrome c3 from Shewanella frigidimarina (fcc3) is a tetrahaem periplasmic protein of 64 kDa with fumarate reductase activity. This work reports the first example of NMR techniques applied to the assignment of the thermodynamic order of oxidation of the four individual haems for such large protein, expanding its applicability to a wide range of proteins. NMR data from partially and fully oxidised samples of fcc3 and a mutated protein with an axial ligand of haem IV replaced by alanine were compared with calculated chemical shifts, allowing the structural assignment of the signals and the unequivocal determination of the order of oxidation of the haems. As oxidation progresses the fcc3 haem domain is polarised, with haems I and II much more oxidised than haems III and IV, haem IV being the most reduced. Thus, during catalysis as an electron is taken by the flavin adenosine dinucleotide from haem IV, haem III is eager to re-reduce haem IV, allowing the transfer of two electrons to the active site.

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.

A computational study of Shewanella oneidensis MR-1: structural prediction and functional inference of hypothetical proteins

The genomes of many organisms have been sequenced in the last 5 years. Typically about 30% of predicted genes from a newly sequenced genome cannot be given functional assignments using sequence comparison methods. In these situations three-dimensional structural predictions combined with a suite of computational tools can suggest possible functions for these hypothetical proteins. Suggesting functions may allow better interpretation of experimental data (e.g., microarray data and mass spectroscopy data) and help experimentalists design new experiments. In this paper, we focus on three hypothetical proteins of Shewanella oneidensis MR-1 that are potentially related to iron transport/metabolism based on microarray experiments. The threading program PROSPECT was used for protein structural predictions and functional annotation, in conjunction with literature search and other computational tools. Computational tools were used to perform transmembrane domain predictions, coiled coil predictions, signal peptide predictions, sub-cellular localization predictions, motif prediction, and operon structure evaluations. Combined computational results from all tools were used to predict roles for the hypothetical proteins. This method, which uses a suite of computational tools that are freely available to academic users, can be used to annotate hypothetical proteins in general.

Characterization of the Shewanella oneidensis MR-1 decaheme cytochrome MtrA: expression in Escherichia coli confers the ability to reduce soluble Fe(III) chelates

Shewanella oneidensis MR-1 has the metabolic capacity to grow anaerobically using Fe(III) as a terminal electron acceptor. Growth under these conditions results in the de novo synthesis of a number of periplasmic c-type cytochromes, many of which are multiheme in nature and are thought to be involved in the Fe(III) respiratory process. To begin a biochemical study of these complex cytochromes, the mtrA gene that encodes an approximate 32-kDa periplasmic decaheme cytochrome has been heterologously expressed in Escherichia coli. Co-expression of mtrA with a plasmid that contains cytochrome c maturation genes leads to the assembly of a correctly targeted holoprotein, which covalently binds ten hemes. The recombinant MtrA protein has been characterized by magnetic circular dichroism, which shows that all ten hemes have bis-histidine axial ligation. EPR spectroscopy detected only eight of these hemes, all of which are low spin and provides evidence for a spin-coupled pair of hemes in the oxidized state. Redox titrations of MtrA have been carried out with optical- and EPR-monitored methods, and the hemes are shown to reduce over the potential range -100 to -400 mV. In intact cells of E. coli, MtrA is shown to obtain electrons from the host electron transport chain and pass these onto host oxidoreductases or a range of soluble Fe(III) species. This demonstrates the promiscuous nature of this decaheme cytochrome and its potential to serve as a soluble Fe(III) reductase in intact cells.

Thermodynamic characterization of a tetrahaem cytochrome isolated from a facultative aerobic bacterium, Shewanella frigidimarina: a putative redox model for flavocytochrome c3

The facultative aerobic bacterium Shewanella frigidimarina produces a small c-type tetrahaem cytochrome (86 residues) under anaerobic growth conditions. This protein is involved in the respiration of iron and shares 42% sequence identity with the N-terminal domain of a soluble flavocytochrome, isolated from the periplasm of the same bacterium, which also contains four c -type haem groups. The thermodynamic properties of the redox centres and of an ionizable centre in the tetrahaem cytochrome were determined using NMR and visible spectroscopy techniques. This is the first detailed thermodynamic study performed on a tetrahaem cytochrome isolated from a facultative aerobic bacterium and reveals that this protein presents unique features. The redox centres have negative and different redox potentials, which are modulated by redox interactions between the four haems (covering a range of 8-56 mV) and by redox-Bohr interactions between the haems and an ionizable centre (-4 to -36 mV) located in close proximity to haem III. All of the interactions between the five centres are clearly dominated by electrostatic effects and the microscopic reduction potential of haem III is the one most affected by the oxidation of the other haems and by the protonation state of the molecule. Altogether, this study indicates that the tetrahaem cytochrome isolated from S. frigidimarina (Sfc) has the thermodynamic properties to work as an electron wire between its redox partners. Considering the high degree of sequence identity between Sfc and the cytochrome domain of flavocytochrome c(3), the structural similarities of the haem core, and that the macroscopic potentials are also identical, the results obtained in this work are rationalized in order to put forward a putative redox model for flavocytochrome c(3).

A directional electron transfer regulator based on heme-chain architecture in the small tetraheme cytochrome c from Shewanella oneidensis

The macroscopic and microscopic redox potentials of the four hemes of the small tetraheme cytochrome c from Shewanella oneidensis were determined. The microscopic redox potentials show that the order of reduction is from hemes in the C-terminal domain (hemes 3 and 4) to the N-terminal domain (heme 1), demonstrating the polarization of the tetraheme chain during reduction. This makes heme 4 the most efficient electron delivery site. Furthermore, multi-step reduction of other redox centers through either heme 4 or heme 3 is shown to be possible. This has provided new insights into the two-electron reduction of the flavin in the homologous flavocytochrome c-fumarate reductase.

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.

Determination of the orientation of the axial ligands and of the magnetic properties of the haems in the tetrahaem ferricytochrome from Shewanella frigidimarina

The unambiguous assignment of the nuclear magnetic resonance (NMR) signals of the alpha-substituents of the haems in the tetrahaem cytochrome isolated from Shewanella frigidimarina NCIMB400, was made using a combination of homonuclear and heteronuclear experiments. The paramagnetic (13)C shifts of the nuclei directly bound to the porphyrin of each haem group were analysed in the framework of a model for the haem electronic structure. The analysis yields g-tensors for each haem, which allowed the assignment of some electron paramagnetic resonance (EPR) signals to specific haems, and the orientation of the magnetic axes relative to each haem to be established. The orientation of the axial ligands of the haems was determined semi-empirically from the NMR data, and the structural results were compared with those of the homologous tetrahaem cytochrome from Shewanella oneidensis MR-1 showing significant similarities between the two proteins.

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.

Crystal structures at atomic resolution reveal the novel concept of "electron-harvesting" as a role for the small tetraheme cytochrome c

The genus Shewanella produces a unique small tetraheme cytochrome c that is implicated in the iron oxide respiration pathway. It is similar in heme content and redox potential to the well known cytochromes c(3) but related in structure to the cytochrome c domain of soluble fumarate reductases from Shewanella sp. We report the crystal structure of the small tetraheme cytochrome c from Shewanella oneidensis MR-1 in two crystal forms and two redox states. The overall fold and heme core are surprisingly different from the soluble fumarate reductase structures. The high resolution obtained for an oxidized orthorhombic crystal (0.97 A) revealed several flexible regions. Comparison of the six monomers in the oxidized monoclinic space group (1.55 A) indicates flexibility in the C-terminal region containing heme IV. The reduced orthorhombic crystal structure (1.02 A) revealed subtle differences in the position of several residues, resulting in decreased solvent accessibility of hemes and the withdrawal of a positive charge from the molecular surface. The packing between monomers indicates that intermolecular electron transfer between any heme pair is possible. This suggests there is no unique site of electron transfer on the surface of the protein and that electron transfer partners may interact with any of the hemes, a process termed "electron-harvesting." This optimizes the efficiency of intermolecular electron transfer by maximizing chances of productive collision with redox partners.

Effect of hydrogen-bond networks in controlling reduction potentials in Desulfovibrio vulgaris (Hildenborough) cytochrome C3 probed by site-specific mutagenesis

Cytochromes C3 isolated from Desulfovibrio spp. are periplasmic proteins that play a central role in energy transduction by coupling the transfer of electrons and protons from hydrogenase. Comparison between the oxidized and reduced structures of cytochrome C3 isolated from Desulfovibrio vulgaris (Hildenborough) show that the residue threonine 24, located in the vicinity of heme III, reorients between these two states [Messias, A. C., Kastrau, D. H. W., Costa, H. S., LeGall, J., Turner, D. L., Santos, H., and Xavier, A. V. (1998) J. Mol. Biol. 281, 719-739]. Threonine 24 was replaced with valine by site-directed mutagenesis to elucidate its effect on the redox properties of the protein. The NMR spectra of the mutated protein are very similar to those of the wild type, showing that the general folding and heme core architecture are not affected by the mutation. However, thermodynamic analysis of the mutated cytochrome reveals a large alteration in the microscopic reduction potential of heme III (75 and 106 mV for the protonated forms of the fully reduced and oxidized states, respectively). The redox interactions involving this heme are also modified, while the remaining heme-heme interactions and the redox-Bohr interactions are less strongly affected. Hence, the order of oxidation of the hemes in the mutated cytochrome is different from that in the wild type, and it has a higher overall affinity for electrons. This is consistent with the replacement of threonine 24 by valine preventing the formation of a network of hydrogen bonds, which stabilizes the oxidized state. The mutated protein is unable to perform a concerted two-electron step between the intermediate oxidation stages, 1 and 3, which can occur in the wild-type protein. Thus, replacing a single residue unbalances the global network of cooperativities tuned to control thermodynamically the directionality of the stepwise electron transfer and may affect the functionality of the protein.