Structure of dimerized assimilatory NADPH-dependent sulfite reductase reveals the minimal interface for diflavin reductase binding

preprint OA: closed
Full text JSON View at publisher

Abstract

Abstract Escherichia coli NADPH-dependent assimilatory sulfite reductase (SiR) fixes sulfur for incorporation into sulfur-containing biomolecules. SiR is composed of two subunits: an NADPH, FMN, and FAD-binding diflavin reductase and an iron siroheme/Fe4S4 cluster-containing oxidase. How they interact has been unknown for over 50 years because SiR is highly flexible, thus has been intransigent for traditional X-ray or cryo-EM structural analysis. A combination of the chameleon plunging system with a fluorinated lipid overcame the challenge of preserving a dimer between the subunits for high-resolution (2.84 Å) cryo-EM analysis. Here, we report the first structure of the reductase/oxidase complex, revealing how they interact in a minimal interface. Further, we determined the structural elements that discriminate between pairing a siroheme-containing oxidase with a diflavin reductase or a ferredoxin partner to channel the six electrons that reduce sulfite to sulfide.
Full text 134,317 characters · extracted from preprint-html · click to expand
Structure of dimerized assimilatory NADPH-dependent sulfite reductase reveals the minimal interface for diflavin reductase binding | Research Square window.SnipcartSettings = { analytics: { enabled: false } }; (function() { var accessVector = localStorage.getItem('access_vector') || ''; window.dataLayer = window.dataLayer || []; if (accessVector) { window.dataLayer.push({ user: { profile: { profileInfo: { snid: accessVector } } } }); } })(); (function(w,d,s,l,i){w[l]=w[l]||[];w[l].push({'gtm.start':new Date().getTime(),event:'gtm.js'});var f=d.getElementsByTagName(s)[0],j=d.createElement(s),dl=l!='dataLayer'?'&l='+l:'';j.async=true;j.src='https://www.googletagmanager.com/gtm.js?id='+i+dl;f.parentNode.insertBefore(j,f);})(window,document,'script','dataLayer','GTM-K279D39R'); Browse Preprints In Review Journals COVID-19 Preprints AJE Video Bytes Research Tools Research Promotion AJE Professional Editing AJE Rubriq About Preprint Platform In Review Editorial Policies Our Team Advisory Board Help Center Sign In Submit a Preprint Cite Share Download PDF Article Structure of dimerized assimilatory NADPH-dependent sulfite reductase reveals the minimal interface for diflavin reductase binding M. Elizabeth Stroupe, Behrouz Ghazi Esfahani, Nidhi Walia, Kasahun Neselu, and 4 more This is a preprint; it has not been peer reviewed by a journal. https://doi.org/ 10.21203/rs.3.rs-4758050/v1 This work is licensed under a CC BY 4.0 License Status: Published Journal Publication published 26 Mar, 2025 Read the published version in Nature Communications → Version 1 posted You are reading this latest preprint version Abstract Escherichia coli NADPH-dependent assimilatory sulfite reductase (SiR) fixes sulfur for incorporation into sulfur-containing biomolecules. SiR is composed of two subunits: an NADPH, FMN, and FAD-binding diflavin reductase and an iron siroheme/Fe 4 S 4 cluster-containing oxidase. How they interact has been unknown for over 50 years because SiR is highly flexible, thus has been intransigent for traditional X-ray or cryo-EM structural analysis. A combination of the chameleon plunging system with a fluorinated lipid overcame the challenge of preserving a dimer between the subunits for high-resolution (2.84 Å) cryo-EM analysis. Here, we report the first structure of the reductase/oxidase complex, revealing how they interact in a minimal interface. Further, we determined the structural elements that discriminate between pairing a siroheme-containing oxidase with a diflavin reductase or a ferredoxin partner to channel the six electrons that reduce sulfite to sulfide. Biological sciences/Structural biology Biological sciences/Structural biology/Electron microscopy/Cryoelectron microscopy NADPH-dependent assimilatory sulfite reductase hemoflavoprotein oxidoreductase cryo-EM Figures Figure 1 Figure 2 Figure 3 Figure 4 Introduction Assimilatory sulfite reduction by NADPH-dependent sulfite reductase (SiR) is essential to produce sulfide for incorporation into sulfur-containing biomolecules. In g-proteobacteria like Escherichia coli , SiR is a multimeric oxidoreductase composed of an octameric diflavin reductase (SiRFP) and four independently binding subunits of a siroheme-containing hemoprotein (SiRHP) 1-3 . The E. coli SiRFP subunit is homologous to cytochrome P450 (CYP) reductase (CPR) 4 , the reductase domain of the bacterial CYP/CPR fusion CYP102A1/CYPBM3 5 , the reductase domain of nitric oxide synthase (NOSr) 6,7 , and methionine synthase reductase (MSR) 8 . One of the hallmarks of this diflavin reductase family is that they are exceptionally conformationally malleable, which makes structural analysis challenging. For example, to date there are no high-resolution structures of the full-length NOS homodimer, the complex between methionine synthase and MSR, or the SiR heterododecameric holoenzyme (for simplicity, here referred to as a dodecamer). The structure of the CYP/CPR heterodimer and CYPBM3 are known 9-11 , but CYP/CPR form a 1:1 heterodimer, whereas the oxidase and reductase domains are fused in CYPBM3, so little can be inferred about other homologs that function as higher-order protein complexes. Despite its homology to other well-studied diflavin reductases, SiRFP is unique because it assembles into an octamer through its N-terminal 52 amino acids 12 . Removing those amino acids results in a 60 kDa monomer (SiRFP-60), which binds SiRHP as a 1:1 heterodimer with reduced activity 3 . Further, removing the complete N-terminal FMN-binding flavodoxin (Fld) domain results in a 43 kDa monomer that contains just the NADPH- and FAD-binding NADP + ferredoxin reductase (FNR) domain (SiRFP-43), which also binds SiRHP as a 1:1 heterodimer but is inactive for electron transfer 13 (Abbreviations and theoretical molecular weights are summarized in Table S1). SiRHP has few known homologs because of its unique siroheme/Fe 4 S 4 cluster assembly that form the sulfite-binding active site 14 . Assimilatory SiRs or siroheme-dependent nitrite reductases (NiRs) from other bacterial species or plants have a similar hemoprotein but use a transiently bound ferredoxin as their electron source 15-17 . SiRs that are responsible for energy conversion, dissimilatory sulfite reductases (DSRs), share a common siroheme binding fold but are heterotetrameric and are often fused to auxiliary domains 18 . Their electron donors are poorly understood. Here, we show the first high-resolution cryogenic electron microscopy (cryo-EM) structure of the minimal SiRFP/SiRHP dimer, which elucidates their binding interface to understand how SiRHP tightly binds SiRFP’s FNR domain. Results The SiRHP-SiRFP interaction is highly sensitive to cryogenic TEM (cryo-EM) preparation We determined the structure of the SiRFP/SiRHP dimer from three modified, minimal dimers, each of which is named by the change to SiRFP and its resulting molecular weight (Table S1 ). First, we truncated SiRFP to remove the N-terminal Fld domain (SiRFP-43/SiRHP). This is an inactive dimer, as the Fld domain is required for electron transfer. Nevertheless, the two subunits bind tightly and is the most simplified complex between SiRFP and SiRHP 19 . Second, we truncated both SiRFP’s N-terminal octamerization domain as well as the linker between the Fld and FNR domains to create a monomeric SiRFP that can be locked in an open position (SiRFP-60D/SiRHP) 20 . Third, we generated a variant of monomeric SiRFP-60 lacking reactive cysteines into which we engineered a disulfide bond between the Fld and FNR domains (SiRFP-60X/SiRHP) 21 . Each variant is highly sensitive to traditional blotting plunge-freezing methods for cryo-EM preservation. To overcome this sensitivity, we combined the protection of a high critical micelle concentration, fluorinated lipid, fos-choline-8 (FF8, Creative Biolabs, Shirley, NY, USA ) 22 , with the blot-free, rapid plunging afforded by the chameleon system (SPT Labtech, Melbourn, UK) 23 . This cryo-EM sample preparation helped us to retain each intact complex within near ideal ice thickness and avoid denaturation at the air water interface (Fig. S1 ). The smallest complex (SiRFP-43/SiRHP) showed well-aligned 2D class averages, however the 3D structure revealed structural anisotropy, either due to its small size/asymmetric geometry or from a preferred orientation, that limited high-resolution analysis despite the absence of mobile elements (Figs. S1A and S2). SiRFP-60D/SiRHP showed moderate-resolution density (3.54 Å) for the SiRFP FNR domain and SiRHP, however the N-terminal Fld domain was not visible (Figs. S1B and S3A). The 2.84 Å-resolution structure of SiRFP-60X/SiRHP revealed the most detail for SiRFP’s Fld and FNR domains, despite a lack of density for the linker between SiRFP’s domains in the highest-resolution reconstruction (Figs. 1 A, S1C, and S3B). High resolution features for each of the cofactors in both subunits supported this reported resolution (Fig. S4). Therefore, we analyzed the SiRFP-SiRHP interface for this construct in detail. SiRFP-SiRHP binding The SiR dodecameric holoenzyme is composed of oligomers of the dimers discussed here and is about 800 kDa in mass. Despite this large mass, the binding interface between the minimal SiRFP/SiRHP dimer is small, reaching the surface area of 1,138 Å 2 relative to the overall surface of 43,610 Å 2 . SiRHP alone has a solvent-exposed surface of 25,680 Å 2 . SiRFP-60X alone has a solvent-exposed surface of 21,930 Å 2 . That is, for a large complex only about 2.6% of the solvent-exposed surface is buried upon subunit binding. This is consistent with hydrogen-deuterium exchange data on the complex that reveals single, short peptides from each subunit that become occluded upon binding 3 , the sequences of which predicted the interface would be dominated by hydrophobic interactions (Fig. 1 B). The interface is governed by the N-terminus of SiRHP. The structure of this region is previously uncharacterized as it is proteolytically removed in the X-ray crystal structure of E. coli SiRHP 24 . These 80 amino acids follow the topology helix 1 - loop - helix 2 - turn - helix 3 - b-strand 1 - helix 4 - loop - b-strand 2 (Fig. 1 C). Only amino acids from the turn, helix 4, and surrounding loops directly interact with SiRFP. The regions that are N-terminal to the interface interact with domain 1 or the N-terminal half of the parachute domain ( i.e. the first sulfite or nitrite reductase repeat (S/NiRR) 24,25 ), breaking SiRHP’s pseudo two-fold symmetry (Figs. 1 A and S5A). An extension to the parachute domain that is not ordered in the original crystal structure, from amino acids 184–209, helps to hold the N-terminus in place (Fig. S5A). Those that are N-terminal to the interface approach the distal active site, but do not contribute significantly to anion or siroheme binding as they are held back by interactions to SiRFP, discussed below. Moving N-terminally along the peptide, it then turns back to form the loop that binds a pocket in SiRFP before moving away from SiRFP. The N-terminal most amino acids reach all the way to the other side of SiRHP, interacting with the N-terminus of the a-helix (h11) that precedes the linker that joins the two S/NiRRs and mimics the siroheme binding site (Fig. S5B) 24 . The central interaction that pegs the subunits together is a p-cation interaction between h-Lys27 from SiRHP (for simplicity, amino acids from SiRHP will be designated with the prefix “h-“) and f-His258 from SiRFP (similarly, amino acids from SiRFP will be designated with the prefix “f-“) (Figs. 1 D and E). This interaction is buttressed by f-Phe496 and f-Val500, which have previously been shown essential for SiRFP-SiRHP binding 13 . Further hydrophobic and p-stacking interactions dominate the interface. For example, h-Leu40 inserts into a pocket in SiRFP formed between f-b-sheet 17 and f-a-helix 18, which includes f-Phe496. h-Ile65 Cg2 sits 3.3 Å from the plane of f-Arg250’s guanidinium group (Figs. 1 D and F), which is rotated 90 o from its position in free SiRFP 20 . h-Gln72Cg also packs into a pocket formed by the backbone atoms of f-Ile247, f-Thr248, and f-Gly249, pinned in place by the h-Ile65/f-Arg250 and h-Lys73/f-His258 interactions. Farther from the interface, there is another stacking interaction between the guanidinium group from h-Arg63 and the h-Phe437 aromatic ring that stabilized the deformed helix that includes h-Ile65 (Figs. 1 D and G). h-Phe437 is adjacent, through h-A443, to SiRHP’s iron-sulfur cluster. In this way, the binding interface between SiRFP and SiRHP reaches to the SiRHP active site through a network of hydrophobic interactions (Fig. S6). Neither ionic interactions nor hydrogen bonds play a direct role in the interface, but rather stabilize the amino acids and structural elements that mediate the interface (Figs. 1 C, D and S6). For example, a hydrogen bond network from f-Thr404Og through f-Tyr498OH and finally to f-His258Nd1 position its imidazole ring for the interaction with h-Lys73. An ionic interaction between h-Lys127Nz and h-Asp38Od2 reach across the loop, presenting h-Leu40 to project into SiRFP’s pocket. An additional ionic interaction between h-Asp61d1 and h-Arg66NHd also stabilize the deformed helix that contains h-Ile65 and turn it towards f-Arg250. The SiRHP interaction with SiRFP differs from SIR interaction with ferredoxin g-proteobacteria couple a diflavin reductase, SiRFP, to SiRHP. In contrast, other organisms like Zea mays and Mycobacteria tuberculosis use a ferredoxin (Fd) as their reductase partner 15 , 17 . In Z. mays SIR, Fd bridges SIR’s C-terminal domain 2 to a loop between the first two b-strands (amino acids Asp110 to Asn118), positioning the Fd iron-sulfur cluster near the SIR metal sites 26 . The interaction is bolstered by en face stacking between Fd Tyr37 and SIR Arg324 (Fig. 2 A). In SiRHP, the equivalent element between the structurally conserved b-strands 1 and 2 is considerably longer, stretching from h-Asp62 to h-Arg77 and containing a short helix from h-Arg63 to h-Glu71 (Fig. 2 B). This loop contains the critical residues h-Gln72 and h-Lys73 that anchor the interaction with SiRFP’s FNR domain, which faces away from where Fd binds to the Z. mays homolog (Fig. 2 C). Further, the arginine in SIR that stacks with Fd Tyr37 is not conserved in SiRHP – the equivalent position is h-Gly262, despite the structural conservation of the loop between b-strands 7 and 8 that contributes to siroheme binding (Fig. 2 D). The SiRHP active site loop is locked in its closed conformation When bound to SiRFP, SiRHP’s anion binding loop (h-Asn149 to h-Arg153) is in its closed position, held in place by a long, through-space interaction between h-Arg53 and h-Asn149 (Fig. 3 A). h-Arg53 is, in turn, held in place by stacking between its guanidinium group and h-Tyr58. The ring of h-Tyr58 subsequently sits over the methyl group on the siroheme pyrroline A ring. The only other new protein/siroheme interaction identified in this now complete structure of SiRHP is an ionic bond between h-Gln60 and the propionyl group from siroheme pyrroline ring B. The siroheme is saddle-shaped, as in free SiRHP and unlike in dissimilatory SIR 27 – 30 . h-Arg153 is flipped away from the bound phosphate. The other three anion binding amino acids, h-Arg83, h-Lys215, and h-Lys217, remain largely unchanged from free SiRHP (Figs. 3 A-C) 24 . This conformation differs from the various redox and anion-bound structures of free SiRHP 24 , 31 – 33 . In the phosphate-bound, free SiRHP structure that lacks the N-terminal 80 amino acid extension 24 , the loop is flipped open such that h-Ala146-h-Ala148 are disordered (Fig. 3 B). Upon reduction and sulfite binding, h-Arg153 flips over to interact with the smaller anion and the loop becomes ordered 32 . h-Asn149 points away from the active site. In this way, SiRFP binding to SiRHP, with the ordering of SiRHP’s N-terminus, induces an intermediate structure with elements of both the oxidized, phosphate bound and the reduced, sulfite bound conformations (Fig. 3 D). In the original SiRHP X-ray crystal structure, the siroheme iron is significantly domed above the siroheme nitrogens, indicative of an oxidized Fe 3+ (Fig. 3 B) 24 . Subsequent chemical reduction experiments show the doming flattens upon conversion to Fe 2+ , commensurate with release of the phosphate to allow substrate binding (Fig. 3 B) 32 . Contemporary X-ray diffraction experiments show that this reduction is beam-induced, but within the constraints of the crystal the phosphate remains bound to the siroheme iron in the active site 33 . In this cryo-EM structure, the siroheme iron appears to be in the plane of the siroheme nitrogens, suggesting that it has also been reduced by the electron beam. The central density for the phosphate is 3.5 Å from the siroheme iron and its pyramidal shape is rotated such that the oxygen-iron bond is broken (Figs. 3 A and S4C). SiRFP is highly mobile Although the 2D class averages in all three datasets appeared to show little orientation preference with discernable features, further analysis revealed each to have unique properties related to the SiRFP variant used to generate the dimer (Table S1 , Figs. S1-3). SiRFP-43/SiRHP (107 kDa in mass): The 2D class averages for SiRFP-43/SiRHP appeared to show high-resolution features, however the initial 3D models were poorly aligned, likely due to a combination of small mass, preferred orientation, and limitations in grid preparation, so the refined volume did not achieve high resolution despite its simplified form. 3D variability analysis in CryoSPARC 34 did not reveal significant conformational mobility, however orientation analysis and calculation of the 3DFSC 35 showed the particles harbored a preferred orientation (Fig. S2). Nevertheless, the absence of SiRFP’s Fld domain did not alter the SiRFP-SiRHP interaction. SiRFP-60D/SiRHP (123 kDa): As with SiRFP-43/SiRHP, this dimer showed 2D class averages with high-resolution features for the core of the dimer (Figs. S1B and S3A). Nevertheless, the initial models were of inconsistent structure, so we checked for heterogeneity using “3D variability” in CryoSPARC, which revealed a distinctive degree of movement in the Fld domain. To better understand the degree of flexibility for the Fld domain from the SiRFP-60D/SiRHP dimer, we performed both “3D flex” in CryoSPARC as well as “heterogenous refinement” in cryoDRGN 36 with the particles from refinement on the main heterodimer body. This analysis, anchored on the SiRFP-FNR/SiRHP dimer, identified a dramatic movement of SiRFP’s Fld domain relative to the FNR domain, swinging 20 o between the most compact and most open forms (Video S1 and Fig. 4 A). The density for the linker joining the domains is not visible. In its most compact conformation, the Fld domain reaches the canonical “closed” conformation in which the Fld domain tucks into a cavity in the FNR domain (Fig. 4 B). In the most open conformation, the Fld domain assumes a different position to that seen in the “open” conformation determined by X-ray crystallography of the same monomeric variant and the average solution envelope determined from SANS, intermediate between the fully opened and closed conformations (Fig. 4 C) 19 , 20 . SiRFP-60X/SiRHP (124 kDa): The Fld and FNR domains are covalently linked in this dimer, thus restricting the swinging domain motion of the Fld domain seen in the SiRFP-60D/SiRHP dimer. No preferred orientation was identified either in 2D or 3D analysis (Fig. S3B). Nevertheless, there is no density for the linker between the Fld and FNR domains in the highest-resolution reconstruction, which suggests that the 36 amino acid linker is exceptionally dynamic but allows SiRFP’s two domains to come into close contact to catch them in a crosslink. Locking the Fld and FNR domains in the closed conformation resulted in the most stable variant for high-resolution analysis (2.84 Å), discussed above. With further classification of the SiRFP-60X/SiRHP structure, a sub-class of ~ 20,000 particles (10% of the particles used for highest resolution reconstruction for this dataset) emerged that revealed the linker between the two domains. The linker runs from the C-terminal end of the Fld domain (amino acid f-Ser207) to the N-terminal end of the FNR domain (amino acid f-Ile226) (Fig. 4 D). There is neither defined secondary structure nor contact with either of the SiRFP domains, thus the linker is not constrained in a single conformation across the entire ensemble of particles even when the domains are crosslinked. The highly mobile nature of these 36 amino acids explains how the Fld domain can reposition for contact with the oxidase partner, SiRHP, presumably to mediate electron transfer. Discussion The interface between the SiRHP and SiRFP subunits of the SiR holoenzyme is surprisingly minimal, driven by p-cation, hydrophobic, and planar stacking between aromatic amino acids. The non-canonical planar stacking interactions with strategic arginine residues are not as uncommon as would be expected from the canonical role that arginine plays in electrostatic interactions and are often found at important structural elements of the protein or enzyme 37 . In the case of the SiRFP/SiRHP dimer, these interactions either mediate the dimer interface itself or establish the structural platform for other amino acids to mediate the interface. Polar or electrostatic interactions play supporting roles to align the non-polar interactions because the arginine amino groups are free to make independent, more canonical, interactions with other side chains. The curious combination of non-ionic interactions that holds the complex together explains why cryo-EM analysis has been elusive – the hydrophobic air-water interface in a thin film for cryogenic plunging is destructive for such interactions 38 – 41 . We performed extensive trial-and-error probing grid type, grid substrate, hole size, method for hydrophilizing the grid, detergents, and plunging methods to identify what seems to be a singular condition that allowed cryogenic preservation for cryo-EM analysis: the combination of FF8 and the chameleon blot-free, rapid plunging. Neither FF8 alone with traditional plunging nor chameleon alone without FF8 was sufficient. Such low-throughput screening to find the one combination of variables that allows for high-quality sample preparation highlights the need for further technology development on the front-end of cryo-EM structure determination. High-throughput data collection with direct electron detectors is not sufficient to overcome cryo-EM sample preparation challenges. Surprisingly minimal structural elements determine the evolutionary pressure for siroheme-dependent SiR oxidases to partner with a diflavin reductase like SiRFP or a simpler Fd reductase. For example, in comparing the only known structures of monomeric assimilatory SiR oxidases with their reductase partners – E. coli SiRHP (this work) and Z. mays SIR 26 – a single loop extension between a core two-strand b-sheet discriminates between Fd binding in the case of Z. mays SIR and the binding of SiRHP to the FNR domain of SiRFP in the case of E. coli . E. coli has only one fdx gene that does not fulfill the reductase function for SiRHP, given that strains of E. coli lacking the gene encoding SiRFP ( cysJ ) cannot grow in a low sulfur environment 14 . The SiRFP/SiRHP dimeric structure reported here shows that when SiRHP binds SiRFP, the anion-binding active site loop is locked in its closed conformation, constrained from disorder. Nevertheless, h-Arg153 is in its phosphate-binding conformation. The position of the loop is, thus, impacted by the presence of SiRHP’s N-terminus that is absent in the X-ray crystal structure. The importance of this observation is understood in comparison to a series of X-ray crystal structures in various redox states and bound to various anions 31 , 32 . This series of structures shows how the active site re-orients with the changing state, namely that the loop between h-Asn149 to h-Arg153 is disordered in the oxidized state. Upon reduction the loop becomes ordered and h-Arg153 flips to bind the substrate. The question that remains is how SiRFP’s Fld domain approaches SiRHP’s metallic active site for productive electron transfer as it does not mediate the structural interaction between the subunits. Clearly there is dramatic flexibility within SiRFP, even when the Fld and FNR domains are covalently attached by an engineered crosslink – we cannot visualize the linker between the domains without extensive classification. Within this minimal dimer, the crosslink is intramolecular but in the SiR dodecameric holoenzyme, we have observed both intermolecular and intramolecular crosslinks 21 . Further, despite our extensive efforts at determining the structure of the full complex, the dodecamer has proven to be highly heterogeneous. Together, these observations support the hypothesis that within the holoenzyme, there is not a singular inter-subunit interaction between the Fld and FNR domains of any given SiRFP subunit within the octamer. The basis for SiRFP’s extreme flexibility is that the linker between the Fld and FNR domains is 36 amino acids long, compared to the 12 amino acid long linker in CPR 42 . Heterogeneity analysis of SiRFP-D in complex with SiRHP identified a dominant vector of motion parallel to the major axis of the FNR domain. This is in contrast to the highly extended conformation seen in the X-ray crystal structure of free SiRFP-D 20 . In this way, there is not a single axis along which the Fld domain moves relative to the FNR domain such that the Fld domain makes a productive interaction with SiRHP within a SiRHP/SiRFP dimer. Nevertheless, prior biochemical experimentation shows that cross-subunit interactions occur within the SiRFP octamer 21 . Thus, it is likely that a unique, single electron transfer-pathway does not exist within the full, dodecameric holoenzyme complex. Conclusions The combined technological developments in cryo-EM over the past decade all played a role in determining the structure of this elusive complex, from cryogenic sample preservation to data collection to image analysis. In doing so, we show that the subunits of NADPH-dependent assimilatory sulfite reductase bind through an interface governed by the N-terminus of SiRHP and the FNR domain of SiRFP. The interaction is surprisingly minimal, governed by a set of hydrophobic and stacking interactions. Structural pairing between a siroheme-dependent oxidase and diflavin reductase or ferredoxin appears to fall to a short loop between a conserved 2-stranded b sheet. The high mobility of SiRFP’s Fld domain relative to its FNR domain likely explains how a minimal dimer maintains redox-dependent functionality and provides a mechanism for redundant electron transfer pathways within the dodecamer that efficiently performs the six-electron reduction of sulfite to sulfide. Materials and methods Protein Expression and Purification Each SiRFP/SiRHP dimer variant was generated and purified as previously described 2 , 3 , 13 , 19 , 21 . Briefly, untagged SiRHP was co-purified with the following SiRFP variants: 1) SiRFP-43, a 43 kDa monomeric SiRFP variant lacking the N-terminal Fld domain and linker 19 ; 2) SiRFP-60D, a 60 kDa monomeric variant of SiRFP generated by truncating its first 52 amino acids with an additional internal truncation of six amino acid (D-AAPSQS) in the linker that joins the Fld and FNR domains 20 ; and 3) SiRFP-60X, the 60 kDa monomer with an engineered disulfide crosslink between the Fld and FNR domains, as has previously been analyzed in octameric SiRFP 21 . In SiRFP-60X, two background cysteines (C162T and C552S) were altered to avoid unwanted crosslinking and two cysteines were added (E121C and N556C) to form a disulfide bond between Fld and FNR domains of SiRFP, previously confirmed by mass spectroscopy analysis 21 . DNA sequencing confirmed the presence of all deletions and mutations in the variants. E. coli LMG194 cells (Invitrogen, Carlsbad, CA, USA) were transformed with the pBAD plasmid containing the genes encoding SiRFP-60D, SiRFP-60X, SiRFP-43 or SiRHP. All variants were expressed independently. Dimers were formed as followed: for SiRFP-43/SiRHP and SiRFP-60D/SiRHP, cells were mixed before lysis and the dimers co-purified. The SiRFP-60X/SiRHP dimer was assembled by reconstituting SiRFP-60X with 2 Eq SiRHP for 1 h on ice, as before 2 , 19 , 21 . A six-histidine tag was present in all SiRFP variants whereas SiRHP was untagged. Each recombinant E. coli strain was grown and induced at 25°C with 0.05% L-arabinose. SiRFP-43, SiRFP-60D and SiRHP cells were expressed for 4 hr whereas SiRFP-60X was expressed overnight. All variants were purified using a combination of Ni-NTA affinity chromatography (Cytiva, Marlborough, MA, USA), anion exchange HiTrap-Q HP chromatography (Cytiva, Marlborough, MA, USA) and Sephacryl S300-HR size exclusion chromatography (Cytiva, Marlborough, MA, USA) using previously optimized SPG buffers (17 mM succinic acid, 574 mM sodium dihydrogen phosphate, pH 6.8, 374 mM glycine, 200 mM NaCl) 2 , 13 , 19 , 21 . All variants have been extensively characterized with UV-Vis spectroscopy, size exclusion chromatography, SDS PAGE analysis for correct stoichiometry and with SANS for solution monodispersity (Fig. S7) 2 , 13 , 19 , 21 . Cryo-EM sample preparation SiRFP-43/SiRHP, SiRFP-60D/SiRHP and SiRFP-60X/SiRHP were all plunged using the chameleon blotless plunging system (SPT Labtech, Melbourn, UK) at 10 mg/mL protein. chameleon self-wicking grids 23 were backed with 18-gauge gold (Ted Pella, Redding CA, USA) on an Auto 306 vacuum coater (BOC Edwards, West Sussex, UK). The following glow discharge (GD)/wicking times (WTs) were used: SiRFP-43/SiRHP: 30 s GD, 130 ms WT; SiRFP-60D/SiRHP: 80 s GD, 195 ms WTs; SiRFP-60X/SiRHP 45 s GD, 175 ms WTs. All samples were premixed with FC-8 detergent 22 at 2 mM final concentration. The plunged grids were clipped and then screened for high quality with a Titan Krios operating the Leginon software package 43 . Data collection and processing SiRFP-43/SiRHP: 14,600 movies were collected at 300 KV on a Titan Krios using a GATAN K3 camera with 0.844 Å/pixel. After motion correction with Motioncor2 44 in the Relion GUI 45 , CTF estimation was performed using CTFFIND4 46 . Particles were picked by “blob picker” followed by “template picker” in CryoSPARC 34 . Initial 2D classification, followed by multiple rounds of 2D class selection/classification, identified 1,500,000 particles that were used for initial map building and non-uniform refinement in CryoSPARC. This process achieved a reported 3.6 Å-resolution map of the minimal dimer. However, the map features did not reflect the reported resolution. “Orientation diagnostic” job from CryoSPARC was performed to confirm this interpretation. To further confirm and diminish the preferred orientation issue, “Rebalance 2D” was performed to put particles into 7 super classes and limit the total particles in each superclass down to total of 350,000 and 105,000 particles. Non-uniform refinement 47 followed by orientation diagnostics performed in CryoSPARC produced a more accurate resolution (4.31 Å, 4.74 Å respectively) and better quality map. deepEMhancer 48 was used to sharpen these maps (Fig. S2). Particle picking performed by the TOPAZ algorithm 49 gave the same result. SiRFP-60D/SiRHP: 25,488 movies were collected at 300 KV on a Titan Krios using a GATAN K3 camera with 0.844 Å/pixel. Motion correction, CTF estimation, and particle picking were performed as for SiRFP-43/SiRHP. After multiple rounds of 2D classification, ~ 550,000 particles were used for initial map building and non-uniform refinement in CryoSPARC to achieve the final 3.49 Å-resolution structure of the dimer, masked to omit SiRFP’s Fld domain. Multiple refinements with different masking were performed. Masking to include the whole complex, including SiRFP’s Fld domain, resulted in a low resolution reconstruction for the Fld domain. Particles were then down sampled and imported to cryoDRGN to perform heterogenous refinement, giving a series of volumes showing the Fld domain movement. CryoSPARC 3D Flex 50 gave the same result as cryoDRGN. SiRFP-60X/SiRHP: ~10,000 movies were collected at 300 KV on a Titan Krios using a DE Apollo camera with 0.765 Å/pixel. Motion correction, CTF estimation, and particle picking were performed as above. Multiple rounds of 2D classifications identified ~ 185,000 particles that were used for initial map building and non-uniform refinement in CryoSPARC to achieve a 2.84 Å-resolution structure of the entire dimer, including SiRFP’s Fld domain. Further 3D variability was performed to track the Fld linker. After multiple rounds of 3D classification 10% of the particles (~ 20,000) were used to perform a non-uniform refinement, resulted in a 3.61 Å overall resolution structure, sharpened by deepEMhancer including the linker. Model building and refinement. Model building was initiated using the “fit in map” algorithm in Chimera 51 using the atomic model of SiRFP from PDB ID 6EFV 20 , fitting each of the Fld or FNR domains independently or the atomic model of SiRHP generated in AlphaFold 52 , 53 to capture its N-terminal 80 amino acids. Iterative real-space refinement in PHENIX 54 with manual fitting in Coot 55 yielded a model with a correlation of 0.85 (Table S3). The model deposited in PDB as 9C91and the map deposited in EMDB as EMD-45359. A stack from the final 3D reconstruction particles was deposited in the EMPIAR database under the accession number EMPIAR-12180. Declarations Acknowledgements and funding We kindly thank Dr.s Scott Stagg and Christopher Stroupe for helpful discussions. Some of this work was performed at the Simons Electron Microscopy Center at the New York Structural Biology Center, with major support from the Simons Foundation (SF349247). Florida State University supports cryo-EM in the Biological Imaging Resource Center, which houses the following equipment used in this study: a Gatan Solaris Plasma Cleaner (NIH grant S10 RR024564), a Hitachi HT7800 (NSF grant 2017869), a ThermoFisher Vitrobot Mark IV (NIH grant S10 RR024564), an SPI chameleon plunging system (NIH grant R24 GM145964), a ThermoFisher Titan Krios (NIH grant S10 RR025080), and a DE Apollo direct electron detector (NIH grant R35 GM139616). This work was further supported by National Science Foundation grants MCB1856502 and CHE1904612 to M.E.S. Author contributions BGE participated in data collection and performed the data analysis. NW and BGE prepared the protein specimen. KN and MA prepared cryogenic specimen and participated in data collection. IA optimized protein specimen preparation. AW and JHM supervised cryogenic specimen preparation and data collection. MES conceived of the study and oversaw experimentation. References Siegel LM, Davis PS (1974) Reduced nicotinamide adenine dinucleotide phosphate-sulfite reductase of enterobacteria. IV. The Escherichia coli hemoflavoprotein: subunit structure and dissociation into hemoprotein and flavoprotein components. J Biol Chem 249:1587–1598 Murray DT et al (2022) Neutron scattering maps the higher-order assembly of NADPH-dependent assimilatory sulfite reductase. Biophys J 121:1799–1812. https://doi.org/10.1016/j.bpj.2022.04.021 Askenasy I et al (2015) The N-terminal Domain of Escherichia coli Assimilatory NADPH-Sulfite Reductase Hemoprotein Is an Oligomerization Domain That Mediates Holoenzyme Assembly. J Biol Chem 290:19319–19333. https://doi.org/10.1074/jbc.M115.662379 Wang M et al (1997) Three-dimensional structure of NADPH-cytochrome P450 reductase: prototype for FMN- and FAD-containing enzymes. Proc Natl Acad Sci U S A 94:8411–8416 Girvan HM et al (2006) Flavocytochrome P450 BM3 and the origin of CYP102 fusion species. Biochem Soc Trans 34:1173–1177. https://doi.org/10.1042/BST0341173 Garcin ED et al (2004) Structural basis for isozyme-specific regulation of electron transfer in nitric-oxide synthase. J Biol Chem 279:37918–37927. https://doi.org/10.1074/jbc.M406204200 Zhang J et al (2001) Crystal structure of the FAD/NADPH-binding domain of rat neuronal nitric-oxide synthase. Comparisons with NADPH-cytochrome P450 oxidoreductase. J Biol Chem 276:37506–37513. https://doi.org/10.1074/jbc.M105503200 Olteanu H, Banerjee R (2001) Human methionine synthase reductase, a soluble P-450 reductase-like dual flavoprotein, is sufficient for NADPH-dependent methionine synthase activation. J Biol Chem 276:35558–35563. https://doi.org/10.1074/jbc.M103707200 Freeman SL et al (2018) Solution structure of the cytochrome P450 reductase-cytochrome. J Biol Chem 293:5210–5219. https://doi.org/10.1074/jbc.RA118.001941 Zhang L et al (2020) Structural insight into the electron transfer pathway of a self-sufficient P450 monooxygenase. Nat Commun 11:2676. https://doi.org/10.1038/s41467-020-16500-5 Zhang H et al (2018) The full-length cytochrome P450 enzyme CYP102A1 dimerizes at its reductase domains and has flexible heme domains for efficient catalysis. J Biol Chem 293:7727–7736. https://doi.org/10.1074/jbc.RA117.000600 Zeghouf M, Fontecave M, Coves J (2000) A simplifed functional version of the Escherichia coli sulfite reductase. J Biol Chem 275:37651–37656. https://doi.org/10.1074/jbc.M005619200 M005619200 [pii] Askenasy I et al (2018) Structure-Function Relationships in the Oligomeric NADPH-Dependent Assimilatory Sulfite Reductase. Biochemistry 57:3764–3772. https://doi.org/10.1021/acs.biochem.8b00446 Siegel LM, Murphy MJ, Kamin H (1973) Reduced nicotinamide adenine dinucleotide phosphate-sulfite reductase of enterobacteria. I. The Escherichia coli hemoflavoprotein: molecular parameters and prosthetic groups. J Biol Chem 248:251–264 Nakayama M, Akashi T, Hase T (2000) Plant sulfite reductase: molecular structure, catalytic function and interaction with ferredoxin. J Inorg Biochem 82:27–32 Schnell R, Sandalova T, Hellman U, Lindqvist Y, Schneider G (2005) Siroheme- and [Fe4-S4]-dependent NirA from Mycobacterium tuberculosis is a sulfite reductase with a covalent Cys-Tyr bond in the active site. J Biol Chem 280:27319–27328. https://doi.org/M502560200 Pinto R, Harrison JS, Hsu T, Jacobs WR Jr., Leyh TS (2007) Sulfite reduction in mycobacteria. J Bacteriol 189:6714–6722. https://doi.org/10.1128/jb.00487-07 Jespersen M, Pierik AJ, Wagner T (2023) Structures of the sulfite detoxifying F 420-dependent enzyme from Methanococcales. Nat Chem Biol 19:695–702. https://doi.org/10.1038/s41589-022-01232-y Murray DT, Weiss KL, Stanley CB, Nagy G, Stroupe ME (2021) Small-angle neutron scattering solution structures of NADPH-dependent sulfite reductase. J Struct Biol 213:107724. https://doi.org/10.1016/j.jsb.2021.107724 Tavolieri AM et al (2019) NADPH-dependent sulfite reductase flavoprotein adopts an extended conformation unique to this diflavin reductase. J Struct Biol. https://doi.org/10.1016/j.jsb.2019.01.001 Walia N et al (2023) Domain crossover in the reductase subunit of NADPH-dependent assimilatory sulfite reductase. J Struct Biol 215:108028. https://doi.org/10.1016/j.jsb.2023.108028 Kampjut D, Steiner J, Sazanov LA (2021) Cryo-EM grid optimization for membrane proteins. iScience 24:102139. https://doi.org/10.1016/j.isci.2021.102139 Levitz TS et al (2022) Approaches to Using the Chameleon: Robust, Automated, Fast-Plunge cryoEM Specimen Preparation. Front Mol Biosci 9:903148. https://doi.org/10.3389/fmolb.2022.903148 Crane BR, Siegel LM, Getzoff ED (1995) Sulfite reductase structure at 1.6 A: evolution and catalysis for reduction of inorganic anions. Science 270:59–67 Crane BR, Getzoff ED (1996) The relationship between structure and function for the sulfite reductases. Curr Opin Struct Biol 6:744–756. https://doi.org/S0959-440X (96)80003-0 [pii] Kim JY, Nakayama M, Toyota H, Kurisu G, Hase T (2016) Structural and mutational studies of an electron transfer complex of maize sulfite reductase and ferredoxin. J Biochem. https://doi.org/10.1093/jb/mvw016 Oliveira TF et al (2008) The crystal structure of Desulfovibrio vulgaris dissimilatory sulfite reductase bound to DsrC provides novel insights into the mechanism of sulfate respiration. J Biol Chem 283:34141–34149. https://doi.org/M805643200 Oliveira TF et al (2011) Structural insights into dissimilatory sulfite reductases: structure of desulforubidin from desulfomicrobium norvegicum. Front Microbiol 2:71. https://doi.org/10.3389/fmicb.2011.00071 Parey K, Warkentin E, Kroneck PM, Ermler U (2010) Reaction cycle of the dissimilatory sulfite reductase from Archaeoglobus fulgidus. Biochemistry 49:8912–8921. https://doi.org/10.1021/bi100781f Schiffer A et al (2008) Structure of the dissimilatory sulfite reductase from the hyperthermophilic archaeon Archaeoglobus fulgidus. J Mol Biol 379, 1063–1074 https://doi.org/S0022-2836(08)00456-7 [pii] 10.1016/j.jmb.2008.04.027 Crane BR, Siegel LM, Getzoff ED (1997) Probing the catalytic mechanism of sulfite reductase by X-ray crystallography: structures of the Escherichia coli hemoprotein in complex with substrates, inhibitors, intermediates, and products. Biochemistry 36:12120–12137. https://doi.org/10.1021/bi971066i bi971066i [pii] Crane BR, Siegel LM, Getzoff ED (1997) Structures of the siroheme- and Fe4S4-containing active center of sulfite reductase in different states of oxidation: heme activation via reduction-gated exogenous ligand exchange. Biochemistry 36:12101–12119. https://doi.org/10.1021/bi971065q bi971065q [pii] Smith KW, Stroupe ME (2012) Mutational analysis of sulfite reductase hemoprotein reveals the mechanism for coordinated electron and proton transfer. Biochemistry 51:9857–9868. https://doi.org/10.1021/bi300947a Punjani A, Rubinstein JL, Fleet DJ, Brubaker MA (2017) cryoSPARC: algorithms for rapid unsupervised cryo-EM structure determination. Nat Methods 14:290–296. https://doi.org/10.1038/nmeth.4169 Aiyer S, Zhang C, Baldwin PR, Lyumkis D (2021) Evaluating Local and Directional Resolution of Cryo-EM Density Maps. Methods Mol Biol 2215:161–187. https://doi.org/10.1007/978-1-0716-0966-8_8 Zhong ED, Bepler T, Berger B, Davis JH (2021) CryoDRGN: reconstruction of heterogeneous cryo-EM structures using neural networks. Nat Methods 18:176–185. https://doi.org/10.1038/s41592-020-01049-4 Flocco MM, Mowbray SL (1994) Planar stacking interactions of arginine and aromatic side-chains in proteins. J Mol Biol 235:709–717. https://doi.org/10.1006/jmbi.1994.1022 D'Imprima E et al (2019) Protein denaturation at the air-water interface and how to prevent it. Elife 8. https://doi.org/10.7554/eLife.42747 Glaeser RM et al (2016) Factors that Influence the Formation and Stability of Thin, Cryo-EM Specimens. Biophys J 110:749–755. https://doi.org/10.1016/j.bpj.2015.07.050 Glaeser RM, Han BG (2017) Opinion: hazards faced by macromolecules when confined to thin aqueous films. Biophys Rep 3:1–7. https://doi.org/10.1007/s41048-016-0026-3 Glaeser RM, PROTEINS, INTERFACES, AND, CRYO-EM GRIDS (2018) Curr Opin Colloid Interface Sci 34, 1–8 https://doi.org/10.1016/j.cocis.2017.12.009 Xia C et al (2011) Conformational changes of NADPH-cytochrome P450 oxidoreductase are essential for catalysis and cofactor binding. J Biol Chem 286:16246–16260. https://doi.org/10.1074/jbc.M111.230532 Suloway C et al (2005) Automated molecular microscopy: the new Leginon system. J Struct Biol 151:41–60. https://doi.org/S1047-8477( 05)00072-9 [pii] 10.1016/j.jsb.2005.03.010 Zheng SQ et al (2017) MotionCor2: anisotropic correction of beam-induced motion for improved cryo-electron microscopy. Nat Methods 14:331–332. https://doi.org/10.1038/nmeth.4193 Scheres SH (2012) RELION: implementation of a Bayesian approach to cryo-EM structure determination. J Struct Biol 180:519–530. https://doi.org/10.1016/j.jsb.2012.09.006 Rohou A, Grigorieff N (2015) CTFFIND4: Fast and accurate defocus estimation from electron micrographs. J Struct Biol 192:216–221 Punjani A, Zhang H, Fleet DJ (2020) Non-uniform refinement: adaptive regularization improves single-particle cryo-EM reconstruction. Nat Methods 17:1214–1221. https://doi.org/10.1038/s41592-020-00990-8 Sanchez-Garcia R et al (2021) DeepEMhancer: a deep learning solution for cryo-EM volume post-processing. Commun Biol 4:874. https://doi.org/10.1038/s42003-021-02399-1 Bepler T et al (2019) Positive-unlabeled convolutional neural networks for particle picking in cryo-electron micrographs. Nat Methods 16:1153–1160. https://doi.org/10.1038/s41592-019-0575-8 Punjani A, Fleet D (2023) 3D Flexible Refinement: Determining Structure and Motion of Flexible Proteins from Cryo-EM. Microsc Microanal 29:1024. https://doi.org/10.1093/micmic/ozad067.518 Pettersen EF et al (2004) UCSF Chimera - A visualization system for exploratory research and analysis. J Comput Chem 25:1605–1612. https://doi.org/10.1002/jcc.20084 Jumper J et al (2021) Highly accurate protein structure prediction with AlphaFold. Nature 596:583–589. https://doi.org/10.1038/s41586-021-03819-2 Roney JP, Ovchinnikov S (2022) State-of-the-Art Estimation of Protein Model Accuracy Using AlphaFold. Phys Rev Lett 129:238101. https://doi.org/10.1103/PhysRevLett.129.238101 Adams PD et al (2010) PHENIX: a comprehensive Python-based system for macromolecular structure solution. Acta Crystallogr D Biol Crystallogr 66:213–221. https://doi.org/S0907444909052925 [pii] 10.1107/S0907444909052925 Emsley P, Lohkamp B, Scott WG, Cowtan K (2010) Features and development of Coot. Acta Crystallogr D Biol Crystallogr 66:486–501. https://doi.org/S0907444910007493 [pii] 10.1107/S0907444910007493 Additional Declarations There is NO Competing Interest. Supplementary Files VideoS1.mp4 Supplemental Video 1 supplementalNC.docx Cite Share Download PDF Status: Published Journal Publication published 26 Mar, 2025 Read the published version in Nature Communications → Version 1 posted You are reading this latest preprint version Research Square lets you share your work early, gain feedback from the community, and start making changes to your manuscript prior to peer review in a journal. As a division of Research Square Company, we’re committed to making research communication faster, fairer, and more useful. We do this by developing innovative software and high quality services for the global research community. Our growing team is made up of researchers and industry professionals working together to solve the most critical problems facing scientific publishing. Also discoverable on Platform About Our Team In Review Editorial Policies Advisory Board Help Center Resources Author Services Accessibility API Access RSS feed Manage Cookie Preferences © Research Square 2026 | ISSN 2693-5015 (online) Privacy Policy Terms of Service Do Not Sell My Personal Information {"props":{"pageProps":{"initialData":{"identity":"rs-4758050","acceptedTermsAndConditions":true,"allowDirectSubmit":false,"archivedVersions":[],"articleType":"Article","associatedPublications":[],"authors":[{"id":332140072,"identity":"b6851b1a-e1c2-41b7-a9ae-8ae1c3ff0dc5","order_by":0,"name":"M. Elizabeth Stroupe","email":"data:image/png;base64,iVBORw0KGgoAAAANSUhEUgAAAZAAAAAyAQMAAABI0h/eAAAABlBMVEX///8AAABVwtN+AAAACXBIWXMAAA7EAAAOxAGVKw4bAAAA50lEQVRIiWNgGAWjYFAD9h6GAwwMcmC2BAG1jA1giucMSIsxKVokckAkEVp029ufP66oOJzHL/n24IEPDAb2BgeYD97mwaPF7MwZw8YzZw4XS87OSzg4g8EgccMBtmRrvFpu5DA2NrYdTtxwO8fgMA/DnwSDAzxm0ni13H/+EKxl/80zIC0gh/F/w6/lBoMhxBYJHrAWxg0HeNjwazmTYziz4Ux6scQZkF8MDBJnHmYztpyDT8vx4w8+NlRY5/G3nz384UOFgT3f8eaHN97g0QIDCRDKAIiZiVCOpGUUjIJRMApGARYAACXWUlZGzZTgAAAAAElFTkSuQmCC","orcid":"https://orcid.org/0000-0002-4393-1951","institution":"Florida State University","correspondingAuthor":true,"prefix":"","firstName":"M.","middleName":"Elizabeth","lastName":"Stroupe","suffix":""},{"id":332140073,"identity":"29c3fb5c-9d71-4e91-ba2a-3161d60fb740","order_by":1,"name":"Behrouz Ghazi Esfahani","email":"","orcid":"","institution":"Florida State University","correspondingAuthor":false,"prefix":"","firstName":"Behrouz","middleName":"Ghazi","lastName":"Esfahani","suffix":""},{"id":332140074,"identity":"e7927fc8-19af-4be2-a4bd-6a2dfccf1c6f","order_by":2,"name":"Nidhi Walia","email":"","orcid":"","institution":"Florida State University","correspondingAuthor":false,"prefix":"","firstName":"Nidhi","middleName":"","lastName":"Walia","suffix":""},{"id":332140075,"identity":"f2ffc63d-c053-4b60-aace-75726cb87b00","order_by":3,"name":"Kasahun Neselu","email":"","orcid":"","institution":"New York Structural Biology Center","correspondingAuthor":false,"prefix":"","firstName":"Kasahun","middleName":"","lastName":"Neselu","suffix":""},{"id":332140076,"identity":"7fc656cf-8e3b-432a-ac28-de762984c731","order_by":4,"name":"Mahira Aragon","email":"","orcid":"","institution":"New York Structural Biology Center","correspondingAuthor":false,"prefix":"","firstName":"Mahira","middleName":"","lastName":"Aragon","suffix":""},{"id":332140077,"identity":"71125be0-345d-4ae2-bd80-189dd3707dc0","order_by":5,"name":"Isabel Askenasy","email":"","orcid":"https://orcid.org/0000-0002-2775-2590","institution":"University of Cambridge","correspondingAuthor":false,"prefix":"","firstName":"Isabel","middleName":"","lastName":"Askenasy","suffix":""},{"id":332140078,"identity":"2ce5d0ed-bd35-4b7a-b30a-1953b9d95d8a","order_by":6,"name":"Hui Wei","email":"","orcid":"","institution":"Rutgers University","correspondingAuthor":false,"prefix":"","firstName":"Hui","middleName":"","lastName":"Wei","suffix":""},{"id":332140079,"identity":"d0ebdbb6-ab3c-4cc6-aaef-7e1cbd275ed8","order_by":7,"name":"Joshua Mendez","email":"","orcid":"https://orcid.org/0000-0002-3317-5654","institution":"New York Structural Biology Center","correspondingAuthor":false,"prefix":"","firstName":"Joshua","middleName":"","lastName":"Mendez","suffix":""}],"badges":[],"createdAt":"2024-07-17 17:50:23","currentVersionCode":1,"declarations":"","doi":"10.21203/rs.3.rs-4758050/v1","doiUrl":"https://doi.org/10.21203/rs.3.rs-4758050/v1","draftVersion":[],"editorialEvents":[{"content":"https://doi.org/10.1038/s41467-025-58037-5","type":"published","date":"2025-03-26T04:00:00+00:00"}],"editorialNote":"","failedWorkflow":false,"files":[{"id":61284928,"identity":"de031f19-26b0-42f7-bf1d-7357695fc2db","added_by":"auto","created_at":"2024-07-29 06:12:15","extension":"png","order_by":1,"title":"Figure 1","display":"","copyAsset":false,"role":"figure","size":848787,"visible":true,"origin":"","legend":"\u003cp\u003eSiRFP and SiRHP interact through SiRFP’s FNR domain and SiRHP’s N-terminus. \u003cstrong\u003eA)\u003c/strong\u003e The SiRFP/SiRHP interface is minimal, governed by amino acids from three loops from SiRHP that fit into pockets on SiRFP’s FNR domain. SiRFP’s Fld domain is pink, SiRFP’s FNR domain is light blue. SiRHP’s N-terminal 80 amino acids are dark blue. SiRHP’s core two S/NiRRs are green. \u003cstrong\u003eB)\u003c/strong\u003e The interface between the SiR subunits is governed primarily by hydrophobic interactions, in teal. Charged amino acids are gold and neutral amino acids are white. The subunits are displayed such that they are rotated 90\u003csup\u003eo\u003c/sup\u003e away from one another to show their interacting faces. \u003cstrong\u003eC)\u003c/strong\u003e The topology diagram of SiRHP, highlighting the N-terminal 80 amino acids whose structure was previously unknown but that govern the interactions with SiRFP. The amino acids that interact with SiRFP are localized to three loops. Hydrophobic amino acids are green, cationic amino acids are purple, and polar amino acids are yellow. The first S/NiRR is dark green whereas the second S/NiRR is light green, with the connecting h11 and linker in yellow and cyan, respectively (“Created with\u0026nbsp;BioRender.com.”). \u003cstrong\u003eD)\u003c/strong\u003e The interface between SiRFP and SiRHP is minimal. The SiRFP-FNR is light blue with light purple highlighted residues that are important for forming the SiRHP binding site. SiRHP’S N-terminal 80 amino acids are blue and the S/NiRR repeats are green with mauve highlighted residues that mediate the interaction with SiRFP. \u003cstrong\u003eE-G)\u003c/strong\u003e Select densities for important side chain interactions that mediate the interface between the SiRFP and SiRHP interface, colored as in \u003cstrong\u003eD\u003c/strong\u003e.\u003c/p\u003e","description":"","filename":"image1.png","url":"https://assets-eu.researchsquare.com/files/rs-4758050/v1/e42d204067886178622e21ba.png"},{"id":61284927,"identity":"c8f826c5-4fde-4e77-9446-bdb58da70c14","added_by":"auto","created_at":"2024-07-29 06:12:15","extension":"png","order_by":2,"title":"Figure 2","display":"","copyAsset":false,"role":"figure","size":619974,"visible":true,"origin":"","legend":"\u003cp\u003e\u003cem\u003eZea mays\u003c/em\u003e SIR and \u003cem\u003eEscherichia coli\u003c/em\u003e SiRHP interact with their reductase partners in different ways. \u003cstrong\u003eA)\u003c/strong\u003e SIR uses a ferredoxin (Fd, pink) as its reductase partner, which binds such that its Fe\u003csub\u003e2\u003c/sub\u003eS\u003csub\u003e2\u003c/sub\u003e cluster is close to SIR’s Fe\u003csub\u003e4\u003c/sub\u003eS\u003csub\u003e4\u003c/sub\u003e cluster, positioned between SIR’s N-terminus (yellow) and domain 2 (green). (PDB 5H92\u003csup\u003e26\u003c/sup\u003e) \u003cstrong\u003eB)\u003c/strong\u003e SiRHP h4 (mauve) is the sole secondary structural feature that discriminates between SiRFP binding and Fd binding. SiRHP’s N-terminal 80 amino acids are dark blue. SIR’s N-terminal amino acids are yellow with the element equivalent to SiRHP h4 in medium blue. \u003cstrong\u003eC)\u003c/strong\u003e SiRFP (light blue) and Fd (pink) bind on different faces of SiRHP (dark blue) or SIR (yellow). Only the core S/NiRRs from SiRHP (green) are shown for clarity. Other elements are colored as in \u003cstrong\u003eB\u003c/strong\u003e. \u003cstrong\u003eD)\u003c/strong\u003e SIR R324 interacts with Fd Y37. That interaction is prevented in SiRHP because the position equivalent to SIR R324 is h-G262 (mauve).\u003c/p\u003e","description":"","filename":"image2.png","url":"https://assets-eu.researchsquare.com/files/rs-4758050/v1/76ac4b86bd8dfa2a3890e6ce.png"},{"id":61284930,"identity":"49105a73-97e4-4c94-b961-0a96b316fda4","added_by":"auto","created_at":"2024-07-29 06:12:16","extension":"png","order_by":3,"title":"Figure 3","display":"","copyAsset":false,"role":"figure","size":851017,"visible":true,"origin":"","legend":"\u003cp\u003eSiRHP’s active site loop that forms the complete substrate binding site is in its closed conformation with the intact N-terminus, bound to SiRFP. Panels are labeled by their binding partner and oxidation state. \u003cstrong\u003eA)\u003c/strong\u003e When SiRFP binds oxidized SiRHP with the bound, inhibitory phosphate, the anion-binding loop containing h-N149 is ordered but with h-R153 in a similar orientation as in free SiRHP. Amino acids important for anion binding or loop position are dark pink. The siroheme methyl group from ring A is a dark green sphere and the interaction between h-Q60 and the ring B propionyl group is marked by a gray dash. \u003cstrong\u003eB)\u003c/strong\u003e In the absence of its N-terminal 80 amino acids, SiRFP, or the sulfite substrate, SiRHP’s anion binding loop is disordered (pink ribbons). The amino acids that are important for anion binding in response to oxidation state or oxidase binding are dark purple (PDB 1AOP\u003csup\u003e24\u003c/sup\u003e). Upon reduction and with bound sulfite, the active site loop becomes ordered and h-R153 flips over to bind the substrate (light gray ribbons). The amino acids that are important for anion binding in response to oxidation state or oxidase binding are light purple (PDB 2GEP\u003csup\u003e31,32\u003c/sup\u003e).\u003cstrong\u003e C)\u003c/strong\u003e Superimposition of the SiRFP-bound SiRHP (green) and the reduced, sulfite-bound, SiRHP (light gray) shows the anion-binding loop to be in similar position as in reduced SiRHP, primed for sulfite binding. \u003cstrong\u003eD) \u003c/strong\u003eSuperimposition\u003cstrong\u003e \u003c/strong\u003eof the three SiRHP structures bound to different anions or binding partner highlights the intermediate conformation induced by SiRFP binding. This intermediate conformation shows ordering of SiRHP’s N-terminus and the active site loop in a phosphate-bound, oxidized state. Models are colored as in \u003cstrong\u003eA-C\u003c/strong\u003e.\u003c/p\u003e","description":"","filename":"image3.png","url":"https://assets-eu.researchsquare.com/files/rs-4758050/v1/d8f74892ebb01b9dcada4425.png"},{"id":61284931,"identity":"5658a7ac-0ddb-493c-ac75-48dbc399ebd2","added_by":"auto","created_at":"2024-07-29 06:12:16","extension":"png","order_by":4,"title":"Figure 4","display":"","copyAsset":false,"role":"figure","size":677075,"visible":true,"origin":"","legend":"\u003cp\u003eSiRFP’s Fld domain is highly mobile relative to its FNR domain. \u003cstrong\u003eA)\u003c/strong\u003e One of the latent variable dimensions of motion demonstrates an opening of the Fld domain, even in SiRFP-60D, with the motion parallel to the midline of the FNR domain (light blue). The most-closed position is dark pink, the most open is light pink. SiRHP is green, with its N-terminus in dark blue. The angle is measured between the lines connecting the last traceable amino acid in the FNR domain with the open and closed Fld volume center of mass. \u003cstrong\u003eB)\u003c/strong\u003e The map of the most closed conformation of the Fld domain (dark pink) corresponds to the position of the Fld domain in the dimer formed with SiRFP-60X, far from where it is in the extended conformation of SiRFP-60D in crystals (teal, PDB 6EFV,\u003csup\u003e20\u003c/sup\u003e). \u003cstrong\u003eC)\u003c/strong\u003e The map of the most opened conformation (light pink) shows that the Fld domain lands in an intermediary position between a highly extended position seen in the X-ray crystal structure (teal, PDB 6EFV\u003csup\u003e20\u003c/sup\u003e) and the closed conformation seen in \u003cstrong\u003eB\u003c/strong\u003e. \u003cstrong\u003eD)\u003c/strong\u003e Even with the Fld and FNR domains constrained by a crosslink, the 36 amino acid-long linker between the domains is largely disordered, visible only in a small (~20,000) subset of particles (purple volume). \u0026nbsp;\u0026nbsp;\u003c/p\u003e","description":"","filename":"image4.png","url":"https://assets-eu.researchsquare.com/files/rs-4758050/v1/ca33391d10c97c6543cace4f.png"},{"id":79332903,"identity":"ef754a60-d024-4c28-839c-65a9b0ed6193","added_by":"auto","created_at":"2025-03-27 07:05:40","extension":"pdf","order_by":0,"title":"","display":"","copyAsset":false,"role":"manuscript-pdf","size":4071744,"visible":true,"origin":"","legend":"","description":"","filename":"manuscript.pdf","url":"https://assets-eu.researchsquare.com/files/rs-4758050/v1/e16d13c5-52ad-4300-bfc8-d2e0e2d2b9a3.pdf"},{"id":61285447,"identity":"922bc906-ff67-4d02-9841-a758253968f0","added_by":"auto","created_at":"2024-07-29 06:20:16","extension":"mp4","order_by":1,"title":"","display":"","copyAsset":false,"role":"supplement","size":373654,"visible":true,"origin":"","legend":"Supplemental Video 1","description":"","filename":"VideoS1.mp4","url":"https://assets-eu.researchsquare.com/files/rs-4758050/v1/26aefb6e190491a43e85a999.mp4"},{"id":61284932,"identity":"75ba612c-a847-40c3-9d9e-aa80a8461649","added_by":"auto","created_at":"2024-07-29 06:12:16","extension":"docx","order_by":2,"title":"","display":"","copyAsset":false,"role":"supplement","size":24523425,"visible":true,"origin":"","legend":"","description":"","filename":"supplementalNC.docx","url":"https://assets-eu.researchsquare.com/files/rs-4758050/v1/c39b52b0141fdda895ebf999.docx"}],"financialInterests":"There is \u003cb\u003eNO\u003c/b\u003e Competing Interest.","formattedTitle":"Structure of dimerized assimilatory NADPH-dependent sulfite reductase reveals the minimal interface for diflavin reductase binding","fulltext":[{"header":"Introduction","content":"\u003cp\u003eAssimilatory sulfite reduction by NADPH-dependent sulfite reductase (SiR) is essential to produce sulfide for incorporation into sulfur-containing biomolecules. In g-proteobacteria like \u003cem\u003eEscherichia coli\u003c/em\u003e, SiR is a multimeric oxidoreductase composed of an octameric diflavin reductase (SiRFP) and four independently binding subunits of a siroheme-containing hemoprotein (SiRHP)\u003csup\u003e1-3\u003c/sup\u003e.\u003c/p\u003e\n\u003cp\u003eThe \u003cem\u003eE. coli\u003c/em\u003e SiRFP subunit is homologous to cytochrome P450 (CYP) reductase (CPR)\u003csup\u003e4\u003c/sup\u003e, the reductase domain of the bacterial CYP/CPR fusion CYP102A1/CYPBM3\u003csup\u003e5\u003c/sup\u003e, the reductase domain of nitric oxide synthase (NOSr)\u003csup\u003e6,7\u003c/sup\u003e, and methionine synthase reductase (MSR)\u003csup\u003e8\u003c/sup\u003e. One of the hallmarks of this diflavin reductase family is that they are exceptionally conformationally malleable, which makes structural analysis challenging. For example, to date there are no high-resolution structures of the full-length NOS homodimer, the complex between methionine synthase and MSR, or the SiR heterododecameric holoenzyme (for simplicity, here referred to as a dodecamer). The structure of the CYP/CPR heterodimer and CYPBM3 are known\u003csup\u003e9-11\u003c/sup\u003e, but CYP/CPR form a 1:1 heterodimer, whereas the oxidase and reductase domains are fused in CYPBM3, so little can be inferred about other homologs that function as higher-order protein complexes.\u003c/p\u003e\n\u003cp\u003eDespite its homology to other well-studied diflavin reductases, SiRFP is unique because it assembles into an octamer through its N-terminal 52 amino acids\u003csup\u003e12\u003c/sup\u003e. Removing those amino acids results in a 60 kDa monomer (SiRFP-60), which binds SiRHP as a 1:1 heterodimer with reduced activity\u003csup\u003e3\u003c/sup\u003e. Further, removing the complete N-terminal FMN-binding flavodoxin (Fld) domain results in a 43 kDa monomer that contains just the NADPH- and FAD-binding NADP\u003csup\u003e+\u003c/sup\u003e ferredoxin reductase (FNR) domain (SiRFP-43), which also binds SiRHP as a 1:1 heterodimer but is inactive for electron transfer\u003csup\u003e13\u003c/sup\u003e (Abbreviations and theoretical molecular weights are summarized in Table S1).\u003c/p\u003e\n\u003cp\u003eSiRHP has few known homologs because of its unique siroheme/Fe\u003csub\u003e4\u003c/sub\u003eS\u003csub\u003e4\u003c/sub\u003e cluster assembly that form the sulfite-binding active site\u003csup\u003e14\u003c/sup\u003e. Assimilatory SiRs or siroheme-dependent nitrite reductases (NiRs) from other bacterial species or plants have a similar hemoprotein but use a transiently bound ferredoxin as their electron source\u003csup\u003e15-17\u003c/sup\u003e. SiRs that are responsible for energy conversion, dissimilatory sulfite reductases (DSRs), share a common siroheme binding fold but are heterotetrameric and are often fused to auxiliary domains\u003csup\u003e18\u003c/sup\u003e. Their electron donors are poorly understood. Here, we show the first high-resolution cryogenic electron microscopy (cryo-EM) structure of the minimal SiRFP/SiRHP dimer, which elucidates their binding interface to understand how SiRHP tightly binds SiRFP\u0026rsquo;s FNR domain.\u003c/p\u003e"},{"header":"Results","content":"\u003cdiv id=\"Sec2\" class=\"Section2\"\u003e \u003ch2\u003eThe SiRHP-SiRFP interaction is highly sensitive to cryogenic TEM (cryo-EM) preparation\u003c/h2\u003e \u003cp\u003eWe determined the structure of the SiRFP/SiRHP dimer from three modified, minimal dimers, each of which is named by the change to SiRFP and its resulting molecular weight (Table \u003cspan refid=\"MOESM1\" class=\"InternalRef\"\u003eS1\u003c/span\u003e). First, we truncated SiRFP to remove the N-terminal Fld domain (SiRFP-43/SiRHP). This is an inactive dimer, as the Fld domain is required for electron transfer. Nevertheless, the two subunits bind tightly and is the most simplified complex between SiRFP and SiRHP\u003csup\u003e\u003cspan citationid=\"CR19\" class=\"CitationRef\"\u003e19\u003c/span\u003e\u003c/sup\u003e. Second, we truncated both SiRFP\u0026rsquo;s N-terminal octamerization domain as well as the linker between the Fld and FNR domains to create a monomeric SiRFP that can be locked in an open position (SiRFP-60D/SiRHP)\u003csup\u003e20\u003c/sup\u003e. Third, we generated a variant of monomeric SiRFP-60 lacking reactive cysteines into which we engineered a disulfide bond between the Fld and FNR domains (SiRFP-60X/SiRHP)\u003csup\u003e21\u003c/sup\u003e.\u003c/p\u003e \u003cp\u003eEach variant is highly sensitive to traditional blotting plunge-freezing methods for cryo-EM preservation. To overcome this sensitivity, we combined the protection of a high critical micelle concentration, fluorinated lipid, fos-choline-8 (FF8, Creative Biolabs, Shirley, NY, USA )\u003csup\u003e\u003cspan citationid=\"CR22\" class=\"CitationRef\"\u003e22\u003c/span\u003e\u003c/sup\u003e, with the blot-free, rapid plunging afforded by the chameleon system (SPT Labtech, Melbourn, UK)\u003csup\u003e\u003cspan citationid=\"CR23\" class=\"CitationRef\"\u003e23\u003c/span\u003e\u003c/sup\u003e. This cryo-EM sample preparation helped us to retain each intact complex within near ideal ice thickness and avoid denaturation at the air water interface (Fig. \u003cspan refid=\"MOESM1\" class=\"InternalRef\"\u003eS1\u003c/span\u003e). The smallest complex (SiRFP-43/SiRHP) showed well-aligned 2D class averages, however the 3D structure revealed structural anisotropy, either due to its small size/asymmetric geometry or from a preferred orientation, that limited high-resolution analysis despite the absence of mobile elements (Figs. S1A and S2). SiRFP-60D/SiRHP showed moderate-resolution density (3.54 \u0026Aring;) for the SiRFP FNR domain and SiRHP, however the N-terminal Fld domain was not visible (Figs. S1B and S3A). The 2.84 \u0026Aring;-resolution structure of SiRFP-60X/SiRHP revealed the most detail for SiRFP\u0026rsquo;s Fld and FNR domains, despite a lack of density for the linker between SiRFP\u0026rsquo;s domains in the highest-resolution reconstruction (Figs.\u0026nbsp;\u003cspan refid=\"Fig1\" class=\"InternalRef\"\u003e1\u003c/span\u003eA, S1C, and S3B). High resolution features for each of the cofactors in both subunits supported this reported resolution (Fig. S4). Therefore, we analyzed the SiRFP-SiRHP interface for this construct in detail.\u003c/p\u003e \u003cp\u003e \u003c/p\u003e \u003c/div\u003e\n\u003ch3\u003eSiRFP-SiRHP binding\u003c/h3\u003e\n\u003cp\u003eThe SiR dodecameric holoenzyme is composed of oligomers of the dimers discussed here and is about 800 kDa in mass. Despite this large mass, the binding interface between the minimal SiRFP/SiRHP dimer is small, reaching the surface area of 1,138 \u0026Aring;\u003csup\u003e\u003cspan citationid=\"CR2\" class=\"CitationRef\"\u003e2\u003c/span\u003e\u003c/sup\u003e relative to the overall surface of 43,610 \u0026Aring;\u003csup\u003e\u003cspan citationid=\"CR2\" class=\"CitationRef\"\u003e2\u003c/span\u003e\u003c/sup\u003e. SiRHP alone has a solvent-exposed surface of 25,680 \u0026Aring;\u003csup\u003e\u003cspan citationid=\"CR2\" class=\"CitationRef\"\u003e2\u003c/span\u003e\u003c/sup\u003e. SiRFP-60X alone has a solvent-exposed surface of 21,930 \u0026Aring;\u003csup\u003e\u003cspan citationid=\"CR2\" class=\"CitationRef\"\u003e2\u003c/span\u003e\u003c/sup\u003e. That is, for a large complex only about 2.6% of the solvent-exposed surface is buried upon subunit binding. This is consistent with hydrogen-deuterium exchange data on the complex that reveals single, short peptides from each subunit that become occluded upon binding\u003csup\u003e\u003cspan citationid=\"CR3\" class=\"CitationRef\"\u003e3\u003c/span\u003e\u003c/sup\u003e, the sequences of which predicted the interface would be dominated by hydrophobic interactions (Fig.\u0026nbsp;\u003cspan refid=\"Fig1\" class=\"InternalRef\"\u003e1\u003c/span\u003eB).\u003c/p\u003e \u003cp\u003eThe interface is governed by the N-terminus of SiRHP. The structure of this region is previously uncharacterized as it is proteolytically removed in the X-ray crystal structure of \u003cem\u003eE. coli\u003c/em\u003e SiRHP\u003csup\u003e\u003cspan citationid=\"CR24\" class=\"CitationRef\"\u003e24\u003c/span\u003e\u003c/sup\u003e. These 80 amino acids follow the topology helix 1 - loop - helix 2 - turn - helix 3 - b-strand 1 - helix 4 - loop - b-strand 2 (Fig.\u0026nbsp;\u003cspan refid=\"Fig1\" class=\"InternalRef\"\u003e1\u003c/span\u003eC). Only amino acids from the turn, helix 4, and surrounding loops directly interact with SiRFP. The regions that are N-terminal to the interface interact with domain 1 or the N-terminal half of the parachute domain (\u003cem\u003ei.e.\u003c/em\u003e the first sulfite or nitrite reductase repeat (S/NiRR)\u003csup\u003e24,25\u003c/sup\u003e), breaking SiRHP\u0026rsquo;s pseudo two-fold symmetry (Figs.\u0026nbsp;\u003cspan refid=\"Fig1\" class=\"InternalRef\"\u003e1\u003c/span\u003eA and S5A). An extension to the parachute domain that is not ordered in the original crystal structure, from amino acids 184\u0026ndash;209, helps to hold the N-terminus in place (Fig. S5A). Those that are N-terminal to the interface approach the distal active site, but do not contribute significantly to anion or siroheme binding as they are held back by interactions to SiRFP, discussed below. Moving N-terminally along the peptide, it then turns back to form the loop that binds a pocket in SiRFP before moving away from SiRFP. The N-terminal most amino acids reach all the way to the other side of SiRHP, interacting with the N-terminus of the a-helix (h11) that precedes the linker that joins the two S/NiRRs and mimics the siroheme binding site (Fig. S5B)\u003csup\u003e\u003cspan citationid=\"CR24\" class=\"CitationRef\"\u003e24\u003c/span\u003e\u003c/sup\u003e.\u003c/p\u003e \u003cp\u003eThe central interaction that pegs the subunits together is a p-cation interaction between h-Lys27 from SiRHP (for simplicity, amino acids from SiRHP will be designated with the prefix \u0026ldquo;h-\u0026ldquo;) and f-His258 from SiRFP (similarly, amino acids from SiRFP will be designated with the prefix \u0026ldquo;f-\u0026ldquo;) (Figs.\u0026nbsp;\u003cspan refid=\"Fig1\" class=\"InternalRef\"\u003e1\u003c/span\u003eD and E). This interaction is buttressed by f-Phe496 and f-Val500, which have previously been shown essential for SiRFP-SiRHP binding\u003csup\u003e\u003cspan citationid=\"CR13\" class=\"CitationRef\"\u003e13\u003c/span\u003e\u003c/sup\u003e. Further hydrophobic and p-stacking interactions dominate the interface. For example, h-Leu40 inserts into a pocket in SiRFP formed between f-b-sheet 17 and f-a-helix 18, which includes f-Phe496. h-Ile65 Cg2 sits 3.3 \u0026Aring; from the plane of f-Arg250\u0026rsquo;s guanidinium group (Figs.\u0026nbsp;\u003cspan refid=\"Fig1\" class=\"InternalRef\"\u003e1\u003c/span\u003eD and F), which is rotated 90\u003csup\u003eo\u003c/sup\u003e from its position in free SiRFP\u003csup\u003e\u003cspan citationid=\"CR20\" class=\"CitationRef\"\u003e20\u003c/span\u003e\u003c/sup\u003e. h-Gln72Cg also packs into a pocket formed by the backbone atoms of f-Ile247, f-Thr248, and f-Gly249, pinned in place by the h-Ile65/f-Arg250 and h-Lys73/f-His258 interactions. Farther from the interface, there is another stacking interaction between the guanidinium group from h-Arg63 and the h-Phe437 aromatic ring that stabilized the deformed helix that includes h-Ile65 (Figs.\u0026nbsp;\u003cspan refid=\"Fig1\" class=\"InternalRef\"\u003e1\u003c/span\u003eD and G). h-Phe437 is adjacent, through h-A443, to SiRHP\u0026rsquo;s iron-sulfur cluster. In this way, the binding interface between SiRFP and SiRHP reaches to the SiRHP active site through a network of hydrophobic interactions (Fig. S6).\u003c/p\u003e \u003cp\u003eNeither ionic interactions nor hydrogen bonds play a direct role in the interface, but rather stabilize the amino acids and structural elements that mediate the interface (Figs.\u0026nbsp;\u003cspan refid=\"Fig1\" class=\"InternalRef\"\u003e1\u003c/span\u003eC, D and S6). For example, a hydrogen bond network from f-Thr404Og through f-Tyr498OH and finally to f-His258Nd1 position its imidazole ring for the interaction with h-Lys73. An ionic interaction between h-Lys127Nz and h-Asp38Od2 reach across the loop, presenting h-Leu40 to project into SiRFP\u0026rsquo;s pocket. An additional ionic interaction between h-Asp61d1 and h-Arg66NHd also stabilize the deformed helix that contains h-Ile65 and turn it towards f-Arg250.\u003c/p\u003e\n\u003ch3\u003eThe SiRHP interaction with SiRFP differs from SIR interaction with ferredoxin\u003c/h3\u003e\n\u003cp\u003eg-proteobacteria couple a diflavin reductase, SiRFP, to SiRHP. In contrast, other organisms like \u003cem\u003eZea mays\u003c/em\u003e and \u003cem\u003eMycobacteria tuberculosis\u003c/em\u003e use a ferredoxin (Fd) as their reductase partner\u003csup\u003e\u003cspan citationid=\"CR15\" class=\"CitationRef\"\u003e15\u003c/span\u003e,\u003cspan citationid=\"CR17\" class=\"CitationRef\"\u003e17\u003c/span\u003e\u003c/sup\u003e. In \u003cem\u003eZ. mays\u003c/em\u003e SIR, Fd bridges SIR\u0026rsquo;s C-terminal domain 2 to a loop between the first two b-strands (amino acids Asp110 to Asn118), positioning the Fd iron-sulfur cluster near the SIR metal sites\u003csup\u003e\u003cspan citationid=\"CR26\" class=\"CitationRef\"\u003e26\u003c/span\u003e\u003c/sup\u003e. The interaction is bolstered by \u003cem\u003een face\u003c/em\u003e stacking between Fd Tyr37 and SIR Arg324 (Fig.\u0026nbsp;\u003cspan refid=\"Fig2\" class=\"InternalRef\"\u003e2\u003c/span\u003eA).\u003c/p\u003e \u003cp\u003e \u003c/p\u003e \u003cp\u003eIn SiRHP, the equivalent element between the structurally conserved b-strands 1 and 2 is considerably longer, stretching from h-Asp62 to h-Arg77 and containing a short helix from h-Arg63 to h-Glu71 (Fig.\u0026nbsp;\u003cspan refid=\"Fig2\" class=\"InternalRef\"\u003e2\u003c/span\u003eB). This loop contains the critical residues h-Gln72 and h-Lys73 that anchor the interaction with SiRFP\u0026rsquo;s FNR domain, which faces away from where Fd binds to the \u003cem\u003eZ. mays\u003c/em\u003e homolog (Fig.\u0026nbsp;\u003cspan refid=\"Fig2\" class=\"InternalRef\"\u003e2\u003c/span\u003eC). Further, the arginine in SIR that stacks with Fd Tyr37 is not conserved in SiRHP \u0026ndash; the equivalent position is h-Gly262, despite the structural conservation of the loop between b-strands 7 and 8 that contributes to siroheme binding (Fig.\u0026nbsp;\u003cspan refid=\"Fig2\" class=\"InternalRef\"\u003e2\u003c/span\u003eD).\u003c/p\u003e\n\u003ch3\u003eThe SiRHP active site loop is locked in its closed conformation\u003c/h3\u003e\n\u003cp\u003eWhen bound to SiRFP, SiRHP\u0026rsquo;s anion binding loop (h-Asn149 to h-Arg153) is in its closed position, held in place by a long, through-space interaction between h-Arg53 and h-Asn149 (Fig.\u0026nbsp;\u003cspan refid=\"Fig3\" class=\"InternalRef\"\u003e3\u003c/span\u003eA). h-Arg53 is, in turn, held in place by stacking between its guanidinium group and h-Tyr58. The ring of h-Tyr58 subsequently sits over the methyl group on the siroheme pyrroline A ring. The only other new protein/siroheme interaction identified in this now complete structure of SiRHP is an ionic bond between h-Gln60 and the propionyl group from siroheme pyrroline ring B. The siroheme is saddle-shaped, as in free SiRHP and unlike in dissimilatory SIR\u003csup\u003e\u003cspan additionalcitationids=\"CR28 CR29\" citationid=\"CR27\" class=\"CitationRef\"\u003e27\u003c/span\u003e\u0026ndash;\u003cspan citationid=\"CR30\" class=\"CitationRef\"\u003e30\u003c/span\u003e\u003c/sup\u003e. h-Arg153 is flipped away from the bound phosphate. The other three anion binding amino acids, h-Arg83, h-Lys215, and h-Lys217, remain largely unchanged from free SiRHP (Figs.\u0026nbsp;\u003cspan refid=\"Fig3\" class=\"InternalRef\"\u003e3\u003c/span\u003eA-C)\u003csup\u003e\u003cspan citationid=\"CR24\" class=\"CitationRef\"\u003e24\u003c/span\u003e\u003c/sup\u003e.\u003c/p\u003e \u003cp\u003e \u003c/p\u003e \u003cp\u003eThis conformation differs from the various redox and anion-bound structures of free SiRHP\u003csup\u003e\u003cspan citationid=\"CR24\" class=\"CitationRef\"\u003e24\u003c/span\u003e,\u003cspan additionalcitationids=\"CR32\" citationid=\"CR31\" class=\"CitationRef\"\u003e31\u003c/span\u003e\u0026ndash;\u003cspan citationid=\"CR33\" class=\"CitationRef\"\u003e33\u003c/span\u003e\u003c/sup\u003e. In the phosphate-bound, free SiRHP structure that lacks the N-terminal 80 amino acid extension\u003csup\u003e\u003cspan citationid=\"CR24\" class=\"CitationRef\"\u003e24\u003c/span\u003e\u003c/sup\u003e, the loop is flipped open such that h-Ala146-h-Ala148 are disordered (Fig.\u0026nbsp;\u003cspan refid=\"Fig3\" class=\"InternalRef\"\u003e3\u003c/span\u003eB). Upon reduction and sulfite binding, h-Arg153 flips over to interact with the smaller anion and the loop becomes ordered\u003csup\u003e\u003cspan citationid=\"CR32\" class=\"CitationRef\"\u003e32\u003c/span\u003e\u003c/sup\u003e. h-Asn149 points away from the active site. In this way, SiRFP binding to SiRHP, with the ordering of SiRHP\u0026rsquo;s N-terminus, induces an intermediate structure with elements of both the oxidized, phosphate bound and the reduced, sulfite bound conformations (Fig.\u0026nbsp;\u003cspan refid=\"Fig3\" class=\"InternalRef\"\u003e3\u003c/span\u003eD).\u003c/p\u003e \u003cp\u003eIn the original SiRHP X-ray crystal structure, the siroheme iron is significantly domed above the siroheme nitrogens, indicative of an oxidized Fe\u003csup\u003e3+\u003c/sup\u003e (Fig.\u0026nbsp;\u003cspan refid=\"Fig3\" class=\"InternalRef\"\u003e3\u003c/span\u003eB)\u003csup\u003e\u003cspan citationid=\"CR24\" class=\"CitationRef\"\u003e24\u003c/span\u003e\u003c/sup\u003e. Subsequent chemical reduction experiments show the doming flattens upon conversion to Fe\u003csup\u003e2+\u003c/sup\u003e, commensurate with release of the phosphate to allow substrate binding (Fig.\u0026nbsp;\u003cspan refid=\"Fig3\" class=\"InternalRef\"\u003e3\u003c/span\u003eB)\u003csup\u003e\u003cspan citationid=\"CR32\" class=\"CitationRef\"\u003e32\u003c/span\u003e\u003c/sup\u003e. Contemporary X-ray diffraction experiments show that this reduction is beam-induced, but within the constraints of the crystal the phosphate remains bound to the siroheme iron in the active site\u003csup\u003e\u003cspan citationid=\"CR33\" class=\"CitationRef\"\u003e33\u003c/span\u003e\u003c/sup\u003e. In this cryo-EM structure, the siroheme iron appears to be in the plane of the siroheme nitrogens, suggesting that it has also been reduced by the electron beam. The central density for the phosphate is 3.5 \u0026Aring; from the siroheme iron and its pyramidal shape is rotated such that the oxygen-iron bond is broken (Figs.\u0026nbsp;\u003cspan refid=\"Fig3\" class=\"InternalRef\"\u003e3\u003c/span\u003eA and S4C).\u003c/p\u003e\n\u003ch3\u003eSiRFP is highly mobile\u003c/h3\u003e\n\u003cp\u003eAlthough the 2D class averages in all three datasets appeared to show little orientation preference with discernable features, further analysis revealed each to have unique properties related to the SiRFP variant used to generate the dimer (Table \u003cspan refid=\"MOESM1\" class=\"InternalRef\"\u003eS1\u003c/span\u003e, Figs. S1-3).\u003c/p\u003e \u003cp\u003eSiRFP-43/SiRHP (107 kDa in mass): The 2D class averages for SiRFP-43/SiRHP appeared to show high-resolution features, however the initial 3D models were poorly aligned, likely due to a combination of small mass, preferred orientation, and limitations in grid preparation, so the refined volume did not achieve high resolution despite its simplified form. 3D variability analysis in CryoSPARC\u003csup\u003e\u003cspan citationid=\"CR34\" class=\"CitationRef\"\u003e34\u003c/span\u003e\u003c/sup\u003e did not reveal significant conformational mobility, however orientation analysis and calculation of the 3DFSC\u003csup\u003e\u003cspan citationid=\"CR35\" class=\"CitationRef\"\u003e35\u003c/span\u003e\u003c/sup\u003e showed the particles harbored a preferred orientation (Fig. S2). Nevertheless, the absence of SiRFP\u0026rsquo;s Fld domain did not alter the SiRFP-SiRHP interaction.\u003c/p\u003e \u003cp\u003eSiRFP-60D/SiRHP (123 kDa): As with SiRFP-43/SiRHP, this dimer showed 2D class averages with high-resolution features for the core of the dimer (Figs. S1B and S3A). Nevertheless, the initial models were of inconsistent structure, so we checked for heterogeneity using \u0026ldquo;3D variability\u0026rdquo; in CryoSPARC, which revealed a distinctive degree of movement in the Fld domain. To better understand the degree of flexibility for the Fld domain from the SiRFP-60D/SiRHP dimer, we performed both \u0026ldquo;3D flex\u0026rdquo; in CryoSPARC as well as \u0026ldquo;heterogenous refinement\u0026rdquo; in cryoDRGN\u003csup\u003e\u003cspan citationid=\"CR36\" class=\"CitationRef\"\u003e36\u003c/span\u003e\u003c/sup\u003e with the particles from refinement on the main heterodimer body. This analysis, anchored on the SiRFP-FNR/SiRHP dimer, identified a dramatic movement of SiRFP\u0026rsquo;s Fld domain relative to the FNR domain, swinging 20\u003csup\u003eo\u003c/sup\u003e between the most compact and most open forms (Video S1 and Fig.\u0026nbsp;\u003cspan refid=\"Fig4\" class=\"InternalRef\"\u003e4\u003c/span\u003eA). The density for the linker joining the domains is not visible. In its most compact conformation, the Fld domain reaches the canonical \u0026ldquo;closed\u0026rdquo; conformation in which the Fld domain tucks into a cavity in the FNR domain (Fig.\u0026nbsp;\u003cspan refid=\"Fig4\" class=\"InternalRef\"\u003e4\u003c/span\u003eB). In the most open conformation, the Fld domain assumes a different position to that seen in the \u0026ldquo;open\u0026rdquo; conformation determined by X-ray crystallography of the same monomeric variant and the average solution envelope determined from SANS, intermediate between the fully opened and closed conformations (Fig.\u0026nbsp;\u003cspan refid=\"Fig4\" class=\"InternalRef\"\u003e4\u003c/span\u003eC)\u003csup\u003e\u003cspan citationid=\"CR19\" class=\"CitationRef\"\u003e19\u003c/span\u003e,\u003cspan citationid=\"CR20\" class=\"CitationRef\"\u003e20\u003c/span\u003e\u003c/sup\u003e.\u003c/p\u003e \u003cp\u003e \u003c/p\u003e \u003cp\u003eSiRFP-60X/SiRHP (124 kDa): The Fld and FNR domains are covalently linked in this dimer, thus restricting the swinging domain motion of the Fld domain seen in the SiRFP-60D/SiRHP dimer. No preferred orientation was identified either in 2D or 3D analysis (Fig. S3B). Nevertheless, there is no density for the linker between the Fld and FNR domains in the highest-resolution reconstruction, which suggests that the 36 amino acid linker is exceptionally dynamic but allows SiRFP\u0026rsquo;s two domains to come into close contact to catch them in a crosslink. Locking the Fld and FNR domains in the closed conformation resulted in the most stable variant for high-resolution analysis (2.84 \u0026Aring;), discussed above.\u003c/p\u003e \u003cp\u003eWith further classification of the SiRFP-60X/SiRHP structure, a sub-class of ~\u0026thinsp;20,000 particles (10% of the particles used for highest resolution reconstruction for this dataset) emerged that revealed the linker between the two domains. The linker runs from the C-terminal end of the Fld domain (amino acid f-Ser207) to the N-terminal end of the FNR domain (amino acid f-Ile226) (Fig.\u0026nbsp;\u003cspan refid=\"Fig4\" class=\"InternalRef\"\u003e4\u003c/span\u003eD). There is neither defined secondary structure nor contact with either of the SiRFP domains, thus the linker is not constrained in a single conformation across the entire ensemble of particles even when the domains are crosslinked. The highly mobile nature of these 36 amino acids explains how the Fld domain can reposition for contact with the oxidase partner, SiRHP, presumably to mediate electron transfer.\u003c/p\u003e"},{"header":"Discussion","content":"\u003cp\u003eThe interface between the SiRHP and SiRFP subunits of the SiR holoenzyme is surprisingly minimal, driven by p-cation, hydrophobic, and planar stacking between aromatic amino acids. The non-canonical planar stacking interactions with strategic arginine residues are not as uncommon as would be expected from the canonical role that arginine plays in electrostatic interactions and are often found at important structural elements of the protein or enzyme\u003csup\u003e\u003cspan citationid=\"CR37\" class=\"CitationRef\"\u003e37\u003c/span\u003e\u003c/sup\u003e. In the case of the SiRFP/SiRHP dimer, these interactions either mediate the dimer interface itself or establish the structural platform for other amino acids to mediate the interface. Polar or electrostatic interactions play supporting roles to align the non-polar interactions because the arginine amino groups are free to make independent, more canonical, interactions with other side chains.\u003c/p\u003e \u003cp\u003eThe curious combination of non-ionic interactions that holds the complex together explains why cryo-EM analysis has been elusive \u0026ndash; the hydrophobic air-water interface in a thin film for cryogenic plunging is destructive for such interactions\u003csup\u003e\u003cspan additionalcitationids=\"CR39 CR40\" citationid=\"CR38\" class=\"CitationRef\"\u003e38\u003c/span\u003e\u0026ndash;\u003cspan citationid=\"CR41\" class=\"CitationRef\"\u003e41\u003c/span\u003e\u003c/sup\u003e. We performed extensive trial-and-error probing grid type, grid substrate, hole size, method for hydrophilizing the grid, detergents, and plunging methods to identify what seems to be a singular condition that allowed cryogenic preservation for cryo-EM analysis: the combination of FF8 and the chameleon blot-free, rapid plunging. Neither FF8 alone with traditional plunging nor chameleon alone without FF8 was sufficient. Such low-throughput screening to find the one combination of variables that allows for high-quality sample preparation highlights the need for further technology development on the front-end of cryo-EM structure determination. High-throughput data collection with direct electron detectors is not sufficient to overcome cryo-EM sample preparation challenges.\u003c/p\u003e \u003cp\u003eSurprisingly minimal structural elements determine the evolutionary pressure for siroheme-dependent SiR oxidases to partner with a diflavin reductase like SiRFP or a simpler Fd reductase. For example, in comparing the only known structures of monomeric assimilatory SiR oxidases with their reductase partners \u0026ndash; \u003cem\u003eE. coli\u003c/em\u003e SiRHP (this work) and \u003cem\u003eZ. mays\u003c/em\u003e SIR\u003csup\u003e\u003cspan citationid=\"CR26\" class=\"CitationRef\"\u003e26\u003c/span\u003e\u003c/sup\u003e \u0026ndash; a single loop extension between a core two-strand b-sheet discriminates between Fd binding in the case of \u003cem\u003eZ. mays\u003c/em\u003e SIR and the binding of SiRHP to the FNR domain of SiRFP in the case of \u003cem\u003eE. coli\u003c/em\u003e. \u003cem\u003eE. coli\u003c/em\u003e has only one \u003cem\u003efdx\u003c/em\u003e gene that does not fulfill the reductase function for SiRHP, given that strains of \u003cem\u003eE. coli\u003c/em\u003e lacking the gene encoding SiRFP (\u003cem\u003ecysJ\u003c/em\u003e) cannot grow in a low sulfur environment\u003csup\u003e\u003cspan citationid=\"CR14\" class=\"CitationRef\"\u003e14\u003c/span\u003e\u003c/sup\u003e.\u003c/p\u003e \u003cp\u003eThe SiRFP/SiRHP dimeric structure reported here shows that when SiRHP binds SiRFP, the anion-binding active site loop is locked in its closed conformation, constrained from disorder. Nevertheless, h-Arg153 is in its phosphate-binding conformation. The position of the loop is, thus, impacted by the presence of SiRHP\u0026rsquo;s N-terminus that is absent in the X-ray crystal structure. The importance of this observation is understood in comparison to a series of X-ray crystal structures in various redox states and bound to various anions\u003csup\u003e\u003cspan citationid=\"CR31\" class=\"CitationRef\"\u003e31\u003c/span\u003e,\u003cspan citationid=\"CR32\" class=\"CitationRef\"\u003e32\u003c/span\u003e\u003c/sup\u003e. This series of structures shows how the active site re-orients with the changing state, namely that the loop between h-Asn149 to h-Arg153 is disordered in the oxidized state. Upon reduction the loop becomes ordered and h-Arg153 flips to bind the substrate.\u003c/p\u003e \u003cp\u003eThe question that remains is how SiRFP\u0026rsquo;s Fld domain approaches SiRHP\u0026rsquo;s metallic active site for productive electron transfer as it does not mediate the structural interaction between the subunits. Clearly there is dramatic flexibility within SiRFP, even when the Fld and FNR domains are covalently attached by an engineered crosslink \u0026ndash; we cannot visualize the linker between the domains without extensive classification. Within this minimal dimer, the crosslink is intramolecular but in the SiR dodecameric holoenzyme, we have observed both intermolecular and intramolecular crosslinks\u003csup\u003e\u003cspan citationid=\"CR21\" class=\"CitationRef\"\u003e21\u003c/span\u003e\u003c/sup\u003e. Further, despite our extensive efforts at determining the structure of the full complex, the dodecamer has proven to be highly heterogeneous. Together, these observations support the hypothesis that within the holoenzyme, there is not a singular inter-subunit interaction between the Fld and FNR domains of any given SiRFP subunit within the octamer.\u003c/p\u003e \u003cp\u003eThe basis for SiRFP\u0026rsquo;s extreme flexibility is that the linker between the Fld and FNR domains is 36 amino acids long, compared to the 12 amino acid long linker in CPR\u003csup\u003e\u003cspan citationid=\"CR42\" class=\"CitationRef\"\u003e42\u003c/span\u003e\u003c/sup\u003e. Heterogeneity analysis of SiRFP-D in complex with SiRHP identified a dominant vector of motion parallel to the major axis of the FNR domain. This is in contrast to the highly extended conformation seen in the X-ray crystal structure of free SiRFP-D\u003csup\u003e\u003cspan citationid=\"CR20\" class=\"CitationRef\"\u003e20\u003c/span\u003e\u003c/sup\u003e. In this way, there is not a single axis along which the Fld domain moves relative to the FNR domain such that the Fld domain makes a productive interaction with SiRHP within a SiRHP/SiRFP dimer. Nevertheless, prior biochemical experimentation shows that cross-subunit interactions occur within the SiRFP octamer\u003csup\u003e\u003cspan citationid=\"CR21\" class=\"CitationRef\"\u003e21\u003c/span\u003e\u003c/sup\u003e. Thus, it is likely that a unique, single electron transfer-pathway does not exist within the full, dodecameric holoenzyme complex.\u003c/p\u003e"},{"header":"Conclusions","content":"\u003cp\u003eThe combined technological developments in cryo-EM over the past decade all played a role in determining the structure of this elusive complex, from cryogenic sample preservation to data collection to image analysis. In doing so, we show that the subunits of NADPH-dependent assimilatory sulfite reductase bind through an interface governed by the N-terminus of SiRHP and the FNR domain of SiRFP. The interaction is surprisingly minimal, governed by a set of hydrophobic and stacking interactions. Structural pairing between a siroheme-dependent oxidase and diflavin reductase or ferredoxin appears to fall to a short loop between a conserved 2-stranded b sheet. The high mobility of SiRFP\u0026rsquo;s Fld domain relative to its FNR domain likely explains how a minimal dimer maintains redox-dependent functionality and provides a mechanism for redundant electron transfer pathways within the dodecamer that efficiently performs the six-electron reduction of sulfite to sulfide.\u003c/p\u003e"},{"header":"Materials and methods","content":"\u003cdiv id=\"Sec10\" class=\"Section2\"\u003e \u003ch2\u003eProtein Expression and Purification\u003c/h2\u003e \u003cp\u003eEach SiRFP/SiRHP dimer variant was generated and purified as previously described\u003csup\u003e\u003cspan citationid=\"CR2\" class=\"CitationRef\"\u003e2\u003c/span\u003e,\u003cspan citationid=\"CR3\" class=\"CitationRef\"\u003e3\u003c/span\u003e,\u003cspan citationid=\"CR13\" class=\"CitationRef\"\u003e13\u003c/span\u003e,\u003cspan citationid=\"CR19\" class=\"CitationRef\"\u003e19\u003c/span\u003e,\u003cspan citationid=\"CR21\" class=\"CitationRef\"\u003e21\u003c/span\u003e\u003c/sup\u003e. Briefly, untagged SiRHP was co-purified with the following SiRFP variants: 1) SiRFP-43, a 43 kDa monomeric SiRFP variant lacking the N-terminal Fld domain and linker\u003csup\u003e\u003cspan citationid=\"CR19\" class=\"CitationRef\"\u003e19\u003c/span\u003e\u003c/sup\u003e; 2) SiRFP-60D, a 60 kDa monomeric variant of SiRFP generated by truncating its first 52 amino acids with an additional internal truncation of six amino acid (D-AAPSQS) in the linker that joins the Fld and FNR domains\u003csup\u003e\u003cspan citationid=\"CR20\" class=\"CitationRef\"\u003e20\u003c/span\u003e\u003c/sup\u003e; and 3) SiRFP-60X, the 60 kDa monomer with an engineered disulfide crosslink between the Fld and FNR domains, as has previously been analyzed in octameric SiRFP\u003csup\u003e\u003cspan citationid=\"CR21\" class=\"CitationRef\"\u003e21\u003c/span\u003e\u003c/sup\u003e. In SiRFP-60X, two background cysteines (C162T and C552S) were altered to avoid unwanted crosslinking and two cysteines were added (E121C and N556C) to form a disulfide bond between Fld and FNR domains of SiRFP, previously confirmed by mass spectroscopy analysis\u003csup\u003e\u003cspan citationid=\"CR21\" class=\"CitationRef\"\u003e21\u003c/span\u003e\u003c/sup\u003e. DNA sequencing confirmed the presence of all deletions and mutations in the variants.\u003c/p\u003e \u003cp\u003e \u003cem\u003eE. coli\u003c/em\u003e LMG194 cells (Invitrogen, Carlsbad, CA, USA) were transformed with the pBAD plasmid containing the genes encoding SiRFP-60D, SiRFP-60X, SiRFP-43 or SiRHP. All variants were expressed independently. Dimers were formed as followed: for SiRFP-43/SiRHP and SiRFP-60D/SiRHP, cells were mixed before lysis and the dimers co-purified. The SiRFP-60X/SiRHP dimer was assembled by reconstituting SiRFP-60X with 2 Eq SiRHP for 1 h on ice, as before\u003csup\u003e\u003cspan citationid=\"CR2\" class=\"CitationRef\"\u003e2\u003c/span\u003e,\u003cspan citationid=\"CR19\" class=\"CitationRef\"\u003e19\u003c/span\u003e,\u003cspan citationid=\"CR21\" class=\"CitationRef\"\u003e21\u003c/span\u003e\u003c/sup\u003e. A six-histidine tag was present in all SiRFP variants whereas SiRHP was untagged. Each recombinant \u003cem\u003eE. coli\u003c/em\u003e strain was grown and induced at 25\u0026deg;C with 0.05% L-arabinose. SiRFP-43, SiRFP-60D and SiRHP cells were expressed for 4 hr whereas SiRFP-60X was expressed overnight. All variants were purified using a combination of Ni-NTA affinity chromatography (Cytiva, Marlborough, MA, USA), anion exchange HiTrap-Q HP chromatography (Cytiva, Marlborough, MA, USA) and Sephacryl S300-HR size exclusion chromatography (Cytiva, Marlborough, MA, USA) using previously optimized SPG buffers (17 mM succinic acid, 574 mM sodium dihydrogen phosphate, pH 6.8, 374 mM glycine, 200 mM NaCl)\u003csup\u003e\u003cspan citationid=\"CR2\" class=\"CitationRef\"\u003e2\u003c/span\u003e,\u003cspan citationid=\"CR13\" class=\"CitationRef\"\u003e13\u003c/span\u003e,\u003cspan citationid=\"CR19\" class=\"CitationRef\"\u003e19\u003c/span\u003e,\u003cspan citationid=\"CR21\" class=\"CitationRef\"\u003e21\u003c/span\u003e\u003c/sup\u003e. All variants have been extensively characterized with UV-Vis spectroscopy, size exclusion chromatography, SDS PAGE analysis for correct stoichiometry and with SANS for solution monodispersity (Fig. S7)\u003csup\u003e\u003cspan citationid=\"CR2\" class=\"CitationRef\"\u003e2\u003c/span\u003e,\u003cspan citationid=\"CR13\" class=\"CitationRef\"\u003e13\u003c/span\u003e,\u003cspan citationid=\"CR19\" class=\"CitationRef\"\u003e19\u003c/span\u003e,\u003cspan citationid=\"CR21\" class=\"CitationRef\"\u003e21\u003c/span\u003e\u003c/sup\u003e.\u003c/p\u003e \u003c/div\u003e \u003cdiv id=\"Sec11\" class=\"Section2\"\u003e \u003ch2\u003eCryo-EM sample preparation\u003c/h2\u003e \u003cp\u003eSiRFP-43/SiRHP, SiRFP-60D/SiRHP and SiRFP-60X/SiRHP were all plunged using the chameleon blotless plunging system (SPT Labtech, Melbourn, UK) at 10 mg/mL protein. chameleon self-wicking grids\u003csup\u003e\u003cspan citationid=\"CR23\" class=\"CitationRef\"\u003e23\u003c/span\u003e\u003c/sup\u003e were backed with 18-gauge gold (Ted Pella, Redding CA, USA) on an Auto 306 vacuum coater (BOC Edwards, West Sussex, UK). The following glow discharge (GD)/wicking times (WTs) were used: SiRFP-43/SiRHP: 30 s GD, 130 ms WT; SiRFP-60D/SiRHP: 80 s GD, 195 ms WTs; SiRFP-60X/SiRHP 45 s GD, 175 ms WTs. All samples were premixed with FC-8 detergent\u003csup\u003e\u003cspan citationid=\"CR22\" class=\"CitationRef\"\u003e22\u003c/span\u003e\u003c/sup\u003e at 2 mM final concentration. The plunged grids were clipped and then screened for high quality with a Titan Krios operating the Leginon software package\u003csup\u003e\u003cspan citationid=\"CR43\" class=\"CitationRef\"\u003e43\u003c/span\u003e\u003c/sup\u003e.\u003c/p\u003e \u003cp\u003e \u003cem\u003eData collection and processing\u003c/em\u003e \u003c/p\u003e \u003cp\u003e \u003col\u003e \u003cspan\u003e \u003cli\u003e \u003cp\u003eSiRFP-43/SiRHP: 14,600 movies were collected at 300 KV on a Titan Krios using a GATAN K3 camera with 0.844 \u0026Aring;/pixel. After motion correction with Motioncor2\u003csup\u003e44\u003c/sup\u003e in the Relion GUI\u003csup\u003e\u003cspan citationid=\"CR45\" class=\"CitationRef\"\u003e45\u003c/span\u003e\u003c/sup\u003e, CTF estimation was performed using CTFFIND4\u003csup\u003e46\u003c/sup\u003e. Particles were picked by \u0026ldquo;blob picker\u0026rdquo; followed by \u0026ldquo;template picker\u0026rdquo; in CryoSPARC\u003csup\u003e\u003cspan citationid=\"CR34\" class=\"CitationRef\"\u003e34\u003c/span\u003e\u003c/sup\u003e. Initial 2D classification, followed by multiple rounds of 2D class selection/classification, identified 1,500,000 particles that were used for initial map building and non-uniform refinement in CryoSPARC. This process achieved a reported 3.6 \u0026Aring;-resolution map of the minimal dimer. However, the map features did not reflect the reported resolution. \u0026ldquo;Orientation diagnostic\u0026rdquo; job from CryoSPARC was performed to confirm this interpretation. To further confirm and diminish the preferred orientation issue, \u0026ldquo;Rebalance 2D\u0026rdquo; was performed to put particles into 7 super classes and limit the total particles in each superclass down to total of 350,000 and 105,000 particles. Non-uniform refinement\u003csup\u003e\u003cspan citationid=\"CR47\" class=\"CitationRef\"\u003e47\u003c/span\u003e\u003c/sup\u003e followed by orientation diagnostics performed in CryoSPARC produced a more accurate resolution (4.31 \u0026Aring;, 4.74 \u0026Aring; respectively) and better quality map. deepEMhancer\u003csup\u003e\u003cspan citationid=\"CR48\" class=\"CitationRef\"\u003e48\u003c/span\u003e\u003c/sup\u003e was used to sharpen these maps (Fig. S2). Particle picking performed by the TOPAZ algorithm\u003csup\u003e\u003cspan citationid=\"CR49\" class=\"CitationRef\"\u003e49\u003c/span\u003e\u003c/sup\u003e gave the same result.\u003c/p\u003e \u003c/li\u003e \u003c/span\u003e \u003cspan\u003e \u003cli\u003e \u003cp\u003eSiRFP-60D/SiRHP: 25,488 movies were collected at 300 KV on a Titan Krios using a GATAN K3 camera with 0.844 \u0026Aring;/pixel. Motion correction, CTF estimation, and particle picking were performed as for SiRFP-43/SiRHP. After multiple rounds of 2D classification, ~\u0026thinsp;550,000 particles were used for initial map building and non-uniform refinement in CryoSPARC to achieve the final 3.49 \u0026Aring;-resolution structure of the dimer, masked to omit SiRFP\u0026rsquo;s Fld domain. Multiple refinements with different masking were performed. Masking to include the whole complex, including SiRFP\u0026rsquo;s Fld domain, resulted in a low resolution reconstruction for the Fld domain. Particles were then down sampled and imported to cryoDRGN to perform heterogenous refinement, giving a series of volumes showing the Fld domain movement. CryoSPARC 3D Flex\u003csup\u003e\u003cspan citationid=\"CR50\" class=\"CitationRef\"\u003e50\u003c/span\u003e\u003c/sup\u003e gave the same result as cryoDRGN.\u003c/p\u003e \u003c/li\u003e \u003c/span\u003e \u003cspan\u003e \u003cli\u003e \u003cp\u003eSiRFP-60X/SiRHP: ~10,000 movies were collected at 300 KV on a Titan Krios using a DE Apollo camera with 0.765 \u0026Aring;/pixel. Motion correction, CTF estimation, and particle picking were performed as above. Multiple rounds of 2D classifications identified\u0026thinsp;~\u0026thinsp;185,000 particles that were used for initial map building and non-uniform refinement in CryoSPARC to achieve a 2.84 \u0026Aring;-resolution structure of the entire dimer, including SiRFP\u0026rsquo;s Fld domain. Further 3D variability was performed to track the Fld linker. After multiple rounds of 3D classification 10% of the particles (~\u0026thinsp;20,000) were used to perform a non-uniform refinement, resulted in a 3.61 \u0026Aring; overall resolution structure, sharpened by deepEMhancer including the linker.\u003c/p\u003e \u003c/li\u003e \u003c/span\u003e \u003c/ol\u003e \u003c/p\u003e \u003cp\u003e \u003cem\u003eModel building and refinement.\u003c/em\u003e \u003c/p\u003e \u003cp\u003eModel building was initiated using the \u0026ldquo;fit in map\u0026rdquo; algorithm in Chimera\u003csup\u003e\u003cspan citationid=\"CR51\" class=\"CitationRef\"\u003e51\u003c/span\u003e\u003c/sup\u003e using the atomic model of SiRFP from PDB ID 6EFV\u003csup\u003e\u003cspan citationid=\"CR20\" class=\"CitationRef\"\u003e20\u003c/span\u003e\u003c/sup\u003e, fitting each of the Fld or FNR domains independently or the atomic model of SiRHP generated in AlphaFold\u003csup\u003e\u003cspan citationid=\"CR52\" class=\"CitationRef\"\u003e52\u003c/span\u003e,\u003cspan citationid=\"CR53\" class=\"CitationRef\"\u003e53\u003c/span\u003e\u003c/sup\u003e to capture its N-terminal 80 amino acids. Iterative real-space refinement in PHENIX\u003csup\u003e\u003cspan citationid=\"CR54\" class=\"CitationRef\"\u003e54\u003c/span\u003e\u003c/sup\u003e with manual fitting in Coot\u003csup\u003e\u003cspan citationid=\"CR55\" class=\"CitationRef\"\u003e55\u003c/span\u003e\u003c/sup\u003e yielded a model with a correlation of 0.85 (Table S3). The model deposited in PDB as 9C91and the map deposited in EMDB as EMD-45359. A stack from the final 3D reconstruction particles was deposited in the EMPIAR database under the accession number EMPIAR-12180.\u003c/p\u003e \u003c/div\u003e"},{"header":"Declarations","content":"\u003ch2\u003eAcknowledgements and funding\u003c/h2\u003e \u003cp\u003eWe kindly thank Dr.s Scott Stagg and Christopher Stroupe for helpful discussions. Some of this work was performed at the Simons Electron Microscopy Center at the New York Structural Biology Center, with major support from the Simons Foundation (SF349247). Florida State University supports cryo-EM in the Biological Imaging Resource Center, which houses the following equipment used in this study: a Gatan Solaris Plasma Cleaner (NIH grant S10 RR024564), a Hitachi HT7800 (NSF grant 2017869), a ThermoFisher Vitrobot Mark IV (NIH grant S10 RR024564), an SPI chameleon plunging system (NIH grant R24 GM145964), a ThermoFisher Titan Krios (NIH grant S10 RR025080), and a DE Apollo direct electron detector (NIH grant R35 GM139616). This work was further supported by National Science Foundation grants MCB1856502 and CHE1904612 to M.E.S.\u003c/p\u003e \u003ch2\u003eAuthor contributions\u003c/h2\u003e \u003cp\u003eBGE participated in data collection and performed the data analysis. NW and BGE prepared the protein specimen. KN and MA prepared cryogenic specimen and participated in data collection. IA optimized protein specimen preparation. AW and JHM supervised cryogenic specimen preparation and data collection. MES conceived of the study and oversaw experimentation.\u003c/p\u003e"},{"header":"References","content":"\u003col\u003e\u003cli\u003e\u003cspan\u003eSiegel LM, Davis PS (1974) Reduced nicotinamide adenine dinucleotide phosphate-sulfite reductase of enterobacteria. IV. The Escherichia coli hemoflavoprotein: subunit structure and dissociation into hemoprotein and flavoprotein components. J Biol Chem 249:1587\u0026ndash;1598\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eMurray DT et al (2022) Neutron scattering maps the higher-order assembly of NADPH-dependent assimilatory sulfite reductase. Biophys J 121:1799\u0026ndash;1812. \u003cspan class=\"ExternalRef\"\u003e\u003cspan class=\"RefSource\"\u003ehttps://doi.org/10.1016/j.bpj.2022.04.021\u003c/span\u003e\u003cspan address=\"10.1016/j.bpj.2022.04.021\" targettype=\"DOI\" class=\"RefTarget\"\u003e\u003c/span\u003e\u003c/span\u003e\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eAskenasy I et al (2015) The N-terminal Domain of Escherichia coli Assimilatory NADPH-Sulfite Reductase Hemoprotein Is an Oligomerization Domain That Mediates Holoenzyme Assembly. J Biol Chem 290:19319\u0026ndash;19333. \u003cspan class=\"ExternalRef\"\u003e\u003cspan class=\"RefSource\"\u003ehttps://doi.org/10.1074/jbc.M115.662379\u003c/span\u003e\u003cspan address=\"10.1074/jbc.M115.662379\" targettype=\"DOI\" class=\"RefTarget\"\u003e\u003c/span\u003e\u003c/span\u003e\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eWang M et al (1997) Three-dimensional structure of NADPH-cytochrome P450 reductase: prototype for FMN- and FAD-containing enzymes. Proc Natl Acad Sci U S A 94:8411\u0026ndash;8416\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eGirvan HM et al (2006) Flavocytochrome P450 BM3 and the origin of CYP102 fusion species. Biochem Soc Trans 34:1173\u0026ndash;1177. \u003cspan class=\"ExternalRef\"\u003e\u003cspan class=\"RefSource\"\u003ehttps://doi.org/10.1042/BST0341173\u003c/span\u003e\u003cspan address=\"10.1042/BST0341173\" targettype=\"DOI\" class=\"RefTarget\"\u003e\u003c/span\u003e\u003c/span\u003e\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eGarcin ED et al (2004) Structural basis for isozyme-specific regulation of electron transfer in nitric-oxide synthase. J Biol Chem 279:37918\u0026ndash;37927. \u003cspan class=\"ExternalRef\"\u003e\u003cspan class=\"RefSource\"\u003ehttps://doi.org/10.1074/jbc.M406204200\u003c/span\u003e\u003cspan address=\"10.1074/jbc.M406204200\" targettype=\"DOI\" class=\"RefTarget\"\u003e\u003c/span\u003e\u003c/span\u003e\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eZhang J et al (2001) Crystal structure of the FAD/NADPH-binding domain of rat neuronal nitric-oxide synthase. Comparisons with NADPH-cytochrome P450 oxidoreductase. J Biol Chem 276:37506\u0026ndash;37513. \u003cspan class=\"ExternalRef\"\u003e\u003cspan class=\"RefSource\"\u003ehttps://doi.org/10.1074/jbc.M105503200\u003c/span\u003e\u003cspan address=\"10.1074/jbc.M105503200\" targettype=\"DOI\" class=\"RefTarget\"\u003e\u003c/span\u003e\u003c/span\u003e\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eOlteanu H, Banerjee R (2001) Human methionine synthase reductase, a soluble P-450 reductase-like dual flavoprotein, is sufficient for NADPH-dependent methionine synthase activation. J Biol Chem 276:35558\u0026ndash;35563. \u003cspan class=\"ExternalRef\"\u003e\u003cspan class=\"RefSource\"\u003ehttps://doi.org/10.1074/jbc.M103707200\u003c/span\u003e\u003cspan address=\"10.1074/jbc.M103707200\" targettype=\"DOI\" class=\"RefTarget\"\u003e\u003c/span\u003e\u003c/span\u003e\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eFreeman SL et al (2018) Solution structure of the cytochrome P450 reductase-cytochrome. J Biol Chem 293:5210\u0026ndash;5219. \u003cspan class=\"ExternalRef\"\u003e\u003cspan class=\"RefSource\"\u003ehttps://doi.org/10.1074/jbc.RA118.001941\u003c/span\u003e\u003cspan address=\"10.1074/jbc.RA118.001941\" targettype=\"DOI\" class=\"RefTarget\"\u003e\u003c/span\u003e\u003c/span\u003e\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eZhang L et al (2020) Structural insight into the electron transfer pathway of a self-sufficient P450 monooxygenase. Nat Commun 11:2676. \u003cspan class=\"ExternalRef\"\u003e\u003cspan class=\"RefSource\"\u003ehttps://doi.org/10.1038/s41467-020-16500-5\u003c/span\u003e\u003cspan address=\"10.1038/s41467-020-16500-5\" targettype=\"DOI\" class=\"RefTarget\"\u003e\u003c/span\u003e\u003c/span\u003e\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eZhang H et al (2018) The full-length cytochrome P450 enzyme CYP102A1 dimerizes at its reductase domains and has flexible heme domains for efficient catalysis. J Biol Chem 293:7727\u0026ndash;7736. \u003cspan class=\"ExternalRef\"\u003e\u003cspan class=\"RefSource\"\u003ehttps://doi.org/10.1074/jbc.RA117.000600\u003c/span\u003e\u003cspan address=\"10.1074/jbc.RA117.000600\" targettype=\"DOI\" class=\"RefTarget\"\u003e\u003c/span\u003e\u003c/span\u003e\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eZeghouf M, Fontecave M, Coves J (2000) A simplifed functional version of the Escherichia coli sulfite reductase. J Biol Chem 275:37651\u0026ndash;37656. \u003cspan class=\"ExternalRef\"\u003e\u003cspan class=\"RefSource\"\u003ehttps://doi.org/10.1074/jbc.M005619200\u003c/span\u003e\u003cspan address=\"10.1074/jbc.M005619200\" targettype=\"DOI\" class=\"RefTarget\"\u003e\u003c/span\u003e\u003c/span\u003eM005619200 [pii]\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eAskenasy I et al (2018) Structure-Function Relationships in the Oligomeric NADPH-Dependent Assimilatory Sulfite Reductase. Biochemistry 57:3764\u0026ndash;3772. \u003cspan class=\"ExternalRef\"\u003e\u003cspan class=\"RefSource\"\u003ehttps://doi.org/10.1021/acs.biochem.8b00446\u003c/span\u003e\u003cspan address=\"10.1021/acs.biochem.8b00446\" targettype=\"DOI\" class=\"RefTarget\"\u003e\u003c/span\u003e\u003c/span\u003e\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eSiegel LM, Murphy MJ, Kamin H (1973) Reduced nicotinamide adenine dinucleotide phosphate-sulfite reductase of enterobacteria. I. The Escherichia coli hemoflavoprotein: molecular parameters and prosthetic groups. J Biol Chem 248:251\u0026ndash;264\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eNakayama M, Akashi T, Hase T (2000) Plant sulfite reductase: molecular structure, catalytic function and interaction with ferredoxin. J Inorg Biochem 82:27\u0026ndash;32\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eSchnell R, Sandalova T, Hellman U, Lindqvist Y, Schneider G (2005) Siroheme- and [Fe4-S4]-dependent NirA from Mycobacterium tuberculosis is a sulfite reductase with a covalent Cys-Tyr bond in the active site. J Biol Chem 280:27319\u0026ndash;27328. https://doi.org/M502560200\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003ePinto R, Harrison JS, Hsu T, Jacobs WR Jr., Leyh TS (2007) Sulfite reduction in mycobacteria. J Bacteriol 189:6714\u0026ndash;6722. \u003cspan class=\"ExternalRef\"\u003e\u003cspan class=\"RefSource\"\u003ehttps://doi.org/10.1128/jb.00487-07\u003c/span\u003e\u003cspan address=\"10.1128/jb.00487-07\" targettype=\"DOI\" class=\"RefTarget\"\u003e\u003c/span\u003e\u003c/span\u003e\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eJespersen M, Pierik AJ, Wagner T (2023) Structures of the sulfite detoxifying F 420-dependent enzyme from Methanococcales. Nat Chem Biol 19:695\u0026ndash;702. \u003cspan class=\"ExternalRef\"\u003e\u003cspan class=\"RefSource\"\u003ehttps://doi.org/10.1038/s41589-022-01232-y\u003c/span\u003e\u003cspan address=\"10.1038/s41589-022-01232-y\" targettype=\"DOI\" class=\"RefTarget\"\u003e\u003c/span\u003e\u003c/span\u003e\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eMurray DT, Weiss KL, Stanley CB, Nagy G, Stroupe ME (2021) Small-angle neutron scattering solution structures of NADPH-dependent sulfite reductase. J Struct Biol 213:107724. \u003cspan class=\"ExternalRef\"\u003e\u003cspan class=\"RefSource\"\u003ehttps://doi.org/10.1016/j.jsb.2021.107724\u003c/span\u003e\u003cspan address=\"10.1016/j.jsb.2021.107724\" targettype=\"DOI\" class=\"RefTarget\"\u003e\u003c/span\u003e\u003c/span\u003e\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eTavolieri AM et al (2019) NADPH-dependent sulfite reductase flavoprotein adopts an extended conformation unique to this diflavin reductase. J Struct Biol. \u003cspan class=\"ExternalRef\"\u003e\u003cspan class=\"RefSource\"\u003ehttps://doi.org/10.1016/j.jsb.2019.01.001\u003c/span\u003e\u003cspan address=\"10.1016/j.jsb.2019.01.001\" targettype=\"DOI\" class=\"RefTarget\"\u003e\u003c/span\u003e\u003c/span\u003e\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eWalia N et al (2023) Domain crossover in the reductase subunit of NADPH-dependent assimilatory sulfite reductase. J Struct Biol 215:108028. \u003cspan class=\"ExternalRef\"\u003e\u003cspan class=\"RefSource\"\u003ehttps://doi.org/10.1016/j.jsb.2023.108028\u003c/span\u003e\u003cspan address=\"10.1016/j.jsb.2023.108028\" targettype=\"DOI\" class=\"RefTarget\"\u003e\u003c/span\u003e\u003c/span\u003e\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eKampjut D, Steiner J, Sazanov LA (2021) Cryo-EM grid optimization for membrane proteins. iScience 24:102139. \u003cspan class=\"ExternalRef\"\u003e\u003cspan class=\"RefSource\"\u003ehttps://doi.org/10.1016/j.isci.2021.102139\u003c/span\u003e\u003cspan address=\"10.1016/j.isci.2021.102139\" targettype=\"DOI\" class=\"RefTarget\"\u003e\u003c/span\u003e\u003c/span\u003e\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eLevitz TS et al (2022) Approaches to Using the Chameleon: Robust, Automated, Fast-Plunge cryoEM Specimen Preparation. Front Mol Biosci 9:903148. \u003cspan class=\"ExternalRef\"\u003e\u003cspan class=\"RefSource\"\u003ehttps://doi.org/10.3389/fmolb.2022.903148\u003c/span\u003e\u003cspan address=\"10.3389/fmolb.2022.903148\" targettype=\"DOI\" class=\"RefTarget\"\u003e\u003c/span\u003e\u003c/span\u003e\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eCrane BR, Siegel LM, Getzoff ED (1995) Sulfite reductase structure at 1.6 A: evolution and catalysis for reduction of inorganic anions. Science 270:59\u0026ndash;67\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eCrane BR, Getzoff ED (1996) The relationship between structure and function for the sulfite reductases. Curr Opin Struct Biol 6:744\u0026ndash;756. \u003cspan class=\"ExternalRef\"\u003e\u003cspan class=\"RefSource\"\u003ehttps://doi.org/S0959-440X\u003c/span\u003e\u003cspan address=\"https://doi.org/S0959-440X\" targettype=\"URL\" class=\"RefTarget\"\u003e\u003c/span\u003e\u003c/span\u003e (96)80003-0 [pii]\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eKim JY, Nakayama M, Toyota H, Kurisu G, Hase T (2016) Structural and mutational studies of an electron transfer complex of maize sulfite reductase and ferredoxin. J Biochem. \u003cspan class=\"ExternalRef\"\u003e\u003cspan class=\"RefSource\"\u003ehttps://doi.org/10.1093/jb/mvw016\u003c/span\u003e\u003cspan address=\"10.1093/jb/mvw016\" targettype=\"DOI\" class=\"RefTarget\"\u003e\u003c/span\u003e\u003c/span\u003e\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eOliveira TF et al (2008) The crystal structure of Desulfovibrio vulgaris dissimilatory sulfite reductase bound to DsrC provides novel insights into the mechanism of sulfate respiration. J Biol Chem 283:34141\u0026ndash;34149. https://doi.org/M805643200\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eOliveira TF et al (2011) Structural insights into dissimilatory sulfite reductases: structure of desulforubidin from desulfomicrobium norvegicum. Front Microbiol 2:71. \u003cspan class=\"ExternalRef\"\u003e\u003cspan class=\"RefSource\"\u003ehttps://doi.org/10.3389/fmicb.2011.00071\u003c/span\u003e\u003cspan address=\"10.3389/fmicb.2011.00071\" targettype=\"DOI\" class=\"RefTarget\"\u003e\u003c/span\u003e\u003c/span\u003e\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eParey K, Warkentin E, Kroneck PM, Ermler U (2010) Reaction cycle of the dissimilatory sulfite reductase from Archaeoglobus fulgidus. Biochemistry 49:8912\u0026ndash;8921. \u003cspan class=\"ExternalRef\"\u003e\u003cspan class=\"RefSource\"\u003ehttps://doi.org/10.1021/bi100781f\u003c/span\u003e\u003cspan address=\"10.1021/bi100781f\" targettype=\"DOI\" class=\"RefTarget\"\u003e\u003c/span\u003e\u003c/span\u003e\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eSchiffer A et al (2008) Structure of the dissimilatory sulfite reductase from the hyperthermophilic archaeon Archaeoglobus fulgidus. \u003cem\u003eJ Mol Biol\u003c/em\u003e 379, 1063\u0026ndash;1074 \u003cspan class=\"ExternalRef\"\u003e\u003cspan class=\"RefSource\"\u003ehttps://doi.org/S0022-2836(08)00456-7\u003c/span\u003e\u003cspan address=\"https://doi.org/S0022-2836(08)00456-7\" targettype=\"URL\" class=\"RefTarget\"\u003e\u003c/span\u003e\u003c/span\u003e [pii] \u003cspan class=\"ExternalRef\"\u003e\u003cspan class=\"RefSource\"\u003e10.1016/j.jmb.2008.04.027\u003c/span\u003e\u003cspan address=\"10.1016/j.jmb.2008.04.027\" targettype=\"DOI\" class=\"RefTarget\"\u003e\u003c/span\u003e\u003c/span\u003e\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eCrane BR, Siegel LM, Getzoff ED (1997) Probing the catalytic mechanism of sulfite reductase by X-ray crystallography: structures of the Escherichia coli hemoprotein in complex with substrates, inhibitors, intermediates, and products. Biochemistry 36:12120\u0026ndash;12137. \u003cspan class=\"ExternalRef\"\u003e\u003cspan class=\"RefSource\"\u003ehttps://doi.org/10.1021/bi971066i bi971066i\u003c/span\u003e\u003cspan address=\"10.1021/bi971066i bi971066i\" targettype=\"DOI\" class=\"RefTarget\"\u003e\u003c/span\u003e\u003c/span\u003e[pii]\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eCrane BR, Siegel LM, Getzoff ED (1997) Structures of the siroheme- and Fe4S4-containing active center of sulfite reductase in different states of oxidation: heme activation via reduction-gated exogenous ligand exchange. Biochemistry 36:12101\u0026ndash;12119. \u003cspan class=\"ExternalRef\"\u003e\u003cspan class=\"RefSource\"\u003ehttps://doi.org/10.1021/bi971065q bi971065q\u003c/span\u003e\u003cspan address=\"10.1021/bi971065q bi971065q\" targettype=\"DOI\" class=\"RefTarget\"\u003e\u003c/span\u003e\u003c/span\u003e[pii]\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eSmith KW, Stroupe ME (2012) Mutational analysis of sulfite reductase hemoprotein reveals the mechanism for coordinated electron and proton transfer. Biochemistry 51:9857\u0026ndash;9868. \u003cspan class=\"ExternalRef\"\u003e\u003cspan class=\"RefSource\"\u003ehttps://doi.org/10.1021/bi300947a\u003c/span\u003e\u003cspan address=\"10.1021/bi300947a\" targettype=\"DOI\" class=\"RefTarget\"\u003e\u003c/span\u003e\u003c/span\u003e\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003ePunjani A, Rubinstein JL, Fleet DJ, Brubaker MA (2017) cryoSPARC: algorithms for rapid unsupervised cryo-EM structure determination. Nat Methods 14:290\u0026ndash;296. \u003cspan class=\"ExternalRef\"\u003e\u003cspan class=\"RefSource\"\u003ehttps://doi.org/10.1038/nmeth.4169\u003c/span\u003e\u003cspan address=\"10.1038/nmeth.4169\" targettype=\"DOI\" class=\"RefTarget\"\u003e\u003c/span\u003e\u003c/span\u003e\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eAiyer S, Zhang C, Baldwin PR, Lyumkis D (2021) Evaluating Local and Directional Resolution of Cryo-EM Density Maps. Methods Mol Biol 2215:161\u0026ndash;187. \u003cspan class=\"ExternalRef\"\u003e\u003cspan class=\"RefSource\"\u003ehttps://doi.org/10.1007/978-1-0716-0966-8_8\u003c/span\u003e\u003cspan address=\"10.1007/978-1-0716-0966-8_8\" targettype=\"DOI\" class=\"RefTarget\"\u003e\u003c/span\u003e\u003c/span\u003e\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eZhong ED, Bepler T, Berger B, Davis JH (2021) CryoDRGN: reconstruction of heterogeneous cryo-EM structures using neural networks. Nat Methods 18:176\u0026ndash;185. \u003cspan class=\"ExternalRef\"\u003e\u003cspan class=\"RefSource\"\u003ehttps://doi.org/10.1038/s41592-020-01049-4\u003c/span\u003e\u003cspan address=\"10.1038/s41592-020-01049-4\" targettype=\"DOI\" class=\"RefTarget\"\u003e\u003c/span\u003e\u003c/span\u003e\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eFlocco MM, Mowbray SL (1994) Planar stacking interactions of arginine and aromatic side-chains in proteins. J Mol Biol 235:709\u0026ndash;717. \u003cspan class=\"ExternalRef\"\u003e\u003cspan class=\"RefSource\"\u003ehttps://doi.org/10.1006/jmbi.1994.1022\u003c/span\u003e\u003cspan address=\"10.1006/jmbi.1994.1022\" targettype=\"DOI\" class=\"RefTarget\"\u003e\u003c/span\u003e\u003c/span\u003e\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eD'Imprima E et al (2019) Protein denaturation at the air-water interface and how to prevent it. Elife 8. \u003cspan class=\"ExternalRef\"\u003e\u003cspan class=\"RefSource\"\u003ehttps://doi.org/10.7554/eLife.42747\u003c/span\u003e\u003cspan address=\"10.7554/eLife.42747\" targettype=\"DOI\" class=\"RefTarget\"\u003e\u003c/span\u003e\u003c/span\u003e\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eGlaeser RM et al (2016) Factors that Influence the Formation and Stability of Thin, Cryo-EM Specimens. Biophys J 110:749\u0026ndash;755. \u003cspan class=\"ExternalRef\"\u003e\u003cspan class=\"RefSource\"\u003ehttps://doi.org/10.1016/j.bpj.2015.07.050\u003c/span\u003e\u003cspan address=\"10.1016/j.bpj.2015.07.050\" targettype=\"DOI\" class=\"RefTarget\"\u003e\u003c/span\u003e\u003c/span\u003e\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eGlaeser RM, Han BG (2017) Opinion: hazards faced by macromolecules when confined to thin aqueous films. Biophys Rep 3:1\u0026ndash;7. \u003cspan class=\"ExternalRef\"\u003e\u003cspan class=\"RefSource\"\u003ehttps://doi.org/10.1007/s41048-016-0026-3\u003c/span\u003e\u003cspan address=\"10.1007/s41048-016-0026-3\" targettype=\"DOI\" class=\"RefTarget\"\u003e\u003c/span\u003e\u003c/span\u003e\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eGlaeser RM, PROTEINS, INTERFACES, AND, CRYO-EM GRIDS (2018) \u003cem\u003eCurr Opin Colloid Interface Sci\u003c/em\u003e 34, 1\u0026ndash;8 \u003cspan class=\"ExternalRef\"\u003e\u003cspan class=\"RefSource\"\u003ehttps://doi.org/10.1016/j.cocis.2017.12.009\u003c/span\u003e\u003cspan address=\"10.1016/j.cocis.2017.12.009\" targettype=\"DOI\" class=\"RefTarget\"\u003e\u003c/span\u003e\u003c/span\u003e\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eXia C et al (2011) Conformational changes of NADPH-cytochrome P450 oxidoreductase are essential for catalysis and cofactor binding. J Biol Chem 286:16246\u0026ndash;16260. \u003cspan class=\"ExternalRef\"\u003e\u003cspan class=\"RefSource\"\u003ehttps://doi.org/10.1074/jbc.M111.230532\u003c/span\u003e\u003cspan address=\"10.1074/jbc.M111.230532\" targettype=\"DOI\" class=\"RefTarget\"\u003e\u003c/span\u003e\u003c/span\u003e\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eSuloway C et al (2005) Automated molecular microscopy: the new Leginon system. J Struct Biol 151:41\u0026ndash;60. \u003cspan class=\"ExternalRef\"\u003e\u003cspan class=\"RefSource\"\u003ehttps://doi.org/S1047-8477(\u003c/span\u003e\u003cspan address=\"https://doi.org/S1047-8477(\" targettype=\"URL\" class=\"RefTarget\"\u003e\u003c/span\u003e\u003c/span\u003e 05)00072-9 [pii] \u003cspan class=\"ExternalRef\"\u003e\u003cspan class=\"RefSource\"\u003e10.1016/j.jsb.2005.03.010\u003c/span\u003e\u003cspan address=\"10.1016/j.jsb.2005.03.010\" targettype=\"DOI\" class=\"RefTarget\"\u003e\u003c/span\u003e\u003c/span\u003e\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eZheng SQ et al (2017) MotionCor2: anisotropic correction of beam-induced motion for improved cryo-electron microscopy. Nat Methods 14:331\u0026ndash;332. \u003cspan class=\"ExternalRef\"\u003e\u003cspan class=\"RefSource\"\u003ehttps://doi.org/10.1038/nmeth.4193\u003c/span\u003e\u003cspan address=\"10.1038/nmeth.4193\" targettype=\"DOI\" class=\"RefTarget\"\u003e\u003c/span\u003e\u003c/span\u003e\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eScheres SH (2012) RELION: implementation of a Bayesian approach to cryo-EM structure determination. J Struct Biol 180:519\u0026ndash;530. \u003cspan class=\"ExternalRef\"\u003e\u003cspan class=\"RefSource\"\u003ehttps://doi.org/10.1016/j.jsb.2012.09.006\u003c/span\u003e\u003cspan address=\"10.1016/j.jsb.2012.09.006\" targettype=\"DOI\" class=\"RefTarget\"\u003e\u003c/span\u003e\u003c/span\u003e\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eRohou A, Grigorieff N (2015) CTFFIND4: Fast and accurate defocus estimation from electron micrographs. J Struct Biol 192:216\u0026ndash;221\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003ePunjani A, Zhang H, Fleet DJ (2020) Non-uniform refinement: adaptive regularization improves single-particle cryo-EM reconstruction. Nat Methods 17:1214\u0026ndash;1221. \u003cspan class=\"ExternalRef\"\u003e\u003cspan class=\"RefSource\"\u003ehttps://doi.org/10.1038/s41592-020-00990-8\u003c/span\u003e\u003cspan address=\"10.1038/s41592-020-00990-8\" targettype=\"DOI\" class=\"RefTarget\"\u003e\u003c/span\u003e\u003c/span\u003e\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eSanchez-Garcia R et al (2021) DeepEMhancer: a deep learning solution for cryo-EM volume post-processing. Commun Biol 4:874. \u003cspan class=\"ExternalRef\"\u003e\u003cspan class=\"RefSource\"\u003ehttps://doi.org/10.1038/s42003-021-02399-1\u003c/span\u003e\u003cspan address=\"10.1038/s42003-021-02399-1\" targettype=\"DOI\" class=\"RefTarget\"\u003e\u003c/span\u003e\u003c/span\u003e\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eBepler T et al (2019) Positive-unlabeled convolutional neural networks for particle picking in cryo-electron micrographs. Nat Methods 16:1153\u0026ndash;1160. \u003cspan class=\"ExternalRef\"\u003e\u003cspan class=\"RefSource\"\u003ehttps://doi.org/10.1038/s41592-019-0575-8\u003c/span\u003e\u003cspan address=\"10.1038/s41592-019-0575-8\" targettype=\"DOI\" class=\"RefTarget\"\u003e\u003c/span\u003e\u003c/span\u003e\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003ePunjani A, Fleet D (2023) 3D Flexible Refinement: Determining Structure and Motion of Flexible Proteins from Cryo-EM. Microsc Microanal 29:1024. \u003cspan class=\"ExternalRef\"\u003e\u003cspan class=\"RefSource\"\u003ehttps://doi.org/10.1093/micmic/ozad067.518\u003c/span\u003e\u003cspan address=\"10.1093/micmic/ozad067.518\" targettype=\"DOI\" class=\"RefTarget\"\u003e\u003c/span\u003e\u003c/span\u003e\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003ePettersen EF et al (2004) UCSF Chimera - A visualization system for exploratory research and analysis. J Comput Chem 25:1605\u0026ndash;1612. \u003cspan class=\"ExternalRef\"\u003e\u003cspan class=\"RefSource\"\u003ehttps://doi.org/10.1002/jcc.20084\u003c/span\u003e\u003cspan address=\"10.1002/jcc.20084\" targettype=\"DOI\" class=\"RefTarget\"\u003e\u003c/span\u003e\u003c/span\u003e\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eJumper J et al (2021) Highly accurate protein structure prediction with AlphaFold. Nature 596:583\u0026ndash;589. \u003cspan class=\"ExternalRef\"\u003e\u003cspan class=\"RefSource\"\u003ehttps://doi.org/10.1038/s41586-021-03819-2\u003c/span\u003e\u003cspan address=\"10.1038/s41586-021-03819-2\" targettype=\"DOI\" class=\"RefTarget\"\u003e\u003c/span\u003e\u003c/span\u003e\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eRoney JP, Ovchinnikov S (2022) State-of-the-Art Estimation of Protein Model Accuracy Using AlphaFold. Phys Rev Lett 129:238101. \u003cspan class=\"ExternalRef\"\u003e\u003cspan class=\"RefSource\"\u003ehttps://doi.org/10.1103/PhysRevLett.129.238101\u003c/span\u003e\u003cspan address=\"10.1103/PhysRevLett.129.238101\" targettype=\"DOI\" class=\"RefTarget\"\u003e\u003c/span\u003e\u003c/span\u003e\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eAdams PD et al (2010) PHENIX: a comprehensive Python-based system for macromolecular structure solution. Acta Crystallogr D Biol Crystallogr 66:213\u0026ndash;221. \u003cspan class=\"ExternalRef\"\u003e\u003cspan class=\"RefSource\"\u003ehttps://doi.org/S0907444909052925\u003c/span\u003e\u003cspan address=\"https://doi.org/S0907444909052925\" targettype=\"URL\" class=\"RefTarget\"\u003e\u003c/span\u003e\u003c/span\u003e [pii] \u003cspan class=\"ExternalRef\"\u003e\u003cspan class=\"RefSource\"\u003e10.1107/S0907444909052925\u003c/span\u003e\u003cspan address=\"10.1107/S0907444909052925\" targettype=\"DOI\" class=\"RefTarget\"\u003e\u003c/span\u003e\u003c/span\u003e\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eEmsley P, Lohkamp B, Scott WG, Cowtan K (2010) Features and development of Coot. Acta Crystallogr D Biol Crystallogr 66:486\u0026ndash;501. \u003cspan class=\"ExternalRef\"\u003e\u003cspan class=\"RefSource\"\u003ehttps://doi.org/S0907444910007493\u003c/span\u003e\u003cspan address=\"https://doi.org/S0907444910007493\" targettype=\"URL\" class=\"RefTarget\"\u003e\u003c/span\u003e\u003c/span\u003e \u003cspan class=\"ExternalRef\"\u003e\u003cspan class=\"RefSource\"\u003e[pii] 10.1107/S0907444910007493\u003c/span\u003e\u003cspan address=\"[pii] 10.1107/S0907444910007493\" targettype=\"DOI\" class=\"RefTarget\"\u003e\u003c/span\u003e\u003c/span\u003e\u003c/span\u003e\u003c/li\u003e\u003c/ol\u003e"}],"fulltextSource":"","fullText":"","funders":[],"hasAdminPriorityOnWorkflow":false,"hasManuscriptDocX":true,"hasOptedInToPreprint":true,"hasPassedJournalQc":"","hasAnyPriority":true,"hideJournal":false,"highlight":"","institution":"","isAcceptedByJournal":true,"isAuthorSuppliedPdf":false,"isDeskRejected":"","isHiddenFromSearch":false,"isInQc":false,"isInWorkflow":false,"isPdf":false,"isPdfUpToDate":true,"isWithdrawnOrRetracted":false,"journal":{"display":true,"email":"[email protected]","identity":"nature-portfolio","isNatureJournal":true,"hasQc":false,"allowDirectSubmit":false,"externalIdentity":"","sideBox":"","snPcode":"","submissionUrl":"","title":"Nature Portfolio","twitterHandle":"","acdcEnabled":false,"dfaEnabled":false,"editorialSystem":"ejp","reportingPortfolio":"","inReviewEnabled":true,"inReviewRevisionsEnabled":false},"keywords":"NADPH-dependent assimilatory sulfite reductase, hemoflavoprotein, oxidoreductase, cryo-EM","lastPublishedDoi":"10.21203/rs.3.rs-4758050/v1","lastPublishedDoiUrl":"https://doi.org/10.21203/rs.3.rs-4758050/v1","license":{"name":"CC BY 4.0","url":"https://creativecommons.org/licenses/by/4.0/"},"manuscriptAbstract":"\u003cp\u003e\u003cem\u003eEscherichia coli\u003c/em\u003e NADPH-dependent assimilatory sulfite reductase (SiR) fixes sulfur for incorporation into sulfur-containing biomolecules. SiR is composed of two subunits: an NADPH, FMN, and FAD-binding diflavin reductase and an iron siroheme/Fe\u003csub\u003e4\u003c/sub\u003eS\u003csub\u003e4\u003c/sub\u003e cluster-containing oxidase. How they interact has been unknown for over 50 years because SiR is highly flexible, thus has been intransigent for traditional X-ray or cryo-EM structural analysis. A combination of the chameleon plunging system with a fluorinated lipid overcame the challenge of preserving a dimer between the subunits for high-resolution (2.84 Å) cryo-EM analysis. Here, we report the first structure of the reductase/oxidase complex, revealing how they interact in a minimal interface. Further, we determined the structural elements that discriminate between pairing a siroheme-containing oxidase with a diflavin reductase or a ferredoxin partner to channel the six electrons that reduce sulfite to sulfide.\u003c/p\u003e","manuscriptTitle":"Structure of dimerized assimilatory NADPH-dependent sulfite reductase reveals the minimal interface for diflavin reductase binding","msid":"","msnumber":"","nonDraftVersions":[{"code":1,"date":"2024-07-29 06:12:11","doi":"10.21203/rs.3.rs-4758050/v1","editorialEvents":[],"status":"published","journal":{"display":true,"email":"[email protected]","identity":"nature-communications","isNatureJournal":true,"hasQc":false,"allowDirectSubmit":false,"externalIdentity":"NCOMMS","sideBox":"Learn more about [Nature Communications](http://www.nature.com/ncomms/)","snPcode":"","submissionUrl":"https://mts-ncomms.nature.com/","title":"Nature Communications","twitterHandle":"","acdcEnabled":true,"dfaEnabled":true,"editorialSystem":"ejp","reportingPortfolio":"Nature Communications","inReviewEnabled":true,"inReviewRevisionsEnabled":false}}],"origin":"","ownerIdentity":"8beaebd7-f67d-411e-a3a5-50557009de1c","owner":[],"postedDate":"July 29th, 2024","published":true,"recentEditorialEvents":[],"rejectedJournal":[],"revision":"","amendment":"","status":"published-in-journal","subjectAreas":[{"id":35168207,"name":"Biological sciences/Structural biology"},{"id":35168208,"name":"Biological sciences/Structural biology/Electron microscopy/Cryoelectron microscopy"}],"tags":[],"updatedAt":"2025-03-27T07:05:34+00:00","versionOfRecord":{"articleIdentity":"rs-4758050","link":"https://doi.org/10.1038/s41467-025-58037-5","journal":{"identity":"nature-communications","isVorOnly":false,"title":"Nature Communications"},"publishedOn":"2025-03-26 04:00:00","publishedOnDateReadable":"March 26th, 2025"},"versionCreatedAt":"2024-07-29 06:12:11","video":"","vorDoi":"10.1038/s41467-025-58037-5","vorDoiUrl":"https://doi.org/10.1038/s41467-025-58037-5","workflowStages":[]},"version":"v1","identity":"rs-4758050","journalConfig":"researchsquare"},"__N_SSP":true},"page":"/article/[identity]/[[...version]]","query":{"redirect":"/article/rs-4758050","identity":"rs-4758050","version":["v1"]},"buildId":"qtupq5eGEP_6zYnWcrvyt","isFallback":false,"isExperimentalCompile":false,"dynamicIds":[84888],"gssp":true,"scriptLoader":[]}

Text is read by the "Ask this paper" AI Q&A widget below. Extraction quality varies by source — PMC NXML preserves structure cleanly, OA-HTML may include some navigation residue, and OA-PDF can have broken hyphenation. The publisher copy (via DOI) is the canonical version.

My notes (saved in your browser only)

Ask this paper AI returns verbatim quotes from the full text · source: preprint-html

Answers must be backed by verbatim quotes from this paper's full text. Hallucinated quotes are dropped automatically; if no verbatim passage answers the question, we say so. How this works

Citation neighborhood (no data yet)

We don't have any in-corpus citations linked to this paper yet. This is a recent paper (2024) — citers typically take a year or two to land, and the OpenAlex reference graph may still be filling in.

Source provenance

europepmc
last seen: 2026-05-20T01:45:00.602351+00:00