Search

Secondary coordination sphere (SCS) interactions have been shown to play important roles in tuning reduction potentials and electron transfer (ET) properties of the Type 1 copper proteins, but the precise roles of these interactions are not fully understood. In this work, we examined the influence of F114P, F114N, and N47S mutations in the SCS on the electronic structure of the T1 copper center in azurin (Az) by studying the hyperfine couplings of (i) histidine remote N ε nitrogens and (ii) the amide N p using the two-dimensional (2D) pulsed electron paramagnetic resonance (EPR) technique HYSCORE (hyperfine sublevel correlation) combined with quantum mechanics/molecular mechanics (QM/MM) and DLPNO-CCSD calculations. Our data show that some components of hyperfine tensor and isotropic coupling in N47SAz and F114PAz (but not F114NAz) deviate by up to ∼±20% from WTAz, indicating that these mutations significantly influence the spin density distribution between the Cu II site and coordinating ligands. Furthermore, our calculations support the assignment of N p to the backbone amide of residue 47 (both in Asn and Ser variants). Since the spin density distributions play an important role in tuning the covalency of the Cu-Scys bond of Type 1 copper center that has been shown to be crucial in controlling the reduction potentials, this study provides additional insights into the electron spin factor in tuning the reduction potentials and ET properties.

The "Histidine-brace" (His-brace) copper-binding site, composed of Cu(His) 2 with a backbone amine, is found in metalloproteins with diverse functions. A primary example is lytic polysaccharide monooxygenase (LPMO), a class of enzymes that catalyze the oxidative depolymerization of polysaccharides, providing not only an energy source for native microorganisms but also a route to more effective industrial biomass conversion. Despite its importance, how the Cu His-brace site performs this unique and challenging oxidative depolymerization reaction remains to be understood. To answer this question, we have designed a biosynthetic model of LPMO by incorporating the Cu His-brace motif into azurin, an electron transfer protein. Spectroscopic studies, including ultraviolet-visible (UV-Vis) absorption and electron paramagnetic resonance, confirm copper binding at the designed His-brace site. Moreover, the designed protein is catalytically active towards both cellulose and starch, the native substrates of LPMO, generating degraded oligosaccharides with multiturnovers by C1 oxidation. It also performs oxidative cleavage of the model substrate 4-nitrophenyl-D-glucopyranoside, achieving a turnover number ~9% of that of a native LPMO assayed under identical conditions. This work presents a rationally designed artificial metalloenzyme that acts as a structural and functional mimic of LPMO, which provides a promising system for understanding the role of the Cu His-brace site in LPMO activity and potential application in polysaccharide degradation.

Reaction centers from Rhodobacter sphaeroides with residue M265 mutated from isoleucine to threonine, serine, and asparagine (M265IT, M265IS, and M265IN, respectively) in the Q A -· state are studied by high-resolution electron spin echo envelope modulation (ESEEM) and electron nuclear double resonance spectroscopy methods to investigate the structural characteristics of these mutants influencing the redox properties of the Q A site. All three mutants decrease the redox midpoint potential ( E m ) of Q A by ∼0.1 V, yet the mechanism for this drop in E m is unclear. In this work, we examine (i) the hydrogen bonding interactions between Q A - · and residues histidine M219 and alanine M260, (ii) the electron spin density distribution of the semiquinone, and (iii) the orientations of the ubiquinone methoxy substituents. 13 C measurements show no significant contribution of methoxy dihedral angles to the observed decrease in E m for the Q A mutants. Instead, 14 N three-pulse ESEEM data suggest that electrostatic or hydrogen bond formation between the mutated M265 side chain and His-M219 N δ may be involved in the observed lowering of the Q A midpoint potential. For mutant M265IN, analysis of the proton hyperfine couplings reveals a weakened hydrogen bond network, resulting in an altered Q A - · spin density distribution. The magnetic resonance study presented here is most consistent with an electrostatic or structural perturbation of the His-M219 N δ hydrogen bond in these mutants as a mechanism for the ∼0.1 V decrease in Q A E m .

A set of C43(DE3) and BL21(DE3) Escherichia coli host strains that are auxotrophic for various amino acids is briefly reviewed. These strains require the addition of a defined set of one or more amino acids in the growth medium, and have been specifically designed for overproduction of membrane or water-soluble proteins selectively labelled with stable isotopes, such as 2H, 13C and 15N. The strains described here are available for use and have been deposited into public strain banks. Although they cannot fully eliminate the possibility of isotope dilution and mixing, metabolic scrambling of the different amino acid types can be minimized through a careful consideration of the bacterial metabolic pathways. The use of a suitable auxotrophic expression host strain with an appropriately isotopically labelled growth medium ensures high levels of isotope labelling efficiency as well as selectivity for providing deeper insight into protein structure-function relationships., (© The Author(s) 2020. Published by Oxford University Press on behalf of the Japanese Biochemical Society. All rights reserved.)

The primary and secondary coordination spheres of metal binding sites in metalloproteins have been investigated extensively, leading to the creation of high-performing functional metalloproteins; however, the impact of the overall structure of the protein scaffold on the unique properties of metalloproteins has rarely been studied. A primary example is the binuclear Cu A center, an electron transfer cupredoxin domain of photosynthetic and respiratory complexes and, recently, a protein coregulated with particulate methane and ammonia monooxygenases. The redox potential, Cu-Cu spectroscopic features, and a valence delocalized state of Cu A are difficult to reproduce in synthetic models, and every artificial protein Cu A center to-date has used a modified cupredoxin. Here, we present a fully functional Cu A center designed in a structurally nonhomologous protein, cytochrome c peroxidase (C c P), by only two mutations (Cu A C c P). We demonstrate with UV-visible absorption, resonance Raman, and magnetic circular dichroism spectroscopy that Cu A C c P is valence delocalized. Continuous wave and pulsed (HYSCORE) X-band EPR show it has a highly compact g z area and small A z hyperfine principal value with g and A tensors that resemble axially perturbed Cu A . Stopped-flow kinetics found that Cu A formation proceeds through a single T2Cu intermediate. The reduction potential of Cu A C c P is comparable to native Cu A and can transfer electrons to a physiological redox partner. We built a structural model of the designed Cu binding site from extended X-ray absorption fine structure spectroscopy and validated it by mutation of coordinating Cys and His residues, revealing that a triad of residues (R48C, W51C, and His52) rigidly arranged on one α-helix is responsible for chelating the first Cu(II) and that His175 stabilizes the binuclear complex by rearrangement of the C c P heme-coordinating helix. This design is a demonstration that a highly conserved protein fold is not uniquely necessary to induce certain characteristic physical and chemical properties in a metal redox center.

