# Redox-controlled reorganization and flavin strain within the ribonucleotide reductase R2b-NrdI complex monitored by serial femtosecond crystallography

Juliane John, Oskar Aurelius, Vivek Srinivas, Patricia Saura, In-Sik Kim, Asmit Bhowmick, Philipp S Simon, Medhanjali Dasgupta, Cindy Pham, Sheraz Gul, et al.

► To cite this version:

Juliane John, Oskar Aurelius, Vivek Srinivas, Patricia Saura, In-Sik Kim, et al. Redox-controlled reorganization and flavin strain within the ribonucleotide reductase R2b-NrdI complex monitored by serial femtosecond crystallography. eLife, 2022, 11, pp.e79226. 10.7554/eLife.79226. hal-04907084

HAL Id: hal-04907084

https://hal.science/hal-04907084v1

Submitted on 22 Jan 2025

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_1.png" /></figure>
<figure><img src="image_2.png" /></figure>

# Redox-controlled reorganization and flavin strain within the ribonucleotide reductase R2b-Nrdl complex monitored by serial femtosecond crystallography

Juliane John¹, Oskar Aurelius¹,², Vivek Srinivas¹, Patricia Saura¹, In-Sik Kim³, Asmit Bhowmick³, Philipp S Simon³, Medhanjali Dasgupta³, Cindy Pham³, Sheraz Gul³, Kyle D Sutherlin³, Pierre Aller⁴,⁵, Agata Butryn⁴,⁵†, Allen M Orville⁴,⁵, Mun Hon Cheah⁶, Shigeki Owada⁷,⁸, Kensuke Tono⁷,⁸, Franklin D Fuller⁹, Alexander Batyuk⁹, Aaron S Brewster³, Nicholas K Sauter³, Vittal K Yachandra³, Junko Yano³, Ville RI Kaila¹, Jan Kern³*, Hugo Lebrette¹*, Martin Högbom¹*

¹Department of Biochemistry and Biophysics, Arrhenius Laboratories for Natural Sciences, Stockholm University, Stockholm, Sweden; ²MAX IV Laboratory, Lund University, Lund, Sweden; ³Molecular Biophysics and Integrated Bioimaging Division, Lawrence Berkeley National Laboratory, Berkeley, United States; ⁴Diamond Light Source Ltd, Harwell Science and Innovation Campus, Didcot, United Kingdom; ⁵Research Complex at Harwell, Harwell Science and Innovation Campus, Didcot, United Kingdom; ⁶Department of Chemistry - Ångström, Molecular Biomimetics, Uppsala University, Uppsala, Sweden; ⁷Japan Synchrotron Radiation Research Institute, Sayo-gun, Japan; ⁸RIKEN SPring-8 Center, Sayo-gun, Japan; ⁹LCLS, SLAC National Accelerator Laboratory, Menlo Park, United States

**Abstract** Redox reactions are central to biochemistry and are both controlled by and induce protein structural changes. Here, we describe structural rearrangements and crosstalk within the *Bacillus cereus* ribonucleotide reductase R2b-Nrdl complex, a di-metal carboxylate-flavoprotein system, as part of the mechanism generating the essential catalytic free radical of the enzyme. Femtosecond crystallography at an X-ray free electron laser was utilized to obtain structures at room temperature in defined redox states without suffering photoreduction. Together with density functional theory calculations, we show that the flavin is under steric strain in the R2b-Nrdl protein complex, likely tuning its redox properties to promote superoxide generation. Moreover, a binding site in close vicinity to the expected flavin O₂ interaction site is observed to be controlled by the redox state of the flavin and linked to the channel proposed to funnel the produced superoxide species from Nrdl to the di-manganese site in protein R2b. These specific features are coupled to further structural changes around the R2b-Nrdl interaction surface. The mechanistic implications for the control of reactive oxygen species and radical generation in protein R2b are discussed.

## Editor's evaluation

This paper reports a fundamental set of new results that are obtained using compelling methods in protein crystallography and related fields to investigate and visualize the complex mechanism of an enzyme. The paper will be of interest to a broad audience in structural biology, biochemistry, and enzymology, providing a detailed mechanism of an important biological system and demonstrating

*For correspondence:
jfkern@lbl.gov (JK);
hugo.lebrette@dbb.su.se (HL);
hogbom@dbb.su.se (MH)

**Present address:**
†Macromolecular Machines
Laboratory, The Francis Crick
Institute, London, United
Kingdom; ‡Laboratoire de
Microbiologie et Génétique
Moléculaires (LMGM), Centre de
Biologie Intégrative (CBI), CNRS,
UPS, Université de Toulouse,
Toulouse, France

**Competing interest:** The authors declare that no competing interests exist.

**Funding:** See page 18

**Preprinted:** 14 April 2022

**Received:** 04 April 2022

**Accepted:** 22 August 2022

**Published:** 09 September 2022

**Reviewing Editor:** Donald
Hamelberg, Georgia State
University, United States

© Copyright John et al. This
article is distributed under the
terms of the Creative Commons
Attribution License, which
permits unrestricted use and
redistribution provided that the
original author and source are
credited.
useful tools. The work is timely and has implications for future investigations of complex biochemical processes.

Introduction

Ribonucleotide reductases (RNRs) are essential metalloenzymes that employ sophisticated radical chemistry to reduce the 2'-OH group of ribonucleotides and thereby produce deoxyribonucleotides (dNTPs), the building blocks of DNA. This reaction is the only known pathway for the *de novo* synthesis of dNTPs. Three classes of RNRs are differentiated based on their structural features and radical generating mechanism. Class I consists of a radical generating subunit, R2 and a catalytic subunit, R1. Oxygen is required to generate a radical in the R2 subunit, which is reversibly shuttled to R1 to initiate the ribonucleotide reduction. Class I, found in eubacteria and all eukaryotes, is to date divided into five subclasses, Ia-le, based on the type of metal cofactor, metal ligands, and radical storage in R2 (Högbom et al., 2020). The R2 subunit is characterized by a ferritin-like fold housing two metal ions coordinated by six conserved residues (Nordlund and Reichard, 2006), that is two histidines and four carboxylates (with the exception of subclass Ie, Blaesi et al., 2018; Srinivas et al., 2018). The class Ib R2 subunit (R2b) metal site can bind two manganese ions or two iron ions. The two metal ions are oxidized from the M(II)/M(II) state to a short-lived M(III)/M(IV) intermediate which decays to M(III)/M(III) while producing a radical species on an adjacent tyrosyl residue (Tyr-), where it is also stored (Cotruvo et al., 2013; Cotruvo and Stubbe, 2010). The di-manganese R2b was shown to be the physiologically relevant form of R2b (Cotruvo and Stubbe, 2010; Cox et al., 2010), it accumulates higher amounts of radical and has higher enzymatic activity than the di-iron R2b. Molecular oxygen (O₂) can directly activate the di-iron R2b (Huque et al., 2000). In contrast, the di-manganese R2b form cannot react with O₂ but requires a superoxide radical (O₂⁻) provided by NrdI, a flavoprotein which is generally encoded in the same operon as R2b (Figure 1; Berggren et al., 2014; Cotruvo and Stubbe, 2010; Roca et al., 2008).

NrdI is a small, globular protein that binds the redox-active flavin mononucleotide (FMN) cofactor. NrdI stands out in comparison to other flavodoxins, which perform single one- or two-electron transfers, by being able to perform two successive one-electron reductions (*Cotruvo and Stubbe, 2008*). Fully reduced NrdI with hydroquinone FMN (NrdIhq) can reduce O₂ to O₂⁻ while being oxidized to semiquinone NrdI (NrdIsq) (*Cotruvo et al., 2013*). NrdIsq can in turn produce a second O₂⁻ from another O₂ molecule and become fully oxidized (NrdIox) in the process (*Berggren et al., 2014*). The oxidation state of NrdI can be followed by ultraviolet-visible light absorption (UV-vis) spectroscopy. NrdI is faint yellow in the hydroquinone state, dark blue in the semiquinone state, and bright orange in the oxidized state (*Røhr et al., 2010; Figure 1, Figure 2—figure supplement 1*). NrdI binds tightly to R2b in a 1:1 ratio with FMN at the protein-protein interface, forming a dimer of heterodimers. Previously obtained crystal structures of the R2b-NrdI complex from *Escherichia coli* and *Bacillus cereus* show that a conserved channel connects FMN and the metal site in R2b (Boal et al., 2010;

<figure><img src="image_3.png" /><figcaption>**Figure 1.** Activation of the di-manganese centre in ribonucleotide reductase class Ib R2. Both hydroquinone (hq) and semiquinone (sq) flavin mononucleotide (FMN) of NrdI can reduce molecular oxygen to superoxide, which is shuttled to the metal site in R2b and activates the di-manganese cofactor. FMN is oxidized either to semiquinone or the fully oxidized form (ox) in the process.</figcaption></figure>
Hammerstad et al., 2014). This channel is proposed to shuttle O₂⁻ produced by NrdI in complex with R2b to the metal site, where the radical is formed. Here, we investigate the structural details and reaction mechanism of the activation of NrdI by O₂ as well as the subsequent shuttling of O₂⁻ to the metal site for the generation of Tyr· in R2b.

Studying structural changes in redox-active proteins is challenging from an experimental point of view. X-ray crystallography is a well-suited method for investigating both high-resolution structural details in the redox-active centres and overall reorganization of proteins. However, both metalloenzymes and flavodoxins are very sensitive to photoreduction during the exposure to the high-energy X-ray radiation of synchrotron sources. Synchrotron crystal structures of RNR and NrdI proteins invariably suffer from photoreduction, which complicates the determination of the oxidation state of redox-active centres and has proven to be a problem for obtaining fully oxidized structures (Grave et al., 2019; Johansson et al., 2010; Røhr et al., 2010). While there are well-established experimental workarounds like short exposure times, helical data collection, and serial synchrotron crystallography the time-scales necessary to obtain the diffraction data using these methods cannot fully avoid the reducing effects of the X-ray radiation (Spence, 2017). Serial femtosecond crystallography (SFX) overcomes this problem by the so-called diffraction-before-destruction principle (Doerr, 2011; Nass, 2019; Neutze et al., 2000). Short intense X-ray pulses of the duration of only a few femtoseconds are produced by an X-ray free electron laser (XFEL), illuminating one crystal and giving rise to one diffraction pattern at a time. The crystal is generally destroyed in the process and needs to be replaced by a new one for the next image. This method allows to record the scattering information before atomic displacement has time to occur resulting in a dataset effectively free from the effects of radiation damage (Spence, 2017).

Here, we present the first SFX structures of R2b in complex with NrdI in the oxidized and hydroquinone state. The datasets reveal redox-dependent structural rearrangements both in the FMN-binding pocket of NrdI and at the protein interface. Despite the significant rearrangement in the direct vicinity of the cofactor, FMN itself changes surprisingly little between redox states. This marks an interesting contrast to structures of free NrdI, which display significant conformational change of FMN between different oxidation states. We conclude that the R2b-NrdI complex formation is restricting FMN movement and inhibits this conformational change with possible implications for its redox properties and O₂ reactivity. Moreover, we performed quantum chemical density functional theory (DFT) calculations to computationally estimate the strain energy and additional effects on the FMN redox properties. We also describe the first R2b-NrdI complex with a di-manganese metal centre from B. cereus. The metal coordination in both structures blocks access to the channel that connects FMN and the metal site providing further information on gating and control of catalytic reactive oxygen species.

## Results

### SFX crystal structures of *B. cereus* R2b in complex with oxidized and hydroquinone NrdI collected at ambient temperature

An initial crystallization condition for the R2b-NrdI complex was optimized to yield a sufficient amount of crystals smaller than 100 µm in the longest axis (Figure 2—figure supplement 1), the maximum acceptable crystal size for the experimental setup. Crystals were initially tested at SACLA (SPring-8 Angstrom Compact free electron LAser, Japan) under aerobic conditions with a grease extruder setup (data not shown) (Sugahara et al., 2015; Tono et al., 2013). The crystals diffracted to 2 Å and proved to be stable under room temperature for several days and sturdy enough to handle the physical stress of being manipulated for the experiment. At LCLS (Linac Coherent Light Source at SLAC National Accelerator Laboratory, USA), we obtained two structures of the *B. cereus* R2b-NrdI (*BcR2b-NrdI*) complex in different defined redox states by SFX. The datasets were collected at room temperature under anaerobic conditions using the drop-on-demand sample delivery method (Fuller et al., 2017) (see Materials and methods for details).

