Phenothiazine derivative-accelerated microbial extracellular electron transfer in bioelectrochemical system.

In bioelectrochemical system (BES) the extracellular electron transfer (EET) from bacteria to anode electrode is recognized as a crucial step that governs the anodic reaction efficiency. Here, we report a novel approach to substantially enhance the microbial EET by immobilization of a small active phenothiazine derivative, methylene blue, on electrode surface. A comparison of the currents generated by Shewanella oneidensis MR-1 and its mutants as well as the electrochemical analytical results reveal that the accelerated EET was attributed to enhanced interactions between the bacterial outer-membrane cytochromes and the immobilized methylene blue. A further investigation into the process using in situ Raman spectro-electrochemical method coupled with density functional theory calculations demonstrates that the electron shuttling was achieved through the change of the molecule conformation of phenothiazine in the redox process. These results offer useful information for engineering BES.

In conventional waste treatment processes, not only substantial energy in organic matter is wasted but also extra energy should be invested1. Bioelectrochemical systems (BES) is recognized as an energy-saving alternative because it allows direct recovery of energy from waste streams like sewage12. From filling knowledge gaps in earth's biogeochemical cycles to powering different devices such as nano-sensors and detectors, BES offers a versatile platform for various technologies3456. One common feature of all BES is a microbially catalyzed oxidation of substrate and subsequently transport of the generated electrons to the anode. However, the extracellular electron transfer (EET) from bacteria to anode electrode is usually a limiting step in this process and is recognized to be a crucial step in engineering BES7. In an effort to accelerate the anodic EET, one of the useful approaches is to use nanostructured chemical catalysts8910.

Three distinct EET pathways, i.e., through cell-surface proteins, mediator compounds (such as flavins), and cellular pili, also named nanowires, have been proposed in BESs. For some bacteria such as Shewanella, usually more than one type of EET pathways are involved111213. These EET pathways inspire us in designing more powerful materials for BES in a biomimetic way. In our previous study, we fabricated a network-structured graphene oxide, which acted as nanowires to significantly enhance the anodic EET14. In addition, in light of a positive role of electron shuttles in EET but a unstable nature of these molecule in solution1215, immobilization of redox-active molecule on the electrode has also been adopted as an efficient way to promote microbial EET. For example, immobilization of 9,10-anthraquinone-2,6-disulfate on anode electrode using polyethyleneimine as a binder has been reported to improved microbial current generation16. However, one shortcoming of polyethyleneimine is the poor conductivity, which can to some extent limit the EET process. Furthermore, despite of intensive studies in the role of redox mediators in microbial reduction of electrodes, there has not been a conclusive consensus whether electrode modification with synthetic redox mediators could promote EET.

Phenazines are aromatic compounds with a dibenzo annulated azine structure. In their derivatives, a methyl and/or amino group is attached to the benzene rings, where one N can be substituted by S (phenothiazine) for example in the case of methylene blue (MB)17. These structures make them easy to be electropolymerized on electrode surface18. Like other electron shuttles, most phenothiazine derivatives have sufficiently-low reduction potentials ranging from −180 to 100 mV under usual physiological and environmental conditions. Thus, the redox reactions between reduced phenazines or thiazines and terminal electron acceptors are thermodynamically feasible19. These features make them ideal candidates for electrode modification to enhance the microbial EET at BES anode.

Here, we attempt to accelerate anodic microbial EET in a BES by electropolymerization of small active phenothiazine derivatives on a carbon electrode. To achieve this, Shewanella was chosen as a model anodic microbe, and the anode electrode was electropolymerized with exogenous phenothiazine. The interactions between Shewanella and the as-prepared anode electrode were investigated by recoding the microbial current generation. In addition, the role of immobilized phenothiazine in EET was elucidated by coupling in situ Raman spectral-electrochemical analysis with density functional theory (DFT) to probe the conformation variations at a molecular level.

Results

Electron transfer at different potentials