Cytochrome bo 3 , one of three terminal oxygen reductases in the aerobic respiratory chain of Escherichia coli , has been well characterized as a ubiquinol oxidase. The ability of cytochrome bo 3 to catalyze the two-electron oxidation of ubiquinol-8 requires the enzyme to stabilize the one-electron oxidized ubisemiquinone species that is a transient intermediate in the reaction. Cytochrome bo 3 has been shown recently to also utilize demethylmenaquinol-8 as a substrate that, along with menaquinol-8, replaces ubiquinol-8 when E. coli is grown under microaerobic or anaerobic conditions. In this work, we show that its steady-state turnover with 2,3-dimethyl-1,4-naphthoquinol, a water-soluble menaquinol analogue, is just as efficient as with ubiquinol-1. Using pulsed electron paramagnetic resonance spectroscopy, we demonstrate that the same residues in cytochrome bo 3 that stabilize the semiquinone state of ubiquinone also stabilize the semiquinone state of menaquinone, with the hydrogen bond strengths and the distribution of unpaired spin density accommodated for the different substrate. Catalytic function with menaquinol is more tolerant of mutations at the active site than with ubiquinol. A mutation of one of the stabilizing residues (R71H in subunit I) that eliminates the ubiquinol oxidase activity of cytochrome bo 3 does not abolish activity with soluble menaquinol analogues.

A comparative study of the 1 H and 14 N hyperfine interactions between the Cu A site in an engineered Cu A center in azurin (WT-Cu A Az) and its His120Gly variant (H120G-Cu A Az) using the two-dimensional ESEEM technique, HYSCORE, is reported. HYSCORE spectroscopy has clarified conflicting results in previous electron paramagnetic resonance (EPR) and electron nuclear double resonance (ENDOR) studies and found clear differences between the two Cu A azurins. Specifically, a hyperfine coupling A N⊥ of 15.3 MHz was determined for the first time from the frequencies of double-quantum transitions of 14 N histidine nitrogens coordinated to Cu A in WT-Cu A Az. In contrast, such coupling was not observed in the spectra of H120G-Cu A Az, indicating at least a several megahertz increase in A N⊥ for the coordinated nitrogen in this variant. In addition, 14 N HYSCORE spectra of WT-Cu A Az show interaction with only one type of weakly coupled nitrogen assigned to the remote N ε atom of coordinated imidazole residues based on the quadrupole coupling constant ( e 2 Qq/4 h) of ∼0.4 MHz. The spectrum of H120G-Cu A Az resolves additional features typical for backbone peptide nitrogens with larger e 2 Qq/4 h values of ∼0.7 MHz. Hyperfine couplings with these nitrogens vary between ∼0.4 and 0.7 MHz. In addition, the two resolved cross-peaks from C β protons in H120G-Cu A Az display only ∼1 MHz shifts relative to the corresponding peaks in WT-Cu A Az. These new findings have provided the first experimental evidence of the previous density functional theory analysis that predicted changes in the delocalized electron spin population of ∼0.02-0.03 (i.e., ∼10%) on copper and sulfur atoms of the Cu A center in H120 variants relative to WT-Cu A Az and resolved contradicting results between EPR and ENDOR studies of the valence distribution in Cu A Az and its variants.

A machine learning approach is presented for analyzing complex two-dimensional hyperfine sublevel correlation electron paramagnetic resonance (HYSCORE EPR) spectra with the proficiency of an expert spectroscopist. The computer vision algorithm requires no training on experimental data; rather, all of the spin physics required to interpret the spectra are learned from simulations alone. This approach is therefore applicable even when insufficient experimental data exist to train the algorithm. The neural network is demonstrated to be capable of utilizing the full information content of two-dimensional 14 N HYSCORE spectra to predict the magnetic coupling parameters and their underlying probability distributions that were previously inaccessible. The predicted hyperfine ( a, T) and 14 N quadrupole ( K, η) coupling constants deviate from the previous manual analyses of the experimental spectra on average by 0.11 MHz, 0.09 MHz, 0.19 MHz, and 0.09, respectively.

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 13 C β -cysteine labeling of model [2Fe-2S] proteins to help bypass these problems. Analysis of the 13 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 13 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 1 -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.

Hydrogen bonding between semiquinone (SQ) intermediates and side-chain or backbone nitrogens in protein quinone processing sites (Q-sites) is a common motif. Previous studies on SQs from multiple protein environments have reported specific features in the 15 N HYSCORE spectra not reproducible by a theory based on fixed hyperfine parameters, and the source of these lineshape distortions remained unknown. In this work, using the spectra of the SQ in the Q-sites of wild-type and mutant D75H cytochrome bo 3 ubiquinol oxidase from Escherichia coli, we have explained the observed additional features as originating from a-strain of the isotropic hyperfine coupling. In two-dimensional spectra, the a-strain manifests as well-resolved lineshape distortions of the basic cross-ridges and accompanying lines of low intensity in the opposite quadrant that allow its direct analysis. We have shown that their appearance is regulated by the relative values of the strain width, Δ a, and parameter, δ = |2 a + T| - 4ν 15N . α-strain provides a direct measure of the structural dynamics and heterogeneity of the O···H···N bond in the SQ systems.

Iron-sulfur clusters are one of the most versatile and ancient classes of redox mediators in biology. The roles that these metal centers take on are predominantly determined by the number and types of coordinating ligands (typically cysteine and histidine) that modify the electronic structure of the cluster. Here we map the spin density distribution onto the cysteine ligands for the three major classes of the protein-bound, reduced [2Fe-2S](His) n (Cys) 4-n (n = 0, 1, 2) cluster by selective cysteine- 13 C β isotope labeling. The spin distribution is highly asymmetric in all three systems and delocalizes further along the reduced Fe 2+ ligands than the nonreducible Fe 3+ ligands for all clusters studied. The preferential spin transfer onto the chemically reactive Fe 2+ ligands is consistent with the structural concept that the orientation of the cluster in proteins is not arbitrarily decided, but rather is optimized such that it is likely to facilitate better electronic coupling with redox partners. The resolution of all cysteine- 13 C β hyperfine couplings and their assignments provides a measure of the relative covalencies of the metal-thiolate bonds not readily available to other techniques.

Determining the complete electron spin density distribution for protein-bound radicals, even with advanced pulsed electron paramagnetic resonance (EPR) methods, is a formidable task. Here we present a strategy to overcome this problem combining multifrequency HYSCORE and ENDOR measurements on site-specifically 13 C-labeled samples with DFT calculations on model systems. As a demonstration of this approach, pulsed EPR experiments are performed on the primary Q A and secondary Q B ubisemiquinones of the photosynthetic reaction center from Rhodobacter sphaeroides 13 C-labeled at the ring and tail positions. Despite the large number of nuclei interacting with the unpaired electron in these samples, two-dimensional X- and Q-band HYSCORE and orientation selective Q-band ENDOR resolve and allow for a characterization of the eight expected 13 C resonances from significantly different hyperfine tensors for both semiquinones. From these results we construct, for the first time, the most complete experimentally determined maps of the s- and p π -orbital spin density distributions for any protein organic cofactor radical to date. This work lays a foundation for understanding the relationship between the electronic structure of semiquinones and their functional properties, and introduces new techniques for mapping out the spin density distribution that are readily applicable to other systems.

Lysyl oxidase (LOX) is a copper amine oxidase that cross-links collagens and elastin in connective tissue and plays an important role in fibrosis, cancer development, and formation of the "metastatic niche". Despite its important biological functions, the structure of human LOX remains unknown (unlike that of an unrelated LOX, from Pichia pastoris). Here, we expressed active LOX from Drosophila melanogaster, DmLOXL1, a close homologue of human LOX, and characterized it by MS, UV-vis, activity, and inhibition assays. We then used bioinformatics, electron paramagnetic resonance, electron spin-echo envelope modulation, and hyperfine sublevel-correlation (HYSCORE) spectroscopies to probe Cu-ligand bonding finding direct evidence for pH-dependent Cu-His interactions. At pH = 9.3, the spectroscopic data indicated primarily a single His bound to Cu, but at pH = 7.5, there was evidence for a ∼ 1:1 mixture of species containing 1 and 3 His ligands. We then used HYSCORE to probe possible interactions between the LOX inhibitor BAPN (β-aminopropionitrile; 1-[ 13 C 15 N]cyano-2-aminoethane) and the copper center-finding none. Overall, the results are of interest since they provide new spectroscopic information about the nature of the catalytic site in LOX, an important anticancer drug target.

