Matrix-assisted laser desorption ionization mass spectrometry of membrane proteins: demonstration of a simple method to determine subunit molecular weights of hydrophobic subunits

A Bacterial Stress Response Regulates Respiratory Protein Complexes To Control Envelope Stress Adaptation

The Cpx envelope stress response mediates adaptation to stresses that affect protein folding within the envelope of Gram-negative bacteria. Recent transcriptome analyses revealed that the Cpx response impacts genes that affect multiple cellular functions predominantly associated with the cytoplasmic membrane. In this study, we examined the connection between the Cpx response and the respiratory complexes NADH dehydrogenase I and cytochrome bo₃ in enteropathogenic Escherichia coli We found that the Cpx response directly represses the transcription of the nuo and cyo operons and that Cpx-mediated repression of these complexes confers adaptation to stresses that compromise envelope integrity. Furthermore, we found that the activity of the aerobic electron transport chain is reduced in E. coli lacking a functional Cpx response despite no change in the transcription of either the nuo or the cyo operon. Finally, we show that expression of NADH dehydrogenase I and cytochrome bo₃ contributes to basal Cpx pathway activity and that overproduction of individual subunits can influence pathway activation. Our results demonstrate that the Cpx response gauges and adjusts the expression, and possibly the function, of inner membrane protein complexes to enable adaptation to envelope stress.IMPORTANCE Bacterial stress responses allow microbes to survive environmental transitions and conditions, such as those encountered during infection and colonization, that would otherwise kill them. Enteric microbes that inhabit or infect the gut are exposed to a plethora of stresses, including changes in pH, nutrient composition, and the presence of other bacteria and toxic compounds. Bacteria detect and adapt to many of these conditions by using envelope stress responses that measure the presence of stressors in the outermost compartment of the bacterium by monitoring its physiology. The Cpx envelope stress response plays a role in antibiotic resistance and host colonization, and we have shown that it regulates many functions at the bacterial inner membrane. In this report, we describe a novel role for the Cpx response in sensing and controlling the expression of large, multiprotein respiratory complexes at the cytoplasmic membrane of Escherichia coli The significance of our research is that it will increase our understanding of how these stress responses are involved in antibiotic resistance and the mechanisms used by bacteria to colonize the gut.

Location of the Substrate Binding Site of the Cytochrome bo3 Ubiquinol Oxidase from Escherichia coli

Cytochrome bo₃ is a respiratory proton-pumping oxygen reductase that is a member of the heme-copper superfamily that utilizes ubiquinol-8 (Q₈H₂) as a substrate. The current consensus model has Q₈H₂ oxidized at a low affinity site (Q_L), passing electrons to a tightly bound quinone cofactor at a high affinity site (Q_H site) that stabilizes the one-electron reduced ubisemiquinone, facilitating the transfer of electrons to the redox active metal centers where O₂ is reduced to water. The current work shows that the Q₈ bound to the Q_H site is more dynamic than previously thought. In addition, mutations of residues at the Q_H site that do not abolish activity have been re-examined and shown to have properties expected of mutations at the substrate binding site (Q_L): an increase in the K_M of the substrate ubiquinol-1 (up to 4-fold) and an increase in the apparent K_i of the inhibitor HQNO (up to 8-fold). The data suggest that there is only one binding site for ubiquinol in cyt bo₃ and that site corresponds to the Q_H site.

Purification and biochemical properties of a cytochrome bc complex from the aerobic hyperthermophilic archaeon Aeropyrum pernix

Background

The bioenergetics of Archaea with respect to the evolution of electron transfer systems is very interesting. In contrast to terminal oxidases, a canonical bc₁ complex has not yet been isolated from Archaea. In particular, c-type cytochromes have been reported only for a limited number of species.

Results

Here, we isolated a c-type cytochrome-containing enzyme complex from the membranes of the hyperthermophilic archaeon, Aeropyrum pernix, grown aerobically. The redox spectrum of the isolated c-type cytochrome showed a characteristic α-band peak at 553 nm corresponding to heme C. The pyridine hemochrome spectrum also revealed the presence of heme B. In non-denaturing polyacrylamide gel electrophoresis, the cytochrome migrated as a single band with an apparent molecular mass of 80 kDa, and successive SDS-PAGE separated the 80-kDa band into 3 polypeptides with apparent molecular masses of 40, 30, and 25 kDa. The results of mass spectrometry indicated that the 25-kDa band corresponded to the hypothetical cytochrome c subunit encoded by the ORF APE_1719.1. In addition, the c-type cytochrome-containing polypeptide complex exhibited menaquinone: yeast cytochrome c oxidoreductase activities.

