Search tips
Search criteria 


Logo of ploscompComputational BiologyView this ArticleSubmit to PLoSGet E-mail AlertsContact UsPublic Library of Science (PLoS)
PLoS Comput Biol. Dec 2009; 5(12): e1000619.
Published online Dec 24, 2009. doi:  10.1371/journal.pcbi.1000619
PMCID: PMC2789320
Bistability of Mitochondrial Respiration Underlies Paradoxical Reactive Oxygen Species Generation Induced by Anoxia
Vitaly A. Selivanov,1,2,3 Tatyana V. Votyakova,4* Jennifer A. Zeak,4 Massimo Trucco,4 Josep Roca,2 and Marta Cascante1*
1Departament de Bioquimica i Biologia Molecular, Facultat de Biologia, Institut de Biomedicina at Universitat de Barcelona IBUB and IDIBAPS Hospital Clinic, Barcelona, Catalunya, Spain
2Hospital Clínic, IDIBAPS, CIBERES; Universitat de Barcelona, Barcelona, Catalunya, Spain
3A.N.Belozersky Institute of Physico-Chemical Biology, Moscow State University, Moscow, Russia
4Department of Pediatrics, The University of Pittsburgh School of Medicine and The Children's Hospital of Pittsburgh, Diabetes Institute, Pittsburgh, Pennsylvania, United States of America
G. R. Scott Budinger, Editor
Northwestern University, United States of America
* E-mail: tav2/at/ (TVV); martacascante/at/ (MC)
Performed the experiments: TVV JAZ. Analyzed the data: VAS MT JR MC. Wrote the paper: VAS TVV MC.
Received April 15, 2009; Accepted November 17, 2009.
Increased production of reactive oxygen species (ROS) in mitochondria underlies major systemic diseases, and this clinical problem stimulates a great scientific interest in the mechanism of ROS generation. However, the mechanism of hypoxia-induced change in ROS production is not fully understood. To mathematically analyze this mechanism in details, taking into consideration all the possible redox states formed in the process of electron transport, even for respiratory complex III, a system of hundreds of differential equations must be constructed. Aimed to facilitate such tasks, we developed a new methodology of modeling, which resides in the automated construction of large sets of differential equations. The detailed modeling of electron transport in mitochondria allowed for the identification of two steady state modes of operation (bistability) of respiratory complex III at the same microenvironmental conditions. Various perturbations could induce the transition of respiratory chain from one steady state to another. While normally complex III is in a low ROS producing mode, temporal anoxia could switch it to a high ROS producing state, which persists after the return to normal oxygen supply. This prediction, which we qualitatively validated experimentally, explains the mechanism of anoxia-induced cell damage. Recognition of bistability of complex III operation may enable novel therapeutic strategies for oxidative stress and our method of modeling could be widely used in systems biology studies.
The levels of reactive oxygen species (ROS) that are generated as a side product of mitochondrial respiratory electron transport largely define the extent of oxidative stress in living cells. Free radicals formed in electron transport, such as ubisemiquinone, could pass their non-paired electron directly to oxygen, thus producing superoxide radical that gives rise to a variety of ROS. It is well known in clinical practice that upon recommencing oxygen supply after anoxia a tissue produces much more ROS than before the anoxia, and the state of high ROS production is stable. The mechanism of switching from low to high ROS production by temporal anoxia was unknown, in part because of the lack of detailed mathematical description of hundreds of redox states of respiratory complexes, which are formed in the process of electron transport. A new methodology of automated construction of large systems of differential equations allowed us to describe the system in detail and predicts that the mechanism of paradoxical effect of anoxia-reoxygenation could be defined by the properties of complex III of mitochondrial respiratory chain. Our experiments confirmed that the effect of hypoxia-reoxygenation is confined by intramitochondrial processes since it is observed in isolated mitochondria.
The pathologic consequences of anoxia–reoxygenation, including the oxidative stress associated with increased production of reactive oxygen species (ROS) in mitochondria, form the basis of major diseases, including heart disease, age-related degenerative conditions and ischemic syndrome in reperfusion [1]. The use of novel antioxidants, which addresses the consequences of oxidative stress, has proven to be effective in organ preservation [2], but there is no doubt that a better understanding of the causes of elevated ROS production during the anoxia/reoxygenation would help to introduce novel strategies addressing the primary events of this clinical phenomenon.
Paradoxically, ROS production increases during severe hypoxia despite of decrease of oxygen concentration; this ROS increase acts as a metabolic signal for cell adaptation to oxygen deficiency [3]. Moreover, when cells that have been exposed to anoxia are returned to their normal oxygen supply, the rate of ROS production, instead of returning to the low baseline, greatly increases, often leading to cell death. This phenomenon of reperfusion injury after ischemia has been well known [4],[5], but its mechanism remains unclear and it is analyzed here.
Mitochondrial respiratory electron transport, which is schematically shown in Figure S1, is generally accepted as a process related with ROS production in living cells. In particular, o-site of quinol oxidation in complex III (Qo) is one of the most frequently considered sites of superoxid anion generation [3],[6],[7], and the considered here mechanism for the stimulation of ROS production by anoxia is related to this site. The proposed mechanism is based on our analysis of the well-known Q-cycle (ubiquinone (Q) oxidation/reduction) mechanism of electron and proton transport performed by complex III in the mitochondrial respiratory chain [8][13], as schematically shown in Figure 1. The main components of complex III are cytochrome b, containing two hemes characterized by low (bL) and high (bH) midpoint potentials, Rieske protein containing an iron-sulfur redox center (FeS), and cytochrome c1. The FeS center of the Rieske protein accepts one electron from a bound ubiquinol (QH2) producing a highly reactive anion radical of ubiquinone (Q) known also as semiquinone radical (further referred as SQ) and releasing two H+ to the cytosolic side. The electron accepted by the FeS center then is delivered to c1 and passes further downstream in the electron transport. The semiquinone radical is normally transformed into ubiquinone as it delivers its unpaired electron to cytochrome bL, which then passes it to bH. However, there is a probability that the semiquinone radical delivers its unpaired electron directly to oxygen, producing a superoxide radical [14][16]. Semiquinone oxidation by any of these two mechanisms transforms it into ubiquinone, which then dissociates from complex III, binds at the matrix side and receives two electrons from cytochrome bh that are derived from oxidation of two QH2 molecules as described above. This process subsequently generates semiquinone and QH2 again, taking protons from the matrix. The dissociation of the newly produced QH2 and its subsequent binding at the cytosolic side starts the cycle again.
Figure 1
Figure 1
Scheme of the redox reactions in complex III of mitochondrial respiratory chain.
A great deal of problems in studying ROS production in respiratory chain induced by anoxia arises because of the lack of tools for systematic theoretical analysis of this complicated system. Various combinations of oxidized and reduced states of the electron transporters create 400 possible different redox states for complex III, which are formed in the process of electron transport and which should be taken into account in an analysis of electron transport and the related ROS production. The high levels of complexity has, until now, precluded understanding of complex III functionality. To overcome this problem, we have developed a special software tool that takes into account all possible redox states of complex III. The software automatically constructs a system of differential equations describing their evolution based on several rules defined by the types of reactions between the electron transporters, as outlined in Methods and in more detail in Text S1. Our current approach to simulating mitochondrial respiration has inherited the main principles of automated construction of large equation systems used for stable isotope tracer data analysis [17][19].
Using this tool for the simulations of evolution of all redox states of complex III, we found that complex III has a property of bistability: it can persist in two different steady states at the same set of parameters; evolution to one or another steady state is defined only by the initial state of the system. One such state is characterized by a high ROS production rate. The system can be switched to this latter state either by an increase in succinate supply or by a decrease of oxygen availability, and can persist in it after a return to the initial conditions. This behavior explains the mechanics of paradoxical increase in ROS generation induced by anoxia and its further increase after return to a normal oxygen supply. Experimental data presented here qualitatively confirmed such a bistable behavior.
The phenomenon of bistability in complex III operation
The model constructed, as described in Methods, is based on the generally accepted Q-cycle mechanism of complex III operation and predicts that it can evolve to a one of two different steady states. The direction of complex III evolution depends on the initial state of the system. If initially it is in highly reduced state, characterized by high levels of ubiquinol and semiquinone, it remains in this state, as shows thick gray lines in Figures 2A and 2B. However, if the initial levels of reduction is lower, the system evolves to another steady state, the same for a variety of initial states, as Figures 2A and 2B show. Thus, the Q-cycle mechanism defines bistability of complex III operation: with the same parameters it could function in a mode characterized by either low or high SQ levels (low or high ROS producing states respectively).
Figure 2
Figure 2
The prediction of bistability in complex III operation.
Two branches of steady states in the continuum of succinate concentrations
Various factors could trigger complex III from one mode to another, specifically, the availability of succinate could be a triggering factor. Figure 3A shows the time course of the transition from low to high ROS producing steady state, induced by an increase in succinate concentration (eq (17)). The capacity for ROS production is reflected in the concentration of free radicals, which can directly interact with oxygen, such as the semiquinone radical (SQ) bound to the Qo site of complex III [14],[20]. An increase in succinate supply, and respective reduction of free ubiquinone (Q) to ubiquinol (QH2), switches the system from low ROS producing state, characterized by low levels of SQ, to the high ROS producing state, characterized by high levels of SQ. In accordance with mass conservation, reduction of Q to QH2 results in the deficiency of Q; this decreases electron flow and promotes reduction of cytochrome b. The QH2 bound at Qo site can freely deliver its first electron to the cytochrome c1, but cannot deliver the second electron to cytochrome bL because latter is already reduced. As a result, highly active SQ radicals are accumulated.
Figure 3
Figure 3
Bi-stability of the reaction system catalyzed by complex III.
The continuum of steady state SQ levels and electron flow dependent on succinate supply is shown in Figures 3B and 3C. If the system is initially in a low ROS producing state, it switches to the high ROS production when the succinate concentration reaches a certain threshold value as Figure 3B shows. However, if the succinate concentration decreases back when the system initially is in high ROS producing state, it does not follow the same pattern, as Figure 3C shows. This example demonstrates hysteresis in complex III behavior.
The value of transmembrane potential is essential for such hysteretic behavior. As Figure S2 shows, the region of bistability is clearly distinguishable at transmembrane potential of ~200 mV, which corresponds to the state 4 of mitochondrial respiration. The fall down of transmembrane potential to below 150 mV (as could be when the addition of ADP switches mitochondria to the state 3 of respiration), switches complex III from high to low ROS producing state. At high transmembrane potential the region of bistability persists over a large variation of model parameters as the sensitivity analysis shown in Text S1 and Figures S3 and S4 indicates. A switch from one steady state to another one essentially redistributes the fractions of various redox states of the complex as Figure S5 shows.
The electron flux in high ROS producing steady state, being restricted by the deficiency of Q, is low; maintaining it requires low succinate supply and this steady state persists even if substrate supply decreases until it falls down below the minimal threshold, as displayed in Figure 3C. The reduction of molecular oxygen by semiquinone radicals transforms them into Q, producing an acceptor able to take electrons from bh, and thus to activate the Q-cycle. Thus, high ROS production rate, in a way, is a means to return back to low ROS production. Even a low rate of such an electron leak to oxygen helps the system to revert to low ROS production (see sensitivity analysis presented in Text S1). Normally direct transfer to oxygen insignificantly contributes to the total electron flow, but in the case of extremely high complex II activity, if it reduces all available Q, the electron flow to complex IV would equalize the ROS production rate.
Experimental confirmation of bistability
The predicted phenomenon of bistability explored in Figure 3 could be observed experimentally, as Figure 4A shows. The isolated rat brain mitochondria incubated with succinate (state 4 of respiration) are in a high ROS producing state (blue trace “ros”), however, temporal presence of ADP, which is rapidly transformed into ATP, switches them to a low ROS producing mode (red trace “ros,ADP”). When ADP is present, ATP synthesis lowers membrane potential (shown in reverse direction as measured by quenching of fluorescence, red trace “mp,ADP”). However, after complete transformation of the added ADP into ATP, mitochondria again come back to state 4 of respiration, the membrane potential increases to the same levels as that without ADP (blue trace “mp”), but ROS production remains to be much lower (red trace “ros,ADP”). ADP only switched respiratory chain from high to low ROS production state, which could be maintained under the same conditions of incubation. The addition of ATP to the medium had no effect on ROS production (blue trace “ros”). In this case ROS production was the same as in the absence of ATP (Figure S6).
Figure 4
Figure 4
Experimental observations of bistability in mitochondrial electron transport and their simulation in the model.
Figure 4B shows two branches of steady state ROS production rates measured before addition of ADP (stars) and after its consumption (squares) (as Figure 4A explains). Solid lines generated by computer simulation show that our model has reproduced these experimental data after implementing the Michaelis' dependence of complex II activity on succinate concentration (eq (17)). Thus, isolated mitochondria show the predicted bistable behavior (Figure 3) and the model simulates the particular experimental data (Figure 4B). In these experiments 5 mM of pyruvate were present at variable concentrations of succinate to ensure that overall substrate supply is not a limiting factor. In parallel experiments it was shown that at different succinate concentrations membrane potential was invariable (not shown). It should be noted that the parameters used as characteristics of succinate oxidation (Vm, Km) in complex II, in fact characterize several processes including succinate transport. In our in vitro experiments the amount of succinate must be sufficiently high to maintain the rate of its transport needed for ubiquinone reduction. The intramitochondrial succinate concentration, which in vivo is produced in TCA cycle and which is sufficient to trigger high ROS production, could be much less than that necessary to be added externally. The actual parameters, which govern the oxidation of succinate produced in the matrix in vivo, could be different, and since phenomenologically the mechanism admits the forward and reverse switch between the steady states (Figure 3), it could be expected under some conditions, such as ischemia, which could induce several fold increase in intracellular succinate content and, as a consequence, increase in FADH2 [21]. Another factor, which substitutes succinate in triggering ROS production in vivo, is reduction of ubiquinone by complex I, whose substrate, NADH, is accumulated during anoxia [21].
Anoxia-induced switch to high ROS producing steady state
The lack of oxygen itself can induce a switch to high ROS production if initially respiratory chain is in a low ROS production mode as Figure 4A shows. For this simulation the electron flux from complex III to complex IV was described as Michaelis' function of oxygen concentration (eq (18)). Parameters were chosen so that oxygen availability limits the flux when oxygen concentration falls to about 0.1% O2, (in accordance with [13],[22]). The decrease of oxygen availability at the excess of succinate results in a reduction of ubiquinone and cytochrome b, and, finally, an accumulation of SQ radicals. As Figure 5A shows, the accumulation of SQ radicals proceeds in conjunction with limiting the electron flow by the lack of oxygen. Apparently, the physiological response to mild hypoxia, which could be observed at 5% O2 [23], is related with another mechanisms.
Figure 5
Figure 5
Theoretical predictions and experimental observations of the bistability in mitochondrial electron transport triggered by anoxia.
If after passing through nearly anoxic conditions oxygen availability increases, the system remains in a high ROS producing state. Thus, the functional organization of electron transport in complex III, which has the property of bistability, explains the apparent paradox of increased ROS production induced by anoxia and re-oxygenation.
Figure 5B shows that anoxia induces a switch to a high ROS producing mode in isolated brain mitochondria, which were switched to a low ROS producing mode by temporal presence of ADP. Experiment was started with an open cuvette and normal oxygen supply by an addition of mitochondria to the incubation medium containing 1 mM ADP. Low value of membrane potential is indicative of mitochondria being initially in state 3. Upon conversion of all present ADP into ATP, mitochondria transited into state 4; an increased membrane potential is indicative of that transition. During the first stage of experiment mitochondria released low, but measurable amounts of ROS. An anoxic state was achieved by bubbling of the cuvette with N2 for 2–3 min and then closing it with a Teflon lid. Slow depolarization of mitochondria and their persistence in the depolarized state for the duration of closure indicated an achievement of anoxic state. Notably, the ROS signal stayed flat during that stage. At the third stage of the experiment, the resumption of oxygen supply was achieved by opening the cuvette and bubbling it with air. At this stage mitochondria regained membrane potential up to the initial level thus showing their metabolic integrity. Significantly higher than at initial normoxic stage amounts of ROS were released upon the resumption of oxygen supply.
Although the damaging effect of reoxygenation after anoxia is well known in general, presented here experimental data for the first time directly demonstrate anoxia-induced increase in ROS at the levels of isolated mitochondria.
While at constant oxygen concentration the levels of SQ at Qo site could be considered as a relative rate of ROS production, the change of oxygen concentration must also affect it. Figure 6 shows the prediction of ROS production rate when the oxygen concentration changes as Figure 5A shows. For the calculations of ROS production eq. (16) was used with various KROS. Figure 6 indicates that the increase in ROS production could be detectable despite the decrease of oxygen concentration, if KROS is in the range of Km for respiration. Increase of KROS results in the decreases ROS production before it starts to affect the respiration rate that induces increase of SQ. In this case the change of SQ would not be well detectable. Thus, if experimentally ROS production was found to increase under anoxic conditions [3], it means that KROS must be in the range of Km for respiration or lower.
Figure 6
Figure 6
Theoretical predictions of ROS production under hypoxic conditions.
The algorithms for automated construction of large systems of differential equations, which are presented in Methods, allowed us to simulate in detail the Q-cycle mechanism of complex III operation and to find that the property of bistability is inherited from this mechanism. Due to this property it can operate in a low or high ROS producing mode at the same values of parameters, which explains the difference in ROS generation by respiratory chain before and after an episode of anoxia. The presented experiments qualitatively confirm the existence of two different ROS producing steady states in isolated mitochondria at the same microenvironmental conditions and a possibility to switch between them by the temporal presence of ADP and/or episodes of anoxia.
The following biochemical events were found to underlie the transition from low to high ROS producing state: (i) ubiquinone pool is reduced to ubiquinol, (ii) lack of oxidized form (ubiquinone) restricts cytochrome bH oxidation and results in bH and bL reduction, which, in turn, (iii) restricts the oxidation of semiquinone radicals, thus increasing their concentration. If substrate concentration is then decreased below a threshold, the remaining electron flow is sufficient to maintain this high ROS producing state and therefore low and high ROS producing states could exist at the same set of parameters.
The increase of ROS production in tissues caused by ischemia-reperfusion is a well known clinical phenomenon; it is of great importance in various systemic diseases and organ transplantation. However, it was not known, which system is responsible for the switch from normally low ROS producing state to high ROS production, and how such a switch proceeds. We predicted and experimentally confirmed that complex III could respond to anoxia-reoxygenation in the same manner as it was observed at tissue levels. This indicates that the revealed mechanism could constitute the main part of tissue response to anoxia-reoxygenation. It does not exclude, though, a possible contribution of other parts of respiratory chain and/or other parts of metabolic and signaling pathways as modifiers of the analyzed phenomenon or even as other individual sources of ROS; rather, our findings open a new direction in investigation of the molecular mechanisms of bistable behavior that underlies the signaling effect of anoxia and oxidative stress caused by reoxygenation. In living cells ischemia combines several factors affecting ROS production, such as reduction of electron transport chain, anoxia, converting ATP into ADP and AMP, and deamination of the latter resulting in decrease of cellular adenilate pool [24]. Various combinations of these factors could create various types of behavior. Understanding their underlying mechanisms will improve the clinical care of diseases related to anoxia and also the quality of organ transplantation.
The proposed method for kinetic description of multienzyme complexes inherited from the established method of automated equation construction, which we developed for a different area of investigation, namely the analysis of stable isotopic isomers distribution in the intermediates of central carbohydrate metabolism [17][19]. Thus, it is not restricted by an area of mitochondrial respiration and could be further modified for applications in different areas, for instant for the analysis of metabolic signalling, which requires a description of multiple states of receptor molecules. Using this method will allow further expanding the model, including in it the other parts of electron transport chain and central metabolism, which provides substrates for mitochondrial respiration, ion transport that defines transmembrane potential, etc. Thus, the proposed method of differential equations construction opens a new direction for most detailed mathematical analysis of complex bioenergetic systems.
The model simulates the reactions performed by complex III, which overall expression is:
A mathematical equation, expression, or formula.
 Object name is pcbi.1000619.e001.jpg
