X-ray structure and enzymatic study of a bacterial NADPH oxidase highlight the activation mechanism of eukaryotic NOX
eLife assessment
In this manuscript, the authors investigate the properties of prokaryotic NADPH oxidases (NOX) and discuss the implications for NOX regulation and function. The structure of the S. pneumoniae Nox protein is an important step forward in our understanding of procaryotic NOX enzymes and the characterization and interpretation are convincing. The results will be of interest to structural biologists as well as biochemists focusing on enzymatic functions.
https://doi.org/10.7554/eLife.93759.3.sa0Important: Findings that have theoretical or practical implications beyond a single subfield
- Landmark
- Fundamental
- Important
- Valuable
- Useful
Convincing: Appropriate and validated methodology in line with current state-of-the-art
- Exceptional
- Compelling
- Convincing
- Solid
- Incomplete
- Inadequate
During the peer-review process the editor and reviewers write an eLife Assessment that summarises the significance of the findings reported in the article (on a scale ranging from landmark to useful) and the strength of the evidence (on a scale ranging from exceptional to inadequate). Learn more about eLife Assessments
Abstract
NADPH oxidases (NOX) are transmembrane proteins, widely spread in eukaryotes and prokaryotes, that produce reactive oxygen species (ROS). Eukaryotes use the ROS products for innate immune defense and signaling in critical (patho)physiological processes. Despite the recent structures of human NOX isoforms, the activation of electron transfer remains incompletely understood. SpNOX, a homolog from Streptococcus pneumoniae, can serves as a robust model for exploring electron transfers in the NOX family thanks to its constitutive activity. Crystal structures of SpNOX full-length and dehydrogenase (DH) domain constructs are revealed here. The isolated DH domain acts as a flavin reductase, and both constructs use either NADPH or NADH as substrate. Our findings suggest that hydride transfer from NAD(P)H to FAD is the rate-limiting step in electron transfer. We identify significance of F397 in nicotinamide access to flavin isoalloxazine and confirm flavin binding contributions from both DH and Transmembrane (TM) domains. Comparison with related enzymes suggests that distal access to heme may influence the final electron acceptor, while the relative position of DH and TM does not necessarily correlate with activity, contrary to previous suggestions. It rather suggests requirement of an internal rearrangement, within the DH domain, to switch from a resting to an active state. Thus, SpNOX appears to be a good model of active NOX2, which allows us to propose an explanation for NOX2’s requirement for activation.
Introduction
NADPH oxidases (NOX) are membranous enzymes that ‘professionally’ produce the reactive oxygen species (ROS) superoxide or H2O2 (Vermot et al., 2021); the ‘professional’ production in NOX contrasts with many other sources of cellular ROS which result from metabolic leaking or byproducts. The first-characterized human enzyme of the NOX family, now called NOX2, has been known since 1978 (Segal and Jones, 1978), and provided a paradigm of its ROS product as a cytotoxic agent in defense mechanisms. However, from 1999 to 2001, six other isoforms of NOX were discovered: NOX1/3/4/5 and DUOX1/2 (Bánfi et al., 2001; De Deken et al., 2000; Dupuy et al., 1999; Geiszt et al., 2000; Kikuchi et al., 2000; Suh et al., 1999). These other isoforms of NOX function in many tissues and cell types and the ROS they produce play roles beyond cytotoxic defense, as biochemical and signaling agents. In addition to the identification of NOX2 deficiency as the cause of chronic granulomatous disease (Segal and Jones, 1978; Segal and Jones, 1979), dysregulation of other NOX family members has since been highlighted in pathophysiologies like hyperthyroidism (De Deken et al., 2000; Dupuy et al., 1999) and inflammatory bowel diseases along with well-documented participation in normal physiological process like cardiovascular tone regulation (Rajagopalan et al., 1996), sense of balance (Bánfi et al., 2004), and fertility (Musset et al., 2012) among others. Interestingly, these integrated processes directly link to multicellularity leading to the belief that true NOXes are found only in multicellular eukaryotes (Lalucque and Silar, 2003).
The various isoforms display significant diversity in assembly with accessory membrane proteins (p22phox for NOX1-4, DUOXA1/2 for the corresponding DUOX1/2) or not (NOX5); in phosphorylation-promoted association of cytosolic factors required for activation (NOX 1–3) or not (NOX4-5, DUOX1/2); in the presence of EF hands conferring Ca2+ dependence (NOX5, DUOX1/2) or not (NOX1-4); and in the species of ROS produced, superoxide for NOX1-3,5 and H2O2 for NOX4 and DUOXs enzymes (Vermot et al., 2021).
Despite their variety, all NOX family members share a common catalytic core comprised of two conserved domains. The transmembrane (TM) domain belongs to the Ferric Reductase Domain (FRD) superfamily, which consists of six TM helices containing conserved histidines in TM 3 and 5 that chelate up to two heme b cofactors (Zhang et al., 2013). The cytosolic dehydrogenase (DH) domain belongs to the Ferredoxin NADP +Reductase (FNR) family which has motifs to bind a flavin cofactor and a pyridine dinucleotide NAD(P)H (Taylor et al., 1993). Thus, the NOX core catalytic domain contains the necessary machinery to move electrons from the electron donor – NADPH – through the flavin and heme cofactors to the ultimate electron acceptor, molecular oxygen, to make the ROS product.
The fusion of the TM (N-terminal) and DH (C-terminal) domains constitutes the conserved architecture for the catalytic subunit of the NOX family members. The widespread appearance of separate FRD and DH domains in prokaryotes motivated a search for prokaryotic NOX homologs. First, separate E. coli FRD and FNR representatives MsrQ and Fre were identified and shown to interact, producing a NOX-like flow of electrons from NADPH to flavin and through the hemes to an electron acceptor (Juillan-Binard et al., 2017). Another bioinformatic search identified a large number of prokaryotic sequences with fused FRD and FNR domains (Hajjar et al., 2017). We demonstrated that one such protein from Streptococcus pneumoniae, called SpNOX, is a bona fide NOX that produces ROS and shares the spectral signature and biochemical hallmarks (resistance to cyanide, inhibition by DPI) of eukaryotic NOX (Hajjar et al., 2017), although its biological function in S. pneumoniae has not yet been determined. These results demonstrate definitively that NOXes also exist in prokaryotes. SpNOX purifies as a monomer and is constitutively active (Hajjar et al., 2017), so it represents the minimal functional core shared by all NOX family members. Furthermore, unlike eukaryotic NOX family members, it overexpresses well in E. coli and retains robust activity in detergent solution. These qualities make SpNox well suited to structural determination by several methods.
Published crystal structures of separate TM and DH domains of a NOX5 homolog CsNOX from the cyanobacterium Cylindrospermum stagnale (Magnani et al., 2017) provided novel structural information about the NOX family. A low-resolution structure of full-length SpNOX (Vermot et al., 2020) provided the first look at the fused TM and DH domains in NOX proteins, while cryo-EM structures of mammalian DUOX1 (Sun, 2020; Wu et al., 2021) and NOX2 (Liu et al., 2022; Noreng et al., 2022) furnished important new insights. Still, much remains to be discovered regarding electron transfer, the essential catalytic function of the NOX family. SpNOX thus remains an attractive prototype through which to examine structure/function relationships common to NOX family members.
Here, we obtained X-ray structures of WT and F397W mutant of DH-only, and of a F397W mutant of full-length SpNOX. We combined enzymological and structural characterization to elucidate insights into substrate and cofactor recognition, as well as kinetic and mechanistic aspects of electron transfer for the NOX enzyme family. Comparative analysis of this constitutively active prokaryotic NOX with the previously described human NOX2 and DUOX1 structures provides a molecular explanation NOX2’s requirement for an activation step.
Results and discussion
Production of SpNOX and SpNOXDH domain
SpNOX shares only about 25% identity with the human NOX homologs; nonetheless the NOX family members share several well-conserved motifs identified in eukaryotic NOX by Massari et al., 2022. Figure 1A shows the logos corresponding to the eukaryotic NOX motifs in Massari et al., 2022, but constructed from our alignment containing about half prokaryotic sequences. Other than the motif identified as NADPH3, for which our joint prokaryotic-eukaryotic NOX alignment lacked consensus, SpNOX (and other prokaryotic NOX) recognizably share all the motifs identified by Massari et al., 2022.
The C-terminal motif identified by our alignment includes several residues past the C-terminal end of the motif identified from eukaryotes, indicating the high frequency of additional residues in prokaryotic NOX. A specific acid/base pair (E283R284 human NOX2 numbering) highly conserved in eukaryotic NOX appears directly after the predicted end of the 6th TM helix (Sharpe et al., 2010), at the inter-domain hinge (Figure 1B); SpNOX has the similar pair Q180K181 at this position which prompted us to define K181 as the start of SpNOXDH. A typical SDS-PAGE of purified SpNOX and SpNOXDH is shown (Figure 1—figure supplement 1).
Optimization of the SpNOXDH thermostability and activity
SpNOXDH showed a tendency to aggregate, so we investigated buffer conditions to optimize its stability, solubility and activity (Gekko and Timasheff, 1981; Martins de Oliveira et al., 2018). SpNOXDH thermostability and activity were assessed through wide ranges of salt concentrations (0.05–1 M) and glycerol percentage (0–20%). The temperatures of thermal unfolding were determined by nanoDSF and demonstrated a constant improvement of Tm with NaCl concentration and glycerol percentage increases (Figure 2A and B). SpNOXDH activity was assessed by reduction of cytochrome c. In reactions without NaCl and glycerol, the DH domain efficiently supported cytochrome c reduction (Figure 2C and D). In contrast to thermal stability, however, SpNOXDH exhibits a dramatic loss of activity with increasing NaCl and glycerol. We generated a 3D surface plot of Tm as a function of NaCl and glycerol with specific activity as color scale (Figure 2E) which allowed us to optimize buffer conditions jointly. The region highlighted by (○) in Figure 2E indicates that at 300 mM NaCl and 5% glycerol SpNOXDH exhibits both high activity (5 mol reduced cyt. c/s/mol.SpNOXDH) and a good thermostability (melting temperature of 60 °C).
We then optimized the pH in order to conjointly maximize the thermostability and the activity of SpNOXDH. Direct comparison of the kinetic parameters of SpNOXDH and SpNOX necessitated identifying a pH that would work well for both constructs. The dependences of SpNOXDH activity (Figure 2—figure supplement 1A) and stability (Figure 2—figure supplement 1B), evaluated along a range of pH (respectively 5–7 and 4–9.5) in the presence of 300 mM NaCl 5% glycerol, exhibit a progressive decrease with increasing alkalinity. Like SpNOXDH, the thermostability of SpNOX follows a similar and progressive decrease of Tm with alkalinity (Figure 2—figure supplement 1B), but in contrast the activity of SpNOX shows an optimum at higher alkalinity (Figure 2—figure supplement 1A). The two constructs show similar activity at pH 6.5, so we chose this as the pH. Thus, for all experiments in which we compare activity of the full-length SpNOX to the truncated SpNOXDH, we chose a standard buffer of 50 mM bis TRIS-propane pH 6.5, 300 mM NaCl, 5% glycerol.
Functional characterization of SpNOXDH
The DH domain of NOXes, including SpNOX, as a homolog of the FNR superfamily, should catalyze the transfer of two electrons from the cytosolic NADPH to the FAD coenzyme; in full-length NOX the electrons are then transferred singly to the hemes and across the membrane to molecular O2. The functionality of the DH domain of SpNOX was evaluated by monitoring the oxidation of NADPH substrate and also by monitoring ferricytochrome c reduction which in the absence of the TM domain indicates the capacity of FAD to transfer electrons to soluble O2. Superoxide generation was assessed in aerobic conditions using SOD-inhibitable ferricytochrome c reduction.
Like purified or membrane preparations of both the full-length and DH constructs of other NOX family members (Han et al., 2001; Nisimoto et al., 2004; Nisimoto et al., 2010), SpNOXDH requires supplemental FAD to sustain both superoxide production, which can be observed in the presence of Cyt c (Figure 3A), and NADPH oxidation, which can be observed in the absence of Cyt c (Figure 3B). The dependence of SpNOXDH on flavin in the oxidation of NADPH was assessed by addition of diphenyleneiodonium (DPI), a common specific flavoenzyme inhibitor (Nisimoto et al., 2014), which inhibited the enzyme activity (Figure 3C).
In Figure 3D NADPH oxidation and superoxide generation by SpNOXDH in aerobic conditions was followed simultaneously at both 340 and 550 nm. The addition of superoxide dismutase (SOD) led to a ~40% inhibition of SpNOXDH’s Cyt. c reductase activity, confirming that the reduction of Cyt. c observed in the system is partly due to O2●- generation (Figure 3D and Table 1). However, the remaining~60%Cyt. c reductase activity suggests direct electron transfer from the FADH2 produced by the SpNOXDH to the Cyt. c. Thus, SpNOXDH, separated from the TM domain, retained the redox properties and sensitivity displayed in the context of the whole SpNOX protein (Hajjar et al., 2017). The strict requirement of FAD addition for SpNOXDH activity and its µM level of affinity suggests that the flavin behaves as a co-substrate rather than a prosthetic group. As an isolated domain, SpNOXDH may work as a flavin reductase enzyme (Gaudu et al., 1994; Fieschi et al., 1995; Nivière et al., 1996), making an interesting parallel with the E. coli flavin reductase fre (also in the FNR family) which reduces the MsrQ transmembrane component (in the FRD family, like the NOX TM domain, see Figure 7; Caux et al., 2021; Juillan-Binard et al., 2017).
Comparative flavin recognition properties and activities of SpNOXDH vs SpNOX
To characterize the flavin reductase activity of the DH domain and the relative contributions to FAD binding of the TM and DH domains of SpNOX, we determined the Km and kcat of SpNOX and SpNOXDH cytochrome c reductase activity using flavins FAD, FMN, riboflavin, and lumiflavin (Table 2). We observed a progressive loss of affinity for SpNOX as flavin size diminishes, with 263-fold higher affinity for FAD than for lumiflavin (Table 2). In contrast, the SpNOXDH retains constant affinity equal to that of lumiflavin (~10 μM) regardless of the flavin analog. These results suggest that the isolated DH domain binds to FAD primarily at the isoalloxazine ring of the cofactor, and the addition of ribitol (riboflavin), phosphate groups (FMN), and adenosine nucleotide (FAD) do not further interact significantly with the isolated DH domain. However, the full-length SpNOX appears to interact along the entire length of the FAD cofactor, indicating significant interactions between the FAD and the TM domain. This FAD binding site shared between the TM and DH domains may contribute to the correct placement of the cofactor for electron transfer to the proximal heme. On the other hand, in full-length SpNOX the kcat appears to be independent of the flavin derivative used, suggesting that the binding interactions between TM and FAD play no catalytic role and that only the isoalloxazine ring needs to be taken into account for the electron pathway from NADPH to the first heme. Finally, the kcat of full-length and DH-only constructs are the same order of magnitude, with ratios around 1 (Table 2), showing that the presence (or absence) of the hemes in the TM domain has little effect on the rate. It follows that the reduction of FAD – that is, hydride transfer from the nicotinamide – is the rate limiting step in the electron transfer from substrate to O2.
SpNOX uses NADPH and NADH with similar affinity and activity
The typical Michaelis Menten type kinetics parameters of SpNOX and SpNOXDH were obtained for increasing NAD(P)H concentrations (Table 3). Full length and DH constructs showed the same order of magnitude affinity to each other for NADPH or NADH as substrate (Table 3). Furthermore, the value of 1 for NADPH/NADH affinity for each construct shows that SpNOX does not discriminate significantly between NADPH and NADH; the agreement of two different assays – cytochrome c reduction and NADPH oxidation – strengthens this observation. Thus, in SpNOX, the NAD(P)H binding energy arises from interactions just in the DH domain with little contribution from the TM domain of SpNOX. The lack of discrimination of SpNOX between NADH and NADPH represents a significant difference from eukaryotic NOX which strongly select NADPH as substrate. In these experiments that used a range of NAD(P)H concentration, kcat’s of the same order of magnitude for both constructs again supports hydride transfer as the limiting step in NOX electron flow.
Structural characterization of SpNOX and SpNOXDH : mutation of F397
We obtained large hexagonal crystals of SpNOX but with poor diffraction power beyond 5 Å, and only 6% between 3 and 5 Å, with clear anisotropy. We therefore tried stabilizing mutations based on the work of Deng et al., 1999 in pea FNR, homologous to the NOX DH domain. They showed that the FNR C-terminal Y308 stabilizes the isoalloxazine site of the FAD cofactor and regulates access of the NADP + substrate to this site; other works extended the discussion of the importance of this residue to the NOX DH domain (Zhen et al., 1998; Kean et al., 2017). We reasoned that changing the corresponding residue of SpNOX to W would stabilize FAD binding and inhibit NADPH binding, thereby reducing mobility of the two subdomains relative to each other. Supporting this approach, Magnani et al reported the stabilizing effect in CsNOX of the serendipitous insertion of a Trp that stacks with the FAD isoalloxazine ring (Ceccon et al., 2017; Magnani et al., 2017). Similarly, we reasoned that changing the corresponding F residue to S would increase NADPH access to FAD, potentially facilitating crystals with NADPH bound.
Based on sequence alignment (Figure 1), either of two phenylalanine residues in SpNox, F397, and F399, could correspond to the conserved C-terminal residue in FNR. We mutated both positions to tryptophan in the SpNOXDH construct and monitored activity to test their effects. In our standard buffer used for comparing DH with full-length constructs, F399W mutants showed normal activity (Table 3), while F397W mutants showed no activity (Figure 4). We interpreted this to mean that W in position 397 prevents hydride transfer by immobile stacking with the isoalloxazine ring of FAD, while W in position 399 has little effect on hydride transfer. These results indicate that F397 of SpNOX corresponds to the terminal Y of FNR.
We subsequently explored mutations of F397 in both full length and DH constructs of SpNOX. The flavin content of purified F397W is much higher than that of WT, while the F397S is not able to retain flavin in its binding site (Figure 4A), consistent with our hypothesis about flavin binding. As displayed in Figure 4B, in our standard activity conditions, F397W has little or no activity while F397S shows enhanced activity compared to WT. Increasing the NADPH concentration 10-fold resulted in low but significant activity of the F397W mutant (not shown). These results indicate that the W and S mutations of position 397 do not significantly affect the structural integrity of the SpNOX constructs. We used Michaelis-Menten experiments to test the affinity for both NADPH and NADH of F397S mutants of full-length enzyme and DH constructs; such experiments were not possible with F397W because of its very low activity (Figure 4). F397S mutants show similar results to WT (Table 3) in both constructs indicating that S in this position does not greatly impact affinity for either substrate; this mutation produced an increase in kcat for both constructs, consistent with the hypothesis that S in this position provides increased access of substrate to the catalytic site.
High resolution structure of the SpNOX dehydrogenase domain
The use of the full-length SpNOX F397W construct resulted in about 20% of the crystals produced displaying diffraction limits close to 4 Å, with some crystals showing an anisotropic diffraction between 3 and 3.3 Å in a privileged direction of the crystals, as often observed for membrane proteins (Martin et al., 2023; Robert et al., 2017). Despite the improvement of the crystal quality, the anisotropy precluded resolution without a good model for molecular replacement. However, the soluble SpNOXDH domain represents more than 50% of the full-length enzyme, making it a suitable model for use in molecular replacement. We obtained some crystals of both SpNOXDH WT and F397W, co-crystallized with FAD (Figure 5—figure supplement 1). Unfortunately, crystals of the SpNOX F397S mutant with or without NADPH could not be obtained.
De novo resolution at 1.94 Å of SpNOXDH F397W was possible because we co-crystallized in the presence of sodium bromide and collected a bromide-SAD dataset (see methods). To solve the SpNOXDH WT structure, we used the SpNOXDH F397W structure as a model, and achieved 2.5 Å resolution. As expected, SpNOXDH presents a typical fold of the FNR superfamily of reductase domain containing two sub-domains, the FAD-binding domain (FBD) and an NADPH-binding domain (NBD) (Figure 5A). Both WT and F397W SpNOXDH display 3 molecules per asymmetric unit and the two structures are identical. In the asymmetric unit of SpNOXDH WT, one of the three molecules contains no flavin (Figure 5—figure supplement 2.). In the asymmetric unit of SpNOXDH F397W, all molecules of the unit contain a well-defined FAD (Figure 5—figure supplement 2). These observations agree with flavin contents observed in solution (Figure 4A) and support the hypothesis that the F397W mutant stabilizes FAD binding.
In strong correlation with the affinity data obtained for the flavins shown in Table 2, the DH domain mainly interacts with the isoalloxazine ring through polar contacts with the protein backbone and a few residues (P234, F235, and S23) belonging to the conserved FAD-1 binding motif (Massari et al., 2022; see Figure 1). The DH domain shows almost no contact with the ribitol and adenosine parts of the FAD. With the exception of possible polar contact of the first phosphate group with the backbone of the protein, the rest of the FAD molecule points out into the solvent, in the direction of the missing TM domain (Figure 5D).
In the wild type, the isoalloxazine ring is sandwiched between FBD residues P234 and F235 on its si-face, and the NBD residue F397 on its re-face (Figure 5D). Superposition of the different SpNOXDH WT molecules from the crystal asymmetric unit shows that the F397 aromatic ring imperfectly occupies the position in front of the re-face of the isoalloxazine ring (Figure 5B), with orientation 40° with respect to the plane of the isoalloxazine ring, and distance at 5.3 Å. In contrast (and as predicted from the homologous FNR structure), the tryptophan of the F397W mutation makes an extensive and much tighter parallel stacking interaction at a distance of 3.5 Å from the isoalloxazine ring (Figure 5C). From these observations, WT SpNOXDH could well behave as a flavin reductase, able to use different flavins as substrate through the sole recognition of the isoalloxazine ring. As Figure 4A shows, the F397W mutant may fix the flavin as a cofactor. Figure 5E shows the sidechains of both the F and W residue overlaid, on which is superimposed a putative binding mode of NADPH based on the pea FNR structure (Deng et al., 1999). Trp at position 397 appears to preclude access of the nicotinamide to the isoalloxazine, preventing direct hydride transfer. Phe in this position, however, leaves more space in this pocket and the sidechain has more degrees of freedom, making it feasible for the nicotinamide to displace it and transfer hydride to the isoalloxazine. These observations most likely explain the drastic difference of activity between WT and F397W; they probably also explain the enhanced activity observed in the F397S mutant which provides even more space near the isoalloxazine (Table 3). The increase of the global kcat in the F397S mutant is an additional argument to support hydride transfer as the limiting step in NOX electron transfer.
Structure of full-length SpNOX
The high-resolution model of the DH domain of SpNOX allowed resolution of the anisotropic crystals of full-length SpNOX-F397W to 3.6 Å. Small extracellular loops limiting crystal contacts together with intrinsic interdomain flexibility of SpNOX are probably the source of the crystal anisotropy due to subtle variation of packing within the crystals (Figure 6—figure supplement 1). Like the previously reported structures of the NOX family (Liu et al., 2022; Magnani et al., 2017; Sun, 2020; Wu et al., 2021), the crystallographic model of SpNOX harbors the N-terminal TM domain that encompasses 6 membrane-spanning alpha helices chelating two hemes, and the C-terminal cytosolic DH domain (Figure 6A). The density corresponding to the two heme groups is packed in the channel formed by TM helices 2–5. The inner and outer heme are respectively coordinated by the imidazole rings of pairs of histidines H69 (TM3)/H129 (TM5) and H82 (TM3)/H142 (TM5), holding them perpendicular to the plane of the membrane (Vermot et al., 2020). The metal-to-metal distance of the two hemes is ~21.5 Å, while the shortest interatomic distance is ~9.8 Å, between the vinyl groups.
As observed in previously solved structures of NOX enzymes (Liu et al., 2022; Magnani et al., 2017; Sun, 2020; Wu et al., 2021), a similar electron transfer relay across the membrane was mapped. SpNOX has 9.2 Å edge-to-edge distance between the isoalloxazine ring to the inner heme; a conserved aromatic residue, Y136, at equal distance from both hemes (6.9 Å), or even F107 also on the pathway to the hemes, could assist the electron transfer between the hemes as previously suggested (Magnani et al., 2017).
FAD binding site in the full-length SpNOX: role of the D-Loop region
We observed electronic density corresponding to the FAD cofactor at the interface of the DH and TM domain. FAD binding in the full-length construct preserves the π-stacking interaction with W397 and the interactions with F235 and P234 observed in the DH-only structure. In addition, the full-length SpNOX structure shows that the adenine extremity of FAD is packed in a pocket including the amino acid F399 from the DH domain– providing a potential role for this residue – as well as residues from the D-loop of the TM domain. Indeed, Y122 and R126 could stabilize the adenine ring simultaneously. Y122 interacts through π-stacking with the adenine ring (Figure 6B). R126 is at a reasonable distance to bridge with the adenine ring and simultaneously a carboxylate of the proximal heme. However, the resolution of the structure makes the accurate positioning of side chains difficult, particularly in this region which displays large electron densities.
To check the validity of our side-chain observations we generated alanine mutants of Y122 and R126 and checked the affinity of several flavins for each mutant (Figure 6B and Table 4). We determined affinity for NADPH for all mutants; similar affinity in all mutants suggests no effect of these mutations on global enzyme structure (Table 4). The results in Table 4 shows that only the affinity for FAD is impacted by Y122A and R126A, confirming the role of these residues in the interaction with the FAD’s adenine ring, which is missing in the other flavins tested. Y122A diminishes FAD affinity by a factor of 6 and R126A by a factor of 23. This latter suggests that R126 is crucial for FAD affinity; possibly the simultaneous interaction between R125 and both heme and flavin (proposed above) could hold the adenine ring in a position to be sandwiched between Y122 on one side and F399 on the other. Surprisingly, while the position corresponding to Y122 is very highly conserved as an aromatic residue in the full NOX enzyme family, the position corresponding to SpNOX R126 shows a difference between eukaryotes and prokaryotes. In eukaryotes, the position corresponding to R126 is very highly conserved as an aromatic, while in prokaryotes it is most often a basic residue, although aromatics are not uncommon. The interactions with the TM domain’s D-loop (or lack thereof) provide an explanation for the differences among the affinities of the flavins for the full-length protein (Table 2).
Access to the distal heme: diverse situations from prokaryotic to eukaryotic enzymes
The FRD found in NOXs but also in MsrQ and STEAP proteins, transfers electrons from the cytoplasmic to the extracellular side of the cell. Thus, access of the distal heme is essential to the specificity of the electron acceptor and its product. Figure 7 compares the distal heme access of a panel of FRD-containing proteins from prokaryotes and eukaryotes. From the structures available to date, including the present study, some FRDs (MsrQ, SpNOX, and CsNOX) exhibit small extracellular loops and free access to the distal heme (Figure 7A), while others (NOX2) have long extracellular loops (Figure 7A) and sometimes another extracellular domain (DUOX) (Figure 7B) create a cap that limits access. Homology models of both NOX4 (O’Neill et al., 2018) and NOX1 (Ward et al., 2023) also predict similar caps. In NOX with large caps, the distal hemes are deeply masked, implying that O2 must diffuse into the protein. On the basis of structures, several potential tunnels have been proposed for NOX2 and DUOXs (Liu et al., 2022; Noreng et al., 2022; Sun, 2020; Wu et al., 2021) and for the homology model of NOX1 (Ward et al., 2023). For most of the tunnels identified, bottleneck constriction would require protein dynamics to allow diffusion and the exact pathway for O2 still remains to be defined. However, a potential entry site with direct access to the distal heme from the lateral side can be seen in the NOXs with large caps (Figure 7B). In addition, an internal pocket in close proximity to the distal heme has been proposed as the site of O2 reduction through an outer sphere electron transfer mechanism. This pocket is characterized in eukaryotic enzymes by a conserved His and Arg involved in O2 binding and heme interaction; the residues characteristic of the O2 reducing center in eukaryotic FRD domains of NOX and DUOX enzymes are not conserved in SpNOX. The importance of the conserved His and Arg in the reducing center pocket of eukaryotes are strongly supported by CGD-induced mutations (Cross et al., 1995; Picciocchi et al., 2011; Rae et al., 1998). A pathological NOX1 variant was used to show that an additional Asn residue, separated from the conserved His by two residues, is necessary for efficient O2 reduction Ward et al., 2023; this Asn is highly conserved in animal NOX (see Helix 3 motif in Figure 1) but not prokaryotes. Entry site, tunnel size and reducing center pocket location ensure the selectivity to O2 as unique electron acceptor for NOXs with large caps.
FRD containing proteins with direct access to the distal heme may reflect more versatility with regard to final electron acceptors. MsrQ, to date the only FRD-containing prokaryotic protein with clearly identified physiological function, provides an interesting case. In E. coli, MsrQ transfers its electrons to MsrP, a periplasmic methionine sulfoxide reductase, which repairs proteins damaged by methionine oxidation (Gennaris et al., 2015). The absence of long extracellular loops allows MsrP to dock on MsrQ and to receive electrons directly from the distal heme to its molybdopterin center. Similarly, SpNOX’s exposed heme may confer its ability to directly reduce Cyt. c (Hajjar et al., 2017). The physiological functions of neither SpNOX nor CsNOX are known; in addition to their NOX activity, they may use electron acceptors other than O2.
Domain orientation in SpNOX structure: comparison with other NOX structures
The previous crystal structure of a prokaryotic NOX was in fact two independent structures of the separate TM and DH domains of CsNOX (Magnani et al., 2017). With the SpNOX structure presented here we provide the first entire prokaryotic NOX structure, affording a view of the orientation of the TM domain relative to the DH. The importance of the relative orientation of the domains in NOX proteins has been underlined as a critical point in the activation and the efficiency of electron transfer (Liu et al., 2022; Wu et al., 2021). Comparison of the recent structures of an inactive ‘resting state’ NOX2 and of the presumably active ‘high Ca2+ state’ DUOX1, emphasized the large 49° rotation of the DH domain of one structure relative to the other. This motion leads to a 12.4 Å movement of the FAD, which was suggested as critical for the activation of the electron transfer (Liu et al., 2022). Including SpNOX into this movement analysis is all the more interesting given that SpNOX is constitutively active. Here, Figure 8A shows that in our SpNOX structure, the relative position of TM vs DH domain is intermediate between the inactive (NOX2) and active (DUOX) structures. Figure 8B details the potential intermediate motions of SpNOX: the relative position of SpNOX DH with respect to the inactive state NOX2 DH rotates 13° and translates 5 Å toward the plane of membrane, while the relative position of the high Ca2+ ‘activated’ state of DUOX1 DH rotates almost an additional 49° and translates an additional 7.6 Å. SpNOX’s domain orientation is much closer to that observed in the NOX2 resting state than that of DUOX1.
The relative positions of the cofactors involved in the first steps of electron transfer are highlighted in Figure 8C. The distance between hemes in a given FRD domain is invariant; thus, the efficiency of electron transfer depends on two main factors: the distance between the nicotinamide of NAD(P)H and the isoalloxazine ring of FAD, and the distance between the isoalloxazine ring of FAD and the proximal heme. Surprisingly, there is no evident correlation between DH motion from inactive to active state and contraction of the distance between FAD and heme. In fact, this distance actually increases between inactive NOX2 and constitutively active SpNOX, and from SpNOX to ‘activated’ (high Ca2+) DUOX. Nonetheless, in all cases these distances, ranging from 7.2 to 10.2 from edge to edge (from methyl group of the flavin ring to the carboxylate of the heme), all remain compatible with biologically relevant electron tunneling probability, calculated as <14 Å. (Page et al., 1999).
These observations led us to question the paradigm that switching between resting and activated states requires DH motion. The lack of DH motion between the low and the high Ca2+ state human DUOX1 structures adds further hesitation. The DH domains of these two structures superimpose with a RMSD of 0.246 Å, and the FAD and NADPH of these structures also superimpose perfectly. The distance between the nicotinamide ring of the NADPH and the isoalloxazine ring of the FAD even in the activated DUOX1 structure is surprisingly large for efficient hydride transfer (Figure 8D). In line with conclusions from DH homologs (Deng et al., 1999; Hermoso et al., 2002), our data (Table 3, Figures 4 and 8) suggests that the SpNOX F397W sidechain represents the position that would be occupied by the nicotinamide ring of NADPH. The distance between the Trp397 sidechain and the isoalloxazine ring of FAD in our structure is 3.3 Å, similar to the corresponding distance of these moieties in homologous DH domains and suitable for hydride transfer. In the case of DUOX1, it has been suggested (Wu et al., 2021) that to enhance hydride transfer, the ‘high-Ca2+ activated’ state requires ‘tensing’ of the DH from a relaxed orientation, which can be envisioned as the relative motion of the NBD subdomain (containing NADPH) toward the FBD (containing FAD), rather than a large swing or translation of the entire DH domain relative to the TM domain.
These considerations suggest that differences in relative domain orientations in NOXes might represent isoform specificity rather than intermediate activation states. If so, these differences and those of the cofactor distances could fine tune the electron transfer properties and resulting efficiency of ROS production, linking the isoform with its physiological context. In this context, the impact of the lipid environment on domain orientation is suggested by longstanding observations of lipid effects on NOX activity (Koshkin and Pick, 1993; Shpungin et al., 1989) and by the observation of a lipid in the murine DUOX structure that mediates the interaction between the TM and DH domain through a direct connection with the phospho-ADP-ribose part of the NADPH (Sun, 2020). Structures solved in detergent like ours presented here may be released from interaction with the membrane, increasing flexibility. The question of what happens in NOX2 –whether DH tensing, domain motion, or something else – remains elusive and awaits further structural characterization. In any case the SpNOX domain orientation presented here makes it a rather close model of human NOX2.
A close-up view of NOX’s NAD(P)H binding domains vs the FNR gold standard
The striking difference in distance between NADPH and FAD at the active site (Figure 8D) parallels a structural difference between NOX2 and DUOX1 on one hand, and SpNOX on the other. Here, we use the well-characterized FNR structure as a basis of comparison for DH domains. The NADPH binding groove has characteristic motifs universally conserved in DH domains (Piubelli et al., 2000; motifs NADPH1 and NADPH4 in Figure 1); these motifs provide specific contacts with the NADPH ligand, while the critical aromatic in the C-terminal motif is important for FAD stabilization and apparently displaced by NADPH binding (Deng et al., 1999; Kean et al., 2017). In FNR, the position of the NADPH1 and NADPH4 motifs – representative of the full binding groove for NADPH – are separated from the flavin binding site by a narrow space (Figure 9 top left), while in DUOX1 the corresponding space is much wider (Figure 9, top right). Superposition of the four DH domains makes it clear that in SpNOX the distance between the binding grooves for the two ligands is short, like in FNR (Figure 9 middle left), while the distance between the ligand grooves in NOX2 is longer, like in DUOX1 (Figure 9, middle right). NADPH-binding residues of FNR are mostly conserved in SpNOX (Figure 9—figure supplement 1), while the two Arg residues that interact with the 2’-phosphate of NADPH in DUOX 1 are conserved in NOX2, although the side chain of one of those (R446) is not visible in the CryoEM structure of NOX2 (Figure 9, bottom left).
Although the SpNOX and NOX2, in the resting state, structures lack NADPH, superposition produces an excellent picture of the location of the NADPH in these proteins. Figure 9 (bottom) shows that when the proteins are superimposed, the NADPH from pea FNR fits extremely well into the NADPH binding groove of SpNOX and reproduces the short distance between nicotinamide and isoalloxazine seen in FNR. Similarly, the superimposed DUOX1 NADPH fits very well into the corresponding binding groove of NOX2. Even this crude docking with no energy minimization leaves little doubt about the overall NADPH position in SpNOX and NOX2. Our functional characterization of mutants in position 397 strongly support the picture presented here. Of course, new structures for NOX2 and SpNOX in complex with NAD(P)H would be a welcome addition to confirm these observations.
Taking into account, the similar interdomain organization and comparable electron pathways (Figure 8D) of SpNOX and NOX2, the different topology at the interface of their respective FBD and NBD subdomains would explain the switch from an inactive conformation of NOX2 in its resting state to a constitutively active conformation of SpNOX. This idea is consistent with the idea that electron transfer in the high Ca2+-state DUOX1 structure may require a ‘tensing’ motion that brings the two subdomains, and hence their ligand binding grooves, closer together (Wu et al., 2021). The mechanism by which the inactive resting state of NOX2 changes to active, and how this situation fits with the activated state of the DUOX1 structure, remains to be elucidated.
Conclusion
The important roles of NOX in human, animal, and plant health and disease explain the field’s focus on eukaryotic NOX. However, prokaryotic NOX’s advantages in production and handling offer opportunities to decipher structural and mechanistic details relevant to the whole family (eg. Magnani et al., 2017; Vermot et al., 2020). We present here extensive enzymological characterization of SpNOX activity and also new structural information on both full-length enzyme and the isolated DH. First, the particular effort to optimize conditions for SpNOXDH allowed the investigation of the first step of NOX electron transport, hydride transfer from NAD(P)H to FAD (Figure 2 and Table 1). We showed that NADPH and NADH are equivalent electron donors for SpNOX, constituting a significant difference from eukaryotic NOX. Comparison with the full length SpNOX demonstrates that NAD(P)H binding is largely independent from the TM domain, at least in our conditions in which SpNOX is solubilized in LMNG detergent. Comparison of kcat obtained from both constructs suggests that the hydride transfer from NAD(P)H to FAD is the rate limiting step in the overall NOX reaction (Table 3). Structural and functional characterization of WT and mutant F397W and F397S highlight the critical role of the aromatic at 397: Phe provides flavin isoalloxazine ring interaction and stabilization and also properties that allow entrance of the nicotinamide ring to the active site, similar to the mechanism long accepted in FNR family domains. Parallel analysis of flavin analogs using both DH and full-length SpNOX shows that the DH domain alone has no selectivity for flavin type. SpNOX, on the other hand, displays a marked selectivity for FAD with affinity at the nM level, indicating binding contribution from the TM domain. This is supported by interactions visible in the entire SpNOX structure, and a related mutational analysis, highlighting the role of the D-Loop region in the interaction with the adenine ring of the flavin (Figure 6B).
Structural comparison of NOX structures indicates that accessibility of the distal heme of SpNOX and CsNOX may reflect plasticity in the final electron acceptor, while the inaccessibility of other NOXes may enforce O2 as a substrate (Figure 7). On the other hand, this comparison also generates a new conundrum, how to explain hydride transfer in activated states of NOX. SpNOX shares similar domain organization with NOX2, and the orientations and distances between cofactors in the electron transfer pathway are similar (Figure 9—figure supplement 2 for overlay of both structures). Surprisingly, DUOX1 domain organization and cofactor distances are quite different, raising the possibility that domain organization is isoform specific. Considering the constitutive activity of SpNOX, we believe that the electron pathway as organized in NOX2 is poised for electron transfer, and that the difference in position of the NADPH binding groove could explain the difference between active and inactive states. The ‘high Ca2+’ activated state of DUOX1 also displays the distant binding groove suggesting that an additional step might be required to reach an efficient active state. Rather than a global repositioning of the DH, activation of NOX2 may require an internal rearrangement within the DH that brings the NADPH groove closer to the FAD site; a similar rearrangement may be necessary for electron transfer in activated DUOX1. Conformational change within the DH may agree with a previous study concluding that p67phox, which activates NOX2, regulates hydride transfer from NADPH to the FAD rather than NADPH binding. (Nisimoto et al., 1999). In this view, the assembly of p67phox with the DH could promote the required rearrangement within the DH domain. However, the resting state of NOX2 has been reported to resist initial binding of NADPH (Liu et al., 2022). This could be a second level of control for NOX2.
During the review process of this work, another group published a pre-print presenting several CryoEM structures of SpNOX (Dubach et al., 2023). Remarkably, considering the small size of the protein and its solubilization in detergent, they solved several structures from 2.33 to 2.75 Å resolution. Their high-resolution in the TM domain allowed them to improve the orientation of the adenine part of FAD and of the distal heme. However, apart from these adjustments, their Cryo-EM and our crystal structures are essentially the same (residus Y122 of the D loop is stabilizing the adenine from FAD, the distances between the cofactors are the same as well as the access to the active site or the distal heme). The work of Dubach et al., 2023 includes both functional data and additional structural data that complements ours. For example, they mutated F107, one of the potential aromatic relays for the electron transfer between the hemes. This mutation did not produce any functional impact, supporting the idea of a direct electron transfer from the proximal to the distal heme, although only simultaneous F107 and Y136 mutations would definitively close the interheme ‘relay’ possibility. An important contribution of their Cryo-EM work includes structures of SpNOX with bound NAD(P)H; these fully validate the binding site we anticipated (Figure 9—figure supplement 1), and confirm the critical roles of residue K250, S348 and Y353 in NAD(P)H binding. In addition, these structures highlight the capacity of the active site to accommodate the F397 side chain, the nicotinamide ring and the flavin ring within the groove where hydride transfer should occur. Their structures document the role of F397 in positioning the nicotinamide core within this groove and in controlling hydride transfer, since the F sidechain is located between the flavin and nicotinamide rings. In one of the NADPH-containing structures, they used a F397A mutation, which mimics the situation in which the F sidechain has moved out from between the nicotinamide and flavin rings, to trap a nicotinamide conformation productive for hydride transfer, approaching 3.6 Å from the isoalloxazine ring. These data are fully consistent with our enzymatic and structural characterization of F397 mutants. Altogether both our structures and theirs strongly support the critical role of F397 in the tuning of the hydride transfer mechanism. Taken together, the observations from both the Cryo-EM and our crystal structures of SpNOX reinforce the need for close proximity of the NADPH and FAD binding grooves for hydride transfer, exemplified in this constitutively active enzyme; this in turn supports our hypothesis, presented herein, that the activated state of eukaryotic NOX requires a specific relative motion of the NBD toward the FBD to achieve a productive topology in the DH domain.
In conclusion, SpNOX seems to be a good model of an activated form of NOX2. The molecular tools and structures provided in this work will allow further investigations to understand additional molecular features at work in those fascinating NOX enzymes.
Between the last revision and the proof-reading process of this article, a CryoEM structure of an activated form of neutrophil NADPH oxidase was published in Nature by Liu et al., 2024. It illustrates and magnificently lends credence to all the hypotheses, arising from this work on SpNOX, showing a relative movement of the NBD towards the FBD subdomain and thus allowing productive binding of NADPH. This confirms SpNOX, constitutively active, as a model of the activated form of NOX2 and thus gives all the more importance to the mechanistic characterization presented in this article.
Materials and methods
Structure-based sequence alignment
Request a detailed protocolA large set of 549 NOX homolog sequences distributed about evenly between prokaryotes and eukaryotes, was aligned using MAFFT (Katoh et al., 2019). The sequences were split into transmembrane (TM) and dehydrogenase (DH) domains identified by homology and including 30 overlapping residues between the domains. Each domain was realigned separately with PromalS3D (Pei et al., 2008). The TM alignment was constrained to the crystal structure of CsNOX TM domain (5O0T) while the DH alignment was constrained to the crystal structure of CsNOX DH domain (5O0X) and DH domain homologs 1GJR, 1FNB, 1GVH, 2EIX, 3A1F; all other parameters were at default values. The resulting multiple alignments were inspected in Jalview (Waterhouse et al., 2009); those with missing secondary structural elements or with unusually large insertions were removed and remaining sequences were pruned to 91% (TM domain) or 95% (DH domain) identity. The final sequence sets (172 sequences of TM and 207 sequences of DH domain) were realigned in Promals3D with default parameters and the same structural constraints specified above. Logos of subsequences of the full alignment were created in WebLogo (Crooks et al., 2004).
Cloning of SpNOX, SpNOXDH WT, F397S and F397W
Request a detailed protocolThe synthetic SpNOX gene optimized for expression in E. coli including a polyhistidine tag and thrombin cleavage site at its N-terminus was amplified via PCR and subcloned into pET-30b (Novagen) between the NdeI/ BamHI sites (Vermot et al., 2020).
The expression vector encoding His-tagged SpNOXDH - aa 181–417 of full-length SpNOX - was obtained from SpNOX by site-directed mutagenesis according to the manufacturer’s protocol (QuikChange Lightning Site-Directed Mutagenesis Kit, Agilent), using forward primer 5’gtctggtcccgcgtggcagtaaaattagctttccgtatctggg3’ and reverse primer 5’cccagatacggaaagctaattttactgccacgcgggaccagac3’.
F397S or F397W mutants of SpNOX and SpNOXDH were created using forward primer 5’ACGGAACTGATCTACGAAGGCTCTAAATTCAAATGAGAATTCGAGC3’ and reverse primer 5’GCTCGAATTCTCATTTGAATTTAGAGCCTTCGTAGATCAGTTCCGT3’ for F397S or forward primer 5’CGGAACTGATCTACGAAGGCTGGAAATTCAAATGAGAATTCGAG3’ and reverse primer 5’CTCGAATTCTCATTTGAATTTCCAGCCTTCGTAGATCAGTTCCG3’ for F397W.
We retain the nomenclature ‘F397x’ for these mutations in both SpNOX and SpNOXDH constructs, although in SpNOXDH construct that residue is actually in position 217.
Production of SpNOX, SpNOXDH WT, F397S and F397W
Request a detailed protocolThe expression of SpNOX WT, F397S and F397W was performed as described in Vermot et al., 2020 for the full-length WT SpNOX.
The DH-only domain proteins SpNOXDH WT, -F397S and -F397W were over-expressed in E. coli BL21(DE3) using TB medium supplemented with 50 µg.ml–1 kanamycine. At culture OD600=1, expression was induced by the addition of 0.5 mM IPTG and the culture was grown overnight at 18 °C. Cells were harvested by centrifugation (5000 × g, 20 min, 4 °C) and stored at –20 °C. The purification procedure described below is the final optimized version for SpNOXDH stabilization. The pellet was resuspended in lysis buffer (50 mM TRIS pH7, 1 M NaCl, 10 mM Imidazole and 5% glycerol) supplemented with deoxyribonuclease and a mix of protease inhibitors (chymostatin, leupeptin, antipain, pepstatin at 1 µg.ml–1, aprotinin at 5 µg/ml), disrupted by sonication, and cleared by centrifugation at 4 °C (18,500 rpm, JA-25.50 Beckman, 20 min,). The clarified extract was purified on a 4 ml His60 Ni-IDA Superflow resin (Ozyme) previously equilibrated in high salt buffer (HS buffer; 50 mM TRIS pH7, 1 M NaCl, 10 mM Imidazole and 5% glycerol). After an extensive wash (20 column volumes, HS buffer), the protein was eluted with 50 mM TRIS pH7, 1 M NaCl, 300 mM Imidazole and 5% glycerol and the buffer was exchanged against a stabilization buffer (50 mM Bis- TRIS propane pH6.5, 0.3 M NaCl, 5% glycerol) via gravity using a desalting column (PD10 17-0851-01, GE Healthcare). The selected fractions were pooled and loaded onto a Superdex75 size exclusion chromatography (GE Healthcare) equilibrated in the same buffer. The protein was eluted as a monomer and concentrated by ultrafiltration (10 kDa centrifugal filter unit, Amicon). The protein purity was checked on SDS-PAGE and the concentration was determined by UV-visible spectroscopy using respectively ε280nm = 24870 M−1.cm−1 for SpNOXDH WT and SpNOXDH F397S or ε280nm = 30370 M−1.cm−1 for SpNOXDH F397W respectively.
Thermostability assays for SpNOXDH stability improvement
Request a detailed protocolThermal stability assays were performed using a Prometheus NT.48 nanoDSF instrument (NanoTemper Inc, Germany). Fifty μl of sample per condition were prepared to perform experiments in triplicates using standard-treated glass capillaries. The excitation at 280 nm was optimized to yield emission intensities of intrinsic fluorescence at 330 and 350 nm, in the absence of heat gradient, in the range of 5000–16,000 A.U. The temperature gradient was set to an increase of 1 °C/min in a range from 20 °C to 95 °C. Protein unfolding was measured by detecting the temperature-dependent change in the ratio of fluorescence emission at 350 nm vs 330 nm. Samples used for testing pH impact on stability were done by diluting 10 times concentrated SpNOXDH solution (2 mg.ml–1) to a 0.2 mg.ml–1 protein solution in the following 0.1 M buffers depending on the tested pH: citric acid buffer for pH 4, 5, 5.5 and 6; sodium acetate for pH 4.5; bis TRIS-propane for pH 6.5 and 8, Pipes for pH 7; TRIS for pH 7 and 7.5; Bicine buffer for pH 9 and CAPS buffer for pH 9.5.
Activity assay on purified SpNOX and SpNOXDH domain
Request a detailed protocolThe activity of the purified DH domain was measured under aerobic conditions following simultaneously the oxidation of NAD(P)H (decrease of absorbance at 340 nm, ε=6.22 M−1.cm−1) and the reduction of cytochrome c (increase of absorbance at 550 nm, ε=21.1 M−1.cm−1), as a function of time, using a Cary 50 UV-visible spectrophotometer (Varian).
A typical assay was as follows: 50 mM bis TRIS-propane pH 6.5, 0.003% LMNG supplemented with 200 µM NADPH, 300 mM NaCl, and 5% glycerol, was added to a 2 mm quartz cuvette. Then 0.5 µg.ml–1 SpNOXDH (final) and 10 µM FAD (final) were added, mixing after each addition. For Cyt c reductase assay, 100 µM cytochrome c was added prior to the addition of the protein. Final assay volume was 500 µl. When full length SpNOX was used alone, not in parallel study with SpNOXDH, we used a buffer optimized for activity of the full length enzyme: Tris pH 7, 300 mM NaCl, and 0.003% LMNG.
The inhibition of superoxide production was verified by addition of DPI at a final concentration of 50 µM or Superoxide dismutase (10 U) to the reaction mixture after a sufficient time to allow the slope to be characterized in the presence or absence of the inhibitor.
Data analysis of thermostability and activity assays
Request a detailed protocolThe results from thermostability and activity assays were analyzed simultaneously using the JMP statistical software (Version 13.0. SAS Institute Inc, Cary, NC, 1989–2019) to plot 2D contour plots of Tm and activity as a function of salt concentrations and glycerol. The median of triplicates was used to compute the graphs, and linear interpolations were done between measurement points to fill the graphs. A 2D linear interpolation was performed on thermostability and activity results as a function of pH and salt concentration using Python to produce a 3D figure.
Michaelis-Menten kinetic analysis for SpNOX and SpNOXDH
Request a detailed protocolAll chemicals and reagents were purchased from Sigma-Aldrich. Various NAD(P)H concentrations were obtained by serial dilutions from the saturating NAD(P)H solution in triplicate in a 96-well (flat bottom) Greiner microplate (equilibrated at 25 °C). The reaction was initiated (at 25 °C) by injection of a mixture of the stabilization buffer (50 mM bis TRIS-propane pH6.5, 300 mM NaCl, 10 µM FAD, 5% glycerol and 0.025 mM LMNG) containing SpNoxDH (250 ng) and cytochrome c (100 µM, final concentration) or (50 mM Tris pH 7, 300 mM NaCl, 10 µM FAD and 0.025 mM LMNG) containing SpNOX (250 ng) and cytochrome c (100 µM, final concentration) and the absorbance at 550 nm was recorded within 3 min using a monochromator multimode microplate reader (Clariostar, BMG Labtech) after orbital shaking. Michaelis-Menten saturation curves of SpNOX, initial apparent velocities as a function of NAD(P)H concentration, Km (µM), and kcat (s–1) were determined with Prism software. We determined the kinetic parameters for all flavin analogs using a similar procedure, holding NAD(P)H concentration constant at 200 µM and varying the flavin concentrations; similarly, kinetic parameters of SpNOXDH were determined for NAD(P)H based on NAD(P)H oxidation by monitoring the absorbance at 340 nm within 3 min. Affinity is formally described by Kd, but in the text we use Km as a ‘relevant approximation of affinity’ (Srinivasan, 2022).
Crystallization of SpNOXDH domains
Request a detailed protocolCrystallization experiments were carried out at the High Throughput Crystallization Laboratory (HTX Lab) at EMBL Grenoble using automated protocols (Cornaciu et al., 2021; Dimasi et al., 2007; Dupeux et al., 2011; Mariaule et al., 2014). SpNOXDH and SpNOXDH F397W were concentrated to 15 mg.ml–1. The crystallization experiments were carried out using the sitting-drop vapor-diffusion method with a crystallization robot (Mosquito, SPTLabtech). A total of 0.1 µL of protein solution and 0.1 µL of reservoir were mixed to equilibrate against 45 µL reservoir solution at 20 °C in 96-well CrystalDirect plates (MiTeGen) and automatically imaged in RockImager robot (Formulatrix). Initially, commercially available screening kits from NeXtal (Classic-suite, PEGs, JCSG) and Molecular Dimension (MemTrans, MemMeso, MemGold, PACT, Morpheus) were used to identify initial crystallization conditions that were optimized in subsequent experiments. Optimal crystals of SpNOXDH grew in crystallization drops with precipitant conditions of 35% PEG-MME500, 0.1 M Sodium Citrate pH5. Optimal crystals of SpNOXDHF397W grew in crystallization drops with precipitant condition composed of 17% to 20% w/v PEG3350, 0.1 M Bis-TRIS propane pH 6.5 and 0.2 M Sodium bromide. Crystals of SpNOXDH WT and SpNOXDH F397W typically appeared after 24 hr and reached full size in 7 days with dimensions of respectively 100×20 × 10 μm3 and 80×30 × 20 μm3. Automated high-throughput crystal cryo-cooling and harvesting were performed with CrystalDirect Technology (Cipriani et al., 2012; Márquez and Cipriani, 2014; Zander et al., 2016).
Data collection and structural solution of SpNOXDH F397W
Request a detailed protocolA bromide-SAD dataset was collected at the macromolecular crystallography beamline ID30B of the European Synchrotron radiation facility, equipped with Pilatus3 6 M detector (McCarthy et al., 2018). The wavelength was set above the Br K-edge at 0.9184 Å in order to maximize the anomalous signal, and a highly redundant data set at a resolution of 1.94 Å, including 7200 images with 0.05 degrees oscillation, was collected.
Data processing was performed by the GrenADES parallelproc pipeline (Monaco et al., 2013) that is based on XDS (Kabsch, 2010) and scaled using aimless Evans and Murshudov, 2013 from the CCP4-suite (Winn et al., 2011). Attempts to solve the structure with the molecular replacement method using homologous models failed. Therefore, we experimentally solved the structure by the SAD method, using the bromide anomalous dispersion signal from crystals grown in presence of 0.2 M sodium bromide. Experimental phasing and initial model building was performed using CRANK2 (Skubák and Pannu, 2013) The final model was obtained by successive rounds of manual and automated refinement using the programs REFMAC (Murshudov et al., 1997), and BUSTER (Bricogne et al., 2017). Crystallographic data statistics are presented in Table 5. The PDB deposition code for the SpNOXDH F397W is 8qq1.
Data collection and structural solution of SpNOXDH WT
Request a detailed protocolX-ray diffraction data were collected at ID30A-1/MASSIF-1 beamline at ESRF Grenoble (Bowler et al., 2015). A highly redundant data set at a resolution of 2.4 Å, including 3600 images with 0.1 degrees oscillation, was collected. The structure of SpNOX DH WT was solved at a 2.4 Å resolution by molecular replacement using the coordinates of the chain A of SpNOX DH F397W variant structure (PDB 8qq1) as a search template in Molrep searching for 3 mol/a.u. The initial solution was evaluated using COOT (Emsley et al., 2010) revealing no major clashes and good global packing with contacts between chains and classical size solvent channels. The final model was obtained by successive rounds of refinement using the programs REFMAC (Murshudov et al., 1997) and manual construction using COOT (Emsley et al., 2010). Crystallographic data statistics are presented in Table 5. The PDB deposition code for the SpNOXDH WT is 8qq5.
Crystallization, data collection, and processing of full-length SpNOX
Request a detailed protocolSpNOX F397W at 4.04 mg.ml–1 in 50 mM TRIS pH 7, 300 mM NaCl, 0.025 mM LMNG, 0.01 mM FAD was crystallized in 30.5% PEG 300, 0.15 M Li2SO4, 0.15 M NaCl and 0.1 M MES pH 6, using the vapor diffusion method in a 6 µL protein sample plus 6 µL of well hanging drop at 20 °C. Hexagonal brick-red crystals were obtained, mounted on mesh litho-loops and flash frozen in native conditions in liquid nitrogen. X-ray diffraction data were collected at ID30A-1/MASSIF-1 beamline at ESRF Grenoble (Bowler et al., 2015). A dataset of 1800 images was collected with an oscillation range of 0.2° per image and an exposure time of 1.699 s and a crystal to detector distance refined to 500.9 mm. As observed for many membrane protein crystals (Martin et al., 2023; Robert et al., 2017), the data collected were very anisotropic going to 3.2 Å resolution in the strongest diffracting direction. After data processing, the dataset was finally cut at 3.6 Å to ensure sufficient completeness in the highest resolution shell. Inclusion of further data only resulted in more noise without obvious gain in 2mFo-DFc electron density. Data were initially processed using the program XDS as spherical to the highest resolution possible (3.2 Å) even though spherical statistics were not usable. Staraniso analysis for diffraction anisotropy (staraniso@globalphasing.org) revealed that completeness was 85.4% in the highest resolution shell, therefore revealing that all the data collectable for this crystal had been collected. Data was cut at the diffraction limits suggested by the Staraniso server using only the first 1000 images.
Phases were solved by molecular replacement using PHASER (McCoy et al., 2007), on amplitudes, with data corrected for anisotropy using Staraniso. The search models were the high-resolution DH domain (PDB 8qq5) and a model of the transmembrane domain obtained using PROMALS3D (Pei et al., 2008) and I-TASSER (Yang and Zhang, 2015) on the basis of multiple sequence alignment of NOX homologs and validated by structural alignment on CsNOX (Vermot et al., 2020). The molecular replacement solution was achieved with placing the DH domain first, followed by the TM domain, yielding LLG = 91.8 and TFZ = 10.9. The phasing resulted in a model with one molecule of SpNOX in the asymmetric unit and the space group P6422. Crystal packing generation showed a sensible solution (classic type II packing of membrane protein and no clash), showing good density for the DH domain and clearly identifiable transmembrane helices, representing a typical initial anisotropic low-resolution map for a membrane protein solved by molecular replacement.
The resulting model was not biologically possible with an obvious error in the positioning of the transmembrane domain, which was inverted by 180° around the longitudinal axis. A new transmembrane domain model was generated on the basis of a sequence alignment between CsNOX and SpNOX and the use of CHAINSAW (CCP4). A new molecular replacement with PHASER (McCoy et al., 2007) without modification of the DH domain position resulted in a new orientation of transmembrane domain, which was biologically compatible. Fo-Fc map of the model harbored two positive peaks at 6.1σ and 3.6 σ corresponding to the position of proximal and distal hemes, respectively, sandwiched between two key histidine residues, allowing the sequence assignment of the TM domain. Iterative manual building in COOT (Emsley et al., 2010) and LORESTR pipeline (Kovalevskiy et al., 2016), using TLS (5–183 for TM domain and 184–399 for DH domain) allowed large repositioning of some α–helixes of transmembrane domain. A final comparison of a SpNox model generated with AlphaFold (Jumper et al., 2021) drove the repositioning of some loops and some sequence shifts in transmembrane helixes not associated with hemes. Final refinement was made with the help of ISOLDE (Croll, 2018) and LORESTR, ramachandran and rotamers outliers were corrected, yielding a final model with R=26.2% and Rfree = 32.0% that was deposited in the Protein Data Bank under the PDB accession code 8qq7. Crystallographic data statistics are presented in Table 5. The PDB deposition code for the SpNOX F397W is 8qq7.
Data availability
Diffraction data have been deposited in PDB under the accession codes 8qq1, 8qq5 and 8qq7. All other data generated or analyzed during this study are included in the manuscript and figure supplements.
-
RCSB Protein Data BankID 8QQ1. SpNOX dehydrogenase domain, mutant F397W in complex with Flavin adenine dinucleotide (FAD).
-
RCSB Protein Data BankID 8QQ5. Structure of WT SpNox DH domain: a bacterial NADPH oxidase.
-
RCSB Protein Data BankID 8QQ7. Structure of SpNOX: a Bacterial NADPH oxidase.
References
-
A Ca(2+)-activated NADPH oxidase in testis, spleen, and lymph nodesThe Journal of Biological Chemistry 276:37594–37601.https://doi.org/10.1074/jbc.M103034200
-
NOX3, a superoxide-generating NADPH oxidase of the inner earThe Journal of Biological Chemistry 279:46065–46072.https://doi.org/10.1074/jbc.M403046200
-
MASSIF-1: a beamline dedicated to the fully automatic characterization and data collection from crystals of biological macromoleculesJournal of Synchrotron Radiation 22:1540–1547.https://doi.org/10.1107/S1600577515016604
-
Engineering stability in NADPH oxidases: A common strategy for enzyme productionMolecular Membrane Biology 34:67–76.https://doi.org/10.1080/09687688.2018.1535141
-
CrystalDirect: a new method for automated crystal harvesting based on laser-induced photoablation of thin filmsActa Crystallographica Section D Biological Crystallography 68:1393–1399.https://doi.org/10.1107/S0907444912031459
-
The Automated Crystallography Pipelines at the EMBL HTX Facility in GrenobleJournal of Visualized Experiments 01:62491.https://doi.org/10.3791/62491
-
ISOLDE: a physically realistic environment for model building into low-resolution electron-density mapsActa Crystallographica. Section D, Structural Biology 74:519–530.https://doi.org/10.1107/S2059798318002425
-
A variant X-linked chronic granulomatous disease patient (X91+) with partially functional cytochrome bThe Journal of Biological Chemistry 270:8194–8200.https://doi.org/10.1074/jbc.270.14.8194
-
Cloning of two human thyroid cDNAs encoding new members of the NADPH oxidase familyThe Journal of Biological Chemistry 275:23227–23233.https://doi.org/10.1074/jbc.M000916200
-
Expression, crystallization and X-ray data collection from microcrystals of the extracellular domain of the human inhibitory receptor expressed on myeloid cells IREM-1Acta Crystallographica Section F Structural Biology and Crystallization Communications 63:204–208.https://doi.org/10.1107/S1744309107004903
-
A thermal stability assay can help to estimate the crystallization likelihood of biological samplesActa Crystallographica Section D Biological Crystallography 67:915–919.https://doi.org/10.1107/S0907444911036225
-
Purification of a Novel Flavoprotein Involved in the Thyroid NADPH OxidaseJournal of Biological Chemistry 274:37265–37269.https://doi.org/10.1074/jbc.274.52.37265
-
Features and development of CootActa Crystallographica. Section D, Biological Crystallography 66:486–501.https://doi.org/10.1107/S0907444910007493
-
How good are my data and what is the resolution?Acta Crystallographica. Section D, Biological Crystallography 69:1204–1214.https://doi.org/10.1107/S0907444913000061
-
The mechanism and substrate specificity of the NADPH:flavin oxidoreductase from Escherichia coliThe Journal of Biological Chemistry 270:30392–30400.https://doi.org/10.1074/jbc.270.51.30392
-
The NAD(P)H:flavin oxidoreductase from Escherichia coli as a source of superoxide radicalsThe Journal of Biological Chemistry 269:8182–8188.
-
A Two-component NADPH Oxidase (NOX)-like System in Bacteria Is Involved in the Electron Transfer Chain to the Methionine Sulfoxide Reductase MsrPThe Journal of Biological Chemistry 292:2485–2494.https://doi.org/10.1074/jbc.M116.752014
-
XDSActa Crystallographica. Section D, Biological Crystallography 66:125–132.https://doi.org/10.1107/S0907444909047337
-
MAFFT online service: multiple sequence alignment, interactive sequence choice and visualizationBriefings in Bioinformatics 20:1160–1166.https://doi.org/10.1093/bib/bbx108
-
Automated refinement of macromolecular structures at low resolution using prior informationActa Crystallographica. Section D, Structural Biology 72:1149–1161.https://doi.org/10.1107/S2059798316014534
-
NADPH oxidase: an enzyme for multicellularity?Trends in Microbiology 11:9–12.https://doi.org/10.1016/s0966-842x(02)00007-0
-
Estimation of crystallization likelihood through a fluorimetric thermal stability assayMethods in Molecular Biology 1091:189–195.https://doi.org/10.1007/978-1-62703-691-7_13
-
BookCrystaldirect: A novel approach for automated crystal harvesting based on Photoablation of thin filmsIn: Chen YW, editors. Structural Genomics: General Applications. Humana Press. pp. 197–203.https://doi.org/10.1007/978-1-62703-691-7_14
-
Specific Xray diffraction patterns of membrane proteins caused by secondary structure collinearityBiochimica et Biophysica Acta. Biomembranes 1865:184065.https://doi.org/10.1016/j.bbamem.2022.184065
-
ID30B - a versatile beamline for macromolecular crystallography experiments at the ESRFJournal of Synchrotron Radiation 25:1249–1260.https://doi.org/10.1107/S1600577518007166
-
Phaser crystallographic softwareJournal of Applied Crystallography 40:658–674.https://doi.org/10.1107/S0021889807021206
-
Automatic processing of macromolecular crystallography X-ray diffraction data at the ESRFJournal of Applied Crystallography 46:804–810.https://doi.org/10.1107/S0021889813006195
-
Refinement of macromolecular structures by the maximum-likelihood methodActa Crystallographica Section D Biological Crystallography 53:240–255.https://doi.org/10.1107/S0907444996012255
-
NOX5 in human spermatozoa: Expression, function, and regulationThe Journal of Biological Chemistry 287:9376–9388.https://doi.org/10.1074/jbc.M111.314955
-
The p67(phox) activation domain regulates electron flow from NADPH to flavin in flavocytochrome b(558)The Journal of Biological Chemistry 274:22999–23005.https://doi.org/10.1074/jbc.274.33.22999
-
Nox4: a hydrogen peroxide-generating oxygen sensorBiochemistry 53:5111–5120.https://doi.org/10.1021/bi500331y
-
Is the NAD(P)H:flavin oxidoreductase from Escherichia coli a member of the ferredoxin-NADP+ reductase family?. Evidence for the catalytic role of serine 49 residueThe Journal of Biological Chemistry 271:16656–16661.https://doi.org/10.1074/jbc.271.28.16656
-
Structure of the core human NADPH oxidase NOX2Nature Communications 13:6079.https://doi.org/10.1038/s41467-022-33711-0
-
Quantitative interaction analysis permits molecular insights into functional NOX4 NADPH oxidase heterodimer assemblyThe Journal of Biological Chemistry 293:8750–8760.https://doi.org/10.1074/jbc.RA117.001045
-
PROMALS3D: a tool for multiple protein sequence and structure alignmentsNucleic Acids Research 36:2295–2300.https://doi.org/10.1093/nar/gkn072
-
Role of putative second transmembrane region of Nox2 protein in the structural stability and electron transfer of the phagocytic NADPH oxidaseThe Journal of Biological Chemistry 286:28357–28369.https://doi.org/10.1074/jbc.M111.220418
-
Competition between C-terminal tyrosine and nicotinamide modulates pyridine nucleotide affinity and specificity in plant ferredoxin-NADP(+) reductaseThe Journal of Biological Chemistry 275:10472–10476.https://doi.org/10.1074/jbc.275.14.10472
-
X-Linked Chronic Granulomatous Disease: Mutations in the CYBB Gene Encoding the gp91-phox Component of Respiratory-Burst OxidaseThe American Journal of Human Genetics 62:1320–1331.https://doi.org/10.1086/301874
-
Activation of the Superoxide Forming NADPH Oxidase in a Cell-free System by Sodium Dodecyl SulfateJournal of Biological Chemistry 264:9195–9203.https://doi.org/10.1016/S0021-9258(18)60513-7
-
Automatic protein structure solution from weak X-ray dataNature Communications 4:2777.https://doi.org/10.1038/ncomms3777
-
A guide to the Michaelis-Menten equation: steady state and beyondThe FEBS Journal 289:6086–6098.https://doi.org/10.1111/febs.16124
-
Structures of mouse DUOX1-DUOXA1 provide mechanistic insights into enzyme activation and regulationNature Structural & Molecular Biology 27:1086–1093.https://doi.org/10.1038/s41594-020-0501-x
-
Overview of the CCP4 suite and current developmentsActa Crystallographica. Section D, Biological Crystallography 67:235–242.https://doi.org/10.1107/S0907444910045749
-
I-TASSER server: new development for protein structure and function predictionsNucleic Acids Research 43:W174–W181.https://doi.org/10.1093/nar/gkv342
-
Automated harvesting and processing of protein crystals through laser photoablationActa Crystallographica Section D Structural Biology 72:454–466.https://doi.org/10.1107/S2059798316000954
-
Probing the role of the carboxyl terminus of the gp91phox subunit of neutrophil flavocytochrome b558 using site-directed mutagenesisThe Journal of Biological Chemistry 273:6575–6581.https://doi.org/10.1074/jbc.273.11.6575
Article and author information
Author details
Funding
European Commission (iNEXT Discovery (grant agreement ID: 871037))
- Jose Antonio Marquez
European Commission (Fragment Screen (grant agreement ID: 101094131))
- Jose Antonio Marquez
Agence Nationale de la Recherche (FRISBI (ANR-10- INSB-0502))
- Franck Fieschi
Agence Nationale de la Recherche (GRAL (ANR-10-LABX-49-01))
- Franck Fieschi
Agence Nationale de la Recherche (Bandit (ANR17-CE11-0013))
- Jerome Dupuy
- Franck Fieschi
Université Grenoble Alpes (Emergence program)
- Susan ME Smith
- Franck Fieschi
Université Grenoble Alpes (Invited professorship)
- Susan ME Smith
The funders had no role in study design, data collection and interpretation, or the decision to submit the work for publication.
Acknowledgements
Funding for the automated crystallography pipelines at EMBL Grenoble was provided by the grants iNEXT Discovery (grant agreement ID: 871037) and Fragment Screen (grant agreement ID: 101094131) funded by the European Commission and through Instruct-ERIC. We want to thank the EMBL-ESRF Joint Structural Biology Group (JSBG) and in particular Andrew McCarthy for access and support in the use of MX beam lines at the ESRF. This work used the platforms of the Grenoble Instruct center (UMS 3518 CNRS-CEA-UJF-EMBL) with support from FRISBI (ANR-10- INSB-0502) and GRAL (ANR-10-LABX-49–01) within the Grenoble Partnership for Structural Biology. This work was supported by the French Agence Nationale de la Recherche Bandit Project (ANR17-CE11-0013) to FF and JD. AV was supported through the Emergence program from the Univ. Grenoble Alpes; SMES was supported by an invited professorship from Univ. Grenoble-Alpes and through Emergence partner Kennesaw State University.
Version history
- Preprint posted:
- Sent for peer review:
- Reviewed Preprint version 1:
- Reviewed Preprint version 2:
- Version of Record published:
Cite all versions
You can cite all versions using the DOI https://doi.org/10.7554/eLife.93759. This DOI represents all versions, and will always resolve to the latest one.
Copyright
© 2024, Petit-Hartlein et al.
This article is distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use and redistribution provided that the original author and source are credited.
Metrics
-
- 605
- views
-
- 113
- downloads
-
- 3
- citations
Views, downloads and citations are aggregated across all versions of this paper published by eLife.
Download links
Downloads (link to download the article as PDF)
Open citations (links to open the citations from this article in various online reference manager services)
Cite this article (links to download the citations from this article in formats compatible with various reference manager tools)
Further reading
-
- Biochemistry and Chemical Biology
- Structural Biology and Molecular Biophysics
Pre-mRNA splicing is catalyzed in two steps: 5ʹ splice site (SS) cleavage and exon ligation. A number of proteins transiently associate with spliceosomes to specifically impact these steps (first and second step factors). We recently identified Fyv6 (FAM192A in humans) as a second step factor in Saccharomyces cerevisiae; however, we did not determine how widespread Fyv6’s impact is on the transcriptome. To answer this question, we have used RNA sequencing (RNA-seq) to analyze changes in splicing. These results show that loss of Fyv6 results in activation of non-consensus, branch point (BP) proximal 3ʹ SS transcriptome-wide. To identify the molecular basis of these observations, we determined a high-resolution cryo-electron microscopy (cryo-EM) structure of a yeast product complex spliceosome containing Fyv6 at 2.3 Å. The structure reveals that Fyv6 is the only second step factor that contacts the Prp22 ATPase and that Fyv6 binding is mutually exclusive with that of the first step factor Yju2. We then use this structure to dissect Fyv6 functional domains and interpret results of a genetic screen for fyv6Δ suppressor mutations. The combined transcriptomic, structural, and genetic studies allow us to propose a model in which Yju2/Fyv6 exchange facilitates exon ligation and Fyv6 promotes usage of consensus, BP distal 3ʹ SS.
-
- Biochemistry and Chemical Biology
- Neuroscience
TIMM50, an essential TIM23 complex subunit, is suggested to facilitate the import of ~60% of the mitochondrial proteome. In this study, we characterized a TIMM50 disease-causing mutation in human fibroblasts and noted significant decreases in TIM23 core protein levels (TIMM50, TIMM17A/B, and TIMM23). Strikingly, TIMM50 deficiency had no impact on the steady-state levels of most of its putative substrates, suggesting that even low levels of a functional TIM23 complex are sufficient to maintain the majority of TIM23 complex-dependent mitochondrial proteome. As TIMM50 mutations have been linked to severe neurological phenotypes, we aimed to characterize TIMM50 defects in manipulated mammalian neurons. TIMM50 knockdown in mouse neurons had a minor effect on the steady state level of most of the mitochondrial proteome, supporting the results observed in patient fibroblasts. Amongst the few affected TIM23 substrates, a decrease in the steady state level of components of the intricate oxidative phosphorylation and mitochondrial ribosome complexes was evident. This led to declined respiration rates in fibroblasts and neurons, reduced cellular ATP levels, and defective mitochondrial trafficking in neuronal processes, possibly contributing to the developmental defects observed in patients with TIMM50 disease. Finally, increased electrical activity was observed in TIMM50 deficient mice neuronal cells, which correlated with reduced levels of KCNJ10 and KCNA2 plasma membrane potassium channels, likely underlying the patients’ epileptic phenotype.