Conclusion

In conclusion, we showed that A. pernix, a hyperthemophilic archaeon, has a "full" bc complex that includes a c-type cytochrome, and to the best of our knowledge, A. pernix is the first archaea from which such a bc complex has been identified. However, an electron donor candidates for cytochrome c oxidase, such as a blue copper protein, have not yet been identified in the whole genome data of this archaeon. We are currently trying to identify an authentic substrate between a bc complex and terminal oxidase.

Mutation analysis of the interaction of B-type cytochrome c oxidase with its natural substrate cytochrome c-551

Heme-copper oxidases in the respiratory chain are classified into three subfamilies: A-, B- and C-types. Cytochrome bo(3)-type cytochrome c oxidase from thermophilic Bacillus is a B-type oxidase that is thought to interact with cytochrome c through hydrophobic interactions. This is in contrast to A-type oxidases, which bind cytochrome c molecules primarily through electrostatic forces between acidic residues in the oxidase subunit II and basic residues within cytochromes. In order to investigate the substrate-binding site in cytochrome bo(3), eight acidic residues in subunit II were mutated to corresponding neutral residues and enzymatic activity was measured using cytochrome c-551 from closely related Bacillus PS3. The mutation of E116, located at the interface to subunit I, decreased the k(cat) value most prominently without affecting the K(m) value, indicating that the residue is important for electron transfer. The mutation of D99, located close to the Cu(A) site, largely affected both values, suggesting that it is important for both electron transfer and substrate binding. The mutation of D49 and E84 did not affect enzyme kinetic parameters, but the mutation of E64, E66 and E68 lowered the affinity of cytochrome bo(3) for cytochrome c-551 without affecting the k(cat) value. These three residues are located at the front of the hydrophilic globular domain and distant from the Cu(A) site, suggesting that these amino acids compose an acidic patch for a second substrate-binding site. This is the first report on site-directed mutagenesis experiments of a B-type heme-copper oxidase.

Proteomics of mitochondrial inner and outer membranes

For the proteomic study of mitochondrial membranes, documented high quality mitochondrial preparations are a necessity to ensure proper localization. Despite the state-of-the-art technologies currently in use, there is no single technique that can be used for all studies of mitochondrial membrane proteins. Herein, we use examples to highlight solubilization techniques, different chromatographic methods, and developments in gel electrophoresis for proteomic analysis of mitochondrial membrane proteins. Blue-native gel electrophoresis has been successful not only for dissection of the inner membrane oxidative phosphorylation system, but also for the components of the outer membrane such as those involved in protein import. Identification of PTMs such as phosphorylation, acetylation, and nitration of mitochondrial membrane proteins has been greatly improved by the use of affinity techniques. However, understanding of the biological effect of these modifications is an area for further exploration. The rapid development of proteomic methods for both identification and quantitation, especially for modifications, will greatly impact the understanding of the mitochondrial membrane proteome.

A novel approach to analyze membrane proteins by laser mass spectrometry: from protein subunits to the integral complex

A novel laser-based mass spectrometry method termed LILBID (laser-induced liquid bead ion desorption) is applied to analyze large integral membrane protein complexes and their subunits. In this method the ions are IR-laser desorbed from aqueous microdroplets containing the hydrophobic protein complexes solubilized by detergent. The method is highly sensitive, very efficient in sample handling, relatively tolerant to various buffers, and detects the ions in narrow, mainly low-charge state distributions. The crucial experimental parameter determining whether the integral complex or its subunits are observed is the laser intensity: At very low intensity level corresponding to an ultrasoft desorption, the intact complexes, together with few detergent molecules, are transferred into vacuum. Under these conditions the oligomerization state of the complex (i.e., its quaternary structure) may be analyzed. At higher laser intensity, complexes are thermolyzed into subunits, with any residual detergent being stripped off to yield the true mass of the polypeptides. The model complexes studied are derived from the respiratory chain of the soil bacterium Paracoccus denitrificans and include complexes III (cytochrome bc(1) complex) and IV (cytochrome c oxidase). These are well characterized multi-subunit membrane proteins, with the individual hydrophobic subunits being composed of up to 12 transmembrane helices.