The two electrons taken from QH2 then are transferred to cytochrome c and further to complex IV and oxygen in conjunction with translocation of four more protons. Although the presented model does not consider the details of this latter process, it accounts for it stoichiometry.
In the course of complex III reaction its core (consisted of cytochrome b with its two hemes, bH and bL, cytochrome c1, and Fe-S containing center of Rieske protein) bind and dissosiate ubiquinone/ubiquinol either in the matrix (Qi) or cytosolic (Qo) side of the inner mitochondrial membrane. This binding of quinones produce four different species of the complex:
A mathematical equation, expression, or formula.
 Object name is pcbi.1000619.e002.jpg
A mathematical equation, expression, or formula.
 Object name is pcbi.1000619.e003.jpg
A mathematical equation, expression, or formula.
 Object name is pcbi.1000619.e004.jpg
A mathematical equation, expression, or formula.
 Object name is pcbi.1000619.e005.jpg
The electron transporters consisting the species (1)–(4) could be either in reduced or oxidized state. Various combinations of redox states of transporters produce many forms of each specie. The concentrations of such forms, which are transformed one into another during the complex operation, are the variables of the model. All the concentrations calculated for different redox forms of the same specie are stored in a specific array. To facilitate the references to variables inside the program, the two redox states of each transporter represented either by the digit “0” for oxidized state or by “1” for the reduced one. For instance, the binary “0011” (or decimal “3”) with respect to the array of core (1) refers to the concentration of a specie with oxidized bH and bL and reduced c1 and FeS. On the other hand, these digits together form a binary number, which signifies that the concentration of this redox form occupies third position (starting from 0) in this array. Thus, the program uses such integer number for two purposes: to designate a specific redox form and in the same time to refer to the place of its concentration in the array reserved for the concentrations of all possible redox forms.
The total number of redox forms for a specie is defined by the number of transporters constituting it, being 2n, where n is the number of transporters. In this way, core (1) has 16 different redox states varying from 0000 to 1111; species (2) and (3) have 64 states each, varying from 000000 to 111111, and, respectively, the specie (4) has 256 redox states. The program constructs differential equations for each of these redox states based on the algorithms that reflect the redox reactions between transporters in accordance with the known biochemical mechanism of electron transport in complex III. As is schematically shown in Figure 2, the model takes into account 12 types of reactions, nine of them constitute the proper Q-cycle and three others are just related with complex III activity: superoxid anion formation (10), QH2 supply by succinate dehydrogenase, and cytocrome c1 oxidation by cytochrome c.
The algorithms for automated calculation of reaction rates for each redox form participated in thece 12 types of reactions is described in details next.
Reaction 0 transports the first electron of QH2 bound on the positive (p-, or cytosolic) side to Fe3+ of Rieske protein and releases two protons to the intermembrane space.
A mathematical equation, expression, or formula.
 Object name is pcbi.1000619.e006.jpg
