If you don't remember your password, you can reset it by entering your email address and clicking the Reset Password button. You will then receive an email that contains a secure link for resetting your password
If the address matches a valid account an email will be sent to __email__ with instructions for resetting your password
To whom correspondence should be addressed: School of Biosciences, Centre for Biocatalysis, University of Exeter, Stocker Road, Exeter EX4 4QD, United Kingdom. Tel.: 44-1329-264-675; Fax: 44-1392-263-434.
* This work was funded in part by Biotechnology and Biological Sciences Research Council (BBSRC) Grants P17219 and BBS/B/10110 (to C. S. B. and D. J. R.) and BBSRC/Engineering and Physical Sciences Research Council (life science interface) Grant BB/D00781X/1 (to C. S. B., R. J. L., and D. J. R.). Sequencing of the T. selenatis genome was funded by the Systems Biology Theme, University of Exeter. The on-line version of this article (available at http://www.jbc.org) contains supplemental Figs. S1 and S2. 1 A recipient of a BBSRC Agri-food committee studentship. 2 Supported by a La Trobe University postgraduate award.
Selenate reductase (SER) from Thauera selenatis is a periplasmic enzyme that has been classified as a type II molybdoenzyme. The enzyme comprises three subunits SerABC, where SerC is an unusual b-heme cytochrome. In the present work the spectropotentiometric characterization of the SerC component and the identification of redox partners to SER are reported. The mid-point redox potential of the b-heme was determined by optical titration (Em + 234 ± 10 mV). A profile of periplasmic c-type cytochromes expressed in T. selenatis under selenate respiring conditions was undertaken. Two c-type cytochromes were purified (∼24 and ∼6 kDa), and the 24-kDa protein (cytc-Ts4) was shown to donate electrons to SerABC in vitro. Protein sequence of cytc-Ts4 was obtained by N-terminal sequencing and liquid chromatography-tandem mass spectrometry analysis, and based upon sequence similarities, was assigned as a member of cytochrome c4 family. Redox potentiometry, combined with UV-visible spectroscopy, showed that cytc-Ts4 is a diheme cytochrome with a redox potential of +282 ± 10 mV, and both hemes are predicted to have His-Met ligation. To identify the membrane-bound electron donors to cytc-Ts4, growth of T. selenatis in the presence of respiratory inhibitors was monitored. The specific quinol-cytochrome c oxidoreductase (QCR) inhibitors myxothiazol and antimycin A partially inhibited selenate respiration, demonstrating that some electron flux is via the QCR. Electron transfer via a QCR and a diheme cytochrome c4 is a novel route for a member of the DMSO reductase family of molybdoenzymes.
)) and catalyze either hydride or oxygen transfer as generalized by Reaction 1.
(REACTION 1)
These soluble enzymes consist of three subunits and in addition to the b-heme cytochrome (γ-subunit), they comprise an iron-sulfur protein (β-subunit) coordinating 1 × [3Fe-4S] cluster and 3 × [4Fe-4S] clusters, and a catalytic component (α-subunit) that coordinates a [4Fe-4S] cluster and the active site molybdopterin guanine dinucleotide cofactor (
) (Fig. 1). The reductases play a pivotal function, coupling the reduction of substrates to the generation of the proton-motive force (PMF). Identifying the route by which electrons are transferred to these reductases is vital to understanding their bioenergetics (
). How periplasmic substrate reduction can generate a PMF, which is sufficient to support growth, is of considerable interest. The use of selenate and selenite as bacterial respiratory substrates has been well documented (
). By far the most well studied selenate-respiring bacterium has been T. selenatis; however, the characterization of the selenate respiratory system has not yet addressed the mechanism by which electrons are transferred from the cytoplasmic membrane to the periplasmic selenate reductase (SerABC). Because SerABC appears to be soluble and not associated with the inner or outer membrane (
), it is likely that another protein acts as a shuttle to take electrons from the membrane and deliver them to the b-heme of SerC. Reduced SerC then donates electrons, via the iron-sulfur clusters of SerAB (
FIGURE 1Subunit organization and redox cofactors of periplasmic type II molybdoenzymes. Direction of electron transfer is indicated by the arrows. Reductases include selenate and chlorate reductase. Dehydrogenases include dimethylsulfide and ethylbenzene dehydrogenases.
The b-heme of SerC is presumed to be the site of electron entry to the selenate reductase (SER) enzyme complex but also remains poorly studied. The recently solved crystal structure of EBDH reveals that its γ-subunit also binds an unusual b-type heme, which is coordinated by methionine and lysine ligands (
). This coordination favors the uncharged ferrous heme over the charged ferric form and as such modulates a relatively high redox potential of +254 mV. Similarly, optical potentiometric redox titrations of the heme b in the γ-subunit of dimethylsulfide dehydrogenase (DMSDH) revealed a high redox potential of Em = +324mV (
). By contrast, the cytoplasmic-facing respiratory nitrate reductase NAR from Escherichia coli (another member of the type II molybdoenzymes) possesses a membrane-bound γ-subunit (NarI) that coordinates two bis-His-ligated b-hemes that have much lower redox potentials at +20 and +120 mV (
). In the case of NAR, the low redox potentials of the heme favors electron flow away from NarI, toward the molybdenum guanine dinucleotide cofactor in NarG, and ultimately resulting in the reduction of nitrate to nitrite. Critically, the dehydrogenases (EBDH and DMSDH) function in the reverse direction (Fig. 1); it has been argued that, due to the unusually high potential b-heme in γ-subunits of both EBDH and DMSDH, that the γ-subunit is the site of electron egress from these enzymes to their redox partners (
). Consequently, the high potential b-heme might function to draw electrons away from, rather than toward, the molybdopterin guanine dinucleotide cofactor during catalysis as displayed in the nitrate reductase enzymes. Analysis of the amino acid sequence alignment of EBDH and DMSDH with selenate reductase (SER) shows that the coordinating methionine (Met-138) and lysine (Lys-228) residues are also conserved in SerC. This then raises interesting questions as to what is the redox potential for the b-heme in the reductase members of this distinct type II clade and what are the likely redox partners?
Among the related type II molybdoenzymes, there are a number of different routes of electron transfer to terminal reductases, although the majority catalyze reactions with a lower redox potential than the selenate/selenite couple (+475 mV). For example; the DMSO/DMS couple is +160 mV and trimethylamine N-oxide/trimethylamine couple is +130 mV (
). Within the DMSO reductase group, a number of pathways for electron transport have been established, and a variety of membrane-bound proteins have been implicated. These membrane-bound cytochromes are responsible for extracting electrons from the quinol pool so that they can be transferred to the associated downstream terminal reductases. The existing literature provides several routes by which electrons are transferred. For example, the membrane-bound nitrate reductase NarGH (
) receives its electrons from its dedicated NarI membrane integral subunit, which contains two b-hemes; the DMSO reductase from E. coli has a similar organization of subunits (
) accept electrons from dedicated tetra- or penta-heme quinol dehydrogenases in the membrane (TorC/DorC/NapC). In the case of Nap, there is also an alternative pathway; the nap gene cluster in E. coli for example, also contains napG and napH, which encode iron-sulfur proteins responsible for electron transfer from ubiquinol to NapC (
) postulated that the selenate reductase of T. selenatis could receive electrons from quinol-cytochrome c oxidoreductase (QCR) via a high potential c-type cytochrome, due to the high potential of the selenate/selenite couple. This has not been seen for any of the other type II molybdoenzymes, so showing electron donation from QCR, via a soluble cytochrome c to SerABC, would represent a novel electron transport chain within this group. The aim of the present work was to resolve the electron transfer pathways during selenate respiration in T. selenatis.
EXPERIMENTAL PROCEDURES
Growth of T. selenatis, Production of Periplasmic Fractions, and Purification of SER
T. selenatis was grown anaerobically at 30 °C in mineral salts medium (
) containing yeast extract (0.1%), with either selenate or nitrite (10 mm) as terminal electron acceptors and acetate (10 mm) as the electron donor in 10-liter batch cultures. Cultures were harvested during late log phase (after 16–18 h growth) at A600 0.6–0.7, and spheroplasts were prepared as described previously (
). The spheroplasts were removed by centrifugation (25,000 × g, 20 min), and the supernatant, containing the periplasm, was retained. Periplasms were analyzed by SDS-PAGE, and c-type cytochromes were visualized using the heme stain method. SER was purified as detailed by Schröder et al. (
). Purified enzyme was stored at −80 °C. Purity was determined by SDS-PAGE, and positive identification of the presence of the SerC component was achieved by N-terminal sequence analysis as described previously (
Soluble protein from a 10-liter culture was loaded onto a Q-Sepharose Fast-Flow anion-exchange column (70-ml bed volume, Amersham Biosciences) pre-equilibrated with 3 column volumes of 30 mm Tris, pH 8.5. The column was washed with a further column volume of buffer to remove unbound proteins. A 24-kDa c-type cytochrome was eluted using a 300-ml gradient of 0–300 mm NaCl in 30 mm Tris, pH 8.5. Protein concentration was monitored using absorbance at 280 nm, and the presence of cytochromes was monitored using absorbance at 410 nm. The fractions containing the 24-kDa c-type cytochrome as determined by SDS-PAGE were pooled and concentrated using a 15-ml 10-kDa molecular weight cut-off centrifugal concentrator at 4000 × g. The concentrated protein was loaded onto a Superdex 200 16/60 gel-filtration column (Amersham Biosciences) equilibrated with 2 column volumes of 30 mm Tris, pH 8.5, and eluted in the same buffer. The fractions containing the 24-kDa c-type cytochrome were pooled and loaded onto a 1-ml Mono Q anion-exchange column (Amersham Biosciences) equilibrated with 30 mm Tris, pH 8.5, and eluted with a gradient of 0–1 m NaCl in 30 mm Tris, pH 8.5. As a final purification step, the fractions containing the 24-kDa c-type cytochrome were pooled and transferred to dialysis tubing (6-kDa MWCO, Fisher), and the protein was then dialyzed overnight at 4 °C, with gentle stirring, into 30 mm Tris, pH 8.5, containing 1 m ammonium sulfate. The dialyzed protein was bound to a 25-ml Fast-Flow phenyl-Sepharose column (Amersham Biosciences) pre-equilibrated with 1 m ammonium sulfate in 30 mm Tris, pH 8.5, and subsequently eluted with a decreasing gradient of 1 to 0 m ammonium sulfate. The cytochrome was judged pure by SDS-PAGE, concentrated as before, and stored at 4 °C for short term use, or −80 °C for longer term storage. All chromatography procedures were carried out using an AKTAPrime (Amersham Biosciences) purification system, except for the Mono Q column, which was run on an AKTABasic column (Amersham Biosciences). The soluble fraction from a 2-liter culture was also prepared as described above. This soluble fraction was concentrated using a 15-ml 5-kDa molecular weight cut-off centrifugal concentrator (Millipore) to a final volume of 5 ml, and then loaded onto a Superdex 200 16/60 gel-filtration column pre-equilibrated with 20 mm Tris, pH 8.0. Proteins were eluted in the same buffer. Fractions were monitored for presence of protein using the absorbance at 280 nm for the presence of cytochrome using the absorbance at 410 nm. Fractions containing the 6-kDa cytochrome were identified by SDS-PAGE, pooled, and concentrated in a 5-ml 5-kDa MWCO concentrator.
Electronic Absorption Spectroscopy
UV-visible spectra of protein samples were recorded on a Varian Cary 4E UV-visible spectrophotometer between 350 and 650 nm. Reduced samples were achieved by adding excess sodium dithionite and oxidized samples by the addition of potassium ferricyanide.
Electron Transfer Assay
To determine whether purified cytochromes were able to donate electrons to SerABC in vitro, a solution of ∼5 μm cytochrome in 50 mm phosphate buffer (pH 7.5) was degassed in a sealed cuvette. A weak solution (10 mm) of sodium dithionite was also degassed and then titrated into the cuvette using a 10-μl Hamilton syringe, until the cytochrome was fully reduced, as determined by wavelength scanning UV-visible spectroscopy (350–700 nm). SerABC and selenate were added to final concentrations of 1 μm and 20 mm, respectively, and re-oxidation of the cytochrome was monitored by recording the spectra. Cytochromes were similarly tested for selenate and selenite reductase activity by incubation with the relevant substrate but no SerABC.
Optical Redox Titrations
Optical spectra of the samples were measured on a Hitachi U-3310 spectrophotometer between 350 and 700 nm. Protein samples were added to a stirred cuvette constantly sparged with oxygen-free nitrogen to minimize oxidation from air. Redox mediators (2,3,5,6-tetramethyl-p-phenylendiamine, 1,2-napthoquinone, phenozinemethosulfate, phenylethosulfate, juglone, duroquinone, and menadion) were added to a final concentration of 6 μm each, and the potential was measured using an electrode (
). The electrode was initially calibrated with reference to a saturated solution of quinhydrone (295 mV versus standard hydrogen electrode). A degassed solution of 5 mm sodium dithionite was prepared and added in 1-μl steps, recording the potential, and the spectra at each step. Once the sample was fully reduced, 3 mm potassium ferricyanide was titrated in similarly to re-oxidize the sample. The fraction of protein reduced at each potential was calculated and fitted to a Nernstian equation using the curve-fitting program TableCurve 2D (Systat Software Inc.).
Sequence Analysis of the 24-kDa Cytochrome
N-terminal sequencing of the 24-kDa cytochrome was undertaken at the Pinnacle Proteomic facility, Newcastle University. Native 24-kDa cytochrome c was also digested in solution using Trypsin Gold (Promega, Madison, WI) according to the manufacturer's instructions. Briefly, Trypsin Gold in 50 mm acetic acid was mixed with the cytochrome (in 30 mm Tris, pH 8.5) in a ratio of 1:20 w/w protease:protein. The mixture was incubated at 37 °C for 1 h and analyzed by SDS-PAGE. Protein bands of interest were excised for analysis, which was carried out using liquid chromatography/tandem mass spectrometry at the University of York Proteomics Department.
Sequencing the T. selenatis Genome Using Illumina Sequencing
Genomic DNA was extracted and fragmented to using the Bioruptor apparatus, and fragments of 200 bp were selected using an agarose gel-purification protocol. These fragments were then prepared using the standard Illumina paired-end library-preparation protocols. The resulting DNA was purified using an acrylamide gel, and these were sequenced using an Illumina GA2 sequencer using a single lane loaded at 6 pm. Illumina version 2 reagents were used together with Illumina SCS 2.3 and Illumina Pipeline 1.3 software to base call. This produced 15,422,089 36-bp paired-end sequences. These were subsequently trimmed to remove or shorten reads that contained adaptor sequence. Read 2 contained unusually high error rates toward the 3′-end of the reads. Where appropriate these were shortened prior to any de novo assembly to remove bases with quality scores of <20.
MAQ (available on-line from sourceforge.net) was used to align these sequences to the reference genome Thauera MZ1T (GenBankTM accession CP001281 and CP001282). This was achieved with an average depth of 32.9 across non-gap regions. Unmapped reads were then collated and assembled de novo using Velvet (version 0.7.47) in single-end mode using a hash length of 21 (
). These unmapped reads generated an assembly with 14,870 contigs and an N50 length of 1223 bp. These contigs represent genomic regions present in T. selenatis but which are absent in Thauera MZ1T. Additionally a completely de novo assembly was produced using Velvet in paired-end mode using a parameter sweep of hash length and coverage cut-off. A hash length of 25 and a coverage cut-off of 5 were found to produce the highest N50. This resulted in an assembly of 769 contigs with an N50 value of 135,919 bp (made up of 17 contigs) and a total length of 9.48 Mb.
Micro-culture Growth and Respiratory Inhibition Studies
To monitor the effects of changing culture conditions and inhibitors on the growth of T. selenatis, a method utilizing micro-culture was developed. A 500-ml bottle of mineral salts medium was prepared without the addition of selenate or other electron acceptor. The bottle was sealed with rubber septa and sparged with oxygen-free nitrogen for at least 1 h. Each well of a 96-well microtiter plate had 250 μl of the degassed medium added to it, and electron acceptor was added in varying concentrations from degassed stock solutions. Inhibitors, myxothiazol, antimycin A, and 2-n-heptyl-4-hydroxyquinoline N-oxide (HQNO), at a range of concentrations, were added to wells as required. Wells were inoculated with 1 μl of washed cells from an exponential phase T. selenatis culture (1 ml of cells at A600 nm 0.6–0.7 were harvested by centrifugation at 12,000 × g, 2 min, and resuspended in 1 ml of the degassed medium). Inoculated plates were incubated in a FLUOstar Optima microplate reader (BMG Labtech) attached to a nitrogen cylinder to keep the plate in an anaerobic atmosphere. The plate was maintained at a constant temperature of 30 °C, and A600 nm values were measured in each well every 15 min using the kinetics-measuring program, with 5-s shaking before each reading. Each plate had 6 wells without the addition of bacteria to ensure no contaminants were growing, and to act as a baseline, and 6 wells with non-inhibiting growing conditions for T. selenatis to check the culture was growing typically. The maximum gradient at exponential growth was determined by fitting the growth curves to the Gompertz equation (
). All growth conditions were replicated at a range of substrate concentrations and using non-linear regression; the maximum growth rate (μmax) of the culture was calculated using the Monod equation (Equation 1) (
The c-type cytochromes expressed during selenate respiration in T. selenatis were profiled by heme staining SDS-PAGE resolved periplasmic proteins. Several soluble cytochromes were identified (referred to as cytc-Ts1–Ts7, in size order) (Fig. 2A) and by comparing the expression profile to nitrate grown cells, a number of cytochromes were found to be more highly expressed during selenate respiration. cytc-Ts1 was present under both growth conditions but to a lesser extent when grown on selenate. cytc-Ts2 was poorly resolved but present as a band at ∼40 kDa in selenate-grown cells. cytc-Ts3 (∼33 kDa) was again expressed under both growth conditions. Cytochromes identified as cytc-Ts5 and cytc-Ts6 were resolved only in selenate grown cells. cytc-Ts4 (24 kDa) and cytc-Ts7 (6 kDa) were both significantly up-regulated under selenate growth conditions and selected for further investigation. cytc-Ts4 and cytc-Ts7 were purified (Fig. 2A) and tested as potential electron donors to SerABC.
FIGURE 2Cytochrome c profile from T. selenatis.A, SDS-PAGE gels stained for c-type cytochromes, showing periplasmic cytochromes (cytc-Ts1–7) expressed during selenate and/or nitrate respiration. Lane 1, Invitrogen SeeBlue Plus2 Prestained Standard; lane 2, periplasmic fraction from selenate grown cells (heme stained); and lane 3, periplasmic fraction from nitrate grown cells (heme stained). 7 μg of protein was loaded in each lane. Lane 4, purified cytc-Ts4 (∼24 kDa) stained with Invitrogen SeeBlue Plus2; lane 5, purified cytc-Ts4 (∼24 kDa) stained with heme stain; lane 6, purified cytc-Ts7 (∼6 kDa) stained with Invitrogen SeeBlue Plus2; and lane 7, purified cytc-Ts7 (∼6 kDa) stained with heme stain. B, spectra of 25 μm cytc-Ts4 (∼24-kDa protein) as purified (solid line), reduced with sodium dithionite (dotted line), and oxidized with potassium ferricyanide (dashed line). C, spectra of 25 μm cytc-Ts7 (∼6-kDa protein) as purified (solid line), reduced with sodium dithionite (dotted line), and oxidized with potassium ferricyanide (dashed line).
The reduced UV-visible spectrum of the purified 6-kDa cytochrome (cytc-Ts7) was typical of a low spin c-type cytochrome, with absorbance maxima at 415, 522, and 551 nm (Fig. 2C). The oxidized spectrum shows a shift in the Soret band to 409 nm and a broad feature around 525 nm. To test the ability of this cytochrome to donate electrons to SerABC, it was reduced by titration with a weak solution of sodium dithionite, and mixed with purified selenate reductase and selenate in an anaerobic assay. No re-oxidation of the heme was observed suggesting that the 6-kDa cytochrome is not able to act as an electron donor to SerABC in vitro (data not shown).
Characterization of the 24-kDa Cytochrome (cytc-Ts4)
UV-visible spectroscopy of the 24-kDa cytochrome (cytc-Ts4) was again typical of a low spin heme c (Fig. 2B). It was purified in a partially reduced state, but did not readily become oxidized, even after lengthy exposure to air. Protein reduced with dithionite showed a Soret band at 414 nm, an α-band at 550 nm, and a β-band at 519 nm. The α-band showed asymmetry, with a shoulder at 545 nm suggesting a composite band, possibly due to the presence of more than one heme in the protein. Oxidized with ferricyanide, the cytochrome had a Soret band at 409 nm and a broad feature at 528 nm. As for cytc-Ts7, cytc-Ts4 was tested for its ability to donate electrons to SerABC in an anaerobic assay. In this case, after incubation with enzyme and substrate, the visible absorbance spectrum was shifted (Fig. 3). The Soret band shifted from 414 to 410 nm, and decreased in intensity, and the α/β peaks broadened and decreased in intensity, indicative of re-oxidation of the heme. cytc-Ts4 incubated with SerABC, selenate, or selenite alone did not exhibit any shifts in absorbance (data not shown), suggesting that this cytochrome is capable of donating electrons to SerABC for the reduction of selenate to selenite in vitro, and therefore could be an in vivo electron donor to SerABC. The rate of cytc-Ts4 heme oxidation was determined to be some 50 times slower than the observed rate of SER turnover determined using methyl viologen as the electron donor (
) have shown that NAR from Paracoccus denitrificans has a 50-fold higher Vmax with methyl viologen than with duroquinol.
FIGURE 3Spectroscopic assay of electron donation from cytc-Ts4 to SerABC. cytc-Ts4 (5 μm) reduced with dithionite (dashed line) was mixed with 1 μm SerABC and 20 mm selenate and re-oxidation was monitored by wavelength scanning UV-visible spectroscopy (solid line).
The redox potential of the heme of cytc-Ts4 was measured spectrophotometrically by monitoring the absorbance of the α-band at 551 nm (adjusted by reference to an isosbestic point at 560 nm) during reduction or oxidation of the heme (Fig. 4A). The fraction of protein reduced at each measured potential was calculated, then plotted as a function of redox potential versus the standard hydrogen electrode (Fig. 4B). Nernst curves for n = 1 and n = 2 (where n is the number of electrons transferred) were fitted and are also shown in Fig. 4B. Interestingly, the data fits more closely to an n = 2 curve, suggesting the cytochrome has two hemes. A mid-point potential of +282 ± 10 mV for the cytochrome was calculated, but we were unable to resolve two separate potentials.
FIGURE 4Optical redox titration of cytc-Ts4.A, α and β region of cytc-Ts4 spectrum during titrations. The black arrow indicates increasing redox potential. B, fraction of cytochrome reduced as calculated with reference to absorbance at 551 nm, as a function of redox potential. Nernst curves for n = 1 electrons (dashed line) and n = 2 electrons (solid line) are also shown.
MALDI-TOF analysis of cytc-Ts4 gave a molecular mass of 23,558 Da, not including heme groups. N-terminal sequencing (Pinnacle Proteomic facility, Newcastle University) and liquid chromatography/tandem mass spectrometry of tryptic digest fragments (University of York Proteomics Department) resulted in 60 amino acids of sequence data, which was used to identify the gene for this cytochrome from the draft genome assembly of T. selenatis. Peptides were searched against a custom Blast (
) database containing all contigs generated from the sequence data, which could not be mapped to Thauera MZ1T and separately to another database containing the complete set of contigs generated from all Illumina T. selenatis data. Due to the short lengths of the peptides only two significant hits were obtained, and these contigs were used as targets for further investigation. Open reading frame prediction was carried out using the CLC-Bio software package, and candidate regions were searched via Blast against the non-redundant database of protein sequences. The sequence of cytc-Ts4 (GenBankTM accession GU570563) was used to search the database, and an alignment of similar cytochromes is shown in supplemental Fig. S1. Alignment of the predicted amino acid sequences shows the presence of two conserved CXXCH motifs consistent with the binding of two covalently attached hemes. The presence of two further conserved methionine residues suggests that both hemes are His-Met coordinated.
Determination of the Redox Potential of the b-heme in SerC
Purified SER was used to determine the midpoint potential of the b-heme in SerC, the point of electron entry to the reductase from the cytc redox partner. The redox potential of the heme was determined by optical redox titration by monitoring the absorbance of the heme at 558 nm (adjusted with reference to an isosbestic point at 575 nm) during sequential oxidation and reduction of the heme (Fig. 5A). The slight drift observed in the spectrum was due to the contribution of the [Fe-S] clusters in SerAB. The absorbance was monitored at the α-band wavelength rather than that of the Soret band, because the redox mediators used do not absorb in the α-band region, but might interfere at Soret band wavelengths. The fraction of protein reduced at each measured potential was calculated and then plotted as a function of redox potential (Fig. 5B). The Nernst curve corresponding to n = 1 electrons was fitted to the data giving a midpoint potential of +234 ± 10 mV.
FIGURE 5Optical redox titration of SerC.A, the α and β region of b-heme spectrum during titrations. The black arrow indicates increasing redox potential. Samples were prepared in 30 mm Tris-HCl, pH 7.5. B, fraction of cytochrome b reduced as calculated with reference to absorbance at 558 nm, as a function of redox potential. Nernst curve for n = 1 electron.
Quinol-cytochrome c Oxidoreductase Is an Electron Donor to Selenate Reductase in T. selenatis
Soluble cytochromes acting as electron carriers need to accept electrons from an upstream donor. Commonly, this comes in the form of a QCR, which can contribute to the PMF via the Q-cycle. To investigate whether the QCR is involved in selenate respiration and could therefore be responsible for electron transfer to the di-heme cytochrome (cytc-Ts4), classic inhibitors of respiratory chains were used during selenate and nitrite respiration. The inhibitors myxothiazol and antimycin A specifically inhibit the QCR, but at different sites (
). Myxothiazol binds to the Qo site on the b heme subunit of QCR. The Qo site is where quinol is oxidized to quinone, and the binding of myxothiazol stops electrons being transferred from the Q-pool to downstream electron acceptors. By contrast, antimycin A binds to the Qi site of QCR preventing the reduction of quinone to quinol. HQNO is a less specific quinone analogue, which inhibits not only the QCR (
To culture T. selenatis under more reproducible conditions, and allow multiple replicates of growth curves, a method of culturing T. selenatis in a 96-well microplate was devised. All growth conditions were replicated in 10 wells, allowing averaging of the optical densities over each time point. Before averaging, the optical density of a control set of wells of medium with no bacterial culture was subtracted from the growth curves. To test whether concentrations of inhibitor were high enough to inhibit T. selenatis, the microplate method was used to investigate the effect of myxothiazol on nitrite respiration, to see if complete inhibition could be achieved. Thauera sp. is known to express a cd1 nitrite reductase (
). A range of nitrite concentrations was used, with myxothiazol at 10 μm. Ten wells contained 5 mm nitrite, but no inhibitor. Myxothiazol was seen to inhibit nitrite respiration very effectively, and no growth was observed (Fig. 6A). This confirms the effectiveness of myxothiazol as an inhibitor of the QCR in T. selenatis.
FIGURE 6Inhibition of respiration by myxothiazol, antimycin A, and HQNO.A, optical density of T. selenatis cultures grown with 10 mm selenate as electron acceptor in the absence and presence of inhibitors. Shown are no inhibitor (■), 10 μm myxothiazol (♦), 20 μm HQNO (●), and both 10 μm myxothiazol and 20 μm HQNO (▾). Positive control: T. selenatis grown with 5 mm nitrite as electron acceptor in the presence of 10 μm myxothiazol (○). B, optical density of T. selenatis cultures grown with 10 mm selenate as electron acceptor in the absence and presence of antimycin A. Shown are no inhibitor (■), 10 μm antimycin A (♦), and 20 μm antimycin A (●). C, specific growth rate (μmax) of T. selenatis cultures grown on increasing selenate concentration in the absence (■) and presence (●) of myxothiazol (10 μm). In all cases error bars represent standard deviation of n = 10 cultures.
T. selenatis was cultured under selenate respiratory conditions, with varying concentrations of selenate and 10 μm myxothiazol. Ten wells contained 10 mm selenate but lacked myxothiazol, as a positive control. The resulting growth curves (Fig. 6A) indicate that selenate respiration is partially inhibited by myxothiazol, suggesting that the QCR is involved in this process. The experiments were repeated in the presence of the inhibitor antimycin A (at both 10 and 20 μm), and the results (Fig. 6B) are consistent with the observations made using myxothiazol, strongly indicating the involvement of QCR in electron transport. Again, even at 20 μm antimycin, some growth is detected (final A600 nm ∼ 0.3 unit). Analysis of the specific growth rates (μmax) at increasing selenate concentration shows that Ks is unaffected (∼2.5 mm); however, in the presence of myxothiazol, a 3-fold decrease in μmax is recorded (no inhibitor = 0.32 ± 0.03 h−1; plus myxothiazol = 0.10 ± 0.02 h−1) (Fig. 6C). Therefore, it can be seen, from the growth curve analysis, that both myxothiazol and antimycin A retarded the growth rate of T. selenatis during selenate respiration but did not inhibit it completely.
Inhibition of Selenate Respiration by HQNO
Because myxothiazol has been shown to fully inhibit the QCR in nitrite reduction, but is only partially effective at inhibiting selenate respiration, HQNO was used as an inhibitor that affects a wider range of proteins with quinol dehydrogenase activity. HQNO was shown to inhibit growth completely with no significant change in culture A for ∼14 h post inoculation. Cultures maintained for 14–18 h showed a slight increase in culture A presumably due to the breakdown or evaporation of HQNO from the growth medium (Fig. 6A). By combining both inhibitors, HQNO and myxothiazol, there was no significant difference between the growth curves obtained with HQNO alone. The results (Fig. 6A) show that partial inhibition of selenate respiration is obtained by myxothiazol (and antimycin A), whereas complete inhibition of selenate respiration is obtained by HQNO. The combined data suggest that there is more than one route by which electrons can be transferred to selenate reductase: via the QCR, which is inhibited by myxothiazol and antimycin A, and via another type of quinol dehydrogenase, which is inhibited by HQNO.
DISCUSSION
The bioenergetics of selenate respiration in T. selenatis has remained unresolved since the isolation of T. selenatis by Macy and co-workers in 1993 (
). The respiration of selenate results in the sequential reduction of selenate to selenite to elemental selenium, yet it is only selenate that can function as a respiratory substrate. SER is a soluble molybdoenzyme composed of three subunits and is expressed from an operon of four genes (
), three for the structural components and the forth (SerD) a private chaperone for the functional assembly of SerA. No obvious membrane-bound component has been identified. Although others have speculated as to how SER receives electrons from the Q-pool and how this is coupled to the generation of PMF (
), biochemical evidence for the electron transfer pathway in T. selenatis has remained wanting. In the present study we have demonstrated that SER can receive electrons from a di-heme cytochrome of the cytc4 family and is driven by electron transfer via a branched pathway that includes electron transfer from QCR or QDH. The demonstration of an electron transfer pathway to a periplasmic molybdoezyme that utilizes a high potential di-heme cytochrome c carrier and electron flux through QCR is the first of its kind and unprecedented in other molybdo-oxidoreductase systems.
A number of periplasmic c-type cytochromes were up-regulated when T. selenatis was grown on selenate as the sole electron acceptor. During this study, two of these cytochromes have been purified and tested as potential redox partners to SER. A 24-kDa cytochrome, termed cytc-Ts4, was shown to donate electrons to SER and was characterized further. N-terminal sequence and liquid chromatography/tandem mass spectrometry analyses identified peptides as those resulting from members of the cytochrome c4 family. Attempts to use PCR and sequence the cytc-Ts4 gene using degenerate primers derived from other c4 homologues failed. Consequently, a holistic approach was adopted, and the entire genome of T. selenatis was sequenced using an Illumina sequencer. The complete cytc-Ts4 gene sequence was subsequently identified. A database search with the translated sequence of this gene produced a large number of similar proteins, most of which are classified as cytochromes c4, a group of di-heme cytochromes c characterized by high redox potentials, a low intensity asymmetrical α band around 550 nm and a low α/β peak ratio (
), but the discovery of a cytochrome c4 in T. roseopersicina, a facultative anaerobic photosynthetic bacterium, suggests that their roles are much more varied (
). Recently, a cytc4 has also been found to function as a physiological electron carrier in the photosynthetic electron transport chain in Rubrivivax gelatinosus (
). The two hemes in each case both have His-Met axial ligation, and the methionine residues identified as ligands are conserved throughout the cytochrome c4 group, and in the cytc-Ts4. The sequence similarity and high redox potential of the 24-kDa cytochrome suggest that we can assign cytc-Ts4 to the cytochrome c4 group. The involvement of a c4 cytochrome in anaerobic selenate respiration defines a new role for this class of cytochromes.
Of the other cytochromes up-regulated during growth of T. selenatis on selenate, a 6-kDa cytochrome was also purified (cytc-Ts7). No evidence was obtained to suggest that cytc-Ts7 could transfer electrons to SER. Furthermore, no direct reduction of selenate or selenite was observed. The function of cytc-Ts7 remains to be established. Chlorate reductase from I. dechloratans is also a member of the TAT-translocated type II molybdoenzymes and is closely related to SER. It has been shown recently that a 6-kDa cytochrome from I. dechloratans was readily oxidized in the presence of chlorate and cell homogenate containing chlorate reductase (
), suggesting that a 6-kDa cytochrome could be the redox partner for chlorate reductase. The authors suggest that the 6-kDa cytochrome might accept electrons from QCR, but no biochemical evidence for the involvement of QCR was reported. More recently, a cytochrome c gene located at the gene cluster for chlorate respiration in I. dechloratans has been cloned and overexpressed in E. coli but demonstrated not to function as an electron donor to purified chlorate reductase (
). It is interesting to speculate what the functions of the additional c-type cytochromes are in T. selenatis when grown on selenate. The mechanism by which selenite is reduced remains unknown. What is clear is that the reduction of selenite results in the formation of selenium nanospheres that are found in the surrounding medium. For the elemental selenium to form extracellular particles, it is considered likely that selenite is reduced either in the periplasm or on the outer cell surface, and as a consequence the other cytochromes identified might form an electron conduit to selenite reductase.
The involvement of QCR in the respiratory electron transfer chain of T. selenatis was suggested from the genome analysis (GenBankTM accession EU732596.1) and confirmed by complete inhibition of nitrite respiration by myxothiazol. T. selenatis expresses a cd1 nitrite reductase (also identified in the T. selenatis genome (GenBankTM accession AY078264.1), which commonly receives electrons from QCR via soluble intermediates (cytochromes/azurins). Selenate respiration was only partially inhibited (70% reduction in μmax) by the same concentration of myxothiazol or antimycin A and completely inhibited by HQNO. QCR binds two quinol molecules at the periplasmic face of the membrane and moves two electrons to the cytoplasmic face of the membrane, where it reduces one quinone molecule. This so-called Q-cycle translocates two positive charges per quinol oxidized. Given that myxothiazol and antimycin A are specific inhibitors of QCR, a redox loop coupling a Q-cycle mechanism to selenate reduction would seem plausible. HQNO is a less specific quinone analogue, which inhibits both QCR and other membrane-bound QDHs (
). HQNO alone, or in combination with myxothiazol, resulted in the complete inhibition of growth of T. selenatis. These data suggest that selenate reductase can accept electrons from more that one electron donor. A branched pathway that involves both QCR and QDH is proposed (Fig. 7). A QDH partner has not been identified, but analysis of the c-type cytochromes in membrane fractions from T. selenatis, grown on selenate, reveals that cytochromes of ∼30 and ∼20 kDa are present, but attempts to purify these components have so far proved unsuccessful (data not shown). A napC/nirT homologue (
The ability to utilize a QCR-dependent bioenergetic pathway is facilitated by the high redox potential of the selenate/selenite couple (Em +475 mV). This allows the selenate reductase to function by using a high potential cytochrome c as the electron donor. The interesting feature of SER thus relates to the point of electron entry and the unusual b-heme in SerC. The determination of the redox potential of the SerC heme has allowed us to add to our knowledge of the redox centers of the selenate reductase and is the first periplasmic reductase of the Type II molybdoenzymes to have the potential of its b-heme determined. Supplemental Fig. S2 shows the known midpoint potentials of all redox centers within SER, from the redox potential of the substrate to the potential of b-heme subunit. The value of Em = +234 ± 10 mV for the SerC b-heme is very similar to that determined for EbdC (+254 mV). EBDH functions by withdrawing electrons from ethylbenzene, where the b-heme functions as an electron attractant. Conceptually, to convert from a dehydrogenase function to a reductase function, one might assume that the potential of the corresponding b-heme in the type II molybdoreductases might be retuned to function at a much lower redox-potential, thus facilitating electron donation. This appears not to be the case for selenate reductase. It is of interest to compare SerC/EbdC cytochromes to the membrane-bound cytochrome associated with the respiratory nitrate reductase (NarI). NarI contains two b-hemes, but both have much lower redox potential (+120 and +20 mV). Modeling the three-dimensional structure of SerC, using the x-ray coordinates of EbdC from A. aromaticum, predicts that the b-heme in SerC is also coordinated by methionine (Met-138) and lysine (Lys-228) ligands (Fig. 8). In NarI, both hemes have iron with bis-histidinyl coordination, and this could contribute to the differences in redox potential. The sulfur of methionine residues is a good electron acceptor and could contribute to stabilizing the reduced state of the heme, therefore significantly raising the midpoint potential.
FIGURE 8b-Heme coordination in SerC.A, predicted structure of SerC modeled upon the crystal structure of EBDH γ-subunit. The superposition shows SerC in green and the EbdC in cyan. B, detail of heme pocket showing ligands. Homology modeling was performed using the homology tools in the software package MOETM. The sequence of T. selenatis SerC was aligned and modeled against the structure of EBDH γ-subunit from A. aromaticum (PDB ID 2IVFC). The sequence was modeled from residues 27–239. The modeling was performed using the protein template alone, and additionally using the environment of the template structure with the heme ligand. An ensemble of 30 intermediate models was created, and the best intermediate model as minimized using the CHARM22 force field to a root mean square gradient of 0.01. The quality of the models was assessed by PROCHECK (
Evidently the combination of a high potential b-heme in SerC and the unusual high potential redox couple of the substrate (selenate/selenite) has provided a system that can readily accept electrons from cyt-c and QCR redox partners. By using these redox partners, selenate reduction can be linked to a Q-cycle mechanism. Because selenate reductase functions on the positive side of the cytoplasmic membrane, the reduction of selenate to selenite will consume two pumped protons. If SER was to receive electrons only from a membrane-bound QDH, then only two protons would be released to the periplasm, resulting in no net gain of proton electro-chemical gradient (0q+/2e−). Because the reduction of selenite to selenium does not support growth, selenate reduction via a QDH could only be energy conserving if electron input into the Q-pool was proton-motive (i.e. via an NADH dehydrogenase, a formate dehydrogenase, or hydrogenase). By using the Q-cycle mechanism a net gain of 2q+/2e− of proton electrochemical gradient can be maintained irrespective of the electron donor.
Finally, the Q-cycle-coupling mechanism presented for selenate reduction in T. selenatis, while a novel mechanism for electron transport to molybdoenzymes in bacteria, might also be utilized in the Archaea (
). Sequence analysis has suggested that the clade of TAT-translocated type II molybdoenzymes, of which SER is a member, also contains a number of unusual respiratory nitrate reductases (pNAR). These reductases all possess a TAT leader peptide and are translocated to the positive side of the membrane. The pNAR enzymes are distributed among the Archaea, and examples have been identified in the hyperthermophiles and halophiles. Interestingly, those from halophiles (Haloferax mediterranei and Haloarcula marismortui) lack a NarI-like component so the mechanism of electron transport is not obvious. However, genetic analysis reveals the presence of a sequence that encodes a protein (NarC) similar to the di-heme subunits of QCR. Adjacent to the narC gene is a gene that encodes NarB, which is predicted to encode a Rieske iron-sulfur protein. The combination of NarB and NarC in H. mediterranei has led to the suggestion that nitrate respiration in some organisms might be driven by a Q-cycle mechanism (
). By coupling a Q-cycle mechanism to nitrate reduction would provide a pNAR system that is bioenergetically equivalent to the well characterized NAR (or nNAR) systems in bacteria.
Acknowledgments
We thank Ann Reilly (University of East Anglia), Iona Knight (University of Exeter), and Dr. Carys Watts and Dr. Helen Ridley (University of Newcastle) for technical support and useful discussions. We also thank Dr. Kirsty Line (University of Exeter) for providing structural models.