# The FAD synthetase from the human pathogen Streptococcus pneumoniae: a bifunctional enzyme exhibiting activity-dependent redox requirements

Maria Sebastian, Erandi Lira-Navarrete, Ana Serrano, Carlos Marcuello

Angles, Adrián Velázquez-Campoy, Anabel Lostao, Ramon Hurtado-Guerrero,

Milagros Medina, Marta Martinez-Julvez

► To cite this version:

Maria Sebastian, Erandi Lira-Navarrete, Ana Serrano, Carlos Marcuello Angles, Adrián Velázquez-Campoy, et al.. The FAD synthetase from the human pathogen Streptococcus pneumoniae: a bi-functional enzyme exhibiting activity-dependent redox requirements. Scientific Reports, 2017, 7, 10.1038/s41598-017-07716-5. [hal-02619626](https://hal.inrae.fr/hal-02619626)

HAL Id: hal-02619626

https://hal.inrae.fr/hal-02619626

Submitted on 25 May 2020

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>

Distributed under a Creative Commons Attribution 4.0 International License
SCIENTIFIC REPORTS

OPEN

# The FAD synthetase from the human pathogen *Streptococcus pneumoniae*: a bifunctional enzyme exhibiting activity-dependent redox requirements

María Sebastián<fnref n="1" />,², Erandi Lira-Navarrete²,⁶, Ana Serrano<fnref n="1" />,²,⁷, Carlos Marcuello³,⁸, Adrián Velázquez-Campoy<fnref n="1" />,²,⁴,⁵, Anabel Lostao ID³,⁴, Ramón Hurtado-Guerrero²,⁴, Milagros Medina ID<fnref n="1" />,² & Marta Martínez-Júlvez ID<fnref n="1" />,²

Received: 2 May 2017

Accepted: 28 June 2017

Published online: 08 August 2017

Prokaryotic bifunctional FAD synthetases (FADSs) catalyze the biosynthesis of FMN and FAD, whereas in eukaryotes two enzymes are required for the same purpose. FMN and FAD are key cofactors to maintain the flavoproteome homeostasis in all type of organisms. Here we shed light to the properties of the hitherto unstudied bacterial FADS from the human pathogen *Streptococcus pneumoniae* (*SpnFADS*). As other members of the family, *SpnFADS* catalyzes the three typical activities of prokaryotic FADSs: riboflavin kinase (RFK), ATP:FMN:adenylyltransferase (FMNAT), and FAD pyrophosphorylase (FADpp). However, several *SpnFADS* biophysical properties differ from those of other family members. In particular; i) the RFK activity is not inhibited by the riboflavin (RF) substrate, ii) the FMNAT and FADSpp activities require flavin substrates in the reduced state, iii) binding of adenine nucleotide ligands is required for the binding of flavinic substrates/products and iv) the monomer is the preferred state. Collectively, our results add interesting mechanistic differences among the few prokaryotic bifunctional FADSs already characterized, which might reflect the adaptation of the enzyme to relatively different environments. In a health point of view, differences among FADS family members provide us with a framework to design selective compounds targeting these enzymes for the treatment of diverse infectious diseases.

Prokaryotic FAD synthetases (FADSs) are bifunctional and bimodular enzymes that first catalyze the synthesis of flavin mononucleotide (FMN) from riboflavin (RF, vitamin B₂) through an ATP:riboflavin kinase activity (RFK, EC 2.7.1.26) and subsequently the adenylylation of FMN to produce flavin adenine dinucleotide (FAD) using an ATP:FMN:adenylyltransferase activity (FMNAT, EC 2.7.7.2)<fnref n="1" />⁻³. For some of them, the FMNAT activity has also been shown to be reversible, exhibiting therefore a FAD pyrophosphorylase activity (FADpp)⁴,⁵. In mammals and yeasts two independent enzymes catalyze these reactions⁶⁻⁹. Eukaryotic RFKs exhibit structural homology with the C-terminal module of prokaryotic bifunctional FADSs¹⁰⁻¹². On the contrary eukaryotic FMNATs do not exhibit similarity to the prokaryotic N-terminal modules, which belong to the nucleotidyltransferase

<footnote n="1">¹Departamento de Bioquímica y Biología Molecular y Celular, Facultad de Ciencias, Universidad de Zaragoza, Zaragoza, Spain. ²Instituto de Biocomputación y Física de Sistemas Complejos (BIFI) and GBsC-CSIC and BIFI-CSIC Joint Units, Universidad de Zaragoza, Zaragoza, Spain. ³Laboratorio de Microscopías Avanzadas (LMA), Instituto de Nanociencia de Aragón (INA), and Fundación INA, Universidad de Zaragoza, Zaragoza, Spain. ⁴Fundación ARAID, Diputación General de Aragón, Aragón, Spain. ⁵Aragon Institute for Health Research (IIS Aragon), Zaragoza, 50009, Spain. ⁶Present address: Copenhagen Center for Glycomics, Department of Cellular and Molecular Medicine, School of Dentistry, University of Copenhagen, Copenhagen, DK-2200, Denmark. ⁷Present address: Centro de Investigaciones Biológicas, CSIC, Ramiro de Maeztu 9, E-28040, Madrid, Spain. ⁸Present address: Univ Reims, Lab Rech Nanosci, EA4682, F-51100 Reims and INRA, FARE Lab, F-51100, Reims, France. María Sebastián and Erandi Lira contributed equally to this work. Correspondence and requests for materials should be addressed to M.M. (email: mmedina@unizar.es) or M.M.-J. (email: mmartine@unizar.es)</footnote>
superfamily<sup>13-15</sup>. Since FMN and FAD are essential and versatile cofactors in a plethora of vital redox processes<sup>16-20</sup>, and also participate in signal transduction during apoptosis<sup>21</sup>, embryonic development<sup>22</sup>, chromatin remodeling<sup>23</sup>, nucleotide synthesis<sup>24</sup>, tRNA methylation<sup>25</sup>, protein folding<sup>26</sup>, and light dependent processes in plants as well as in sensory transduction in microorganisms<sup>27-31</sup>, their biosynthesis is worth to be understood.

Structures of prokaryotic FADSs are available only for *Thermotoga maritima* (TmFADS, PDB ID: 1MRZ,¹¹,³²), *Corynebacterium ammoniagenes* (CaFADS, PDB ID: 2X0K¹⁰) and the human pathogen *Streptococcus pneumoniae* (*SpnFADS*, PDB ID: 3OP1). Except for the exhaustive biochemical and biophysical characterization of CaFADS¹,⁵,³³, the rest of family members have been poorly studied⁴, ³⁴, ³⁵. *SpnFADS* shows the typical sequence and structural features of prokaryotic bifunctional FADSs, and as such it would be predicted to contain RFK, FMNAT and FADpp activities (Supplementary Fig. S1; r.m.s.d. of 1.63 Å encompassing 247 residues and 30% sequence identity with CaFADS). The N-terminal module (residues 1–186) is an α/β dinucleotide binding domain that by homology to other family members is expected to catalyze the FMNAT and FADpp activities¹⁰, ¹¹. Likewise, the C-terminal module (residues 187–305), predicted to catalyze the RFK activity¹⁰, ³⁶, consists of a β barrel formed by six antiparallel β-strands and a long terminal α-helix. The different structures together with the significant differences at the catalytic centers between the N-terminal modules of bifunctional FADSs and the mammalian FMNAT enzymes¹⁰, ¹⁴ make prokaryotic FADSs attractive exploitable targets for the treatment of human infectious diseases³⁷. In this context, *SpnFADS* might be considered as a drug target for *S. pneumoniae*, a pathogen almost exclusive for humans, that is the leading cause of invasive bacterial pneumoniae disease in children, in the elderly and in immunodepressed patients³⁸⁻⁴⁰.

We present here the biochemical and biophysical characterization of *SpnFADS*, the only FADS from a pathogenic prokaryote whose structure is available. Our data indicate that *SpnFADS* catalyzes the RFK, FMNAT and FADpp activities. However, it differs from other family members in two major catalytic aspects: its RFK activity does not exhibit inhibition by the RF substrate and the FMNAT activity requires the flavin substrate in its reduced state for substrate binding and catalysis. In addition, the prior binding of adenine nucleotides is required for the recognition of flavins, and the enzyme exits mainly in the monomeric state during catalysis. These results exemplify how the same type of enzymes from different organisms has evolved different strategies to perform the same reaction, presumably due to their unique environment. Our data also point to these bacterial enzymes as promising pharmacological targets because not only differ between themselves but also from the human enzyme.

## Results