A targeted proteomic approach for the analysis of rat liver mitochondrial outer membrane proteins with extensive sequence coverage

Membrane proteins play an important role in cellular function. However, their analysis by mass spectrometry often is hindered by their hydrophobicity and/or low abundance. In this article, we present a method for the mass spectrometric analysis of membrane proteins based on the isolation of the resident membranes, isolation of the proteins by gel electrophoresis, and electroelution followed by enzymatic digestion by both trypsin and proteinase K. With this method, we have achieved 82-99% sequence coverage for the membrane proteins carnitine palmitoyltransferase-I (CPT-I), long-chain acyl-CoA synthetase (LCAS), and voltage-dependent anion channel (VDAC), isolated from rat liver mitochondrial outer membranes, including the transmembrane domains of these integral membrane proteins. This high sequence coverage allowed the identification of the isoforms of the proteins under study. This methodology provides a targeted approach for examining membrane proteins in detail.

Nonacid cleavable detergents applied to MALDI mass spectrometry profiling of whole cells

Although cleavable detergents were first synthesized a number of years ago, they have only recently been successfully applied to problems involving biological molecules. Recent reports have demonstrated that these compounds are useful for applications involving both 2D PAGE and mass spectrometry. However, most cleavable surfactants have utilized acid-labile functional groups to affect cleavage. In applications where extreme pH is required, acid cleavable detergents have limited usefulness. We report the synthesis of fluoride cleavable silane compounds and photolabile cinnamate esters as cleavable detergents having alternative cleavage chemistries than previously reported cleavable detergents. These compounds were applied to whole cell analysis using MALDI mass spectrometry, and it was demonstrated that their use results in an increase in the number of proteins analyzed by increasing protein solubility.

The Role of Extra Fragment at the C-terminal of Cytochrome b (Residues 421–445) in the Cytochrome bc1 Complex from Rhodobacter sphaeroides

Sequence alignment of cytochrome b of the cytochrome bc1 complex from various sources reveals that bacterial cytochrome b contain an extra fragment at the C terminus. To study the role of this fragment in bacterial cytochrome bc1 complex, Rhodobacter sphaeroides mutants expressing His-tagged cytochrome bc1 complexes with progressive deletion from this fragment (residues 421-445) were generated and characterized. The cytbDelta-(433-445) bc1 complex, in which 13 residues from the C-terminal end of this fragment are deleted, has electron transfer activity, subunit composition, and physical properties similar to those of the complement complex, indicating that this region of the extra fragment is not essential. In contrast, the electron transfer activity, binding of cytochrome b, ISP, and subunit IV to cytochrome c1, redox potentials of cytochromes b and c1 in the cytbDelta-(427-445), cytbDelta-(425-445), and cytbDelta-(421-445) mutant complexes, in which 19, 21, or all residues of this fragment are deleted, decrease progressively. EPR spectra of the [2Fe-2S] cluster and the cytochromes b in these three deletion mutant bc1 complexes are also altered; the extent of spectral alteration increases as this extra fragment is shortened. These results indicate that the first 12 residues (residues 421-432) from the N-terminal end of the C-terminal extra fragment of cytochrome b are essential for maintaining structural integrity of the bc1 complex.

Phosphorylation of rat liver mitochondrial carnitine palmitoyltransferase-I: effect on the kinetic properties of the enzyme

Hepatic carnitine palmitoyltransferase-I (CPT-IL) isolated from mitochondrial outer membranes obtained in the presence of protein phosphatase inhibitors is readily recognized by phosphoamino acid antibodies. Mass spectrometric analysis of CPT-IL tryptic digests revealed the presence of three phosphopeptides including one with a protein kinase CKII (CKII) consensus site. Incubation of dephosphorylated outer membranes with protein kinases and [gamma-32P]ATP resulted in radiolabeling of CPT-I only by CKII. Using mass spectrometry, only one region of phosphorylation was detected in CPT-I isolated from CKII-treated mitochondria. The sequence of the peptide and position of phosphorylated amino acids have been determined unequivocally as FpSSPETDpSHRFGK (residues 740-752). Furthermore, incubation of dephosphorylated outer membranes with CKII and unlabeled ATP led to increased catalytic activity and rendered malonyl-CoA inhibition of CPT-I from competitive to uncompetitive. These observations identify a new mechanism for regulation of hepatic CPT-I by phosphorylation.