The program simulates this reaction for all redox states of complex III capable to perform it, i.e. for the forms with attached QH2 at Qo site and oxidized Fe-S center of Rieske protein. In binary terminology described above for species bH-bL-c1-FeS-Q-Q the reaction (5) could be written as follows:
A mathematical equation, expression, or formula.
 Object name is pcbi.1000619.e007.jpg
here x designate any redox state of a respective transporter, 0 or 1. The expressions for reaction rates according to the mass action law are:
A mathematical equation, expression, or formula.
 Object name is pcbi.1000619.e008.jpg
The program passes through the array for species (2), checks if the position of variable satisfy the requirements for substrate given in (5a) (that three last digits are 011), calculates the position of respective product in the array (three last digits are 101 and the digits designated as x are the same as for chosen substrate), then calculates fluxes (6) and adds them to the derivatives for the substrate and product. For the derivatives the program also has an array, where positions are organized in the same way as for concentrations.
For the other species capable to perform this reaction, Qi-Qi-bH-bL-c1-FeS-Qo-Qo, it could be written as
A mathematical equation, expression, or formula.
 Object name is pcbi.1000619.e009.jpg
and reaction rates expressed similar to (6).
The model takes into account the ratio of forward and reverse rate constants, which could be defined from the difference of midpoint redox potentials (Em):
A mathematical equation, expression, or formula.
 Object name is pcbi.1000619.e010.jpg