**SpnFADS is purified free of flavin ligands.** Typically, the yield after purification was 3.5 mg of *SpnFADS* per liter of *E.coli* culture. Enriched *SpnFADS* fractions along purification did not show notorious absorbance at 450 nm, indicating the purified protein does not contain any flavin bound under these conditions. Thus, the final absorption spectrum of purified *SpnFADS* exhibits a single peak at 279 nm (Supplementary Fig. S2) and in contrast to *CaFADS*¹³, *SpnFADS* does not co-purify with FAD, suggesting a much lower affinity for the oxidized flavinic substrates/products. The far-ultraviolet (UV) circular dichroism (CD) spectrum of *SpnFADS* exhibits a negative band at 208 nm typical of secondary α-helix structure and a slightly less intense negative band around 222 nm (Supplementary Fig. S2). Analysis of the far-UV CD spectrum using the method of Raussens⁴¹ predicted 28.3, 19.9, 12.5 and 34.5% of content in α helix, β strand, turn and random, respectively, in agreement with its 3D-structure containing 29% α helix and 28% β strand (Supplementary Fig. S1). Near-UV CD spectrum of *SpnFADS* shows a broad 265–290 nm negative band with minima at 280.5 and 287 nm related with the protein tertiary structure organization (Supplementary Fig. S2). Altogether these data indicate that the purified protein is correctly folded.

**SpnFADS catalyzes the RFK, FMNAT and FADpp activities.** A qualitative assay was initially performed to evaluate whether *SpnFADS* was able to catalyze the RFK, FMNAT and FADpp activities by resolving the flavin reaction products by thin layer chromatography (TLC). As shown in Fig. 1A, *SpnFADS* catalyzes the conversion of RF into FMN, both in the absence and in the presence of a reducing agent. However, reducing conditions (2.5 mM of sodium dithionite) are required for the transformation of FMN into FAD (Fig. 1A), as well as of FAD into FMN (Fig. 1B).

**The reduced state of flavins is critical for the FMNAT activity but not for the RFK one.** Prior to determining the kinetic parameters for the RFK activity of *SpnFADS*, we quantitatively evaluated the influence of the reducing environment (sodium dithionite) as well as of Mg²⁺ on this reaction. Under saturating concentrations of substrates (RF and ATP) rates for RF transformation were independent on the reducing agent concentration (Supplementary Fig. S3), indicating that the *SpnFADS* RFK activity does not depend on its redox environment. On the contrary, Mg²⁺ enhances the RFK activity, being the rate of the process practically constant at cation concentrations over 0.4 mM (Supplementary Fig. S3). Thus, herein a concentration of 0.8 mM MgCl₂ was used. Steady-state rates for the RFK activity showed saturation profiles for both substrates (Fig. 1C), contrary to *CaFADS* that has strong inhibition by excess of RF⁵, ³³. Thus, profiles for the RFK activity of *SpnFADS* were fit to the Michaelis-Menten model, allowing determination of *k*<sub>cat</sub>, *K*<sub>m</sub><sup>RF</sup> and *K*<sub>m</sub><sup>ATP</sup> values of 18.7 min⁻¹, 0.9 µM and 54.0 µM, respectively. These values were only slightly influenced by the presence of the reducing agent (Table 1, Fig. 1C).

The FMNAT activity of *SpnFADS* was determined following the transformation of FMN into FAD. We again evaluated first the influence of Mg²⁺ concentration and of reducing conditions in the FMNAT activity. The presence of sodium dithionite was a requirement for the transformation of FMN into FAD (Supplementary Fig. S3), being the rate of the process constant over 2 mM sodium dithionite. On the contrary, the FMNAT activity was very poor in the presence of milder reducing agents as DTT (not shown), suggesting that strong reducing conditions are required to achieve catalysis. Maximal activity was observed in the range of 0.8–1.2 mM of the divalent
<figure><img src="image_2.png" /><figcaption>**Figure 1.** The RFK, FMNAT and FADpp activities of *SpnFADS*. (A) Thin layer chromatography of the products of the RFK and FMNAT activities. Reaction mixtures contained around 500 nM *SpnFADS*, 50 µM RF, 0.2 mM ATP, and several concentrations of sodium dithionite. First lane to the left corresponds to a standard solution containing RF, FMN and FAD at 50 µM each. (B) Thin layer chromatography of the products of FADpp activity. Reaction mixtures contained around 500 nM *SpnFADS*, 50 µM FAD, 0.2 mM PPI, and several concentrations of sodium dithionite. First lane corresponds to the standard solutions and second line to a free FAD sample with 10 mM MgCl₂. (C) Rates for RFK activity as a function of the RF concentration at saturating ATP (412 µM) without (●) and with 2.6 mM (■) sodium dithionite. (D) Rates for the FMNAT activity as a function of the FMN concentration at saturating ATP (252 µM) and 3 mM sodium dithionite. Unless otherwise stated all measurements were carried out at 25 °C in 20 mM PIPES, pH 7.0, with 0.8 mM MgCl₂.</figcaption></figure>

cation (Supplementary Fig. S3), being the observed rates in its absence up to 5-fold lower. Thus, 0.8 mM MgCl₂ and 3 mM sodium dithionite were selected to further characterize the FMNAT kinetic parameters of *SpnFADS*. Observed rates exhibited Michaelis-Menten saturation profiles (Fig. 1D) with $k_{cat}$, $K_m^{FMN}$ and $K_m^{ATP}$ values of 34.9 min⁻¹, 9.8 µM and 31.6 µM, respectively (Table 1). The need of reducing conditions for the FMNAT activity suggests that either the isoalloxazine ring of FMN needs to be reduced for the reaction to occur, or the reaction depends on a conformational change of the flavin binding site of the FMNAT module of *SpnFADS* induced by the redox environment. To clarify this matter an additional experiment was performed. The FMNAT activity of *SpnFADS* was measured under anaerobic conditions, but in the absence of dithionite, using as substrate the photo-reduced FMN. HPLC analysis of the reaction products revealed that *SpnFADS* is able to transform FMN into FAD in the samples containing photo-reduced FMN in contrast with the negative controls where FMN is oxidized. Therefore, we conclude that FMN needs to be in its reduced form to get adenylated by *SpnFADS*.

**SpnFADS does not bind oxidized flavins by itself.** Isothermal titration calorimetry (ITC) thermograms for the titration of *SpnFADS* with oxidized RF, FMN or FAD, either in presence or absence of Mg²⁺, indicated that under the assayed aerobic conditions none of these flavins binds to the protein with detectable enthalpic contribution (Supplementary Fig. S4). Therefore, we also titrated *SpnFADS* with FMN and FAD in a flavin reducing buffer and in the presence of 0.8 mM Mg²⁺. As shown in Fig. 2A, under such conditions binding of the reduced
<table><caption>**Table 1.** Steady-state kinetic parameters for the RFK and FMNAT activities of SpnFADS. Parameters obtained at 25 °C in 20 mM PIPES, pH 7.0, 0.8 mM MgCl₂ and, when indicated 3 mM sodium dithionite as reductant. Estimated errors in kcat and Km are within ± 10%. $^a$kcat measured at 6 µM of RF. $^b$Data from$^{33}$. Data were obtained at 25 °C in 20 mM PIPES, pH 7.0 and 0.8 mM and 10 mM MgCl₂, respectively, for the RFK and FMNAT activities. $^c$Inhibition by substrate prevented the determination of true parameters and these correspond to apparent constants, $k_{cat,app}$, $K_m^{RF}$, $K_{m,app}^{ATP}$, $k_{cat,app}/K_{m,app}^{RF}$ and $k_{cat,app}/K_{m,app}^{ATP}$. $^d$Data from$^{35}$. Data obtained in 50 mM potassium phosphate, pH 7.5. $^e$kcat measured at 27.7 µM of FMN. $^f$Data from$^{42}$. Data obtained in 50 mM potassium phosphate, pH 7.5. $^g$Data from$^{34}$. Data obtained in 100 mM potassium phosphate, pH 7.5. $^h$FADS from *Escherichia coli*. Data from$^{43}$. Data obtained in 50 mM potassium phosphate, pH 7.5.</caption><thead><tr><th>RFK activity</th><th>Reductant</th><th>k<sub>cat</sub> (min<sup>-1</sup>)</th><th>K<sub>m</sub><sup>RF</sup> (μM)</th><th>k<sub>cat</sub>K<sub>m</sub><sup>RF</sup> (min<sup>-1</sup> μM<sup>-1</sup>)</th><th>k<sub>cat</sub> (min<sup>-1</sup>)</th><th>K<sub>m</sub><sup>ATP</sup> (μM)</th><th>k<sub>cat</sub>K<sub>m</sub><sup>ATP</sup> (min<sup>-1</sup> μM<sup>-1</sup>)</th></tr></thead><tbody><tr><td rowspan="2">SpnFADS</td><td>–</td><td>18.7</td><td>0.9</td><td>21.0</td><td>19.5<sup>a</sup></td><td>54.0</td><td>0.4</td></tr><tr><td>+</td><td>24.1</td><td>0.6</td><td>41.5</td><td>13.4<sup>a</sup></td><td>34.1</td><td>0.4</td></tr><tr><td>CaFADS<sup>b,c</sup></td><td>–</td><td>&lt; 302</td><td>&lt; 13</td><td>23.2</td><td>68</td><td>13.7</td><td>4.9</td></tr><tr><td>LmFADSI<sup>d</sup></td><td>+</td><td>0.40</td><td>6.9</td><td>0.057</td><td></td><td></td><td></td></tr><tr><td>HsRFK<sup>f</sup></td><td>+</td><td>0.50</td><td>36</td><td>0.014</td><td></td><td></td><td></td></tr><tr><td>BsFADS<sup>g</sup></td><td>+</td><td>0.70</td><td>55</td><td>0.013</td><td></td><td></td><td></td></tr><tr><td>SdFADS<sup>g</sup></td><td>+</td><td>0.30</td><td>40</td><td>0.007</td><td></td><td></td><td></td></tr><tr><td>EcFADS<sup>h</sup></td><td>+</td><td>0.39</td><td>2</td><td>0.20</td><td></td><td></td><td></td></tr><tr><th>FMNAT activity</th><th></th><th>k<sub>cat</sub> (min<sup>-1</sup>)</th><th>K<sub>m</sub><sup>FMN</sup> (μM)</th><th>k<sub>cat</sub>K<sub>m</sub><sup>FMN</sup> (min<sup>-1</sup> μM<sup>-1</sup>)</th><th>k<sub>cat</sub> (min<sup>-1</sup>)</th><th>K<sub>m</sub><sup>ATP</sup> (μM)</th><th>k<sub>cat</sub>K<sub>m</sub><sup>ATP</sup> (min<sup>-1</sup> μM<sup>-1</sup>)</th></tr><tr><td>SpnFADS</td><td>+</td><td>34.9</td><td>9.8</td><td>3.6</td><td>23.2<sup>e</sup></td><td>31.6</td><td>0.7</td></tr><tr><td>CaFADS<sup>b</sup></td><td>–</td><td>17.0</td><td>1.2</td><td>14.3</td><td>17.0</td><td>35.7</td><td>0.5</td></tr><tr><td>LmFADSI<sup>d</sup></td><td>+</td><td>4.39</td><td>29.2</td><td>0.15</td><td></td><td></td><td></td></tr><tr><td>LmFADSI<sup>d</sup></td><td>+</td><td>0.51</td><td>12.9</td><td>0.04</td><td></td><td></td><td></td></tr><tr><td>EcFADS<sup>h</sup></td><td>+</td><td>0.06</td><td>4</td><td>0.015</td><td></td><td></td><td></td></tr></tbody></table>

