The loneliness of the electrons in the bc1 complex

Correlated particle transport enables biological free energy transduction

Studies of biological transport frequently neglect the explicit statistical correlations among particle site occupancies (i.e., they use a mean-field approximation). Neglecting correlations sometimes captures biological function, even for out-of-equilibrium and interacting systems. We show that neglecting correlations fails to describe free energy transduction, mistakenly predicting an abundance of slippage and energy dissipation, even for networks that are near reversible and lack interactions among particle sites. Interestingly, linear charge transport chains are well described without including correlations, even for networks that are driven and include site-site interactions typical of biological electron transfer chains. We examine three specific bioenergetic networks: a linear electron transfer chain (as found in bacterial nanowires), a near-reversible electron bifurcation network (as in complex III of respiration and other recently discovered structures), and a redox-coupled proton pump (as in complex IV of respiration).

Kinetic Mathematical Modeling of Oxidative Phosphorylation in Cardiomyocyte Mitochondria

Oxidative phosphorylation (OXPHOS) is an oxygen-dependent process that consumes catabolized nutrients to produce adenosine triphosphate (ATP) to drive energy-dependent biological processes such as excitation-contraction coupling in cardiomyocytes. In addition to in vivo and in vitro experiments, in silico models are valuable for investigating the underlying mechanisms of OXPHOS and predicting its consequences in both physiological and pathological conditions. Here, we compare several prominent kinetic models of OXPHOS in cardiomyocytes. We examine how their mathematical expressions were derived, how their parameters were obtained, the conditions of their experimental counterparts, and the predictions they generated. We aim to explore the general landscape of energy production mechanisms in cardiomyocytes for future in silico models.

Mechanisms and mathematical modelling of ROS production by the mitochondrial electron transport chain

Reactive oxygen species (ROS) are recognised both as damaging molecules and intracellular signalling entities. In addition to its role in ATP generation, the mitochondrial electron transport chain (ETC) constitutes a relevant source of mitochondrial ROS, in particular during pathological conditions. Mitochondrial ROS homeostasis depends on species- and site-dependent ROS production, their bioreactivity, diffusion, and scavenging. However, our quantitative understanding of mitochondrial ROS homeostasis has thus far been hampered by technical limitations, including lack of truly site- and/or ROS-specific reporter molecules. In this context, the use of computational models is of great value to complement and interpret empirical data, as well as to predict variables that are difficult to assess experimentally. During the last decades, various mechanistic models of ETC-mediated ROS production have been developed. Although these often-complex models have generated novel insights, their parameterisation, analysis, and integration with other computational models is not straightforward. In contrast, phenomenological (sometimes termed "minimal") models use a relatively small set of equations to describe empirical relationship(s) between ROS-related and other parameters, and generally aim to explore system behaviour and generate hypotheses for experimental validation. In this review, we first discuss ETC-linked ROS homeostasis and introduce various detailed mechanistic models. Next, we present how bioenergetic parameters (e.g. NADH/NAD⁺ ratio, mitochondrial membrane potential) relate to site-specific ROS production within the ETC and how these relationships can be used to design minimal models of ROS homeostasis. Finally, we illustrate how minimal models have been applied to explore pathophysiological aspects of ROS.

In Silico Modeling of the Mitochondrial Pumping Complexes with Markov State Models

The mechanism of proton pumping by the mitochondrial electron transport chain complexes is still enigmatic after decades of research. Recently, there has been interest in in silico Markov state models to model the mitochondrial pumping complexes at the microscopic level, and this chapter describes the methods of constructing and simulating such models.

The modified Q-cycle: A look back at its development and forward to a functional model

On looking back at a lifetime of research, it is interesting to see, in the light of current progress, how things came to be, and to speculate on how things might be. I am delighted in the context of the Mitchell prize to have that excuse to present this necessarily personal view of developments in areas of my interests. I have focused on the Q-cycle and a few examples showing wider ramifications, since that had been the main interest of the lab in the 20 years since structures became available, - a watershed event in determining our molecular perspective. I have reviewed the evidence for our model for the mechanism of the first electron transfer of the bifurcated reaction at the Q_o-site, which I think is compelling. In reviewing progress in understanding the second electron transfer, I have revisited some controversies to justify important conclusions which appear, from the literature, not to have been taken seriously. I hope this does not come over as nitpicking. The conclusions are important to the final section in which I develop an internally consistent mechanism for turnovers of the complex leading to a state similar to that observed in recent rapid-mix/freeze-quench experiments, reported three years ago. The final model is necessarily speculative but is open to test.

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

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

