Mutations in cytochrome b that affect kinetics of the electron transfer reactions at center N in the yeast cytochrome bc1 complex

We have examined the pre-steady-state kinetics and thermodynamic properties of the b hemes in variants of the yeast cytochrome bc₁ complex that have mutations in the quinone reductase site (center N). Trp-30 is a highly conserved residue, forming a hydrogen bond with the propionate on the high potential b heme (b_H heme). The substitution by a cysteine (W30C) lowers the redox potential of the heme and an apparent consequence is a lower rate of electron transfer between quinol and heme at center N. Leu-198 is also in close proximity to the b_H heme and a L198F mutation alters the spectral properties of the heme but has only minor effects on its redox properties or the electron transfer kinetics at center N. Substitution of Met-221 by glutamine or glutamate results in the loss of a hydrophobic interaction that stabilizes the quinone ligands. Ser-20 and Gln-22 form a hydrogen-bonding network that includes His-202, one of the carbonyl groups of the ubiquinone ring, and an active-site water. A S20T mutation has long-range structural effects on center P and thermodynamic effects on both b hemes. The other mutations (M221E, M221Q, Q22E and Q22T) do not affect the ubiquinol oxidation kinetics at center P, but do modify the electron transfer reactions at center N to various extents. The pre-steady reduction kinetics suggest that these mutations alter the binding of quinone ligands at center N, possibly by widening the binding pocket and thus increasing the distance between the substrate and the b_H heme. These results show that one can distinguish between the contribution of structural and thermodynamic factors to center N function.     

Domain movement of iron sulfur protein in cytochrome bc1 complex is facilitated by the electron transfer from cytochrome b(L) to b(H)

The key step of the "protonmotive Q-cycle" mechanism for cytochrome bc1 complex is the bifurcated oxidation of ubiquinol at the Qp site. ISP is reduced when its head domain is at the b-position and subsequent move to the c1 position, to reduce cytochrome c1, upon protein conformational changes caused by the electron transfer from cytochrome b(L) to b(H). Results of analyses of the inhibitory efficacy and the binding affinity, determined by isothermal titration calorimetry, of Pm and Pf, on different redox states of cytochrome bc1 complexes, confirm this speculation. Pm inhibitor has a higher affinity and better efficacy with the cytochrome b(H) reduced complex and Pf binds better and has a higher efficacy with the ISP reduced complex.     

Thermodynamic and EPR studies of slowly relaxing ubisemiquinone species in the isolated bovine heart complex I

Previously, we investigated ubisemiquinone (SQ) EPR spectra associated with NADH-ubiquinone oxidoreductase (complex I) in the tightly coupled bovine heart submitochondrial particles (SMP). Based upon their widely differing spin relaxation rate, we distinguished SQ spectra arising from three distinct SQ species, namely SQ(Nf) (fast), SQ(Ns) (slow), and SQ(Nx) (very slow). The SQ(Nf) signal was observed only in the presence of the proton electrochemical gradient (deltamu(H)(+)), while SQ(Ns) and SQ(Nx) species did not require the presence of deltamu(H+). We have now succeeded in characterizing the redox and EPR properties of SQ species in the isolated bovine heart complex I. The potentiometric redox titration of the g(z,y,x)=2.00 semiquinone signal gave the redox midpoint potential (E(m)) at pH 7.8 for the first electron transfer step [E(m1)(Q/SQ)] of -45 mV and the second step [E(m2)(SQ/QH(2))] of -63 mV. It can also be expressed as [E(m)(Q/QH(2))] of -54 mV for the overall two electron transfer with a stability constant (K(stab)) of the SQ form as 2.0. These characteristics revealed the existence of a thermodynamically stable intermediate redox state, which allows this protein-associated quinone to function as a converter between n=1 and n=2 electron transfer steps. The EPR spectrum of the SQ species in complex I exhibits a Gaussian-type spectrum with the peak-to-peak line width of approximately 6.1 G at the sample temperature of 173 K. This indicates that the SQ species is in an anionic Q(-) state in the physiological pH range. The spin relaxation rate of the SQ species in isolated complex I is much slower than the SQ counterparts in the complex I in situ in SMP. We tentatively assigned slow relaxing anionic SQ species as SQ(Ns), based on the monophasic power saturation profile and several fold increase of its spin relaxation rate in the presence of reduced cluster N2. The current study also suggests that the very slowly relaxing SQ(Nx) species may not be an intrinsic complex I component. The functional role of SQ(Ns) is further discussed in connection with the SQ(Nf) species defined in SMP in situ.     