forms of FMN and FAD to the enzyme was detected. Both reduced cofactors bind to a single SpnFADS binding site (Fig. 2A and Table 2). $K_d^{FMN}$ is in the low micromolar range, and binding of reduced FAD is 15-fold weaker. Enthalpic contribution to the binding are very similar for both flavins, but the entropic contributions favor FMN binding (Fig. 2C). These results thus indicate that the FMNAT/FADpp sites of SpnFADS only bind the reduced states of FMN or FAD.

Thermograms for the titration of SpnFADS with ATP indicated a single binding site for this nucleotide (presumably at the FMNAT site as in CaFADS), both in the absence and presence of 0.8 mM of Mg$^{2+}$ (Fig. 2B and Table 2), although binding affinities and interaction parameters differ between these two conditions. Thus, ATP and ATP:Mg$^{2+}$ allocate differently into the protein binding site, suggesting that the ion contributes to the nucleotide-protein-bound conformation. Titrations at higher Mg$^{2+}$ concentrations showed two independent ATP binding sites, putatively related with each catalytic site, with similar affinities between them, and also similar to that of the site found at lower Mg$^{2+}$ concentrations. However, enthalpic and entropic contributions to the ATP binding result highly influenced by Mg$^{2+}$; in its absence, and at low concentrations, the nucleotide binding is enthalpically driven with an opposing entropic contribution, while at high cation concentrations relatively mild entropic and enthalpic contributions favor the binding (Fig. 2C). Thermodynamic parameters for the SpnFADS titration with ADP also differ depending on Mg$^{2+}$ concentration. In the absence of the cation, a single high affinity binding site is observed, being the binding driven by a favorable enthalpic contribution. At 10 mM of Mg$^{2+}$, a second binding site appears with lower affinity for the nucleotide, being the averaged $K_d$ ($K_{d,av}$) at least 7-times higher. ADP binding in the presence of 10 mM Mg$^{2+}$ is also enthalpically driven, but in this case the enthalpic contribution is much higher, being the entropic term mightily unfavorable. Thus high Mg$^{2+}$ concentrations induce the generation of a second binding site for the adenine nucleotide ligands, both ATP and ADP, and influence the way the ligands allocate into the active site, as also supported by the different thermodynamic contributions to free energy of the binding.

Thermograms for the titration of mixtures containing SpnFADS plus ADP with oxidized RF or FMN showed that the presence of the adenine nucleotide promotes binding of both of these flavins (Fig. 2 and Table 2). Under these conditions RF in presence of 10 mM MgCl₂, apparently binds at two independent binding sites with moderate affinity, $K_{d,av}$ 51.3 µM, being the averaged binding enthalpically driven and the entropic contribution highly unfavorable. However, in the absence of Mg$^{2+}$ only one binding site with strong affinity is detected. By contrast FMN binds at a single site with stronger affinity than those found for RF, being the binding driven by the enthalpic contribution and the binding parameters just modulated by the presence of Mg$^{2+}$. Finally, binding of FAD to SpnFADS in the presence of PPI could not be detected by enthalpic change, similarly to that occurred when titrating the PPI-free enzyme with oxidized FAD.

**SpnFADS exits mainly as a monomer but other oligomeric states are populated during catalysis.** SpnFADS eluted in a main peak with a small shoulder in the gel filtration column used as a last step in its purification. The presence of the shoulder might be indicative of freshly purified SpnFADS stabilizing quaternary assemblies, as described for CaFADS$^1$. To check for such possibility we passed through the gel filtration
<figure><img src="image_3.png" /><figcaption>B</figcaption></figure>

<figure><img src="image_4.png" /><figcaption>C</figcaption></figure>

<figure><img src="image_5.png" /></figure>

**Figure 2.** Binding of ligands to *SpnFADS*. (A) Calorimetric titrations of *SpnFADS* with reduced FMN and FAD, performed at 20 °C in 20 mM PIPES, pH 7.0, 0.8 mM MgCl₂, 4 mM sodium dithionite, 10 mM glucose, 0.5 U/mL glucose oxidase (degassed). (B) Calorimetric titrations of *SpnFADS* with ATP and of the preformed *SpnFADS*:ADP complex with RF or FMN, at 20 °C in 20 mM PIPES, pH 7.0, and at the indicated concentrations of MgCl₂. In all titration figures upper panels show the thermograms for the interactions and the lower panels show the corresponding binding isotherms with integrated normalized heats. (C) Thermodynamic dissections of the interaction of *SpnFADS* with reduced FAD, reduced FMN and ATP, as well as of the preformed *SpnFADS*:ADP complex with RF and FMN. The binding Gibbs energy (ΔG), enthalpy (ΔH), and entropy (−TΔS) are represented in dark grey, light grey and white bars, respectively. MgCl₂ concentrations in mM used in each experiment are indicated by numbers.

column a freshly purified *SpnFADS* fraction considered as monomeric. This protein eluted as a main broad peak at ~14.5 mL (~96% of the protein) but also produced a shoulder at ~13 mL (Fig. 3A). The peak at ~14.5 mL deconvolutes into three smaller peaks with apparent MWs of 61.2, 63.7 and 70 kDa. These values are slightly higher than the theoretic one for the monomer (34.8 kDa). The elongated structure of *SpnFADS* and the different conformations expected for its external loops could justify the reasonable assignment of those peaks to either different conformers of monomers or dimers. In contrast, the apparent MW of the shoulder at 13 mL is 127.4 kDa,
<table><caption>**Table 2.** Thermodynamic parameters for the interaction of *SpnFADS* with the substrates and products of its activities obtained by ITC. Data obtained at 20 °C in 20 mM PIPES, pH 7.0. Estimated errors in *K*<sub>d</sub> are ±20% and ±0.3 kcal/mol in ΔH and -TΔS. n.d. No enthalpic change was detected along the titration. <sup>a</sup>Measurements performed in the presence of 4 mM sodium dithionite, 10 mM glucose and 0.5 U/ml glucose oxidase. <sup>b</sup>Data corresponding to average values of two similar and independent binding sites (*K*<sub>d,av</sub>, ΔG<sub>av</sub>, ΔH<sub>av</sub>, -TΔS<sub>av</sub>).</caption><thead><tr><th>Receptor</th><th>Ligand</th><th>MgCl<sub>2</sub> (mM)</th><th>K<sub>d</sub> (µM)</th><th>ΔG (kcal/mol)</th><th>ΔH (kcal/mol)</th><th>-TΔS (kcal/mol)</th><th>N (binding stoichiometry)</th></tr></thead><tbody><tr><td rowspan="12">SpnFADS</td><td rowspan="3">ATP</td><td>0</td><td>1.7</td><td>-7.9</td><td>-19.1</td><td>11.2</td><td>~1</td></tr><tr><td>0.8</td><td>9.4</td><td>-6.9</td><td>-27.8</td><td>23.0</td><td>~1</td></tr><tr><td>10</td><td>9.2<sup>b</sup></td><td>-6.9<sup>b</sup></td><td>-3.0<sup>b</sup></td><td>-3.9<sup>b</sup></td><td>~2</td></tr><tr><td rowspan="2">ADP</td><td>0</td><td>4.1</td><td>-7.3</td><td>-12.9</td><td>5.5</td><td>~1</td></tr><tr><td>10</td><td>30.0</td><td>-6.2</td><td>-19.4</td><td>13.3</td><td>~1</td></tr><tr><td rowspan="2">RF, FMN, FAD</td><td>0</td><td>n.d.</td><td>n.d.</td><td>n.d.</td><td>n.d.</td><td>n.d.</td></tr><tr><td>10</td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>FMN<sup>a</sup></td><td>0.8</td><td>0.8</td><td>-8.3</td><td>-5.8</td><td>-2.5</td><td>~1</td></tr><tr><td>FAD<sup>a</sup></td><td>0.8</td><td>12.5</td><td>-6.7</td><td>-5.6</td><td>-1.1</td><td>~1</td></tr><tr><td rowspan="4">SpnFADS:ADP</td><td rowspan="2">RF</td><td>0</td><td>0.3</td><td>-9.0</td><td>-17.2</td><td>8.2</td><td>~1</td></tr><tr><td>10</td><td>51.3<sup>b</sup></td><td>-5.9<sup>b</sup></td><td>-51.6<sup>b</sup></td><td>45.7<sup>b</sup></td><td>~2</td></tr><tr><td rowspan="2">FMN</td><td>0</td><td>1.5</td><td>-7.9</td><td>-12.8</td><td>4.9</td><td>&lt;1</td></tr><tr><td>10</td><td>5.9</td><td>-7.1</td><td>-28.6</td><td>21.5</td><td>&lt;1</td></tr><tr><td>SpnFAD:PP<sub>i</sub></td><td>FAD</td><td>0</td><td>n.d.</td><td>n.d.</td><td>n.d.</td><td>n.d.</td><td>n.d.</td></tr></tbody></table>

