Structures of human dual oxidase 1 complex in low-calcium and high-calcium states

Dual oxidases (DUOXs) produce hydrogen peroxide by transferring electrons from intracellular NADPH to extracellular oxygen. They are involved in many crucial biological processes and human diseases, especially in thyroid diseases. DUOXs are protein complexes co-assembled from the catalytic DUOX subunits and the auxiliary DUOXA subunits and their activities are regulated by intracellular calcium concentrations. Here, we report the cryo-EM structures of human DUOX1-DUOXA1 complex in both high-calcium and low-calcium states. These structures reveal the DUOX1 complex is a symmetric 2:2 hetero-tetramer stabilized by extensive inter-subunit interactions. Substrate NADPH and cofactor FAD are sandwiched between transmembrane domain and the cytosolic dehydrogenase domain of DUOX. In the presence of calcium ions, intracellular EF-hand modules might enhance the catalytic activity of DUOX by stabilizing the dehydrogenase domain in a conformation that allows electron transfer. Dual oxidases (DUOXs), assembled from the catalytic DUOX and the auxiliary DUOXA subunits, produce hydrogen peroxide by transferring electrons from intracellular NADPH to extracellular oxygen in a calcium-activated manner. Here authors report the cryo-EM structures of human DUOX1-DUOXA1 complex in both high-calcium and low-calcium states.

R eactive oxygen species (ROS) are oxygen-containing chemical species that are highly reactive, such as hydrogen peroxide and superoxide anion 1 . They participate in many physiological processes and are implicated in several pathological conditions 1 . ROS can be generated by a class of dedicated enzymes called NADPH oxidase (NOX) in a highly regulated manner. These enzymes are multi-pass transmembrane proteins that catalyze the reduction of extracellular or luminal oxygen by intracellular NADPH to generate superoxide anion or hydrogen peroxide. NOX proteins are involved in many biological processes, including host defense, differentiation, development, cell growth and survival, cytoskeletal reorganization, and modification of the extracellular matrix 2 .
Comprising the human NOX protein family are NOX1-5 and DUOX1-2 (ref. 2 ). NOX2 protein catalyzes the production of superoxide anion during phagocytosis in neutrophils and is essential for host defense 3 . DUOX1-2 proteins are highly expressed in thyroid tissue and they catalyze the production of hydrogen peroxide, which is important for the biosynthesis of thyroid hormones 4 . The function of DUOX protein requires physical interactions with an auxiliary protein called dual oxidase maturation factor (DUOXA) 5 . DUOXA promotes the maturation and proper plasma membrane localization of DUOX 5 . DUOX protein is encoded by two homologous genes in human, namely DUOX1 and DUOX2. Similarly, DUOXA protein is encoded by DUOXA1 and DUOXA2. Loss-of-function mutations of DUOX2 or DUOXA2 in human cause congenital hypothyroidism 6 . Because of the important role of DUOX in thyroid tissue, they are also named thyroid oxidase 4 .
NOX family proteins share a common catalytic core, formed by a heme-coordinating transmembrane domain (TMD) and a cytosolic dehydrogenase (DH) domain 7 . The DH domain binds intracellular substrate NADPH and cofactor FAD, and shares sequence homology to the ferredoxin-NADP + reductase (FNR), which is composed of two subdomains 8 . In addition to the shared TMD-DH catalytic core of NOX, the functional DUOX protein has an additional large N-terminal extracellular peroxidase homology domain (PHD) and a long intracellular loop 0 containing two EF-hand domains, and it requires an auxiliary DUOXA protein for proper function. The activity of DUOX is regulated by intracellular calcium concentration 4 . Prior to our studies, the structures of NOX family members are only available in the form of isolated domains, including the DH domain (PDB ID: 5O0X) 9 and TMD (PDB ID: 5O0T) 9 of NOX5 from the algea Cylindrospermum stagnale (csNOX5) and a subdomain of human NOX2 DH domain (PDB ID: 3A1F). Despite the functional importance of DUOX and other NOX family members, their structures in the context of full-length functional protein complex are still unknown. Several open questions for DUOX remain elusive: How is the DH domain engaged with TMD to perform the catalytic redox reaction? How does DUOXA protein interact and co-assemble with DUOX? How is the activity of DUOX regulated by intracellular calcium? To answer these fundamental questions, we sought to characterize DUOX-DUOXA protein complex both structurally and functionally. Here, we present the cryo-EM structures of human DUOX1-DUOXA1 (hDUOX1-hDUOXA1) complex in both highcalcium and low-calcium states, providing insights into the structure and mechanism of calcium activation for DUOX.