Intermonomer electron transfer between the low-potential b hemes of cytochrome bc₁

Cytochrome (cyt) bc(1) is a structural dimer with its monomers consisting of the Fe-S protein, cyt b, and cyt c(1) subunits. Its three-dimensional architecture depicts it as a symmetrical homodimer, but the mobility of the head domain of the Fe-S protein indicates that the functional enzyme exists in asymmetrical heterodimeric conformations. Here, we report a new genetic system for studying intra- and intermonomer interactions within the cyt bc(1) using the facultative phototrophic bacterium Rhodobacter capsulatus. The system involves two different sets of independently expressed cyt bc(1) structural genes carried by two plasmids that are coharbored by a cell without its endogenous enzyme. Our results indicate that coexpressed cyt bc(1) subunits were matured, assorted, and assembled in vivo into homo- and heterodimeric enzymes that can bear different mutations in each monomer. Using the system, the occurrence of intermonomer electron transfer between the low-potential b hemes of cyt bc(1) was probed by choosing mutations that perturb electron transfer at the hydroquinone oxidation (Q(o)) and quinone reduction (Q(i)) sites of the enzyme. The data demonstrate that active heterodimeric variants, formed of monomers carrying mutations that abolish only one of the two (Q(o) or Q(i)) active sites of each monomer, are produced, and they support photosynthetic growth of R. capsulatus. Detailed analyses of the physicochemical properties of membranes of these mutants, as well as purified homo- and heterodimeric cyt bc(1) preparations, demonstrated that efficient and productive electron transfer occurs between the low-potential b(L) hemes of the monomers in a heterodimeric enzyme. Overall findings are discussed with respect to intra- and intermonomer interactions that take place during the catalytic turnover of cyt bc(1).     

Respiratory cytochrome c oxidase can be efficiently reduced by the photosynthetic redox proteins cytochrome c6 and plastocyanin in cyanobacteria

Plastocyanin and cytochrome c6 are two small soluble electron carriers located in the intrathylacoidal space of cyanobacteria. Although their role as electron shuttle between the cytochrome b6f and photosystem I complexes in the photosynthetic pathway is well established, their participation in the respiratory electron transport chain as donors to the terminal oxidase is still under debate. Here, we present the first time-resolved analysis showing that both cytochrome c6 and plastocyanin can be efficiently oxidized by the aa3 type cytochrome c oxidase in Nostoc sp. PCC 7119. The apparent electron transfer rate constants are ca. 250 and 300 s(-1) for cytochrome c6 and plastocyanin, respectively. These constants are 10 times higher than those obtained for the oxidation of horse cytochrome c by the oxidase, in spite of being a reaction thermodynamically more favourable.     

Conformation-driven and semiquinone-gated proton-pump mechanism in the NADH-ubiquinone oxidoreductase (complex I)

A novel mechanism for proton/electron transfer is proposed for NADH-quinone oxidoreductase (complex I) based on the following findings: (1) EPR signals of the protein-bound fast-relaxing semiquinone anion radicals (abbreviated as Q(Nf)-) are observable only in the presence of proton-transmembrane electrochemical potential; (2) Iron-sulfur cluster N2 and Q(Nf)- are directly spin-coupled; and (3) The projection of the interspin vector extends only 5A along the membrane normal [Yano, T., Dunham, W.R. and Ohnishi, T. (2005) Biochemistry, 44, 1744-1754]. We propose that the proton pump is operated by redox-driven conformational changes of the quinone binding protein. In the input state, semiquinone is reduced to quinol, acquiring two protons from the N (matrix) side of the mitochondrial inner membrane and an electron from the low potential (NADH) side of the respiratory chain. A conformational change brings the protons into position for release at the P (inter-membrane space) side of the membrane via a proton-well. Concomitantly, an electron is donated to the quinone pool at the high potential side of the coupling site. The system then returns to the original state to repeat the cycle. This hypothesis provides a useful frame work for further investigation of the mechanism of proton translocation in complex I.     