<figure><img src="image_6.png" /><figcaption>**Figure 3.** Hydrodynamic properties of *SpnFADS*. (A) Gel-filtration chromatogram of pure *SpnFADS* (20 µM) using a precalibrated Superdex™ 200 10/300 GL column (GE Healthcare) equilibrated with 50 mM Tris/HCl, pH 8.0, 150 mM NaCl, and calibrated with the Gel Filtration Calibration Kit LMW (GE Healthcare). The inset shows the percentage of the fractions obtained from the chromatogram. (B) 12% Native PAGE of pure *SpnFADS*. The two major bands are denoted by number; 1 for predicted oligomers and 2, the most populated, for the putative monomer. (C) 15% SDS-PAGE of *SpnFADS* without crosslinker and after incubation with the</figcaption></figure>

suggesting this small protein fraction contains high MW oligomers. Pure preparations of *SpnFADS* also migrated in two bands on native gel electrophoresis, being one of them highly populated (band 1, Fig. 3B). According to the gel filtration data this band 1 can be putatively assigned to the monomer.

SDS-PAGE resolved the protein in a single band corresponding to the monomer MW, but pre-incubation of the protein with bis(sulfosuccinimidyl) suberate (BS3) crosslinker stabilized an additional band with an apparent MW of 78.6 kDa that can be assigned to a dimer (Fig. 3C). Altogether these data suggest that purified *SpnFADS* can be found in a monomer-dimer equilibrium that is mainly shifted towards the monomeric form, with traces of higher MW assemblies.

We used atomic force microscopy (AFM) to shed light on the *SpnFADS* assembly in solution, as well as on the possible effect of ligands in the monomer-oligomer equilibrium (Fig. 4A-F). AFM topography images of free *SpnFADS* indicate a compact monomer (height ~7 ± 1 nm) (Fig. 4A). Incubation of the enzyme with different combinations of substrates and/or its reaction products produced some dimers and trimers, while ADP also promoted the formation of tetramers (Fig. 4 and Table 3). Two types of images were found for dimers. In one of them the two protomers dispose parallel in the same plane with the tallest module of one interacting protomer with the smallest of the other, leaving the rest of these molecules free of contact (Fig. 4B). In the other dimer, the
<figure><img src="image_7.png" /><figcaption>**Figure 4.** Image analysis of *SpnFADS* assemblies. Topography AFM images of a scanned area of 167 × 167 nm and Z-height profiles associated to the lines on the features of *SpnFADS* samples. (A) Monomer. (B) Dimer on the same plane. (C) Dimer with an overlapping on a monomer on the other. (D) Trimer from (B) showing an overlapped dimer interacting with a monomer. (E) Trimer from (C) showing overlapping of the dimer on a monomer. The different sizes of the original scanned images slightly change the roughness of the surface and the lateral resolution. This affects the size of the features in X-Y but not the Z-height measured in the profiles. (F) Tetramer. (G) Blender simulation of different assemblies according to the presence or absence of ligands and assembly of the three *SpnFADS* molecules in the crystal asymmetric unit (PDB ID: 3OP1).</figcaption></figure>

most frequent, both protomers are on the same plane establishing different angles and slightly overlapping one module of the first onto one of the modules of the other (Fig. 4C). Both dimers appear as intermediate steps in the formation of the two trimeric species visualized in the images (Fig. 4D and E), which situate an additional monomer over one of the monomers of the dimeric structures. Considering the disposition adopted by the three molecules present in the asymmetric unit of the *SpnFADS* crystal (Fig. 4G), where each protomer is situated in a different plane, the trimers in solution differ from such crystallographic association. Finally, tetramers form two layers of two overlapped monomers (Fig. 4F). Incubation of *SpnFADS* with ATP and FMN, ADP and FMN, or FAD and PPᵢ, produced the higher ratios of trimeric species, ~41%, while the RFK activity reactants, ATP and RF, generated 30% of dimeric molecules and only 16% of trimers. ATP by itself produced 37% of trimers, meanwhile
<table><caption>**Table 3.** Percentages of *SpnFADS* assemblies imaged by AFM. Samples of 0.5 µM *SpnFADS* were analyzed in 20 mM PIPES, pH 6.0, 2 mM DTT and 0.8 mM MgCl₂ at room temperature. When indicated, ADP, ATP or PPi were added at 250 µM, while oxidized FAD, FMN and RF were used at 50 µM. Features correspond to image units, while molecules refer to the amount of individual protein monomers in the corresponding image units. Error associated to percentage determination was ± 10%.</caption><thead><tr><th rowspan="2">Ligands</th><th rowspan="2">Units</th><th colspan="4">ASSOCIATION STATE</th></tr><tr><th>Monomers (%)</th><th>Dimers (%)</th><th>Trimers (%)</th><th>Tetramers (%)</th></tr></thead><tbody><tr><td rowspan="2">—</td><td>Features</td><td>97</td><td>3</td><td>—</td><td>—</td></tr><tr><td>Molecules</td><td>94</td><td>6</td><td>—</td><td>—</td></tr><tr><td rowspan="2">ADP</td><td>Features</td><td>60</td><td>7</td><td>13</td><td>20</td></tr><tr><td>Molecules</td><td>31</td><td>7</td><td>20</td><td>47</td></tr><tr><td rowspan="2">ATP</td><td>Features</td><td>70</td><td>12</td><td>18</td><td>—</td></tr><tr><td>Molecules</td><td>47</td><td>16</td><td>37</td><td>—</td></tr><tr><td rowspan="2">ATP, RF</td><td>Features</td><td>73</td><td>20</td><td>7</td><td>—</td></tr><tr><td>Molecules</td><td>54</td><td>30</td><td>16</td><td>—</td></tr><tr><td rowspan="2">ATP, FMN</td><td>Features</td><td>73</td><td>7</td><td>20</td><td>—</td></tr><tr><td>Molecules</td><td>50</td><td>9</td><td>41</td><td>—</td></tr><tr><td rowspan="2">ADP, FMN</td><td>Features</td><td>63</td><td>15</td><td>22</td><td>—</td></tr><tr><td>Molecules</td><td>40</td><td>19</td><td>41</td><td>—</td></tr><tr><td rowspan="2">FAD, PPi</td><td>Features</td><td>73</td><td>7</td><td>20</td><td>—</td></tr><tr><td>Molecules</td><td>50</td><td>9</td><td>41</td><td>—</td></tr></tbody></table>

ADP makes 20% of trimers and 47% of tetramers (Table 3). Incubation with the reactants of the FMNAT and FADpp activities resembles the same oligomeric yields. Since neither oxidized FMN nor FAD are able to bind to *SpnFADS* in the absence of adenine nucleotides, these data suggest that binding of adenine nucleotides (ATP and ADP) promotes the observed quaternary assemblies.

Discussion