The Q-cycle mechanism of the bc 1 complex is now well enough understood to allow application of advanced computational approaches to the study of atomistic processes. In addition to the main features of the mechanism, these include control and gating of the bifurcated reaction at the Q o -site, through which generation of damaging reactive oxygen species is minimized. We report a new molecular dynamics model of the Rhodobacter sphaeroides bc 1 complex implemented in a native membrane, and constructed so as to eliminate blemishes apparent in earlier Rhodobacter models. Unconstrained MD simulations after equilibration with ubiquinol and ubiquinone respectively at Q o - and Q i -sites show that substrate binding configurations at both sites are different in important details from earlier models. We also demonstrate a new Q o -site intermediate, formed in the sub-ms time range, in which semiquinone remains complexed with the reduced iron sulfur protein. We discuss this, and a spring-loaded mechanism for modulating interactions of the iron sulfur protein with occupants of the Q o -site, in the context of control and gating roles. Such atomistic features of the mechanism can usefully be explored through simulation, but we stress the importance of constraints from physical chemistry and biology, both in setting up a simulation and in interpreting results.

The respiratory cytochrome bo 3 ubiquinol oxidase from Escherichia coli has a high-affinity ubiquinone binding site that stabilizes the one-electron reduced ubisemiquinone (SQ H ), which is a transient intermediate during the electron-mediated reduction of O 2 to water. It is known that SQ H is stabilized by two strong hydrogen bonds from R71 and D75 to ubiquinone carbonyl oxygen O1 and weak hydrogen bonds from H98 and Q101 to O4. In this work, SQ H was investigated with orientation-selective Q-band (∼34 GHz) pulsed 1 H electron-nuclear double resonance (ENDOR) spectroscopy on fully deuterated cytochrome (cyt) bo 3 in a H 2 O solvent so that only exchangeable protons contribute to the observed ENDOR spectra. Simulations of the experimental ENDOR spectra provided the principal values and directions of the hyperfine (hfi) tensors for the two strongly coupled H-bond protons (H1 and H2). For H1, the largest principal component of the proton anisotropic hfi tensor T z' = 11.8 MHz, whereas for H2, T z' = 8.6 MHz. Remarkably, the data show that the direction of the H1 H-bond is nearly perpendicular to the quinone plane (∼70° out of plane). The orientation of the second strong hydrogen bond, H2, is out of plane by ∼25°. Equilibrium molecular dynamics simulations on a membrane-embedded model of the cyt bo 3 Q H site show that these H-bond orientations are plausible but do not distinguish which H-bond, from R71 or D75, is nearly perpendicular to the quinone ring. Density functional theory calculations support the idea that the distances and geometries of the H-bonds to the ubiquinone carbonyl oxygens, along with the measured proton anisotropic hfi couplings, are most compatible with an anionic (deprotonated) ubisemiquinone.

Unlike photosystem II (PSII) in higher plants, bacterial photosynthetic reaction centers (bRCs) from Proteobacteria have an additional peripheral membrane subunit "H". The H subunit is necessary for photosynthetic growth, but can be removed chemically in vitro. The remaining LM dimer retains its activity to perform light-induced charge separation. Here we investigate the influence of the H subunit on interactions between the primary semiquinone and the protein matrix, using a combination of site-specific isotope labeling, pulsed electron paramagnetic resonance (EPR), and density functional theory (DFT) calculations. The data reveal substantially weaker binding interactions between the primary semiquinone and the LM dimer than observed for the intact bRC; the amount of electron spin transferred to the nitrogen hydrogen bond donors is significantly reduced, the methoxy groups are more free to rotate, and the spectra indicate a heterogeneous mixture of bound semiquinone states. These results are consistent with a loosening of the primary quinone binding pocket in the absence of the H subunit.

Cytochrome aa3-600 is a terminal oxidase in the electron transport pathway that contributes to the electrochemical membrane potential by actively pumping protons. A notable feature of this enzyme complex is that it uses menaquinol as its electron donor instead of cytochrome c when it reduces dioxygen to water. The enzyme stabilizes a menasemiquinone radical (SQ) at a high affinity site that is important for catalysis. One of the residues that interacts with the semiquinone is Arg70. We have made the R70H mutant and have characterized the menasemiquinone radical by advanced X- and Q-band EPR. The bound SQ of the R70H mutant exhibits a strong isotropic hyperfine coupling (a(14)N ≈ 2.0 MHz) with a hydrogen bonded nitrogen. This nitrogen originates from a histidine side chain, based on its quadrupole coupling constant, e(2)qQ/h = 1.44 MHz, typical for protonated imidazole nitrogens. In the wild-type cyt aa3-600, the SQ is instead hydrogen bonded with Nε from the Arg70 side chain. Analysis of the (1)H 2D electron spin echo envelope modulation (ESEEM) spectra shows that the mutation also changes the number and strength of the hydrogen bonds between the SQ and the surrounding protein. Despite the alterations in the immediate environment of the SQ, the R70H mutant remains catalytically active. These findings are in contrast to the equivalent mutation in the close homologue, cytochrome bo3 ubiquinol oxidase from Escherichia coli, where the R71H mutation eliminates function.

By utilizing a combined pulsed EPR and DFT approach, the high-resolution structure of the QB site semiquinone (SQB) was determined. The development of such a technique is crucial toward an understanding of protein-bound semiquinones on the structural level, as (i) membrane protein crystallography typically results in low resolution structures, and (ii) obtaining protein crystals in the semiquinone form is rarely feasible. The SQB hydrogen bond network was investigated with Q- (∼34 GHz) and X-band (∼9.7 GHz) pulsed EPR spectroscopy on fully deuterated reactions centers from Rhodobacter sphaeroides. Simulations in the SQB g-tensor reference frame provided the principal values and directions of the H-bond proton hyperfine tensors. Three protons were detected, one with an anisotropic tensor component, T = 4.6 MHz, assigned to the histidine NδH of His-L190, and two others with similar anisotropic constants T = 3.2 and 3.0 MHz assigned to the peptide NpH of Gly-L225 and Ile-L224, respectively. Despite the strong similarity in the peptide couplings, all hyperfine tensors were resolved in the Q-band ENDOR spectra. The Euler angles describing the series of rotations that bring the hyperfine tensors into the SQB g-tensor reference frame were obtained by least-squares fitting of the spectral simulations to the ENDOR data. These Euler angles show the locations of the hydrogen bonded protons with respect to the semiquinone. Our geometry optimized model of SQB used in previous DFT work is in strong agreement with the angular constraints from the spectral simulations, providing the foundation for future joint pulsed EPR and DFT semiquinone structural determinations in other proteins.

