Orientation of the g-Tensor Axes of the Rieske Subunit in the Cytochrome bc1 Complex

Catalytic Reactions and Energy Conservation in the Cytochrome bc1 and b6f Complexes of Energy-Transducing Membranes

This review focuses on key components of respiratory and photosynthetic energy-transduction systems: the cytochrome bc1 and b6f (Cytbc1/b6f) membranous multisubunit homodimeric complexes. These remarkable molecular machines catalyze electron transfer from membranous quinones to water-soluble electron carriers (such as cytochromes c or plastocyanin), coupling electron flow to proton translocation across the energy-transducing membrane and contributing to the generation of a transmembrane electrochemical potential gradient, which powers cellular metabolism in the majority of living organisms. Cytsbc1/b6f share many similarities but also have significant differences. While decades of research have provided extensive knowledge on these enzymes, several important aspects of their molecular mechanisms remain to be elucidated. We summarize a broad range of structural, mechanistic, and physiological aspects required for function of Cytbc1/b6f, combining textbook fundamentals with new intriguing concepts that have emerged from more recent studies. The discussion covers but is not limited to (i) mechanisms of energy-conserving bifurcation of electron pathway and energy-wasting superoxide generation at the quinol oxidation site, (ii) the mechanism by which semiquinone is stabilized at the quinone reduction site, (iii) interactions with substrates and specific inhibitors, (iv) intermonomer electron transfer and the role of a dimeric complex, and (v) higher levels of organization and regulation that involve Cytsbc1/b6f. In addressing these topics, we point out existing uncertainties and controversies, which, as suggested, will drive further research in this field.

EPR-Derived Structure of a Paramagnetic Intermediate Generated by Biotin Synthase BioB

Biotin (vitamin B7) is an enzyme cofactor required by organisms from all branches of life but synthesized only in microbes and plants. In the final step of biotin biosynthesis, a radical S-adenosyl-l-methionine (SAM) enzyme, biotin synthase (BioB), converts the substrate dethiobiotin to biotin through the stepwise formation of two C-S bonds. Previous electron paramagnetic resonance (EPR) spectroscopic studies identified a semistable intermediate in the formation of the first C-S bond as 9-mercaptodethiobiotin linked to a paramagnetic [2Fe-2S] cluster through one of its bridging sulfides. Herein, we report orientation-selected pulse EPR spectroscopic results that reveal hyperfine interactions between the [2Fe-2S] cluster and a number of magnetic nuclei (e.g., 57Fe, 15N, 13C, and 2H) introduced in a site-specific manner via biosynthetic methods. Combining these results with quantum chemical modeling gives a structural model of the intermediate showing that C6, the target of the second hydrogen-atom abstraction, is now in close proximity to the nascent thioether sulfur and is ideally positioned for the second C-S bond forming event.

g-Tensor Directions in the Protein Structural Frame of Hyperthermophilic Archaeal Reduced Rieske-Type Ferredoxin Explored by 13C Pulsed Electron Paramagnetic Resonance

Interpretation of magnetic resonance data in the context of structural and chemical biology requires prior knowledge of the g-tensor directions for paramagnetic metallo-cofactors with respect to the protein structural frame. Access to this information is often limited by the strict requirement of suitable protein crystals for single-crystal electron paramagnetic resonance (EPR) measurements or the reliance on protons (with ambiguous locations in crystal structures) near the paramagnetic metal site. Here we develop a novel pulsed EPR approach with selective ¹³C_β-cysteine labeling of model [2Fe-2S] proteins to help bypass these problems. Analysis of the ¹³C_β-cysteine hyperfine tensors reproduces the g-tensor of the Pseudomonas putida ISC-like [2Fe-2S] ferredoxin (FdxB). Its application to the hyperthermophilic archaeal Rieske-type [2Fe-2S] ferredoxin (ARF) from Sulfolobus solfataricus, for which the single-crystal EPR approach was not feasible, supports the best-fit g _x-, g _z-, and g _y-tensor directions of the reduced cluster as nearly along Fe-Fe, S-S, and the cluster plane normal, respectively. These approximate principal directions of the reduced ARF g-tensor, explored by ¹³C pulsed EPR, are less skewed from the cluster molecular axes and are largely consistent with those previously determined by single-crystal EPR for the cytochrome bc₁-associated, reduced Rieske [2Fe-2S] center. This suggests the approximate g-tensor directions are conserved across the phylogenetically and functionally divergent Rieske-type [2Fe-2S] proteins.