Sequence alignments, structural comparisons¹⁰, ¹³, ³⁶ and ligand binding models (Supplementary Fig. S1, and Fig. 5A) suggested that *SpnFADS* contains all the features of prokaryotic FADSs contributing to substrates stabilization and catalysis. Accordingly, we show here that the enzyme catalyzes the RFK, FMNAT and FADpp activities, but we also highlight different requirements for catalysis and binding of ligands relative to other family members. Interestingly while transformation of RF into FMN by *SpnFADS* occurs similarly under both aerobic and reducing conditions, the FMNAT and FADpp activities require reduced flavins (Fig. 1). Reducing conditions hardly influence any of the *CaFADS* activities⁵, ³³ (Table 1), but some bifunctional FADSs as well as the monofunctional *HsRFK* (*Homo sapiens* RFK) are more active in the presence of high concentrations of reducing agents³⁴, ³⁵, ⁴², ⁴³. Moreover, studies of enriched samples of the FADS from *Bacillus subtilis* (*BsFADS*) early suggested specificity for reduced flavins in all its activities⁴.

Despite the differences in the redox requirements for the RFK and FMNAT activities between *SpnFADS* and *CaFADS*, both enzymes showed the highest efficiency for the synthesis of FMN and FAD compared to the other family enzymes reported up to now (Table 1)⁵, ³³, ³⁵, ⁴², with the lack of inhibition by the RF substrate for the RFK activity in *SpnFADS* (Fig. 1C) as a notable difference between them. Our ITC data indicate that *SpnFADS* does not bind any of its oxidized flavinic substrates when adenine nucleotides are not present (Supplementary Fig. S4), contrary to *CaFADS*⁵, ³³. Reducing conditions are not a requirement for the binding of adenine nucleotides (ATP and ADP) (Fig. 2), and the presence of the adenine nucleotide also induces binding of oxidized RF or FMN. Nevertheless, such ternary interactions in *SpnFADS* are slightly weaker than in *CaFADS*¹. Moreover, in the presence of the adenine nucleotide, oxidized FMN only binds at one site, which according to the enzyme capability to transform oxidized RF into FMN must be the RFK site (Table 2). The fact that oxidized RF does not bind to free *SpnFADS*, the lack of inhibition by RF, and the enzyme being catalytically efficient at low RF concentrations, suggest a dissimilar RFK mechanism between this enzyme and *CaFADS*⁴⁴. Furthermore, the enzymatic and ligand binding characterizations of *SpnFADS* altogether confirm that competent flavin binding for catalysis at the FMNAT site requires the flavin ring in its reduced state. In addition, *SpnFADS* needs higher concentrations of substrates than *CaFADS* to achieve maximal efficiency.

The need for reduced flavinic substrates in *BsFADS* was interpreted as the geometry of the active site only allowing the incorporation of isoalloxazine rings in the reduced bended conformations, contrary to the highly stable planar organization adopted by their oxidized states⁴, ¹⁷, ⁴⁵, ⁴⁶. Despite the high conservation grade in overall folding (Supplementary Fig. S1) and key active sites residues, *CaFADS* and *SpnFADS* are dissimilar in some regions involved in ligand binding (Fig. 5B and Supplementary Fig. S1). In the FMNAT module of *SpnFADS* helixes α310n and α5n break loops L4n and L6n, respectively. These regions form the external part of the cavity where substrates for the FMNAT and FADpp activities bind. The sites for the ribityl and isoalloxazine moieties are formed in both enzymes by hydrophobic and aromatic side-chains (Fig. 6A), whose differential spatial position might determine cavity volumes and, therefore, modes for flavin binding. The cavity appears slightly broader and deeper in *SpnFADS* than in *CaFADS*, although in *SpnFADS* the isoalloxazine binding site surface is rougher
<figure><img src="image_8.png" /><figcaption>**Figure 5.** The structure of *SpnFADS*. (A) Model (based on PDB ID: 3OP1) for the interaction of *SpnFADS* with its ligands. The FMNAT module is shown in green and the RFK in yellow. Flavin and adenine nucleotide ligands are in CPK colored sticks with carbons in white. Mg²⁺ ions are shown as grey balls. Modelling performed using PyMol. (B) Superposition of the structures of *SpnFADS* (green, PDB ID: 3OP1) and *CaFADS* (orange, PDB ID: 2X0K). Dissimilar regions are highlighted in blue and pink, respectively, including the numbers of implicated residues in *SpnFADS*.</figcaption></figure>

and shows a more negative surface electrostatic potential (Fig. 6B). The side chains of two Met residues, M54 and M118, form also part of the flavin ring cavity in *SpnFADS*, with M54 particularly prone to contact with the inner isoalloxazine ring (Fig. 6C). These Met residues are not conserved in most family members, which have instead hydrophobic residues¹³. Moreover, M54 forms part of a conserved sequence motif that is included only in few members such as *BsFADS*.

Structural differences are also observed when comparing the *SpnFADS* RFK module with the structures available for the *CaFADS* RFK module free and in complex with ligands (Figs 5 and 7)³⁶: i) D214 at loop L2c of *SpnFADS* salt bridges R190 at β1c; ii) loop L3c is 12 residues shorter; iii) the C-terminal is considerably shorter; and iv) an extra α helix, α1c, breaks loop L1c-Flap I. In addition, conformations of the conserved PTAN motif and of L1c-Flap I in free *SpnFADS* resemble those elements in the *CaFADS*:FMN:ADP-Mg²⁺ complex. Such organization leaves the adenine nucleotide binding site open and the PTAN motif ready to stabilize the ligand (Fig. 7), contrary to free *CaFADS*³⁶. Thus, the weaker affinity of *SpnFADS* for the RFK activity products might be related with the lack of large conformational changes at the active site for substrates association and products dissociation during catalysis. These traits surely also contribute to the lack of inhibition by the RF substrate in *SpnFADS*.

Finally, substrates and products of *SpnFADS* induce formation of dimers, tetramers and particularly trimers with head-to-tail contacts between the RFK module and the FMNAT module of different protomers (Table 3, Fig. 4). None of these organizations resembles the *CaFADS* dimer-of-trimers and in general the *SpnFADS* monomer appears to be majority⁴⁷. The PDB:PISA server predicts dimeric assemblies with head-to-tail fits (Supplementary Fig. S5 and Table S1), but it does not predict any higher oligomer as stable in solution. Differences in oligomerization profiles might relate with differences in protein oligomerization regions. In fact, a long L3c, only expected for *Corynebacterium* and *Mycobacterium* species, is predicted as a determinant to stabilize the dimer-of-trimers⁴⁷. The shorter loop in *SpnFADS* agrees with this hypothesis, being less prone to oligomerization and establishing different contacts among protomers, which might also contribute to the dissimilar binding, kinetic and inhibition properties among species. The dimer-of-trimers in *CaFADS* envisaged a controlled interplay between the regions involved in the formation of the RFK and/or FMNAT catalytic sites of different protomers as well as in flavin delivery to client apoflavoproteins⁴⁸, ⁴⁹. However, such mechanism might not be required for *SpnFADS*, since in this enzyme the main regulation mechanism for the FAD synthesis seems to be the reducing environment.

In conclusion, *SpnFADS* is a prokaryotic bifunctional FADS that folds in two modules, the C-terminal one being responsible for the RFK activity and the N-terminal one accounting for the FMNAT and FADpp activities. Three main characteristics distinguish *SpnFADS* from the previously well characterized *CaFADS*: i) its FMNAT module only binds and transforms flavins in their reduced state, ii) its RFK activity does not exhibit inhibition by excess of RF substrate and iii) only dimers and trimers appear to be formed during catalysis. Considering the low number of members of the prokaryotic FADS family hitherto characterized, a large dissimilarity of optimal
<figure><img src="image_9.png" /><figcaption>**Figure 6.** Comparative analysis of the *SpnFADS* and *CaFADS* active sites in the *FMNAT* module. (A) Site for ligands binding in *SpnFADS* (green) and *CaFADS* (orange). Hydrophobic and aromatic side chains are shown in sticks. (B) Surface electrostatic potential for the isoalloxazine binding cavity in *SpnFADS* (left panel) and *CaFADS* (right panel). Arrows point to differences in sign of surface electrostatic potential at the *SpnFADS* and *CaFADS* isoalloxazine sites. (C) Relative disposition of M54 and M118 regarding the isoalloxazine moiety of the flavin substrate. In all panels, FMN and AMP have been modelled on the crystals structures of *SpnFADS* (PDB ID: 3OP1) and *CaFADS* (PDB ID: 2X0K) and represented in CPK sticks with carbons in grey¹⁰.</figcaption></figure>
<figure><img src="image_10.png" /><figcaption>**Figure 7.** Comparative analysis of the RFK sites of *SpnFADS*, *CaFADS* and the *CaFADS:FMN:ADP-Mg*²⁺ ternary complex. (A) Cartoon representation of the RFK domains of *SpnFADS* (green, PDB ID: 3OP1), *CaFADS* (orange, PDB ID: 2X0K) and ternary complex of the RFK module of *CaFADS* (blue, PDB ID: 5A8A). L1c-FlapI is highlighted in red, pink and yellow, respectively. Residues of the PTAN motif are in CPK sticks. The FMN and ADP ligands bound to the *CaFADS* RFK module are in white CPK and Mg²⁺ is included as a green sphere. (B) Detail of the PTAN motif in *SpnFADS* and superposition of the PTAN motives of the three structures. (C) Superposition of the RFK domain of *SpnFADS* (green) and *CaFADS* (orange). The salt bridge between D214 at L2c and R190 at β1c of *SpnFADS* is in black dashes. Equivalent residues in *CaFADS*, D215 and T190, are drawn in sticks. Loops L3c are colored in blue (*SpnFADS*) and pink (*CaFADS*).</figcaption></figure>