An electron paramagnetic resonance investigation of the electron transfer reactions in the chlorophyll d containing photosystem I of Acaryochloris marina

Electron paramagnetic resonance (EPR) spectroscopy reveals functional and structural similarities between the reaction centres of the chlorophyll d-binding photosystem I (PS I) and chlorophyll a-binding PS I. Continuous wave EPR spectrometry at 12K identifies iron-sulphur centres as terminal electron acceptors of chlorophyll d-binding PS I. A transient light-induced electron spin echo (ESE) signal indicates the presence of a quinone as the secondary electron acceptor (Q) between P(740)(+) and the iron-sulphur centres. The distance between P(740)(+) and Q(-) was estimated within point-dipole approximation as 25.23+/-0.05A, by the analysis of the electron spin echo envelope modulation.     

Rapid electron transfer between monomers when the cytochrome bc1 complex dimer is reduced through center N

We have obtained evidence for electron transfer between cytochrome b subunits of the yeast bc(1) complex dimer by analyzing pre-steady state reduction of cytochrome b in the presence of center P inhibitors. The kinetics and extent of cytochrome b reduced by quinol in the presence of variable concentrations of antimycin decreased non-linearly and could only be fitted to a model in which electrons entering through one center N can equilibrate between the two cytochrome b subunits of the bc(1) complex dimer. The b(H) heme absorbance in a bc(1) complex inhibited at center P and preincubated with substoichiometric concentrations of antimycin showed a red shift upon the addition of substrate, which indicates that electrons from the uninhibited center N in one monomer are able to reach the b(H) heme at the antimycin-blocked site in the other. The extent of cytochrome b reduction by variable concentrations of menaquinol could only be fitted to a kinetic model that assumes electron equilibration between center N sites in the dimer. Kinetic simulations showed that non-rate-limiting electron equilibration between the two b(H) hemes in the dimer through the two b(L) hemes is possible upon reduction through one center N despite the thermodynamically unfavorable b(H) to b(L) electron transfer step. We propose that electron transfer between cytochrome b subunits minimizes the formation of semiquinone-ferrocytochrome b(H) complexes at center N and favors ubiquinol oxidation at center P by increasing the amount of oxidized cytochrome b.     

Nitric oxide reacts with the ferryl-oxo catalytic intermediate of the CuB-lacking cytochrome bd terminal oxidase

Cytochrome bd is a bacterial respiratory oxidase carrying three hemes but no copper. We show that nitric oxide (NO) reacts with the intermediate F of cytochrome bd from Azotobacter vinelandii: (i) with a 1:1 stoichiometry, (ii) rapidly (k=1.2 +/- 0.1 x 10(5)M(-1)s(-1) at 20 degrees C), and (iii) yielding the oxidized enzyme with nitrite bound to heme d at the active site. Unexpectedly, the NO reaction mechanism of this catalytic intermediate in the Cu(B)-lacking cytochrome bd appears similar to that of beef heart cytochrome c oxidase, where Cu(B) was proposed to play a key role.     

Evidence for electron equilibrium between the two hemes bL in the dimeric cytochrome bc1 complex

Structural analysis of the dimeric mitochondrial cytochrome bc1 complex suggests that electron transfer between inter-monomer hemes bL-bL may occur during bc1 catalysis. Such electron transfer may be facilitated by the aromatic pairs present between the two bL hemes in the two symmetry-related monomers. To test this hypothesis, R. sphaeroides mutants expressing His6-tagged bc1 complexes with mutations at three aromatic residues (Phe-195, Tyr-199, and Phe-203), located between two bL hemes, were generated and characterized. All three mutants grew photosynthetically at a rate comparable to that of wild-type cells. The bc1 complexes prepared from mutants F195A, Y199A, and F203A have, respectively, 78%, 100%, and 100% of ubiquinol-cytochrome c reductase activity found in the wild-type complex. Replacing the Phe-195 of cytochrome b with Tyr, His, or Trp results in mutant complexes (F195Y, F195H, or F195W) having the same ubiquinol-cytochrome c reductase activity as the wild-type. These results indicate that the aromatic group at position195 of cytochrome b is involved in electron transfer reactions of the bc1 complex. The rate of superoxide anion (O2*) generation, measured by the chemiluminescence of 2-methyl-6-(p-methoxyphenyl)-3,7-dihydroimidazo[1,2-alpha]pyrazin-3-one hydrochloride-O2* adduct during oxidation of ubiquinol, is 3 times higher in the F195A complex than in the wild-type or mutant complexes Y199A or F203A. This supports the idea that the interruption of electron transfer between the two bL hemes enhances electron leakage to oxygen and thus decreases the ubiquinol-cytochrome c reductase activity.     