Figure 1a shows the microbial current generated at a potential-poised electrode when an active phenothiazine of small molecules (either Toluidine blue O, Thionine, or MB) was used as the electron mediator. A BES without dosing any of the three phenothiazine derivatives was used as the control, which showed very low current density. A higher current density was produced with MB than with Toluidine blue O or Thionine. The current density peaked at 0.81 A/m2 at the 10th hour for the MB-dosed BES, but only reached 0.65 and 0.21 A/m2 for systems with Toluidine blue O and Thionine, respectively. By plotting the average of the peak current density of the three phenothiazine derivatives vs. their respective redox potential at the 10th hour (determined at pH 7.0 vs SHE; Figure 1b)20, it was found that the phenothiazine derivatives with a lower redox potential tended to have a higher peak current density (Figure 1b). MB had a lowest redox potential of 0.011 V, and also possessed the largest peak current density among the three phenothiazine derivatives. Thus, in the sequential experiments, MB was chosen as the model phenothiazine derivative to modify the electrode.

Figure 1

Microbial EET recorded by the potential-poised electrodes in the presence of phenothiazine (a); plot of the maximum current densities of the three phenothiazine derivatives vs. their respective redox potentials (b).

Enhanced microbial current generation by MB film

When a raw carbon paper was used as the control anode, the current density increased to 0.09 A/m2 (Figure 2a). In comparison, with the MB-modified carbon paper, the current density drastically increased and peaked at 1.11 A/m2, indicating a distinctive improvement in the microbial electron transfer (Figure 2a). To validate whether such a result was repeatable, each experiment was repeated. A similar current curve was obtained (Figure 2S). Notably, in our study only planar electrode was used, the modification of MB did not increase the available electrode surface area, and the electrode was poised at a relative low potential (0.1 V vs. Ag/AgCl) compared with the previous reports (from 0.12 to 0.2 V vs. Ag/AgCl; Table S1). Nevertheless, the obtained current density here was comparable to those with porous nano-structured electrodes reported in other studies (1.0–1.8 A/m2)1421.

Figure 2

Microbial EET at the MB-immobilized electrode with Shewanella oneidensis MR-1 wild type.

(a) ΔOmcA/MtrC mutant (b) and ΔMtrB mutant (c).

It has been widely recognized that the c-type cytochromes in bacterial outer membrane play a significant role in direct electron transfer from cell surface to solid minerals or electrodes because of their unique location on cellular surface in the respiratory electron-transfer process. To investigate whether the interactions between MB and c-type cytochromes are responsible for the promoted electron transfer at the bacteria/electrode interface, S. oneidensis MR-1 mutants (OmcA/MtrC) were used. The current densities were 0.01 and 0.11 A/m2 only for the control and the MB-modified electrodes, respectively (Figure 2b). Deletion of specific cytochromes significantly decreased the microbial current, indicating that both MtrC and OmcA played a key role in the electron transfer between the bacteria and MB interface. The experiment with mutant ΔMtrB further confirms the above results (Figure 2c).

It is well known that S. oneidensis MR-1 is able to excrete compounds such as flavins to shuttle electrons between cells and solid electron acceptor12. The CV results show that the BES with MB had a substantially lower redox current response to flavins than the control (Figure 3). The performance of the MB-modified electrode at different flavins concentrations in a range of 0–1.0 μM was further evaluated (Figure S3). The peak current with 0.5-μM flaving exhibited the highest value of 1.29 A/m2, which was slightly higher than that without flavin dosing (1.18 A/m2). This suggests that the MB-modified electrode did not promote the flavin-based electron transfer. This result is inconsistent with the enhanced current using the MB-modified electrode (Figure 2a). Therefore, the above results further demonstrate that the direct contact between MB and the cells was the main mechanism for the enhanced electron transfer.

Figure 3

CV of 5-μM flavin in fresh mineral medium solution with the dose 0.1 M KCl at scanning rate of 0.05 V/s using the bactieria-free glassy carbon electrode and the MB-modified glassy carbon electrode.

Raman spectroelectrochemical analysis on MB-mediated electron transfer

A higher net intensity of spectra was observed at a higher electrode potential (Figure 4), where MB was in its oxidized form. The number and position of the bands, as well as their relative intensities, changed substantially with the electrode potential. Obviously, the MB layer underwent redox changes with the varying electrode potentials.

Figure 4

Potential-dependent Raman spectra of MB on Au electrode.