Results
Structure determination. To express the hDUOX1 protein, we constructed the matured hDUOX1 protein  in frame with N-terminal GFP tag guided by a rat FSHβ signal peptide for efficient secretion 10 . The molecular weight of GFP-tagged hDUOX1 is 219 kDa. To monitor the formation of DUOX1-DUOXA1 complex, we fused the hDUOXA1 protein with a C-terminal MBP-mScarlet tag to increase its molecular weight to 106 kDa. Fluorescence size-exclusion chromatography (FSEC) showed the coexpression of hDUOXA1 effectively shifted the peak of hDUOX1 toward higher molecular weight, suggesting the formation of a stable hDUOX1-hDUOXA1 hetero-complex (Fig. S1a, b). The peak positions indicated hDUOX1 migrated as a monomer, while hDUOX1-hDUOXA1 complex migrated as a heterotetramer (Fig.  S1b). Moreover, we found the co-expression of DUOX1 and DUOXA1 resulted in cell membranes that showed robust calciumactivated, NADPH-dependent hydrogen peroxide production detected by the Amplex Red assay 11 (Fig. 1a-c). In the low-calcium condition, DUOX1-DUOXA1 complex showed low basal activity (Fig. 1b, c). Addition of calcium not only reduced the K m , but also increased the K cat of DUOX1 complex, leading to the overall enhancement of enzymatic activity (Fig. 1b,c).
We solubilized and reconstituted hDUOX1-hDUOXA1 complex into peptidisc using NSPr 12 . Purified peptidisc sample showed a homogenous peak on SEC (Fig. S1c) and two major protein bands on SDS-PAGE, both of which could be trimmed upon PNGase F treatment (Fig. S1d), suggesting both of hDUOX1 and hDUOXA1 were modified by N-linked glycosylation. UV-vis spectrum showed the peptidisc sample has characteristic Soret band with peak at 415 nm ( Fig. S1e), indicating proper Fe (III) heme incorporation. Moreover, the highly purified peptidisc sample recapitulated the calcium-activated NADPH-dependent hydrogen peroxide production observed on membrane (Fig. S1f), confirming that the calcium-dependent activation is a built-in mechanism of hDUOX1-hDUOXA1 protein complex. However, we found the maximum activity of purified peptidisc sample was lower than the activity measured using crude cell membrane ( Fig. 1b and Fig. S1f), suggesting either membrane bilayer or endogenous lipids might play a role on DUOX activity. We prepared cryo-EM grids using the peptidisc sample, either in the presence of 2.5 mM ethylene glycol tetraacetic acid (EGTA; low calcium) or 0.5 mM free calcium (high calcium). Both samples contained 0.1 mM FAD as cofactor and 0.5 mM NADPH as substrate.
Single particle cryo-EM analysis showed the purified protein was homogeneous and showed twofold symmetry (Figs. S2-S4). The overall resolution of cryo-EM maps in the low-calcium and high-calcium states reached 2.7 and 2.6 Å, respectively (Table S1). The extracellular domains and TMD showed better local resolution than the cytosolic domains, suggesting the higher mobility of the cytosolic domains (Figs. S2g and S4g). To further improve the map quality of cytosolic domains, we exploited symmetry expansion 13 and multibody refinement 14 by dividing one protomer into the large body (the extracellular domain and TMD) and the small body (the cytosolic domains; Figs. S2c and S4c). The final resolutions of cytosolic domain reached 3.4 and 3.2 Å for the low-calcium and high-calcium states, respectively (Figs. S2-4 and Table S1). The high map quality and available homology structures allowed us to build the order regions of the complex, which encompassed 88% of DUOX1 and 79% of DUOXA1 (Figs. S5-8 and Table S1). In the following text, we will focus on the high-calcium state structure unless noted otherwise, because of its higher resolution.
The architecture of hDUOX1-hDUOXA1 protein complex. hDUOX1 subunits and hDUOXA1 subunits co-assemble into a 2:2 heterotetrameric protein complex with molecular weight 457 kDa. The complex encompasses 140 Å × 105 Å × 160 Å 3D space and has an overall twofold rotational symmetry (Fig. 1d-f). Vertically, the complex can be divided into three layers: the extracellular layer, the transmembrane layer, and the cytosolic layer (Fig. 1d). In the extracellular layer, the two large N-terminal PHD domains of hDUOX1 pack against each other diagonally and are buttressed by the extracellular domain of DUOXA1 from beneath ( Fig. 1d-f). The transmembrane layer is formed by 24 transmembrane helices and harbors the heme-binding sites that provide the electron transfer pathway across the membrane (Fig. 1g)   surface of this cavity is highly hydrophobic (Fig. S3i) and there are several lipid molecules bound on this surface (Fig. S3i), suggesting this cavity is probably filled by phospholipids on the cell membrane. The cytosolic layer is comprised of the catalytic DH domain and regulatory domains for intracellular calcium sensing (Fig. 1f, i).
Structure of the catalytic hDUOX1 subunit. hDUOX1 is the catalytic subunit of the complex (Fig. 2). On the extracellular side of hDUOX1 resides the large N-terminal PHD domain which shares sequence homology with several peroxidases, such as peroxidase A from Dictyostelium discoideum (DdPoxA, PDB ID: 6ERC) 15 (Fig. S9a). Functional peroxidases utilize histidinecoordinated heme as the cofactor for catalysis. However, key residues for heme binding, such as the heme ligand histidine, are missing in the PHD of hDUOX1. Indeed, we did not observe any heme density in the structure of hDUOX1 PHD, suggesting PHD is probably not enzymatic functional in term of peroxidase activity. Close inspection of the map reveals two putative cation densities in PHD. One cation (cation binding site 1, CBS1) is coordinated by the side chains of D397 and T332, and the main chain carbonyl groups of V399, T332, and R395 (Fig. S9b). The second cation (CBS2) is coordinated by the side chain of D109, D174, S176, and T170, and the carbonyl groups of T170 and W172 (Fig. S9b). We observed strong densities in these two sites in both low-calcium and high-calcium conditions ( Fig. S9b), suggesting the bound cations might be sodium ions which were present in large quantities in our protein sample or calcium ions that bind very tightly. Both CBS1 and CBS2 are evolutionary conserved in DUOX (Fig. S5) and DdPoxA 15 (Fig. S9c), indicating their functional importance. Interestingly, we found both CBS1 mutant (D397A + T332A) and CBS2 mutant (D109A + D174A) of DUOX1 failed to co-assemble with DUOXA1 ( Fig.  S9d). Because CBS1 and CBS2 are away from the subunit interfaces in the DUOX1-DUOXA1 complex, we speculate these mutants probably affect the folding of PHD domain, suggesting the role of CBS1 and CBS2 in protein stability. The PHD packs on top of the TMD of DUOX1 through multiple noncovalent interactions (Fig. 2a). Moreover, a disulfide bond between C118 on PHD and C1165 on loop C of TMD further staples the bottom of PHD onto the top of TMD (Fig. 2a, d). In the TMD, hDUOX1 has an extra bent M0 helix at the periphery of M3 and M4, together with the canonical six TM helices of NOX protein family. M1-M6 of hDUOX1 form two heme-binding sites within the TMD. H1144 on M3 and H1238 on M5 coordinate the outer heme (Fig. 2b). H1130 on M3 and H1225 on M5 coordinate the inner heme (Fig. 2c). These four histidines are absolutely conserved in NOX family proteins (Fig. S6). We observed a spherical density surrounded by the invariant R1087 on M2, H1148 on M3, and outer heme-coordinating residue H1144 ( Fig. 2b and Fig. S3c). Previous studies showed mutations of the csNOX5 residues corresponding to R1087 and H1148 of hDUOX1 affected the reoxidation of dithionite-reduced TMD by oxygen and this site was proposed to be the oxygen substrate binding site, namely oxygen-reducing center (Fig. S10a, b) 9 . Our structure observations in hDUOX1 support the hypothesis.
Preceding the M1 helix of DUOX1 TMD, an amphipathic preM1 helix floats on the inner leaflet of plasma membrane (Fig. 1h). This helix was previously observed in csNOX5 (ref. 9 ) and is probably a shared feature of NOX family proteins. Between M0 and preM1 is a long cytosolic fragment loop 0. Cryo-EM maps reveal that the N-terminal of loop 0 is a domain rich of β sheets (Figs. S3g and S10c). Structural search using DALI server 16 identified the β sheets-rich domain is a crypto pleckstrin homology-like domain (PHLD) that shares little sequence homology, but high structural similarity to the PH domain proteins (Fig. S10c) 17 .
Following the PHLD, two EF-hand type calcium-binding domains (EF1 and EF2) form a compact helical module that is connected to the PHLD through αC (Fig. S6). Residues predicted to be responsible for calcium binding in EF1 and EF2 are evolutionary conserved in DUOX family proteins (Fig. S6). Although we did not observe the strong densities for small calcium ions due to poor local resolution (Figs. S2-3), the structure of EF-hand module closely resembles the small subunit of calcium-dependent protein phosphatase calcineurin in the calcium-bound state (PDB ID: 4IL1) 18 (Fig. S10d), suggesting both EF1 and EF2 are loaded with calcium in the high-calcium state. Based on the homology structure (4IL1), side chains of D828, D830, N832, and E839 and the main chain carbonyl group of Y834 chelate one calcium ion in EF1 (Fig. S10e) and side chains of D864, D866, N868, E875, and the main chain carbonyl group of L870 chelate another calcium ion in EF2 (Fig. S10f). It is reported that mutations of any of these calcium-binding sites abolished calcium activation 19 and E879K mutation in hDUOX2 (E875 in hDUOX1) leads to congenital hypothyroidism 20 , emphasizing their importance in calcium activation.
The C-terminal catalytic DH domain is connected to M6 of DUOX1 TMD via a short linker (Fig. 1h, i). The DH domain of hDUOX1 has a canonical DH fold and its structure is similar to csNOX5 (ref. 9 ; Fig.1i and Fig. S7). We observed strong densities for both FAD cofactor and NADPH substrate, and their binding sites were contributed from not only DH domain, but also TMD (Fig. 2h, i), as described later.
Inter-domain interactions in the high-calcium state. In the high-calcium state, individual domains of DUOX1 in the cytosolic layer are stabilized by multiple inter-domain interactions. The PHLD interacts with adjacent TMD and DH domains (Fig. 2e). The main chain carbonyl group of K653 on PHLD makes hydrogen bond with R1215 on loop D of TMD (Fig. 2e). Side chain of R674 of PHLD interacts with the main chain carbonyl group of E1348 and I1349 on α1 of the DH domain (Fig. 2f). The EF1-EF2 module in the high-calcium state shapes a crevice that embraces α4 and post α4 loop of the DH domain ( Fig. 2g and Fig. S10g). The interactions between the EF module and DH are mainly hydrophobic and involve F768, F772, F807, F819, F840, and F847 of the EF module, L1463, M1467, I1470, F1475, V1478, and F1484 of the DH domain (Fig. 2g). In addition, K814 of the EF module makes electrostatic interaction with E1281 on β2 of DH (Fig. 2g). The interactions between the EF module and the DH domain of hDUOX1 mimic the interactions between calcineurin subunit B and A in the calcium-bound state (PDB ID: 4IL1) 18 (Fig. S10g, h).
The linker between the EF module and preM1 helices binds in a groove on the surface of the DH domain (Fig. 1i). DH docks onto the bottom of TMD via polar interaction between R1270 on M6 and D1367 on β7, and between R1113 on loop B of TMD and N1550 of DH (Fig. 2h, i). It is reported that R1111Q mutation in hDUOX2 (R1113 in hDUOX1) was identified in congenital hypothyroidism patients 20 , highlighting the importance of this inter-domain interaction. Moreover, both the FAD cofactor and NADPH substrate bind at the interface between DH and TMD. R1214 and R1131 in TMD form electrostatic interaction with phosphate of FAD. D1128 makes hydrogen bonding with ribose of FAD (Fig. 2h). E1039 and N1040 in TMD make hydrogen bonding with adenosine ring of NADPH, and R1036 make cation-π interaction with both adenosine ring and electrostatic interaction with phosphate group of NADPH (Fig. 2i). Notably, R1495, R1424, and R1036 all participate in electrostatic interactions with the phosphate group of NADPH ribose, providing structural mechanism to distinguish NADPH from NADH (Fig. 2i). Through structural comparison, we found the NADPH-binding site in the csNOX5 structure was blocked by the artificially engineered C-terminal insertion, which was introduced into previous crystallization construct 9 (Fig. S10i). Moreover, the adenosine group of FAD has a 180°flip compared with structure of the isolated DH of csNOX5 (Fig. S10j). This is probably because D1128 on TMD stabilizes the ribose of FAD in such a conformation to make the connecting adenosine group of FAD in close proximity with inner heme for electron transfer (Fig. 2h). Taken together, the binding of FAD and NADPH at the interface The putative electron transfer pathway. The measured edge-toedge distances between NADPH and FAD, between FAD and inner heme, and between inner heme and outer heme are 8.2, 3.9, and 6.7 Å, respectively (Fig. 3a). Although, it is possible that there are additional protein residues on DUOX1 that rely electrons from NADPH to FAD, such as W378 between two hemes in csNOX5 (ref. 9 ), the distance between NADPH and FAD is larger than that in the canonical FNR protein, such as 3.2 Å in spinach FNR (sFNR, PDB ID: 1QFZ) 21 . Through structural comparison, we found the DH domain of DUOX1 shows a relaxed conformation, in which two subdomains are loosely packed, while both the DH of csNOX5 and sFNR show a tense conformation and their two subdomains are tightly packed against each other to bring FAD and NADPH into close proximity for electron transfer (Fig. 3b-e). Therefore, the electron transfer efficiency in the current structure of DUOX1 is not optimal. Because the DUOX1   complex on cell membrane exhibited higher activity ( Fig. 1b and Fig. S1f), it is possible that lipids on cell membrane or the bilayer environment could somehow affect the structure of DUOX1 to enhance its electron transfer efficiency. At the terminus of electron transfer chain near the extracellular side, the initial product of oxygen-reducing reaction is superoxide anion. We probed the possible pathways for oxygen entrance and for superoxide anion exit with CAVER 22 , using the oxygenreducing center as the starting point. We located four possible tunnels: tunnel A is formed by M1, M2, M5, and M6 and is capped by loop E on top (Fig. 3f); tunnel B is surrounded by M2, loop A, loop C, and loop E (Fig. 3g); tunnel C is embraced by M3, M4, and loop C (Fig. 3h); and tunnel D is enclosed by M3, M4, loop C, and M0 (Fig. 3i). The bottleneck radii of these tunnels arẽ 1 Å (Fig. 3j), which may allow the permeation of small oxygen substrate under dynamic motion of DUOX1 protein. Further analysis showed tunnels B-D are all surrounded by hydrophobic residues (Fig. 3g-i), which are unfavorable for superoxide anion permeation. In contrast, tunnel A is gated by hydrophilic R1087 on M2, R1062 on M1, and R1248 and Q1245 on loop E (Fig. 3f). We speculate the highly positively charged constriction of tunnel A would strongly attract the negatively charged superoxide anions, and this might be essential for the dismutation reaction between two superoxide anions to generate uncharged hydrogen peroxide for diffusion. Therefore, manipulations that may alter the constrictions of tunnels A-D would affect superoxide anion intermediate leakage. Indeed, it is reported that mutations on DUOX1 loop A or on DUOXA1 N-terminus peptide (NTP) which interacts with and stabilizes loop A would change the ratio of superoxide anion and hydrogen peroxide produced, probably by affecting the leakage of superoxide anions through these tunnels [23][24][25] .
Structure of hDUOXA1 and mechanism of complex assembly. DUOXA protein is an essential auxiliary subunit for DUOX enzyme 5 (Fig. S11a) and it has an extracellular N-terminus that is important for hydrogen peroxide generation 24,25 . We observed the NTP of hDUOXA1 extends and packs onto the PHD-TMD junction of the distal hDUOX1 subunit (Fig. 4a-e). Side chains of F8, F10, and Y11 of NTP insert into the hydrophobic groove formed by loop C, loop A, and PHD of hDUOX1 (Fig. 4c). In addition, K15 of DUOXA1 NTP makes electrostatic interactions with D1077 of DUOX1 (Fig. 4d). This agrees with previous data showing DUOXA1 NTP interacts with DUOX1 loop A (ref. 23 ). hDUOXA1 has five transmembrane helices. Lower part of TM1 interacts with preM1 and M1 of hDUOX1 (Fig. S11a). The remaining four helices and associated extracellular loops share structural similarity with claudin superfamily members, such as claudin-9 (PDB ID: 6OV2) 26 (Fig. S11). The extracellular loop between TM2 and TM3 folds into a compact claudin-like domain (CLD) composed of four β strands and two α helices (Figs. S8 and S11). CLD forms extensive interactions with both distal and proximal DUOX1 subunits (Fig. 4a, b), emphasizing its important role in the complex assembly. This agrees with previous studies showing that splicing variants at TM2-TM3 loop have distinct behavior in supporting the activity of DUOX1 (ref. 24 ). Moreover, we found an ordered N-linked glycosylation decoration on N109 of hDUOXA1 and its branched sugar moieties make extensive polar interactions with both DUOXA1 and DUOX1 subunits    Fig. 4a, b). The PHD of two DUOX1 subunits also interact with each other (Fig. 4a, b). Close to the dyad axis, R50 and R507 on one PHD make polar interactions with E41 and F313 on the opposite PHD (Fig. 4e). We further analyzed the effects of interface mutations on the tetramer assembly, and found mutations of R50E, R507E, and R507A all severely affect tetrameric peak formation on FSEC (Fig. 4f). These structural information and biochemical data revealed the detailed inter-subunit interactions that dictate the heterotetramer assembly.
Conformational change of DUOX1 complex upon calcium activation. The consensus map in the low-calcium state showed the cytosolic layer had poor local resolution, which was improved by multibody refinement 14 (Fig. S4). Further molecular flexibility analysis 14 showed the cytosolic domains (small body) in the lowcalcium state were sampling a broad range of orientations relative to the TMD, evidenced by the plateau-shaped distribution on the histogram of the major eigenvector (Fig. S4f). This is in great contrast to the normal distribution in the high-calcium state (Fig.  S2f), suggesting the cytosolic layer in the low-calcium state is more flexible. We compared the structures in the low-calcium state and high-calcium state, and found structural changes in the extracellular layer and transmembrane layer are small (Fig. 5a). However, there are large conformational changes of the regulatory PHLD and EF-hand module in the cytosolic layer (Fig. 5a-c and Movie S1). In the absence of calcium, the EF module switches from an extended shape into a more contracted shape (Fig. 5d, e), which reconfigures the interface between the EF module and α4 of the DH domain, resulting in a loosely packed structure (Fig. 5f). In the low-calcium state, EF2 moves away from the DH domain. The Cα atom of A894 on αJ of EF2 has 40 Å displacement (Fig. 5b). PHLD rotates away from the TMD and DH domains, and αA of PHLD has 17.2°outward rotation (Fig. 5c). As a result, several inter-domain interactions observed in the high-calcium state were disrupted and therefore the docking of DH domain onto TMD is weakened by these structural changes, leading to a higher mobility of DH domain (Fig. S4g). We propose the increased mobility of DH domain negatively correlates with the electron transfer efficiency and thus the catalytic activity of DUOX. In addition, because TMD also contributes to FAD and NADPH binding, the increased mobility of the DH domain would result in the reduced affinity of NADPH as well. This is in agreement with the markedly reduced K cat and moderately increased K m in the low-calcium state, as we observed (Fig. 1c).
During the preparation of this manuscript, another group reported the structures of mouse DUOX1-DUOXA1 complex 27 . Interestingly, they found mouse DUOX1 complex exists in both heterodimeric and heterotetrameric form, and they proposed the activity of DUOX1 complex is regulated by dimer-tetramer assembly 27 . This is in contrast to our observation that majority of hDUOX1 complex is in a homogenous tetrameric form (Fig. S1c). Whether this difference is due to different protein preparation procedure or different species (mouse vs human) remain elusive. Moreover, the intracellular PHLD and EF domains were not resolved in mouse DUOX1 complex structure because of insufficient map quality 27 . The overall structures of resolved parts between mouse and hDUOX1 complex are similar with root mean square deviation of 1.521 and 0.908 Å for DUOX1 and DUOXA1 subunit, respectively (Fig. S11c). However, detailed structural comparison revealed several differences especially in the atomic models of FAD (Fig. S11d) and NADPH (Fig. S11e), probably due to the poor local map quality of mouse DUOX1 complex (EMD-21964).