Advanced paramagnetic resonance spectroscopies of iron-sulfur proteins: Electron nuclear double resonance (ENDOR) and electron spin echo envelope modulation (ESEEM)

The advanced electron paramagnetic resonance (EPR) techniques, electron nuclear double resonance (ENDOR) and electron spin echo envelope modulation (ESEEM) spectroscopies, provide unique insights into the structure, coordination chemistry, and biochemical mechanism of nature's widely distributed iron-sulfur cluster (FeS) proteins. This review describes the ENDOR and ESEEM techniques and then provides a series of case studies on their application to a wide variety of FeS proteins including ferredoxins, nitrogenase, and radical SAM enzymes. This article is part of a Special Issue entitled: Fe/S proteins: Analysis, structure, function, biogenesis and diseases.

Fast proton-coupled electron transfer observed for a high-fidelity structural and functional [2Fe-2S] Rieske model

Rieske cofactors have a [2Fe–2S] cluster with unique {His₂Cys₂} ligation and distinct Fe subsites. The histidine ligands are functionally relevant, since they allow for coupling of electron and proton transfer (PCET) during quinol oxidation in respiratory and photosynthetic ET chains. Here we present the highest fidelity synthetic analogue for the Rieske [2Fe–2S] cluster reported so far. This synthetic analogue 5^x– emulates the heteroleptic {His₂Cys₂} ligation of the [2Fe–2S] core, and it also serves as a functional model that undergoes fast concerted proton and electron transfer (CPET) upon reaction of the mixed-valent (ferrous/ferric) protonated 5H^2– with TEMPO. The thermodynamics of the PCET square scheme for 5^x– have been determined, and three species (diferric 5^2–, protonated diferric 5H^–, and mixed-valent 5^3–) have been characterized by X-ray diffraction. pK_a values for 5H^– and 5H^2– differ by about 4 units, and the reduction potential of 5H^– is shifted anodically by about +230 mV compared to that of 5^2–. While the N–H bond dissociation free energy of 5H^2– (60.2 ± 0.5 kcal mol^–1) and the free energy, ΔG°_CPET, of its reaction with TEMPO (−6.3 kcal mol^–1) are similar to values recently reported for a homoleptic {N₂/N₂}-coordinated [2Fe–2S] cluster, CPET is significantly faster for 5H^2– with biomimetic {N₂/S₂} ligation (k = (9.5 ± 1.2) × 10⁴ M^–1 s^–1, ΔH^‡ = 8.7 ± 1.0 kJ mol^–1, ΔS^‡ = −120 ± 40 J mol^–1 K^–1, and ΔG^‡ = 43.8 ± 0.3 kJ mol^–1 at 293 K). These parameters, and the comparison with homoleptic analogues, provide important information and new perspectives for the mechanistic understanding of the biological Rieske cofactor.

Protein conformational changes involved in the cytochrome bc1 complex catalytic cycle

Early structures of the cytochrome bc1 complex revealed heterogeneity in the position of the soluble portion of the Rieske iron sulfur protein subunit, implicating a movement of this domain during function. Subsequent biochemical and biophysical works have firmly established that the motion of this subunit acts in the capacity of a conformationally assisted electron transfer step during the already complicated catalytic mechanism described within the modified version of Peter Mitchells Q cycle. How the movement of this subunit is initiated or how the frequency of its motion is controlled as a function of other steps during the catalysis remain topics of debate within the active research communities. This review addresses the historical aspects of the discovery and description of this movement, while attempting to provide a context for the involvement of conformational motion in the catalysis and efficiency of the enzyme. This article is part of a Special Issue entitled: Respiratory complex III and related bc complexes.

ISC-like [2Fe-2S] ferredoxin (FdxB) dimer from Pseudomonas putida JCM 20004: structural and electron-nuclear double resonance characterization

The crystal structure of the ISC-like [2Fe-2S] ferredoxin (FdxB), probably involved in the de novo iron-sulfur cluster biosynthesis (ISC) system of Pseudomonas putida JCM 20004, was determined at 1.90-Å resolution and displayed a novel tail-to-tail dimeric form. P. putida FdxB lacks the consensus free cysteine usually present near the cluster of ISC-like ferredoxins, indicating its primarily electron transfer role in the iron-sulfur cluster. Orientation-selective electron-nuclear double resonance spectroscopic analysis of reduced FdxB in conjunction with the crystal structure has identified the innermost Fe2 site with a high positive spin population as the nonreducible iron retaining the Fe(3+) valence and the outermost Fe1 site as the reduced iron with a low negative spin density. The average g (max) direction is skewed, forming an angle of about 27.3° (±4°) with the normal of the [2Fe-2S] plane, whereas the g (int) and g (min) directions are distributed in the cluster plane, presumably tilted by the same angle with respect to this plane. These results are related to those for other [2Fe-2S] proteins in different electron transport chains (e.g. adrenodoxin) and suggest a significant distortion of the electronic structure of the reduced [2Fe-2S] cluster under the influence of the protein environment around each iron site in general.