Expression, purification, and characterization of BsSco, an accessory protein involved in the assembly of cytochrome c oxidase in Bacillus subtilis

The studies described here were performed to characterize further the plasma membrane associated protein BsSco, which is the product of the gene ypmQ, in Bacillus subtilis. BsSco is a member of the Sco family of proteins found in the inner mitochondrial membrane of yeast and humans and implicated as an accessory protein in the assembly of the Cu(A) site of cytochrome c oxidase. We have cloned the gene expressing BsSco, placed a six-histidine tag on its C-terminus, and over-expressed this protein in B. subtilis. Recombinant BsSco with the his-tag has been purified from Triton X-100-solubilized plasma membranes by nickel metal affinity chromatography. Mass spectral analysis of the purified protein is consistent with processing of BsSco by signal peptidase II removing an N-terminal putative transmembrane sequence to leave an acyl-glyceryl moiety at cysteine residue 19. Antibodies, raised against purified, recombinant BsSco, were used to characterize the timing of the level of native BsSco in batch cultures of wild-type B. subtilis. There is a marked lag in the level of native BsSco, but it does appear prior to cytochrome c oxidase, which is expressed in late stage growth. This work supports a role for BsSco in the assembly of the Cu(A) site of cytochrome c oxidase and its functional relationship to the Sco proteins found in eukaryotic cells.

Mass spectrometric detection of protein, lipid and heme components of cytochrome c oxidase from R. sphaeroides and the stabilization of non-covalent complexes from the enzyme

The cytochrome c oxidase enzyme from the Rhodobacter sphaeroides bacteria exists as a complex of four peptide subunits, two hemes, and a variety of lipids and metal ions held together by non-covalent forces. While the native enzyme functions as an associated unit, this complex usually dissociates during MALDI- TOF analysis. Through the use of matrix additives such as sucrose, the complete complex and partial complexes can be stabilized in the MALDI-TOF experiment. The dissociation of the complex allows for the detection of the components of the enzyme. The direct detection of associated lipids from an aqueous solution of the intact enzyme may eliminate the need for enzyme disruption and lipid extraction. The partial dissociation of multisubunit enzymes in such experiments may allow for the determination of subunit-subunit and subunit-lipid interactions

Mass spectrometric characterization of mitochondrial electron transport complexes: subunits of the rat heart ubiquinol-cytochrome c reductase

Complex III of the mitochondrial electron transport chain, ubiquinol-cytochrome c reductase, was isolated by blue native polyacrylamide gel electrophoresis. Ten of the 11 polypeptides present in this complex were detected directly by matrix-assisted laser desorption/ionization time-of-flight mass spectrometry (MALDI-TOFMS) following electroelution of the active complex. Tryptic and chymotryptic digestion of the complex permit the identification of specific peptides from all of the protein subunits with 70% coverage of the 250 kDa complex. The mass of all 11 proteins was confirmed by second dimension Tricine sodium dodecyl sulfate polyacrylamide gel electrophoresis (SDS-PAGE) and elution of the separated polypeptides. Additionally, the identity of the core I, core II, cytochrome c and the Rieske iron-sulfur protein were confirmed by liquid chromatography/tandem mass spectrometry (LC/MS/MS) characterization of the peptides generated by in-gel trypsin digestion of the SDS-PAGE separated proteins. The methodology demonstrated for analyzing this membrane-bound electron transport complex should be applicable to other membrane complexes, particularly the other mitochondrial electron transport complexes. The MS analysis of the peptides obtained by in-gel digestion of the intact complex permits the simultaneous characterization of the native proteins and modifications that contribute to mitochondrial deficits that have been implicated as contributing to pathological conditions.

Full subunit coverage liquid chromatography electrospray ionization mass spectrometry (LCMS+) of an oligomeric membrane protein: cytochrome b(6)f complex from spinach and the cyanobacterium Mastigocladus laminosus

