X-ray structure of the direct electron transfer-type FAD glucose dehydrogenase catalytic subunit complexed with a hitchhiker protein

The X-ray structure of the catalytic subunit of Burkholderia cepacia FAD glucose dehydrogenase complexed with a hitchhiker protein was determined as a representative molecule of direct electron transfer-type FAD-dependent dehydrogenase complexes. The 3Fe–4S cluster is present at the surface of the catalytic subunit and serves in the intramolecular and intermolecular electron transfer from FAD to the electron-transfer subunit.


Introduction
Various sugar oxidoreductases (dehydrogenases) have been reported to be inherently capable of direct electron transfer to electrodes composed of carbon materials or to gold electrodes. These dehydrogenases harbor an electron-transfer domain or subunit, together with a catalytic domain or subunit. The catalytic domains or subunits, which are responsible for catalyzing sugar oxidation, are categorized by their cofactor: flavin adenine dinucleotide (FAD) or pyrroloquinoline quinone (PQQ). The electron-transfer domains or subunits, which are responsible for transferring electrons to the external electron acceptor, are also categorized by the type of heme (heme b or c) that is present in the electron-transfer domain or subunit.
One of the representative groups of direct electron transfertype (DET-type) dehydrogenases consists of cellobiose dehydrogenases (CDHs), which are composed of a catalytic domain harboring FAD and a heme b-type electron-transfer domain. In CDHs, open and closed states were identified by ISSN 2059-7983 the approach of the cytochrome domain to the catalytic domain via a flexible linker (Tan et al., 2015).
The other representative protein group that is capable of direct electron transfer consists of FAD-dependent dehydrogenase complexes. These complexes are composed of a catalytic subunit with FAD, an electron-transfer subunit containing three heme c moieties and a small subunit. The isolation, characterization, bioelectrochemical studies and application of several bacterial FAD-dependent dehydrogenase complexes have been reported, including bacterial glucose dehydrogenase (FADGDH; Inose et al., 2003;Sode et al., 1996;Tsuya et al., 2006;Yamazaki et al., 1999;Yamaoka & Sode, 2007;Yamaoka et al., 2008), fructose dehydrogenase (FDH; Ameyama et al., 1981;Kawai et al., 2013;), 2-ketod-gluconate dehydrogenase (KGDH; Kataoka et al., 2015;Shinagawa et al., 1981) and sorbitol dehydrogenase (Toyama et al., 2005). These FAD-dependent dehydrogenase complexes have the potential to directly transfer electrons to an electrode because of the presence of the heme c subunit. However, no structural information is currently available for any subunits from DET-type FAD-dependent dehydrogenase complexes.
Our research group has been studying a representative FAD-dependent dehydrogenase complex (FADGDH) derived from Burkholderia cepacia SM4 (BcGDH). BcGDH comprises three distinct subunits: the catalytic subunit (-subunit), which contains an FAD cofactor in its redox center, shows catalytic activity and oxidizes the first hydroxyl group of glucose, the small subunit (-subunit), a hitchhiker protein of the bacterial TAT secretion system that is necessary for the proper folding and secretion of the -subunit (Yamaoka et al., 2004), and the membrane-bound subunit with three heme c moieties (-subunit) that is responsible for the transfer of electrons between the active-site cofactor and external electron acceptors. Owing to the presence of the -subunit, BcGDH is capable of transferring electrons directly to an electrode, making it an ideal molecule for glucose sensors and applications in a variety of biomedical devices Yamashita et al., 2018). In addition, the BcGDH complex, which is BcGDH lacking the -subunit, also exhibits dye-mediated glucose dehydrogenase activity (Inose et al., 2003). Recently, based on biochemical analyses and electron paramagnetic resonance spectroscopy, we reported the presence of a 3Fe-4S cluster in the catalytic subunit (Shiota et al., 2016). The 3Fe-4S cluster is located in the cysteine-rich region, which is conserved in the catalytic subunits of previously reported FAD-dependent dehydrogenase complexes. The 3Fe-4S cluster is responsible for electron transfer from FAD (intramolecular) to the multiheme c subunit (intermolecular), which is the key position for understanding the features of this group of enzymes that are capable of direct electron transfer.
Another notable feature of the enzymes in the FADdependent dehydrogenase complex is the presence of a small subunit that is essential for the functional expression of the FAD-harboring catalytic subunit. Considering their primary structure, particularly the signal sequences necessary for secretion, these small subunits are predicted to be hitchhiker proteins that are needed for secretion of the catalytic subunit into the periplasmic space (Yamaoka et al., 2004). However, no structural information is available for any types of hitchhiker proteins or their complexes with targeted proteins.
In this study, we determined the X-ray structure of BcGDH, consisting of the BcGDH catalytic subunit complexed with the small (hitchhiker) subunit. The structure of the FAD-binding catalytic subunit was similar to those of several glucose-methanol-choline (GMC) oxidoreductases. The catalytic sites of BcGDH were conserved compared with other GMC oxidoreductases. In the structure of BcGDH, the 3Fe-4S cluster was located at the surface of the catalytic subunit. The structure of the complex of the catalytic subunit with the small subunit revealed that these two molecules were connected through disulfide bonds and hydrophobic interactions. Site-directed mutagenesis studies were performed to elucidate the role of the disulfide bond. The structural similarities to other FAD-dependent dehydrogenase complexes and to fumarate reductase are also discussed.

