The mechanism of ubihydroquinone oxidation at the Qo-site of the cytochrome bc1 complex

Mechanism of the Primary Charge Transfer Reaction in the Cytochrome bc1 Complex

The bc₁ complex is a critical enzyme for the ATP production in photosynthesis and cellular respiration. Its biochemical function relies on the so-called Q-cycle, which is well established and operates via quinol substrates that bind inside the protein complex. Despite decades of research, the quinol-protein interaction, which initiates the Q-cycle, has not yet been completely described. Furthermore, the initial charge transfer reactions of the Q-cycle lack a physical description. The present investigation utilizes classical molecular dynamics simulations in tandem with quantum density functional theory calculations, to provide a complete and consistent quantitative description of the primary events that occur within the bc₁ complex upon quinol binding. In particular, the electron and proton transfer reactions that trigger the Q-cycle in the bc₁ complex from Rhodobacter capsulatus are studied. The coupled nature of these charge transfer reactions was revealed by obtaining the transition energy path connecting configurations of the Q_o-site prior and after the transfers. The analysis of orbitals and partial charge distribution of the different states of the Q_o-site has further supported the conclusion. Finally, key structural elements of the bc₁ complex that trigger the charge transfer reactions were established, manifesting the importance of the environment in the process, which is furthermore evidenced by free energy calculations.

The cytochrome b Zn binding amino acid residue histidine 291 is essential for ubihydroquinone oxidation at the Qo site of bacterial cytochrome bc1

The ubiquinol:cytochrome (cyt) c oxidoreductase (or cyt bc₁) is an important membrane protein complex in photosynthetic and respiratory energy transduction. In bacteria such as Rhodobacter capsulatus it is constituted of three subunits: the iron-sulfur protein, cyt b and cyt c₁, which form two catalytic domains, the Q_o (hydroquinone (QH₂) oxidation) and Q_i (quinone (Q) reduction) sites. At the Q_o site, the pathways of bifurcated electron transfers emanating from QH₂ oxidation are known, but the associated proton release routes are not well defined. In energy transducing complexes, Zn²⁺ binding amino acid residues often correlate with proton uptake or release pathways. Earlier, using combined EXAFS and structural studies, we identified Zn coordinating residues of mitochondrial and bacterial cyt bc₁. In this work, using the genetically tractable bacterial cyt bc₁, we substituted each of the proposed Zn binding residues with non-protonatable side chains. Among these mutants, only the His291Leu substitution destroyed almost completely the Q_o site catalysis without perturbing significantly the redox properties of the cofactors or the assembly of the complex. In this mutant, which is unable to support photosynthetic growth, the bifurcated electron transfer reactions that result from QH₂ oxidation at the Q_o site, as well as the associated proton(s) release, were dramatically impaired. Based on these findings, on the putative role of His291 in liganding Zn, and on its solvent exposed and highly conserved position, we propose that His291 of cyt b is critical for proton release associated to QH₂ oxidation at the Q_o site of cyt bc₁.

Human Mitochondrial Cytochrome b Variants Studied in Yeast: Not All Are Silent Polymorphisms

Variations in mitochondrial DNA (mtDNA) cytochrome b (mt‐cyb) are frequently found within the healthy population, but also occur within a spectrum of mitochondrial and common diseases. mt‐cyb encodes the core subunit (MT‐CYB) of complex III, a central component of the oxidative phosphorylation system that drives cellular energy production and homeostasis. Despite significant efforts, most mt‐cyb variations identified are not matched with corresponding biochemical data, so their functional and pathogenic consequences in humans remain elusive. While human mtDNA is recalcitrant to genetic manipulation, it is possible to introduce human‐associated point mutations into yeast mtDNA. Using this system, we reveal direct links between human mt‐cyb variations in key catalytic domains of MT‐CYB and significant changes to complex III activity or drug sensitivity. Strikingly, m.15257G>A (p.Asp171Asn) increased the sensitivity of yeast to the antimalarial drug atovaquone, and m.14798T>C (p.Phe18Leu) enhanced the sensitivity of yeast to the antidepressant drug clomipramine. We demonstrate that while a small number of mt‐cyb variations had no functional effect, others have the capacity to alter complex III properties, suggesting they could play a wider role in human health and disease than previously thought. This compendium of new mt‐cyb‐biochemical relationships in yeast provides a resource for future investigations in humans.

New functional and biophysical insights into the mitochondrial Rieske iron-sulfur protein from genetic suppressor analysis in C. elegans