Crystal structure of oxidized cytochrome c(6A) from Arabidopsis thaliana

Compared with algal and cyanobacterial cytochrome c(6), cytochrome c(6A) from higher plants contains an additional loop of 12 amino acid residues. We have determined the first crystal structure of cytochrome c(6A) from Arabidopsis thaliana at 1.5 Angstrom resolution in order to help elucidate its function. The overall structure of cytochrome c(6A) follows the topology of class I c-type cytochromes in which the heme prosthetic group covalently binds to Cys16 and Cys19, and the iron has octahedral coordination with His20 and Met60 as the axial ligands. Two cysteine residues (Cys67 and Cys73) within the characteristic 12 amino acids loop form a disulfide bond, contributing to the structural stability of cytochrome c(6A). Our model provides a chemical basis for the known low redox potential of cytochrome c(6A) which makes it an unsuitable electron carrier between cytochrome b(6)f and PSI.     

Computation of threshold conditions for epidemiological models and global stability of the disease-free equilibrium (DFE)

One goal of this paper is to give an algorithm for computing a threshold condition for epidemiological systems arising from compartmental deterministic modeling. We calculate a threshold condition T(0) of the parameters of the system such that if T(0)<1 the disease-free equilibrium (DFE) is locally asymptotically stable (LAS), and if T(0)>1, the DFE is unstable. The second objective, by adding some reasonable assumptions, is to give, depending on the model, necessary and sufficient conditions for global asymptotic stability (GAS) of the DFE. In many cases, we can prove that a necessary and sufficient condition for the global asymptotic stability of the DFE is R(0)< or =1, where R(0) is the basic reproduction number [O. Diekmann, J.A. Heesterbeek, Mathematical Epidemiology of Infectious Diseases: Model Building, Analysis and Interpretation, Wiley, New York, 2000]. To illustrate our results, we apply our techniques to examples taken from the literature. In these examples we improve the results already obtained for the GAS of the DFE. We show that our algorithm is relevant for high dimensional epidemiological models.     

The antioxidative function of eicosapentaenoic acid in a marine bacterium, Shewanella marinintestina IK-1

When the eicosapentaenoic acid (EPA)-deficient mutant strain IK-1Delta8 of the marine EPA-producing Shewanella marinintestina IK-1 was treated with various concentrations of hydrogen peroxide (H(2)O(2)), its colony-forming ability decreased more than that of the wild type. Protein carbonylation, induced by treating cells with 0.01 mM H(2)O(2) under bacteriostatic conditions, was enhanced only in cells lacking EPA. The amount of cells recovered from the cultures was decreased more significantly by the presence of H(2)O(2) for cells lacking EPA than for those producing EPA. Treatment of the cells with 0.1 mM H(2)O(2) resulted in much lower intracellular concentrations of H(2)O(2) being consistently detected in cells with EPA than in those without EPA. These results suggest that cellular EPA can directly protect cells against oxidative damage by shielding the entry of exogenously added H(2)O(2) in S. marinintestina IK-1.     

Escherichia coli engineered to produce eicosapentaenoic acid becomes resistant against oxidative damages

The colony-forming ability of Escherichia coli genetically engineered to produce eicosapentaenoic acid (EPA) grown in 3mM hydrogen peroxide (H(2)O(2)) was similar to that of untreated cells. It was rapidly lost in the absence of EPA. H(2)O(2)-induced protein carbonylation was enhanced in cells lacking EPA. The fatty acid composition of the transformants was unaffected by H(2)O(2) treatment, but the amount of fatty acids decreased in cultures of cells lacking EPA and increased in cultures of cells producing EPA, suggesting that cellular EPA is stable in the presence of H(2)O(2) in vivo and may protect cells directly against oxidative damage. We discuss the possible role of EPA in partially blocking the penetration of H(2)O(2) into cells through membranes containing EPA.     