Recombinant expression of BcGDHca
The structural genes encoding the catalytic () and small () GDH subunits of B. cepacia sp. SM4 (FERMBP-7306) were subcloned into the high-expression vector pTrc99A with a His tag at the C-terminus of the -subunit. The constructed plasmid (designated pTrc-His) was transformed into the bacterial host Escherichia coli BL21 (DE3) for expression. In the present study, a complex of the -subunit and His-tagged subunit was expressed as recombinant wild-type BcGDH. The -subunit consists of 168 amino acids (18 kDa), including 47 amino acids of the signal peptide at the N-terminal region, and the mature -subunit contains 121 amino acids (13 kDa). The subunit consists of 539 amino acids and was produced as a catalytic domain of 60 kDa.
Transformed E. coli were cultured in 500 ml conical flasks containing 100 ml ZYP-5052 medium (Studier, 2005) in a rotary shaker at 293 K for 48 h.
The cells were harvested by centrifugation and then resuspended in 20 mM sodium phosphate buffer containing 20 mM imidazole and 0.5 M NaCl pH 7.0. After resuspension, the cells were disrupted with a French press. The lysate was centrifuged at 10 000g for 15 min at 277 K to remove the the insoluble fraction consisting of cell debris and inclusion bodies. The resulting supernatant, which was designated the crude extract, was purified by FPLC.

Enzyme purification
The recombinant BcGDH complex was purified using nickel-chelate chromatography and cation-exchange research papers chromatography. The crude extract was loaded onto a HisTrap HP column (1 ml; GE Healthcare Life Sciences, Uppsala, Sweden) that had been equilibrated with 20 mM sodium phosphate buffer containing 20 mM imidazole and 0.5 M NaCl pH 7.0 and was washed with the same buffer. GDH was then eluted with ten column volumes of a stepwise imidazole gradient (70, 380 and 500 mM imidazole in 20 mM sodium phosphate buffer and 0.5 M NaCl pH 7.0) at a rate of 1 ml min À1 for each step. The fractions with the highest activities were pooled and dialyzed overnight in 10 mM potassium phosphate buffer pH 6.0.
The pooled fractions were subsequently loaded onto a Resource S column (5 ml; GE Healthcare, Little Chalfont, England) that had been equilibrated with 10 mM potassium phosphate buffer pH 6.0 and were washed with the same buffer. GDH was eluted with 20 column volumes of a linear NaCl gradient (0-1 M NaCl in 10 mM potassium phosphate buffer pH 6.0) at a rate of 5 ml min À1 . The purified enzyme was concentrated to 8.9 mg ml À1 and the buffer was exchanged to Milli-Q water using Amicon Ultra-15 (nominal molecular-weight limit 3000; Merck Millipore, Carrigtwohill, Ireland). The protein concentrations were measured using a DC Protein Assay kit (Bio-Rad Laboratories, Hercules, California, USA).

