# A Biochemical Approach to Study the Role of the Terminal Oxidases in Aerobic Respiration in Shewanella oneidensis MR-1

Sébastien Le Laz, Arlette Kpebe, Marielle Bauzan, Sabrina Lignon, Marc Rousset, Myriam Brugna

► To cite this version:

Sébastien Le Laz, Arlette Kpebe, Marielle Bauzan, Sabrina Lignon, Marc Rousset, et al. A Biochemical Approach to Study the Role of the Terminal Oxidases in Aerobic Respiration in Shewanella oneidensis MR-1. PLoS ONE, 2014, 9 (1), e0086343. doi:10.1371/journal.pone.0086343

HAL Id: hal-01494488

https://amu.hal.science/hal-01494488v1

Submitted on 5 Oct 2018

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_1.png" /></figure>

Distributed under a Creative Commons CC BY 4.0 - Attribution - International License
# A Biochemical Approach to Study the Role of the Terminal Oxidases in Aerobic Respiration in *Shewanella oneidensis* MR-1

Sébastien Le Laz¹, Arlette Kpebe¹, Marielle Bauzan², Sabrina Lignon³, Marc Rousset¹, Myriam Brugna¹*

1 CNRS, Aix-Marseille Université, Laboratoire de Bioénergétique et Ingénierie des Proteines, UMR 7281, IMM, Marseille, France, 2 CNRS, Aix-Marseille Université, Unité de fermentation, FR3479, IMM, Marseille, France, 3 CNRS, Aix-Marseille Université, Plate-forme Protéomique, FR3479, IMM, MaP IBiSA, Marseille, France

## Abstract

The genome of the facultative anaerobic γ-proteobacterium *Shewanella oneidensis* MR-1 encodes for three terminal oxidases: a *bd*-type quinol oxidase and two heme-copper oxidases, a A-type cytochrome *c* oxidase and a *cbb*₃-type oxidase. In this study, we used a biochemical approach and directly measured oxidase activities coupled to mass-spectrometry analysis to investigate the physiological role of the three terminal oxidases under aerobic and microaerobic conditions. Our data revealed that the *cbb*₃-type oxidase is the major terminal oxidase under aerobic conditions while both *cbb*₃-type and *bd*-type oxidases are involved in respiration at low-O₂ tensions. On the contrary, the low O₂-affinity A-type cytochrome *c* oxidase was not detected in our experimental conditions even under aerobic conditions and would therefore not be required for aerobic respiration in *S. oneidensis* MR-1. In addition, the deduced amino acid sequence suggests that the A-type cytochrome *c* oxidase is a *ccaa*₃-type oxidase since an uncommon extra-C terminal domain contains two *c*-type heme binding motifs. The particularity of the aerobic respiratory pathway and the physiological implication of the presence of a *ccaa*₃-type oxidase in *S. oneidensis* MR-1 are discussed.

**Citation:** Le Laz S, Kpebe A, Bauzan M, Lignon S, Rousset M, et al. (2014) A Biochemical Approach to Study the Role of the Terminal Oxidases in Aerobic Respiration in *Shewanella oneidensis* MR-1. PLoS ONE 9(1): e86343. doi:10.1371/journal.pone.0086343

**Editor:** Alessandro Giuffrè, National Research Council of Italy (CNR), Italy

**Received:** November 13, 2013; **Accepted:** December 11, 2013; **Published:** January 22, 2014

**Copyright:** © 2014 Le Laz et al. This is an open-access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.

**Funding:** This work was supported by research grants from the Agence Nationale de la Recherche (EngineeringH2cyano and Algo-H2 projects), the Centre National de la Recherche Scientifique (CNRS) and the Aix-Marseille Université. The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript.

**Competing Interests:** The authors have declared that no competing interests exist.

* E-mail: mbrugna@imm.cnrs.fr

## Introduction

In the oxygen respiratory systems, electrons of low-redox potential electron donors are transferred through a series of membrane-bound proteins or complexes and finally, the reduction of molecular oxygen to water is catalyzed by enzymes called terminal oxidases. These oxygen reductases are complicated integral membrane multi-subunit complexes grouped into two major superfamilies. Most of them belong to the well-characterized heme-copper oxidases (HCO) superfamily [1]. HCO have been named cytochrome *c* oxidases or quinol oxidases, depending on the nature of their electron donor and are able to pump protons across membrane. Additionally, based on biochemical and structural differences in their catalytic subunits and on phylogenetic analysis, a classification of HCO into three families was proposed [2]: i) type A (mitochondrial-like oxidases or *aa*₃-type), ii) type B (*ba*₃-type oxidases) and iii) type C (*cbb*₃-type oxidases) only detected in bacteria [3]. Cytochrome *bd*-type oxidases, phylogenetically unrelated to the HCO [4], represent the second major oxidases superfamily [5]. Widely distributed among prokaryotes, *bd*-type oxygen reductases function as quinol-oxidases and are bioenergetically less efficient than HCO since they generate a proton motive force by transmembrane charge separation rather than by pumping protons [6–8]. In addition to its role in cell bioenergetics, many studies suggest that cytochrome *bd* oxidase may be implicated in other important physiological functions. In

particular, the enzyme seems to be involved in the bacterial response to a wide variety of stress conditions such as alkalinization of the medium, high temperatures, hydrogen peroxide and nitrosative (NO) stresses [5,9–11]. Cytochrome *bd* oxidase could also play a determinant role in bacterial pathogenicity by protecting bacteria against the NO-mediated host immune response [12].

In contrast to the mitochondrial respiratory systems, most bacteria have branched-respiratory chains terminating in multiple oxidases or use alternative electron acceptors. This enables them to respond to changes in the environment and contributes to their ability to colonize many microoxic and anoxic environments [13–15]. Each oxidase type in bacteria is expected to have a specific affinity for O₂ [16]. Under O₂ limitation, many bacteria induce high O₂-affinity oxidases to respire traces of molecular oxygen. The *cbb*₃-type cytochrome *c* oxidase and the *bd*-type quinol oxidase are generally considered to have a high affinity for oxygen [17,18]. On the other hand, the A-type oxidases have low affinity for oxygen and are predicted to provide the capacity to grow in high-O₂ environments only [16].

*Shewanella* species are facultative anaerobic gram-negative γ-proteobacteria that exhibit extensive respiratory versatility using a broad spectrum of terminal electron acceptors such as fumarate, dimethyl sulfoxide, NO₂⁻, NO₃⁻, Fe(III), As(V), U(VI), Mn(IV), Cr(VI), Tc(VII) elemental sulfur and azo dyes [19–22]. Genome
analysis of *Shewanella oneidensis* MR-1 revealed the presence of genes coding for terminal oxidases: a *bd*-type quinol oxidase and two HCO, a A-type cytochrome *c* oxidase and a *cbb*₃-type oxidase [19]. Latterly, it has been reported that cytochrome *bd* oxidase confers nitrite resistance to *S. oneidensis* MR-1 and plays a significant role in oxygen respiration under microaerobic but not aerobic conditions [23]. Marritt *et al.* [24] hypothesized that, in accordance with their predicted affinity for molecular oxygen, *S. oneidensis* MR-1 A-type and C-type cytochrome *c* oxidases should operate under aerobic and microaerobic conditions respectively. However, the physiological role of the specific terminal oxidases in *S. oneidensis* MR-1 remains to be revealed.

Very recently, while our work was in progress, Zhou *et al.* [25] conducted a relevant study mainly based on gene expression providing a better understanding of the oxygen respiratory metabolism by *S. oneidensis* MR-1. However, because of translational repressor proteins, antisense RNA and posttranslational modifications, mRNA abundance level can be unrelated to active protein abundance level. Thus, many studies suggest that mRNA expression patterns are necessary but insufficient for quantitative description of biological systems [26–29]. For example, comparison of the transcriptome and proteome data for 27 proteins regulated by the ferric uptake regulator in *S. oneidensis* MR-1 revealed that the expression patterns was correlated with the gene expression data for about half of the proteins whereas microarray data and proteome data were not correlated or even inversely correlated for the others proteins [27].

In this work, we used a biochemical approach and directly measured oxidase activities coupled to mass-spectrometry analysis to further characterize the terminal segment of the respiratory chain of *S. oneidensis* MR-1 depending on the growth conditions. Measurements performed on solubilized membranes of wild-type and terminal oxidases deficient strains clearly indicate that the quinol oxidase activity, arising from the *bd*-type oxidase, is extremely higher under microaerobic conditions than under aerobic conditions which is consistent with the fact that *bd*-type oxidases are generally induced under O₂ limited conditions. Surprisingly, *cbb*₃-type oxidase is the only cytochrome *c* oxidase present in *S. oneidensis* MR-1 membranes under both aerobic and microaerobic conditions. The overall results suggest that the A-type cytochrome *c* oxidase, predicted to have a low affinity for O₂, is not present in the bacterium, in our growth conditions, whatever the oxygen tension. A thorough comparison between our results and those obtained by Zhou and colleagues [25] is developed in the discussion section.