Highly active cytochrome b(6)f complexes from spinach and the cyanobacterium Mastigocladus laminosus have been analyzed by liquid chromatography with electrospray ionization mass spectrometry (LCMS+). Both size-exclusion and reverse-phase separations were used to separate protein subunits allowing measurement of their molecular masses to an accuracy exceeding 0.01% (+/-3 Da at 30,000 Da). The products of petA, petB, petC, petD, petG, petL, petM, and petN were detected in complexes from both spinach and M. laminosus, while the spinach complex also contained ferredoxin-NADP(+) oxidoreductase (Zhang, H., Whitelegge, J. P., and Cramer, W. A. (2001) Flavonucleotide:ferredoxin reductase is a subunit of the plant cytochrome b(6)f complex. J. Biol. Chem. 276, 38159-38165). While the measured masses of PetC and PetD (18935.8 and 17311.8 Da, respectively) from spinach are consistent with the published primary structure, the measured masses of cytochrome f (31934.7 Da, PetA) and cytochrome b (24886.9 Da, PetB) modestly deviate from values calculated based upon genomic sequence and known post-translational modifications. The low molecular weight protein subunits have been sequenced using tandem mass spectrometry (MSMS) without prior cleavage. Sequences derived from the MSMS spectra of these intact membrane proteins in the range of 3.2-4.2 kDa were compared with translations of genomic DNA sequence where available. Products of the spinach chloroplast genome, PetG, PetL, and PetN, all retained their initiating formylmethionine, while the nuclear encoded PetM was cleaved after import from the cytoplasm. While the sequences of PetG and PetN revealed no discrepancy with translations of the spinach chloroplast genome, Phe was detected at position 2 of PetL. The spinach chloroplast genome reports a codon for Ser at position 2 implying the presence of a DNA sequencing error or a previously undiscovered RNA editing event. Clearly, complete annotation of genomic data requires detailed expression measurements of primary structure by mass spectrometry. Full subunit coverage of an oligomeric intrinsic membrane protein complex by LCMS+ presents a new facet to intact mass proteomics.

Cytochrome c oxidase contains an extra charged amino acid cluster in a new type of respiratory chain in the amino-acid-producing Gram-positive bacterium Corynebacterium glutamicum

The membranes from Corynebacterium glutamicum cells contain a hydrophobic di-haem C protein as the cytochrome c subunit of the new type of cytochrome bc complex (complex III in the respiratory chain) encoded by the qcrCAB operon [Sone, N., Nagata, K., Kojima, H., Tajima, J., Kodera, Y., Kanamaru, T., Noguchi, S. & Sakamoto, J. (2001). Biochim Biophys Acta 1503, 279-290]. To characterize complex IV, cytochrome c oxidase and its structural genes were isolated. The oxidase is of the cytochrome aa(3) type, but mass spectrometry indicated that the haem is haem As, which contains a geranylgeranyl side-chain instead of a farnesyl group. The enzyme is a SoxM-type haem-copper oxidase composed of three subunits. Edman degradation and mass spectrometry suggested that the N-terminal signal sequence of subunit II is cleaved and that the new N-terminal cysteine residue is diacylglycerated, while neither subunit I nor subunit III is significantly modified. The genes for subunits II (ctaC) and III (ctaE) are located upstream of the qcrCAB operon, while that for subunit I (ctaD) is located separately. The oxidase showed low enzyme activity with extrinsic substrates such as cytochromes c from horse heart or yeast, and has the Cu(A)-binding motif in its subunit II. A prominent structural feature is the insertion of an extra charged amino acid cluster between the beta2 and beta4 strands in the substrate-binding domain of subunit II. The beta2-beta4 loop of this oxidase is about 30 residues longer than that of major cytochrome c oxidases from mitochondria and proteobacteria, and is rich in both acidic and basic residues. These findings suggest that the extra charged cluster may play a role in the interaction of the oxidase with the cytochrome c subunit of the new type of bc complex.

Structural characterization of two aquaporins isolated from native spinach leaf plasma membranes

Two members of the aquaporin family, PM28A and a new one, PM28C, were isolated and shown to be the major constituents of spinach leaf plasma membranes. These two isoforms were identified and characterized by matrix-assisted laser desorption ionization-mass spectrometry. Edman degradation yielded the amino acid sequence of two domains belonging to the new isoform. PM28B, a previously described isoform, was not found in our preparations. Scanning transmission electron microscopy mass analysis revealed both PM28 isoforms to be tetrameric. Two types of particles, a larger and a smaller one, were found by transmission electron microscopy of negatively stained solubilized proteins and by atomic force microscopy of PM28 two-dimensional crystals. The ratio of larger to smaller particles observed by transmission electron microscopy and single particle analysis correlated with the ratio of PM28A to PM28C determined by matrix-assisted laser desorption ionization-mass spectrometry. The absence of PM28B and the ratio of PM28A to PM28C indicate that these plasma membrane intrinsic proteins are differentially expressed in spinach leaves. These findings suggest that differential expression of the various aquaporin isoforms may regulate the water flux across the plasma membrane, in addition to the known mechanism of regulation by phosphorylation.