Site-directed mutagenesis
Site-directed mutagenesis of the target amino acids (Cys213 in the -subunit and Cys152 in the -subunit) was accomplished using the QuikChange mutagenesis kit (Agilent, Santa Clara, California, USA) according to the manufacturer's instructions. All mutations were confirmed by nucleotide sequencing.

Enzyme assay
The activities of crude extracts and the purified recombinant BcGDH complex were determined using methods described in a previous study (Inose et al., 2003) with slight modifications. The enzyme sample was incubated at room temperature with 10 mM potassium phosphate buffer pH 7.0 containing 6 mM 5-methylphenazinium methylsulfate (phenazine methosulfate; PMS), 0.06 mM 2,6-dichlorophenolindophenol (DCIP) and various concentrations of glucose. The activity was determined by monitoring the decrease in the absorbance of DCIP at 600 nm and using the molar absorption coefficient of DCIP (16.3 mM cm À1 at pH 7.0) to calculate the enzyme activity. The molar absorption coefficient of DCIP was determined by measuring the absorbance of fixed concentrations of DCIP at 600 nm in 10 mM potassium phosphate buffer pH 7.0. One unit of enzyme activity is defined as the amount of enzyme that oxidizes 1 mmol glucose per minute.

Crystallization
Initial crystal screening was performed using the sittingdrop vapor-diffusion method with a Mosquito system (TTP Labtech, Hertfordshire, England). The protein concentration of wild-type BcGDH was 8.9 mg ml À1 in Milli-Q water.
After a few days, yellow crystals were observed in a reservoir solution containing 60% Tacsimate pH 7.0. Well diffracting crystals were obtained in a droplet containing a mixture of 1.5 ml protein solution (5.7 mg ml À1 in Milli-Q water) and 0.75 ml reservoir solution (59.9-60.2% Tacsimate pH 7.0) in a well containing 50 ml reservoir solution using the sitting-drop method at 293 K.
The methionines in BcGDH were replaced with selenomethionines (SeMet BcGDH) in order to determine the initial phases for the structure factors of wild-type BcGDH. However, crystals of SeMet BcGDH were not obtained under the same conditions as those of wild-type BcGDH. Since the prepared SeMet BcGDH contained a small amount of nonprocessed -subunit (18 kDa), which was observed on SDS-PAGE gels, the crystallization of SeMet BcGDH was attempted in the presence of proteases using Proti-Ace (Hampton Research, California, USA). Some crystals of SeMet BcGDH appeared in a droplet consisting of 1.0 ml protein solution (10.2 mg ml À1 in Milli-Q water), 0.2 ml of a 0.1 mg ml À1 subtilisin solution and 1.0 ml reservoir solution (60% Tacsimate pH 7.0) in a well containing 50 ml reservoir solution using the sitting-drop method at 293 K.

X-ray crystallography
Single crystals of both wild-type BcGDH and SeMet BcGDH were mounted in cryoloops and directly flashcooled in a stream of nitrogen gas at 100 K. X-ray diffraction data were collected using an ADSC Quantum 270 CCD detector system on the PF-AR NE3A beamline at the High Energy Accelerator Research Organization (KEK), Tsukuba, Japan. Diffraction data were processed using HKL-2000 (Otwinowski & Minor, 1997) and the CCP4 suite . Although X-ray diffraction data were collected from the crystal of BcGDH to 2.2 Å resolution, the data actually used for structure determination were truncated at 2.6 Å resolution owing to an extensively high R merge in the outermost shell. The unit-cell parameters of the crystal were large (a = b = 110.5, c = 524.9 Å ) and the diffraction pattern was characterized by strong anisotropy.
The initial phases of SeMet BcGDH were obtained using the single-wavelength anomalous dispersion (SAD) method with the AutoSol program (Terwilliger, 2004). Since the unitcell parameters of the crystal of wild-type BcGDH were isomorphous to those of the SeMet derivative, the phases were transferred directly to the former and the model was constructed using AutoBuild in the PHENIX system (Adams et al., 2010;Afonine et al., 2012).
The anomalous dispersion of Fe atoms was utilized in order to determine the number and the positions of Fe atoms in the iron-sulfur cluster. A SAD data set from a crystal of wild-type BcGDH was collected to 1.74086 Å resolution on beamline PF-AR NW12A at KEK using an ADSC Quantum 210r CCD detector system. Data-collection and refinement statistics for all data sets are listed in Table 1