Ubiquinone forms an integral part of the electron transport chain in cellular respiration and photosynthesis across a vast number of organisms. Prior experimental results have shown that the photosynthetic reaction center (RC) from Rhodobacter sphaeroides is only fully functional with a limited set of methoxy-bearing quinones, suggesting that specific interactions with this substituent are required to drive electron transport and the formation of quinol. The nature of these interactions has yet to be determined. Through parameterization of a CHARMM-compatible quinone force field and subsequent molecular dynamics simulations of the quinone-bound RC, we have investigated and characterized the interactions of the protein with the quinones in the Q(A) and Q(B) sites using both equilibrium simulation and thermodynamic integration. In particular, we identify a specific interaction between the 2-methoxy group of ubiquinone in the Q(B) site and the amide nitrogen of GlyL225 that we implicate in locking the orientation of the 2-methoxy group, thereby tuning the redox potential difference between the quinones occupying the Q(A) and Q(B) sites. Disruption of this interaction leads to weaker binding in a ubiquinone analogue that lacks a 2-methoxy group, a finding supported by reverse electron transfer electron paramagnetic resonance experiments of the Q(A)⁻Q(B)⁻ biradical and competitive binding assays.

Specific isotopic labeling at the residue or substituent level extends the scope of different spectroscopic approaches to the atomistic level. Here we describe (13)C isotopic labeling of the methyl and methoxy ring substituents of ubiquinone, achieved through construction of a methionine auxotroph in Rhodobacter sphaeroides strain BC17 supplemented with l-methionine with the side chain methyl group (13)C-labeled. Two-dimensional electron spin echo envelope modulation (HYSCORE) was applied to study the (13)C methyl and methoxy hyperfine couplings in the semiquinone generated in situ at the Qi site of the bc1 complex in its membrane environment. The data were used to characterize the distribution of unpaired spin density and the conformations of the methoxy substituents based on density functional theory calculations of (13)C hyperfine tensors in the semiquinone of the geometry-optimized X-ray structure of the bc1 complex (Protein Data Bank entry 1PP9 ) with the highest available resolution. Comparison with other proteins indicates individual orientations of the methoxy groups in each particular case are always different from the methoxy conformations in the anion radical prepared in a frozen alcohol solution. The protocol used in the generation of the methionine auxotroph is more generally applicable and, because it introduces a gene deletion using a suicide plasmid, can be applied repeatedly.

Recent studies have shown that only quinones with a 2-methoxy group can act simultaneously as the primary (Q A ) and secondary (Q B ) electron acceptors in photosynthetic reaction centers from purple bacteria such as Rb. sphaeroides . 13 C HYSCORE measurements of the 2-methoxy group in the semiquinone states, SQ A and SQ B , were compared with DFT calculations of the 13 C hyperfine couplings as a function of the 2-methoxy dihedral angle. X-ray structure comparisons support 2-methoxy dihedral angle assignments corresponding to a redox potential gap (Δ E m ) between Q A and Q B of 175-193 mV. A model having a methyl group substituted for the 2-methoxy group exhibits no electron affinity difference. This is consistent with the failure of a 2-methyl ubiquinone analogue to function as Q B in mutant reaction centers with a Δ E m of ∼160-195 mV. The conclusion reached is that the 2-methoxy group is the principal determinant of electron transfer from Q A to Q B in type II photosynthetic reaction centers with ubiquinone serving as both acceptor quinones.

X- and Q-band pulsed EPR spectroscopy was applied to study the interaction of the QA site semiquinone (SQA) with nitrogens from the local protein environment in natural abundance (14)N and in (15)N uniformly labeled photosynthetic reaction centers of Rhodobacter sphaeroides. The hyperfine and nuclear quadrupole tensors for His-M219 Nδ and Ala-M260 peptide nitrogen (Np) were estimated through simultaneous simulation of the Q-band (15)N Davies ENDOR, X- and Q-band (14,15)N HYSCORE, and X-band (14)N three-pulse ESEEM spectra, with support from DFT calculations. The hyperfine coupling constants were found to be a((14)N) = 2.3 MHz, T = 0.3 MHz for His-M219 Nδ and a((14)N) = 2.6 MHz, T = 0.3 MHz for Ala-M260 Np. Despite that His-M219 Nδ is established as the stronger of the two H-bond donors, Ala-M260 Np is found to have the larger value of a((14)N). The nuclear quadrupole coupling constants were estimated as e(2)Qq/4h = 0.38 MHz, η = 0.97 and e(2)Qq/4h = 0.74 MHz, η = 0.59 for His-M219 Nδ and Ala-M260 Np, respectively. An analysis of the available data on nuclear quadrupole tensors for imidazole nitrogens found in semiquinone-binding proteins and copper complexes reveals these systems share similar electron occupancies of the protonated nitrogen orbitals. By applying the Townes-Dailey model, developed previously for copper complexes, to the semiquinones, we find the asymmetry parameter η to be a sensitive probe of the histidine Nδ-semiquinone hydrogen bond strength. This is supported by a strong correlation observed between η and the isotropic coupling constant a((14)N) and is consistent with previous computational works and our own semiquinone-histidine model calculations. The empirical relationship presented here for a((14)N) and η will provide an important structural characterization tool in future studies of semiquinone-binding proteins.

The secondary quinone anion radical QB(-) (SQB) in reaction centers of Rhodobacter sphaeroides interacts with Nδ of His-L190 and Np (peptide nitrogen) of Gly-L225 involved in hydrogen bonds to the QB carbonyls. In this work, S-band (∼3.6 GHz) ESEEM was used with the aim of obtaining a complete characterization of the nuclear quadrupole interaction (nqi) tensors for both nitrogens by approaching the cancelation condition between the isotropic hyperfine coupling and (14)N Zeeman frequency at lower microwave frequencies than traditional X-band (9.5 GHz). By performing measurements at S-band, we found a dominating contribution of Nδ in the form of a zero-field nqi triplet at 0.55, 0.92, and 1.47 MHz, defining the quadrupole coupling constant K = e(2)qQ/4h = 0.4 MHz and associated asymmetry parameter η = 0.69. Estimates of the hyperfine interaction (hfi) tensors for Nδ and Np were obtained from simulations of 1D and 2D (14,15)N X-band and three-pulse (14)N S-band spectra with all nuclear tensors defined in the SQB g-tensor coordinate system. From simulations, we conclude that the contribution of Np to the S-band spectrum is suppressed by its strong nqi and weak isotropic hfi comparable to the level of hyperfine anisotropy, despite the near-cancelation condition for Np at S-band. The excellent agreement between our EPR simulations and DFT calculations of the nitrogen hfi and nqi tensors to SQB is promising for the future application of powder ESEEM to full tensor characterizations.

Only quinones with a 2-methoxy group can act simultaneously as the primary (QA) and secondary (QB) electron acceptors in photosynthetic reaction centers from Rhodobacter sphaeroides. (13)C hyperfine sublevel correlation measurements of the 2-methoxy in the semiquinone states, SQA and SQB, were compared with quantum mechanics calculations of the (13)C couplings as a function of the dihedral angle. X-ray structures support dihedral angle assignments corresponding to a redox potential gap (ΔEm) between QA and QB of ~180 mV. This is consistent with the failure of a ubiquinone analogue lacking the 2-methoxy to function as QB in mutant reaction centers with a ΔEm of ≈160-195 mV.