Because of the complexity of this process, DFT calculations were performed to unambiguously identify the corresponding surface species. The detailed information about the calculated results is provided in the Supporting Information (SI). Strong Raman bands located within the region of 1609–1339 cm−1 belong to skeleton stretch vibrations of MB in the oxidized form according to the calculated results. An intense band centered at 1620 cm−1 for high electrode potentials corresponds primarily to aromatic ring stretching vibrations of short C3-N15, C12-N16 bonds, and aromatic ring C-C in-plane stretching vibrations coupled with scissor bending vibrations of C-H in the -CH3 group. A strong band at 1444 cm−1, observed for an oxidized form of MB was assigned to ring C-C stretching vibrations or asymmetric C-N stretches22. Calculations have shown that this mode composed primarily from bending deformation vibrations of H28-C17-H29 and H33-C19-H34 bonds, scissor bending vibrations of H30-C18-H32 and H36-C20-H37 bonds, asymmetric stretches vibration of C6-N7-C8 bond, and skeleton stretch vibrations of Rings I and III (Table S2).

A shift of electrode potential to a lower value, especially lower than −0.3 V, resulted in a decrease in the intensity at 1620 cm−1 and 1434 cm−1 band length. According to the DFT calculations, the length of C6-N7, N7-C8, C5-S10, and S10-C9 bonds would increase with a lowering potential (Table 1). This leads to the expansion of Ring II, which is responsible for a decrease in the intensity for 1620 cm−1 band, while the rock vibrations of C-H bonds in Rings I and III coupled with scissor bending vibrations of C-H in the -CH3 group results in a decrease in intensity of 1434 cm−1 bond. With a decrease in potential, the intensity at 1393 cm−1 increased and peaked at around −0.15 V. By referring to Table S3, it could be concluded that this band is associated with rock vibrations of N7-H39 bond, and the scissor bending vibrations of C6-N7-C8 bond. Therefore, the conformation changes of MB in the redox process can be easily probed according to the relative intensity variations of these characteristic bands.

Table 1

Geometry Parameters of Optimized Molecular Structures of MB in Its Oxidized and Reduced Forms

Bond	Length (oxidized form) (Å)	Length (reduced form) (Å)
C3-N15	1.358	1.398
C12-N16	1.358	1.399
C6-N7	1.339	1.408
N7-C8	1.339	1.408
C5-S10	1.747	1.783
S10-C9	1.747	1.783
Bond	Angle (oxidized form) (°)	Angle (reduced form) (°)
C6-N7-C8	123.171	120.983
C5-S10-C9	103.218	99.519

Discussion

In this work, a small-molecular-weight active phenothiazine derivative, MB, was immobilized on electrode surface to substantially enhance the microbial EET, and the interactions between the bacterial outer-membrane cytochromes and such a derivative were probed using in situ Raman spectro-electrochemical method coupled with DFT calculations for the first time.

Phenothiazine have multiple forms with different redox properties23, which influenced the electron transfer capacity. In this study, we attempted to explore the functions of a series of phenothiazine derivatives in the microbial EET of S. oneidensis MR-1 by directly recording the current. The results clearly show that the maximum microbial current generation was closely related to the redox potentials of the phenothiazine derivatives. This is consistent with previous studies on microbial reduction of lepidocrocite and azo dyes24. This also implies that there is a correlation between the rate of mediated redox reactions and the reduction potential of electron transfer mediators2425. With a consideration of different substrates and electron acceptors, mediators are chosen by their tunable redox potentials in the aqueous solution. Thermodynamically, substrates have their standard potential values, e.g., −0.29 V for CO2/acetate, −0.19 V for pyruvate/lactate26. The standard potential of the mediators should be as close to that of the substrate as possible in order to obtain a maximum energy output. Thus, for a specific electron flow route, it is possible to lower the biological energy acquisition and increase the energy conversion efficiency by choosing an appropriate mediator, such as MB in this study. Generally, a bacterium needs a minimum free energy of about −20 kJ/mol to exploit the free energy change in a reaction27. When this thermodynamic requirement is met, phenothiazine derivatives with a lower redox potential tend to have a greater electron mediation capacity in the subsequent electron transfer process from electron shuttle to electron acceptor (i.e., electrode).