Results
3.1. The overall structure of BcGDHca Fig. 1 and Supplementary Fig. S1 show the overall structure of the complex of theand -subunits of BcGDH and the topology of the protein, respectively. The -subunit consists of five -helices (colored red). The overall structure of the -subunit comprises 15 -helices (colored blue) and 17 -strands and adopts an FAD-binding fold. The additional domain contains a six-stranded antiparallel -sheet surrounded by six -helices and a protruding loop including two -helices (11 and 12) facing towards the -subunit. Two distinguishing long loop regions are located between 2 and 3 and between 4 and 8. The former loop (Ala39-Leu83) contains a unique -helix (2) that is located in close proximity to the -subunit at the center of the -subunit complex and above the iron-sulfur cluster shown in the red circle [ Fig. 1(a)], which is described later. The latter loop (Glu197-Asn229) surrounds the iron-sulfur cluster and contains a cysteine cluster (Cys212, Cys213, Cys218 and Cys222); three of the four cysteines are involved in forming the 3Fe-4S cluster.
The structure around FAD in the -subunit is shown in Fig. 2. FAD is surrounded by part of the long loop between Ala98 and Ser109, 1, 15, 2 and 6. Since continuous electron density was observed between FAD and His105 in the -subunit of BcGDH [ Fig. 2(b)], a covalent bond was deemed to form between FAD and His105 in the -subunit of BcGDH.
In the crystal, two pairs of BcGDH complexes were observed in the asymmetric unit ( Supplementary Fig. S2). As shown in the surface model presented in Figs. 1(a) and 1(b), BcGDH forms a heterodimer with a bent form, and the bent complex (-subunits) faces the back side of the other complex in the crystal ( Supplementary Fig. S2). Although the protein seems to form a heterotetramer, 0 0 , according to PISA analysis (http://www.ebi.ac.uk/ pdbe/pisa/) this assembly is unstable and forms two pairs ofand -subunits ( and 0 0 ). A gel-filtration chromatogram revealed the absence of an oligomeric form of BcGDH in solution (data not shown). Therefore, the observed 0 0 heterotetramer in the crystal structure is an artifact of crystallization. Since each molecule of and 0 0 is almost identical, with r.m.s. deviations for C atoms of 0.48 ( and 0 ) and 0.89 ( and 0 ), the structural description concentrates on unless otherwise specified.

Iron-sulfur cluster
The locations of Fe atoms were identified using the anomalous dispersion method and were observed in the expected positions of the iron-sulfur cluster. Three Fe-atom sites were identified in the iron-sulfur cluster of each molecule [Figs. 3(a)-3(d)]. The simulated-annealing OMIT maps of sulfur ions in the iron-sulfur cluster and the disulfide bond indicated that this iron-sulfur cluster is a 3Fe-4S cluster coordinated by Cys212, Cys218 and Cys222 of the -subunit. A neighboring cysteine, Cys213, in the -subunit forms a disulfide bond with Cys152 of the -subunit. The unique cysteine cluster of BcGDH contributes to the formation of the 3Fe-4S cluster for electron transfer and the disulfide bond for stabilization of the structure of the complex.

Interface between the cand a-subunits
The -subunit tightly binds the -subunit and forms a stable heterodimer. The hydrophobic loop region at the C-terminus  of the -subunit (colored green in the black circle in Fig. 4) contacts one of the distinguishing long loops, including 2 in the -subunit at the interface (colored in cyan). A disulfide bond is located at the interface between theand -subunits (Cys152 in the -subunt and Cys213 in the -subunit), as indicated by the red circle in Fig. 4, and is in close proximity to the hydrophobic cluster created by the C-terminal region in the -subunit indicated by the black circle (Leu146, Val147, Ile148, Pro153, Pro156, Gly157, Phe158, Trp159, Ala160 and Pro163) and the surrounding residues of the -subunit (Leu172, Pro173, Leu174, Phe176, Leu333, Trp334, Pro335, Gly336, Gly338, Pro339 and Met342). A 3Fe-4S cluster is located next to the disulfide bond, and the iron-sulfur cluster is located in the hydrophobic environment formed by the above hydrophobic cluster, the hydrophobic loop region (Met219, Pro223, Ile224, Ala226 and Met227, colored cyan) of the -subunit, including four cysteine residues (Cys212, Cys213, Cys218 and Cys222) involved in formation of the 3Fe-4S cluster and the disulfide bond, and two alanine residues (Ala107 and Ala108) located between 3Fe-4S and the isoalloxazine ring of FAD. Furthermore, 4 and 5 of the -subunit make contacts with protruding helices, including the 11 and 12 helices of the -subunit (in the area indicated by the blue circle). These hydrophobic contacts contribute to the tightly bound subunit interface.

Site-directed mutagenesis
The crystal structure of the BcGDH complex revealed the presence of an inter-subunit disulfide bond between the side chains of Cys213 in the -subunit and Cys152 in the -subunit. As reported previously (Shiota et al., 2016), substitution of Cys213 in the -subunit by serine [(Cys213Ser)] only exerts a limited effect on the kinetic parameters of the BcGDH complex at room temperature. Thus, the inter-subunit disulfide bond is not essential for enzyme activity. A mutant BcGDH complex with a Cys152Ser mutation in the -subunit [(Cys152Ser)] was constructed and characterized to further confirm this hypothesis. Table 2  Overall structure of BcGDH. (a) Overall structure of BcGDH shown as cartoon (left) and surface (right) models. In the cartoon model, the -subunit is shown as red -helices and green loop regions, and the -subunit is represented with blue -helices, yellow -strands and gray loop regions. In the surface model, theand -subunits are shown in red and blue, respectively, and are tightly bound to each other. The bound FAD is represented as an orange stick model. The cysteine cluster indicated with a red circle includes the iron-sulfur cluster shown as spheres. (b) A view of the model rotated 180 around the x axis. mutant BcGDH complex at room temperature. In the absence of the electron-transfer subunit (-subunit), the BcGDH complex shows relatively low dye-mediated glucose dehydrogenase activity at the conventionally used concentration of the primary electron accepter (PMS; Yamazaki et al., 1999). Therefore, the enzyme activity of the BcGDH complex was determined using 6 mM PMS, a concentration that is tenfold higher than the condition used for the BcGDH complex containing the electron-transfer subunit ( Supplementary Fig. S4). The V max of the (Cys152Ser) complex is only moderately lower than that of the wild-type complex. The K m of the (Cys152Ser) complex is comparable to the K m of the wild-type complex.
While the substitutions of these cysteine residues did not substantially affect the kinetic parameters of the BcGDH complex at room temperature, disulfide bonds often contribute to the stability of the tertiary and/or quaternary structure of proteins. We therefore studied the enzyme activity of the cysteine-substituted mutants at higher temperatures. The wildtype BcGDH complex showed maximum activity at 343 K owing to its high thermal stability, whereas both the (Cys152Ser) and (Cys213Ser) complexes showed maximum activity at approximately 303-313 K (Fig. 5). A substantial decrease in the optimal reaction temperatures of the (Cys152Ser) and (Cys213Ser) complexes suggested an important role of the inter-subunit disulfide bond in maintaining the thermal stability of the wild-type BcGDH complex.  Table S1). These enzymes are categorized as members of the FAD-containing glucosemethanol-choline oxidoreductase (GMC) family.