Conformation of the c552:aa3 electron transfer complex in Paracoccus denitrificans studied by EPR on oriented samples

The EPR spectral parameters of aa(3) oxidase and cyt c(552) from Paracoccus denitrificans were studied in purified oxidase and enriched cyt c(552). The orientation of the g-tensors of hemes a and c(552) were determined on partially ordered membranes, enriched cyt c(552) and a c(552):aa(3) subcomplex. The known correlation of g-tensor to molecular axes in histidine/methionine ligated hemes permits us to position cyt c(552) with respect to the parent membrane. Taken together with previous data on the interaction surface between aa(3) oxidase and cyt c(552), these results allow us to arrive at a single conformation for the c(552):aa(3) electron transfer complex.     

The electric field generated by photosynthetic reaction center induces rapid reversed electron transfer in the bc1 complex

The cytochrome bc(1) complex is the central enzyme of respiratory and photosynthetic electron-transfer chains. It couples the redox work of quinol oxidation and cytochrome reduction to the generation of a proton gradient needed for ATP synthesis. When the quinone processing Q(i)- and Q(o)-sites of the complex are inhibited by both antimycin and myxothiazol, the flash-induced kinetics of the b-heme chain, which transfers electrons between these sites, are also expected to be inhibited. However, we have observed in Rhodobacter sphaeroides chromatophores, that when a fraction of heme b(H) is reduced, flash excitation induces fast (half-time approximately 0.1 ms) oxidation of heme b(H), even in the presence of antimycin and myxothiazol. The sensitivity of this oxidation to ionophores and uncouplers, and the absence of any delay in the onset of this reaction, indicates that it is due to a reversal of electron transfer between b(L) and b(H) hemes, driven by the electrical field generated by the photosynthetic reaction center. In the presence of antimycin A, but absence of myxothiazol, the second and following flashes induce a similar ( approximately 0.1 ms) transient oxidation of approximately 10% of the cytochrome b(H) reduced on the first flash. From the observed amplitude of the field-induced oxidation of heme b(H), we estimate that the equilibrium constant for sharing one electron between hemes b(L) and b(H) is 10-15 at pH 7. The small value of this equilibrium constant modifies our understanding of the thermodynamics of the Q-cycle, especially in the context of a dimeric structure of bc(1) complex.     

Pressure equilibrium and jump study on unfolding of 23-kDa protein from spinach photosystem II

Pressure-induced unfolding of 23-kDa protein from spinach photosystem II has been systematically investigated at various experimental conditions. Thermodynamic equilibrium studies indicate that the protein is very sensitive to pressure. At 20 degrees C and pH 5.5, 23-kDa protein shows a reversible two-state unfolding transition under pressure with a midpoint near 160 MPa, which is much lower than most natural proteins studied to date. The free energy (DeltaG(u)) and volume change (DeltaV(u)) for the unfolding are 5.9 kcal/mol and -160 ml/mol, respectively. It was found that NaCl and sucrose significantly stabilize the protein from unfolding and the stabilization is associated not only with an increase in DeltaG(u) but also with a decrease in DeltaV(u). The pressure-jump studies of 23-kDa protein reveal a negative activation volume for unfolding (-66.2 ml/mol) and a positive activation volume for refolding (84.1 ml/mol), indicating that, in terms of system volume, the protein transition state lies between the folded and unfolded states. Examination of the temperature effect on the unfolding kinetics indicates that the thermal expansibility of the transition state and the unfolded state of 23-kDa protein are closer to each other and they are larger than that of the native state. The diverse pressure-refolding pathways of 23-kDa protein in some conditions were revealed in pressure-jump kinetics.     

Photo-induced cyclic electron transfer involving cytochrome bc1 complex and reaction center in the obligate aerobic phototroph Roseobacter denitrificans.