*BcR2b* protein was produced metal free to allow full control over metal loading during complex reconstitution. The metal content of the protein was determined by total reflection X-ray fluorescence (TXRF) and only trace amounts of metals could be detected. The iron and manganese content per R2b monomer corresponded to metal-to-protein molar ratios of 0.27% ± 0.04% and 0.07% ± 0.04%, respectively. The *BcR2b-NrdI* complex was reconstituted *in vitro* by mixing both proteins in a molar
<table><caption>**Table 1.** Data collection and refinement statistics.
Values in parenthesis are for the highest resolution shell.</caption><thead><tr><th></th><th>BcR2b<sub>MnMn</sub>-NrdI<sub>ox</sub></th><th>BcR2b<sub>MnMn</sub>-NrdI<sub>hq</sub></th></tr></thead><tbody><tr><td>PDB ID</td><td>7Z3D</td><td>7Z3E</td></tr><tr><td colspan="3"><strong>Data collection statistics</strong></td></tr><tr><td>XFEL source</td><td>LCLS MFX</td><td>LCLS MFX</td></tr><tr><td>Wavelength (Å)</td><td>1.30</td><td>1.30</td></tr><tr><td>Space group</td><td>C222<sub>1</sub></td><td>C222<sub>1</sub></td></tr><tr><td>Unit cell dimensions a, b, c (Å)</td><td>61.4, 125.6, 145.0</td><td>61.2, 125.8, 144.8</td></tr><tr><td>Unit cell angles α, β, γ (°)</td><td>90, 90, 90</td><td>90, 90, 90</td></tr><tr><td>Resolution range (Å)</td><td>51.59–2.0 (2.034–2.0)</td><td>51.47–2.0 (2.034–2.0)</td></tr><tr><td>Unique reflections</td><td>38,319 (1881)</td><td>38,188 (1880)</td></tr><tr><td>Multiplicity</td><td>66.94 (25.42)</td><td>59.92 (26.46)</td></tr><tr><td>Merged lattices</td><td>14,357</td><td>13,509</td></tr><tr><td>Completeness (%)</td><td>99.91 (100)</td><td>99.91 (100)</td></tr><tr><td>Mean I/sigma (I)</td><td>3.247 (0.654)</td><td>3.28 (0.791)</td></tr><tr><td>Wilson B-factor (Å<sup>2</sup>)</td><td>35.45</td><td>33.22</td></tr><tr><td>R<sub>split</sub> (%)</td><td>11.1 (86.3)</td><td>11.1 (76.5)</td></tr><tr><td>CC<sub>1/2</sub></td><td>0.989 (0.32)</td><td>0.988 (0.39)</td></tr><tr><td colspan="3"><strong>Refinement statistics</strong></td></tr><tr><td>Resolution range used in refinement (Å)</td><td>24.17–2.0 (2.07–2.0)</td><td>23.91–2.0 (2.07–2.0)</td></tr><tr><td>Reflections used in refinement</td><td>38,253 (3745)</td><td>38,119 (3729)</td></tr><tr><td>Reflections used for R<sub>free</sub></td><td>1909 (206)</td><td>1897 (203)</td></tr><tr><td>R<sub>work</sub> (%)</td><td>15.91 (30.40)</td><td>15.16 (28.44)</td></tr><tr><td>R<sub>free</sub> (%)</td><td>19.56 (32.27)</td><td>18.37 (31.18)</td></tr><tr><td>RMSD, bond distances (Å)</td><td>0.007</td><td>0.007</td></tr><tr><td>RMSD, bond angles (°)</td><td>0.79</td><td>0.81</td></tr><tr><td>Ramachandran favoured (%)</td><td>99.27</td><td>98.79</td></tr><tr><td>Ramachandran allowed (%)</td><td>0.73</td><td>1.21</td></tr><tr><td>Ramachandran outliers (%)</td><td>0</td><td>0</td></tr><tr><td>Rotamer outliers (%)</td><td>0.79</td><td>1.05</td></tr><tr><td>Clashscore</td><td>1.86</td><td>2.58</td></tr><tr><td>Protein residues (R2b + NrdI)</td><td>415 (299 + 116)</td><td>416 (298 + 118)</td></tr><tr><td>Average B-factor (Å<sup>2</sup>)</td><td>47.27</td><td>43.66</td></tr><tr><td>Macromolecules</td><td>47.24</td><td>43.64</td></tr><tr><td>Ligands</td><td>34.66</td><td>30.38</td></tr><tr><td>Solvent</td><td>50.57</td><td>47.09</td></tr><tr><td>Number of non-H atoms</td><td>3655</td><td>3652</td></tr><tr><td>Macromolecules</td><td>3462</td><td>3477</td></tr><tr><td>Ligands</td><td>33</td><td>33</td></tr><tr><td>Solvent</td><td>160</td><td>142</td></tr></tbody></table>
ratio of 1:1 and set up for batch crystallization. Manganese was present both during the complex reconstitution and in the crystallization condition. Crystals for two different datasets were prepared to investigate different oxidation states of Nrdl. The crystals for the first dataset were grown under aerobic conditions, yielding bright orange crystals, indicating that Nrdl was fully oxidized (Figure 2—figure supplement 1). The structure from these crystals is later referred to as *BcR2b<sub>MnMn</sub>-Nrdl<sub>ox</sub>* (PDB ID: 7Z3D). The second dataset was obtained by reducing *BcR2b<sub>MnMn</sub>-Nrdl<sub>ox</sub>* crystals chemically with sodium dithionite in an anaerobic environment. Consequently, the crystal colour changed from bright orange to faint yellow, indicating that Nrdl underwent a two-electron reduction to the hydroquinone state (Figure 2—figure supplement 1). The crystals were subsequently kept under anaerobic conditions, preventing reoxidation of Nrdl. The corresponding structure is denoted *BcR2b<sub>MnMn</sub>-Nrdl<sub>hq</sub>* (PDB ID: 7Z3E).

## Overall structure