Binding of histidine in the (Cys)3(His)1-coordinated [2Fe-2S] cluster of human mitoNEET

Human mitoNEET is a homodimeric iron-sulfur protein located in the outer mitochondrial membrane with unknown function, but which is known to interact with thiazolidinedione diabetes drugs. Each monomer houses a [2Fe-2S] cluster with an unusual (Cys)(3)(His)(1) ligation. The His ligand is important for enabling cluster release and for tuning the redox potential. We use multifrequency (X-, Ka-, and Q-band) and multitechnique (continuous-wave, electron spin-echo envelope modulation (ESEEM), pulsed electron-nuclear double resonance (ENDOR), and hyperfine sublevel correlation (HYSCORE)) electron paramagnetic resonance spectroscopy to investigate the cluster in its paramagnetic reduced [Fe(2+)Fe(3+)] (S = 1/2) state. It has a rhombic g tensor (2.007, 1.937, 1.897) with an average g value of 1.947 that falls between those of Rieske-type and ferredoxin-type [2Fe-2S] clusters. Simulation and least-squares fitting of orientation-selective Ka- and Q-band ENDOR, 1D ESEEM, and HYSCORE spectra of (14)N and (15)N-labeled mitoNEET yield the principal values and orientations of both the hyperfine tensor ((14)N, A(iso) = -6.25 MHz, T = -0.94 MHz) and the quadrupolar tensor (e(2)Qq/h = -2.47 MHz, eta = 0.38) of the ligating histidine nitrogen N(delta). From these, we can infer the absolute g tensor orientation with respect to the cluster: The g(2) axis is close to perpendicular to the [2Fe-2S] plane, and g(1) and g(3) are in-plane, but skewed from the Fe-Fe and S-S axes. In X-band ENDOR and ESEEM spectra, a weakly coupled nitrogen is visible, most likely the N(epsilon) of the histidine in the protonated state. We find that the cluster is in a valence-localized state, where Fe(2+) is His-bound. The field-sweep spectra show evidence of intercluster dipolar coupling that can be simulated using an uncoupled spin model for each cluster (S(Fe(2+)) = 2, S(Fe(3+)) = 5/2). The parameters determined in this work can function as reporters on how the cluster structure is altered upon pH changes and drug binding.

Continuous-wave and pulsed EPR characterization of the [2Fe-2S](Cys)3(His)1 cluster in rat MitoNEET

CW EPR spectra of reduced [2Fe-2S](Cys)(3)(His)(1) clusters of mammalian mitoNEET soluble domain appear to produce features resulting from the interaction of the electron spins of the two adjacent clusters, which can be explained by employing the local spin model. This model favors the reduction of the outermost iron with His87 and Cys83 ligands, which is supported by orientation-selected hyperfine sublevel correlation (HYSCORE) characterization of the uniformly (15)N-labeled mitoNEET showing one strongly coupled nitrogen from the His87 N(delta) ligand with hyperfine coupling (15)a = 8 MHz. The (14)N and (15)N HYSCORE spectra also exhibit at least two different cross-peaks located near diagonal in the (++) quadrant, with frequencies approximately 2.8 and 2.4 MHz (N2), and the other approximately 4.0 and 3.5 MHz (N1), but did not show any of the larger splitting approximately 1.1-1.4 MHz previously seen with Rieske proteins. Further analysis with partially (15)N(3)-His-labeled protein indicates that His87 N(epsilon) cross-peaks produce resolved features (N2) in the (14)N spectrum but contribute much less than weakly coupled peptide nitrogen species to the (++) quadrant in the (15)N spectrum. It is suggested that these quantitative data may be used in future functional and theoretical studies on the mammalian mitoNEET [2Fe-2S] cluster system.

The reduced [2Fe-2S] clusters in adrenodoxin and Arthrospira platensis ferredoxin share spin density with protein nitrogens, probed using 2D ESEEM