Several intragenic mutations suppress the C. elegans isp-1(qm150) allele of the mitochondrial Rieske iron-sulfur protein (ISP), a catalytic subunit of Complex III of the respiratory chain. These mutations were located in a helical region of the "tether" span of ISP-1, distant from the primary mutation in the extrinsic head, and suppressed all pleiotropic phenotypes associated with the qm150 allele. Analysis of these suppressors revealed control of electron transfer into Complex III through a "spring-loaded" mechanism involving a binding force for formation of enzyme-substrate complex, counter balanced by forces (a chemical "spring") favoring helix formation in the tether. The primary P→S mutation results in inhibition of electron flow into the Q-cycle by decreasing the binding force, and the tether mutations relieve this inhibition by weakening the "spring." In this commentary we discuss additional control features, and relate the primary inhibition to outcomes at the organismal level. In particular, the sensitivity to hyperoxia and the elevated reactive oxygen species (ROS) seen in isp-1(qm150), likely reflect over-reduction of the quinone pool, which is upstream of the inhibited site; at high O2, this would lead to increased ROS production through complex I. We speculate that alternative NADH:ubiquinone oxidoreductase activity in C. elegans from the worm apoptosis inducing factor (AIF) homolog (WAH-1) might also be involved, and that WAH-1 might have a "canary" function in detection of this adverse state (high O2/reduced pool), and a role in protection of the organism by transformation to AIF-like products, and apoptotic recycling of defective cells.

Interplay between the hinge region of iron sulphur protein and the Qo site in the bc1 complex - Analysis of Plasmodium-like mutations in the yeast enzyme

The respiratory chain bc1 complex is central to mitochondrial bioenergetics and the target of antiprotozoals. We characterized a modified yeast bc1 complex that more closely resemble Plasmodium falciparum enzyme. The mutant version was generated by replacing ten cytochrome b Qo site residues by P. falciparum equivalents. The Plasmodium-like changes caused a major dysfunction of the catalytic mechanism of the bc1 complex resulting in superoxide overproduction and respiratory growth defect. The defect was corrected by substitution of the conserved residue Y279 by a phenylalanine, or by mutations in or in the vicinity of the hinge domain of the iron-sulphur protein. It thus appears that side-reactions can be prevented by the substitution Y279F or the modification of the iron-sulphur protein hinge region. Interestingly, P. falciparum - and all the apicomplexan - contains an unusual hinge region. We replaced the yeast hinge region by the Plasmodium version and combined it with the Plasmodium-like version of the Qo site. This combination restored the respiratory growth competence. It could be suggested that, in the apicomplexan, the hinge region and the cytochrome b Qo site have co-evolved to maintain catalytic efficiency of the bc1 complex Qo site.

Mitochondrial Disease-related Mutation G167P in Cytochrome b of Rhodobacter capsulatus Cytochrome bc1 (S151P in Human) Affects the Equilibrium Distribution of [2Fe-2S] Cluster and Generation of Superoxide

Cytochrome bc₁ is one of the key enzymes of many bioenergetic systems. Its operation involves a large scale movement of a head domain of iron-sulfur protein (ISP-HD), which functionally connects the catalytic quinol oxidation Q_o site in cytochrome b with cytochrome c₁. The Q_o site under certain conditions can generate reactive oxygen species in the reaction scheme depending on the actual position of ISP-HD in respect to the Q_o site. Here, using a bacterial system, we show that mutation G167P in cytochrome b shifts the equilibrium distribution of ISP-HD toward positions remote from the Q_o site. This renders cytochrome bc₁ non-functional in vivo. This effect is remediated by addition of alanine insertions (1Ala and 2Ala) in the neck region of the ISP subunit. These insertions, which on their own shift the equilibrium distribution of ISP-HD in the opposite direction (i.e. toward the Q_o site), also act in this manner in the presence of G167P. Changes in the equilibrium distribution of ISP-HD in G167P lead to an increased propensity of cytochrome bc₁ to generate superoxide, which becomes evident when the concentration of quinone increases. This result corroborates the recently proposed model in which “semireverse” electron transfer back to the Q_o site, occurring when ISP-HD is remote from the site, favors reactive oxygen species production. G167P suggests possible molecular effects of S151P (corresponding in sequence to G167P) identified as a mitochondrial disease-related mutation in human cytochrome b. These effects may be valid for other human mutations that change the equilibrium distribution of ISP-HD in a manner similar to G167P.