where n is the number of electrons transported, F = 96500 c/mol is the Faraday constant, R = 8.3 J/(mol×K) is the gas constant, T = 298 K is temperature. The derivation of this equation is based on the basic thermodynamic principles as is described in details in Text S1. The reported values for Em at pH 7 vary as much as 312 mV, [25] or 280 mV, [26] for Fe3+/Fe2+ and 200 mV, [27] or 300 mV, [28] for Q-/QH2; this variability is reflected in the variability of value expressed by (6a) from 0.6 to 18.
Reaction 1 transports the electron received by Fe2+ of Rieske protein from QH2 further to c1:
A mathematical equation, expression, or formula.
 Object name is pcbi.1000619.e011.jpg
All the species (1)–(4) with reduced FeS center of Rieske protein can perform this reaction:
A mathematical equation, expression, or formula.
 Object name is pcbi.1000619.e012.jpg
Reaction rates for core species:
A mathematical equation, expression, or formula.
 Object name is pcbi.1000619.e013.jpg
The ratio of forward and reverse rate constants is defined by an equation similar to (6a):
A mathematical equation, expression, or formula.
 Object name is pcbi.1000619.e014.jpg
Midpoint potentials, for Fe3+/Fe2+ of Riske protein is 312 mV [25] and for c1 it is 341 mV [29]. This difference (ΔEm = 29 mV) results in three-time difference in forward and reverse rate constants.
Reaction 2 transports the second electron of Qo to bL:
A mathematical equation, expression, or formula.
 Object name is pcbi.1000619.e015.jpg
Only the species with ocupied Qo site perform this reaction:
A mathematical equation, expression, or formula.
 Object name is pcbi.1000619.e016.jpg
Reaction rates for bH-bL-c1-FeS-Qo-Qo:
A mathematical equation, expression, or formula.
 Object name is pcbi.1000619.e017.jpg
