Energy transduction through electron/proton cooperativity is at the heart of the metabolism of every living organism Nonetheless, the search for the structural bases sustaining these phenomena has been hindered by the fact that they are usually associated with complex transmembrane proteins of high molecular weight.

Nuclear Magnetic Resonance (NMR) spectroscopy is extremely powerful to study distinct biological systems ranging from biomolecules to specific metabolites. This chapter presents the basic concepts of the technique and illustrates its potential to study such systems. Similarly, to other spectroscopic techniques, the theoretical background of NMR is sustained by detailed mathematics and physical chemistry concepts, which were kept to the minimum. The intent is to introduce the fundamentals of the technique to science students from different backgrounds. The basic concepts of NMR spectroscopy are briefly presented in the first section, and the following sections describe applications in the biosciences field, using electron transfer proteins as model, particularly cytochromes. The heme groups endow cytochromes with particular features making them excellent examples to illustrate the high versatility of NMR spectroscopy. The main methodologies underlying protein solution structure determination are discussed in the second section. This is followed by a description of the main experiments explored to structurally map protein-protein or protein-ligand interface regions in molecular complexes. Finally, it is shown how NMR spectroscopy can assist in the functional characterization of multiheme cytochromes.

The thermoacidophilic archaeon Acidianus ambivalens contains a monomeric 47 kDa type-II NADH dehydrogenase (NDH), which contains a covalently bound flavin. In this work, by a combination of several methods, namely 31P-nuclear magnetic resonance and fluorescence spectroscopies, it is proven that this enzyme contains covalent FMN, a novelty among this family of enzymes, which were so far thought to mainly have the flavin dinucleotide form. Discrimination between several possible covalent flavin linkages was achieved by spectral and fluorescence experiments, which identified an 8α-N(1)-histidylflavin-type of linkage. Analysis of the gene-deduced amino acid sequence of type-II NDH showed no transmembranar helices and allowed the definition of putative dinucleotide and quinone binding motifs. Further, it is suggested that membrane anchoring can be achieved via amphipatic helices.

Using 2D-NMR the four haems of Desulfovibrio vulgaris (Hildenborough) cytochromes, within the X-ray structure were fully cross-assigned according to their redox potential. The strategy used was based on a complete network of chemical exchange connectivities between the NMR signals obtained for all oxidation levels to the corresponding ones in the fully reduced spectrum [1992, Eur. J. Biochem., in press]. This unequivocal cross-assignment disagrees within earlier results obtained for the similar protein from Desulfovibrio vulgaris (Miyazaki F.) [1991, FEBS Lett. 285, 149–151]

Microbial electrosynthesis is an emerging green technology that explores the capability of a particular group of microorganisms to drive their metabolism toward the production of hydrogen or value-added chemicals from electrons supplied by electrode surfaces. The cytochrome PccH showed the largest increase in transcription when electrons are supplied to Geobacter sulfurreducens biofilms. Gene knock-out experiments have shown that the electron transfer toward G. sulfurreducens cells was completely inhibited by the deletion of the gene encoding for cytochrome PccH. This identifies a crucial role for this protein in G. sulfurreducens microbial electrosynthesis mechanisms, which are currently unknown. In this work, we present the backbone (1H, 13C and 15N) and heme assignment for PccH in the oxidized state. The data obtained paves the way to identify and structurally map the molecular interaction regions between the cytochrome PccH and its physiological redox partners.

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.

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.

Three-quarters of the carbon-13 resonances of nuclei attached to the four haems of Desulfovibrio vulgaris ferricytochrome c3 are assigned. Preliminary analysis of their Fermi contact interactions shows that the shifts are directly related to the orientation of both of the axial histidine ligands in each case and the approach can therefore be used to obtain structural information in other cytochromes with bis-histidinyl coordination. The implications for the control of redox potential in cytochromes are discussed.