Induction events and short-term regulation of electron transport in chloroplasts: an overview

Regulation of photosynthetic electron transport at different levels of structural and functional organization of photosynthetic apparatus provides efficient performance of oxygenic photosynthesis in plants. This review begins with a brief overview of the chloroplast electron transport chain. Then two noninvasive biophysical methods (measurements of slow induction of chlorophyll a fluorescence and EPR signals of oxidized P700 centers) are exemplified to illustrate the possibility of monitoring induction events in chloroplasts in vivo and in situ. Induction events in chloroplasts are considered and briefly discussed in the context of short-term mechanisms of the following regulatory processes: (i) pH-dependent control of the intersystem electron transport; (ii) the light-induced activation of the Calvin-Benson cycle; (iii) optimization of electron transport due to fitting alternative pathways of electron flow and partitioning light energy between photosystems I and II; and (iv) the light-induced remodeling of photosynthetic apparatus and thylakoid membranes.

Atomistic and Coarse Grain Topologies for the Cofactors Associated with the Photosystem II Core Complex

Electron transfers within and between protein complexes are core processes of the electron transport chains occurring in thylakoid (chloroplast), mitochondrial, and bacterial membranes. These electron transfers involve a number of cofactors. Here we describe the derivation of molecular mechanics parameters for the cofactors associated with the function of the photosystem II core complex: plastoquinone, plastoquinol, heme b, chlorophyll A, pheophytin, and β-carotene. Parameters were also obtained for ubiquinol and ubiquinone, related cofactors involved in the respiratory chain. Parameters were derived at both atomistic and coarse grain (CG) resolutions, compatible with the building blocks of the GROMOS united-atom and Martini CG force fields, respectively. Structural and thermodynamic properties of the cofactors were compared to experimental values when available. The topologies were further tested in molecular dynamics simulations of the cofactors in their physiological environment, e.g., either in a lipid membrane environment or in complex with the heme binding protein bacterioferritin.

Computational modeling analysis of mitochondrial superoxide production under varying substrate conditions and upon inhibition of different segments of the electron transport chain

A computational mechanistic model of superoxide (O2•-) formation in the mitochondrial electron transport chain (ETC) was developed to facilitate the quantitative analysis of factors controlling mitochondrial O2•- production and assist in the interpretation of experimental studies. The model takes into account all individual electron transfer reactions in Complexes I and III. The model accounts for multiple, often seemingly contradictory observations on the effects of ΔΨ and ΔpH, and for the effects of multiple substrate and inhibitor conditions, including differential effects of Complex III inhibitors antimycin A, myxothiazol and stigmatellin. Simulation results confirm that, in addition to O2•- formation in Complex III and at the flavin site of Complex I, the quinone binding site of Complex I is an additional superoxide generating site that accounts for experimental observations on O2•- production during reverse electron transfer. However, our simulation results predict that, when cytochrome c oxidase is inhibited during oxidation of succinate, ROS production at this site is eliminated and almost all superoxide in Complex I is generated by reduced FMN, even when the redox pressure for reverse electron transfer from succinate is strong. In addition, the model indicates that conflicting literature data on the kinetics of electron transfer in Complex III involving the iron-sulfur protein-cytochrome bL complex can be resolved in favor of a dissociation of the protein only after electron transfer to cytochrome bH. The model predictions can be helpful in understanding factors driving mitochondrial superoxide formation in intact cells and tissues.

Saccharomyces cerevisiae-based mutational analysis of the bc1 complex Qo site residue 279 to study the trade-off between atovaquone resistance and function

The bc1 complex is central to mitochondrial bioenergetics and the target of the antimalarial drug atovaquone that binds in the quinol oxidation (Qo) site of the complex. Structural analysis has shown that the Qo site residue Y279 (Y268 in Plasmodium falciparum) is key for atovaquone binding. Consequently, atovaquone resistance can be acquired by mutation of that residue. In addition to the probability of amino acid substitution, the level of atovaquone resistance and the loss of bc1 complex activity that are associated with the novel amino acid would restrict the nature of resistance-driven mutations occurring on atovaquone exposure in native parasite populations. Using the yeast model, we characterized the effect of all the amino acid replacements resulting from a single nucleotide substitution at codon 279: Y279C, Y279D, Y279F, Y279H, Y279N, and Y279S (Y279C, D, F, H, N, and S). Two residue changes that required a double nucleotide substitution, Y279A and W, were added to the series. We found that mutations Y279A, C, and S conferred high atovaquone resistance but decreased the catalytic activity. Y279F had wild-type enzymatic activity and sensitivity to atovaquone, while the other substitutions caused a dramatic respiratory defect. The results obtained with the yeast model were examined in regard to atomic structure and compared to the reported data on the evolution of acquired atovaquone resistance in P. falciparum.