Electron transfer mediators could assist in microbial energy metabolism by facilitating electron transfer between microbes and electron-donating or accepting substances. Our study reports, for the first time, the direct interactions between c-type cytochromes of electricity-generating bacteria and phenothiazine derivative molecules immobilized on a carbon electrode. The conformation shift of phenothiazine in the redox reaction process was revealed by a combined use of Raman spectroelectrochemical approach and DFT calculation. The results clearly show a distinct improvement of microbial current generation at the electrode resulted from the well communication between c-type cytochromes of S. oneidensis MR-1 and MB.

From a molecular level point of view, elucidation of the redox process of phenothiazine when it interacts with bacteria in the anode is one important fundamental challenge for the mediated electron transfer. The oxidation-reduction reaction of MB used in this study is a two-electron reaction23. The uptake and release of electrons lead to the changes in conformation of MB molecule. The structure variations of MB at the electrode/liquid interface probed by the in-situ Raman approach coupled with DFT calculation in electrochemical oxidation and reduction, reveals its role in the electron shuttling between bacteria and solid electron acceptors. In the oxidation process, two electrons are released to electrode, resulting in the shrink of active Ring II in MB molecule; when it is reduced by bacteria, MB uptakes two electrons and one proton to form a new N-H bond (Figure S4). This is responsible for the expansion of Ring II of MB (Table 1).

For S. oneidensis MR-1, the metal respiratory pathway (Mtr) is recognized as the most important anaerobic respiratory pathway to transfer electron from the cytoplasmic membrane to extracellular electron acceptor28. This pathway is also associated with the electron transfer from cell surface to MFC electrode through direct transfer by cytochromes or via indirect transfer by flavins29. Both MtrC and OmcA are located on the cellular outer membrane surface, and they can be directly involved in electron transfer to an extracellular solid electron acceptor, such as Fe(III) minerals30, or to extracellular electron shuttles12. ΔOmcA/MtrC mutant is severely impaired in the sustained electron transfer, but it could attach to electrodes in a similar manner31. Porphyrin iron is the structure centre of c-type cytochromes, which play an important role in the redox process (Figure 5). Therefore, the electron transfer from c-type cytochromes to phenothiazine electron shuttles was reasonably substituted to the interactions between porphyrin iron and MB molecule. The proposed mechanism for the electron transfer process is illustrated in Figure 5a and b. Porphyrin iron (Fe2+) as the electron donor approaches to the electron acceptor (MB oxide), and the distance between Fe2+ of porphyrin iron and S atom of Ring II in MB oxide is 8.147 Å, while the value in the porphyrin iron (Fe3+)-MB reduction system decreases to 5.175 Å. This indicates that the redox reaction occurs in the surrounding of porphyrin with the range of 5–8 Å, which is the core structure of c-type cytochromes.

Figure 5

Optimized molecular structures of porphyrin iron in reduced form and MB in oxidized form (a) front view and (c) side view, and porphyrin iron in oxidized form and MB in reduced form (b) front view and (d) side view: gray: carbon atom; white: hydrogen atom; blue: nitrogen atom; yellow: sulfur atom; red: oxygen atom; purple: iron atom.

The negative energy change (ΔE = −0.594 eV) of the electron transfer (Table S4) reveals that the porphyrin iron (Fe3+)-MB reduction system (Figure 5b and d) is more stable than the porphyrin iron (Fe2+)-MB oxide system (Figure 5a and c) in the aqueous solution. Therefore, the oxidized form of the MB molecule could easily take electrons from the c-type cytochromes, and transfer them to electron acceptors with a higher potential.

The redox mediators excreted by S. oneidensis MR-1 are usually found at a level of μM or nM concentration32. It could accelerate electron transfer and catalyze electron transfer at a lower applied potential31. However, the redox peak current at the MB-modified electrode was lower than that at the control electrode (Figure 3). This implies that the direct contact between MB and the biofilm should be the main mechanism for the enhanced microbial EET at bacteria/solid electron acceptor interface.