We describe a 1D and 2D ESEEM investigation of VO 2+ adsorbed on hydroxyapatite (HA) at different concentrations and compare with VO 2+ -triphosphate (TPH) complexes studied previously in detail, in an effort to provide more insight into the structure of VO 2+ coordination in bone. Structures of this interaction are important because of the role of bone in the long-term storage of administered vanadium, and the likely role of bone in the steady-state release of vanadium leading to the chronic insulin-enhancing anti-diabetic effects of vanadyl complexes. Three similar sets of cross-peaks from phosphorus nuclei observed in the 31 P HYSCORE spectra of VO 2+ -HA, VO 2+ -TPH, and VO 2+ -bone suggest a common tridentate binding motif for triphosphate moieties to the vanadyl ion. The similarities between the systems present the possibility that in vivo vanadyl coordination in bone is relatively uniform. Experiments with HA samples containing different amounts of adsorbed VO 2+ demonstrate additional peculiarities of the ion-adsorbent interaction which can be expected in vivo . HYSCORE spectra of HA samples show varying relative intensities of 31 P lines from phosphate ligands and 1 H lines, especially lines from protons of coordinated water molecules. This result suggests that the number of equatorial phosphate ligands in HA could be different depending on the water content of the sample and the VO 2+ concentration; complexes of different structure probably contribute to the spectra of VO 2+ -HA. Similar behavior can be also expected in vivo during VO 2+ accumulation in bones.

Ubiquinone is an almost universal, membrane-associated redox mediator. Its ability to accept either one or two electrons allows it to function in critical roles in biological electron transport. The redox properties of ubiquinone in vivo are determined by its environment in the binding sites of proteins and by the dihedral angle of each methoxy group relative to the ring plane. This is an attribute unique to ubiquinone among natural quinones and could account for its widespread function with many different redox complexes. In this work, we use the photosynthetic reaction center as a model system for understanding the role of methoxy conformations in determining the redox potential of the ubiquinone/semiquinone couple. Despite the abundance of X-ray crystal structures for the reaction center, quinone site resolution has thus far been too low to provide a reliable measure of the methoxy dihedral angles of the primary and secondary quinones, QA and QB. We performed 2D ESEEM (HYSCORE) on isolated reaction centers with ubiquinones (13)C-labeled at the headgroup methyl and methoxy substituents, and have measured the (13)C isotropic and anisotropic components of the hyperfine tensors. Hyperfine couplings were compared to those derived by DFT calculations as a function of methoxy torsional angle allowing estimation of the methoxy dihedral angles for the semiquinones in the QA and QB sites. Based on this analysis, the orientation of the 2-methoxy groups are distinct in the two sites, with QB more out of plane by 20-25°. This corresponds to an ≈50 meV larger electron affinity for the QB quinone, indicating a substantial contribution to the experimental difference in redox potentials (60-75 mV) of the two quinones. The methods developed here can be readily extended to ubiquinone-binding sites in other protein complexes.

There are two homologous membrane-bound enzymes in Escherichia coli that catalyze reversible conversion between succinate/fumarate and quinone/quinol. Succinate:ubiquinone reductase (SQR) is a component of aerobic respiratory chains, whereas quinol:fumarate reductase (QFR) utilizes menaquinol to reduce fumarate in a final step of anaerobic respiration. Although, both protein complexes are capable of supporting bacterial growth on either minimal succinate or fumarate media, the enzymes are more proficient in their physiological directions. Here we evaluate factors that may underlie this catalytic bias. This article is part of a Special Issue entitled: Respiratory complex II: Role in cellular physiology and disease., (Copyright © 2013 Elsevier B.V. All rights reserved.)

The in vitro evolution and engineering of quorum-quenching lactonases with enhanced reactivities was achieved using a thermostable GKL enzyme as a template, yielding the E101G/R230C GKL mutant with increased catalytic activity and a broadened substrate range [Chow, J. Y., Xue, B., Lee, K. H., Tung, A., Wu, L., Robinson, R. C., and Yew, W. S. (2010) J. Biol. Chem. 285, 40911-40920]. This enzyme possesses the (β/α)8-barrel fold and is a member of the PLL (phosphotriesterase-like lactonase) group of enzymes within the amidohydrolase superfamily that hydrolyze N-acyl-homoserine lactones, which mediate the quorum-sensing pathways of bacteria. The structure of the evolved N-butyryl-l-homoserine lactone (substrate)-bound E101G/R230C GKL enzyme was determined, in the presence of the inactivating D266N mutation, to a resolution of 2.2 Å to provide an explanation for the observed rate enhancements. In addition, the substrate-bound structure of the catalytically inactive E101N/D266N mutant of the manganese-reconstituted enzyme was determined to a resolution of 2.1 Å and the structure of the ligand-free, manganese-reconstituted E101N mutant to a resolution of 2.6 Å, and the structures of ligand-free zinc-reconstituted wild-type, E101N, R230D, and E101G/R230C mutants of GKL were determined to resolutions of 2.1, 2.1, 1.9, and 2.0 Å, respectively. In particular, the structure of the evolved E101G/R230C mutant of GKL provides evidence of a catalytically productive active site architecture that contributes to the observed enhancement of catalysis. At high concentrations, wild-type and mutant GKL enzymes are differentially colored, with absorbance maxima in the range of 512-553 nm. The structures of the wild-type and mutant GKL provide a tractable link between the origins of the coloration and the charge-transfer complex between the α-cation and Tyr99 within the enzyme active site. Taken together, this study provides evidence of the modulability of enzymatic catalysis through subtle changes in enzyme active site architecture.

The electronic structure and geometry of redox-active metal cofactors in proteins are tuned by the pattern of hydrogen bonding with the backbone peptide matrix. In this study we developed a method for selective amino acid labeling of a hyperthermophilic archaeal metalloprotein with engineered Escherichia coli auxotroph strains, and we applied this to resolve the hydrogen bond interactions with the reduced Rieske-type [2Fe-2S] cluster by two-dimensional pulsed electron spin resonance technique. Because deep electron spin-echo envelope modulation of two histidine (14)N(δ) ligands of the cluster decreased non-coordinating (15)N signal intensities via the cross-suppression effect, an inverse labeling strategy was employed in which (14)N amino acid-labeled archaeal Rieske-type ferredoxin samples were examined in an (15)N-protein background. This has directly identified Lys45 N(α) as providing the major pathway for the transfer of unpaired electron spin density from the reduced cluster by a "through-bond" mechanism. All other backbone peptide nitrogens interact more weakly with the reduced cluster. The extension of this approach will allow visualizing the three-dimensional landscape of preferred pathways for the transfer of unpaired spin density from a paramagnetic metal center onto the protein frame, and will discriminate specific interactions by a "through-bond" mechanism from interactions which are "through-space" in various metalloproteins.