Discussion
In this study, we provided the structures of hDUOX1-hDUOXA1 as a peptidisc-stabilized heterotetrameric protein complex in both high-calcium and low-calcium states. The structure of hDUOX1 complex in the high-calcium state reveals multiple inter-domain interactions that orientate DH and TMD for electron transfer, and thus redox reaction. Removal of calcium ions results in the  Fig. 5 Conformational change of hDUOX1 complex during calcium activation. a Structural comparison of hDUOX-hDUOXA1 complex between the highcalcium state (colored) and the low-calcium states (gray). Protein is shown as cartoon. Regions with large conformational changes are boxed by dashed lines. b Close-up view of the conformational changes of EF-hand module. Cα atom of A894 on αJ helix is used as marker to measure the movement of EF2. c Close-up view of the conformational change of PHLD. The angle between αA helices in the high-calcium and low-calcium states was measured. d, e Conformational differences of EF-hand module between the high-calcium state and the low-calcium state. f Reconfiguration of the interface between EF-hand module and α4 helix of DH domain. Arrows denote movements from high-calcium state into the low-calcium state.
reconfiguration of cytosolic inter-domain interactions which in turn mobilizes the DH domain and lowers the electron transfer efficiency (Fig. 6). These structures provide mechanistic insights into the structure and mechanism of DUOX and other NOX enzymes.