The role of thermodynamic features on the functional activity of electron bifurcating enzymes

Electron bifurcation is a biological mechanism to drive a thermodynamically unfavorable redox reaction through direct coupling with an exergonic reaction. This process allows microorganisms to generate high energy reducing equivalents in order to sustain life and is often found in anaerobic metabolism, where the energy economy of the cell is poor. Recent work has revealed details of the redox energy landscapes for a variety of electron bifurcating enzymes, greatly expanding the understanding of how energy is transformed by this unique mechanism. Here we highlight the plasticity of these emerging landscapes, what is known regarding their mechanistic underpinnings, and provide a context for interpreting their biochemical activity within the physiological framework. We conclude with an outlook for propelling the field toward an integrative understanding of the impact of electron bifurcation.

The proton pumping mechanism of the bc1 complex

The bc₁ complex is a proton pump of the mitochondrial electron transport chain which transfers electrons from ubiquinol to cytochrome c. It operates via the modified Q cycle in which the two electrons from oxidation of ubiquinol at the Q_o center are bifurcated such that the first electron is passed to Cytc via an iron sulfur center and c₁ whereas the second electron is passed across the membrane by b_L and b_H to reduce ubiquinone at the Q_i center. Proton pumping occurs because oxidation of ubiquinol at the Q_o center releases protons to the P-side and reduction of ubiquinone at the Q_i center takes up protons from the N-side. However, the mechanisms which prevent the thermodynamically more favorable short circuit reactions and so ensure precise bifurcation and proton pumping are not known. Here we use statistical thermodynamics to show that reaction steps that originate from high energy states cannot support high flux even when they have large rate constants. We show how the chemistry of ubiquinol oxidation and the structure of the Q_o site can result in free energy profiles that naturally suppress flux through the short circuit pathways while allowing high rates of bifurcation. These predictions are confirmed through in-silico simulations using a Markov state model.

Universal free-energy landscape produces efficient and reversible electron bifurcation

For decades, it was unknown how electron-bifurcating systems in nature prevented energy-wasting short-circuiting reactions that have large driving forces, so synthetic electron-bifurcating molecular machines could not be designed and built. The underpinning free-energy landscapes for electron bifurcation were also enigmatic. We predict that a simple and universal free-energy landscape enables electron bifurcation, and we show that it enables high-efficiency bifurcation with limited short-circuiting (the EB scheme). The landscape relies on steep free-energy slopes in the two redox branches to insulate against short-circuiting using an electron occupancy blockade effect, without relying on nuanced changes in the microscopic rate constants for the short-circuiting reactions. The EB scheme thus unifies a body of observations on biological catalysis and energy conversion, and the scheme provides a blueprint to guide future campaigns to establish synthetic electron bifurcation machines.

Modelling mitochondrial ROS production by the respiratory chain

ROS (superoxide and oxygen peroxide in this paper) play a dual role as signalling molecules and strong oxidizing agents leading to oxidative stress. Their production mainly occurs in mitochondria although they may have other locations (such as NADPH oxidase in particular cell types). Mitochondrial ROS production depends in an interweaving way upon many factors such as the membrane potential, the cell type and the respiratory substrates. Moreover, it is experimentally difficult to quantitatively assess the contribution of each potential site in the respiratory chain. To overcome these difficulties, mathematical models have been developed with different degrees of complexity in order to analyse different physiological questions ranging from a simple reproduction/simulation of experimental results to a detailed model of the possible mechanisms leading to ROS production. Here, we analyse experimental results concerning ROS production including results still under discussion. We then critically review the three models of ROS production in the whole respiratory chain available in the literature and propose some direction for future modelling work.

Optimal efficiency of the Q-cycle mechanism around physiological temperatures from an open quantum systems approach