The molecular evolution of the Qo motif

Quinol oxidation in the catalytic quinol oxidation site (Q_o site) of cytochrome (cyt) bc₁ complexes is the key step of the Q cycle mechanism, which laid the ground for Mitchell’s chemiosmotic theory of energy conversion. Bifurcated electron transfer upon quinol oxidation enables proton uptake and release on opposite membrane sides, thus generating a proton gradient that fuels ATP synthesis in cellular respiration and photosynthesis. The Q_o site architecture formed by cyt b and Rieske iron–sulfur protein (ISP) impedes harmful bypass reactions. Catalytic importance is assigned to four residues of cyt b formerly described as PEWY motif in the context of mitochondrial complexes, which we now denominate Q_o motif as comprehensive evolutionary sequence analysis of cyt b shows substantial natural variance of the motif with phylogenetically specific patterns. In particular, the Q_o motif is identified as PEWY in mitochondria, α- and ε-Proteobacteria, Aquificae, Chlorobi, Cyanobacteria, and chloroplasts. PDWY is present in Gram-positive bacteria, Deinococcus–Thermus and haloarchaea, and PVWY in β- and γ-Proteobacteria. PPWF only exists in Archaea. Distinct patterns for acidophilic organisms indicate environment-specific adaptations. Importantly, the presence of PDWY and PEWY is correlated with the redox potential of Rieske ISP and quinone species. We propose that during evolution from low to high potential electron-transfer systems in the emerging oxygenic atmosphere, cyt bc₁ complexes with PEWY as Q_o motif prevailed to efficiently use high potential ubiquinone as substrate, whereas cyt b with PDWY operate best with low potential Rieske ISP and menaquinone, with the latter being the likely composition of the ancestral cyt bc₁ complex.

The evolution of respiratory O2/NO reductases: an out-of-the-phylogenetic-box perspective

Complex life on our planet crucially depends on strong redox disequilibria afforded by the almost ubiquitous presence of highly oxidizing molecular oxygen. However, the history of O2-levels in the atmosphere is complex and prior to the Great Oxidation Event some 2.3 billion years ago, the amount of O2 in the biosphere is considered to have been extremely low as compared with present-day values. Therefore the evolutionary histories of life and of O2-levels are likely intricately intertwined. The obvious biological proxy for inferring the impact of changing O2-levels on life is the evolutionary history of the enzyme allowing organisms to tap into the redox power of molecular oxygen, i.e. the bioenergetic O2 reductases, alias the cytochrome and quinol oxidases. Consequently, molecular phylogenies reconstructed for this enzyme superfamily have been exploited over the last two decades in attempts to elucidate the interlocking between O2 levels in the environment and the evolution of respiratory bioenergetic processes. Although based on strictly identical datasets, these phylogenetic approaches have led to diametrically opposite scenarios with respect to the history of both the enzyme superfamily and molecular oxygen on the Earth. In an effort to overcome the deadlock of molecular phylogeny, we here review presently available structural, functional, palaeogeochemical and thermodynamic information pertinent to the evolution of the superfamily (which notably also encompasses the subfamily of nitric oxide reductases). The scenario which, in our eyes, most closely fits the ensemble of these non-phylogenetic data, sees the low O2-affinity SoxM- (or A-) type enzymes as the most recent evolutionary innovation and the high-affinity O2 reductases (SoxB or B and cbb3 or C) as arising independently from NO-reducing precursor enzymes.

Identification of ubiquinol binding motifs at the Qo-site of the cytochrome bc1 complex

Enzymes of the bc₁ complex family power the biosphere through their central role in respiration and photosynthesis. These enzymes couple the oxidation of quinol molecules by cytochrome c to the transfer of protons across the membrane, to generate a proton-motive force that drives ATP synthesis. Key for the function of the bc₁ complex is the initial redox process that involves a bifurcated electron transfer in which the two electrons from a quinol substrate are passed to different electron acceptors in the bc₁ complex. The electron transfer is coupled to proton transfer. The overall mechanism of quinol oxidation by the bc₁ complex is well enough characterized to allow exploration at the atomistic level, but details are still highly controversial. The controversy stems from the uncertain binding motifs of quinol at the so-called Q_o active site of the bc₁ complex. Here we employ a combination of classical all atom molecular dynamics and quantum chemical calculations to reveal the binding modes of quinol at the Q_o-site of the bc₁ complex from Rhodobacter capsulatus. The calculations suggest a novel configuration of amino acid residues responsible for quinol binding and support a mechanism for proton-coupled electron transfer from quinol to iron–sulfur cluster through a bridging hydrogen bond from histidine that stabilizes the reaction complex.

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.