We have used X-band ESEEM to study the reduced [2Fe-2S] cluster in adrenodoxin and Arthrospira platensis ferredoxin. By use of a 2D approach (HYSCORE), we have shown that the cluster is involved in weak magnetic interactions with several nitrogens in each protein. Despite substantial differences in the shape and orientational dependence of individual cross-peaks, the major spectral features in both proteins are attributable to two peptide nitrogens (N1 and N2) with similar hyperfine couplings approximately 1.1 and approximately 0.70 MHz. The couplings determined represent a small fraction (0.0003-0.0005) of the unpaired spin density of the reduced cluster transferred to these nitrogens over H-bond bridges or the covalent bonds of cysteine ligands. Simulation of the HYSCORE spectra has allowed us to estimate the orientation of the nuclear quadrupole tensors of N1 and N2 in the g-tensor coordinate system. The most likely candidates for the role of N1 and N2 have been identified in the protein environment by comparing magnetic-resonance data with crystallographic structures of the oxidized proteins. A possible influence of redox-linked structural changes on ESEEM data is analyzed using available structures for related proteins in two redox states.

Magnetic interactions sense changes in distance between heme b(L) and the iron-sulfur cluster in cytochrome bc(1)

During the operation of cytochrome bc₁, a key enzyme of biological energy conversion, the iron−sulfur head domain of one of the subunits of the catalytic core undergoes a large-scale movement from the catalytic quinone oxidation Q_o site to cytochrome c₁. This changes a distance between the two iron−two sulfur (FeS) cluster and other cofactors of the redox chains. Although the role and the mechanism of this movement have been intensely studied, they both remain poorly understood, partly because the movement itself is not easily traceable experimentally. Here, we take advantage of magnetic interactions between the reduced FeS cluster and oxidized heme b_L to use dipolar enhancement of phase relaxation of the FeS cluster as a spectroscopic parameter which with a unique clarity and specificity senses changes in the distance between those two cofactors. The dipolar relaxation curves measured by EPR at Q-band in a glass state of frozen solution (i.e., under the conditions trapping a dynamic distribution of FeS positions that existed in a liquid phase) of isolated cytochrome bc₁ were compared with the curves calculated for the FeS cluster occupying distinct positions in various crystals of cytochrome bc₁. This comparison revealed the existence of a broad distribution of the FeS positions in noninhibited cytochrome bc₁ and demonstrated that the average equilibrium position is modifiable by inhibitors or mutations. To explain the results, we assume that changes in the equilibrium distribution of the FeS positions are the result of modifications of the orienting potential gradient in which the diffusion of the FeS head domain takes place. The measured changes in the phase relaxation enhancement provide the first direct experimental description of changes in the strength of dipolar coupling between the FeS cluster and heme b_L.

Across membrane communication between the Q(o) and Q(i) active sites of cytochrome bc(1)

The ubihydroquinone:cytochrome c oxidoreductase (cyt bc(1)) contains two catalytically active domains, termed the hydroquinone oxidation (Q(o)) and quinone reduction (Q(i)) sites, which are distant from each other by over 30 A. Previously, we have reported that binding of inhibitors to the Q(i) site on one (n) side of the energy-transducing membrane changes the local environment of the iron-sulfur (Fe/S) protein subunit residing in the Q(o) site on the other (p) side of the lipid bilayer [Cooley, J. W., Ohnishi, T., and Daldal, F. (2005) Biochemistry 44, 10520-10532]. These findings best fit a model whereby the Q(o) and Q(i) sites of the cyt bc(1) are actively coupled in spite of their distant locations. Because the Fe/S protein of the cyt bc(1) undergoes a large-scale (macro) domain movement during catalysis, we examined various macromobility-defective Fe/S subunit mutants to assess the role of this motion on the coupling of the active sites and also during the multiple turnovers of the enzyme. By monitoring the changing environments of the Fe/S protein [2Fe-2S] cluster upon addition of Q(i) site inhibitors in selected mutants, we found that the Q(o)-Q(i) site interactions manifest differently depending on the ability of the Fe/S protein to move between the cytochrome b and cytochrome c(1) subunits of the enzyme. In the presence of antimycin A, an immobile Fe/S protein mutant exhibited no changes in its EPR spectra. In contrast, mobility-restricted mutants showed striking alterations in the EPR line shapes and revealed two discrete subpopulations in respect to the [2Fe-2S] cluster environments at the Q(o) site. These findings led us to conclude that the mobility of the Fe/S protein is involved in its response to the occupancy of the Q(i) site by different molecules. We propose that the heterogeneity seen might reflect the distinct responses of the two Fe/S proteins at the Q(o) sites of the dimeric enzyme upon the occupancy of the Q(i) sites and discuss it in terms of the function of the dimeric cyt bc(1) during its multiple turnovers.