In the Q(B) site of the Rhodobacter sphaeroides photosynthetic reaction center, the donation of a hydrogen bond from the hydroxyl group of Ser-L223 to the ubisemiquinone formed after the first flash is debatable. In this study, we use a combination of spectroscopy and quantum mechanics/molecular mechanics (QM/MM) calculations to comprehensively explore this topic. We show that ENDOR, ESEEM, and HYSCORE spectroscopic differences between mutant L223SA and the wild-type sample (WT) are negligible, indicating only minor perturbations in the ubisemiquinone spin density for the mutant sample. Qualitatively, this suggests that a strong hydrogen bond does not exist in the WT between the Ser-L223 hydroxyl group and the semiquinone O(1) atom, as removal of this hydrogen bond in the mutant should cause a significant redistribution of spin density in the semiquinone. We show quantitatively, using QM/MM calculations, that a WT model in which the Ser-L223 hydroxyl group is rotated to prevent hydrogen bond formation with the O(1) atom of the semiquinone predicts negligible change for the L223SA mutant. This, together with the better agreement between key QM/MM calculated and experimental hyperfine couplings for the non-hydrogen-bonded model, leads us to conclude that no strong hydrogen bond is formed between the Ser-L223 hydroxyl group and the semiquinone O(1) atom after the first flash. The implications of this finding for quinone reduction in photosynthetic reaction centers are discussed.

Selective (15)N isotope labeling of the cytochrome bo(3) ubiquinol oxidase from Escherichia coli with auxotrophs was used to characterize the hyperfine couplings with the side-chain nitrogens from residues R71, H98, and Q101 and peptide nitrogens from residues R71 and H98 around the semiquinone (SQ) at the high-affinity Q(H) site. The two-dimensional ESEEM (HYSCORE) data have directly identified N(ε) of R71 as an H-bond donor carrying the largest amount of unpaired spin density. In addition, weaker hyperfine couplings with the side-chain nitrogens from all residues around the SQ were determined. These hyperfine couplings reflect a distribution of the unpaired spin density over the protein in the SQ state of the Q(H) site and the strength of interaction with different residues. The approach was extended to the virtually inactive D75H mutant, where the intermediate SQ is also stabilized. We found that N(ε) of a histidine residue, presumably H75, carries most of the unpaired spin density instead of N(ε) of R71, as in wild-type bo(3). However, the detailed characterization of the weakly coupled (15)N atoms from selective labeling of R71 and Q101 in D75H was precluded by overlap of the (15)N lines with the much stronger ~1.6 MHz line from the quadrupole triplet of the strongly coupled (14)N(ε) atom of H75. Therefore, a reverse labeling approach, in which the enzyme was uniformly labeled except for selected amino acid types, was applied to probe the contribution of R71 and Q101 to the (15)N signals. Such labeling has shown only weak coupling with all nitrogens of R71 and Q101. We utilize density functional theory-based calculations to model the available information about (1)H, (15)N, and (13)C hyperfine couplings for the Q(H) site and to describe the protein-substrate interactions in both enzymes. In particular, we identify the factors responsible for the asymmetric distribution of the unpaired spin density and ponder the significance of this asymmetry to the quinone's electron transfer function.

Amino-acid selective isotope labeling of proteins offers numerous advantages in mechanistic studies by revealing structural and functional information unattainable from a crystallographic approach. However, efficient labeling of proteins with selected amino acids necessitates auxotrophic hosts, which are often not available. We have constructed a set of auxotrophs in a commonly used Escherichia coli expression strain C43(DE3), a derivative of E. coli BL21(DE3), which can be used for isotopic labeling of individual amino acids or sets of amino acids. These strains have general applicability to either soluble or membrane proteins that can be expressed in E. coli. We present examples in which proteins are selectively labeled with (13)C- and (15)N-amino acids and studied using magic-angle spinning solid-state NMR and pulsed EPR, demonstrating the utility of these strains for biophysical characterization of membrane proteins, radical-generating enzymes and metalloproteins., (Copyright © 2011 Elsevier Inc. All rights reserved.)

When the superoxide radical O(2)(•-) is generated on reaction of KO(2) with water in dimethyl sulfoxide, the decay of the radical is dramatically accelerated by inclusion of quinones in the reaction mix. For quinones with no or short hydrophobic tails, the radical product is a semiquinone at much lower yield, likely indicating reduction of quinone by superoxide and loss of most of the semiquinone product by disproportionation. In the presence of ubiquinone-10, a different species (I) is generated, which has the EPR spectrum of superoxide radical. However, pulsed EPR shows spin interaction with protons in fully deuterated solvent, indicating close proximity to the ubinquinone-10. We discuss the nature of species I, and possible roles in the physiological reactions through which ubisemiquinone generates superoxide by reduction of O(2) through bypass reactions in electron transfer chains.

Enzymatic N(2) reduction proceeds along a reaction pathway composed of a sequence of intermediate states generated as a dinitrogen bound to the active-site iron-molybdenum cofactor (FeMo-co) of the nitrogenase MoFe protein undergoes six steps of hydrogenation (e(-)/H(+) delivery). There are two competing proposals for the reaction pathway, and they invoke different intermediates. In the 'Distal' (D) pathway, a single N of N(2) is hydrogenated in three steps until the first NH(3) is liberated, and then the remaining nitrido-N is hydrogenated three more times to yield the second NH(3). In the 'Alternating' (A) pathway, the two N's instead are hydrogenated alternately, with a hydrazine-bound intermediate formed after four steps of hydrogenation and the first NH(3) liberated only during the fifth step. A recent combination of X/Q-band EPR and (15)N, (1,2)H ENDOR measurements suggested that states trapped during turnover of the α-70(Ala)/α-195(Gln) MoFe protein with diazene or hydrazine as substrate correspond to a common intermediate (here denoted I) in which FeMo-co binds a substrate-derived [N(x)H(y)] moiety, and measurements reported here show that turnover with methyldiazene generates the same intermediate. In the present report we describe X/Q-band EPR and (14/15)N, (1,2)H ENDOR/HYSCORE/ESEEM measurements that characterize the N-atom(s) and proton(s) associated with this moiety. The experiments establish that turnover with N(2)H(2), CH(3)N(2)H, and N(2)H(4) in fact generates a common intermediate, I, and show that the N-N bond of substrate has been cleaved in I. Analysis of this finding leads us to conclude that nitrogenase reduces N(2)H(2), CH(3)N(2)H, and N(2)H(4) via a common A reaction pathway, and that the same is true for N(2) itself, with Fe ion(s) providing the site of reaction.

In the photosynthetic reaction center from Rhodobacter sphaeroides, the primary (Q(A)) and secondary (Q(B)) electron acceptors are both ubiquinone-10, but with very different properties and functions. To investigate the protein environment that imparts these functional differences, we have applied X-band HYSCORE, a 2D pulsed EPR technique, to characterize the exchangeable protons around the semiquinone (SQ) in the Q(A) and Q(B) sites, using samples of (15)N-labeled reaction centers, with the native high spin Fe(2+) exchanged for diamagnetic Zn(2+), prepared in (1)H(2)O and (2)H(2)O solvent. The powder HYSCORE method is first validated against the orientation-selected Q-band ENDOR study of the Q(A) SQ by Flores et al. (Biophys. J.2007, 92, 671-682), with good agreement for two exchangeable protons with anisotropic hyperfine tensor components, T, both in the range 4.6-5.4 MHz. HYSCORE was then applied to the Q(B) SQ where we found proton lines corresponding to T ≈ 5.2, 3.7 MHz and T ≈ 1.9 MHz. Density functional-based quantum mechanics/molecular mechanics (QM/MM) calculations, employing a model of the Q(B) site, were used to assign the observed couplings to specific hydrogen bonding interactions with the Q(B) SQ. These calculations allow us to assign the T = 5.2 MHz proton to the His-L190 N(δ)H···O(4) (carbonyl) hydrogen bonding interaction. The T = 3.7 MHz spectral feature most likely results from hydrogen bonding interactions of O1 (carbonyl) with both Gly-L225 peptide NH and Ser-L223 hydroxyl OH, which possess calculated couplings very close to this value. The smaller 1.9 MHz coupling is assigned to a weakly bound peptide NH proton of Ile-L224. The calculations performed with this structural model of the Q(B) site show less asymmetric distribution of unpaired spin density over the SQ than seen for the Q(A) site, consistent with available experimental data for (13)C and (17)O carbonyl hyperfine couplings. The implications of these interactions for Q(B) function and comparisons with the Q(A) site are discussed., (© 2011 American Chemical Society)