The Q-cycle mechanism entering the electron and proton transport chain in oxygenic photosynthesis is an example of how biological processes can be efficiently investigated with elementary microscopic models. Here we address the problem of energy transport across the cellular membrane from an open quantum system theoretical perspective. We model the cytochrome \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${b}_{6}\,f$$\end{document}b6f protein complex under cyclic electron flow conditions starting from a simplified kinetic model, which is hereby revisited in terms of a Markovian quantum master equation formulation and spin-boson Hamiltonian treatment. We apply this model to theoretically demonstrate an optimal thermodynamic efficiency of the Q-cycle around ambient and physiologically relevant temperature conditions. Furthermore, we determine the quantum yield of this complex biochemical process after setting the electrochemical potentials to values well established in the literature. The present work suggests that the theory of quantum open systems can successfully push forward our theoretical understanding of complex biological systems working close to the quantum/classical boundary.

Separation of photo-induced radical pair in cryptochrome to a functionally critical distance

Cryptochrome is a blue light receptor that acts as a sensor for the geomagnetic field and assists many animals in long-range navigation. The magnetoreceptor function arises from light-induced formation of a radical pair through electron transfer between a flavin cofactor (FAD) and a triad of tryptophan residues. Here, this electron transfer is investigated by quantum chemical and classical molecular dynamics calculations. The results reveal how sequential electron transfer, assisted by rearrangement of polar side groups in the cryptochrome interior, can yield a FAD-Trp radical pair state with the FAD and Trp partners separated beyond a critical distance. The large radical pair separation reached establishes cryptochrome's sensitivity to the geomagnetic field through weakening of distance-dependent exchange and dipole-dipole interactions. It is estimated that the key secondary electron transfer step can overcome in speed both recombination (electron back-transfer) and proton transfer involving the radical pair reached after primary electron transfer.

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

1. Recent results suggest that the major flux is carried by a monomeric function, not by an intermonomer electron flow. 2. The bifurcated reaction at the Qo-site involves sequential partial processes, - a rate limiting first electron transfer generating a semiquinone (SQ) intermediate, and a rapid second electron transfer in which the SQ is oxidized by the low potential chain. 3. The rate constant for the first step in a strongly endergonic, proton-first-then-electron mechanism, is given by a Marcus-Brønsted treatment in which a rapid electron transfer is convoluted with a weak occupancy of the proton configuration needed for electron transfer. 4. A rapid second electron transfer pulls the overall reaction over. Mutation of Glu-295 of cyt b shows it to be a key player. 5. In more crippled mutants, electron transfer is severely inhibited and the bell-shaped pH dependence of wildtype is replaced by a dependence on a single pK at ~8.5 favoring electron transfer. Loss of a pK ~6.5 is explained by a change in the rate limiting step from the first to the second electron transfer; the pK ~8.5 may reflect dissociation of QH. 6. A rate constant (<10(3)s(-1)) for oxidation of SQ in the distal domain by heme bL has been determined, which precludes mechanisms for normal flux in which SQ is constrained there. 7. Glu-295 catalyzes proton exit through H(+) transfer from QH, and rotational displacement to deliver the H(+) to exit channel(s). This opens a volume into which Q(-) can move closer to the heme to speed electron transfer. 8. A kinetic model accounts well for the observations, but leaves open the question of gating mechanisms. For the first step we suggest a molecular "escapement"; for the second a molecular ballet choreographed through coulombic interactions. This article is part of a Special Issue entitled: Respiratory complex III and related bc complexes.

A unifying kinetic framework for modeling oxidoreductase-catalyzed reactions

MOTIVATION

Oxidoreductases are a fundamental class of enzymes responsible for the catalysis of oxidation-reduction reactions, crucial in most bioenergetic metabolic pathways. From their common root in the ancient prebiotic environment, oxidoreductases have evolved into diverse and elaborate protein structures with specific kinetic properties and mechanisms adapted to their individual functional roles and environmental conditions. While accurate kinetic modeling of oxidoreductases is thus important, current models suffer from limitations to the steady-state domain, lack empirical validation or are too specialized to a single system or set of conditions.

RESULTS

To address these limitations, we introduce a novel unifying modeling framework for kinetic descriptions of oxidoreductases. The framework is based on a set of seven elementary reactions that (i) form the basis for 69 pairs of enzyme state transitions for encoding various specific microscopic intra-enzyme reaction networks (micro-models), and (ii) lead to various specific macroscopic steady-state kinetic equations (macro-models) via thermodynamic assumptions. Thus, a synergistic bridge between the micro and macro kinetics can be achieved, enabling us to extract unitary rate constants, simulate reaction variance and validate the micro-models using steady-state empirical data. To help facilitate the application of this framework, we make available RedoxMech: a Mathematica™ software package that automates the generation and customization of micro-models.