# Experimental Procedures

## Bacterial Growth Conditions

The bacterial strains used in this study are listed in Table 1. All the strains were maintained on Luria-Bertani (LB) agar plates containing the appropriate antibiotic. Overnight cultures of the *Escherichia coli* strains CC118 λpir and 1047/pRK2013 used for conjugations were carried out in LB medium at 37°C and 250 rpm. *S. oneidensis* MR-1 and terminal oxidase-deficient strains were grown in LB medium at 30°C either in aerobic conditions or microaerobic conditions. A 30-L BIOSTAT® Cplus-C30-3 fermentor (Sartorius BBI Systems, Germany), controlled by a micro-DCU system and equipped with pH and pO₂ sensors, was used for 30-L batch cultures. Bioreactor was inoculated at OD600 nm~0.2 with 1 L of an overnight flask culture of *S. oneidensis* MR-1 in LB. The pH of the LB medium was regulated at 7.2 with 5 M H₃PO₄. The oxygen partial pressure sensor (InPro® 6820, Mettler-Toledo, Switzerland) was calibrated at 100% in air-

saturated LB at 30°C which corresponds to an estimated oxygen concentration of 8 mg.L⁻¹ (250 µM). For aerobic cultures, the pO₂ was regulated at a minimum of 40% (100 µM O₂) whereas for microaerobic cultures the pO₂ was regulated at a maximum of 6.5% (16 µM O₂). During growth, samples were taken from the cultures at the exponential phase and at the stationary phase, cells were harvested 10 min at 11,400×g and pellets were frozen at -80°C. For growth parameters measurements, strains were inoculated at a OD600 nm~0.05 in 500 mL flask with baffled base containing 100 mL LB and were aerobically grown at 30°C and 250 rpm. Samples were periodically taken in sterile conditions.

## Construction of the Deletion Mutants

*S. oneidensis* deletion mutants were constructed by double crossover events using the suicide plasmid pKNG101 as previously described [30]. Briefly, 500 bp regions upstream and downstream the gene to be deleted were amplified by overlapping PCR with Pfu DNA polymerase (Promega) and *S. oneidensis* MR-1 genomic DNA as template. Primers used for generating PCR products are listed in Table S1. The purified PCR product was cloned into the pGEM®-T Easy vector (Promega) and subcloned into pKNG101 using restriction sites SpeI and Sall. The resulting plasmid was transferred into the plasmid donor strain *E. coli* CC118λpir and mobilized into the appropriate *S. oneidensis* strain by conjugation with the helper strain *E. coli* 1047/pRK2013. Integration of the mutagenesis construct into the chromosome was selected by streptomycin resistance. Clones in which the double recombination events occurred were selected on LB agar plates containing 6% sucrose. Colonies that were resistant to sucrose and susceptible to streptomycin were screened for the deletion event by PCR and confirmed by DNA sequencing. Following this strategy, the SO2364, SO4607 and SO3286 genes encoding respectively the catalytic subunit of the annotated *cbb*₃-type cytochrome *c* oxidase, cytochrome *c* oxidase and cytochrome *d* ubiquinol oxidase were deleted. All the single and double deletion strains are listed in Table 1.

## Preparation and Solubilization of Membranes

About 2 g of cells were resuspended in 20 mM Tris-HCl pH 7.6 supplemented with DNase I and protease inhibitors (Roche) and disrupted twice with an ultra high pressure cell homogenizer (Stansted Fluid Power). Debris and unbroken cells were removed by centrifugation for 10 min at 5,000×g and the supernatant was ultracentrifuged for 1 h at 230,000×g at 4°C to pellet the membranes. Membranes were gently resuspended in 20 mM Tris-HCl pH 7.6 with 5% glycerol to a protein concentration of 10 mg.mL⁻¹ and *n*-dodecyl β-D-maltoside (DDM) was added to a final concentration of 1% (w/v). Solubilization was performed on a rotary tube mixer for 2 h at 4°C. After ultracentrifugation for 1 h at 230,000x g at 4°C, the supernatant containing the DDM-solubilized membranes was retained and protein concentration was determined with a bicinchoninic acid protein assay kit (Sigma) using bovine serum albumin as protein standard.

## Spectral Analysis on Solubilized Membranes

Reduced minus oxidized difference absorbance spectra were recorded on a Lambda 25 UV/VIS spectrophotometer (PerkinElmer) at room temperature. 1 mM potassium ferricyanide was used as the oxidizing agent and a few grains of sodium dithionite was used as the reducing agent.
<table><caption>**Table 1.** Strains used in this study.</caption><thead><tr><th>Strain</th><th>Description</th><th>Reference or source</th></tr></thead><tbody><tr><td><strong><em>E. coli</em></strong></td><td></td><td></td></tr><tr><td>DH5α</td><td>Competent cells for cloning</td><td>New England Biolabs</td></tr><tr><td>1047/pRK2013</td><td>Helper stain for conjugation</td><td>[76]</td></tr><tr><td>CC118λ.pir</td><td>Host strain for pKNG101 replication and donor strain for conjugation</td><td>[77]</td></tr><tr><td><strong><em>S. oneidensis</em></strong><sup>a</sup></td><td></td><td></td></tr><tr><td>MR-1</td><td>Wild-type</td><td>Laboratory collection</td></tr><tr><td>SLL01</td><td>SO2364 deletion mutant/ΔccoN</td><td>This study</td></tr><tr><td>SLL02</td><td>SO3286 deletion mutant/ΔcydA</td><td>This study</td></tr><tr><td>SLL03</td><td>SO4607 deletion mutant/ΔcoxA</td><td>This study</td></tr><tr><td>SLL05</td><td>SO4607 and SO3286 deletion mutant/ΔcoxAΔcydA</td><td>This study</td></tr><tr><td>SLL06</td><td>SO4607 and SO2364 deletion mutant/ΔcoxAΔccoN</td><td>This study</td></tr></tbody></table>

## Cytochrome c Oxidase Activity and O₂ Uptake Measurements

Cytochrome *c* from equine heart was reduced with 75 mM sodium ascorbate in 50 mM sodium phosphate buffer pH 7 with 5% glycerol (v/v). Excess of reductant was removed using a PD10 column (GE Healthcare) and reduced cytochrome *c* was stored at -20°C. Cytochrome *c* oxidase activity was spectrophotometrically measured, on a Lambda 25 UV/VIS spectrophotometer (PerkinElmer) at 550 nm at 30°C in a reaction mixture containing 50 mM sodium phosphate buffer pH 7.4 and 75 µM reduced cytochrome *c*. The reaction was initiated by the addition of solubilized membranes. O₂ reductase activity on solubilized membranes was polarographically measured with a Clark-type oxygen electrode (Hansatech Instruments Oxygraph) in a stirred volume of 1 mL of 50 mM sodium phosphate buffer pH 7.4 using either ubiquinol-1 or NₙNₙN'N'-tetramethyl-p-phenylenediamine (TMPD) as electron donor. Ubiquinol-1 was prepared from ubiquinone-1 (Sigma) by reduction with sodium dithionite and sodium borohydride as previously described [31]. Quinol oxidase activity was determined in the presence of 0.9 mM ubiquinol-1 and 10 mM DTT whereas TMPD oxidase activity assays were carried out with 100 µM TMPD and 10 mM sodium ascorbate. Reactions were initiated by the addition of DDM-solubilized membranes and the consumption of dioxygen was recorded at 30°C. Data were analyzed using the O₂ View software (version 1.02; Hansatech Instruments Ltd.). For inhibition studies, membranes were pre-incubated with 50 µM KCN for 10 min and added to the reaction mixture containing 50 µM KCN.

## Separation by Blue-native gel (BN-gel) Electrophoresis

BN-gels were performed using a mini VE cell electrophoresis apparatus (Amersham Pharmacia) according to the method of Schägger [32,33], as described by Guiral *et al.* [34] except that the intensity was 10 mA/gel and the voltage 250 V. Running times were 3 h. The gels were run with the blue cathode buffer (50 mM Tricine, 7.5 mM imidazole pH 7.0, 0.02% Coomassie Blue G 250 (Serva Blue G)) until the front line had entered into half of the gels and then the cathode buffer was replaced by the one that contained only 0.002% Serva Blue G. The anode buffer was constituted of 25 mM imidazole pH 7.0. The apparent molecular mass of proteins was estimated using the native molecular mass markers NativeMark (Invitrogen).

## In-gel Enzymatic Activities

Cytochrome *c* oxidase activity was revealed on the BN-gel at 30°C using 3,3'-diaminobenzidine and cytochrome *c* from equine heart as previously described [34]. TMPD oxidase activity was revealed on the BN-gel at 30°C using 1 mM TMPD as electron donor in 50 mM sodium phosphate buffer pH 7. Protein bands were cut out from the gel and stored at -20°C before mass spectrometry analysis.