The cytochrome bo(3) ubiquinol oxidase from Escherichia coli resides in the bacterial cytoplasmic membrane and catalyzes the two-electron oxidation of ubiquinol-8 and four-electron reduction of O(2) to water. The one-electron reduced semiquinone forms transiently during the reaction, and the enzyme has been demonstrated to stabilize the semiquinone. The semiquinone is also formed in the D75E mutant, where the mutation has little influence on the catalytic activity, and in the D75H mutant, which is virtually inactive. In this work, wild-type cytochrome bo(3) as well as the D75E and D75H mutant proteins were prepared with ubiquinone-8 (13)C-labeled selectively at the methyl and two methoxy groups. This was accomplished by expressing the proteins in a methionine auxotroph in the presence of l-methionine with the side chain methyl group (13)C-labeled. The (13)C-labeled quinone isolated from cytochrome bo(3) was also used for the generation of model anion radicals in alcohol. Two-dimensional pulsed EPR and ENDOR were used for the study of the (13)C methyl and methoxy hyperfine couplings in the semiquinone generated in the three proteins indicated above and in the model system. The data were used to characterize the transferred unpaired spin densities on the methyl and methoxy substituents and the conformations of the methoxy groups. In the wild type and D75E mutant, the constraints on the configurations of the methoxy side chains are similar, but the D75H mutant appears to have altered methoxy configurations, which could be related to the perturbed electron distribution in the semiquinone and the loss of enzymatic activity.

Cytochrome bo(3) is the major respiratory oxidase located in the cytoplasmic membrane of Escherichia coli when grown under high oxygen tension. The enzyme catalyzes the 2-electron oxidation of ubiquinol-8 and the 4-electron reduction of dioxygen to water. When solubilized and isolated using dodecylmaltoside, the enzyme contains one equivalent of ubiquinone-8, bound at a high affinity site (Q(H)). The quinone bound at the Q(H) site can form a stable semiquinone, and the amino acid residues which hydrogen bond to the semiquinone have been identified. In the current work, it is shown that the tightly bound ubiquinone-8 at the Q(H) site is not displaced by ubiquinol-1 even during enzyme turnover. Furthermore, the presence of high affinity inhibitors, HQNO and aurachin C1-10, does not displace ubiquinone-8 from the Q(H) site. The data clearly support the existence of a second binding site for ubiquinone, the Q(L) site, which can rapidly exchange with the substrate pool. HQNO is shown to bind to a single site on the enzyme and to prevent formation of the stable ubisemiquinone, though without displacing the bound quinone. Inhibition of the steady state kinetics of the enzyme indicates that aurachin C1-10 may compete for binding with quinol at the Q(L) site while, at the same time, preventing formation of the ubisemiquinone at the Q(H) site. It is suggested that the two quinone binding sites may be adjacent to each other or partially overlap., (Copyright © 2010 Elsevier B.V. All rights reserved.)

Photosynthetic reaction centers from Rhodobacter sphaeroides have identical ubiquinone-10 molecules functioning as primary (Q(A)) and secondary (Q(B)) electron acceptors. X-band 2D pulsed EPR spectroscopy, called HYSCORE, was applied to study the interaction of the Q(B) site semiquinone with nitrogens from the local protein environment in natural and (15)N uniformly labeled reactions centers. (14)N and (15)N HYSCORE spectra of the Q(B) semiquinone show the interaction with two nitrogens carrying transferred unpaired spin density. Quadrupole coupling constants estimated from (14)N HYSCORE spectra indicate them to be a protonated nitrogen of an imidazole residue and amide nitrogen of a peptide group. (15)N HYSCORE spectra allowed estimation of the isotropic and anisotropic couplings with these nitrogens. From these data, we calculated the unpaired spin density transferred onto 2s and 2p orbitals of nitrogen and analyzed the contribution of different factors to the anisotropic hyperfine tensors. The hyperfine coupling of other protein nitrogens with the semiquinone is weak (<0.1 MHz). These results clearly indicate that the Q(B) semiquinone forms hydrogen bonds with two nitrogens and provide quantitative characteristics of the hyperfine couplings with these nitrogens, which can be used in theoretical modeling of the Q(B) site. On the basis of the quadrupole coupling constant, one nitrogen can only be assigned to N(delta) of His-L190, consistent with all existing structures. However, we cannot specify between two candidates the residue corresponding to the second nitrogen. Further work employing multifrequency spectroscopic approaches or selective isotope labeling would be desirable for unambiguous assignment of this nitrogen.

Cytochrome aa(3)-600 is one of the principle respiratory oxidases from Bacillus subtilis and is a member of the heme-copper superfamily of oxygen reductases. This enzyme catalyzes the two-electron oxidation of menaquinol and the four-electron reduction of O(2) to 2H(2)O. Cytochrome aa(3)-600 is of interest because it is a very close homologue of the cytochrome bo(3) ubiquinol oxidase from Escherichia coli, except that it uses menaquinol instead of ubiquinol as a substrate. One question of interest is how the proteins differ in response to the differences in structure and electrochemical properties between ubiquinol and menaquinol. Cytochrome bo(3) has a high affinity binding site for ubiquinol that stabilizes a ubi-semiquinone. This has permitted the use of pulsed EPR techniques to investigate the protein interaction with the ubiquinone. The current work initiates studies to characterize the equivalent site in cytochrome aa(3)-600. Cytochrome aa(3)-600 has been cloned and expressed in a His-tagged form in B. subtilis. After isolation of the enzyme in dodecylmaltoside, it is shown that the pure enzyme contains 1 eq of menaquinone-7 and that the enzyme stabilizes a mena-semiquinone. Pulsed EPR studies have shown that there are both similarities as well as significant differences in the interactions of the mena-semiquinone with cytochrome aa(3)-600 in comparison with the ubi-semiquinone in cytochrome bo(3). Our data indicate weaker hydrogen bonds of the menaquinone in cytochrome aa(3)-600 in comparison with ubiquinone in cytochrome bo(3). In addition, the electronic structure of the semiquinone cyt aa(3)-600 is more shifted toward the anionic form from the neutral state in cyt bo(3).