The two structures were solved in space group C222₁ and are of very similar quality with a resolution of 2.0 Å and similar unit cell dimensions (Table 1). The asymmetric unit contains one monomer of the 1:1 *BcR2b-Nrdl* complex. The physiological dimer can be generated by applying crystal symmetry (Figure 2A). Clear electron density maps allowed us to model residues 1–299 of 322 for R2b and 1–116 of 118 for Nrdl in *BcR2b<sub>MnMn</sub>-Nrdl<sub>ox</sub>* and 1–298 of R2b and all 118 residues of Nrdl in *BcR2b<sub>MnMn</sub>-Nrdl<sub>hq</sub>*. Overall *BcR2b<sub>MnMn</sub>-Nrdl<sub>ox</sub>* and *BcR2b<sub>MnMn</sub>-Nrdl<sub>hq</sub>* are similar, as indicated by the Cα root-mean-square deviation (RMSD) value of 0.16 Å (Figure 2B). Refinement of both structures was conducted independently from each other following the same protocol (see Materials and methods for details) with final *R<sub>work</sub>*/*R<sub>free</sub>* of 0.16/0.20 for *BcR2b<sub>MnMn</sub>-Nrdl<sub>ox</sub>* and 0.15/0.18 for *BcR2b<sub>MnMn</sub>-Nrdl<sub>hq</sub>* (Table 1). Hammerstad et al. previously reported two crystal structures of the *BcR2b-Nrdl* complex (Hammerstad et al., 2014), obtained by single-crystal X-ray diffraction at a synchrotron source and under cryogenic conditions. Here, however, R2b harbours a di-iron-active site. These structures will be referred to as *BcR2b<sub>FeFe</sub>-Nrdl-1* and *BcR2b<sub>FeFe</sub>-Nrdl-2* (PDB ID: 4BMO and 4BMP, respectively). Structural alignment between the published synchrotron and our XFEL structures shows that the overall fold of the complex is similar with Cα-RMSD values between 0.42 and 0.47 Å (Figure 2B). Compared to the

<figure><img src="image_4.png" /><figcaption>**Figure 2.** Structure of the *BcR2b-Nrdl* complex. (A) Structure of the *BcR2b<sub>MnMn</sub>-Nrdl<sub>ox</sub>* complex (PDB ID: 7Z3D). In each 1:1 *R2b-Nrdl* dimer, *R2b* and *Nrdl* are coloured in lighter and darker grey, respectively. The second dimer is generated by crystal symmetry. Manganese ions are represented as purple spheres and flavin mononucleotide as yellow sticks at the *R2b-Nrdl* interface. (B) Superimposition of all structures of the *BcR2b-Nrdl* dimer with their corresponding Cα-RMSD in the table. Individual structures are colour coded identically in the figure and the table. Visible N-termini are marked with an arrow (►) and C-termini with a square (■).</figcaption></figure>

The online version of this article includes the following figure supplement(s) for figure 2:

**Figure supplement 1.** Crystallization of the di-manganese *BcR2b-Nrdl* complex.

**Figure supplement 2.** Extended C-terminus of *R2b* in *BcR2b<sub>MnMn</sub>-Nrdl<sub>ox</sub>*.
$BcR2b_{FeFe}$–Nrdl structures we observe a 7–8 residue extended ordered C-terminus for R2b in our structures. The C-terminus in $BcR2b_{MnMn}$–Nrdl$_{ox}$ continues in the same orientation as the C-terminal helix without keeping a helical conformation. It interacts with a groove between R2b and Nrdl extending the R2b–Nrdl-binding area by forming hydrogen bonds to two other helices of R2b and a loop in Nrdl (Figure 2—figure supplement 2).

## R2b complex formation prevents butterfly bend of FMN in Nrdl

Three redox states of FMN are physiologically relevant for Nrdl: oxidized (FMN$_{ox}$), neutral semiquinone (FMN$_{sq}$), and anionic hydroquinone (FMN$_{hq}$) (Cotruvo et al., 2013; Røhr et al., 2010; Figure 3). Reduction of free FMN causes the isoalloxazine to bend along the N5–N10 axis leading to a ‘butterfly bend’ where the isoalloxazine moiety deviates from planarity (Figure 4; Zheng and Ornstein, 1996).

Nrdl binds one molecule of FMN noncovalently. The pyrimidine ring of the 7,8-dimethyl-isoalloxazine ring is buried and tightly bound to the protein by a combination of hydrogen bonds and π-stacking interactions while the benzene ring and part of the phosphate tail are solvent exposed (Figure 4C). The crystal structure of $BcNrdl$ has been described previously by Røhr et al., 2010 and shows the FMN-binding pocket on the protein surface. The same study also presents quantum mechanics/molecular mechanics (QM/MM) calculations of the theoretical butterfly bend for FMN bound to Nrdl. The calculated angles are shallow with -3.1° for FMN$_{ox}$ and -3.9° for FMN$_{sq}$ and more pronounced with 14.1° for FMN$_{hq}$ (Figure 4B). The authors investigated the influence of photoreduction on the butterfly bend in FMN bound to $BcNrdl$. The crystals for the first structure (PDB ID: 2X2O, 1.1 Å resolution) were produced from oxidized protein. The measured butterfly bend in this structure is 4.6° and does thus not correspond well to the calculated FMN$_{ox}$ or FMN$_{sq}$ angle. It does however resemble the calculated theoretical butterfly bend of 3.3° of the physiologically not relevant anionic semiquinone state of FMN (FMN$_{sq}$-) (Figure 4B). The second Nrdl structure (PDB ID: 2X2P, 1.2 Å resolution) structure was produced from protein in the semiquinone state and the butterfly bend was 11° after data collection, similar to the calculated FMN$_{hq}$ angle (Figure 4B). The authors could confirm the reducing effect of the synchrotron radiation on FMN by comparing Raman spectra of the crystals before and after data collection. These structures will be referred to as $BcNrdl_{ox+e}$ and $BcNrdl_{sq+e}$ to emphasize the photoreduction. Johansson et al. observed the same discrepancy between the butterfly bend for a synchrotron structure of initially oxidized Nrdl from *Bacillus anthracis* (BaNrdl), a protein with 99% sequence identity to $BcNrdl$ and identical FMN protein environment (Johansson et al., 2010). BaNrdl$_{ox+e}$ (PDB ID: 2XOD, 1.0 Å resolution) is even more bent than $BcNrdl_{ox+e}$ with 5.7°, indicating significant photoreduction of FMN during data collection (Figure 4—figure supplement 1). Figure 4B lists the FMN angles of the calculated and measured structures.

The use of SFX allowed us, for the first time, to investigate the conformation of oxidized FMN in Nrdl. In $BcR2b_{MnMn}$–Nrdl$_{ox}$, the isoalloxazine moiety of FMN$_{ox}$ is almost planar, with a butterfly bend of 0.2°. This value is comparable to the angle of -3.1° for FMN$_{ox}$ in Nrdl calculated by QM/MM (Røhr et al., 2010). Unexpectedly, the bend of FMN$_{hq}$ in the $BcR2b_{MnMn}$–Nrdl$_{hq}$ structure is minimal with only 0.1° and thus similar to FMN$_{ox}$ in the $BcR2b_{MnMn}$–Nrdl$_{ox}$ structure (Figure 4A, B). Additionally, we calculated an isomorphous difference (Fo(ox)–Fo(hq)) map of $BcR2b_{MnMn}$–Nrdl$_{ox}$ and $BcR2b_{MnMn}$–Nrdl$_{hq}$ using the phases of the oxidized dataset. These maps reduce model bias by directly comparing the experimental data and are very sensitive to subtle changes of atom positions in different datasets (Rould and Carter, 2003). The Fo(ox)–Fo(hq) map shows a slight movement of the oxygen on C4 of FMN indicating a small twist of the pyrimidine ring between both structures but no further bending (Figure 5). This angle of $BcR2b_{MnMn}$–Nrdl$_{hq}$ does not correspond to either the calculated bend of

<figure><img src="image_5.png" /><figcaption>**Figure 3.** The physiologically relevant oxidation states of flavin mononucleotide (FMN) in Nrdl. The red dashed line marks the virtual axis between N5 and N10. The ribityl phosphate group is denoted as R.</figcaption></figure>
<figure><img src="image_6.png" /><figcaption>**Figure 4.** The butterfly bending conformations of flavin mononucleotide (FMN) in different redox states. (A) Overlay of FMN cofactors of *B. cereus* NrdI in different crystal structures: *BcR2b*<sub>MnMn</sub>-NrdI<sub>ox</sub> (dark blue), *BcR2b*<sub>MnMn</sub>-NrdI<sub>ox</sub> (cyan), *BcNrdI*<sub>sq+e</sub> (dark grey, PDB ID: 2X2P), *BcNrdI*<sub>ox+e</sub> (white, PDB ID: 2X2O). (B) Comparison of experimentally measured and computationally calculated FMN bending angles for *BcNrdI* alone and in complex with *BcR2b*. 2Fo-Fc maps are contoured at 2 σ for synchrotron and serial femtosecond crystallography (SFX) structures, the models used are the same as in panel A and color-coded identically. Theoretical butterfly angles obtained by quantum mechanics/molecular mechanics (QM/MM) calculations by *Rohr et al., 2010* and compared *Figure 4 continued on next page*</figcaption></figure>
Figure 4 continued

to angles obtained by density functional theory (DFT) calculations in this work. (C) The benzene ring of FMN in $BcNrdl_{sq+e}$ (PDB ID: 2X2P) is solvent exposed. (D) In $BcR2b_{MnMn}-Nrdl_{hq}$, the FMN-binding pocket is closed by residues contributed by R2b. Relevant residues are shown in sticks. The surface of the FMN-binding pocket in Nrdl in panels (C) and (D) is shown in grey.

The online version of this article includes the following figure supplement(s) for figure 4:

**Figure supplement 1.** Flavin mononucleotide (FMN) in *Bacillus anthracis* $Nrdl_{ox+e}$.

**Figure supplement 2.** Density functional theory (DFT) calculations.

FMN of 14.1° nor the experimental bending angle of $BcNrdl_{sq+e}$ of 11°. In the R2b-Nrdl complex, the solvent exposed part of FMN is covered by R2b, which closes the binding pocket with residues Met16, Phe17, Gln20, Ile197, and Val200 (Figure 4D). In particular, nonpolar contacts formed between the outermost methyl groups of FMN and the side chains of Phe17 and Gln20 are sterically hindering the isoalloxazine moiety from bending. Both side chains are in turn held in place through a network of interactions involving a second shell of residues from R2b and Nrdl that prevent Phe17 and Gln20 from changing conformation should FMN bend (Supplementary file 1).

In order to quantify the strain effects induced by R2b on FMN, we performed quantum chemical DFT calculations based on the R2b-Nrdl complex structures or structures of Nrdl alone (PDB IDs: 2X2O, 2X2P), including all first and second sphere interactions between FMN and the surrounding protein and water molecules, with around 170–250 atoms modelled quantum mechanically (Figure 4—figure supplement 2). The DFT calculations suggest that the isoalloxazine ring of FMN is bent in Nrdl, with bending angles of –2° to 11° depending on the redox state, which compares well with the angles calculated at the QM/MM level by Røhr et al. (Figure 4B, Supplementary file 2). In contrast, we find that the binding of R2b to Nrdl results in a more planar isoalloxazine ring for all oxidation states, with a bending angle of only –3° to +4°, which agrees with the angles obtained experimentally by SFX (Figure 4B) although not fully reproducing the planarity observed in the SFX crystal structures. This geometric distortion introduces molecular strain on FMN with strain energies of around 70–150 meV higher in the R2b-Nrdl complex relative to Nrdl alone (Supplementary file 2). We find that these effects tune the redox potential of FMN by 118 mV and 54 mV in the R2b-Nrdl complex relative to Nrdl alone for the oxidized/semiquinone ($FMN_{sq}/FMN_{ox}$) and semiquinone/hydroquinone ($FMN_{hq}/FMN_{sq}$) redox couples, respectively (Supplementary file 3). These rather large redox tuning effects arise from the combination of strain and electrostatic effects.

## Reorganization of FMN environment and binding position controlled by FMN redox state

The change of redox state of FMN causes movement around the cofactor in both proteins. Alignment of the Nrdl backbone of the $BcR2b_{MnMn}-Nrdl_{ox}$ and the $BcR2b_{MnMn}-Nrdl_{hq}$ structures shows a clear reorganization in the 40 s loop close to the isoalloxazine ring of FMN (Figure 5A, Figure 5—figure supplement 1). The biggest change can be seen for Cα of Gly44 and a smaller shift for Cα of Thr43 and Phe45. The reorganization of the 40 s loop of Nrdl after reduction is also apparent in the superposition of both structures and in the Fo(ox)-Fo(hq) map (Figure 5B, C). Notably, Gly44 of Nrdl is flipped by 180° and forms a hydrogen bond with the hydrogen of N5 on FMN. The hydrogen bond between O4 of FMN and Gly46 shifts slightly. Thr43 in turn is shifted towards FMN and its side chain is rotated by 180°; Phe45 is also slightly rotated (Figure 5C). This redox-dependent conformational change of residues 43–45 was previously also observed in crystal structures of $BcNrdl$ and $BaNrdl$ in the absence of the R2b subunit (PDB ID: 2X2O, 2X2P, 2XOD, 2XOE) (Johansson et al., 2010; Røhr et al., 2010).

Redox-induced changes are also visible on the side of the isoalloxazine ring facing R2b. An unknown molecule is bound within 5 Å from the reactive C4a of FMN in $BcR2b_{MnMn}-Nrdl_{ox}$, which binds O2 to reduce it to $O_2^{\cdot-}$ (Ghislà and Massey, 1989). The molecule is 3-coordinated to Trp74 (Nrdl), Ser196 (R2b), and a water (p-ox in Figure 5D, E). A similar density was previously observed in the di-iron $BcR2b-Nrdl$ synchrotron structures where it was modelled as a chloride ion. In our structure, placing a water in this position does not sufficiently explain the density while modelling it as a chloride ion leaves no residual density (Figure 5—figure supplement 2). Chloride is abundant in the crystallization condition; however, serine and tryptophan are not typical chloride ligands (Carugo, 2014).
<figure><img src="image_7.png" /><figcaption>**Figure 5.** Flavin mononucleotide (FMN) environment in the oxidized and hydroquinone *BcR2b*-Nrdl complex. (A) Structural alignment of Nrdl in *BcR2b*<sub>MnMn</sub>-Nrdl<sub>ox</sub> compared to *BcR2b*<sub>MnMn</sub>-Nrdl<sub>hq</sub>, Cα-RMSD in Å are shown for each residue of Nrdl between the two models. The secondary structure assignment corresponding to the residues of Nrdl is represented in the cartoon below. Nrdl interactions with FMN are marked with light grey background (B) Fo(ox)-Fo(hq) map contoured at 4.5 σ for the *BcR2b*-Nrdl complex with positive density in green and negative density in red. FMNs are shown in sticks, manganese ions as purple spheres, R2b as transparent cartoon in white and Nrdl in transparent black. Differences between both datasets cluster around the FMN at the R2b/Nrdl interface. The distance between metal site and FMN is around 16 Å, marked with a dashed line. (C-E) Closeups of the difference density around the FMN from different angles including superposition of *BcR2b*<sub>MnMn</sub>-Nrdl<sub>hq</sub> in dark blue and *BcR2b*<sub>MnMn</sub>-Nrdl<sub>ox</sub> in cyan. Fo(ox)-Fo(hq) map contoured at 4.5 σ for all panels. Adjacent secondary structure elements are shown as transparent cartoon in white for R2b and black for Nrdl. (C) Rearrangement of the 40 s loop of Nrdl. Gly44 flips 180° in the hydroquinone state and forms a hydrogen bond with N5 of FMN (shown as dashed line). (D) Rearrangements of side chains on the R2b facing side of FMN. Residues of Nrdl are labelled in black, Ser196 from R2b in light grey. An unexplained density, larger than water, moves between the structures, named p-ox and p-hq and is marked as transparent big spheres; waters are represented as small opaque spheres. (E) Same side chains as in (D) are shown at a different angle facing the FMN. FMN is shown transparent in the background for clarity.</figcaption></figure>

The online version of this article includes the following figure supplement(s) for figure 5:

Figure 5 continued on next page
Figure 5 continued

**Figure supplement 1.** Electron density around the 40 s loop of NrdI in the oxidized and reduced complex structures.

**Figure supplement 2.** Different modelling results for unidentified density close to flavin mononucleotide (FMN).

**Figure supplement 3.** Channel between R2b and NrdI.

The exact nature of the molecule could not be determined and was thus designated as an unknown atom in the final model. Interestingly, the binding position undergoes a redox-dependent switch: In the reduced structure the density is instead found at a position between the oxygen on C2 of FMN<sub>hq</sub> and Asn71 of NrdI, 7 Å from the reactive C4a. Two adjacent water molecules make the unknown molecule 4-coordinated (p-hq in **Figure 5D, E**). To accommodate for the change of position residues Met77, Trp74, Asn71, and Glu98 from NrdI and Ser196 from R2b move in a plane parallel to the FMN isoalloxazine ring (**Figure 5D, E**). The unidentified molecule exchanges position with two different water molecules of the well-ordered solvent network connecting FMN with the active site in R2b. The solvent network is housed by a channel between both proteins lined by the side chains of Ser162, Tyr166, Lys263, and Asn267 and the mainchain of Glu195 and Ser196 (*Hammerstad et al., 2014*). The position of the density in $BcR2b_{MnMn}-NrdI_{ox}$ is close to the channel entrance on the R2b surface marked by Lys263, Asn267, and Ser196 (**Figure 5—figure supplement 3**). Despite the large structural rearrangements around the FMN, the Fo(ox)-Fo(hq) map clearly shows that no redox-induced movement protrudes further down the channel and that the coordination of the metal site is unaffected by the change of NrdI redox state (**Figure 5B and Figure 5—figure supplement 3**).

## Access to the metal site is not gated by FMN oxidation state, complex formation, or manganese binding

Both SFX $BcR2b-NrdI$ structures display a di-manganese metal centre with similar coordination. In $BcR2b_{MnMn}-NrdI_{ox}$, the two metal ions of R2b are refined at full occupancy and the coordination sphere is clearly defined in the electron density map (**Figure 6A**). Mn1 is 4-coordinated by His96, Glu93, Asp62, and Glu195. Mn2 is 5-coordinated by His198, Glu93, Glu195, and Glu161. Notably, Glu195 and Glu161 exhibit two alternative conformations each and alternately coordinate Mn2 in a monodentate or bidentate fashion. Importantly, while Glu195 bridges the two ions, Glu161 coordinates only Mn2. No coordinating waters could be identified (**Figure 6B**). The protein complex used in the crystallization for the $BcR2b_{MnMn}-NrdI_{ox}$ dataset has not undergone any oxidation of the metal site and both manganese ions are in the Mn(II)/Mn(II) oxidation state (see Materials and methods for details). Notably, the reduction of FMN within NrdI occurs at about 16 Å from the metal site and does not affect the metal coordination in $BcR2b-NrdI$ (**Figure 5A**). The Fo(ox)-Fo(hq) map shows slightly lower occupancy for both manganese ions in $BcR2b_{MnMn}-NrdI_{hq}$ (**Figure 6—figure supplement 1**). The loss of metals can be explained by the reduction treatment of the crystals for this dataset and Mn1 was modelled with 80%, Mn2 with 90% occupancy in the reduced structure. Both structures show the same metal-metal distance of 3.8–3.9 Å within experimental error for this resolution (**Figure 6B**).

The metal-ligand Glu161 is likely a key player in the activation of the di-manganese centre as it is proposed to gate the access to the metal site for the oxidant produced by NrdI (*Boal et al., 2010*). Two main conformations have been observed for this ligand: In the ‘closed’ conformation the glutamate interacts only with Mn2 and is proposed to prevent $O_2^{\cdot-}$ from reaching the metal site by obstructing the channel. In the ‘open’ conformation, on the other hand, Glu161 bridges both manganese ions leaving space for three water molecules to connect Mn2 to the water network in the channel between FMN and the metal site. The open conformation of the equivalent glutamate has been observed in structures of di-manganese R2b alone, for example from *B. cereus*, *E. coli*, and *Streptococcus sanguinis* (PDB ID: 4BMU, 3N37, 4N83) (**Figure 6C**) but also in the R2b-NrdI complex of *E. coli* (PDB ID: 3N3A) (**Figure 6D; Boal et al., 2010; Hammerstad et al., 2014; Makhlynets et al., 2014**). Thus, NrdI binding to R2b is not responsible for triggering the open conformation of Glu161 (Glu158 in *E. coli* and Glu157 in *S. sanguinis*). It has been hypothesized that the presence of manganese ions in the active site could cause the glutamate to shift to the open conformation as the closed conformation is typically observed in R2b structures containing a di-iron site, for example in the $BcR2b_{FeFe}-NrdI$ complex (PDB ID: 4BMP) (**Figure 6E; Hammerstad et al., 2014**). Indeed, since the di-iron form of R2b
<figure><img src="image_8.png" /><figcaption>**Figure 6.** Comparison of active site architecture in different R2b structures including the radical harbouring tyrosine. Manganese ions are shown as purple, iron ions as orange, and waters as red spheres. The direction of the connecting channel between R2b and NrdI is indicated by a black arrow and 'C'. (A) BcR2b<sub>MnMn</sub>–NrdI<sub>ox</sub> active site. The 2Fo–Fc map is contoured at 2 σ shown as grey mesh and Fo–Fc is contoured at 4 σ with negative density in red (no positive density is present). (B) Superposition of BcR2b<sub>MnMn</sub>–NrdI<sub>hq</sub> and BcR2b<sub>MnMn</sub>–NrdI<sub>ox</sub>. The metal coordination is identical for both structures. (C) Active site of BcR2b<sub>MnMn</sub> (PDB ID: 4BMU). The metal site of chain A is shown with Glu161 in both the open (opaque) and closed (transparent) conformation. (D) Active site of EcR2b<sub>MnMn</sub>–NrdI<sub>hq</sub> (PDB ID: 3N3A). Glu158 is in the open conformation and two waters are in the position of the Glu161 of BcR2b. (E) Active site of the di-iron BcR2b–NrdI complex (PDB ID: 4BMP). The metal coordination is similar to the manganese containing structures with Glu161 in the closed conformation. An additional water is hydrogen bonded to Asp62 and Tyr100 but is not in interaction with iron. Metal–metal distances are shown as dashed and metal ligands as solid lines; residues in alternative conformations are shown transparent for panels (B)–(E).</figcaption></figure>

The online version of this article includes the following figure supplement(s) for figure 6:

**Figure supplement 1.** Comparison of the active site in di-manganese BcR2b<sub>MnMn</sub>–NrdI structures.

is oxidized by O₂ via a NrdI-independent pathway, the movement of the glutamate is not required for radical generation. However, our SFX BcR2b<sub>MnMn</sub>–NrdI structures harbour a di-manganese site and exhibit a Glu161 in the closed conformation, thereby demonstrating that manganese in the active site is not sufficient to induce the glutamate shift. Finally, our radiation damage free structures also show that the opening of the channel towards the metal site is independent from the FMN oxidation state, as Glu161 is in the closed conformation when di-manganese R2b is in complex with either NrdI<sub>ox</sub> or NrdI<sub>hq</sub> (*Figure 6B*). Altogether, our data show that the glutamate shift of Glu161 is not caused by the presence of manganese ions, R2b–NrdI complex formation, a specific NrdI redox state or a combination of these factors.
Discussion

In this study, we present two structures of the manganese containing R2b-Nrdl complex of *B. cereus* with Nrdl in two different oxidation states. These structures were obtained by room temperature SFX data collection using XFEL radiation in contrast to all previously studied R2b or Nrdl structures, which were obtained by classical single-crystal synchrotron data collection at cryogenic temperatures. Determination of the oxidation state of redox-active enzymes via synchrotron radiation is problematic because exposure to X-rays exerts photoreduction on redox centres. Investigation of artefact-free oxidized enzymes is therefore exceedingly challenging with synchrotron radiation. In the context of class Ib RNRs this affects both the redox-active metal centre of R2b and FMN cofactor of Nrdl. SFX eliminates the effects of photoreduction observed during synchrotron-based data collection. We describe SFX structures of the R2b-Nrdl complex with Nrdl in the oxidized and hydroquinone oxidation state. Both 2.0 Å structures are of similar quality as the previously published synchrotron structures of the same complex with a di-iron centre (Hammerstad et al., 2014) showing that the change of methodology does not affect the quality. This allowed detailed examination of structural reorganization induced by changes in FMN redox state. Notably, the FMN conformation per se changes little between the oxidized and hydroquinone oxidation state in R2b-Nrdl. This is surprising since it was previously shown that FMN is bent in reduced Nrdl when not bound to R2b. The complexation of R2b and Nrdl thus prevents FMN bending and exerts strain on the isoalloxazine ring in the hydroquinone state, which was also reproduced in DFT calculations. Interestingly, the calculations revealed that the oxidized FMN is also under strain in the R2b-Nrdl complex in comparison to Nrdl alone. It could also be confirmed by the calculations that R2b causes a shift in redox potential for both redox couples (FMN<sub>sq</sub>/FMN<sub>ox</sub> and FMN<sub>hq</sub>/FMN<sub>sq</sub>). Other flavoproteins have been shown to tune the redox potential of FMN by forcing it into a specific binding angle (Senda et al., 2009; Walsh and Miller, 2003). A recent study by Sorigué et al. presented an SFX structure of the flavoenzyme fatty acid photodecarboxylase (FAP) in complex with an oxidized flavin cofactor exhibiting a butterfly bending angle of 14° (Sorigué et al., 2021). Using time resolved SFX and complementary approaches, the study showed that during the enzymatic cycle of FAP the conformation of the flavin retains the butterfly bend in the different redox states. The authors conclude that the butterfly bend is enforced by the protein scaffold to promote flavin reduction. Following the same line of reasoning we propose that R2b binding to Nrdl restricts FMN bending and thus, opposite to FAP, changes the FMN redox potential of the hydroquinone to favour FMN oxidation and reduction of molecular oxygen to superoxide. This notion is supported by a study by Cotruvo et al. demonstrating that superoxide production by the *in vitro* R2b-Nrdl<sub>hq</sub> complex of *Bacillus subtilis* is about 40 times faster than the production of superoxide by free Nrdl<sub>hq</sub> under the same experimental conditions (Cotruvo et al., 2013). Our calculations support that the effect of strain on the flavin can contribute to downshift the FMN redox potential, which would thermodynamically favour electron transfer. Interestingly, however, the electrostatic effects in the R2b-Nrdl complex are predicted to counteract this effect. From a kinetic perspective, on the other hand, the strain imposed in the FMN<sub>hq</sub> state by the R2b-Nrdl complex could lower the reorganization energy of electron transfer and thus kinetically favour superoxide generation. The full effects of R2b-Nrdl complex formation on the FMN redox properties and their physico-chemical background are complex and warrant further investigation.

Controlling the superoxide production by complex formation could serve as a mechanism to protect
the cell from production of superoxide by free Nrdl. The binding of R2b may also facilitate access of
O₂ to FMN by forming a hydrophobic binding pocket around the benzene ring of FMN. O₂ is reduced
to superoxide at the C4a atom of the isoalloxazine ring of FMN and we observe a molecule bigger
than a water bound about 5 Å away from the C4a atom in the oxidized R2b-Nrdl complex. From the
experimental setup it seems unlikely that the unidentified molecule is a superoxide ion because FMN
never underwent a redox cycle in this structure. With a concentration of 100 mM in the crystallization
condition chloride is a more likely candidate for the molecule, as also suggested by model refinement
(Figure 5—figure supplement 2). However, refinement with a superoxide ion leaves little residual
density, so the binding could fit a molecule of its size (Figure 5—figure supplement 2). Given the
similar negative charge and the comparable electron density for the superoxide and chloride anions,
the observed density could thus represent a potential binding position of superoxide after its gener-
ation. In addition, the density undergoes a redox-dependent switch of position moving further away
from the reactive C4a (7 Å) towards the channel entrance in the reduced R2b-Nrdl complex. This move
<figure><img src="image_9.png" /><figcaption>**Figure 7.** Proposed mechanism for the assembly of the R2b-Nrdl complex. (A) Free Nrdl<sub>hq</sub> and molecular oxygen react slowly to produce superoxide and Nrdl<sub>sq</sub>. (B) R2b and Nrdl form a complex which imposes strain on the flavin mononucleotide (FMN) bend and tunes its redox properties to favour oxygen reduction. The complex formation also generates a redox state-controlled potential superoxide-binding site close to the reactive carbon of FMN, indicated by the dashed and solid circles. Upon exposure to molecular oxygen, superoxide is generated in the complex and shuttled towards R2b.</figcaption></figure>

is accompanied by a change of the putative coordination sphere. Importantly, in both oxidized and reduced states, the unknown molecule is integrated into the conserved hydrogen-bonding network connecting FMN and the metal site and switches position with a water. Even though the exact nature of the observed electron density is unknown we conclude that the binding properties of these key positions are controlled by the FMN redox state and could mark the route superoxide would travel after being produced by the FMN.

Taken together, our results suggest that the production of superoxide by Nrdl and the radical generation in R2b is elegantly orchestrated by the formation of the R2b-Nrdl complex as well as redox state control of binding positions (Figure 7). R2b-Nrdl binding induces conformational strain in the flavin as well as other effects influencing its redox properties to promote superoxide generation. The interaction surface also provides binding positions that are controlled by the redox state of the flavin, presumably involved in gating of channel and metal site access.

# Materials and methods

## Protein expression

The plasmids containing the genes for both *BcR2b* (pET22b-*Bcr2b*) and *BcNrdl* (pET22b-*Bcnrdl*) were kindly provided by Marta Hammerstad and Kristoffer Andersson (University of Oslo). The protein expression of *BcNrdl* was adapted from *Røhr et al., 2010*. *E. coli* BL21(DE3) (New England Biolabs, Frankfurt am Main, Germany) cells were transformed with the pET22b-*Bcnrdl* plasmid. A preculture of Lysogeny Broth medium (Formedium, Norfolk, UK) was inoculated with a single colony, grown over night at 37°C and 200 rpm. The next day large-scale cultures of 1.6 l Terrific Broth medium (Formedium) per glass bottle supplemented with 100 µg/ml carbenicillin (Alfa Aesar, Kandel, Germany) and 1:10,000 (vol/vol) antifoam 204 (Merck, Darmstadt, Germany) were inoculated with 0.5% (vol/vol) of the preculture. The cells were incubated in a LEX bioreactor (Epiphyte3, Toronto, Canada) at 37°C until an optical density at 600 nm of about 0.8 was reached. The cultures were cooled down to 20°C, and the protein expression was induced with 0.8 mM isopropyl β-D-1-thiogalactopyranoside (IPTG) (Formedium). The cells, harvested after 12–16 hr of expression by centrifugation at 4000 × g for 20 min, formed a dark grey pellet. A volume of 1 l culture formed 12–19 g wet cell pellet. Expression of *BcR2b* was adapted from *Tomter et al., 2008* and the protein was expressed in a similar way as
BcNrdl. To express BcR2b metal-free, ethylenediaminetetraacetic acid (EDTA) (PanReac AppliChem, Darmstadt, Germany) was added to the large-scale cultures shortly before induction with IPTG to a final concentration of 1 mM. A volume of 1 l culture gave about 8 g wet cell pellet. The pellets were flash frozen in liquid nitrogen and stored at -20°C until further use.

## Purification

### BcNrdl

The purification of Nrdl was adapted from Rohr et al., 2010. About 20 g of bacterial cell pellet was resuspended in lysis buffer (100 mM Tris-HCI pH 7.5) supplemented with a tablet of EDTA-free cOmplete Protease Inhibitor S2 Cocktail Tablet (PIC) (Roche, Solna, Sweden) and DNase (PanReac Applichem) was added. Cells were lysed with a sonicator, Sonics VCX130 (Sonics, Newtown, CT), and the soluble fraction separated from cell debris by centrifugation at 40,000 × g for 30 min at 4°C. The Nrdl protein was precipitated by slow addition of ammonium sulphate (NH₄SO₄) to the lysate to a final concentration of 60% (wt/vol) (0.37 g/ml) while stirring at 4°C. The precipitate was pelleted by centrifugation at 20,000 × g for 20 min and 4°C and subsequently solubilised with a minimal volume of size-exclusion chromatography (SEC) buffer (50 mM Tris pH 7.5). The protein was desalted by dialysis overnight at 4°C against dialysis buffer (10 mM Tris pH 7.5). The desalted protein was filtered through a 0.45 µm filter to remove precipitate and loaded onto a Q Sepharose High Performance 5 ml anion exchange column (Cytiva, Uppsala, Sweden), washed with SEC buffer and eluted with a gradient from 0% to 50% elution buffer (50 mM Tris-HCI pH 7.5, 1 M KCl). The fractions containing the desired protein could be identified by the bright orange colour of Nrdl (Figure 2—figure supplement 1). Relevant fractions were pooled, concentrated with a Vivaspin 20 centrifugal concentrator with a 30,000 Da molecular weight cut-off polyethersulfone membrane (Vivaspin 30 k concentrator) (Sartorius, Göttingen, Germany) and injected onto a HiLoad Superdex 75 prep grade size-exclusion column (Cytiva). Orange fractions were analysed for purity by sodium dodecyl sulphate-polyacrylamide gel electrophoresis (SDS-PAGE) and fractions with a 95% or higher purity were pooled. The theoretical molecular weight of Nrdl was calculated to be 13,449 Da with ProtParam (Gasteiger et al., 2005) and the extinction coefficient ε<sub>447</sub> = 10.8 mM<sup>-1</sup> cm<sup>-1</sup> (Berggren et al., 2014) was used to determine the protein concentration using UV-vis spectroscopy. The protein was concentrated to 25 mg/ml, flash frozen in liquid nitrogen and stored at -80°C until further use.

### BcR2b

The R2b purification protocol was adapted from Tomter et al., 2008. The protein was produced metal free; for that purpose, EDTA was included in the lysis buffer to inhibit metal uptake. About 20 g of bacterial cell pellet was resuspended in lysis buffer (100 mM Tris-HCI pH 7.5, 1 mM EDTA) with a tablet of PIC and lysed by sonication. Streptomycin sulphate was added to a total of 2.5% (wt/vol) and incubated for 10 min at 20°C for DNA precipitation. The crude lysate was cleared by centrifugation for 30 min at 40,000 × g. The supernatant was cleared from contaminants by adding 40% (0.24 g/ml) NH₄SO₄ at 20°C followed by centrifugation at 20,000 × g for 20 min. The NH₄SO₄ concentration for the remaining supernatant was increased to 50% (0.31 g/ml) at 20° to precipitate R2b. After the second centrifugation at 20,000 × g for 20 min the pellet was dissolved with a minimal volume of lysis buffer. The protein was diluted 4× with buffer A (50 mM Tris pH 7.5, 1.5 M NH₄SO₄) to increase the salt concentration sufficiently. The diluted protein was filtered through a 0.45-µm membrane filter and loaded onto a column packed with 20 ml Phenyl Sepharose High Performance resin (Cytiva). The protein was washed with buffer A and eluted with a gradient with SEC buffer (50 mM Tris pH 7.5). Several column volumes were needed to elute the protein. The elution fractions were analysed by SDS-PAGE. Protein containing fractions were pooled, concentrated by a Vivaspin 50 k concentrator and injected onto a HiLoad Superdex 200 prep grade column (Cytiva). Elution fractions were analysed by SDS-PAGE and pure fractions pooled. The theoretical molecular weight and the extinction coefficient of R2b were calculated with ProtParam to be 37,017 Da and ε<sub>280</sub> = 48.36 mM<sup>-1</sup> cm<sup>-1</sup>, respectively (Gasteiger et al., 2005). R2b was concentrated to 50 mg/ml with a Vivaspin 50 kDa concentrator, flash frozen and stored at -80°C until further use. A total of 20 g bacterial cell pellet yielded about 100-150 mg pure protein.
## Quantification of R2b metal content with TXRF

To control the metal content of *BcR2b* the protein was measured with TXRF spectroscopy using a Bruker PicoFox S2 spectrometer (Bruker, Billerica, MA). Three independently prepared replicates of the concentrated protein at 1.7 mM were mixed 1:1 with a gallium standard at 20 mg/l and dried on top of a siliconized quartz sample carrier. The discs were individually measured before use to avoid external contamination of the sample. The results were analysed with the Bruker Spectra software version 7.8.2.0 provided with the instrument.

## Crystallization of the *BcR2b*-Nrdl complex

The *BcR2b*-Nrdl protein complex was prepared for crystallization by first metal loading *BcR2b* with manganese(II), yielding the reduced R2b<sub>MnMn</sub> and then forming the complex by adding oxidized Nrdl: 0.25 mM *BcR2b* in 50 mM Tris-HCl pH 7.5 was incubated with 12 molar equivalents of Mn(II)Cl<sub>2</sub> for 10 min at 20° C, then 0.25 mM *BcNrdl* was added and the mixture incubated for another 15 min at 20°C to ensure complex formation. New crystallization conditions were screened since the crystallization condition published in *Hammerstad et al., 2014* could not be reproduced. Initial hits in conditions C5 of the JSCG+ crystallisation screen (0.1 M sodium HEPES pH 7.5, 0.8 M sodium phosphate monobasic monohydrate, 0.8 M potassium phosphate monobasic) and the A7 condition of the PACT premier crystallization screen (0.1 M sodium acetate pH 5.0, 0.2 M NaCl, 20% (wt/vol) PEG 6000) (both Molecular Dimensions, Sheffield, UK) were further optimized to yield crystals between 10 and 50 µm in the longest axis (Figure 2—figure supplement 1). The final crystallization protocol was established as described: Crystals spontaneously grew in a hanging drop vapour diffusion experiment after one to two days at 20°C. Crystallization condition A (0.1 M HEPES pH 7.0, 0.6–0.85 M sodium phosphate monobasic monohydrate, 0.6–0.85 M potassium phosphate monobasic) was manually mixed with the protein complex solution in a ratio of 1:1 (1 + 1 µl). The crystals grown in this experiment were used to produce microseeds for crystallization with the batch method. Two crystal containing drops were transferred into 50 µl of crystallization condition A in the bottom of a 1.5 ml microcentrifuge tube. A seed bead (Saint Gobain, Aachen, Germany) was placed into the solution, the crystal crushed by vigorous shaking and the mixture used as seed stock for the following crystallization experiments. Batch crystallization was set up in PCR tubes at 20°C. A volume of 40 µl of crystallization condition B (0.1 M Mn(II)Cl<sub>2</sub>, 0.1 M sodium acetate 5.0, 5% PEG 6000) were pipetted on top of 40 µl of protein complex solution. A volume of 8 µl of microseed solution was added to the tube and everything mixed by pipetting. Orange rhombohedron shaped crystals typically sized between 20 and 50 µm in the longest axis grew over night, forming an orange pellet at the bottom of the tubes. The crystals were resuspended and pooled in 2-ml microcentrifuge tubes (Figure 2—figure supplement 1). Crystals used to collect the dataset for *BcR2b*<sub>MnMn</sub>–Nrdl<sub>ox</sub> could be directly loaded into a Gastight SampleLock (Hamilton, Bonaduz, Switzerland) syringe for sample delivery. The crystals used for the *BcR2b*<sub>MnMn</sub>–Nrdl<sub>hq</sub> dataset had to be chemically reduced first as described below.

## Reduction protocol

Crystals for the *BcR2b*<sub>MnMn</sub>–Nrdl<sub>hq</sub> dataset were chemically reduced before data collection. The reduction affects only FMN in Nrdl, since the active site already contains Mn(II)/Mn(II) due to metal loading with Mn(II)Cl<sub>2</sub>. All following steps were conducted in an anaerobic glove box with O<sub>2</sub> below 10 ppm at room temperature. A volume of 900 µl of pooled crystal slurry was gently centrifuged, forming a dense crystal pellet. Of the supernatant 800 µl were carefully removed, collected separately, and supplemented with freshly prepared, anaerobic sodium dithionite (DT) to a final concentration of 20 mM. The DT-containing supernatant was gently mixed with the crystal pellet and the colour change of the crystals from a bright orange to a faint, light yellow was observed in a matter of minutes. The DT was subsequently washed out by gently spinning down the crystals, removing the supernatant without disturbing the crystal pellet, replacing it with anaerobic wash buffer (100 mM MnCl<sub>2</sub>, 50 mM sodium acetate pH 5.0, 2.5% (wt/vol) PEG 6000, 10% (vol/vol) glycerol) and gently resuspending the crystals. This washing step was repeated three times, adding wash buffer in the last step to a volume of 800 µl. The final concentration of DT in the crystal slurry was below 4 µM. All of the crystal slurry was loaded into a 1 ml gastight SampleLock Hamilton syringe and stored in the anaerobic box until sample injection.
## Data collection

The *BcR2b*–Nrdl crystals were initially tested for stability and diffraction quality at SACLA in Japan during experiment 2017B8085. The sample was delivered with a grease extruder setup installed on site (Sugahara et al., 2015; Tono et al., 2013). Hydroxyethyl cellulose matrix (Sugahara et al., 2017) and the protein crystals were mixed together in a volume ratio matrix to crystal pellet of 9:1 and the mixture ejected with a HPLC pump through a 150 µm nozzle with a flow rate of 1–1.5 µl/min, delivering the sample into the XFEL beam. X-rays were delivered as <10 fs long pulses at 10.9 keV with 30 Hz repetition rate and a typical pulse energy of around 0.32 mJ with a beam size of 2 × 2 µm² (fwhm). The diffraction of the sample was recorded 100 mm downstream of the interaction point on an Octal MPCCD detector. The initial data collection showed stability of the crystals over days and a diffraction quality to about 2 Å.

The datasets used in this paper were collected at LCLS, California during experiment LU50. X-ray pulses at 9.5 keV with a pulse energy of 4 mJ, 30 Hz repetition rate and a duration of around 35 fs were generated and used for X-ray diffraction in the macromolecular femtosecond crystallography (MFX) experimental hutch (Sierra et al., 2019). The diffraction was recorded on a Rayonix MX340 (Rayonix L.L.C, Evanston, USA) detector. The sample was delivered into the X-ray interaction point with the drop-on-tape method; the detailed method was described by Fuller et al., 2017. Briefly, crystal slurry in a 1-ml SampleLock syringe was pumped with a syringe pump (KD scientific, Holliston, MA) at a flow rate of 8–9 µl/min through a silica capillary into a 6 µl sample reservoir. An acoustic transducer transferred crystal containing droplets of 2.5–4 nl volume onto a Kapton tape, which transported the droplets into the X-ray beam at 28°C and 27% relative humidity with a speed of 300 mm/s, which resulted in the crystals being exposed to the He environment for about 0.8 s. The enclosure of the setup was filled with a He atmosphere with an O₂ level below 0.1% for the *BcR2b*<sub>MnMn</sub>–Nrdl<sub>hq</sub> data collection.

## Processing, structure determination, and refinement

The datasets for *BcR2b*<sub>MnMn</sub>–Nrdl<sub>hq</sub> and *BcR2b*<sub>MnMn</sub>–Nrdl<sub>ox</sub> were processed with cctbx.xfel (Brewster et al., 2019b; Hattne et al., 2014; Sauter, 2015) and DIALS (Brewster et al., 2018; Winter et al., 2018). We performed joint refinement of the crystal models against the detector position for each batch to account for small time-dependent variations in detector position and also corrected for the Kapton tape shadow (Fuller et al., 2017). Data were scaled and merged to 2.0 Å resolution using cctbx.xfel.merge with errors determined by the ev11 method (Brewster et al., 2019a). The final resolution cut-off was based on values of CC<sub>1/2</sub> (where it does not monotonically decrease as a function of resolution anymore) and multiplicity (where it falls off the 10-fold threshold), as well as on *R*-factors after initial refinement. Data statistics are available in Table 1. Both structures were solved by molecular replacement and refined independently with the PHENIX Suite (Liebschner et al., 2019). The phases were solved with phenix.phaser (McCoy et al., 2007). A *BcR2b*–Nrdl complex structure (PDB ID: 4BMO; Hammerstad et al., 2014) was modified by manually removing waters and alternate conformations and used as a starting model. The suggested solutions were in the same space group (C222<sub>1</sub>) as the starting model with similar unit cell dimensions (Table 1) with one complex in the asymmetric unit. The *R*-free set of the 4BMO model corresponding to 5% of reflections was assigned to all datasets with phenix.reflection Tools. Restraints of the entry ‘FMN’ of the ligand database provided by phenix were used for the oxidized FMN in *BcR2b*<sub>MnMn</sub>–Nrdl<sub>ox</sub>. Restraints for the hydroquinone FMN in *BcR2b*<sub>MnMn</sub>–Nrdl<sub>hq</sub> were generated with phenix.eLBOW (Moriarty et al., 2009). The datasets were iteratively refined with phenix.refine (Afonine et al., 2012), examined and built in coot (Emsley et al., 2010) and validated with MolProbity (Williams et al., 2018). Refinement of all atoms included isotropic *B*-factors, TLS parameters, occupancy and reciprocal space refinement with a high-resolution cut-off of 2.0 Å. Waters were initially added using phenix.refine and in later refinements corrected manually. The metal occupancy for the reduced structure was fixed manually after several cycles of refinement. The unknown species found in the vicinity of the FMN was modelled as unknown atom (PDB ligand ID: UNX) in the two structures. An overview over refinement and model quality statistics can be found in Table 1, created with phenix.table_one. The refined structures were compared and RMSD values calculated with SSM superimpose (Krissinel and Henrick, 2004). The contacts of FMN in both structures listed in Supplementary file 1 were analysed with the help of LIGPLOT (Wallace et al., 1995).
## DFT calculations

Quantum chemical DFT models of the FMN-binding site in R2b-Nrdl and Nrdl systems were built based on the coordinates of the crystal structures of R2b-Nrdl<sub>ox</sub> and R2b-Nrdl<sub>hq</sub> (this study), and Nrdl<sub>ox</sub> (PDB ID: 2X2O) and Nrdl<sub>sq</sub> (PDB ID: 2X2P) (*Rehr et al., 2010*). The models of the Nrdl system comprised the FMN isoalloxazine and ribityl moieties, sidechains of residues Trp74, Thr43, Thr42, Ser69, and backbone of Thr42, Thr43, Gly44, Phe45, Gly46, Ser69, Gly70, Asn71, Met77, Phe78, Gly79, and by including five water molecules resolved in the crystal structure, comprising 169-170 atoms (**Figure 4—figure supplement 2A**). When applicable, free N-H or C=O backbone groups were removed from the selection. Model systems of the R2b-Nrdl complex included the same FMN and Nrdl groups as in the Nrdl-only systems, as well as Nrdl Met8 sidechain, and sidechain of R2b residues Met16, Phe17, Gln20, Ile197, Val200, and five water molecules resolved in the structure, yielding a system with 245-246 atoms (**Figure 4—figure supplement 2B**). Protein sidechains were cut at the Cα-Cβ, and saturated with hydrogens. During the geometry optimizations, the Cβ positions, as well as Cα positions of the protein backbone, and terminal parts of the ribityl tail of FMN were fixed to their experimental positions. Additional restraints were also introduced on residues from the R2b subunit to account for the interaction of the binding site. Geometry optimizations were performed at the B3LYP-D3/def2-SVP/ε = 4 level of theory, and single point energies were re-evaluated at the B3LYP-D3/def2-TZVP/ε = 4 level (*Becke, 1993; Grimme et al., 2010; Klamt and Schüürmann, 1993; Lee et al., 1988; Schäfer et al., 1992*). The bending angles of the flavin were defined as the dihedral angle between the plane formed by the pyrimidine moiety and N5, N10 atoms, with the plane formed by the benzene moiety and N5, N10 atoms and measured with UCSF Chimera (*Pettersen et al., 2004*). Strain energies, computed at the B3LYP-D3/def2-TZVP level, were evaluated as the difference between the isolated flavin group optimized in its protein environment and the isolated flavin, reoptimized without the protein environment. The systems were modelled in oxidized state (FMN<sub>ox</sub>), the neutral semiquinone state (FMN<sub>sq</sub>), and anionic hydroquinone state (FMN<sub>hq</sub>). All DFT calculations were performed with TURBOMOLE v. 7.5 (*Balasubramani et al., 2020*).

## Figures

Molecular figures were prepared with the PyMOL Molecular Graphics System, Version 2.4.2 Schrödinger, LLC. The surface representation of the channel in **Figure 5—figure supplement 3** was generated with the 'Cavities and Pockets' function of PyMOL after removing all waters from the model and with the surface cavity radius set to -2 and the surface cavity cut-off to -6. **Figure 3** was designed with ChemDraw (PerkinElmer Informatics). The bar graph in **Figure 5A** was generated with GraphPad Prism, Version 9.2.0 for iOS (GraphPad Software, San Diego, CA, USA). All figures except **Figure 3** were designed and assembled in Affinity Designer (Serif (Europe) Ltd, Nottingham, UK).

## Data availability

## Acknowledgements

We thank Robert Bolotovsky, Iris D Young, and Lee J O'Riordan for help with data processing and the staff at LCLS and SACLA. We thank Kristine Grave for input and helpful discussions during data analysis and manuscript writing. AMO, PA, and ABu were supported by Diamond Light Source, the UK Science and Technology Facilities Council (STFC), a jointly funded strategic award from the Wellcome Trust and the Biotechnology and Biological Sciences Research Council (102593 to James Naismith). The DOT instrument used in this research was funded by Department of Energy (DOE), Office of Science, Office of Basic Energy Sciences (BES), Division of Chemical Sciences, Geosciences, and Biosciences (to JK, JY, and VKY). XFEL data were collected under proposal LU50 at LCLS, SLAC, Stanford, USA, and under proposal 2017B8085 at BL2 of SACLA, Japan. The Rayonix detector used at LCLS was supported by the NIH grant S10 OD023453. Use of the LCLS, SLAC National Accelerator Laboratory, is supported by the U.S. DOE, Office of Science, BES, under contract no. DE-AC02-76SF00515. Computational resources were provided by the Swedish National Infrastructure for Computing (SNIC 2021/1-40, SNIC 2022/1-29).
## Additional information

### Funding

<table><thead><tr><td>Funder</td><td>Grant reference number</td><td>Author</td></tr></thead><tbody><tr><td>Vetenskapsrådet</td><td>2017-04018</td><td>Martin Högbom</td></tr><tr><td>Vetenskapsrådet</td><td>2021-03992</td><td>Martin Högbom</td></tr><tr><td>European Research Council</td><td>HIGH-GEAR 724394</td><td>Martin Högbom</td></tr><tr><td>Knut och Alice Wallenbergs Stiftelse</td><td>2017.0275</td><td>Martin Högbom</td></tr><tr><td>Knut och Alice Wallenbergs Stiftelse</td><td>2019.0436</td><td>Martin Högbom</td></tr><tr><td>Knut och Alice Wallenbergs Stiftelse</td><td>2019.0251</td><td>Ville RI Kaila</td></tr><tr><td>National Institutes of Health</td><td>GM133081</td><td>Kyle D Sutherlin</td></tr><tr><td>National Institutes of Health</td><td>GM117126</td><td>Nicholas K Sauter</td></tr><tr><td>National Institutes of Health</td><td>GM55302</td><td>Vittal K Yachandra</td></tr><tr><td>National Institutes of Health</td><td>GM110501</td><td>Junko Yano</td></tr><tr><td>National Institutes of Health</td><td>GM126289</td><td>Jan Kern</td></tr><tr><td>Wellcome Trust</td><td>Investigator Award 210734/Z/18/Z</td><td>Allen M Orville</td></tr><tr><td>Royal Society</td><td>Wolfson Fellowship RSWF\ R2\182017</td><td>Allen M Orville</td></tr></tbody></table>

The funders had no role in study design, data collection, and interpretation, or the decision to submit the work for publication. For the purpose of Open Access, the authors have applied a CC BY public copyright license to any Author Accepted Manuscript version arising from this submission.

### Author contributions

Juliane John, Conceptualization, Data curation, Formal analysis, Investigation, Visualization, Methodology, Writing – original draft, Writing – review and editing, Developed the study, Produced and purified proteins, Developed the crystallization protocol, Performed LCLS experiment, Refined, analysed and interpreted crystal structures, Performed TXRF measurements; Oskar Aurelius, Conceptualization, Investigation, Methodology, Writing – review and editing, Developed the study, Performed LCLS experiment; Vivek Srinivas, Conceptualization, Investigation, Methodology, Developed the study, Performed SACLA experiment; Patricia Saura, Formal analysis, Investigation, Writing – review and editing, Performed and analysed the DFT calculations; In-Sik Kim, Investigation, Methodology, Developed, tested and ran the sample delivery system, Performed LCLS experiment; Asmit Bhowmick, Formal analysis, Investigation, Writing – review and editing, Processed and analysed XFEL data, Performed LCLS experiment; Philipp S Simon, Investigation, Methodology, Writing – review and editing, Developed, tested and ran the sample delivery system, Performed LCLS experiment; Medhunjali Dasgupta, Formal analysis, Writing – review and editing, Processed and analysed XFEL data; Cindy Pham, Investigation, Methodology, Developed, tested and ran the sample delivery system, Performed LCLS experiment; Sheraz Gul, Investigation, Methodology, Developed, tested and ran the sample delivery system, Performed LCLS experiment; Kyle D Sutherlin, Formal analysis, Funding acquisition, Processed and analysed XFEL data, Performed LCLS experiment; Pierre Aller, Investigation, Writing – review and editing, Performed LCLS and SACLA experiments; Agata Butryn, Formal analysis, Investigation, Writing – review and editing, Processed and analysed XFEL data, Performed LCLS experiment;
Allen M Orville, Funding acquisition, Investigation, Methodology, Writing – review and editing, Developed, tested and ran the sample delivery system, Performed LCLS and SACLA experiments; Mun Hon Cheah, Investigation, Writing – review and editing, Performed SACLA experiment; Shigeki Owada, Investigation, Performed SACLA experiment; Kensuke Tono, Investigation, Performed SACLA experiment; Franklin D Fuller, Investigation, Methodology, Developed, tested and ran the sample delivery system, Operated the MFX instrument at LCLS, Performed LCLS and SACLA experiments; Alexander Batyuk, Investigation, Operated the MFX instrument at LCLS, Performed LCLS experiment; Aaron S Brewster, Formal analysis, Investigation, Processed and analysed XFEL data, Performed LCLS experiment; Nicholas K Sauter, Formal analysis, Funding acquisition, Writing – review and editing, Processed and analysed XFEL data, Performed LCLS experiment; Vittal K Yachandra, Resources, Funding acquisition, Investigation, Performed LCLS and SACLA experiments; Junko Yano, Resources, Funding acquisition, Investigation, Writing – review and editing, Performed LCLS and SACLA experiments; Ville RI Kaila, Formal analysis, Funding acquisition, Investigation, Writing – review and editing, Performed and analysed the DFT calculations; Jan Kern, Resources, Funding acquisition, Investigation, Methodology, Project administration, Writing – review and editing, Developed, tested and ran the sample delivery system, Performed LCLS and SACLA experiments; Hugo Lebrette, Conceptualization, Formal analysis, Supervision, Investigation, Methodology, Writing – original draft, Project administration, Writing – review and editing, Developed the study, Produced and purified proteins, Developed the crystallization protocol, Performed SACLA and LCLS experiments, Interpreted crystal structures; Martin Högbom, Conceptualization, Resources, Supervision, Funding acquisition, Investigation, Methodology, Project administration, Writing – review and editing, Developed the study, Performed LCLS and SACLA experiments, Interpreted crystal structures

## Author ORCIDs

Juliane John http://orcid.org/0000-0001-9626-3670

Vivek Srinivas http://orcid.org/0000-0002-0265-1873

Patricia Saura http://orcid.org/0000-0003-2575-9913

Philipp S Simon http://orcid.org/0000-0002-2859-4475

Agata Butryn http://orcid.org/0000-0002-5227-4770

Allen M Orville http://orcid.org/0000-0002-7803-1777

Mun Hon Cheah http://orcid.org/0000-0001-5732-1524

Shigeki Owada http://orcid.org/0000-0002-1451-7612

Kensuke Tono http://orcid.org/0000-0003-1218-3759

Alexander Batyuk http://orcid.org/0000-0002-9393-2880

Aaron S Brewster http://orcid.org/0000-0002-0908-7822

Nicholas K Sauter http://orcid.org/0000-0003-2786-6552

Junko Yano http://orcid.org/0000-0001-6308-9071

Ville RI Kaila http://orcid.org/0000-0003-4464-6324

Jan Kern http://orcid.org/0000-0002-7272-1603

Hugo Lebrette http://orcid.org/0000-0002-8081-181X

Martin Högbom http://orcid.org/0000-0001-5574-9383

## Decision letter and Author response

Decision letter https://doi.org/10.7554/eLife.79226.sa1

Author response https://doi.org/10.7554/eLife.79226.sa2

## Additional files

### Supplementary files

* Supplementary file 1. List of contacts that prevent flavin mononucleotide (FMN) from bending due to R2b-Nrdl complex formation. The contacts of FMN to Nrdl are not listed since those are also formed in the absence of R2b. The only additional contacts formed with between FMN and R2b are with the side chains (sc) of Phe17 and Gln20. The contacts formed by those side chains are also given. Contacts in both the oxidized and reduced structures are listed.

* Supplementary file 2. Bending angle (in degrees) and strain energies ($E_{strain}$, in eV) of the isoalloxazine ring of flavin mononucleotide (FMN) in BcNrdl based on density functional theory (DFT) calculations. The strain energies are computed relative to the same redox state of the flavin, optimized without the protein surroundings. See Materials and methods for definition of the
bending angle.

* Supplementary file 3. Redox potential shifts of Nrdl and R2b-Nrdl. The relative redox potentials, ΔE<sub>tot</sub>, are computed at the B3LYP-D3/def2-TZVP/ε = 4 level relative to the same redox transition of the isolated flavin in ε = 4. The shifts in strain energy (at ε = 4), ΔE<sub>strain</sub>, refer to the same redox transition relative to the isolated flavin, optimized without surroundings. The protein electrostatic shift, ΔE<sub>el</sub>, are estimated based on the difference between ΔE<sub>tot</sub> and ΔE<sub>strain</sub> (ΔE<sub>tot</sub> = ΔE<sub>strain</sub> + ΔE<sub>el</sub>). ΔΔE<sub>tot</sub> are the total shifts in redox potentials of R2b-Nrdl relative to Nrdl alone.

* MDAR checklist

Data availability

The atomic coordinates and crystallographic data have been deposited in the Protein Data Bank under the following accession codes: 7Z3D and 7Z3E.

The following datasets were generated:

<table><thead><tr><td>Author(s)</td><td>Year</td><td>Dataset title</td><td>Dataset URL</td><td>Database and Identifier</td></tr></thead><tbody><tr><td>John J, Lebrette H, Aurelius O, Högbom M</td><td>2022</td><td>XFEL structure of Class Ib ribonucleotide reductase dimanganese(II) NrdF in complex with oxidized Nrdl from Bacillus cereus</td><td>https://www.rcsb.org/structure/7Z3D</td><td>RCSB Protein Data Bank, 7Z3D</td></tr><tr><td>John J, Lebrette H, Aurelius O, Högbom M</td><td>2022</td><td>XFEL structure of Class Ib ribonucleotide reductase dimanganese(II) NrdF in complex with hydroquinone Nrdl from Bacillus cereus</td><td>https://www.rcsb.org/structure/7Z3E</td><td>RCSB Protein Data Bank, 7Z3E</td></tr></tbody></table>

References

* Afonine PV, Grosse-Kunstleve RW, Echols N, Headd JJ, Moriarty NW, Mustyakimov M, Terwilliger TC, Urzhumtsev A, Zwart PH, Adams PD. 2012. Towards automated crystallographic structure refinement with phenix.refine. *Acta Crystallographica. Section D, Biological Crystallography* **68**:352–367. DOI: https://doi.org/10.1107/S0907444912001308, PMID: 22505256

* Balasubramani SG, Chen GP, Coriani S, Diedenhofen M, Frank MS, Franzke YJ, Furche F, Grotjahn R, Harding ME, Hättig C, Hellweg A, Helmich-Paris B, Holzer C, Huniar U, Kaupp M, Marefat Khah A, Karbalaei Khani S, Müller T, Mack F, Nguyen BD, et al. 2020. TURBOMOLE: modular program suite for ab initio quantum-chemical and condensed-matter simulations. *J Chem Phys* **152**:184107. DOI: https://doi.org/10.1063/5.0004635, PMID: 32414256

* Becke AD. 1993. Density-functional thermochemistry. III. the role of exact exchange. *J Chem Phys* **98**:5648–5652. DOI: https://doi.org/10.1063/1.464913

* Berggren G, Duraffourg N, Sahlin M, Sjöberg BM. 2014. Semiquinone-induced maturation of bacillus anthracis ribonucleotide reductase by a superoxide intermediate. *The Journal of Biological Chemistry* **289**:31940–31949. DOI: https://doi.org/10.1074/jbc.M114.592535, PMID: 25262022

* Blaesi EJ, Palowitch GM, Hu K, Kim AJ, Rose HR, Alapati R, Lougee MG, Kim HJ, Taguchi AT, Tan KO, Laremore TN, Griffin RG, Krebs C, Matthews ML, Silakov A, Bollinger JM, Allen BD, Boal AK. 2018. Metal-free class ie ribonucleotide reductase from pathogens initiates catalysis with a tyrosine-derived dihydroxyphenylalanine radical. *PNAS* **115**:10022–10027. DOI: https://doi.org/10.1073/pnas.1811993115, PMID: 30224458

* Boal AK, Cotruvo JA, Stubbe JA, Rosenzweig AC. 2010. Structural basis for activation of class ib ribonucleotide reductase. *Science* **329**:1526–1530. DOI: https://doi.org/10.1126/science.1190187, PMID: 20688982

* Brewster AS, Waterman DG, Parkhurst JM, Gildea RJ, Young ID, O'Riordan LJ, Yano J, Winter G, Evans G, Sauter NK. 2018. Improving signal strength in serial crystallography with DIALS geometry refinement. *Acta Crystallographica. Section D, Structural Biology* **74**:877–894. DOI: https://doi.org/10.1107/S2059798318009191, PMID: 30198898

* Brewster AS, Bhowmick A, Bolotovsky R, Mendez D, Zwart PH, Sauter NK. 2019a. SAD phasing of XFEL data depends critically on the error model. *Acta Crystallographica. Section D, Structural Biology* **75**:959–968. DOI: https://doi.org/10.1107/S2059798319012877, PMID: 31692470

* Brewster AS, Young ID, Lyubimov A, Bhowmick A, Sauter NK. 2019b. Processing serial crystallographic data from xfels or synchrotrons using the cctbx.xfel GUI. *Computational Crystallography Newsletter* **10**:22–39.

* Carugo O. 2014. Buried chloride stereochemistry in the protein data bank. *BMC Structural Biology* **14**:1–7. DOI: https://doi.org/10.1186/s12900-014-0019-8, PMID: 25928393
Cotruvo JA, Stubbe J. 2008. Nrdl, a flavodoxin involved in maintenance of the diferric-tyrosyl radical cofactor in *Escherichia coli* class ib ribonucleotide reductase. *PNAS* **105**:14383-14388. DOI: https://doi.org/10.1073/pnas.0807348105, PMID: 18799738

Cotruvo JA, Stubbe J. 2010. An active dimanganese(III)-tyrosyl radical cofactor in *Escherichia coli* class ib ribonucleotide reductase. *Biochemistry* **49**:1297-1309. DOI: https://doi.org/10.1021/bi902106n, PMID: 20070127

Cotruvo JA, Stich TA, Britt RD, Stubbe J. 2013. Mechanism of assembly of the dimanganese-tyrosyl radical cofactor of class ib ribonucleotide reductase: enzymatic generation of superoxide is required for tyrosine oxidation via a mn(III)mn(IV) intermediate. *Journal of the American Chemical Society* **135**:4027-4039. DOI: https://doi.org/10.1021/ja312457t, PMID: 23402532

Cox N, Ogata H, Stolle P, Reijerse E, Auling G, Lubitz W. 2010. A tyrosyl-dimanganese coupled spin system is the native metalloradical cofactor of the R2F subunit of the ribonucleotide reductase of corynebacterium ammoniagenes. *Journal of the American Chemical Society* **132**:11197-11213. DOI: https://doi.org/10.1021/ja1036995, PMID: 20698687

Doerr A. 2011. Diffraction before destruction. *Nature Methods* **8**:283. DOI: https://doi.org/10.1038/nmeth0411-283, PMID: 21574275

Emsley P, Lohkamp B, Scott WG, Cowtan K. 2010. Features and development of coot. *Acta Crystallographica. Section D, Biological Crystallography* **66**:486-501. DOI: https://doi.org/10.1107/S0907444910007493, PMID: 20383002

Fuller FD, Gul S, Chatterjee R, Burgie ES, Young ID, Lebrette H, Srinivas V, Brewster AS, Michels-Clark T, Clinger JA, Andi B, Ibrahim M, Pastor E, de Lichtenberg C, Hussein R, Pollock CJ, Zhang M, Stan CA, Kroll T, Fransson T, et al. 2017. Drop-on-demand sample delivery for studying biocatalysts in action at X-ray free-electron lasers. *Nature Methods* **14**:443-449. DOI: https://doi.org/10.1038/nmeth.4195, PMID: 28250468

Gasteiger E, Hoogland C, Gattiker A, Duvaud S, Wilkins MR, Appel RD, Bairoch A. 2005. Protein identification and analysis tools on the expasy server. Walker JM (Ed). *In The Proteomics Protocols Handbook*. Humana Press. p. 571-607. DOI: https://doi.org/10.1385/1-59259-890-0-571

Ghisla S, Massey V. 1989. Mechanisms of flavoprotein-catalyzed reactions. *European Journal of Biochemistry* **181**:1-17. DOI: https://doi.org/10.1111/j.1432-1033.1989.tb14688.x, PMID: 2653819

Grave K, Lambert W, Berggren G, Griese JJ, Bennett MD, Logan DT, Högbom M. 2019. Redox-induced structural changes in the di-iron and di-manganese forms of bacillus anthracis ribonucleotide reductase subunit nrdf suggest a mechanism for gating of radical access. *Journal of Biological Inorganic Chemistry* **24**:849-861. DOI: https://doi.org/10.1007/s00775-019-01703-z, PMID: 31410573

Grimme S, Antony J, Ehrlich S, Krieg H. 2010. A consistent and accurate ab initio parametrization of density functional dispersion correction (DFT-D) for the 94 elements H-pu. *J Chem Phys* **132**:154104. DOI: https://doi.org/10.1063/1.3382344, PMID: 20423165

Hammerstad M, Hersleth HP, Tomter AB, Røhr AK, Andersson KK. 2014. Crystal structure of bacillus cereus class ib ribonucleotide reductase di-iron nrdf in complex with nrdi. *ACS Chemical Biology* **9**:526-537. DOI: https://doi.org/10.1021/cb400757h, PMID: 24295378

Hattne J, Echols N, Tran R, Kern J, Gildea RJ, Brewster AS, Alonso-Mori R, Glöckner C, Hellmich J, Laksmono H, Sierra RG, Lassalle-Kaiser B, Lampe A, Han G, Gul S, DiFiore D, Milathianaki D, Fry AR, Miahnahri A, White WE, et al. 2014. Accurate macromolecular structures using minimal measurements from X-ray free-electron lasers. *Nature Methods* **11**:545-548. DOI: https://doi.org/10.1038/nmeth.2887, PMID: 24633409

Högbom M, Sjöberg B, Berggren G. 2020. Radical enzymes. *ELS* **1**:375-393. DOI: https://doi.org/10.1002/9780470015902.a0029205

Huque Y, Fieschi F, Torrents E, Gibert I, Eliasson R, Reichard P, Sahlin M, Sjoberg BM. 2000. The active form of the R2F protein of class ib ribonucleotide reductase from corynebacterium ammoniagenes is a diferric protein. *The Journal of Biological Chemistry* **275**:25365-25371. DOI: https://doi.org/10.1074/jbc.M002751200, PMID: 10801858

Johansson R, Torrents E, Lundin D, Sprenger J, Sahlin M, Sjöberg BM, Logan DT. 2010. High-resolution crystal structures of the flavoprotein nrdi in oxidized and reduced states--an unusual flavodoxin: structural biology. *The FEBS Journal* **277**:4265-4277. DOI: https://doi.org/10.1111/j.1742-4658.2010.07815.x, PMID: 20831589

Klamt A, Schüürmann G. 1993. COSMO: A new approach to dielectric screening in solvents with explicit expressions for the screening energy and its gradient. *Journal of the Chemical Society, Perkin Transactions* **2**:799-805. DOI: https://doi.org/10.1039/P29930000799

Krissinel E, Henrick K. 2004. Secondary-structure matching (SSM), a new tool for fast protein structure alignment in three dimensions. *Acta Crystallographica. Section D, Biological Crystallography* **60**:2256-2268. DOI: https://doi.org/10.1107/S0907444904026460, PMID: 15572779

Lee C, Yang W, Parr RG. 1988. Development of the colle-salvetti correlation-energy formula into a functional of the electron density. *Physical Review. B, Condensed Matter* **37**:785-789. DOI: https://doi.org/10.1103/physrevb.37.785, PMID: 9944570

Liebschner D, Afonine PV, Baker ML, Bunkóczi G, Chen VB, Croll TI, Hintze B, Hung LW, Jain S, McCoy AJ, Moriarty NW, Oeffner RD, Poon BK, Prisant MG, Read RJ, Richardson JS, Richardson DC, Sammito MD, Sobolev OV, Stockwell DH, et al. 2019. Macromolecular structure determination using X-rays, neutrons and electrons: recent developments in phenix. *Acta Crystallographica. Section D, Structural Biology* **75**:861-877. DOI: https://doi.org/10.1107/S2059798319011471, PMID: 31588918

Makhlynets O, Boal AK, Rhodes DV, Kitten T, Rosenzweig AC, Stubbe J. 2014. Streptococcus sanguinis class ib ribonucleotide reductase: high activity with both iron and manganese cofactors and structural insights. *The*
Journal of Biological Chemistry **289**:6259-6272. DOI: https://doi.org/10.1074/jbc.M113.533554, PMID: 24381172

McCoy AJ, Grosse-Kunstleve RW, Adams PD, Winn MD, Storoni LC, Read RJ. 2007. Phaser crystallographic software. *Journal of Applied Crystallography* **40**:658–674. DOI: https://doi.org/10.1107/S0021889807021206, PMID: 19461840

Moriarty NW, Grosse-Kunstleve RW, Adams PD. 2009. Electronic ligand builder and optimization workbench (elbow): a tool for ligand coordinate and restraint generation. *Acta Crystallographica. Section D, Biological Crystallography* **65**:1074–1080. DOI: https://doi.org/10.1107/S0907444909029436, PMID: 19770504

Nass K. 2019. Radiation damage in protein crystallography at X-ray free-electron lasers. *Acta Crystallographica*. *Section D, Structural Biology* **75**:211–218. DOI: https://doi.org/10.1107/S2059798319000317, PMID: 30821709

Neutze R, Wouts R, van der Spoel D, Weckert E, Hajdu J. 2000. Potential for biomolecular imaging with femtosecond X-ray pulses. *Nature* **406**:752–757. DOI: https://doi.org/10.1038/35021099, PMID: 10963603

Nordlund P, Reichard P. 2006. Ribonucleotide reductases. *Annual Review of Biochemistry* **75**:681–706. DOI: https://doi.org/10.1146/annurev.biochem.75.103004.142443, PMID: 16756507

Pettersen EF, Goddard TD, Huang CC, Couch GS, Greenblatt DM, Meng EC, Ferrin TE. 2004. UCSF chimera--a visualization system for exploratory research and analysis. *Journal of Computational Chemistry* **25**:1605–1612. DOI: https://doi.org/10.1002/jcc.20084, PMID: 15264254

Roca I, Torrents E, Sahlin M, Gibert I, Sjöberg BMM. 2008. Nrdl essentiality for class ib ribonucleotide reduction in *Streptococcus pyogenes*. *Journal of Bacteriology* **190**:4849–4858. DOI: https://doi.org/10.1128/JB.00185-08, PMID: 18502861

Røhr AK, Hersleth HP, Andersson KK. 2010. Tracking flavin conformations in protein crystal structures with raman spectroscopy and QM/MM calculations. *Angewandte Chemie* **49**:2324–2327. DOI: https://doi.org/10.1002/anie.200907143, PMID: 20187055

Rould MA, Carter CW. 2003. Isomorphous difference methods. *Methods in Enzymology* **374**:145–163. DOI: https://doi.org/10.1016/S0076-6879(03)74007-5, PMID: 14696372

Sauter NK. 2015. XFEL diffraction: developing processing methods to optimize data quality. *Journal of Synchrotron Radiation* **22**:239–248. DOI: https://doi.org/10.1107/S1600577514028203, PMID: 25723925

Schäfer A, Horn H, Ahlrichs R. 1992. Fully optimized contracted gaussian basis sets for atoms li to kr. *J Chem Phys* **97**:2571–2577. DOI: https://doi.org/10.1063/1.463096

Senda T, Senda M, Kimura S, Ishida T. 2009. Redox control of protein conformation in flavoproteins. *Antioxidants & Redox Signaling* **11**:1741–1766. DOI: https://doi.org/10.1089/ARS.2008.2348, PMID: 19243237

Sierra RG, Batyuk A, Sun Z, Aquila A, Hunter MS, Lane TJ, Liang M, Yoon CH, Alonso-Mori R, Armenta R, Castagna JC, Hollenbeck M, Osier TO, Hayes M, Aldrich J, Curtis R, Koglin JE, Rendahl T, Rodriguez E, Carbajo S, et al. 2019. The macromolecular femtosecond crystallography instrument at the linac coherent light source. *Journal of Synchrotron Radiation* **26**:346–357. DOI: https://doi.org/10.1107/S1600577519001577, PMID: 30855242

Sorigué D, Hadjidemetriou K, Blangy S, Gotthard G, Bonvalet A, Coquelle N, Samire P, Aleksandrov A, Antonucci L, Benachir A, Boutet S, Byrdin M, Cammarata M, Carbajo S, Cuiné S, Doak RB, Foucar L, Gorel A, Grünbein M, Hartmann E, et al. 2021. Mechanism and dynamics of fatty acid photodecarboxylase. *Science* **372**:6538. DOI: https://doi.org/10.1126/science.abd5687, PMID: 33833098

Spence JCH. 2017. XFELs for structure and dynamics in biology. *IUCR J* **4**:322–339. DOI: https://doi.org/10.1107/S2052252517005760, PMID: 28875020

Srinivas V, Lebrette H, Lundin D, Kutin Y, Sahlin M, Lerche M, Eirich J, Branca RMM, Cox N, Sjöberg BM, Högbom M. 2018. Metal-free ribonucleotide reduction powered by a DOPA radical in mycoplasma pathogens. *Nature* **563**:416–420. DOI: https://doi.org/10.1038/s41586-018-0653-6, PMID: 30429545

Sugahara M, Mizohata E, Nango E, Suzuki M, Tanaka T, Masuda T, Tanaka R, Shimamura T, Tanaka Y, Suno C, Ihara K, Pan D, Kakinouchi K, Sugiyama S, Murata M, Inoue T, Tono K, Song C, Park J, Kameshima T, et al. 2015. Grease matrix as a versatile carrier of proteins for serial crystallography. *Nature Methods* **12**:61–63. DOI: https://doi.org/10.1038/nmeth.3172, PMID: 25384243

Sugahara M, Nakane T, Masuda T, Suzuki M, Inoue S, Song C, Tanaka R, Nakatsu T, Mizohata E, Yumoto F, Tono K, Joti Y, Kameshima T, Hatsui T, Yabashi M, Nureki O, Numata K, Nango E, Iwata S. 2017. Hydroxyethyl cellulose matrix applied to serial crystallography. *Scientific Reports* **7**:1–9. DOI: https://doi.org/10.1038/s41598-017-00761-0, PMID: 28386083

Tomter AB, Bell CB, Røhr AK, Andersson KK, Solomon EI. 2008. Circular dichroism and magnetic circular dichroism studies of the biferrous site of the class ib ribonucleotide reductase from bacillus cereus: comparison to the class ia enzymes. *Biochemistry* **47**:11300–11309. DOI: https://doi.org/10.1021/bi801212f, PMID: 18831534

Tono K, Togashi T, Inubushi Y, Sato T, Katayama T, Ogawa K, Ohashi H, Kimura H, Takahashi S, Takeshita K, Tomizawa H, Goto S, Ishikawa T, Yabashi M. 2013. Beamline, experimental stations and photon beam diagnostics for the hard x-ray free electron laser of SACLA. *New Journal of Physics* **15**:083035. DOI: https://doi.org/10.1088/1367-2630/15/8/083035

Wallace AC, Laskowski RA, Thornton JM. 1995. LIGPLOT: a program to generate schematic diagrams of protein-ligand interactions. *Protein Engineering, Design and Selection* **8**:127–134. DOI: https://doi.org/10.1093/protein/8.2.127, PMID: 7630882

Walsh JD, Miller AF. 2003. Flavin reduction potential tuning by substitution and bending. *Journal of Molecular Structure* **623**:185–195. DOI: https://doi.org/10.1016/S0166-1280(02)00719-4, PMID: 12603753
Williams CJ, Headd JJ, Moriarty NW, Prisant MG, Videau LL, Deis LN, Verma V, Keedy DA, Hintze BJ, Chen VB, Jain S, Lewis SM, Arendall WB, Snoeyink J, Adams PD, Lovell SC, Richardson JS, Richardson DC. 2018. MolProbity: more and better reference data for improved all-atom structure validation. *Protein Science* **27**:293-315. DOI: https://doi.org/10.1002/pro.3330, PMID: 29067766

Winter G, Waterman DG, Parkhurst JM, Brewster AS, Gildea RJ, Gerstel M, Fuentes-Montero L, Vollmar M, Michels-Clark T, Young ID, Sauter NK, Evans G. 2018. DIALS: implementation and evaluation of a new integration package. *Acta Crystallographica Section D Structural Biology* **74**:85-97. DOI: https://doi.org/10.1107/S2059798317017235, PMID: 29533234

Zheng YJ, Ornstein RL. 1996. A theoretical study of the structures of flavin in different oxidation and protonation states. *Journal of the American Chemical Society* **118**:9402-9408. DOI: https://doi.org/10.1021/ja9608151