The reported Em (−60 mV for bL [30] and −140 mV for Q/Q [28]) giving ΔEm = 80 mV according to the equation similar to (7a) define the ratio of forward to reverse rate constant as 25.
Reaction 3 transports electron from bL to bH
A mathematical equation, expression, or formula.
 Object name is pcbi.1000619.e018.jpg
All species with reduced bL and oxidized bH perform this reaction:
A mathematical equation, expression, or formula.
 Object name is pcbi.1000619.e019.jpg
Reaction rates for core species:
A mathematical equation, expression, or formula.
 Object name is pcbi.1000619.e020.jpg
Since in this reaction electron moves in the direction from positive outer side to negative matrix side of the membrane, the rate constants depend of electric transmembrane potential [31].
A mathematical equation, expression, or formula.
 Object name is pcbi.1000619.e021.jpg
where ΔΨ is transmembrane potential.
Here the ratio of forward and reverse k0 is defined by midpoint redox potential difference (similar to (7a)). The reported values are −58 mV for bL [30] and 61 mV for bH [30]. This difference (ΔEm = 119 mV) defines the forward rate constant to be two orders of magnitude higher than the reverse rate constant.
Reaction 4 transports first electron from bH to Q on the n-side:
A mathematical equation, expression, or formula.
 Object name is pcbi.1000619.e022.jpg
This reaction is performed by the species with ubiquinone bound at Qi site:
A mathematical equation, expression, or formula.
 Object name is pcbi.1000619.e023.jpg
Reaction rates:
A mathematical equation, expression, or formula.
 Object name is pcbi.1000619.e024.jpg
For bH we accepted Em = 61 mV, reported in [30]; for Q/Q, the reported Em vary from 90 mV [30] to 45 mV [32]. To start, we accepted the value which is close to Em for bH, in this case the forward and reverse rate constants are equal.
Reaction 5 transports second electron from bH to Q on the n-side:
A mathematical equation, expression, or formula.
 Object name is pcbi.1000619.e025.jpg
This reaction can be performed by species that contain semiquinone bound on Qi site:
A mathematical equation, expression, or formula.
 Object name is pcbi.1000619.e026.jpg
Reaction rates:
A mathematical equation, expression, or formula.
 Object name is pcbi.1000619.e027.jpg
For bH Em = 61 mV [30]; for Q/QH2 data in literature vary from 16.5 mV [30] to 150 mV [32]. To start, for this transition we also accepted the value close to Em for bH, in this case the forward and reverse rate constants are equal.
Reaction 6. Ubiquinol binding to and dissociation from the complex III at the cytosolic side (Qo) of the inner mitochondrial membrane. This is a way of transition between species: core (1) and species with occupied Qo site (2), and also (3) and (4):
A mathematical equation, expression, or formula.
 Object name is pcbi.1000619.e028.jpg
The expressions for forward and reverse reaction rates are:
A mathematical equation, expression, or formula.
 Object name is pcbi.1000619.e029.jpg
Simulating this and other binding/dissociation reactions the program works with two arrays, since substrate and product are belong to different species. This does not change much the principles of automated calculation described for the reaction 0. In the same way the program checks the requirements for substrate position (in this case all redox form satisfy it) calculates the position of product specific for each redox form of substrate, calculates the rates (12a) and adds them to the respective derivatives.
For the reaction (12)
A mathematical equation, expression, or formula.
 Object name is pcbi.1000619.e030.jpg
As a starting point we assumed that Cxxxx = Cxxxx11 when QH2 = Q. In this case Kd = 2 nmol/mg of protein.
Reaction 7. Ubiquinone binding and dissociation to the complex III at the matrix side of the inner mitochondrial membrane:
A mathematical equation, expression, or formula.
 Object name is pcbi.1000619.e031.jpg
The expressions for reaction rates:
A mathematical equation, expression, or formula.
 Object name is pcbi.1000619.e032.jpg
In simulations we used the same value:
A mathematical equation, expression, or formula.
 Object name is pcbi.1000619.e033.jpg
Reaction 8. Ubiquinone dissociation and binding to the complex III at the cytosolic side of the inner mitochondrial membrane:
A mathematical equation, expression, or formula.
 Object name is pcbi.1000619.e034.jpg
The expressions for reaction rates:
A mathematical equation, expression, or formula.
 Object name is pcbi.1000619.e035.jpg
Reaction 9. Ubiquinol dissociation and binding to the complex III at the matrix side of the inner mitochondrial membrane:
A mathematical equation, expression, or formula.
 Object name is pcbi.1000619.e036.jpg
The expressions for reaction rates:
A mathematical equation, expression, or formula.
 Object name is pcbi.1000619.e037.jpg
Reaction 10. ROS production in the model.
Semiquinone bound to the Qo site of complex III could be oxidized by molecular oxygen, thus producing superoxide radical:
A mathematical equation, expression, or formula.
 Object name is pcbi.1000619.e038.jpg
This reaction is considered as unidirectional:
A mathematical equation, expression, or formula.
 Object name is pcbi.1000619.e039.jpg
For relative ROS production kROS was taken to be 1 and KROS is specified in the legend to Figure 6.
Reaction 11. Succinate oxidation in complex II resulted in ubiquinone reduction (Q to QH2), which rate V is described by the Michaelis' equation:
A mathematical equation, expression, or formula.
 Object name is pcbi.1000619.e040.jpg
where KQ = 0.5 nmol/mg of protein, KS = 0.25 mM of protein, Vm = 1.88 nmol/mg of protein/s, and S is external succinate concentration. Thus, equation (17) represents succinate dehydrogenase activity in a generalized form not going into details as in above reactions.
Reaction 12. Electron transport from complex III to cytochrome c and further to complex IV.
All species (1)–(4) contribute to this rection:
A mathematical equation, expression, or formula.
 Object name is pcbi.1000619.e041.jpg
The expression for this reaction rate, which in fact integrates all the steps of electron transition from cytochrome c1 to oxygen, takes into account oxygen concentration:
A mathematical equation, expression, or formula.
 Object name is pcbi.1000619.e042.jpg