Electron harvesting bacteria are key targets to develop microbial electrosynthesis technologies, which are valid alternatives for the production of value-added compounds without utilization of fossil fuels. Geobacter sulfurreducens, that is capable of donating and accepting electrons from electrodes, is one of the most promising electroactive bacteria. Its electron transfer mechanisms to electrodes have been progressively elucidated, however the electron harvesting pathways are still poorly understood. Previous studies showed that the periplasmic cytochromes PccH and GSU2515 are overexpressed in current-consuming G. sulfurreducens biofilms. PccH was characterized, though no putative partners have been identified. In this work, GSU2515 was characterized by complementary biophysical techniques and in silico simulations using the AlphaFold neural network. GSU2515 is a low-spin monoheme cytochrome with a disordered N-terminal region and an α-helical C-terminal domain harboring the heme group. The cytochrome undergoes a redox-linked heme axial ligand switch, with Met91 and His94 as distal axial ligand in the reduced and oxidized state, respectively. The reduction potential of the cytochrome is negative and is modulated by the pH in the physiological range: -78 mV at pH 6 and -113 mV at pH 7. Such pH-dependence coupled to the redox-linked switch of the axial ligand allows the cytochrome to drive a proton-coupled electron transfer step that is crucial to confer directionality to the respiratory chain. Biomolecular interactions and electron transfer experiments indicated that GSU2515 and PccH form a redox complex. Overall, the data obtained highlights for the first time how periplasmic proteins bridge the electron transfer between the outer and inner membrane in the electron harvesting pathways of G. sulfurreducens.

The rising interest in the use of Geobacter bacteria for biotechnological applications demands a deep understanding of how these bacteria are able to thrive in a variety of environments and perform extracellular electron transfer. The Geobacter metallireducens bacterium can couple the oxidation of a wide range of compounds to the reduction of several extracellular acceptors, including heavy metals, toxic organic compounds or electrode surfaces. The periplasmic c-type cytochrome PpcA from this bacterium is a member of a family composed of five periplasmic triheme cytochromes, which are important to bridge the electron transfer between the cytoplasm and the extracellular environment. To better understand the functional mechanism of PpcA it is essential to obtain structural data for this cytochrome. In this work, the geometry of the heme axial ligands, as well as the magnetic properties of the hemes were determined for the oxidized form of the cytochrome, using the 13C NMR chemical shifts of the heme α-substituents. The results were further compared with those previously obtained for the homologous cytochrome from Geobacter sulfurreducens. The orientations of the axial histidine planes and the magnetic properties of the hemes are conserved in both proteins. Overall, the results obtained allowed the definition of the orientation of the magnetic axes of PpcA from G. metallireducens, which will be used as constraints to assist the solution structure determination of the cytochrome in the oxidized form.

The unambiguous assignment of the nuclear magnetic resonance (NMR) signals of the α-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 13C 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.

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.

NMR spectroscopy has been applied with great success to study electron transfer proteins
with multiple redox centers. This study aimed to elucidate the redox behavior the enzyme fumarate
reductase from Shewanella frigidimarina and particularly to reveal the electron transfer mechanism
from the N-terminal domain to the active center. We developed a new strategy encompassing the
acquisition of 1H-EXSY bidimensional spectra for observation of chemical exchange connectivities in
partially oxidized samples of fcc3, estimation of the paramagnetic chemical shifts expected for the
heme substituents and their comparison with NMR spectra obtained in the fully oxidized protein. This
study allowed obtaining the order of oxidation of the different groups (II-I-III, IV) and gave insights of
the functional mechanisms that allow fcc3 to efficiently transfer electrons from the N-terminal domain
to the active center.

The cytochrome PccH from Geobacter sulfurreducens (Gs) plays a crucial role in current-consuming fumarate-reducing biofilms. Deletion of pccH gene inhibited completely electron transfer from electrodes toward Gs cells. The pccH gene was cloned and the protein heterologously expressed in Escherichia coli. Complementary biophysical techniques including CD, UV-visible and NMR spectroscopy were used to characterize PccH. This cytochrome contains one low-spin c-type heme with His-Met axial coordination and unusual low-reduction potential. This reduction potential is pH-dependent, within the Gs physiological pH range, and is discussed within the context of the electron transfer mechanisms from electrodes to Gs cells.