conditions for catalysis is envisaged. FADSs have a major functional role providing flavin cofactors to the cellular flavoproteome while also promoting flavins and flavoproteins homeostasis in all living organisms. In this context substrate inhibition, the use of reduced substrates or the formation of oligomeric assemblies might reflect different regulation mechanisms used by the enzyme to withstand relatively different environments. Such, differences among family members provide with a framework to the design of selective compounds targeting FADS for the treatment of diverse infectious diseases.

## Methods

**Cloning, Overexpression and Purification of SpnFADS.** The DNA sequence encoding *SpnFADS* (accession number ADM84672, SPAP_1083), containing the *NdeI* and *XhoI* restriction sites and codons optimized for its expression in *E. coli*, was synthetically produced by Gen-Script. This DNA sequence was cloned into a modified pET-15b vector (*Novagen*) that contains an N-terminal His₆-Tag sequence followed by a PreScission Protease site upstream of the *NdeI* and *XhoI* cloning sites. The final construct was fully sequenced at *Sistemas Genómicos* (www.sistemasgenomicos.com) and transformed into the *E. coli* strain BL21 Star™ (DE3) (*Invitrogen*). Transformed cells were grown at 37°C in 2xTY medium (1.6% (w/v) tryptone, 1% (w/v) yeast extract, and 0.5% (w/v) NaCl) with 100 µg/mL of ampicillin until cells reached an OD₆₀₀ of 0.6. Then the culture was cooled down to 18°C and protein expression was induced by addition of 1 mM isopropyl β-D-1-thiogalactopyranoside. After 20 h of incubation at 18°C, cells were harvested by centrifugation and suspended in 20 mM sodium phosphate, pH 7.4, 500 mM NaCl and 10 mM imidazole, containing lysozyme (1 mg/mL), DNase (0.1 mg/mL) and protease inhibitors (0.2 mM Phenylmethanesulfonyl fluoride (PMSF) and 10 mM benzamidine). After 30 min of incubation at 37°C, cells were broken by sonication and centrifuged. The supernatant containing the soluble protein was loaded into a His-Trap affinity column (HisTrap HP, GE Healthcare) and the washing flow was discarded. Protein was eluted applying a gradient from 10 to 500 mM imidazole in 20 mM sodium phosphate, pH 7.4, 500 mM
NaCl. Buffer was exchanged to 25 mM Tris/HCl, pH 7.4, 150 mM NaCl, using a HiPrep Desalting Column (GE Healthcare). The His₆-Tag was removed by 48 h-incubation with PreScission protease (GE Healthcare) at 4 °C, in ratio 1:5 (w/w), and then the protein was loaded into the HisTrap HP and GSTrap 4B connected columns (GE Healthcare) to eliminate both the remained His₆-Tagged SpnFADS and the PreScission protease (a GST fused product). The unbound fraction was recovered and further purified by size exclusion chromatography using a Superdex™ 200 10/300 GL column (GE Healthcare) previously equilibrated with 20 mM PIPES, pH 7.0, 150 mM NaCl. Protein purity was assessed by 15% SDS-PAGE and NaCl was removed by dialysis in 25 mM Tris/HCl, pH 7.5 or 20 mM PIPES, pH 7.0 depending on later use. Pure protein solution was stored at -80 °C.

**Spectroscopic Analysis.** UV-visible spectra were recorded in a Cary 100 spectrophotometer (Varian) and SpnFADS was quantified using the theoretical ε<sub>279nm</sub> = 28.88 mM⁻¹·cm⁻¹ and a molecular weight of 34,521 Da (ProtParam). Circular dichroism (CD) spectra were recorded with a Chirascan spectropolarimeter (Applied Photophysics Ltd.) at 25 °C. Samples containing ~5 and ~20 µM SpnFADS in 25 mM Tris/HCl, pH 7.5 were used in the far-UV (cuvette path length, 0.1 cm) and near-UV CD (cuvette path length, 0.4 cm), respectively.

**Qualitative and quantitative analysis of the activities of SpnFADS.** RFK, FMNAT and FADpp activities of SpnFADS were qualitatively assayed by separating the different flavins present in the reaction mixtures using thin layer chromatography on Silica Gel SIL-G-25 plates as previously described¹³. Reaction mixtures containing 50 µM flavin (either RF, FMN or FAD), 0.2 mM ATP or PP<sub>i</sub>, 0.8 or 10 mM MgCl<sub>2</sub>, and 0.5 µM SpnFADS in 150 µL of 50 mM Tris/HCl, pH 8.0, were incubated in the presence of different sodium dithionite concentrations (0–16 mM) at 25 °C during 30 min. The reaction was stopped by boiling the samples at 100 °C for 5 min. Flavin mixtures containing 50 µM of each flavin, RF, FMN and FAD, were included as standards. Flavin TLC spots were examined by monitoring their fluorescence under UV light.

RFK and FMNAT activities of SpnFADS were quantitatively measured at 25 °C in 20 mM PIPES, pH 7.0, 0–10 mM MgCl<sub>2</sub> as previously described⁵, ³³, and when indicated 0–10 mM sodium dithionite was also added. Reaction mixtures (total volume 500 µL) containing 0.5–15 µM RF or 1–55 µM FMN, and 10–400 µM ATP were pre-incubated at 25 °C until reaction was initiated by addition of ~25 nM enzyme. In all cases, reaction was stopped after 1 min of incubation through protein denaturation by boiling the sample for 5 min. Precipitated protein was removed from the solution by centrifugation at 18,000 x g during 10 min. The flavin composition of supernatants was analyzed using an Alliance HPLC system (Waters) equipped with a 2707 autosampler and an HSST3 column (4.6 × 150 mm, 3.5 µm, Waters) preceded by a precolumn (4.6 × 20 mm, 3.5 µm, Waters) as previously described⁵, ³³. The produced amounts of FMN or FAD were quantified using their corresponding standard curves acquired under the same conditions. The kinetic data obtained for one substrate at saturating concentrations of the second substrate (as nmol of flavin transformed *per* min) were interpreted using the Michaelis-Menten kinetic model, obtaining *k*<sub>cat</sub> and *K*<sub>m</sub> values with errors of ±10%. Controls for transformations of flavins in the absence of the enzyme under all reaction conditions showed only negligible processes.

Additional experiments were performed to elucidate whether reduced flavins are required for the FMNAT activity, or if the reducing environment is a requirement for the enzyme. Reaction mixtures containing 500 µM ATP, 15 µM FMN, 2 µM 5-deazaRF (5-dRF) and 0.8 mM MgCl₂ in a final volume of 1 mL were deoxygenated during 1 h through Argon bubbling in the dark (to avoid FMN degradation). Once oxygen was displaced from the samples, previously deoxygenated SpnFADS was added at a final concentration of 100 nM. Then samples were illuminated during 15 s to photo-reduce FMN (the triplet state of 5-dRF, obtained by light-excitation, quickly evolves to a radical state that returns to the fundamental state by donating an electron to FMN). After 3 min of additional incubation samples were boiled for 5 min (to stop the reaction) and then centrifuged. Flavins present in the supernatant were identified by HPLC as we indicated above.

**High-Sensitivity Isothermal Titration Calorimetry (ITC).** ITC measurements were performed using an AutoITC200 calorimeter (MicroCal-Malvern) thermostated at 20 °C following the procedures previously reported⁵, ³³. Typically, 300 µM nucleotide (ATP/ADP) or 150–200 µM flavin (RF/FMN/FAD) solutions were used to titrate ~20 µM SpnFADS solutions. When appropriate, ternary titrations were performed by adding flavin solutions to a mixture of the protein pre-incubated with 400 µM ADP. Ligands and SpnFADS were dissolved in 20 mM PIPES pH 7.0, with 0, 0.8 or 10 mM MgCl₂, and degassed prior to titration. Up to 19 injections of 2 µL were programmed with enough time spacing for the signal to recover the baseline. Similar protein concentrations in the calorimetric cell were employed in all experiments to guarantee that the oligomerization state of the protein is the same at any time. The association constant (*K*<sub>a</sub>), the enthalpy change (ΔH) and the stoichiometry (N) were estimated through non-linear regression of the experimental data using a model for one or two independent binding sites; the regression was implemented in Origin 7.0 (OriginLab). The dissociation constant (*K*<sub>d</sub>), the free energy change (ΔG), and the entropy change (ΔS) were obtained from basic thermodynamic relationships. Errors considered in the measured parameters (±20% in *K*<sub>d</sub> and ±0.3 kcal/mol in ΔH and -TΔS) were taken larger than the standard deviation between replicates and the numerical error after the fitting analysis.