In conclusion, the immobilized phenothiazine on solid electrode was found to enhance microbial EET in BES. The in situ Raman spectro-electrochemical analysis clearly demonstrates that the mediated EET process was achieved through changes in molecular conformation of phenothiazine. These findings are beneficial for designing new electrode interface to improve the current density of BES. In addition, they are also useful to better understand the microbial EET mechanism, which is crucial to develop and apply BES, biosensors or related bioelectronic devices. These findings may also have general implications in the biogeochemical cycles, biocorrosion and bioremediation.

Methods

Bacterium and cultivation

The wild type S. oneidensis MR-1 and its mutant strains, kindly provided by Prof. K. H. Nealson at the University of Southern California33, were grown from a frozen stock, first aerobically in 100-ml Luria-Bertani broth, a solution of 10 g/l tryptone, 5 g/l NaCl, and 5 g/l yeast extract, at 30°C for 12 h. Subsequently, the culture of 0.5 ml was anaerobically incubated in a mineral medium (40 ml) containing (per liter) 10 mM Hepes, 0.46 g of NH4Cl, 0.225 g of K2HPO4, 0.225 g of KH2PO4, 0.117 g of MgSO4·7H2O, and 0.225 g of (NH4)2SO4. Prior to autoclaving, 10 ml of a mineral mix (containing per liter 1.5 g of nitrilotriacetic acid, 0.1 g of MnCl2·4H2O, 0.3 g of FeSO4·7H2O, 0.17 g of CoCl2·6H2O, 0.1 g of ZnCl2, 0.04 g of CuSO4·5H2O, 0.005 g of AlK(SO4)2·12H2O, 0.005 g of H3BO3, 0.09 g of Na2MoO4, 0.12 g of NiCl2, 0.02 g of NaWO4·2H2O, and 0.10 g of Na2SeO4) were added31. Medium was sparged with N2 gas, adjusted to pH 7.2, and autoclaved.

All solutions were prepared from reagent grade chemicals without further purification. Toluidine blue O, Thionine, and MB were purchased from Aladdin Reagent Co., China.

Electropolymerization of MB

The electropolymerization of MB on carbon paper electrode (2 × 3 cm2, 190 μm thickness, Toray Co., Japan) was carried out using cyclic voltammetry (CV) in a potential range of −1.0 to 1.2 V at a sweep rate of 0.1 V/s in Britton-Robinson (B-R) buffer solution34 (pH 7.0) with 30 cycles (Figure S1). The MB concentration was usually 0.5 mM. After the MB film was formed on the carbon paper surface, the electrode was rinsed thoroughly with distilled water.

Three-electrode experiments

A single-chamber, three-electrode system was used, in which a Ag/AgCl (3 M KCl) and a Pt wire were used as the reference electrode and counter electrode respectively. For the selection of the phenothiazine derivatives, they were chosen according to their tunable redox potentials in the aqueous solution. Three different mediators were chosen: Toluidine blue O (redox potential, E0 = 0.027 V vs. Standard Hydrogen Electrode, SHE, pH = 7.0), Thionine (E0 = 0.060 V), and MB (E0 = 0.011 V)20. The above mediator solutions were individually added to the mineral medium containing 10 mM lactate after filtration sterilization to make the final concentration of 1 μM. Carbon paper (2 × 3 cm2, 190 μm thickness, Toray Co., Japan) was used as the working electrode.

To evaluate the enhancement of microbial EET by the immobilized MB film, carbon paper electropolymerized with MB was used as the working electrode, while raw carbon paper was used in the control tests. The freshly prepared cell suspension was inoculated into the reactor at a constant voltage of 0.1 V (vs. Ag/AgCl, 3 M KCl) held with a model CHI 1030A electrochemical workstation (Chenhua Instrument Co., China). After the working electrodes were poised at 0.1 V for at least 12 h and the biofilm was formed on the electrodes, 1 mL of concentrated lactate solution was injected to the three-electrode system to make a concentration 10 mM. Each experiment was conducted at least two times.

Electrochemical characterization

CV measurements of flavin at the MB-modified electrode were conducted using a conventional three-electrode cell with a platinum wire as the counter electrode and a Ag/AgCl (3 M KCl) as reference in a CHI 660C Electrochemical Workstation (Chenhua Instrument Co., China). MB-modified glassy carbon electrodes (3-mm diameter) were used as the working electrodes. Before use, the carbon electrodes were carefully polished to a mirror finish with 1.0, 0.3, and 0.05 μm alumina slurries successively. The electrolyte was 5 μM flavin in mineral medium solution with 0.1 M KCl added. The scanning rate was 0.05 V/s. The solution was de-aerated by passing highly pure nitrogen for 20 min before the electrochemical experiments, and a continuous flow of nitrogen was maintained over the sample solution during experiments.