here O2 is a fraction of dissolved oxygen with respect to its content at equilibrium with atmospheric oxygen at normal pressure, k is reaction rate constant, Km is Michaelis constant for interaction with oxygen. These parameters were adjusted so that the model describes limitation of electron flow by oxygen, which was found experimentally. According to the experimental data [22],[13], the electron flow is not limited by O2 availability until the oxygen concentration falls to about 1 µM (<0.1% O2). The model reproduces these data if Km = 0.03, and k = 25 s−1. These parameters were used to calculate the oxygen dependence of semiquinone radical content presented in Figure 5A of the main text.
Complex III content is 0.4 nmol/mg of protein [33], and total ubiquinone content is 4 nmol/mg of protein [33].
Automated construction of differential equations
The system of differential equations describes the evolution of concentrations (C) of all redox forms, which the program references as is described above. Derivative of each redox state is calculated as a sum of all the reaction rates of influx to or efflux from the considered state, which are calculated in accordance with the equations (6–18). For each type of reaction there is a special function, which, once called for a specie, passes from redox state 0 to the last one checking if the redox conditions permitting the reaction are satisfied and, if yes, calculates the flux and adds it (with respective sign) to the derivatives of concentrations of reagents. In this way, to simulate the operation of respiratory complex it is necessary just to check whether the functions simulating all the reactions performed by the species (1)–(4) are called, and these functions automatically would calculate the derivatives for all 400 redox states of complex III.
For the numerical solutions of differential equations the program implements 5th order Runge-Kutta, Bulirch-Stoer, implicit Bulirch-Stoer, Rosenberg methods [34]. The program and brief user guide is free available in
Experimental procedures
Similar technique is described in more detail elsewhere [35],[36].
Isolation of rat brain mitochondria - All procedures involving animals were approved by the Children's Hospital of Pittsburgh and were in compliance with “Principles of Laboratory Animal Care” and the current laws of the United States. Rat brain mitochondria were isolated from the cortex of adult Wistar rats. After removal, tissue was minced and homogenized in ice-cold isolation buffer I (IB I) which contained: 225 mM mannitol, 75 mM sucrose, 5 mM HEPES buffer (pH adjusted to 7.3 with KOH), 0.1 mg/ml fatty acid free BSA, 1 mM tetrapotassium EDTA and 12% Percoll. The homogenate thus obtained was carefully layered on the top of a discontinuous gradient of Percoll (24% and 42%) prepared using the same buffer. The preparation was then centrifuged at 31,000×g for 10 min. The fraction containing the mitochondria located between 42% and 24% Percoll was carefully withdrawn by a syringe and washed from Percoll twice by pelleting in IB I. The resulting mitochondrial suspension was diluted in isolation buffer II (IB II), which was same as IB I, except for the concentration of EDTA (0.1 mM) and lack of albumin, and spun down at 12,000×g for 10 min. The deposit of mitochondria was homogenized in IB II at a final protein concentration of ~20 mg/ml and stored on ice until use. The protein concentration in the mitochondrial samples was determined using a Protein Assay kit (Pierce Chemical Company, Rockford IL) according to the manufacture's instructions. Mitochondria prepared in this way were active for at least 4–5 hours, as determined by their ability to maintain a stable transmembrane potential in the presence of oxidizable substrates.
Hydrogen peroxide measurements were performed in a stirred cuvette mounted in a Shimatzu RF-5301 spectrofluorimeter maintained at 37°C. Mitochondria (0.2 mg/ml of protein) were added to the incubation medium that contained: 125 mM KCl; 2 mM KH2 PO4; 2 mM MgCl2; 10 mM Tris;10 mM HEPES (pH 7.0); 100 µM EGTA; 5 mM pyruvate and particular concentrations of succinate (0.05–5 mM) as the oxidizable substrates. Reporting fluorescent dye for H2O2 was Amplex red (2 µM) which increased fluorescence upon oxidation to resorufin in the presence of 1 U/ml of horseradish peroxidase (HRP) as previously described [36]. Measurements were carried out at excitation/emission wavelengths of 560 nm (slit 1.5 nm)/590 nm (slit 3nm), respectively. Amounts of H2O2 released by mitochondria were estimated by constructing calibration curves using known H2O2 concentrations in the standard incubation buffer together with Amplex red and HRP, but without mitochondria. Normoxia experiments were performed in a conventional square (10×10 mm), open cuvette. Anoxic experiments were performed in a cuvette with a narrow rectangular neck (10×3 mM). Anoxic conditions were made by bubbling the cuvette with N2 and, once the reaction media was depleted of O2, closing it with a Teflon lid. Resumption of oxygen supply was achieved by opening the cuvette and bubbling it with air.
Mitochondrial transmembrane potential, ΔΨm - was estimated using fluorescence quenching of the cationic dye safranine O. Since polarized mitochondria have a negative charge inside, positively charged molecules of safranine O are accumulated inside the matrix; increase in dye concentration inside the matrix leads to fluorescence quenching , thus a decrease in fluorescence corresponds to an increase of membrane potential. The excitation wavelength was 495 nm (slit 3nm) and emission 586 nm (slit 5nm), and the dye concentration used was 2.5 µM [36].
Ethics Statement
All procedures involving animals were approved by the Children's Hospital of Pittsburgh and were in compliance with “Principles of Laboratory Animal Care” and the current laws of the United States.
Text S1
Text file describing the algorithms, parameters, experiments referred in the main text of the manuscript
(0.16 MB PDF)
Figure S1
Scheme of the processes in mitochondrial electron transport chain and their coupling with central energetic metabolism. Respiratory complexes (CI–CIV) deliver electrons from NADH oxidizing it to NAD+ to oxygen reducing it to H2O. This delivery is subdivided into many steps, some of which are coupled with proton transport from matrix to cytosol. Redox reactions in electron transport chain and proton translocation are fed by glycolysis and TCA cycle, which provide NADH, complex 1 substrate, and succinate, complex II substrate. Ubiquinone/ubiquinol (Q/QH2) cycle plays essential role in proton transport and communication between CI–CIII. The scheme of electron transport coupled with proton translocation in the complex III (Q-cycle mechanism) is shown in more detail in the right panel. Oxidizes ubiquinol reducing cytochrome C and releasing two protons to the cytosolic side (positive or p-side; this is reflected in the index of released H+). In addition it translocates two protons from matrix (negative or n-side) to cytosol. Ubiquinol (QH2) delivers its first electron to Fe(3+), releasing 2 protons at the p-side of the inner mitochondrial membrane and producing semiquinone (SQ) radical. Then the latter gives its unpaired electron to cytochrome bl and produced ubiquinone (Q) dissociates from the complex. Free Q is bound at the n-side, receives two electrons from cytochrome bh resulted from oxidation of two QH2 molecules, thus producing subsequently SQ and QH2, taking protons from n-side. Dissociation of the produced QH2 accomplishes a round of the cycle.
(0.06 MB TIF)
Figure S2
Validation of the model by qualitative simulation of the observed experimentally decrease of ROS production and increase of maximal flow rate induced by ADP, and increase of ROS production induced by hypoxia. (a) Decrease of steady state content of semiquinone bound at Qo site (ROS production rate) in conjunction with a decrease of transmembrane potential. Blue squares indicate the steady state levels of SQ bound at Qo site reached from initial states with high SQ levels at various relative succinate concentrations. Orange diamonds indicate the steady state SQ at Qo levels reached from initial states with low SQ levels. These two curves, which indicate bistable behavior of complex III operation, were obtained for transmembrane potential 200 mV that simulates state 4 of respiration. If transmembrane potential falls down to below 150 mV (as could be after a transition to state 3 induced by addition of ADP that stimulates ATP synthesis), the area of bistability is shifted to the right (as indicated by yellow and green triangles) and the same substrate supply corresponds now to low semiquinone content. Thus, if the system operates in high ROS producing mode at the relative succinate levels around 0.8–0.9, addition of ADP would switch it to the low ROS producing mode. (b) Increase of electron flow through the respiratory chain in conjunction with a decrease of transmembrane potential. The labels of curves correspond to the same conditions as in Figure S2a. If the system persists in state 4 of respiration in a high ROS producing mode (orange diamonds), the addition of ADP (transition to state 3) would switch the electron flow to the values indicated by green triangles, and if the substrate supply is around 0.8–0.9, electron flow would increase from 0.4 to 1.4 ng.atom of O/mg/s. (c) Initially mitochondria incubated with succinate are in state 4 and in high ROS production state (orange curve in Figure S2a). The addition of ADP, simulated as decrease of transmembrane potential, induces transition to low ROS production state. After ADP is consumed transmembrane potential increased and ROS production changes to the state indicated by blue curve in Figure S2a. Hypoxia, simulated as further increase of substrate supply and decrease of rate constant of cytochrome c1 oxidation, provoked switch to back to the state indicated by orange curve.
(0.08 MB TIF)
Figure S3
Analysis of sensitivity of the bistability region to the change of model parameters. The basic model parameters, which are listed above and in the main text. To show the sensitivity of bistability region the parameters were changed, which induced the change of area of bistability. Blue and orange lines indicate the bistable region, which corresponds to the indicated above set of parameters. Yellow and green lines illustrate how this region changes in response to the change of a parameter as indicated in the figures.
(0.11 MB TIF)
Figure S4
Analysis of sensitivity of the bistability region to the change of model parameters. Designations are the same as in Figure S3
(0.12 MB TIF)
Figure S5
Distribution of redox states of complex III, which corresponds to low (upper panel), high (lower panel) ROS production states or intermediate non-steady state reached in 3 second after the switch of substrate supply which induced the transition from low to high ROS producing steady state (middle panel). The redox states are positioned as follows: the numbers 0–15 correspond to 16 states of bh-bl-c1-FeS, numbers 16–79 correspond to 64 states of bh-bl-c1-FeS-Q-Q, numbers 80–143 correspond to 64 states of Q-Q-bh-bl-c1-FeS, numbers 144–399 correspond to 256 states of Q-Q-bh-bl-c1-FeS-Q-Q.
(0.08 MB TIF)
Figure S6
(a) Temporary stimulation by ADP switches the electron transport from high to low ROS production mode. Rat brain mitochondria were added to the incubation medium, with subsequent addition of 1 mM ADP (red curves), or with 1 mM of ADP present initially (blue curves). Without ADP, ROS production is much higher than in its presence, however, after ADP addition it decreases. During the transformation of ADP into ATP the membrane potential remains low but, when all added ADP is transformed, the potential increases (curves MP). However, after the completion of transformation, ROS production increased, but not to the original level. Since the produced ATP itself does not affect ROS production, this data indicate that ADP plays a role of trigger, which switches the respiratory chain from high to low ROS production. (b) ATP does not affect ROS production, but adventitious ADP in the ATP preparation was sufficient to switch the electron transport to low ROS producing mode. Blue curves are shown for comparison; they are taken from an experiment, similar to that shown in (a). Black curve shows that ATP added with oligomycin does not change ROS production. Green curves show that a presence of traces of ADP in the preparation of ATP switches a part of mitochondria to low ROS producing mode, and that the addition of oligomycin after ATP does not bring the ROS rate to the original level.
(0.05 MB TIF)
The authors have declared that no competing interests exist.
This work was supported by the European Commission (FP6, BioBridge LSHG-CT-2006-037939, (FP7, Diaprepp Health-F2-2008-202013, Etherpath KBBE-grant agreement n°222639; the Spanish Government and the European Union FEDER funds ( SAF2008-00164 and ISCIII-RTICC (RD06/0020/0046); Generalitat de Catalunya, (2005SGR00204; 2006ITT-10007); Foundation Marato TV3-042010 ( and the National Institutes of Health (, AG 20899 and U01 DK61058-01; and the Department of Defense W81XWH-06-1-0317. The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript.
1. Li C, Jackson RM. Reactive species mechanisms of cellular hypoxia–reoxygenation injury. Am J Physiol. 2002;282:C227–C241. [PubMed]
2. Bottino R, Balamurugan AN, Tse H, Thirunavukkarasu C, Ge X, et al. Response of human islets to isolation stress and the effect of antioxidant treatment. Diabetes. 2004;53:2559–2568. [PubMed]
3. Guzy RD, Hoyos B, Robin E, Chen H, Liu L, et al. Mitochondrial complex III is required for hypoxia-induced ROS production and cellular oxygen sensing. Cell Metab. 2005;1:401–408. [PubMed]
4. Cino M, Del Maestro RF. Generation of hydrogen peroxide by brain mitochondria: the effect of reoxygenation following postdecapitative ischemia. Arch Biochem Biophys. 1989;269:623–638. [PubMed]
5. Cai J, Jones DP. Mitochondrial redox signaling during apoptosis. J Bioenerg Biomembr. 1999;31:327–334. [PubMed]
6. Turrens JF, Alexandre A, Lehninger AL. Ubisemiquinone is the electron donor for superoxide formation by complex III of heart mitochondria. Arch Biochem Biophys. 1985;237:408–414. [PubMed]
7. Ksenzenko M, Konstantinov AA, Khomutov GB, Tikhonov AN, Ruuge EK. Effect of electron transfer inhibitors on superoxide generation in the cytochrome bc1 site of the mitochondrial respiratory chain. FEBS Lett. 1983;155:19–24. [PubMed]
8. Vinogradov AD. Catalytic properties of the mitochondrial NADH–ubiquinone oxidoreductase (Complex I) and the pseudo-reversible active/inactive enzyme transition. Biochim Biophys Acta. 1998;1364:169–185. [PubMed]
9. Ohnishi T, Johnson JE, Jr, Yano T, LoBruttoc R, Widgerb WR. Thermodynamic and EPR studies of slowly relaxing ubisemiquinone species in the isolated bovine heart complex I. FEBS Lett. 2005;579:500–506. [PubMed]
10. Moser CC, Farid TA, Chobot SE, Dutton PL. Electron tunneling chains of mitochondria. Biochim Biophys Acta. 2006;1757:1096–1109. [PubMed]
11. Mulkidjanian AY. Ubiquinol oxidation in the cytochrome bc 1 complex: Reaction mechanism and prevention of short-circuiting. Biochim Biophys Acta. 2005;1709:5–34. [PubMed]
12. Cooley JW, Ohnishi T, Fevzi D. Binding dynamics at the quinone reduction (Qi) site influence the equilibrium interactions of the iron sulfur protein and hydroquinone oxidation (Qo) site of the cytochrome bc1 complex. Biochemistry. 2005;44:10520–10532. [PMC free article] [PubMed]
13. Guzy RD, Schumacker PT. Oxygen sensing by mitochondria at complex III: the paradox of increased reactive oxygen species during hypoxia. Exp Physiol. 2006;91:807–819. [PubMed]
14. Ohnishi T, Trumpower BL. Differential effects of antimycin on ubisemiquinone bound in different environments in isolated succinate-cytochrome c reductase complex. J Biol Chem. 1980;255:3278–3284. [PubMed]
15. Chandel NS, Maltepe E, Goldwasser E, Mathieu CE, Simon MC, et al. Mitochondrial reactive oxygen species trigger hypoxia-induced transcription. Proc Natl Acad Sci U S A. 1998;95:11715–11720. [PubMed]
16. Bell EL, Cai L, Johnstone BH, Traktuev DO, March KL, et al. Mitochondrial reactive oxygen species trigger hypoxia-inducible factor-dependent extension of the replicative life span during hypoxia. Mol Cell Biol. 2007;27:5737–5745. [PMC free article] [PubMed]
17. Selivanov VA, Puigjaner J, Sillero A, Centelles JJ, Ramos-Montoya A, et al. An optimized algorithm for flux estimation from isotopomer distribution in glucose metabolites. Bioinformatics. 2004;20:3387–3397. [PubMed]
18. Selivanov VA, Meshalkina LE, Solovjeva ON, Kuchel PW, Ramos-Montoya A, et al. Rapid simulation and analysis of isotopomer distributions using constraints based on enzyme mechanisms: an example from HT29 cancer cells. Bioinformatics. 2005;21:3558–3564. [PubMed]
19. Selivanov VA, Marin S, Lee PWN, Cascante M. Software for dynamic analysis of tracer-based metabolomic data: estimation of metabolic fluxes and their statistical analysis. Bioinformatics. 2006;22:2806–2812. [PubMed]
20. Ohnishi ST, Ohnishi T, Muranaka S, Fujita H, Kimura H, et al. A Possible Site of Superoxide Generation in the Complex I Segment of Rat Heart Mitochondria. J Bioenerg Biomembr. 2005;37:1–15. [PubMed]
21. Pisarenko O, Studneva I, Khlopkov V, Solomatina E, Ruuge E. An assessment of anaerobic metabolism during ischemia and reperfusion in isolated guinea pig heart. Biochim Biophys Acta. 1988;934:55–63. [PubMed]
22. Wilson DF, Rumsey WL, Green TJ, Vanderkooi JM. The oxygen dependence of mitochondrial oxidative phosphorylation measured by a new optical method for measuring oxygen concentration. J Biol Chem. 1988;263:2712–2718. [PubMed]
23. Jiang BH, Semenza GL, Bauer C, Marti HH. Hypoxia-inducible factor 1 levels vary exponentially over a physiologically relevant range of O2 tension. Am J Physiol. 1996;271:C1172–C1180. [PubMed]
24. Schwartz LM, Reimer KA, Crago MS, Jennings RB. Pharmacological preconditioning with diazoxide slows energy metabolism during sustained ischemia. Exp Clin Cardiol. 2007;12:139–147. [PMC free article] [PubMed]
25. Link TA, Hagen WR, Pierik AJ, Assmann C, von Jagow G. Determination of the redox properties of the Rieske [2Fe-2S] cluster of bovine heart bc1 complex by direct electrochemistry of a water-soluble fragment. Eur J Biochem. 1992;208:685–691. [PubMed]
26. Trumpower BL. Function of the iron-sulfur protein of the cytochrome b-c1 segment in electron-transfer and energy-conserving reactions of the mitochondrial respiratory chain. Biochim Biophys Acta. 1981;639:129–155. [PubMed]
27. Roginsky VA, Barsukova TK, Stegmann HB. Kinetics of redox interaction between substituted quinones and ascorbate under aerobic conditions. Chem Biol Interact. 1999;121:177–197. [PubMed]
28. Snyder CH, Gutierrez-Cirlos EB, Trumpower BL. Evidence for a concerted mechanism of ubiquinol oxidation by the cytochrome bc1 complex. J Biol Chem. 2000;275:13535–13541. [PubMed]
29. Leguijt T, Engels PW, Crielaard W, Albracht SP, Hellingwerf KJ. Abundance, subunit composition, redox properties, and catalytic activity of the cytochrome bc1 complex from alkaliphilic and halophilic, photosynthetic members of the family Ectothiorhodospiraceae. J Bacteriol. 1993;175:1629–1636. [PMC free article] [PubMed]
30. Covian R, Zwicker K, Rotsaert FA, Trumpower BL. Asymmetric and redox-specific binding of quinone and quinol at center N of the dimeric yeast cytochrome bc1 complex. Consequences for semiquinone stabilization. J Biol Chem. 2007;282:24198–24208. [PubMed]
31. Reynolds IA. Incorporation of membrane potential into theoretical analysis of electrogenic ion pumps. Proc Natl Acad Sci USA. 1985;82:6869–6873. [PubMed]
32. Rich PR. Electron and proton transfers through quinones and cytochrome bc complexes. Biochim Biophys Acta. 1984;768:53–79. [PubMed]
33. Demin OV, Kholodenko BN, Skulachev VP. A model of O2 generation in the complex III of the electron transport chain. Mol Cell Biochem. 1998;184:21–33. [PubMed]
34. Press WH, Flannery BP, Teukolsky SA, Vetterling WT. Numerical Recipes in C: The Art of Scientific Computing. NY: Cambridge University Press; 2002.
35. Selivanov VA, Zeak JA, Roca J, Cascante M, Trucco M. The role of external and matrix pH in mitochondrial reactive oxygen species generation. J Biol Chem. 2008;283:29292–29300. [PMC free article] [PubMed]
36. Votyakova TV, Reynolds IJ. DeltaPsi(m)-Dependent and -independent production of reactive oxygen species by rat brain mitochondria. J Neurochem. 2001;79:266–277. [PubMed]
Articles from PLoS Computational Biology are provided here courtesy of
Public Library of Science