The rate-determining step in the overall turnover of the bc(1) complex is electron transfer from ubiquinol to the Rieske iron-sulfur protein (ISP) at the Q(o)-site. Structures of the ISP from Rhodobacter sphaeroides show that serine 154 and tyrosine 156 form H-bonds to S-1 of the [2Fe-2S] cluster and to the sulfur atom of the cysteine liganding Fe-1 of the cluster, respectively. These are responsible in part for the high potential (E(m)(,7) approximately 300 mV) and low pK(a) (7.6) of the ISP, which determine the overall reaction rate of the bc(1) complex. We have made site-directed mutations at these residues, measured thermodynamic properties using protein film voltammetry to evaluate the E(m) and pK(a) values of ISPs, explored the local proton environment through two-dimensional electron spin echo envelope modulation, and characterized function in strains S154T, S154C, S154A, Y156F, and Y156W. Alterations in reaction rate were investigated under conditions in which concentration of one substrate (ubiquinol or ISP(ox)) was saturating and the other was varied, allowing calculation of kinetic terms and relative affinities. These studies confirm that H-bonds to the cluster or its ligands are important determinants of the electrochemical characteristics of the ISP, likely through electron affinity of the interacting atom and the geometry of the H-bonding neighborhood. The calculated parameters were used in a detailed Marcus-Brønsted analysis of the dependence of rate on driving force and pH. The proton-first-then-electron model proposed accounts naturally for the effects of mutation on the overall reaction.

Two-dimensional electron spin-echo envelope modulation (ESEEM) analysis of the uniformly (15)N-labeled archaeal Rieske-type [2Fe-2S] ferredoxin (ARF) from Sulfolobus solfataricus P1 has been conducted in comparison with the previously characterized high-potential protein homologs. Major differences among these proteins were found in the hyperfine sublevel correlation (HYSCORE) lineshapes and intensities of the signals in the (++) quadrant, which are contributed from weakly coupled (non-coordinated) peptide nitrogens near the reduced clusters. They are less pronounced in the HYSCORE spectra of ARF than those of the high-potential protein homologs, and may account for the tuning of Rieske-type clusters in various redox systems.

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.

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.

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.

The selective (15)N isotope labeling was used for the identification of the nitrogen involved in a hydrogen bond formation with the semiquinone in the high-affinity Q(H) site in the cytochrome bo(3) ubiquinol oxidase. This nitrogen produces dominating contribution to X-Band (14)N ESEEM spectra. The 2D ESEEM (HYSCORE) experiments with the Q(H) site SQ in the series of selectively (15)N labeled bo(3) oxidase proteins have directly identified the N(epsilon) of R71 as an H-bond donor. In addition, selective (15)N labeling has allowed us for the first time to determine weak hyperfine couplings with the side-chain nitrogens from all residues around the SQ. Those are reflecting a distribution of the unpaired spin density over the protein in the SQ state of the quinone processing site.

Recent progress in understanding the Q-cycle mechanism of the bc(1) complex is reviewed. The data strongly support a mechanism in which the Q(o)-site operates through a reaction in which the first electron transfer from ubiquinol to the oxidized iron-sulfur protein is the rate-determining step for the overall process. The reaction involves a proton-coupled electron transfer down a hydrogen bond between the ubiquinol and a histidine ligand of the [2Fe-2S] cluster, in which the unfavorable protonic configuration contributes a substantial part of the activation barrier. The reaction is endergonic, and the products are an unstable ubisemiquinone at the Q(o)-site, and the reduced iron-sulfur protein, the extrinsic mobile domain of which is now free to dissociate and move away from the site to deliver an electron to cyt c(1) and liberate the H(+). When oxidation of the semiquinone is prevented, it participates in bypass reactions, including superoxide generation if O(2) is available. When the b-heme chain is available as an acceptor, the semiquinone is oxidized in a process in which the proton is passed to the glutamate of the conserved -PEWY- sequence, and the semiquinone anion passes its electron to heme b(L) to form the product ubiquinone. The rate is rapid compared to the limiting reaction, and would require movement of the semiquinone closer to heme b(L) to enhance the rate constant. The acceptor reactions at the Q(i)-site are still controversial, but likely involve a "two-electron gate" in which a stable semiquinone stores an electron. Possible mechanisms to explain the cyt b(150) phenomenon are discussed, and the information from pulsed-EPR studies about the structure of the intermediate state is reviewed. The mechanism discussed is applicable to a monomeric bc(1) complex. We discuss evidence in the literature that has been interpreted as shown that the dimeric structure participates in a more complicated mechanism involving electron transfer across the dimer interface. We show from myxothiazol titrations and mutational analysis of Tyr-199, which is at the interface between monomers, that no such inter-monomer electron transfer is detected at the level of the b(L) hemes. We show from analysis of strains with mutations at Asn-221 that there are coulombic interactions between the b-hemes in a monomer. The data can also be interpreted as showing similar coulombic interaction across the dimer interface, and we discuss mechanistic implications.

The ubisemiquinone stabilized at the Qi-site of the bc1 complex of Rhodobacter sphaeroides forms a hydrogen bond with a nitrogen from the local protein environment, tentatively identified as ring N from His-217. The interactions of 14N and 15N have been studied by X-band (approximately 9.7 GHz) and S-band (3.4 GHz) pulsed EPR spectroscopy. The application of S-band spectroscopy has allowed us to determine the complete nuclear quadrupole tensor of the 14N involved in H-bond formation and to assign it unambiguously to the Nepsilon of His-217. This tensor has distinct characteristics in comparison with H-bonds between semiquinones and Ndelta in other quinone-processing sites. The experiments with 15N showed that the Nepsilon of His-217 was the only nitrogen carrying any considerable unpaired spin density in the ubiquinone environment, and allowed calculation of the isotropic and anisotropic couplings with the Nepsilon of His-217. From these data, we could estimate the unpaired spin density transferred onto 2s and 2p orbitals of nitrogen and the distance from the nitrogen to the carbonyl oxygen of 2.38+/-0.13A. The hyperfine coupling of other protein nitrogens with semiquinone is <0.1 MHz. This did not exclude the nitrogen of the Asn-221 as a possible hydrogen bond donor to the methoxy oxygen of the semiquinone. A mechanistic role for this residue is supported by kinetic experiments with mutant strains N221T, N221H, N221I, N221S, N221P, and N221D, all of which showed some inhibition but retained partial turnover.

The cytochrome bo3 ubiquinol oxidase catalyzes the two-electron oxidation of ubiquinol in the cytoplasmic membrane of Escherichia coli, and reduces O2 to water. This enzyme has a high affinity quinone binding site (QH), and the quinone bound to this site acts as a cofactor, necessary for rapid electron transfer from substrate ubiquinol, which binds at a separate site (QL), to heme b. Previous pulsed EPR studies have shown that a semiquinone at the QH site formed during the catalytic cycle is a neutral species, with two strong hydrogen bonds to Asp-75 and either Arg-71 or Gln-101. In the current work, pulsed EPR studies have been extended to two mutants at the QH site. The D75E mutation has little influence on the catalytic activity, and the pattern of hydrogen bonding is similar to the wild type. In contrast, the D75H mutant is virtually inactive. Pulsed EPR revealed significant structural changes in this mutant. The hydrogen bond to Arg-71 or Gln-101 that is present in both the wild type and D75E mutant oxidases is missing in the D75H mutant. Instead, the D75H has a single, strong hydrogen bond to a histidine, likely His-75. The D75H mutant stabilizes an anionic form of the semiquinone as a result of the altered hydrogen bond network. Either the redistribution of charge density in the semiquinone species, or the altered hydrogen bonding network is responsible for the loss of catalytic function.

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.