Raman in-situ spectro-electrochemical experiments were conducted in a closed cylindrical three-electrode cell with a platinum wire counter electrode, and an Ag/AgCl reference electrode. The cell was filled with the B-R buffer solution. To obtain enhanced Raman signal, Au disc electrode (2 mm in diameter) electropolymerized with MB was used as the working electrode. Raman spectra were obtained using a LABRAM-HR Raman spectrometer with an excitation wavelength of 514.5 nm generated by an Ar+ laser. The beam was focused to a spot of ca. 1 mm2 area on the electrode surface. The experiments were carried out in 90° scattering geometry. The integration time was 10 s. In order to reduce photo- or thermo-effects and a possible degradation of MB film by the incident light, the spectroelectrochemical cell and the working electrode were periodically moved with respect to the laser beam. All potential values given are referred to Ag/AgCl reference electrode.

DFT simulations

The minimum-energy geometry structures of the oxidized and reduced methylene blue were determined by DFT computation. For the calculations, an all-electron method, which is within the Perdew-Wang 91 (PW91) form of generalized gradient approximation3536 for the exchange-correlation term, is used as implemented in the DMol3 code3738. This method adopts double precision numerical basis sets that take into account p polarization (DNP). The sizes of the DNP basis sets are comparable with the 6–31G** basis and are regarded as more accurate than Gaussian basis sets of the same size39. The DNP basis sets were used to describe the valence electrons, while all-electron core treatment was utilized to describe the core electrons. A spin-polarized scheme was used to deal with the open-shell systems. The energy in each geometry optimization cycle was converged to within 1 × 10−5 Hartree with a maximum displacement and force of 5 × 10−3 Å and 2 × 10−3 Hartree/Å, respectively. Frequency calculations were performed at the same level of theory to characterize the stationary points on the potential surface and to obtain Gibbs free energies, which were calculated at a standard pressure of 1 atm. The calculated vibrational frequencies were compared with the experimental data.

Calculations on the interactions between porphyrin iron and MB

The interactions between porphyrin iron and MB in the process of electron transfer were carried out for all structures using the GULP code4041. The GULP geometry optimization is based on reducing the magnitude of calculated forces and stresses until they become smaller than defined convergence tolerances. The energy tolerance and gradient norm tolerance were smaller than 0.00001 eV and 0.001 eV/Å, respectively. The geometry structure of electron transfer systems was refined to obtain a stable structure with the minimized total energy. Then, the interaction energy of electron transfer and the distance between the active sites of porphyrin iron and MB could be analyzed. The Dreiding forcefield42 is a purely diagonal forcefield with harmonic valence terms and a cosine-Fourier expansion torsion term. The umbrella functional form was used for inversions, which are defined according to the Wilson out-of-plane definition. The van der Waals interactions are described by the Lennard-Jones potential. Electrostatic interactions are described by atomic monopoles and a screened (distance-dependent) Coulombic term. The COSMO implementation in GULP41 is a continuum solvation model, in which the solute molecule forms a cavity within the dielectric continuum of permittivity. The charge distribution of porphyrin iron and MB polarizes the dielectric medium. The response of the dielectric medium was described by the generation of the polarization charges on the cavity surface. In the COSMO solvation model, ‘water' was selected as the electron transfer environment for the calculations.

Author Contributions

X.W.L., X.F.S., G.P.S. and H.Q.Y. designed the experiments; J.J.C. conducted the DFT calculations and molecular dynamics simulations; X.W.L., X.F.S., Y.X.H., J.F.X., W.W.L., Y.Y.Z. and R.L. conducted the experiments; X.W.L., G.P.S., F. Z. and H.Q.Y. contributed to the planning and coordination of the project; X.W.L., G.P.S., F. Z., J.J.C. and H.Q.Y. wrote and edited the manuscript. All authors contributed to discussions about the results and the manuscript.