Interaction of membrane proteins and lipids with solubilizing detergents

Detergents are indispensable in the isolation of integral membrane proteins from biological membranes to study their intrinsic structural and functional properties. Solubilization involves a number of intermediary states that can be studied by a variety of physicochemical and kinetic methods; it usually starts by destabilization of the lipid component of the membranes, a process that is accompanied by a transition of detergent binding by the membrane from a noncooperative to a cooperative interaction already below the critical micellar concentration (CMC). This leads to the formation of membrane fragments of proteins and lipids with detergent-shielded edges. In the final stage of solubilization membrane proteins are present as protomers, with the membrane inserted sectors covered by detergent. We consider in detail the nature of this interaction and conclude that in general binding as a monolayer ring, rather than as a micelle, is the most probable mechanism. This mode of interaction is supported by neutron diffraction investigations on the disposition of detergent in 3-D crystals of membrane proteins. Finally, we briefly discuss the use of techniques such as analytical ultracentrifugation, size exclusion chromatography, and mass spectrometry relevant for the structural investigation of detergent solubilized membrane proteins.

Detergent-solubilized Escherichia coli cytochrome bo3 ubiquinol oxidase: a monomeric, not a dimeric complex

The protein molecular weight, M(r), and hydrodynamic radius, R(s), of Triton X-100-solubilized Escherichia coli cytochrome bo3 were evaluated by computer fitting of sedimentation velocity data with finite element solutions to the Lamm equation. Detergent-solubilized cytochrome bo3 sediments as a homogeneous species with an S20,w of 6.75 s and a D20,w of 3.71 x 10(-7) cm2/s, corresponding to a R(s) of 5.8 nm and a M(r) of 144,000 +/- 3500. The protein molecular weight agrees very well with the value of 143,929 calculated from the four known subunit sequences and the value of 143,025 measured by MALDI mass spectrometry for the histidine-tagged enzyme. We conclude that detergent-solubilized E. coli ubiquinol oxidase is a monomeric complex of the four known subunits.

Sequence analysis of cytochrome bd oxidase suggests a revised topology for subunit I

Numerous sequences of the cytochrome bd quinol oxidase (cytochrome bd) have recently become available for analysis. The analysis has revealed a small number of conserved residues, a new topology for subunit I and a phylogenetic tree involving extensive horizontal gene transfer. There are 20 conserved residues in subunit I and two in subunit II. Algorithms utilizing multiple sequence alignments predicted a revised topology for cytochrome bd, adding two transmembrane helices to subunit I to the seven that were previously indicated by the analysis of the sequence of the oxidase from E. coli. This revised topology has the effect of relocating the N-terminus and C-terminus to the periplasmic and cytoplasmic sides of the membrane, respectively. The new topology repositions I-H19, the putative ligand for heme b595, close to the periplasmic edge of the membrane, which suggests that the heme b595/heme d active site of the oxidase is located near the outer (periplasmic) surface of the membrane. The most highly conserved region of the sequence of subunit I contains the sequence GRQPW and is located in a predicted periplasmic loop connecting the eighth and ninth transmembrane helices. The potential importance of this region of the protein was previously unsuspected, and it may participate in the binding of either quinol or heme d. There are two very highly conserved glutamates in subunit I, E99 and E107, within the third transmembrane helix (E. coli cytochrome bd-I numbering). It is speculated that these glutamates may be part of a proton channel leading from the cytoplasmic side of the membrane to the heme d oxygen-reactive site, now placed near the periplasmic surface. The revised topology and newly revealed conserved residues provide a clear basis for further experimental tests of these hypotheses. Phylogenetic analysis of the new sequences of cytochrome bd reveals considerable deviation from the 16sRNA tree, suggesting that a large amount of horizontal gene transfer has occurred in the evolution of cytochrome bd.