Exploring by pulsed EPR the electronic structure of ubisemiquinone bound at the QH site of cytochrome bo3 from Escherichia coli with in vivo 13C-labeled methyl and methoxy substituents

Evolution of quinol oxidation within the heme‑copper oxidoreductase superfamily

The heme‑copper oxidoreductase (HCO) superfamily is a large superfamily of terminal respiratory enzymes that are widely distributed across the three domains of life. The superfamily includes biochemically diverse oxygen reductases and nitric oxide reductases that are pivotal in the pathways of aerobic respiration and denitrification. The adaptation of HCOs to use quinol as the electron donor instead of cytochrome c has significant implication for the respiratory flexibility and energetic efficiency of the respiratory chains that include them. In this work, we explore the adaptation of this scaffold to two different electron donors, cytochromes c and quinols, with extensive sequence analysis of these enzymes from publicly available datasets. Our work shows that quinol oxidation evolved independently within the HCO superfamily at least seven times. Enzymes from only two of these independently evolved clades have been biochemically well-characterized. Combining structural modeling with sequence analysis, we identify putative quinol binding sites in each of the novel quinol oxidases. Our analysis of experimental and modeling data suggests that the quinol binding site appears to have evolved at the same structural position within the scaffold more than once.

Structure of the cytochrome aa 3 -600 heme-copper menaquinol oxidase bound to inhibitor HQNO shows TM0 is part of the quinol binding site

Virtually all proton-pumping terminal respiratory oxygen reductases are members of the heme-copper oxidoreductase superfamily. Most of these enzymes use reduced cytochrome c as a source of electrons, but a group of enzymes have evolved to directly oxidize membrane-bound quinols, usually menaquinol or ubiquinol. All of the quinol oxidases have an additional transmembrane helix (TM0) in subunit I that is not present in the related cytochrome c oxidases. The current work reports the 3.6-Å-resolution X-ray structure of the cytochrome aa ₃ -600 menaquinol oxidase from Bacillus subtilis containing 1 equivalent of menaquinone. The structure shows that TM0 forms part of a cleft to accommodate the menaquinol-7 substrate. Crystals which have been soaked with the quinol-analog inhibitor HQNO (N-oxo-2-heptyl-4-hydroxyquinoline) or 3-iodo-HQNO reveal a single binding site where the inhibitor forms hydrogen bonds to amino acid residues shown previously by spectroscopic methods to interact with the semiquinone state of menaquinone, a catalytic intermediate.

The Ubiquinol Binding Site of Cytochrome bo3 from Escherichia coli Accommodates Menaquinone and Stabilizes a Functional Menasemiquinone

Cytochrome bo₃, 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₃ 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₃ 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₃ 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₃ does not abolish activity with soluble menaquinol analogues.

Location of the Substrate Binding Site of the Cytochrome bo3 Ubiquinol Oxidase from Escherichia coli

Cytochrome bo₃ is a respiratory proton-pumping oxygen reductase that is a member of the heme-copper superfamily that utilizes ubiquinol-8 (Q₈H₂) as a substrate. The current consensus model has Q₈H₂ oxidized at a low affinity site (Q_L), passing electrons to a tightly bound quinone cofactor at a high affinity site (Q_H site) that stabilizes the one-electron reduced ubisemiquinone, facilitating the transfer of electrons to the redox active metal centers where O₂ is reduced to water. The current work shows that the Q₈ bound to the Q_H site is more dynamic than previously thought. In addition, mutations of residues at the Q_H site that do not abolish activity have been re-examined and shown to have properties expected of mutations at the substrate binding site (Q_L): an increase in the K_M of the substrate ubiquinol-1 (up to 4-fold) and an increase in the apparent K_i of the inhibitor HQNO (up to 8-fold). The data suggest that there is only one binding site for ubiquinol in cyt bo₃ and that site corresponds to the Q_H site.

Searching for the low affinity ubiquinone binding site in cytochrome bo3 from Escherichia coli