Methods
Cell culture. HEK293F suspension cells (Thermo Fisher Scientific) were cultured in Freestyle 293 medium (Thermo Fisher Scientific) supplemented with 1% FBS at 37°C with 6% CO 2 and 70% humidity. Sf9 insect cells (Thermo Fisher Scientific) were cultured in SIM SF (Sino Biological) at 27°C. The cell lines were routinely checked to be negative for mycoplasma contamination.
Protein expression and purification. We constructed a modified BacMam vector 28,29 with N-terminal GFP tag guided by rat FSHβ signal peptide 10 and cloned hDUOX1 cDNA into this vector (please see Supplementary Table 2 for primer list). The hDUOX1 cDNA, we obtained from Prof. Han, has the same protein sequence as NP_059130.2 except the L1178F mutation, a SNP previously observed in AAI14939.1. The cDNAs of hDUOXA1 were cloned into a non-tagged BacMam vector or a modified BacMam vector with C-terminal MBP-mScarlet tag (please see Supplementary Table 2 for primer list) 28,29 . The hDUOX1 mutants were generated by Quick Change methods (please see Supplementary Table 2 for primer list). The two expression cassette were further merged into one bicistronic vector by the LINK sequence on the modified vector 28,30 . The baculoviruses were produced using the Bac-to-Bac system and amplified in Sf9 cells. For protein expression, HEK293F cells cultured in Freestyle 293 medium at density of 2.8 × 10 6 ml −1 were infected with 15% volume of P2 virus. A total of 10 mM sodium butyrate was added to the culture 12 h post infection and transferred to a 30°C incubator for another 36 h before harvesting. Cells were collected by centrifugation at 3999 × g (JLA 8.1000, Beckman) for 10 min, and washed with 20 mM Tris (pH 8.0 at 4°C), 150 mM NaCl, 2 mM EGTA, 2 μg ml −1 aprotinin, 2 μg ml −1 pepstatin, 2 μg ml −1 leupeptin, flash-frozen, and storage at −80°C.
Enzymatic assay. The membrane fractions of DUOX1 for enzymatic assay were prepared as previously reported with minor modification 31 . Briefly, cells were washed with buffer A. After centrifuging at 3999 × g for 10 min at 4°C, the cell pellets were broken using a needle for 12 times in 1 ml of 20 mM Tris pH 8.0 at 4°C containing 0.1 mM dithiothreitol, 10 mM EGTA (pH 8.0), and the mixture of protease inhibitors. The pellet was removed by centrifuging at 3615 × g for 15 min, and the supernatant was collected and then centrifuged at 228,600 × g (TLA100.3, Beckman) for 1 h. The membrane pellet was resuspended in HBS containing 1 mM EGTA. Meanwhile, 10 μl membrane was solubilized by 100 μl TBS + 1% digitonin with the mixture of protease inhibitors for 1 h at 4°C for FSEC. The protein concentrations in the membrane were estimated by comparing their GFP fluorescence signal to that of a purified GFP-tagged DUOX1 complex.
The H 2 O 2 -generating activity of DUOX1 complex was determined using the amplex red assay 32 . The concentrations of H 2 O 2 solution were determined by measuring UV-vis absorbance at 240 nm with spectrophotometer (Pultton) and calculated using molar extinction coefficient of 43.6 M −1 cm −1 . The concentration of H 2 O 2 solution was further validated by reacting with Amplex red to generate resorufin which has ε 571 = 69,000 M −1 cm −1 (ref. 32 ). Then the H 2 O 2 solution with known concentration was used to calibrate the resorufin fluorescence curve (excitation, 530 nm; emission, 590 nm) measured using a Microplate Reader (BioTek Synergy HT) at 37°C.
The H 2 O 2 -generating reaction of the membrane fraction containing DUOX1 complex was performed at 37°C in 0.15 ml of HBS with 1 mM EGTA, 10 μM FAD, 100 μM NADPH, 50 μM amplex red, 0.067 mg ml −1 horseradish peroxidase, and 0.0576 mg ml −1 SOD. Ca 2+ concentrations were determined using fluorescent indicators fura-2 or fluo3-FF. The K cat and K m values of the membrane fraction containing DUOX1 complex were determined at 37°C with different concentrations of NADPH in the presence or absence of 1.4 mM CaCl 2 . The H 2 O 2generating reaction of the purified DUOX1 complex in peptidisc was performed at 27°C in 0.15 ml of HBS + 1 mM EGTA, 10 μM FAD, 100 μM NADPH, 50 μM amplex red, 0.067 mg ml −1 horseradish peroxidase, and 0.0576 mg ml −1 SOD in the presence or absence of 1.1 mM CaCl 2 . Progress of the reactions was monitored continuously by following the increase of the resorufin fluorescence, and the initial reaction rates were obtained by fitting the curve with linear equation. The activity of DUOX1 complex was determined by subtracting the background of the corresponding buffer without enzyme. The data were processed with Microsoft Excel-2013, SigmaPlot-12.0, and GraphPad Prism 6.
Cryo-EM sample preparation and data acquisition. The peptidisc sample was supplemented with 2.5 mM EGTA (low calcium) or 0.5 mM free calcium (high calcium) for cryo-EM analysis, respectively. Both samples contain 100 μM FAD as the cofactor and 500 μM NADPH as the substrate. To overcome the preferred orientation problem, 0.5 mM non-solubilizing detergent fluorinated octyl-maltoside was added to the sample before cyro-EM sample preparation. Aliquots of 1.5 μL protein sample were placed on graphene oxide-coated grids, as previously reported 33 . Grids were blotted for 3 s at 100% humidity and flash-frozen in liquid ethane cooled by liquid nitrogen using Vitrobot Mark I (FEI). Grids were then transferred to a Titan Krios (FEI) electron microscope that was equipped with a Gatan GIF Quantum energy filter and operated at 300 kV accelerating voltage. Image stacks were recorded on a Gatan K2 Summit direct detector in super-resolution counting mode using Serial EM at a nominal magnification of 130,000× (calibrated pixel size of 1.045 Å pixel −1 ), with a defocus ranging from −1.5 to −2.0 μm. Each stack of 32 frames was exposed for 7.12 s, with a total dose~50 e − Å −2 and a dose rate of 8 e − pixel −1 s −1 on detector.
Image processing. The image processing workflow is illustrated in Figs. S2 and S4. A total of 7076 super-resolution movie stacks of the high-calcium state sample and  Fig. 6 Activation mechanism of DUOX1 complex by calcium. Two DUOX1 and one DUOXA1 subunit are shown as cartoon, and colored the same as Fig. 1d. Calcium ions are presented as green spheres. Electron transfer pathways are indicted with gray arrows. NATURE COMMUNICATIONS | https://doi.org/10.1038/s41467-020-20466-9 ARTICLE NATURE COMMUNICATIONS | (2021) 12:155 | https://doi.org/10.1038/s41467-020-20466-9 | www.nature.com/naturecommunications 2076 stacks of the low-calcium state sample were collected using Serial EM, and motion-corrected, dose weighted, and twofold binned to a pixel size of 1.045 Å using MotionCor2 (ref. 34 ). Contrast transfer function (CTF) parameters were estimated with Gctf 35 . Micrographs with ice or ethane contamination, and empty carbon were removed manually. Autopicking were performed using Gautomatch (kindly provided by Kai Zhang). All subsequent classification and reconstruction was performed in Relion 3.1 (ref. 36 ) unless otherwise stated. Reference-free 2D classification was performed to remove contaminants. Initial model was generated using cryoSPARC 37 . Particles were subjected to multi-reference 3D classification 38,39 and random-phase 3D classification 38,39 . Phase-randomized models were generated from the model obtained from previous refinement using randomize software (from the lab of Nikolaus Grigorieff). Further CTF refinement was then performed with Relion 3.1 using C2 symmetry. The particles were then re-extracted, re-centered, and re-boxed from 256 pixels to 320 pixels for consensus refinement in Relion 3.1 (ref. 36 ) and cryoSPARC 37 . To improve the density of cytosolic layer, particles were symmetry expended 13 for multibody refinement 14 . One soft mask (the large body) that covers the extracellular domain together with TMD of one protomer was generated from the consensus map, using UCSF Chimera and Relion 3.1 (ref. 36 ). The other soft mask (the small body) covers the cytosolic domains of the same protomer. 3D multibody refinements 14 were performed using the two soft masks and the parameters determined from previous consensus refinement. The motions of the bodies were analyzed by relion_-flex_analyse in Relion 3.1 (ref. 36 ). The two half-maps of each body generated by 3D multibody refinement were subjected to post-processing in Relion 3.1 (ref. 36 ). The masked and sharpened maps of each body were aligned to the consensus map using UCSF Chimera 40 , zoned to isolated nonoverlapping regions and summed using Relion 3.1 (ref. 36 ) to generate the composite maps for visualization and model building. All of the resolution estimations were based on a Fourier shell correlation of 0.143 cutoff after correction of the masking effect. B-factor used for map sharpening was automatically determined by the post-processing procedure in Relion 3.1 (ref. 36 ). The local resolution was estimated with Relion 3.1 (ref. 36 ).
Model building. The composite maps derived from multibody refinement were used for model building. The structures of PHD, TMD, EF1-2, and DH domains of hDUOX1 were generated using phyre2 server 41 based on PDB ID: 6ERC, 5O0T, 4IL1, and 5O0X, and manually docked into the cryo-EM maps using Chimera 40 . Initial models of PHLD were generated by Rosetta Web Server using ab initio mode 42 , manually selected according to the distances calculated by RaptorX Contact Prediction server 43 , and validated by the fitting between model and cryo-EM densities, especially the location of bulky aromatic residues. The partial model of hDUOXA1 were generated using EM builder 44 . The initial models were iteratively built using Coot 45 and refined using Phenix in real space 46 . Figures were prepared using UCSF chimera 40 , Chimera X 47 , and Pymol.