The thermodynamic parameters which govern the homotropic (e−/e−) and heterotropic (e−/H+) cooperativity in the tetrahaem cytochrome c3 isolated from Desulfovibrio vulgaris (Hildenborough) were determined, using the paramagnetic shifts of haem methyl groups in the NMR spectra of intermediate oxidized states at different pH levels. A model is put forward to explain how the network of positive and negative cooperativities between the four haems and acid/base group(s) enables the protein to achieve a proton-assisted 2e− step.

The monoheme outer membrane cytochrome F (OmcF) from Geobacter sulfurreducens plays an important role in Fe(III) reduction and electric current production. The electrochemical characterization of this cytochrome has shown that its redox potential is modulated by the solution pH (redox-Bohr effect) endowing the protein with the necessary properties to couple electron and proton transfer in the physiological range. The analysis of the OmcF structures in the reduced and oxidized states showed that with the exception of the side chain of histidine 47 (His⁴⁷), all other residues with protonatable side chains are distant from the heme iron and, therefore, are unlikely to affect the redox potential of the protein. The protonatable site at the imidazole ring of His⁴⁷ is in the close proximity to the heme and, therefore, this residue was suggested as the redox-Bohr center. In the present work, we tested this hypothesis by replacing the His⁴⁷ with non-protonatable residues (isoleucine – OmcFH47I and phenylalanine – OmcFH47F). The structure of the mutant OmcFH47I was determined by X-ray crystallography to 1.13 Å resolution and showed only minimal changes at the site of the mutation. Both mutants were ¹⁵N-labeled and their overall folding was confirmed to be the same as the wild-type by NMR spectroscopy. The pH dependence of the redox potential of the mutants was measured by cyclic voltammetry. Compared to the wild-type protein, the magnitude of the redox-Bohr effect in the mutants was smaller, but not fully abolished, confirming the role of His⁴⁷ on the pH modulation of OmcF’s redox potential. However, the pH effect on the heme substituents’ NMR chemical shifts suggested that the heme propionate P₁₃ also contributes to the overall redox-Bohr effect in OmcF. In physiological terms, the contribution of two independent acid–base centers to the observed redox-Bohr effect confers OmcF a higher versatility to environmental changes by coupling electron/proton transfer within a wider pH range.

Electroactive Geobacter bacteria can perform extracellular electron transfer and present a wide metabolic versatility. These bacteria reduce organic, toxic and radioactive compounds, and produce electric current while interacting with electrodes, making them interesting targets for numerous biotechnological applications. Their global electrochemical responses rely on an efficient interface between the inside and the cell's exterior, which is driven by the highly abundant periplasmic triheme PpcA-family cytochromes. The functional features of these cytochromes have been studied in G. sulfurreducens and G. metallireducens, and although they share a high degree of structural homology and sequence identity, their properties are quite distinct. In this work, the heme axial ligand geometries and the magnetic properties of PpcF from G. metallireducens were determined. The data obtained constitute important constraints for the determination of its solution structure in the oxidized state and indicate that the (i) heme core architecture; (ii) axial ligands geometries and (iii) magnetic properties of the cytochrome are conserved compared to the other members of the PpcA-families. Furthermore, the results also indicate that the heme arrangement is crucial to maintain an intrinsic regulation of the protein's redox properties and hence its electron transfer efficiency and functionality.

Flavocytochrome c3 is a periplasmic fumarate reductase with Mr 63.8 kDa, isolated from Shewanella frigidimarina NCIMB400. NMR spectroscopy was tested for its potential to elucidate the oxidation profile of each of the four haem groups in the enzyme, using the strategy developed previously to perform the thermodynamic characterization of small tetrahaem cytochromes (FEBS Lett. 314 (1992) 155). This work shows that, despite the large size of the protein, 2D-NMR NOESY experiments can be used to obtain the network of chemical exchange connectivities, between the signals of specific haem groups in sequential oxidation stages.