AVAILABILITY

The Mathematica™ source code for RedoxMech, the documentation and the experimental datasets are all available from: http://www.igb.uci.edu/tools/sb/metabolic-modeling.

CONTACT

pfbaldi@ics.uci.edu

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Mammalian complex I pumps 4 protons per 2 electrons at high and physiological proton motive force in living cells

Mitochondrial complex I couples electron transfer between matrix NADH and inner-membrane ubiquinone to the pumping of protons against a proton motive force. The accepted proton pumping stoichiometry was 4 protons per 2 electrons transferred (4H(+)/2e(-)) but it has been suggested that stoichiometry may be 3H(+)/2e(-) based on the identification of only 3 proton pumping units in the crystal structure and a revision of the previous experimental data. Measurement of proton pumping stoichiometry is challenging because, even in isolated mitochondria, it is difficult to measure the proton motive force while simultaneously measuring the redox potentials of the NADH/NAD(+) and ubiquinol/ubiquinone pools. Here we employ a new method to quantify the proton motive force in living cells from the redox poise of the bc(1) complex measured using multiwavelength cell spectroscopy and show that the correct stoichiometry for complex I is 4H(+)/2e(-) in mouse and human cells at high and physiological proton motive force.

Role of the -PEWY-glutamate in catalysis at the Q(o)-site of the Cyt bc(1) complex

We re-examine the pH dependence of partial processes of ubihydroquinone (QH(2)) turnover in Glu-295 mutants in Rhodobacter sphaeroides to clarify the mechanistic role. In more crippled mutants, the bell-shaped pH profile of wildtype was replaced by dependence on a single pK at ~8.5 favoring electron transfer. Loss of the pK at 6.5 reflects a change in the rate-limiting step from the first to the second electron transfer. Over the range of pH 6-8, no major pH dependence of formation of the initial reaction complex was seen, and the rates of bypass reactions were similar to the wildtype. Occupancy of the Q(o)-site by semiquinone (SQ) was similar in the wildtype and the Glu→Trp mutant. Since heme b(L) is initially oxidized in the latter, the bifurcated reaction can still occur, allowing estimation of an empirical rate constant <10(3)s(-1) for reduction of heme b(L) by SQ from the domain distal from heme b(L), a value 1000-fold smaller than that expected from distance. If the pK ~8.5 in mutant strains is due to deprotonation of the neutral semiquinone, with Q(•-) as electron donor to heme b(L), then in wildtype this low value would preclude mechanisms for normal flux in which semiquinone is constrained to this domain. A kinetic model in which Glu-295 catalyzes H(+) transfer from QH•, and delivery of the H(+) to exit channel(s) by rotational displacement, and facilitates rapid electron transfer from SQ to heme b(L) by allowing Q(•-) to move closer to the heme, accounts well for the observations.

Measurement of the mitochondrial membrane potential and pH gradient from the redox poise of the hemes of the bc1 complex

The redox potentials of the hemes of the mitochondrial bc(1) complex are dependent on the proton-motive force due to the energy transduction. This allows the membrane potential and pH gradient components to be calculated from the oxidation state of the hemes measured with multi-wavelength cell spectroscopy. Oxidation states were measured in living RAW 264.7 cells under varying electron flux and membrane potential obtained by a combination of oligomycin and titration with a proton ionophore. A stochastic model of bc(1) turnover was used to confirm that the membrane potential and redox potential of the ubiquinone pool could be measured from the redox poise of the b-hemes under physiological conditions assuming the redox couples are in equilibrium. The pH gradient was then calculated from the difference in redox potentials of cytochrome c and ubiquinone pool using the stochastic model to evaluate the ΔG of the bc(1) complex. The technique allows absolute quantification of the membrane potential, pH gradient, and proton-motive force without the need for genetic manipulation or exogenous compounds.

Spectral components of the α-band of cytochrome oxidase