## Protein Identification by In-gel Digestion and Mass Spectrometry

Tryptic digestion experiments: a robotic workstation (Freedom EVO 100, TECAN) was used to perform automated sample preparation, including multiple steps: washes, reduction and alkylation, digestion by trypsin (PROMEGA, proteomics grade with 0.025% ProteasMax) and extraction of tryptic peptides.

Electrospray quadrupole time of flight (ESI-Q-ToF) analyses were performed on a Synapt G1 mass spectrometer (Waters, Manchester) equipped with a NanoLockSpray ion source and coupled to a nano flow UPLC nanoAcquity (Waters, Manchester). Tryptic peptides were dissolved in 3% acetonitrile/0.1% TFA in water and desalted on a C18 nano trap, before on line elution onto a C18 column (BEH 130 C18, 100 µm x 10 cm, 1.7 µm, Waters). Peptides were eluted with a linear gradient from 3% to 50% of mobile phase B (100% acetonitrile/0.1% formic acid) in A (0.1% formic acid in water) for 30 min. The peptides were detected into the mass spectrometer in a positive ion mode using the MS<sup>c</sup> mode. The doubly charged ion of GluFibrinopeptide (785,84 Da) was used as lock mass. Processing of the spectra and protein search were made by Protein Lynx Global Server 2.5.2 (Waters) using the following parameters: *Shewanella oneidensis* database (8 674 entries, with keratins and trypsin), trypsin enzyme, two miscleavages, carbamidomethylation of cysteine as a fixed modification, oxidation of methionine as an optional modification and a mass tolerance set as automatic. Proteins were considered as identified by passing these two filters: minimum 5 fragment consecutive ions from the b/y series matches per peptide and minimum 2 peptide matches per protein.

## Sequence Analysis

Membrane-spanning helices were searched for using the program TMpred [35]. Multiples sequences were aligned using

<footnote>ªAll *S. oneidensis* strains derived from the parental strain *S. oneidensis* MR-1.</footnote>