The cytochrome bo₃ ubiquinol oxidase is one of three respiratory oxygen reductases in the aerobic respiratory chain of Escherichia coli. The generally accepted model of catalysis assumes that cyt bo₃ contains two distinct ubiquinol binding sites: (i) a low affinity (Q_L) site which is the traditional substrate binding site; and (ii) a high affinity (Q_H) site where a "permanently" bound quinone acts as a cofactor, taking two electrons from the substrate quinol and passing them one-by-one to the heme b component of the enzyme which, in turn, transfers them to the heme o₃/Cu_B active site. Whereas the residues at the Q_H site are well defined, the location of the Q_L site remains unknown. The published X-ray structure does not contain quinone, and substantial amounts of the protein are missing as well. A recent bioinformatics study by Bossis et al. [Biochem J. (2014) 461, 305-314] identified a sequence motif G¹⁶³EFX₃GWX₂Y¹⁷³ as the likely Q_L site in the family of related quinol oxidases. In the current work, this was tested by site-directed mutagenesis. The results show that these residues are not important for catalytic function and do not define the Q_L substrate binding site.

Q-Band Electron-Nuclear Double Resonance Reveals Out-of-Plane Hydrogen Bonds Stabilize an Anionic Ubisemiquinone in Cytochrome bo3 from Escherichia coli

The respiratory cytochrome bo₃ 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₂ 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 ¹H electron-nuclear double resonance (ENDOR) spectroscopy on fully deuterated cytochrome (cyt) bo₃ in a H₂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₃ 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.

Oxygen as Acceptor

Like most bacteria, Escherichia coli has a flexible and branched respiratory chain that enables the prokaryote to live under a variety of environmental conditions, from highly aerobic to completely anaerobic. In general, the bacterial respiratory chain is composed of dehydrogenases, a quinone pool, and reductases. Substrate-specific dehydrogenases transfer reducing equivalents from various donor substrates (NADH, succinate, glycerophosphate, formate, hydrogen, pyruvate, and lactate) to a quinone pool (menaquinone, ubiquinone, and dimethylmenoquinone). Then electrons from reduced quinones (quinols) are transferred by terminal reductases to different electron acceptors. Under aerobic growth conditions, the terminal electron acceptor is molecular oxygen. A transfer of electrons from quinol to O₂ is served by two major oxidoreductases (oxidases), cytochrome bo₃ encoded by cyoABCDE and cytochrome bd encoded by cydABX. Terminal oxidases of aerobic respiratory chains of bacteria, which use O₂ as the final electron acceptor, can oxidize one of two alternative electron donors, either cytochrome c or quinol. This review compares the effects of different inhibitors on the respiratory activities of cytochrome bo₃ and cytochrome bd in E. coli. It also presents a discussion on the genetics and the prosthetic groups of cytochrome bo₃ and cytochrome bd. The E. coli membrane contains three types of quinones that all have an octaprenyl side chain (C₄₀). It has been proposed that the bo₃ oxidase can have two ubiquinone-binding sites with different affinities. "WHAT'S NEW" IN THE REVISED ARTICLE: The revised article comprises additional information about subunit composition of cytochrome bd and its role in bacterial resistance to nitrosative and oxidative stresses. Also, we present the novel data on the electrogenic function of appBCX-encoded cytochrome bd-II, a second bd-type oxidase that had been thought not to contribute to generation of a proton motive force in E. coli, although its spectral properties closely resemble those of cydABX-encoded cytochrome bd.

Regulation of the primary quinone binding conformation by the H subunit in reaction centers from Rhodobacter sphaeroides

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.

Plasticity in the High Affinity Menaquinone Binding Site of the Cytochrome aa3-600 Menaquinol Oxidase from Bacillus subtilis

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.

The semiquinone at the Qi site of the bc1 complex explored using HYSCORE spectroscopy and specific isotopic labeling of ubiquinone in Rhodobacter sphaeroides via (13)C methionine and construction of a methionine auxotroph

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.

Conformational differences between the methoxy groups of QA and QB site ubisemiquinones in bacterial reaction centers: a key role for methoxy group orientation in modulating ubiquinone redox potential

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.

Interactions of intermediate semiquinone with surrounding protein residues at the Q(H) site of wild-type and D75H mutant cytochrome bo3 from Escherichia coli

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.

A rapid and robust method for selective isotope labeling of proteins

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.