Proton environment of reduced Rieske iron-sulfur cluster probed by two-dimensional ESEEM spectroscopy

The proton environment of the reduced [2Fe-2S] cluster in the water-soluble head domain of the Rieske iron-sulfur protein (ISF) from the cytochrome bc(1) complex of Rhodobacter sphaeroides has been studied by orientation-selected X-band 2D ESEEM. The 2D spectra show multiple cross-peaks from protons, with considerable overlap. Samples in which (1)H(2)O water was replaced by (2)H(2)O were used to determine which of the observed peaks belong to exchangeable protons, likely involved in hydrogen bonds in the neighborhood of the cluster. By correlating the cross-peaks from 2D spectra recorded at different parts of the EPR spectrum, lines from nine distinct proton signals were identified. Assignment of the proton signals was based on a point-dipole model for interaction with electrons of Fe(III) and Fe(II) ions, using the high-resolution structure of ISF from Rb. sphaeroides. Analysis of experimental and calculated tensors has led us to conclude that even 2D spectra do not completely resolve all contributions from nearby protons. Particularly, the seven resolved signals from nonexchangeable protons could be produced by at least 13 protons. The contributions from exchangeable protons were resolved by difference spectra ((1)H(2)O minus (2)H(2)O), and assigned to two groups of protons with distinct anisotropic hyperfine values. The largest measured coupling exceeded any calculated value. This discrepancy could result from limitations of the point dipole approximation in dealing with the distribution of spin density over the sulfur atoms of the cluster and the cysteine ligands, or from differences between the structure in solution and the crystallographic structure. The approach demonstrated here provides a paradigm for a wide range of studies in which hydrogen-bonding interactions with metallic centers has a crucial role in understanding the function.

Asymmetric Binding of Stigmatellin to the Dimeric Paracoccus denitrificans bc1 Complex

We have investigated the mechanism responsible for half-of-the-sites activity in the dimeric cytochrome bc(1) complex from Paracoccus denitrificans by characterizing the kinetics of inhibitor binding to the ubiquinol oxidation site at center P. Both myxothiazol and stigmatellin induced a 2-3 nm shift of the visible absorbance spectrum of the b(L) heme. The shift generated by myxothiazol was symmetric, with monophasic kinetics that indicate equal binding of this inhibitor to both center P sites. In contrast, stigmatellin generated an asymmetric shift in the b(L) spectrum, with biphasic kinetics in which each phase contributed approximately half of the total magnitude of the spectral change. The faster binding phase corresponded to a more symmetrical shift of the b(L) spectrum relative to the slower binding phase, indicating that approximately half of the center P sites bound stigmatellin more slowly and in a different position relative to the b(L) heme, generating a different effect on its electronic environment. Significantly, the slow stigmatellin binding phase was lost as the inhibitor concentration was increased. This implies that a conformational change is transmitted from one center P site in the dimer to the other upon stigmatellin binding to one monomer, rendering the second site less accessible to the inhibitor. Because the position that stigmatellin occupies at center P is considered to be analogous to that of the quinol substrate at the moment of electron transfer, these results indicate that the productive enzyme-substrate configuration is prevented from occurring in both monomers simultaneously.

Plasmon waveguide resonance spectroscopic evidence for differential binding of oxidized and reduced Rhodobacter capsulatus cytochrome c2 to the cytochrome bc1 complex mediated by the conformation of the Rieske iron-sulfur protein