Titrations performed to elucidate whether SpnFADS was able to bind reduced FMN and FAD were performed using a high precision VP-ITC system (MicroCal LLC), and experimental sets were similar to those for the experiments in the AutoITC200 calorimeter. Thus, 25 µM SpnFADS were titrated at 20 °C with FMN ≈ 200 µM or FAD ≈ 300 µM. Both the protein and the ligands were dissolved in 20 mM PIPES, pH 7.0, 0.8 mM MgCl₂, 4 mM sodium dithionite, 10 mM glucose, 0.5 U/mL glucose oxidase, and degassed. Each titration was initiated by a 4 µL injection followed by 25–28 stepwise injections of 10 µL. Glucose and glucose oxidase were added to avoid oxidation of sodium dithionite by molecular oxygen in order to avoid unspecific oxidation heat that could mask the one corresponding to the flavin:protein interaction.
**Evaluation of the hydrodynamic and quaternary organization properties.** Stabilization of the quaternary organizations of *SpnFADS* was initially evaluated by gel filtration onto a Superdex™ 200 10/300 GL column (*GE Healthcare*), equilibrated with 50 mM Tris/HCl, pH 8.0, 150 mM NaCl, and calibrated with the Gel Filtration Calibration Kit LMW (*GE Healthcare*)$^{47}$. Oligomer formation was also analyzed by incubation of a pure fraction of *SpnFADS* in 50 mM KHPO$_4$, pH 8.0 with BS3 crosslinker (*Thermo Scientific*) for 2 h at 4 °C. The cross reaction was stopped adding the quenching solution (Tris/HCl 500 mM (10x), pH 8.0) up to 50 mM and the different oligomeric states of the protein were visualized on a 15% SDS-PAGE gel.

For single molecule visualization of *SpnFADS* preparations, atomic force microscopy imaging was performed using a Multimode 8 AFM system (Bruker). Images were taken using the Tapping Mode with calibrated soft V-shaped silicon nitride cantilevers carrying pyramidal 2 nm ultrasharp tips characterized by spring constants ranged from 0.01 to 0.03 N/m and nominal resonance frequencies of 7 to 15 kHz (Bruker Probes). For samples containing ligands, the mixtures were pre-incubated at 4 °C for 10 min under mild stirring to favor complex formation prior to immobilization$^{47}$. ADP, ATP and $PP_1$ were added at 250 µM, while FAD, FMN and RF were used at 50 µM. *SpnFADS* solutions (~0.5 µM) were immobilized on small fresh exfoliated muscovite mica pieces (*Electron Microscopy Sciences*) for 10 min at room temperature. The sample and holder were introduced in a liquid cell previously cleaned with 20% isopropanol and Millipore ultrapure water. AFM measurements were conducted in 20 mM PIPES, pH 6.0, freshly prepared 2 mM dithiothreitol (DTT, Sigma-Aldrich) and 0.8 mM $MgCl_2$ at room temperature. Raw images were processed using the WSxM freeware$^{50}$. Percentages of different features were calculated taking several sample areas with several scan sizes, zooming with the WSxM software functions without losing relevant information and discarding artifacts$^{51}$. A considerable amount of features were collected for the statistical population study. Total protein percentages were also calculated considering the number of molecules involved in each type of association. Errors were calculated from the dispersion of results in the analysis of different AFM images corresponding to different areas of the sample.

**Structural prediction of enzyme-ligand interactions.** Models for ligand binding at the RFK and FMNAT sites of *SpnFADS* were produced comparing its structure (chain B of the crystal structure PDB ID: 3OP1) with the previously reported *CaFADS* crystal structures (for the RFK module) or models (FMNAT module) in complex with their oxidized substrates$^{10,36}$. The PyMol 0.99$^{52}$ software was used for complex modelling, structural analysis and figures production.

**Data availability statement.** All data generated or analyzed during this study are included in this published article (and in its Supplementary Information files).

**References**

1. Frago, S., Velázquez-Campoy, A. & Medina, M. The puzzle of ligand binding to Corynebacterium ammoniagenes FAD synthetase. *J Biol Chem* **284**, 6610–6619, doi: M808142200 (2009).

2. Mack, M., van Loon, A. P. & Hohmann, H. P. Regulation of riboflavin biosynthesis in Bacillus subtilis is affected by the activity of the flavokinase/flavin adenine dinucleotide synthetase encoded by ribC. *J Bacteriol* **180**, 950–955 (1998).

3. Manstein, D. J. & Pai, E. F. Purification and characterization of FAD synthetase from Brevibacterium ammoniagenes. *J Biol Chem* **261**, 16169–16173 (1986).

4. Kearney, E. B., Goldenberg, J., Lipsick, J. & Perl, M. Flavokinase and FAD synthetase from Bacillus subtilis specific for reduced flavins. *J Biol Chem* **254**, 9551–9557 (1979).

5. Serrano, A., Frago, S., Velázquez-Campoy, A. & Medina, M. Role of key residues at the flavin mononucleotide (FMN):adenylyltransferase catalytic site of the bifunctional riboflavin kinase/flavin adenine dinucleotide (FAD) Synthetase from Corynebacterium ammoniagenes. *Int J Mol Sci* **13**, 14492–14517, doi: 10.3390/ijms131114492 (2012).

6. Brizio, C., Otto, A., Brandsch, R., Passarella, S. & Barile, M. A protein factor of rat liver mitochondrial matrix involved in flavinylation of dimethylglycine dehydrogenase. *Eur J Biochem* **267**, 4346–4354 (2000).

7. Barile, M., Brizio, C., Valenti, D., De Virgilio, C. & Passarella, S. The riboflavin/FAD cycle in rat liver mitochondria. *Eur J Biochem* **267**, 4888–4900 (2000).

8. Yruela, I., Arilla-Luna, S., Medina, M. & Contreras-Moreira, B. Evolutionary divergence of chloroplasts FAD synthetase proteins. *BMC Evol Biol* **10**, 311 (2010).

9. Mashhadi, Z., Xu, H., Grochowski, L. L. & White, R. H. Archaeal RibL: a new FAD synthetase that is air sensitive. *Biochemistry* **49**, 8748–8755, doi: 10.1021/bi100817q (2010).

10. Herguedas, B., Martinez-Julvez, M., Frago, S., Medina, M. & Hermoso, J. A. Oligomeric state in the crystal structure of modular FAD synthetase provides insights into its sequential catalysis in prokaryotes. *J Mol Biol* **400**, 218–230, doi: S0022-2836(10)00497-3 (2010).

11. Wang, W., Kim, R., Yokota, H. & Kim, S. H. Crystal structure of flavin binding to FAD synthetase of Thermotoga maritima. *Proteins* **58**, 246–248, doi: 10.1002/prot.20207 (2005).

12. Karthikeyan, S. et al. Crystal structure of human riboflavin kinase reveals a beta barrel fold and a novel active site arch. *Structure* **11**, 265–273, doi: S0969212603000248 (2003).

13. Frago, S., Martínez-Julvez, M., Serrano, A. & Medina, M. Structural analysis of FAD synthetase from Corynebacterium ammoniagenes. *BMC Microbiol* **8**, 160, doi: 1471-2180-8-160 (2008).

14. Huerta, C., Borek, D., Machius, M., Grishin, N. V. & Zhang, H. Structure and mechanism of a eukaryotic FMN adenylyltransferase. *J Mol Biol* **389**, 388–400, doi: S0022-2836(09)00452-5 (2009).

15. Krupa, A., Sandhya, K., Srinivasan, N. & Jonnalagadda, S. A conserved domain in prokaryotic bifunctional FAD synthetases can potentially catalyze nucleotide transfer. *Trends Biochem Sci* **28**, 9–12, doi: S0968000402000099 (2003).

16. Macheroux, P., Kappes, B. & Ealick, S. E. Flavogenomics-a genomic and structural view of flavin-dependent proteins. *Fems J* **278**, 2625–2634, doi: 10.1111/j.1742-4658.2011.08202.x (2011).

17. Massey, V. The chemical and biological versatility of riboflavin. *Biochem Soc Trans* **28**, 283–296 (2000).

18. Parsons, H. G. & Dias, V. C. Intramitochondrial fatty acid metabolism: riboflavin deficiency and energy production. *Biochem Cell Biol* **69**, 490–497 (1991).

19. Daff, S. An appraisal of multiple NADPH binding-site models proposed for cytochrome P450 reductase, NO synthase, and related diflavin reductase systems. *Biochemistry* **43**, 3929–3932, doi: 10.1021/bi0360408 (2004).

20. Wang, M. et al. Three-dimensional structure of NADPH-cytochrome P450 reductase: prototype for FMN- and FAD-containing enzymes. *Proc Natl Acad Sci USA* **94**, 8411–8416 (1997).

21. Susin, S. A. et al. Molecular characterization of mitochondrial apoptosis-inducing factor. *Nature* **397**, 441–446 (1999).
22. White, H. B. & Merrill, A. H. Riboflavin-binding proteins. *Annu Rev Nutr* **8**, 279–299, doi:10.1146/annurev.nu.08.070188.001431 (1988).