The tetrahaem cytochrome isolated during anaerobic growth of Shewanella frigidimarina NCIMB400 is a small protein (86 residues) involved in electron transfer to Fe(III), which can be used as a terminal respiratory oxidant by this bacterium. A 3D solution structure model of the reduced form of the cytochrome has been determined using NMR data in order to determine the relative orientation of the haems. The haem core architecture of S. frigidimarina tetrahaem cytochrome differs from that found in all small tetrahaem cytochromes c3 so far isolated from strict anaerobes, but has some similarity to the N-terminal cytochrome domain of flavocytochrome c3 isolated from the same bacterium. NMR signals obtained for the four haems of S. frigidimarina tetrahaem cytochrome at all stages of oxidation were cross-assigned to the solution structure using the complete network of chemical exchange connectivities. Thus, the order in which each haem in the structure becomes oxidised was determined.

The thermodynamic properties of the Desulfovibrio vulgaris (Hildenborough) tetrahaem cytochrome c3 (Dvc3) are rationalised by a model which involves both homotropic (e−/e−) and heterotropic (e−/H+) cooperativity. The paramagnetic shifts of a methyl group from each haem of the DVc3 have been determined in each stage of oxidation at several pH values by means of two-dimensional exchange NMR. The thermodynamic parameters are obtained by fitting the model to the NMR data and to redox titrations followed by visible spectroscopy. They show significant positive cooperativity between two of the haems whereas the remaining interactions appear to be largely electrostatic in origin. These parameters imply that the protein undergoes a proton-assisted two-electron transfer which can be used for energy transduction. Comparison with the crystal structure together with measurement of the kinetics of proton exchange suggest that the pH dependence is mediated by a charged residue(s) readily acessible to the solvent and close to haem I.

The periplasmic sensor domains encoded by genes gsu0582 and gsu0935 are part of methyl accepting chemotaxis proteins in the bacterium Geobacter sulfurreducens (Gs). The sensor domains of these proteins contain a heme-c prosthetic group and a PAS-like fold as revealed by their crystal structures. Biophysical studies of the two domains showed that nitric oxide (NO) binds to the heme in both the ferric and ferrous forms, whereas carbon monoxide (CO) binds only to the reduced form. In order to address these exogenous molecules as possible physiological ligands, binding studies and resonance Raman (RR) spectroscopic characterization of the respective CO and NO adducts were performed in this work. In the absence of exogenous ligands, typical RR frequencies of five-coordinated (5c) high-spin and six-coordinated (6c) low-spin species were observed in the oxidized form. In the reduced state, only frequencies corresponding to the latter were detected. In both sensors, CO binding yields 6c low-spin adducts by replacing the endogenous distal ligand. The binding of NO by the two proteins causes partial disruption of the proximal Fe-His bond, as revealed by the RR fingerprint features of 5cFe-NO and 6cNO-Fe-His species. The measured CO and NO dissociation constants of ferrous GSU0582 and GSU0935 sensors reveal that both proteins have high and similar affinity toward these molecules (Kd ≈ 0.04−0.08 μM). On the contrary, in the ferric form, sensor GSU0582 showed a much higher affinity for NO (Kd ≈ 0.3 μM for GSU0582 versus 17 μM for GSU0935). Molecular dynamics calculations revealed a more open heme pocket in GSU0935, which could account for the different affinities for NO. Taken together, spectroscopic data and MD calculations revealed subtle differences in the binding properties and structural features of formed CO and NO adducts, but also indicated a possibility that a (5c) high-spin/(6c) low-spin redox-linked equilibrium could drive the physiological sensing of Gs cells.