Oxidative redox titrations of the mitochondrial cytochromes were performed in near-anoxic RAW 264.7 cells by inhibiting complex I. Cytochrome oxidation changes were measured with multi-wavelength spectroscopy and the ambient redox potential was calculated from the oxidation state of endogenous cytochrome c. Two spectral components were separated in the α-band range of cytochrome oxidase and they were identified as the difference spectrum of heme a when it has a high (a(H)) or low (a(L)) midpoint potential (E(m)) by comparing their occupancy during redox titrations carried out when the membrane potential (ΔΨ) was dissipated with a protonophore to that predicted by the neoclassical model of redox cooperativity. The difference spectrum of a(L) has a maximum at 605nm whereas the spectrum of a(H) has a maximum at 602nm. The ΔΨ-dependent shift in the E(m) of a(H) was too great to be accounted for by electron transfer from cytochrome c to heme a against ΔΨ but was consistent with a model in which a(H) is formed after proton uptake against ΔΨ suggesting that the spectral changes are the result of protonation. A stochastic simulation was implemented to model oxidation states, proton uptake and E(m) changes during redox titrations. The redox anti-cooperativity between heme a and heme a(3), and proton binding, could be simulated with a model where the pump proton interacted with heme a and the substrate proton interacted with heme a(3) with anti-cooperativity between proton binding sites, but not with a single proton binding site coupled to both hemes.

[The cytochrome bc1 complex in the mitochondrial respiratory chain functions according to the Q cycle hypothesis of Mitchell: the proof using a stochastic approach?]

The bc1 complex is a central complex in the mitochondrial respiratory chain. It links the electrons transfer from ubiquinol (or coenzyme Q) to cytochrome c and proton translocation across the inner mitochondrial membrane. It is widely agreed that the "Q-cycle mechanism" proposed by Mitchell correctly describes the bc1 complex working. It is based on an unexpected separation of the two electrons coming from the coenzyme Q bound at the Q0 site of the bc1 complex. Using the stochastic approach of Gillespie and the known spatial structure of bc1 complexes with the kinetic parameters described by Moser and Dutton we demonstrated the natural emergence of the Q-cycle mechanism and the quasi absence of short-circuits in the functional dimer of bc1 complex without the necessity to invoke any additional mechanism. This approach gives a framework which is well adapted to the modelling of all oxido-reduction reactions of the respiratory chain complexes, normal or mutant.

The Q cycle of cytochrome bc complexes: a structure perspective

Aspects of the crystal structures of the hetero-oligomeric cytochrome bc(1) and b(6)f ("bc") complexes relevant to their electron/proton transfer function and the associated redox reactions of the lipophilic quinones are discussed. Differences between the b(6)f and bc(1) complexes are emphasized. The cytochrome bc(1) and b(6)f dimeric complexes diverge in structure from a core of subunits that coordinate redox groups consisting of two bis-histidine coordinated hemes, a heme b(n) and b(p) on the electrochemically negative (n) and positive (p) sides of the complex, the high potential [2Fe-2S] cluster and c-type heme at the p-side aqueous interface and aqueous phase, respectively, and quinone/quinol binding sites on the n- and p-sides of the complex. The bc(1) and b(6)f complexes diverge in subunit composition and structure away from this core. b(6)f Also contains additional prosthetic groups including a c-type heme c(n) on the n-side, and a chlorophyll a and β-carotene. Common structure aspects; functions of the symmetric dimer. (I) Quinone exchange with the bilayer. An inter-monomer protein-free cavity of approximately 30Å along the membrane normal×25Å (central inter-monomer distance)×15Å (depth in the center), is common to both bc(1) and b(6)f complexes, providing a niche in which the lipophilic quinone/quinol (Q/QH(2)) can be exchanged with the membrane bilayer. (II) Electron transfer. The dimeric structure and the proximity of the two hemes b(p) on the electrochemically positive side of the complex in the two monomer units allow the possibility of two alternate routes of electron transfer across the complex from heme b(p) to b(n): intra-monomer and inter-monomer involving electron cross-over between the two hemes b(p). A structure-based summary of inter-heme distances in seven bc complexes, representing mitochondrial, chromatophore, cyanobacterial, and algal sources, indicates that, based on the distance parameter, the intra-monomer pathway would be favored kinetically. (III) Separation of quinone binding sites. A consequence of the dimer structure and the position of the Q/QH(2) binding sites is that the p-side QH(2) oxidation and n-side Q reduction sites are each well separated. Therefore, in the event of an overlap in residence time by QH(2) or Q molecules at the two oxidation or reduction sites, their spatial separation would result in minimal steric interference between extended Q or QH(2) isoprenoid chains. (IV) Trans-membrane QH(2)/Q transfer. (i) n/p-side QH(2)/Q transfer may be hindered by lipid acyl chains; (ii) the shorter less hindered inter-monomer pathway across the complex would not pass through the center of the cavity, as inferred from the n-side antimycin site on one monomer and the p-side stigmatellin site on the other residing on the same surface of the complex. (V) Narrow p-side portal for QH(2)/Q passage. The [2Fe-2S] cluster that serves as oxidant, and whose histidine ligand serves as a H(+) acceptor in the oxidation of QH(2), is connected to the inter-monomer cavity by a narrow extended portal, which is also occupied in the b(6)f complex by the 20 carbon phytyl chain of the bound chlorophyll.