Flash-induced redox changes of b-type and c-type cytochromes have been studied in chromatophores from the aerobic photosynthetic bacterium Roseobacter denitrificans under redox-controlled conditions. The flash-oxidized primary donor P+ of the reaction center (RC) is rapidly re-reduced by heme H1 (Em,7 = 290 mV), heme H2 (Em,7 = 240 mV) or low-potential hemes L1/L2 (Em,7 = 90 mV) of the RC-bound tetraheme, depending on their redox state before photoexcitation. By titrating the extent of flash-induced low-potential heme oxidation, a midpoint potential equal to -50 mV has been determined for the primary quinone acceptor QA. Only the photo-oxidized heme H2 is re-reduced in tens of milliseconds, in a reaction sensitive to inhibitors of the bc1 complex, leading to the concomitant oxidation of a cytochrome c spectrally distinct from the RC-bound hemes. This reaction involves cytochrome c551 in a diffusional process. Participation of the bc1 complex in a cyclic electron transfer chain has been demonstrated by detection of flash-induced reduction of cytochrome b561, stimulated by antimycin and inhibited by myxothiazol. Cytochrome b561, reduced upon flash excitation, is re-oxidized slowly even in the absence of antimycin. The rate of reduction of cytochrome b561 in the presence of antimycin increases upon lowering the ambient redox potential, most likely reflecting the progressive prereduction of the ubiquinone pool. Chromatophores contain approximately 20 ubiquinone-10 molecules per RC. At the optimal redox poise, approximately 0.3 cytochrome b molecules per RC are reduced following flash excitation. Cytochrome b reduction titrates out at Eh < 100 mV, when low-potential heme(s) rapidly re-reduce P+ preventing cyclic electron transfer. Results can be rationalized in the framework of a Q-cycle-type model.     

Mutations in cytochrome b that affect kinetics of the electron transfer reactions at center N in the yeast cytochrome bc1 complex

We have examined the pre-steady-state kinetics and thermodynamic properties of the b hemes in variants of the yeast cytochrome bc1 complex that have mutations in the quinone reductase site (center N). Trp-30 is a highly conserved residue, forming a hydrogen bond with the propionate on the high potential b heme (bH heme). The substitution by a cysteine (W30C) lowers the redox potential of the heme and an apparent consequence is a lower rate of electron transfer between quinol and heme at center N. Leu-198 is also in close proximity to the b(H) heme and a L198F mutation alters the spectral properties of the heme but has only minor effects on its redox properties or the electron transfer kinetics at center N. Substitution of Met-221 by glutamine or glutamate results in the loss of a hydrophobic interaction that stabilizes the quinone ligands. Ser-20 and Gln-22 form a hydrogen-bonding network that includes His-202, one of the carbonyl groups of the ubiquinone ring, and an active-site water. A S20T mutation has long-range structural effects on center P and thermodynamic effects on both b hemes. The other mutations (M221E, M221Q, Q22E and Q22T) do not affect the ubiquinol oxidation kinetics at center P, but do modify the electron transfer reactions at center N to various extents. The pre-steady reduction kinetics suggest that these mutations alter the binding of quinone ligands at center N, possibly by widening the binding pocket and thus increasing the distance between the substrate and the bH heme. These results show that one can distinguish between the contribution of structural and thermodynamic factors to center N function.     

The evolution of n-player cooperation-threshold games and ESS bifurcations

An evolutionary game of individuals cooperating to obtain a collective benefit is here modelled as an n-player Prisoner's Dilemma game. With reference to biological situations, such as group foraging, we introduce a threshold condition in the number of cooperators required to obtain the collective benefit. In the simplest version, a three-player game, complex behaviour appears as the replicator dynamics exhibits a catastrophic event separating a parameter region allowing for coexistence of cooperators and defectors and a region of pure defection. Cooperation emerges through an ESS bifurcation, and cooperators only thrive beyond a critical point in cost-benefit space. Moreover, a repelling fixed point of the dynamics acts as a barrier to the introduction of cooperation in defecting populations. The results illustrate the qualitative difference between two-player games and multiple player games and thus the limitations to the generality of conclusions from two-player games. We present a procedure to find the evolutionarily stable strategies in any n-player game with cost and benefit depending on the number of cooperators. This was previously done by Motro [1991. Co-operation and defection: playing the field and the ESS. J. Theor. Biol. 151, 145-154] in the special cases of convex and concave benefit functions and constant cost.