Discussion
Among the identified enzymes in the GMC family, pyranose 2-oxidase (P2Ox) shows the highest similarity to the BcGDH -subunit [Supplementary Figs. S5(b) and S5(d)]. The structure of cholesterol oxidase (ChOx) is also similar to the structure of the BcGDH -subunit [Supplementary Figs. S5(c) and S5(e)].
The residues responsible for the catalytic reaction have previously been identified in a variety of GMC oxidoreductases. A His/His pair in GOx, fungal FADGDH, pyranose dehydrogenase (PDH) and aryl-alcohol oxidase (AAOx), and a His/Asn pair in cholesterol oxidase (ChOx), cellobiose dehydrogenase (CDH), choline oxidase (COx) and P2Ox are expected to function as catalytic pairs based on the crystal structures, site-directed mutagenesis, pH-dependence studies or theoretical calculations. The catalytic pair in BcGDH is likely to be His/Asn, represented by His476 and Asn519, and corresponds to His689 and Asn732 in cellobiose dehydrogenase.
The structure of the active site of BcGDH was compared with Phanerochaete chrysosporium cellobiose dehydrogenase (PcCDH) bound to 6-hydroxy-FAD and the inhibitor cello

Figure 2
FAD-binding site and FAD covalently bound to His105 of the -subunit.
(a) FAD-binding site with the surrounding environment. A long loop of -subunit including 2 is colored cyan. FAD is shown as an orange stick model. Other colors are the same as in Fig. 1. ( (Yamashita et al., 2013). In contrast, proteins with a similar structure to the -subunit were not identified in the DALI search. To date, the -subunit has been considered to be a hitchhiker protein that promotes the secretion of the catalytic subunit to the periplasm. In the bacterial twin-arginine translocation (TAT) pathway, folded proteins are transported across the bacterial cytoplasmic membrane through the recognition of N-terminal signal peptides containing the twin-arginine motif. In some representative twin-arginine signal peptides, -helical regions were predicted using the PSIPRED secondary-structure Cysteine cluster of BcGDH. A unique cysteine cluster is formed by one cysteine from the -subunit and four cysteines from the -subunit. Cys152 and Cys213 form a disulfide bond between theand -subunits. Cys212, Cys218 and Cys222 are involved in the iron-sulfur cluster (3Fe-4S). The identified positions of the 3Fe-4S cluster and the disulfide bond between theand -subunits are shown. Simulated-annealing OMIT maps of F3S (the 3Fe-4S iron-sulfur cluster) and five cysteine residues (Cys152, Cys212, Cys213, Cys218 and Cys222) in wild-type BcGDH contoured at 5 are shown in light blue in (a)-(d). The 2F o À F c electron-density maps using the merged data sets of wild-type BcGDH collected at wavelengths of 1.0 and 1.74086 Å (after the initial refinement and before including F3S in refinement) are shown in blue and contoured at 5 in (b), (c) and (d). S and Fe atoms are shown in yellow and orange, respectively. The identified 3Fe-4S cluster consists of three Fe atoms and four S atoms represented as spheres.
prediction method (Palmer et al., 2005). The -subunit of BcGDH contains the twin-arginine motif in the N-terminal signal peptide and was thought to belong to the Tat protein family. Indeed, the structure of the -subunit of BcGDH contains five -helices.
The DALI search revealed some proteins with limited homology (Supplementary Table S2 and Fig. S7), including the N-terminal domain (NTD) of Salmonella typhimurium chemotaxis receptor methyltransferase (CheR) in complex with S-adenosyl-l-homocysteine (SAH), which had the highest Z-score (4.7). These domains of the enzymes were reported to be essential for catalytic activity, although they are located far from the catalytic domain and no role as a hitchhiker protein was reported.
The X-ray structure of BcGDH also revealed the first structure of a hitchhiker protein in complex with the target protein. According to the results of site-directed mutagenesis studies, the formation of disulfide bonds is required to stabilize the catalytic subunit. In other words, the disulfide bond may prevent denaturation of the 3Fe-4S cluster, thereby maintaining the stability of this enzyme even at temperatures greater than 323 K. Thus, the hitchhiker protein may protect the iron-sulfur cluster until the formation of a complex with the electron-transfer subunit after the catalytic subunit has been secreted into the periplasmic space. The electrontransfer subunit is folded in the periplasmic space and forms a quaternary structure with the catalytic subunit complexed with the hitchhiker protein. However, a mutant catalytic subunit or a mutant hitchhiker protein in the enzyme complex was expressed and functional. These results support the lack of a requirement for disulfide bonds in the functional expression and secretion of the complex into the periplasmic space. Indeed, alignments of the primary structures of the catalytic subunits [ Supplementary Fig. S8(a)] and hitchhiker proteins [ Supplementary Fig. S8(b)] of FAD-dependent dehydrogenase complexes reveal that the cysteine residues are not conserved, indicating that the formation of the disulfide bond is not necessary for a complex to form between catalytic subunits and hitchhiker proteins. Therefore, the hitchhiker protein of FAD-dependent dehydrogenase complexes may mainly interact with and recognize the catalytic subunit throughout these interfacial interactions.
Our previous report revealed the presence of a 3Fe-4S cluster in BcGDH (Shiota et al., 2016). The X-ray structure of the catalytic subunit clearly indicated the position of the 3Fe-4S cluster, which is located on the surface of the catalytic subunit. The distance between N5 of FAD and the 3Fe-4S cluster is about 12-13 Å , which is an adequate distance for electron transfer. These results support our hypothesis that the 3Fe-4S cluster functions in the intramolecular electron transfer from FAD and mediates intermolecular electron transfer from the 3Fe-4S cluster to the electron-transfer subunit. The 3Fe-4S cluster is responsible for electron transfer and interacts with the multi-heme c electron-transfer subunit.
As seen in the electron-density map, a covalent bond might form between C8M of FAD and His105 in the -subunit of BcGDH [ Fig. 2(b)]. The position of His105 is conserved in the GMC oxidoreducatase family, such as in P2Oxs (Bannwarth et al., 2004;Halada et al., 2003;Hallberg et al., 2004;Hassan et al., 2013;Spadiut et al., 2010;Tan et al., 2013), COxs (Quaye et al., 2008) and fumarate reductases B (Iverson et al., 1999(Iverson et al., , 2003   Tightly bound interface between theand -subunits. The highly hydrophobic C-terminal tail of the -subunit (colored green) including Cys152 of the -subunit covers the contact areas like a lid at the interface. The disulfide bond between subunits is indicated by a red circle. One of the long loops of the -subunit including 2 is colored cyan and contacts the C-terminal loop region of the -subunit. The other colors are the same as in Fig. 2. 4 and 5 of the -subunit also make contacts with protruding helices (11 and 12) of the -subunit in the area indicated with a blue circle.  The structures of (a) BcGDH, (b) W. succinogenes membrane-bound fumarate reductase (PDB entry 2bs2; Madej et al., 2006) and the flavoprotein shown in light green in (b) were superimposed onto the structure of BcGDH in (c). The colors of BcGDH are the same as in Fig. 2. The 3Fe-4S cluster is indicated by a red dotted circle. Soluble S. putrefaciens fumarate reductase was first superimposed onto the structure of the -subunit of BcGDH with 3.0 r.m.s.d., and the membrane-bound fumarate reductase from W. succinogenes was superimposed onto the structure of soluble fumarate reductase from S. putrefaciens (with 1.8 r.m.s.d.). Consequently, the structure of the membrane-bound fumarate reductase was directly compared with the structure of BcGDH by removing the structure of the soluble fumarate reductase. (b) The two-molecule form of fumarate reductase is shown in the crystal structure (Lancaster et al., 1999). A molecule of fumarate reductase is represented in light pink (right). In the other molecule, subunit A (flavoprotein containing FAD covalently bound to His43), subunit B (iron-sulfur protein including 2Fe-2S, 4Fe-4S and 3Fe-4S clusters) and subunit C (transmembrane-spanning protein including diheme cytochrome b molecules) are colored light green, light orange and dark green, respectively. The bound FAD molecules in the structures are shown as orange stick models. The bound FAD and heme b molecules are shown as orange lines and magenta stick models, respectively. The iron-sulfur clusters are shown as sphere models. The other molecule of fumarate reductase colored light pink was deleted in (c) to clarify the positions of FAD, iron-sulfur clusters and heme b molecules, which are related to electron transfer. bP is the proximal heme and bD is the distal heme. Lancaster et al., 2001;Madej et al., 2006), and this residue forms covalent bonds with FAD.
Next, we attempted to predict the position of the electrontransfer subunit by comparing of the structures of CDH and fumarate reductase, considering the previously elucidated intramolecular and intermolecular electron-transfer pathways (Shiota et al., 2016;Yamashita et al., 2018), as well as the nature of the electron-transfer subunit (Okuda-Shimazaki et al., 2018).
In Supplementary Fig. S9, the structure of the catalytic domain of CDH (closed state) was superimposed onto the structure of the -subunit of BcGDH [ Supplementary Fig.  S9(d)]. The 3Fe-4S cluster between theand -subunits of BcGDH (red dotted circle) is located on the left side, which is the opposite side to the heme b-type (magenta stick) electrontransfer domain of CDH. Considering that the primary electron acceptor of FAD is the 3Fe-4S cluster, and that in the next step intermolecular electron transfer occurs between the 3Fe-4S cluster and the electron-transfer subunit, the position of the electron-transfer subunit of BcGDH would be opposite to that of the heme b-type electron-transfer domain of CDH in the closed state.
Soluble flavocytochrome c fumarate reductase from Shewanella putrefaciens is a periplasmic tetraheme flavocytochrome c that consists of an N-terminal tetraheme cytochrome c domain and a catalytic region that contains the three C-terminal domains. The N-terminal domain containing the tetraheme moiety is connected by an -helical linker to the FAD-binding catalytic domain with noncovalently bound FAD. Membrane-bound diheme-containing quinol:fumarate reductase (QFR) from Wolinella succinogenes is composed of three subunits (A, B and C), in which subunit A contains the catalytic site of fumarate reduction and an FAD covalently bound to His43. Subunit B contains three iron-sulfur clusters and subunit C is a diheme cytochrome b (Lancaster et al., 2001;Madej et al., 2006).
Focusing on the homology between the soluble flavocytochrome c fumarate reductase from S. putrefaciens and the membrane-bound fumarate reductase from W. succinogenes, as well as the homology between soluble flavocytochrome c fumarate reductase from S. putrefaciens and the -subunit of BcGDH, the structure of membrane-bound fumarate reductase from W. succinogenes was superimposed onto the structure of the -subunit of BcGDH (Fig. 6). As shown in the membrane-bound fumarate reductase, the electron is transferred from FAD to the transmembrane protein containing heme bP and heme bD via 2Fe-2S, 4Fe-S and 3Fe-4S clusters. Interestingly, the positions of FAD and the first Fe-S cluster, 2Fe-2S, of fumarate reductase, with a distance of 12.3 Å , are comparable to those of FAD and the 3Fe-4S cluster of BcGDH, with a distance of around 12-13 Å . The distances between Fe-S clusters are 11.0 Å (2Fe-2S and 4Fe-4S) and 9.1 Å (4Fe-4S and 3Fe-4S), whereas the distances from the Fe-S clusters to the heme domains are 17.6 Å (3Fe-4S and cytochrome bP) and 15.6 Å (cytochrome bP and cytochrome bD). Although membrane-bound fumarate reductase contains three Fe-S clusters and two heme domains, the -subunit of BcGDH contains three heme c moieties in its electron-transfer subunit. The intact BcGDH complex including the membranebound -subunit containing three heme c moieties may form a similar overall structure to fumarate reductases for effective electron transfer. Interestingly, our previous study on the -subunit suggested that the electron from the Fe-S cluster is initially transferred to the third heme in the -subunit (the C-terminal heme domain of the -subunit), is then transferred to the second heme, is further transferred to the first heme (the N-terminal heme domain of the -subunit) and is finally transferred to an external artificial electron acceptor. However, when the intact BcGDH is immobilized on the electrode, the electron is transferred from the third heme to the second heme, and is then directly transferred to the electrode. The third heme of BcGDH may correspond to the position between 4Fe-4S and 3Fe-4S of fumarate reductase, the second heme corresponds to heme bP and the first heme corresponds to heme bD.
In conclusion, this study reports the first X-ray structure of a representative DET-type FAD-dependent dehydrogenase complex: the BcGDH catalytic subunit complexed with a hitchhiker protein. The structure of BcGDH revealed a conserved GMC oxidoreductase-type scaffold and a His/Asn catalytic pair, with a unique structure of the 3Fe-4S cluster, which serves as the electron acceptor of FAD and simultaneously serves as the electron donor for electron transfer of the multiheme c subunit. These findings will be essential for improving our understanding of intramolecular and intermolecular electron transfer by DET-type FAD-dependent dehydrogenase complexes, as well as for engineering DETtype enzymes for the development of future bioelectrochemical devices.

Related literature
The following references relate to PDB entries that are mentioned in the supporting information to this article: Batra et al. (2016), Djordjevic & Stock (1997, 1998