23. Forneris, F., Binda, C., Vanoni, M. A., Battaglioli, E. & Mattevi, A. Human histone demethylase LSD1 reads the histone code. *J Biol Chem* **280**, 41360–41365, doi:10.1074/jbc.M509549200 (2005).

24. Myllykallio, H. et al. An alternative flavin-dependent mechanism for thymidylate synthesis. *Science* **297**, 105–107, doi:10.1126/science.1072113 (2002).

25. Nishimasu, H. et al. Atomic structure of a folate/FAD-dependent tRNA T54 methyltransferase. *Proc Natl Acad Sci USA* **106**, 8180–8185, doi:10.1073/pnas.0901330106 (2009).

26. Gross, E., Kastner, D. B., Kaiser, C. A. & Fass, D. Structure of Ero1p, source of disulfide bonds for oxidative protein folding in the cell. *Cell* **117**, 601–610 (2004).

27. Demarsy, E. & Fankhauser, C. Higher plants use LOV to perceive blue light. *Curr Opin Plant Biol* **12**, 69–74, doi:10.1016/j.pbi.2008.09.002 (2009).

28. Kraiselburd, I. et al. A LOV protein modulates the physiological attributes of Xanthomonas axonopodis pv. citri relevant for host plant colonization. *PLoS One* **7**, e38226, doi:10.1371/journal.pone.0038226 (2012).

29. Masuda, S. Light Detection and Signal Transduction in the BLUF Photoreceptors. *Plant Cell Physiol* **55**, 1858, doi:10.1093/pcp/pcu103 (2014).

30. Mehlhorn, J. et al. Light-Induced Rearrangement of the β5 Strand in the BLUF Photoreceptor SyPixD (Slr1694). *J Phys Chem Lett* **6**, 4749–4753, doi:10.1021/acs.jpclett.5b02245 (2015).

31. Tanwar, M., Nahar, S., Gulati, S., Veetil, S. K. & Kateriya, S. Molecular determinant modulates thermal recovery kinetics and structural integrity of the bacterial BLUF photoreceptor. *FEBS Lett* **590**, 2146–2157, doi:10.1002/1873-3468.12227 (2016).

32. Wang, W., Kim, R., Jancarik, J., Yokota, H. & Kim, S. H. Crystal structure of a flavin-binding protein from Thermotoga maritima. *Proteins* **52**, 633–635, doi:10.1002/prot.10353 (2003).

33. Serrano, A. et al. Key residues at the riboflavin kinase catalytic site of the bifunctional riboflavin kinase/FMN adenyllyltransferase from Corynebacterium ammoniagenes. *Cell Biochem Biophys* **65**, 57–68, doi:10.1007/s12013-012-9403-9 (2013).

34. Grill, S., Busenbender, S., Pfeiffer, M., Kohler, U. & Mack, M. The bifunctional flavokinase/flavin adenine dinucleotide synthetase from Streptomyces davawensis produces inactive flavin cofactors and is not involved in resistance to the antibiotic roseoflavin. *J Bacteriol* **190**, 1546–1553, doi:JB.01586-07 (2008).

35. Matern, A., Pedrolli, D., Großhennig, S., Johansson, J. & Mack, M. Uptake and Metabolism of Antibiotics Roseoflavin and 8-Demethyl-8-Aminoriboflavin in Riboflavin-Auxotrophic Listeria monocytogenes. *J Bacteriol* **198**, 3233–3243, doi:10.1128/JB.00388-16 (2016).

36. Herguedas, B. et al. Structural insights into the synthesis of FMN in prokaryotic organisms. *Acta Crystallogr D Biol Crystallogr* **71**, 2526–2542, doi:10.1107/S1399004715019641 (2015).

37. Serrano, A., Ferreira, P., Martinez-Júlvez, M. & Medina, M. The prokaryotic FAD synthetase family: a potential drug target. *Curr Pharm Des* **19**, 2637–2648 (2013).

38. O'Brien, K. L. et al. Burden of disease caused by Streptococcus pneumoniae in children younger than 5 years: global estimates. *Lancet* **374**, 893–902, doi:10.1016/S0140-6736(09)61204-6 (2009).

39. Mitchell, A. M. & Mitchell, T. J. Streptococcus pneumoniae: virulence factors and variation. *Clin Microbiol Infect* **16**, 411–418, doi:10.1111/j.1469-0691.2010.03183.x (2010).

40. Musher, D. M. Infections caused by Streptococcus pneumoniae: clinical spectrum, pathogenesis, immunity, and treatment. *Clin Infect Dis* **14**, 801–807 (1992).

41. Raussens, V., Ruysschaert, J. M. & Goormaghtigh, E. Protein concentration is not an absolute prerequisite for the determination of secondary structure from circular dichroism spectra: a new scaling method. *Anal Biochem* **319**, 114–121 (2003).

42. Pedrolli, D. B. et al. The antibiotics roseoflavin and 8-demethyl-8-amino-riboflavin from Streptomyces davawensis are metabolized by human flavokinase and human FAD synthetase. *Biochem Pharmacol* **82**, 1853–1859, doi:S0006-2952(11)00680-0 (2011).

43. Langer, S., Hashimoto, M., Hobl, B., Mathes, T. & Mack, M. Flavoproteins are potential targets for the antibiotic roseoflavin in Escherichia coli. *J Bacteriol* **195**, 4037–4045, doi:10.1128/JB.00646-13 (2013).

44. Efimov, I., Kuusk, V., Zhang, X. & McIntire, W. S. Proposed steady-state kinetic mechanism for Corynebacterium ammoniagenes FAD synthetase produced by Escherichia coli. *Biochemistry* **37**, 9716–9723, doi:10.1021/bi972817j (1998).

45. Massey, V. Introduction: flavoprotein structure and mechanism. *FASEB J* **9**, 473–475 (1995).

46. Heelis, P. F. The photophysical and photochemical properties of flavins (isoalloxazines). *Chem Soc Rev* **11**, 15–39 (1982).

47. Marcuello, C., Arilla-Luna, S., Medina, M. & Lostao, A. Detection of a quaternary organization into dimer of trimers of Corynebacterium ammoniagenes FAD synthetase at the single-molecule level and at the in cell level. *Biochim Biophys Acta* **1834**, 665–676, doi:S1570-9639(12)00293-2 (2013).

48. Serrano, A. et al. Quaternary organization in a bifunctional prokaryotic FAD synthetase: Involvement of an arginine at its adenyllyltransferase module on the riboflavin kinase activity. *Biochim Biophys Acta* **1854**, 897–906, doi:10.1016/j.bbpap.2015.03.005 (2015).

49. Giancaspero, T. A. et al. Remaining challenges in cellular flavin cofactor homeostasis and flavoprotein biogenesis. *Front Chem* **3**, 30, doi:10.3389/fchem.2015.00030 (2015).

50. Horcas, I. et al. WSXM: A software for scanning probe microscopy and a tool for nanotechnology. *Review of Scientific Instruments* **78**, doi:Artn 013705 doi:10.1063/1.2432410 (2007).

51. Lostao, A., Peleato, M. L., Gomez-Moreno, C. & Fillat, M. F. Oligomerization properties of FurA from the cyanobacterium *Anabaena* sp. PCC 7120: direct visualization by *in situ* atomic force microscopy under different redox conditions. *Biochim Biophys Acta* **1804**, 1723–1729, doi:S1570-9639(10)00100-7 (2010).

52. Delano, W. L. They PyMOL molecular graphics system. DeLano Scientific, San Carlos, CA, USA, http://www.pymol.org (2002).

## Acknowledgements

This work has been supported by the Spanish Ministry of Economy and Competitiveness (MINECO) [BIO2013-42978-P and BIO2016-75183-P to M.M.; CTQ2013-44367-C2-2-P and BFU2016-75633-P to R.H.-G.], and the Government of Aragón-FEDER [B18 and B89]. A.L., A.V.-C. and R.H-G. thank ARAID for financial support. E.L. received a fellowship from Bruker, and M.S. and C.M. from DGA. The authors thank also J.L Díez and I. Echániz for technical support.

## Author Contributions

M.S., E.L.-N. and A.S. designed and performed the experiments. R.H.G. helped in some experiment designs and characterization. A.V.-C. supervised the ITC experiments and their analysis. A.L. designed, analysed and supervised the AFM experiments and C.M. preformed them. M. M.-J. and M.M. conceived the study, its global design and coordination, and supervised the experimental work and data analysis. M.S., A.S., M. M.-J and M.M. wrote the main manuscript. All authors read and approved the final manuscript.
**Additional Information**

**Supplementary information** accompanies this paper at doi:10.1038/s41598-017-07716-5

**Competing Interests:** The authors declare that they have no competing interests.

**Publisher's note:** Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

<figure><img src="image_11.png" /></figure>

**Open Access** This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons license, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons license and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this license, visit http://creativecommons.org/licenses/by/4.0/.

© The Author(s) 2017