Electrogenic microorganisms possess unique redox biological features, being capable of transferring electrons to the cell exterior and converting highly toxic compounds into nonhazardous forms. These microorganisms have led to the development of Microbial Electrochemical Technologies (METs), which include applications in the fields of bioremediation and bioenergy production. The optimization of these technologies involves efforts from several different disciplines, ranging from microbiology to materials science. Geobacter bacteria have served as a model for understanding the mechanisms underlying the phenomenon of extracellular electron transfer, which is highly dependent on a multitude of multiheme cytochromes (MCs). MCs are, therefore, logical targets for rational protein engineering to improve the extracellular electron transfer rates of these bacteria. However, the presence of several heme groups complicates the detailed redox characterization of MCs. In this Review, the main characteristics of electroactive Geobacter bacteria, their potential to develop microbial electrochemical technologies and the main features of MCs are initially highlighted. This is followed by a detailed description of the current methodologies that assist the characterization of the functional redox networks in MCs. Finally, it is discussed how this information can be explored to design optimal Geobacter-mutated strains with improved capabilities in METs.

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.

Accepted Manuscript online November 14, 2016.The periplasmic triheme cytochrome PpcA from Geobacter sulfurreducens is highly abundant; it is the likely reservoir of electrons to the outer surface to assist the reduction of extracellular terminal acceptors; these include insoluble metal oxides in natural habitats and electrode surfaces from which electricity can be harvested. A detailed thermodynamic characterization of PpcA showed that it has an important redox-Bohr effect that might implicate the protein in e-/H+ coupling mechanisms to sustain cellular growth. This functional mechanism requires control of both the redox state and the protonation state. In the present study, isotope-labeled PpcA was produced and the three-dimensional structure of PpcA in the oxidized form was determined by NMR. This is the first solution structure of a G. sulfurreducens cytochrome in the oxidized state. The comparison of oxidized and reduced structures revealed that the heme I axial ligand geometry changed and there were other significant changes in the segments near heme I. The pH-linked conformational rearrangements observed in the vicinity of the redox-Bohr center, both in the oxidized and reduced structures, constitute the structural basis for the differences observed in the pKa values of the redox-Bohr center, providing insights into the e-/H+ coupling molecular mechanisms driven by PpcA in G. sulfurreducens.EET, extracellular electron transfer; IM, inner membrane; IPTG, isopropyl β-d-thiogalactoside; MFCs, microbial fuel cells; NOE, Nuclear Overhauser effect; OM, outer membrane; rmsd, root mean square deviation.

The thermodynamic model of five interacting charge centres (four haems and an ionisable centre), which was used in the characterisation of the thermodynamic properties of Desulfovibrio vulgaris (Hildenborough) cytochrome c3 (c3DvH), is now used to reevaluate the thermodynamic properties in Desulfovibrio vulgaris (Miyazaki F) cytochrome c3 (c3DvM) on the basis of published data (Park, J.-S., Ohmura, T., Kano, K., Sagara, T., Niki, K., Kyogoku, Y. and Akutsu, H. (1996) Biochim. Biophys. Acta 1293, 45–54). Contrary to the assertion of Park et al. (1996), the pH dependence of the proton chemical shifts of haem methyls in c3DvM in several stages of oxidation is well described by the model, which involves both homotropic (e–/e–) and heterotropic (e–/H+) cooperativity. This shows that the pH dependence observed for c3DvM is not significantly more complicated than that observed for c3DvH. Since the parameters which we now obtain for c3DvM are generated with the same model as those from c3DvH, albeit using less precise data, it is possible to make a preliminary comparison of the thermodynamic properties of these two proteins and of their role in energy transduction.
The extrinsic dipolar shifts generated for each methyl group by each of the four haems in c3DvM are also determined. A novel method for approximating the magnetic susceptibility tensors is used: the orientations of the principal axes of the tensors have been shown to be closely related to the geometry of the axial ligands, which is available from the X-ray structure of c3DvM, and the components of the tensors are extrapolated from EPR g values. The inclusion of the calculated haem extrinsic contributions clearly describes the pH dependence of the haem methyls in the core of the protein, close to other haems. This description is most remarkable in the case of the haem methyl 21CH3 II I, for which the "unusual pH dependence" commented on by Park et al. (1996) is easily explained using the thermodynamic parameters determined by our model together with the calculated extrinsic dipolar shifts, thus providing a test of the analysis.