The dissociation constants for the binding of Rhodobacter capsulatus cytochrome c2 and its K93P mutant to the cytochrome bc1 complex embedded in a phospholipid bilayer were measured by plasmon waveguide resonance spectroscopy in the presence and absence of the inhibitor stigmatellin. The reduced form of cytochrome c2 strongly binds to reduced cytochrome bc1 (Kd = 0.02 microM) but binds much more weakly to the oxidized form (Kd = 3.1 microM). In contrast, oxidized cytochrome c2 binds to oxidized cytochrome bc1 in a biphasic fashion with Kd values of 0.11 and 0.58 microM. Such a biphasic interaction is consistent with binding to two separate sites or conformations of oxidized cytochrome c2 and/or cytochrome bc1. However, in the presence of stigmatellin, we find that oxidized cytochrome c2 binds to oxidized cytochrome bc1 in a monophasic fashion with high affinity (Kd = 0.06 microM) and reduced cytochrome c2 binds less strongly (Kd = 0.11 microM) but approximately 30-fold more tightly than in the absence of stigmatellin. Structural studies with cytochrome bc1, with and without the inhibitor stigmatellin, have led to the proposal that the Rieske protein is mobile, moving between the cytochrome b and cytochrome c1 components during turnover. In one conformation, the Rieske protein binds near the heme of cytochrome c1, while the cytochrome c2 binding site is also near the cytochrome c1 heme but on the opposite side from the Rieske site, where cytochrome c2 cannot directly interact with Rieske. However, the inhibitor, stigmatellin, freezes the Rieske protein iron-sulfur cluster in a conformation proximal to cytochrome b and distal to cytochrome c1. We conclude from this that the dual conformation of the Rieske protein is primarily responsible for biphasic binding of oxidized cytochrome c2 to cytochrome c1. This optimizes turnover by maximizing binding of the substrate, oxidized cytochrome c2, when the iron-sulfur cluster is proximal to cytochrome b and minimizing binding of the product, reduced cytochrome c2, when it is proximal to cytochrome c1.

Identification of Hydrogen Bonds to the Rieske Cluster through the Weakly Coupled Nitrogens Detected by Electron Spin Echo Envelope Modulation Spectroscopy

The interaction of the reduced[2Fe-2S] cluster of isolated Rieske fragment from the bc1 complex of Rhodobacter sphaeroides with nitrogens (14N and 15N) from the local protein environment has been studied by X- and S-band pulsed EPR spectroscopy. The two-dimensional electron spin echo envelope modulation spectra of uniformly 15N-labeled protein show two well resolved cross-peaks with weak couplings of approximately 0.3-0.4 and 1.1 MHz in addition to couplings in the range of 6-8 MHz from two coordinating Ndelta of histidine ligands. The quadrupole coupling constants for weakly coupled nitrogens determined from S-band electron spin echo envelope modulation spectra identify them as Nepsilon of histidine ligands and peptide nitrogen (Np), respectively. Analysis of the line intensities in orientation-selected S-band spectra indicated that Np is the backbone N-atom of Leu-132 residue. The hyperfine couplings from Nepsilon and Np demonstrate the predominantly isotropic character resulting from the transfer of unpaired spin density onto the 2s orbitals of the nitrogens. Spectra also show that other peptide nitrogens in the protein environment must carry a 5-10 times smaller amount of spin density than the Np of Leu-132 residue. The appearance of the excess unpaired spin density on the Np of Leu-132 residue indicates its involvement in hydrogen bond formation with the bridging sulfur of the Rieske cluster. The configuration of the hydrogen bond therefore provides a preferred path for spin density transfer. Observation of similar splittings in the 15N spectra of other Rieske-type proteins and [2Fe-2S] ferredoxins suggests that a hydrogen bond between the bridging sulfur and peptide nitrogen is a common structural feature of [2Fe-2S] clusters.

A comparative, two-dimensional 14N ESEEM characterization of reduced [2Fe-2S] clusters in hyperthermophilic archaeal high- and low-potential Rieske-type proteins

Proteins of the Rieske and Rieske-type family contain a [2Fe-2S] cluster with mixed ligation by two histidines and two cysteines, and play important roles in various biological electron transfer reactions. We report here the comparative orientation-selected ESEEM and HYSCORE studies of the reduced clusters from two hyperthermophilic Rieske-type proteins; a high-potential, archaeal Rieske protein called sulredoxin (SDX) from Sulfolobus tokodaii with weak homology to the cytochrome bc-associated Rieske proteins, and a low-potential, archaeal homolog of an oxygenase-associated Rieske-type ferredoxin (ARF) from Sulfolobus solfataricus. (14)N ESEEM and HYSCORE spectra of SDX and ARF show well-defined variations, which are primarily determined by changes of quadrupole couplings (up to 50% depending on the selected orientation) of the two coordinated nitrogens. These are due to variations in coordination geometry of the histidine imidazole ligands rather than to variations of hyperfine couplings of these nitrogens, which do not exceed 8-10%. The measured quadrupole couplings and their differences in the two proteins are consistent with those calculated using the reported crystal structures of high- and low-potential Rieske proteins. These results suggest that exploration of quadrupole tensors might provide a more accurate method for characterization of the histidine coordination in different proteins and mutants than hyperfine tensors, and might have potential applications in a wider range of biological systems.