The cytochrome b6f complex at the crossroad of photosynthetic electron transport pathways

Regulation of photosynthetic electron transport at the level of the cytochrome b6f complex provides efficient performance of the chloroplast electron transport chain (ETC). In this review, after brief overview of the structural organization of the chloroplast ETC, the consideration of the problem of electron transport control is focused on the plastoquinone (PQ) turnover and its interaction with the b6f complex. The data available show that the rates of plastoquinol (PQH2) formation in PSII and its diffusion to the b6f complex do not limit the overall rate of electron transfer between photosystem II (PSII) and photosystem I (PSI). Analysis of experimental and theoretical data demonstrates that the rate-limiting step in the intersystem chain of electron transport is determined by PQH2 oxidation at the Qo-site of the b6f complex, which is accompanied by the proton release into the thylakoid lumen. The acidification of the lumen causes deceleration of PQH2 oxidation, thus impeding the intersystem electron transport. Two other mechanisms of regulation of the intersystem electron transport have been considered: (i) "state transitions" associated with the light-induced redistribution of solar energy between PSI and PSII, and (ii) redistribution of electron fluxes between alternative pathways (noncyclic electron transport and cyclic electron flow around PSI).

pH-dependent regulation of electron transport and ATP synthesis in chloroplasts

This review is focused on pH-dependent mechanisms of regulation of photosynthetic electron transport and ATP synthesis in chloroplasts. The light-induced acidification of the thylakoid lumen is known to decelerate the plastoquinol oxidation by the cytochrome b 6 f complex, thus impeding the electron flow between photosystem II and photosystem I. Acidification of the lumen also triggers the dissipation of excess energy in the light-harvesting antenna of photosystem II, thereby protecting the photosynthetic apparatus against a solar stress. After brief description of structural and functional organization of the chloroplast electron transport chain, our attention is focused on the nature of the rate-limiting step of electron transfer between photosystem II and photosystem I. In the context of pH-dependent mechanism of photosynthetic control in chloroplasts, the mechanisms of plastoquinol oxidation by the cytochrome b 6 f complex have been considered. The light-induced alkalization of stroma is another factor of pH-dependent regulation of electron transport in chloroplasts. Alkalization of stroma induces activation of the Bassham-Benson-Calvin cycle reactions, thereby promoting efflux of electrons from photosystem I to NADP(+). The mechanisms of the light-induced activation of ATP synthase are briefly considered.

Plasticity of the quinone-binding site of the complex II homolog quinol:fumarate reductase

Respiratory processes often use quinone oxidoreduction to generate a transmembrane proton gradient, making the 2H(+)/2e(-) quinone chemistry important for ATP synthesis. There are a variety of quinones used as electron carriers between bioenergetic proteins, and some respiratory proteins can functionally interact with more than one quinone type. In the case of complex II homologs, which couple quinone chemistry to the interconversion of succinate and fumarate, the redox potentials of the biologically available ubiquinone and menaquinone aid in driving the chemical reaction in one direction. In the complex II homolog quinol:fumarate reductase, it has been demonstrated that menaquinol oxidation requires at least one proton shuttle, but many of the remaining mechanistic details of menaquinol oxidation are not fully understood, and little is known about ubiquinone reduction. In the current study, structural and computational studies suggest that the sequential removal of the two menaquinol protons may be accompanied by a rotation of the naphthoquinone ring to optimize the interaction with a second proton shuttling pathway. However, kinetic measurements of site-specific mutations of quinol:fumarate reductase variants show that ubiquinone reduction does not use the same pathway. Computational docking of ubiquinone followed by mutagenesis instead suggested redundant proton shuttles lining the ubiquinone-binding site or from direct transfer from solvent. These data show that the quinone-binding site provides an environment that allows multiple amino acid residues to participate in quinone oxidoreduction. This suggests that the quinone-binding site in complex II is inherently plastic and can robustly interact with different types of quinones.