The flitting of electrons in complex I: a stochastic approach

A stochastic approach based on the Gillespie algorithm is particularly well adapted to describe the time course of the redox reactions that occur inside the respiratory chain complexes because they involve the motion of single electrons between the individual unique redox centres of a given complex. We use this approach to describe the molecular functioning of the peripheral arm of complex I based on its known crystallographic structure and the rate constants of electron tunnelling derived from the Moser and Dutton phenomenological equations. There are several possible electrons pathways but we show that most of them take the route defined by the successive sites and redox centres: NADH+ site-FMN-N3-N1b-N4-N5-N6a-N6b-N2-Q site. However, the electrons do not go directly from NADH towards the ubiquinone molecule. They frequently jump back and forth between neighbouring redox centres with the result that the net flux of electrons through complex I (i.e. net number of electrons reducing a ubiquinone) is far smaller than the number of redox reactions which actually occur. While most of the redox centres are reduced in our simulations the degree of reduction can vary according to the individual midpoint potentials. The high turnover number observed in our simulation seems to indicate that, in the whole complex I, one or several slower step(s) follow(s) the redox reactions involved in the peripheral arm. It also appears that the residence time of FMNH* and SQ* (possible producers of ROS) is low (around 4% and between 1.6% and 5% respectively according to the values of the midpoint potentials). We did not find any evidence for a role of N7 which remains mainly reduced in our simulations. The role of N1a is complex and depends upon its midpoint potential. In all cases its presence slightly decreases the life time of the flavosemiquinone species. These simulations demonstrate the interest of this type of model which links the molecular physico-chemistry of the individual redox reactions to the more global level of the reaction, as is observed experimentally.

Visualizing changes in electron distribution in coupled chains of cytochrome bc(1) by modifying barrier for electron transfer between the FeS cluster and heme c(1)

Cytochrome c₁ of Rhodobacter (Rba.) species provides a series of mutants which change barriers for electron transfer through the cofactor chains of cytochrome bc₁ by modifying heme c₁ redox midpoint potential. Analysis of post-flash electron distribution in such systems can provide useful information about the contribution of individual reactions to the overall electron flow. In Rba. capsulatus, the non-functional low-potential forms of cytochrome c₁ which are devoid of the disulfide bond naturally present in this protein revert spontaneously by introducing a second-site suppression (mutation A181T) that brings the potential of heme c₁ back to the functionally high levels, yet maintains it some 100 mV lower from the native value. Here we report that the disulfide and the mutation A181T can coexist in one protein but the mutation exerts a dominant effect on the redox properties of heme c₁ and the potential remains at the same lower value as in the disulfide-free form. This establishes effective means to modify a barrier for electron transfer between the FeS cluster and heme c₁ without breaking disulfide. A comparison of the flash-induced electron transfers in native and mutated cytochrome bc₁ revealed significant differences in the post-flash equilibrium distribution of electrons only when the connection of the chains with the quinone pool was interrupted at the level of either of the catalytic sites by the use of specific inhibitors, antimycin or myxothiazol. In the non-inhibited system no such differences were observed. We explain the results using a kinetic model in which a shift in the equilibrium of one reaction influences the equilibrium of all remaining reactions in the cofactor chains. It follows a rather simple description in which the direction of electron flow through the coupled chains of cytochrome bc₁ exclusively depends on the rates of all reversible partial reactions, including the Q/QH₂ exchange rate to/from the catalytic sites.