<footnote>doi:10.1371/journal.pone.0086343.t001</footnote>
the program CLUSTAL W [36] and the HCO phylogenetic tree including various *Shewanella* A-type cytochrome *c* oxidases was built with MEGA5 (http://www.megasoftware.net) [37].

## Results

### Sequence Analysis of the *S. oneidensis* MR-1 Terminal Oxidases

**The cbb₃-type oxidase (Fig. 1A).** The structural subunits of the cbb₃-oxidase are encoded by four genes (SO2364-SO2361) organized in the putative *ccaNOQP* operon. Sequence analysis of *ccaN*, encoding the catalytic subunit of 53.5 kDa, indicates 12 predicted transmembrane helices which harbor the conserved histidines required for ligating the heme *b* and the heme *b*₃/CuB binuclear center [38]. Inspection of the derived amino acid sequences revealed that *ccoO* and *ccoP* encode a monoheme *c*-type cytochrome containing a predicted N-terminal transmembrane helix and a diheme *c*-type cytochrome containing two predicted N-terminal helices respectively. *ccoQ* is predicted to encode a small protein of 6.5 kDa containing a transmembrane helix. The *ccoGHIS* (or *fixGHIS*) gene cluster whose expression is required for the assembly of a functional cbb₃ oxidase is found close to the *ccaNOPQ* operon in most bacteria expressing a cbb₃-type oxidase [39–41]. Interestingly, *ccoG*, encoding a putative polyferredoxin, is a monocistronic predicted gene (SO4737), not genomically adjacent to the *cco* genes. Moreover, an additional gene annotated disulfide bond oxidoreductase *dsbD* (SO2357), is found downstream of *ccoS*. Its predicted amino acid sequence exhibits significant similarity (77%) with the protein VC1435 from *Vibrio cholerae* located at the same position in the *cco* gene cluster and required for assembly of the cbb₃ oxidase [42].

**The bd-type quinol oxidase (Fig. 1B).** The gene cluster SO3286–SO3284 comprises the genes coding for the two subunits CydA (58 kDa) and CydB (42 kDa). The deduced amino acid sequences reveal a high similarity to CydA and CydB from *E. coli* (71% and 66.5% of sequence identity, respectively). CydA is predicted to have nine transmembrane helices and sequence alignment clearly indicates that the residues ligating the low-spin heme *b*₅₅ and the high-spin heme *b*₅₅₅ are conserved in *S. oneidensis* MR-1. CydB is predicted to contain eight membrane-spanning helices like its counterpart in *E. coli* [5]. A third gene

located in the 3' region of *cydAB* and annotated *ybgT* in the *S. oneidensis* MR-1 genome encodes a very small protein of 38 aa predicted to contain a membrane spanning helix. This protein shares 58% sequence identity with the *E. coli* YbgT homologue designated as CydX which has recently been described as a subunit of the cytochrome *bd* oxidase complex required for complex activity [43]. A second gene cluster (SO3780–SO3779) is composed of *cydC* and *cydD*, two genes essential for the assembly of cytochrome *bd*-I in *E. coli* [44–47].

**The A-type cytochrome *c* oxidase (Fig. 1C).** The *cox* gene cluster (SO4606–SO4609) is predicted to encode the three subunits of a third terminal oxidase (*Cox*) annotated *aa*₃-type cytochrome *c* oxidase in the *S. oneidensis* MR-1 genome. The predicted subunit *CoxA* (59 kDa) comprises 12 transmembrane helices and presents 76% of sequence similarity with the catalytic subunit of the *aa*₃ oxidase from *Paracoccus denitrificans*. The six histidines binding the heme-copper binuclear center and the low-spin heme are conserved. The conserved residues forming the D and K-channels and the GHPE260V motif allow to classify *Cox* in A-type oxidases in the A1 subfamily. The predicted subunit II *CoxB* (55 kDa) contains the five conserved ligands of the binuclear CuA center typical of cytochrome *c* oxidase. Sequence alignment with the subunit II of the *aa*₃ cytochrome *c* oxidase from *P. denitrificans* revealed an extension of about 230 amino acid residues in the C-terminal part of *CoxB* containing two *c*-type heme binding consensus sequences. This extra domain likely binds two *c*-type hemes, an uncommon feature among HCO. The predicted subunit III *CoxC* (37 kDa) comprises seven predicted transmembrane helices and does not contain any cofactor.

The *cox* gene cluster includes genes involved in the maturation and assembly of the oxidase. *SenC* is known to bind copper and to be involved in the cytochrome *c* oxidase assembly in yeast or bacteria [48–50]. *ctaA* (SO4613) and *ctaB* (SO4614) encode putative heme *a* and heme *o* synthase respectively. *Surf1* is involved in heme *a* incorporation during bacterial cytochrome *c* oxidase biogenesis [51] and *CtaG* is generally considered as the assembly factor that inserts copper into the CuB center [52,53].

As *Cox* from *S. oneidensis* MR-1 has never been biochemically or spectroscopically investigated, the nature of its heme cofactors is unknown. In subunit I of cytochrome *c* oxidase from *P. denitrificans*, Arginine 54 is in interaction with the formyl group of heme *a*.

<figure><img src="image_2.png" /></figure>

**Figure 1.** Gene clusters encoding the respiratory oxidases of *S. oneidensis* MR-1. Genes encoding the structural subunits are shown in dark grey, genes involved in maturation and assembly of the oxidases are shown in light grey and genes encoding hypothetical proteins (*hyp*) or other proteins are shown in white. Upstream each gene cluster are indicated the positions of the putative Fnr boxes. The number indicated for each gene corresponds to the locus tag. Gene clusters encode the cbb₃ oxidase (A), the bd-type quinol oxidase (B) and the A-type cytochrome *c* oxidase (*Cox*). The scale bar represents 1 kb.
doi:10.1371/journal.pone.0086343.g001
[54,55]. This residue which is not present in oxidases known to incorporate a different type of heme at the "low-spin site" [2] is conserved in *S. oneidensis* MR-1 CoxA. The presence of this residue as well as the *ctaA* and *ctaB* genes required for the heme *a* synthesis in the *cox* gene cluster lead us to believe that an *a*-type heme is incorporated at the "low-spin site" in the subunit I of Cox.

Finally, the presence of the *fnr* (fumarate nitrate regulator) gene directly downstream of the *cco* genes and the presence of predicted Fnr-binding sites in the upstream promoter regions of the different gene clusters suggest a direct regulation by the Fnr transcriptional regulator.

## Spectral Properties of *S. oneidensis* MR-1 Solubilized Membranes

To detect the different terminal oxidases, solubilized membranes from *S. oneidensis* MR-1 aerobically or microaerobically grown were studied by light absorption spectroscopy (Figure 2). Spectra show the presence of cytochrome *c* with peaks at 522 nm and 552 nm and cytochrome *b* with characteristic shoulders at 530 nm and 562 nm. These peaks are due to the contribution from diverse membrane cytochromes. The heme *c/heme *b* ratio is higher (about twofold) in membranes isolated from cells grown under microaerobic conditions (Figure 2 A, B) than under aerobic conditions (Figure 2 C, D). The peak at 632 nm indicates the presence of a heme *d*, cofactor of the *bd*-type quinol oxidase and the trough around 650 nm most likely originates from a stable oxygenated cytochrome *d* species [56]. Furthermore, no *a*-type heme was detected on the spectra of membranes isolated from bacterial cells grown in microaerobic or aerobic conditions. Indeed, membranes isolated from this strain show no signal in the 440 nm-region (data not shown) and in the 600 nm-region. These results indicate the presence of the *bd*-type quinol oxidase in the bacterial membranes and the possible presence of the *cbb3* oxidase in all tested conditions. However, the fact that we did not spectroscopically detect heme *a* leads us to believe that Cox is not present in *S. oneidensis* MR-1 membranes even in aerobic conditions. Nevertheless, we cannot rule out the possibility that Cox does not contain any heme *a* as a cofactor. The possible lack of the A-type cytochrome *c* oxidase in the membranes must be confirmed with other techniques.

## TMPD, Cytochrome *c* and Quinol Oxidase Activities in Solubilized Membranes of *S. oneidensis* MR-1

To further investigate the role of the different terminal oxidases in aerobic respiration in *S. oneidensis* MR-1, catalytic activity of cytochrome *c* oxidases and *bd*-type quinol oxidase were spectro-photometrically and polarographically measured with solubilized membranes prepared from cells grown in LB medium under aerobic and microaerobic conditions (Table 2).

Cytochrome *c* oxidase activity was assessed by following oxidation of horse heart reduced cytochrome *c* by spectrophotometry and by monitoring O₂ consumption in the presence of TMPD, an artificial electron donor capable of reducing cytochrome *c*. The values obtained by both methods indicate that cytochrome *c* oxidase activity is relatively constant, regardless of the tested culture conditions and growth phase.

Quinol oxidase activity of the cytochrome *bd*-type oxidase was determined by measuring O₂ uptake with ubiquinol-1 as the electron donor. While the values obtained with membranes from cells grown in aerobic conditions were low and near to the detection limit, quinol-dependent O₂ consumption was significantly higher with membranes from microaerobic cultures with a 20–40 fold increase depending on the growth phase. In addition,

<figure><img src="image_3.png" /><figcaption>**Figure 2.** Reduced minus oxidized difference absorbance spectra of *S. oneidensis* MR-1 solubilized membranes. Solubilized membranes in 20 mM Tris-HCl pH 7.6, 5% glycerol, 1% DDM were oxidized using 1 mM potassium ferricyanide and reduced by adding few grains of sodium dithionite. Membranes were isolated from strain MR-1 grown under microaerobic conditions (spectra A and B) or aerobic conditions (spectra C and D). Bacterial cells were harvested in the exponential (spectra A and C) or stationary phase (spectra B and D) of growth. The concentration of proteins was 4.8 mg.mL⁻¹. The spectra were recorded at room temperature. The vertical bar indicates the absorption scale.

*doi:10.1371/journal.pone.0086343.g002*</figcaption></figure>

the oxygen consumption in the presence of ubiquinol-1 was lower with membranes prepared from cells harvested during the exponential phase of growth compared to the stationary phase under both aerobic and microaerobic conditions. Our data are consistent with the fact that *bd*-type oxygen reductases are generally expressed under O₂-limited conditions i.e. in microaerobic conditions or when cultures enter the stationary phase, the decrease in oxygen tension resulting from the increase in cell density [14,57].

Furthermore, the addition of 50 µM KCN strongly inhibited TMPD-dependent O₂ consumption and cytochrome *c* oxidase activity (90–100%) whereas the quinol-dependent O₂ uptake was slightly affected with a 10–15% inhibition (data not shown). This observation is in accord with the fact that a specific feature of cytochrome *bd*-type oxidase is its much lower sensitivity to cyanide than heme-copper oxygen reductases [5].

## Identification of Cytochrome *c* Oxidases by Mass Spectrometry

In order to identify enzymes responsible for the measured cytochrome *c* oxidase activity, solubilized *S. oneidensis* MR-1 membrane proteins were separated on a BN-gel electrophoresis followed by an in-gel detection of the cytochrome *c* oxidase activity (Fig. 3A). For this purpose, *S. oneidensis* MR-1 cells were grown aerobically or microaerobically and harvested during the exponential or the stationary phase of growth to determine the influence of the oxygen level and the growth phase on the respiratory cytochrome *c* oxidase content in the membranes. BN-gel of the solubilized membranes shows one major brown band of cytochrome *c* oxidase activity, of molecular mass around 240 kDa (Fig. 3A, band I). Another faint band of activity could be detected, of molecular mass around 450 kDa (Fig. 3A, band II). The pattern
<table><caption>**Table 2.** TMPD, cytochrome *c* and quinol oxidase activities in solubilized membranes of *S. oneidensis* MR-1.</caption><thead><tr><th rowspan="2"></th><th colspan="2">Aerobic conditions</th><th colspan="2">Microaerobic conditions</th></tr><tr><th>EP<sup>a</sup></th><th>SP<sup>b</sup></th><th>EP</th><th>SP</th></tr></thead><tbody><tr><td>TMPD oxidase activity (nmol O<sub>2</sub>.min<sup>-1</sup>.mg protein<sup>-1</sup>)</td><td>90+/-6</td><td>95+/-6</td><td>94+/-4</td><td>84+/-3</td></tr><tr><td>Cyt *c* oxidase activity (nmol cyt c.min<sup>-1</sup>.mg protein<sup>-1</sup>)<sup>c</sup></td><td>6.2+/-0.2</td><td>8.6+/-0.5</td><td>7.3+/-0.5</td><td>5.2+/-0.6</td></tr><tr><td>Ubiquinol-1 oxidase activity (nmol O<sub>2</sub>.min<sup>-1</sup>.mg protein<sup>-1</sup>)</td><td>9+/-1</td><td>20+/-1</td><td>364+/-18</td><td>456+/-17</td></tr></tbody></table>

<sup>a</sup>EP: exponential phase of growth.

<sup>b</sup>SP: stationary phase of growth.

<sup>c</sup>cyt: cytochrome.

Listed values are averages of at least three separate experiments.

doi:10.1371/journal.pone.0086343.t002

of cytochrome *c* oxidase activity bands was identical whatever the oxygen tension in the medium and the growth phase. When TMPD was used instead of horse heart cytochrome *c* to reveal cytochrome *c* oxidase activity, the same pattern was also observed (data not shown). With both electron donors, no band was revealed when the in-gel detection was performed in the presence of 50 µM potassium cyanide (data not shown).

To identify the cytochrome *c* oxidases, the protein band I and band II were directly cut out from the gel, for the four growth conditions, digested with trypsin and analyzed by ESI-Q-ToF mass spectrometry (Fig. 3B). Band I and band II were found to contain, among others proteins, CcoO and CcoP, the monohemic and dihemic cytochrome *c* subunits of the *cbb*<sub>3</sub> oxidase respectively. In addition, the CydA subunit of the *bd*-type quinol oxidase was also identified in band I and band II indicating that the *cbb*<sub>3</sub> oxidase and the *bd*-type quinol oxidase co-migrate in the BN gel. However, none of the three subunits of the Cox enzyme was identified in solubilized membranes of *S. oneidensis* MR-1 grown microaerobically or aerobically. These results strongly suggest that, against all expectation, the high O<sub>2</sub>-affinity *cbb*<sub>3</sub> oxidase is responsible for the cytochrome *c* oxidase activity measured in the

membranes of microaerobically but also aerobically grown *S. oneidensis* MR-1 whereas the A-type cytochrome *c* oxidase Cox remains undetectable even under aerobic conditions.

## Growth Parameters and TMPD Oxidase Activity of *S. oneidensis* MR-1 Wild Type and Oxidase Deletion Mutants

The first part of the study suggests that the *bd*-type and the *cbb*<sub>3</sub>-type high-affinity oxidases are important under aerobic conditions while the A-type cytochrome *c* oxidase would be surprisingly dispensable despite its low affinity for O<sub>2</sub>. To confirm these data and determine the physiological significance of the different oxygen reductases under aerobic conditions, deletion mutants were constructed and characterized during the exponential phase of growth under aerobic conditions by measuring growth parameters and TMPD oxidase activities (Table 3). While the strain lacking the *cbb*<sub>3</sub>-type terminal oxidase (SLL01) grew significantly more slowly and reached the stationary phase at lower cell density compared to the wild type, the single oxidase mutants lacking the *bd*-type terminal oxidase (SLL02) or the A-type cytochrome *c* oxidase Cox (SLL03) grew as well as the wild

<figure><img src="image_4.png" /><figcaption>**Figure 3.** Identification of cytochrome *c* oxidase(s) in *S. oneidensis* MR-1 membranes under aerobic and microaerobic conditions. A. In-gel detection of cytochrome *c* oxidase activity in solubilized membranes on BN-gel. Membrane proteins were prepared from *S. oneidensis* MR-1 cells grown aerobically (lanes 1 and 2) or microaerobically (lanes 3 and 4) and harvested during the exponential (lanes 1 and 3) or the stationary (lanes 2 and 4) phase of growth. 130 µg of total proteins were loaded on a 5–15% polyacrylamide gel. The major brown band (band I) and the faint brown band (band II) of activity are indicated by arrows. B. Terminal oxidases subunits identified by ESI-Q-ToF mass spectrometry. Data correspond to identified proteins in bands I and II from exponentially grown cells under aerobic conditions (panel A, lane 1). Table heading: Band, roman figures refer to the protein bands from the BN gel shown in panel A; Protein name, name in NCBI database; NCBI entry, accession number; Peptides, number of unique peptides detected; Coverage, protein sequence coverage by the matching peptides (in %). Cyt: cytochrome.
doi:10.1371/journal.pone.0086343.g003</figcaption></figure>
<table><caption>**Table 3.** Growth parameters and TMPD-dependent oxidase activity of oxidase deletion mutants cultivated in aerobic conditions.</caption><thead><tr><th>Strain</th><th>Relevant genotype</th><th>Doubling time<sup>a</sup></th><th>Relative yield<sup>b</sup></th><th>TMPD oxidase activity<sup>c</sup></th></tr></thead><tbody><tr><td>MR-1</td><td>wild-type</td><td>54.2+/-0.2</td><td>1.00</td><td>90+/-6</td></tr><tr><td>SLL01</td><td>ΔccoN</td><td>78.0+/-0.7</td><td>0.77</td><td>0</td></tr><tr><td>SLL02</td><td>ΔcydA</td><td>53.4+/-2.1</td><td>1.00</td><td>91+/-9</td></tr><tr><td>SLL03</td><td>ΔcoxA</td><td>59.0+/-1.1</td><td>0.96</td><td>97+/-6</td></tr><tr><td>SLL05</td><td>ΔcydAΔcoxA</td><td>56.5+/-2.7</td><td>0.97</td><td>101+/-9</td></tr><tr><td>SLL06</td><td>ΔccoNΔcoxA</td><td>79.9+/-1.2</td><td>0.78</td><td>0</td></tr></tbody></table>

type. Firstly, the results indicate that none of the terminal oxidases are essential for aerobic growth in our conditions. Secondly, the reduced growth rate of SLL01 suggests that the *cbb*₃-type terminal oxidase is the predominant oxidase in exponentially growing *S. oneidensis* MR-1 cells under aerobic conditions.

The strains lacking two terminal oxidases were also characterized. Growth parameters of the strain containing only the *cbb*₃-type oxidase (SLL05) did not really differ from those of the wild-type indicating that the *cbb*₃-type oxidase is sufficient to support maximum growth rate in aerobic conditions. Conversely, doubling time of the strain containing only the *bd*-type oxidase (SLL06) was significantly longer and the growth yield lower than those of the wild-type strain, with values similar to those of the single oxidase mutants lacking *cbb*₃-type oxidase (SLL01). Interestingly, it was not possible, despite numerous efforts, to generate a double oxidase mutant lacking both *bd*-type and *cbb*₃-type oxidases suggesting that this mutant was not viable under aerobic conditions.

Taken together, these results demonstrate that the *cbb*₃-type or the *bd*-type oxidase is required for aerobic growth of *S. oneidensis* MR-1 and that the *cbb*₃-type terminal oxidase is the major oxidase in exponentially growing cells. These findings are consistent with the results obtained in this work with the wild type and also suggest that the low O₂-affinity cytochrome *c* oxidase Cox is unexpectedly not involved in aerobic respiration in our experimental conditions.

Furthermore, TMPD-dependent oxidase activity was determined with *S. oneidensis* MR-1 wild type strain and oxidase deletion mutants cultivated under aerobic conditions. First, the strain SLL06 exhibited no activity indicating that TMPD was not oxidized by the *bd*-type oxidase. This result contrasts with the cytochrome *bd* oxidase from *E. coli* which can oxidize TMPD with a quite high turnover number [58]. However, several previously reported cytochrome *bd*-oxidases have been found to exhibit low TMPD oxidase activity as in *Azotobacter vinelandii* and *Geobacillus thermodenitrificans* K1041 (formerly *Bacillus stearothermophilus* K1041), or no TMPD oxidase activity as in *Bacillus firmus* OF4 [59–62]. Besides, TMPD oxidase activity was not affected in the single oxidase mutants lacking *cydA* (SLL02) or *coxA* (SLL03) compared to the wild-type strain whereas no activity was measured with the strain lacking *ccoN* (SLL01). Finally, TMPD oxidase activity did not differ from that of the wild type in the strain containing only the *cbb*₃-type oxidase (SLL05). These data indicate that TMPD oxidase activity arose from the *cbb*₃-type oxidase only suggesting that the *cbb*₃-type oxidase is the only HCO present in the aerobic membranes of *S. oneidensis* MR-1. Moreover, this could also explain the low values obtained from measurements of cytochrome *c* oxidase activity in Table 2 since a previous study revealed that

horse heart cytochrome *c* was a poor substrate for the *cbb*₃-type oxidase in *Vibrio cholerae*, with a rate of oxygen reduction 10–20 fold lower than with the diheme cytochrome *c*₄ identified as a natural electron donor [63].

## Spectral Properties of *S. oneidensis* MR-1 Oxidase Deletion Mutants Solubilized Membranes

Solubilized membranes from aerobically grown *S. oneidensis* MR-1 and strains lacking terminal oxidases (single and double mutants) were studied by light absorption spectroscopy (Fig. 4). Dithionite-reduced minus ferricyanide-oxidized spectra demonstrate the presence of *c*-type and *b*-type cytochromes (peaks at 522, 552, and shoulders at 530 and 562 nm respectively) in membranes from all strains. The study of the mutants confirmed that the peak at 632 nm and the trough at 650 nm arise from the heme *d* of the *bd*-quinol oxidase since they are absent on spectra from SLL02 and SLL05, both strains lacking this quinol oxidase. Additionally, spectral signals for *a*-type heme (around 440 and 600 nm) were not detected in any of the strains confirming the absence of Cox in the membrane of *S. oneidensis* MR-1 in aerobic conditions.

## Identification of Cytochrome *c* Oxidases in *S. oneidensis* MR-1 Oxidase Mutants Membranes by Mass Spectrometry

Membrane proteins from *S. oneidensis* MR-1 and oxidase mutant strains grown under aerobic conditions were separated on BN-gel electrophoresis and the cytochrome *c* oxidase activity was revealed (Fig. 5). One band of activity, of molecular mass around 240 kDa, was observed in all strains except for the strain lacking the *cbb*₃ oxidase (SLL01) and the strain lacking both the *cbb*₃ and the Cox oxidases (SLL06). If the BN-gel is incubated a longer time for the detection of the cytochrome *c* oxidase activity, a faint band of activity is detected, of molecular mass around 450 kDa (not shown) in all strains except for SLL01 and SLL06. These results strongly suggest that the *cbb*₃-oxidase is responsible for the bands of activity in the solubilized membranes of the bacterium. Furthermore, the analysis of these bands of activity by ESI-Q-ToF mass spectrometry (data not shown) confirmed that the cytochrome *c* oxidase responsible for in-gel activity is the *cbb*₃ oxidase since both CcoO and CcoP, the monohemic and dihemic cytochrome *c* subunits of the *cbb*₃ oxidase respectively are identified. In contrast, none of the *A*-type Cox subunits was identified. Collectively, these results undoubtedly demonstrate that under aerobic conditions, the cytochrome *c* oxidase activity detected in the membranes from *S. oneidensis* MR-1 arises from the high-affinity *cbb*₃ oxidase and

<footnote><sup>a</sup>Expressed in min.</footnote>

<footnote><sup>b</sup>Relative yield is defined as the highest optical density in early stationary phase in mutant strains cultures relative to that in the wild-type strain cultures.</footnote>

<footnote><sup>c</sup>Expressed in nmol O₂·min⁻¹·mg protein⁻¹.</footnote>

<footnote>Data were obtained from triplicate separated experiments.</footnote>

<footnote>doi:10.1371/journal.pone.0086343.t003</footnote>
<figure><img src="image_5.png" /><figcaption>**Figure 4.** Reduced minus oxidized difference absorbance spectra 0 of membranes from aerobically grown *S. oneidensis* MR-1 oxidase mutants. Solubilized membranes in 20 mM Tris-HCl pH 7.6, 5% glycerol, 1% DDM were oxidized using 1 mM potassium ferricyanide and reduced by adding few grains of sodium dithionite. Membranes were isolated from aerobically grown cells harvested in the exponential phase of growth. The concentration of proteins was 11 mg·mL⁻¹. The spectra were recorded at room temperature. The vertical bar indicates the absorption scale.
doi:10.1371/journal.pone.0086343.g004</figcaption></figure>

that, on the contrary, the low-affinity A-type cytochrome *c* oxidase is surprisingly not required for aerobic respiration.

<figure><img src="image_6.png" /><figcaption>**Figure 5.** Detection of cytochrome *c* oxidase activity on BN-gel in aerobically grown *S. oneidensis* MR-1 oxidase mutants. Membrane proteins were prepared from wild-type or oxidase mutants *S. oneidensis* MR-1 cells grown aerobically and harvested during the exponential phase of growth. 130 µg of total proteins were loaded on a 5–15% polyacrylamide gel. The band of activity is indicated by an arrow. Lane 1: *S. oneidensis* MR-1. Lane 2: strain SLL01. Lane 3: strain SLL02. Lane 4: strain SLL03. Lane 5: strain SLL05. Lane 6: strain SLL06.
doi:10.1371/journal.pone.0086343.g005</figcaption></figure>

Discussion

In this report, we investigated the aerobic respiratory pathway
in *S. oneidensis* MR-1 to determine the physiological role of the
three terminal oxidases in low and in high-O₂ environments. Our
results suggest that the high-affinity *bd*-type oxidase is weakly
expressed in aerobic conditions and significantly induced under
microaerobic conditions which is consistent with previous studies
indicating that *bd* oxidases are typically used by bacteria for
aerobic respiration under O₂-limited conditions [5,64,65]. Taken
together, our data also reveal that the *cbb₃*-type oxidase is the most
important terminal oxidase under aerobic conditions and has a
significant role under microaerobic conditions whereas the low
affinity A-type cytochrome *c* oxidase Cox was not detected in the
tested conditions even under aerobic conditions. In addition,
characterization of the terminal oxidase deletion mutants confirms
that the *cbb₃*-type or the *bd*-type oxidase is required for aerobic
growth and that Cox is unexpectedly not involved in aerobic
respiration in these experimental conditions.

Our results are overall consistent with those recently reported
by Zhou *et al*. [25] who revealed, by measuring the mRNA
abundance and the promoter activity in *S. oneidensis* MR-1, that
*cbb*₃-type oxidase and *bd*-type oxidase are important under
microaerobic conditions and the *cbb*₃-type is the major oxidase
under aerobic conditions while Cox has no physiological
significance. However, it is noteworthy that, in their study, the
expression level of the *cco* genes encoding the *cbb*₃ oxidase was
lower under microaerobic conditions than under aerobic condi-
tions and directly proportional to the O₂ level in aerobic cultures,
suggesting that expression of *cco* is favored in O₂-rich environ-
ments. On the contrary, we clearly present evidence that the
oxidase activity arising from the *cbb*₃ oxidase was constant under
the tested conditions and thus did not respond to growth
conditions. Furthermore, if Zhou *et al*. [25] found that the *cyd*
operon was slightly induced at low O₂ concentration compared to
high O₂ conditions (<5-fold increase in mRNA abundance), the
difference observed between aerobic and microaerobic conditions
was much more significant at the protein level with a 40-fold
increase observed by measuring quinol-dependent oxidase activity.
Probably due to translational and/or post-translational regulation
mechanisms, these discrepancies emphasize the importance of
biochemical analysis at the protein level to understand the global
metabolic and regulatory processes in addition to studies of
expression at a transcriptional level.

These findings are intriguing since *cbb3* oxidases, given their
high affinity for O₂, are usually repressed at high oxygen tension in
many bacteria [66–69]. Furthermore, in organisms carrying an
*aa3*-type and a *cbb3*-type oxidase, the low affinity oxidase plays a
dominant role under high oxygen conditions whereas the *cbb3*
oxidase is induced only under low O₂ concentrations [67]. In *S.
oneidensis* MR-1, the uncommon expression pattern of the
cytochrome *c* oxidases is reminiscent of the one found in
*Pseudomonas aeruginosa*. In this bacterium, the *cox* genes encoding
an *aa3*-type oxidase are expressed at a very low level under O₂-rich
growth conditions in LB-medium while the expression level of the
genes encoding the *cbb3*-1 oxidase was high in the same conditions.

Interestingly, the *cox* genes were found to be significantly
induced in *P. aeruginosa* under starvation of iron, carbon or
nitrogen indicating that the *aa3*-type oxidase might be involved in
energy conservation under low nutrient conditions [70,71]. The
induction of the *cox* genes under nutrient starvation or other
environmental stress conditions should thus be considered in *S.
oneidensis* MR-1. Although a functional loss of the A-type Cox
oxidase in *S. oneidensis* MR-1, as proposed by Zhou *et al*. [25],
cannot be excluded, the presence of the *cox* genes as well as of genes involved in the maturation and assembly of the oxidase in all genome-sequenced *Shewanella* strains leads us to believe that the Cox oxidase has a significant role in some not yet defined conditions. Moreover, in the phylogenetic tree built with the sequences of the catalytic subunits of A, B and C-type oxidases from various Bacteria and Archaea, the A-type Cox from different *Shewanella* species, among them the one from *S. oneidensis* MR-1, clustered with bacterial A-type cytochrome *c* oxidases (data not shown). For *Shewanella* genes, the observed branch lengths did not indicate a high evolutionary rate characteristic of pseudogenes, which is again not in favour of the hypothesis of a functional loss of the A-type Cox.

In addition to the particularity of the membrane terminal oxidase content in *S. oneidensis* MR-1, examination of the derived amino acid sequences emphasized the singularity of the Cox oxidase. Indeed, if sequence analysis allowed to classify Cox in the A1 subfamily and to predict the incorporation of *a*-type heme at least at the low-spin site of CoxA, an uncommon extra C-terminal domain carrying two *c*-type heme binding consensus sequences was identified in the predicted subunit CoxB. Subunit II with an extra domain carrying only one *c*-type heme can be found in bacterial terminal oxidases such as in the *caa*<sub>3</sub> oxidase from *Thermus thermophilus* or *Bacillus subtilis* [72,73]. Triguingly, the C-terminal extension of CoxB from *S. oneidensis* MR-1 likely binds two *c*-type hemes suggesting that this enzyme is a *caa*<sub>3</sub>-type HCO. Such an unusual feature was so far described only in *Desulfovibrio* species [74,75] but we discover that the two *c*-type heme binding motifs were also conserved in some γ-proteobacteria such as *Shewanella* (except in *S. violacea* and *S. benthica* KT99) and some species of the

genus *Psychromonas*, *Colwellia* and *Methylosarcina* (data not shown). As in *Desulfovibrio* species, the two heme *c* domains in *S. oneidensis* MR-1 are highly similar to each other (61% sequence identity) suggesting that the diheme domain evolved from a gene duplication event [74].

A study is currently underway to determine experimental conditions under which the *ccaa*<sub>3</sub>-type oxidase Cox would be expressed. It should provide information about the physiological significance of this oxidase in *S. oneidensis* MR-1.

**Supporting Information**

<caption>**Table S1** Oligonucleotides used in this study.
(DOCX)</caption>

**Acknowledgments**

We thank Cécile JOURLIN-CASTELLI, Joshua ARMITANO and Vincent MEJEAN (CNRS, LCB, Marseilles, France) for expert technical assistance in bacterial conjugation experiments, for kindly providing the *Shewanella oneidensis* MR-1 strain and for helpful discussions. We acknowledge Wolfgang NITSCHKE (CNRS, BIP, Marseilles, France) for his help in phylogenetic analysis and Marianne GUIRAL for critical reading of the manuscript. The Proteomic platform is a member of Marseille Protéomique (http://map.univmed.fr/index.html) labeled IBiSA.

**Author Contributions**

Conceived and designed the experiments: SLL AK MR M. Brugna.
Performed the experiments: SLL AK M. Bauzan SL M. Brugna. Analyzed
the data: SLL AK M. Brugna. Wrote the paper: SLL AK M. Brugna.

**References**

1. Garcia-Horsman JA, Barquera B, Rumbley J, Ma J, Gennis RB (1994) The superfamily of heme-copper respiratory oxidases. J Bacteriol 176: 5587-5600.

2. Pereira MM, Santana M, Teixeira M (2001) A novel scenario for the evolution of haem-copper oxygen reductases. Biochim Biophys Acta 1505: 185-208.

3. Brochier-Armanet C, Talla E, Gribaldo S (2009) The multiple evolutionary histories of dioxygen reductases: Implications for the origin and evolution of aerobic respiration. Mol Biol Evol 26: 285-297.

4. Junemann S (1997) Cytochrome *bd* terminal oxidase. Biochim Biophys Acta 1321: 107-127.

5. Borisov VB, Gennis RB, Hemp J, Verkhovsky MI (2011) The cytochrome *bd* respiratory oxygen reductases. Biochim Biophys Acta 1807: 1398-1413.

6. Puustinen A, Finel M, Haltia T, Gennis RB, Wikstrom M (1991) Properties of the two terminal oxidases of *Escherichia coli*. Biochemistry 30: 3936-3942.

7. Belevich I, Borisov VB, Zhang J, Yang K, Konstantinov AA, et al. (2005) Time-resolved electrometric and optical studies on cytochrome *bd* suggest a mechanism of electron-proton coupling in the di-heme active site. Proc Natl Acad Sci U S A 102: 3657-3662.

8. Borisov VB, Murali R, Verkhovskaya ML, Bloch DA, Han H, et al. (2011) Aerobic respiratory chain of *Escherichia coli* is not allowed to work in fully uncoupled mode. Proc Natl Acad Sci U S A 108: 17320-17324.

9. Borisov VB, Forte E, Davletshin A, Mastronicola D, Sarti P, et al. (2013) Cytochrome *bd* oxidase from *Escherichia coli* displays high catalase activity: an additional defense against oxidative stress. FEBS Lett 587: 2214-2218.

10. Borisov VB, Forte E, Konstantinov AA, Poole RK, Sarti P, et al. (2004) Interaction of the bacterial terminal oxidase cytochrome *bd* with nitric oxide. FEBS Lett 576: 201-204.

11. Borisov VB, Forte E, Sarti P, Brunori M, Konstantinov AA, et al. (2007) Redox control of fast ligand dissociation from *Escherichia coli* cytochrome *bd*. Biochem Biophys Res Commun 355: 97-102.

12. Giuffre A, Borisov VB, Mastronicola D, Sarti P, Forte E (2012) Cytochrome *bd* oxidase and nitric oxide: from reaction mechanisms to bacterial physiology. FEBS Lett 586: 622-629.

13. Richardson DJ (2000) Bacterial respiration: a flexible process for a changing environment. Microbiology 146: 551-571.

14. Poole RK, Cook GM (2000) Redundancy of aerobic respiratory chains in bacteria? Routes, reasons and regulation. Adv Microb Physiol 43: 165-224.

15. Bueno E, Mesa S, Bedmar EJ, Richardson DJ, Delgado MJ (2012) Bacterial adaptation of respiration from oxic to microoxic and anoxic conditions: redox control. Antioxid Redox Signal 16: 819-852.

16. Morris RL, Schmidt TM (2013) Shallow breathing: bacterial life at low O₂. Nat Rev Microbiol 11: 205-212.

17. Belevich I, Borisov VB, Bloch DA, Konstantinov AA, Verkhovsky MI (2007) Cytochrome *bd* from *Azotobacter vinelandii*: evidence for high-affinity oxygen binding. Biochemistry 46: 11177-11184.

18. Belevich I, Borisov VB, Konstantinov AA, Verkhovsky MI (2005) Oxygenated complex of cytochrome *bd* from *Escherichia coli*: stability and photolability. FEBS Lett 579: 4567-4570.

19. Heidelberg JF, Paulsen IT, Nelson KE, Gaidos EJ, Nelson WC, et al. (2002) Genome sequence of the dissimilatory metal ion-reducing bacterium *Shewanella oneidensis*. Nat Biotechnol 20: 1118-1123.

20. Liu C, Gorby YA, Zachara JM, Fredrickson JK, Brown CF (2002) Reduction kinetics of Fe(III), Co(III), U(VI), Cr(VI), and Tc(VII) in cultures of dissimilatory metal-reducing bacteria. Biotechnol Bioeng 80: 637-649.

21. Gralnick JA, Vali H, Lies DP, Newman DK (2006) Extracellular respiration of dimethyl sulfoxide by *Shewanella oneidensis* strain MR-1. Proc Natl Acad Sci U S A 103: 4669-4674.

22. Le Laz S, Kpebe A, Lorquin J, Brugna M, Rousset M (2013) H₂-dependent azoreduction by *Shewanella oneidensis* MR-1: involvement of secreted flavins and both [Ni-Fe] and [Fe-Fe] hydrogenases. Appl Microbiol Biotechnol. In Press. 10.1007/s00253-013-5208-z.

23. Fu H, Chen H, Wang J, Zhou G, Zhang H, et al. (2013) Crp-dependent cytochrome *bd* oxidase confers nitrite resistance to *Shewanella oneidensis*. Environ Microbiol 15: 2198-2212.

24. Marritt SJ, McMillan DG, Shi L, Fredrickson JK, Zachara JM, et al. (2012) The roles of CymA in support of the respiratory flexibility of *Shewanella oneidensis* MR-1. Biochem Soc Trans 40: 1217-1221.

25. Zhou G, Yin J, Chen H, Hua Y, Sun L, et al. (2013) Combined effect of loss of the *caa*<sub>3</sub> oxidase and Crp regulation drives Shewanella to thrive in redox-stratified environments. ISME J 7: 1752-1763.

26. Gygi SP, Rochon Y, Franza BR, Aebersold R (1999) Correlation between protein and mRNA abundance in yeast. Mol Cell Biol 19: 1720-1730.

27. Wan XF, Verberkmoes NC, McCue LA, Stanek D, Connelly H, et al. (2004) Transcriptomic and proteomic characterization of the Fur modulon in the metal-reducing bacterium *Shewanella oneidensis*. J Bacteriol 186: 8385-8400.

28. Nie L, Wu G, Zhang W (2006) Correlation of mRNA expression and protein abundance affected by multiple sequence features related to translational efficiency in *Desulfovibrio vulgaris*: a quantitative analysis. Genetics 174: 2229-2243.

29. Nie L, Wu G, Zhang W (2006) Correlation between mRNA and protein abundance in *Desulfovibrio vulgaris*: a multiple regression to identify sources of variations. Biochem Biophys Res Commun 339: 603-610.
30. Kaniga K, Delor I, Cornelis GR (1991) A wide-host-range suicide vector for improving reverse genetics in gram-negative bacteria: inactivation of the blaA gene of *Yersinia enterocolitica*. *Gene* 109: 137–141.

31. Rieske JS (1967) Preparation and properties of reduced coenzyme Q-cytochrome *c* reductase (complex III of the respiratory chain). Methods Enzymol 10: 239–245.

32. Schagger H, von Jagow G (1991) Blue native electrophoresis for isolation of membrane protein complexes in enzymatically active form. *Anal Biochem* 199: 223–231.

33. Schagger H, Cramer WA, von Jagow G (1994) Analysis of molecular masses and oligomeric states of protein complexes by blue native electrophoresis and isolation of membrane protein complexes by two-dimensional native electrophoresis. *Anal Biochem* 217: 220–230.

34. Guiral M, Prunetti L, Lignon S, Lebrun R, Moinier D, et al. (2009) New insights into the respiratory chains of the chemolithoautotrophic and hyperthermophilic bacterium *Aquifex aeolicus*. *J Proteome Res* 8: 1717–1730.

35. Claros MG, von Heijne G (1994) TopPred II: an improved software for membrane protein structure predictions. *Comput Appl Biosci* 10: 685–686.

36. Thompson JD, Higgins DG, Gibson TJ (1994) CLUSTAL W: improving the sensitivity of progressive multiple sequence alignment through sequence weighting, position-specific gap penalties and weight matrix choice. Nucleic Acids Res 22: 4673–4680.

37. Tamura K, Peterson D, Peterson N, Stecher G, Nei M, et al. (2011) MEGA5: molecular evolutionary genetics analysis using maximum likelihood, evolutionary distance, and maximum parsimony methods. *Mol Biol Evol* 28: 2731–2739.

38. Buschmann S, Warkentin E, Xie H, Langer JD, Ermler U, et al. (2010) The structure of *cbb3* cytochrome oxidase provides insights into proton pumping. *Science* 329: 327–330.

39. Preisig O, Zufferey R, Hennecke H (1996) The *Bradyrhizobium japonicum fixGH5* genes are required for the formation of the high-affinity *cbb3*-type cytochrome oxidase. *Arch Microbiol* 165: 297–305.

40. Koch HG, Winterstein C, Saribas AS, Alben JO, Daldal F (2000) Roles of the *ccoGH5* gene products in the biogenesis of the *cbb3*-type cytochrome *c* oxidase. *J Mol Biol* 297: 49–65.

41. Kulajta C, Thumfart JO, Haid S, Daldal F, Koch HG (2006) Multi-step assembly pathway of the *cbb3*-type cytochrome *c* oxidase complex. *J Mol Biol* 355: 989–1004.

42. Braun M, Thony-Meyer L (2005) Cytochrome *c* maturation and the physiological role of *c*-type cytochromes in *Vibrio cholerae*. *J Bacteriol* 187: 5996–6004.

43. VanOrsdel CE, Bhatt S, Allen RJ, Brenner EP, Hobson JJ, et al. (2013) The *Escherichia coli* *CydX* protein is a member of the *CydAB* cytochrome *bd* oxidase complex and is required for cytochrome *bd* oxidase activity. *J Bacteriol* 195: 3640–3650.

44. Georgiou CD, Fang H, Gennis RB (1987) Identification of the *cydC* locus required for expression of the functional form of the cytochrome *d* terminal oxidase complex in *Escherichia coli*. *J Bacteriol* 169: 2107–2112.

45. Poole RK, Hatch L, Cleeter MW, Gibson F, Cox GB, et al. (1993) Cytochrome *bd* biosynthesis in *Escherichia coli*: the sequences of the *cydC* and *cydD* genes suggest that they encode the components of an ABC membrane transporter. *Mol Microbiol* 10: 421–430.

46. Bebbington KJ, Williams HD (1993) Investigation of the role of the *cydD* gene product in production of a functional cytochrome *d* oxidase in *Escherichia coli*. *FEMS Microbiol Lett* 112: 19–24.

47. Poole RK, Gibson F, Wu G (1994) The *cydD* gene product, component of a heterodimeric ABC transporter, is required for assembly of periplasmic cytochrome *c* and of cytochrome *bd* in *Escherichia coli*. *FEMS Microbiol Lett* 117: 217–223.

48. Krummeck G, Rodel G (1990) Yeast *SCO1* protein is required for a post-translational step in the accumulation of mitochondrial cytochrome *c* oxidase subunits I and II. *Curr Genet* 18: 13–15.

49. Nittis T, George GN, Winge DR (2001) Yeast *Sco1*, a protein essential for cytochrome *c* oxidase function is a *Cup*-binding protein. *J Biol Chem* 276: 42520–42526.

50. Swem DL, Swem LR, Setterdahl A, Bauer CE (2005) Involvement of *SenC* in assembly of cytochrome *c* oxidase in *Rhodobacter capsulatus*. *J Bacteriol* 187: 8081–8087.

51. Hannappel A, Bundschuh FA, Ludwig B (2012) Role of *Surfl* in heme recruitment for bacterial COX biogenesis. *Biochim Biophys Acta* 1817: 928–937.

52. Carr HS, George GN, Winge DR (2002) Yeast *Cox11*, a protein essential for cytochrome *c* oxidase assembly, is a *Cup*-binding protein. *J Biol Chem* 277: 31237–31242.

53. Hiser L, Di Valentin M, Hamer AG, Hosler JP (2000) *Cox11p* is required for stable formation of the *Cup* and magnesium centers of cytochrome *c* oxidase. *J Biol Chem* 275: 619–623.

54. Iwata S, Ostermeier C, Ludwig B, Michel H (1995) Structure at 2.8 Å resolution of cytochrome *c* oxidase from *Paracoccus denitrificans*. *Nature* 376: 660–669.

55. Ostermeier C, Harrenga A, Ermler U, Michel H (1997) Structure at 2.7 Å resolution of the *Paracoccus denitrificans* two-subunit cytochrome *c* oxidase complexed with an antibody FV fragment. *Proc Natl Acad Sci U S A* 94: 10547–10553.

56. Poole RK, Kumar C., Salmon L, Chance B. (1983) The 650 and chromophore in *Escherichia coli* is an “oxy-” or oxygenated compound, not the oxidized form of cytochrome oxidase *d*: an hypothesis. *Journal of General Microbiology* 129: 1335–1344.

57. Georgiou CD, Dueweke TJ, Gennis RB (1988) Regulation of expression of the cytochrome *d* terminal oxidase in *Escherichia coli* is transcriptional. *J Bacteriol* 170: 961–966.

58. Yang K, Borisov VB, Konstantinov AA, Gennis RB (2008) The fully oxidized form of the cytochrome *bd* quinol oxidase from *E. coli* does not participate in the catalytic cycle: direct evidence from rapid kinetics studies. *FEBS Lett* 582: 3705–3709.

59. Sone N, Tsukita S, Sakamoto J (1999) Direct correlation between proton translocation and growth yield: an analysis of the respiratory chain of *Bacillus stearothermophilus*. *J Biosci Bioeng* 87: 495–499.

60. Sakamoto J, Matsumoto A, Oobuchi K, Sone N (1996) Cytochrome *bd*-type quinol oxidase in a mutant of *Bacillus stearothermophilus* deficient in *caa3*-type cytochrome *c* oxidase. *FEMS Microbiol Lett* 143: 151–158.

61. Ng TC, Laheri AN, Maier RJ (1995) Cloning, sequencing, and mutagenesis of the cytochrome *c4* gene from *Azotobacter vinelandii*: characterization of the mutant strain and a proposed new branch in the respiratory chain. *Biochim Biophys Acta* 1230: 119–129.

62. Gilmour R, Krulwich TA (1997) Construction and characterization of a mutant of alkaliphilic *Bacillus firmus* OF4 with a disrupted *eta* operon and purification of a novel cytochrome *bd*. *J Bacteriol* 179: 863–870.

63. Chang HY, Ahn Y, Pace LA, Lin MT, Lin YH, et al. (2010) The diheme cytochrome *c4* from *Vibrio cholerae* is a natural electron donor to the respiratory *cbb3* oxygen reductase. *Biochemistry* 49: 7494–7503.

64. Tseng CP, Albrecht J, Gunsalus RP (1996) Effect of microaerophilic cell growth conditions on expression of the aerobic (*cyoABCDE* and *cydAB*) and anaerobic (*narGHJ1*, *frdABCD*, and *dmsABC*) respiratory pathway genes in *Escherichia coli*. *J Bacteriol* 178: 1094–1098.

65. Winstedt L, Yoshida K, Fujita Y, von Wachenfeldt G (1998) Cytochrome *bd* biosynthesis in *Bacillus subtilis*: characterization of the *cydABCD* operon. *J Bacteriol* 180: 6571–6580.

66. Preisig O, Zufferey R, Thony-Meyer L, Appleby CA, Hennecke H (1996) A high-affinity *cbb3*-type cytochrome oxidase terminates the symbiosis-specific respiratory chain of *Bradyrhizobium japonicum*. *J Bacteriol* 178: 1532–1538.

67. Mouncey NJ, Kaplan S (1998) Oxygen regulation of the *ccaN* gene encoding a component of the *cbb3* oxidase in *Rhodobacter sphaeroides* 2.4.1T: involvement of the *FnrL* protein. *J Bacteriol* 180: 2228–2231.

68. Baker SC, Ferguson SJ, Ludwig B, Page MD, Richter OM, et al. (1998) Molecular genetics of the genus *Paracoccus*: metabolically versatile bacteria with bioenergetic flexibility. *Microbiol Mol Biol Rev* 62: 1046–1078.

69. Pitcher RS, Watmough NJ (2004) The bacterial cytochrome *cbb3* oxidases. *Biochim Biophys Acta* 1655: 388–399.

70. Kawakami T, Kuroki M, Ishii M, Igarashi Y, Arai H (2010) Differential expression of multiple terminal oxidases for aerobic respiration in *Pseudomonas aeruginosa*. *Environ Microbiol* 12: 1399–1412.

71. Arai H (2011) Regulation and Function of Versatile Aerobic and Anaerobic Respiratory Metabolism in *Pseudomonas aeruginosa*. *Front Microbiol* 2: 103. 10.3389/fmicb.2011.00103.

72. Lauraeus M, Haltia T, Saraste M, Wikstrom M (1991) *Bacillus subtilis* expresses two kinds of haem-A-containing terminal oxidases. *Eur J Biochem* 197: 699–705.

73. Mather MW, Springer P, Fee JA (1991) Cytochrome oxidase genes from *Thermus thermophilus*. Nucleotide sequence and analysis of the deduced primary structure of subunit IIc of cytochrome *caas*. *J Biol Chem* 266: 5025–5035.

74. Lobo SA, Almeida CC, Carita JN, Teixeira M, Saraiva LM (2008) The haem-copper oxygen reductase of *Desulfovibrio vulgaris* contains a dihaem cytochrome *c* in subunit II. *Biochim Biophys Acta* 1777: 1528–1534.

75. Lamrabet O, Pieulle L, Aubert C, Mouhamar F, Stocker P, et al. (2011) Oxygen reduction in the strict anaerobe *Desulfovibrio vulgaris* Hildenborough: characterization of two membrane-bound oxygen reductases. *Microbiology* 157: 2720–2732.

76. Figurski DH, Helinski DR (1979) Replication of an origin-containing derivative of plasmid RK2 dependent on a plasmid function provided in trans. *Proc Natl Acad Sci U S A* 76: 1648–1652.

77. Herrero M, de Lorenzo V, Timmis KN (1990) Transposon vectors containing non-antibiotic resistance selection markers for cloning and stable chromosomal insertion of foreign genes in gram-negative bacteria. *J Bacteriol* 172: 6557–6567.