# Molecular characterization of the missing electron pathways for butanol synthesis in Clostridium acetobutylicum

Céline Foulquier, Antoine Rivière, Mathieu Heulot, Suzanna dos Reis, Caroline Perdu, Laurence Girbal, Mailys Pinault, Simon Dusséaux, Minyeong Yoo, Philippe Soucaille, et al.

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

## To cite this version:

Céline Foulquier, Antoine Rivière, Mathieu Heulot, Suzanna dos Reis, Caroline Perdu, et al. Molec- ular characterization of the missing electron pathways for butanol synthesis in Clostridium aceto- butylicum. Nature Communications, 2022, 13, 10.1038/s41467-022-32269-1 . hal-03811719

HAL Id: hal-03811719

https://hal.science/hal-03811719

Submitted on 12 Oct 2022

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_2.png" /></figure>
<figure><img src="image_3.png" /></figure>

# Molecular characterization of the missing electron pathways for butanol synthesis in *Clostridium acetobutylicum*

Received: 6 September 2021

Accepted: 22 July 2022

Published online: 10 August 2022

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

Check for updates

Céline Foulquier<fnref n="1" />,³, Antoine Rivière<fnref n="1" />,³, Mathieu Heulot<fnref n="1" />, Suzanna Dos Reis<fnref n="1" />, Caroline Perdu<fnref n="1" />, Laurence Girbal<fnref n="1" />, Mailys Pinault<fnref n="1" />, Simon Dusséaux<fnref n="1" />, Minyeong Yoo<fnref n="1" />,², Philippe Soucaille<fnref n="1" />,²✉ & Isabelle Meynial-Salles<fnref n="1" />

*Clostridium acetobutylicum* is a promising biocatalyst for the renewable production of *n*-butanol. Several metabolic strategies have already been developed to increase butanol yields, most often based on carbon pathway redirection. However, it has previously demonstrated that the activities of both ferredoxin-NADP⁺ reductase and ferredoxin-NAD⁺ reductase, whose encoding genes remain unknown, are necessary to produce the NADPH and the extra NADH needed for butanol synthesis under solventogenic conditions. Here, we purify, identify and partially characterize the proteins responsible for both activities and demonstrate the involvement of the identified enzymes in butanol synthesis through a reverse genetic approach. We further demonstrate the yield of butanol formation is limited by the level of expression of *CA_C0764*, the ferredoxin-NADP⁺ reductase encoding gene and the *bcd* operon, encoding a ferredoxin-NAD⁺ reductase. The integration of these enzymes into metabolic engineering strategies introduces opportunities for developing a homobutanologenic *C. acetobutylicum* strain.

*Clostridium acetobutylicum* is a Gram-positive, spore-forming anaerobic bacterium capable of converting various sugars and polysaccharides to organic acids (acetate and butyrate) and solvents (acetone, butanol, and ethanol). Due to its importance in the industrial production of the bulk chemicals acetone and butanol<fnref n="1" />⁻³ and its potential use in the production of *n*-butanol, a promising biobased liquid fuel with several advantages over ethanol⁴,⁵, much research has focused on (i) understanding the regulation of solvent formation⁶⁻¹⁵ and (ii) metabolically engineering this microorganism to produce high yields of alcohols¹⁶⁻¹⁸.

Using a global system biology approach to the characterization of the solventogenic metabolism of a phosphate-limited chemostat culture of *C. acetobutylicum*, the six steps involved in the conversion of acetyl-CoA to butanol (Fig. 1) were previously characterized¹⁹: the main enzyme responsible for crotonyl-CoA reduction to butyryl-CoA was

shown to be the BCD complex (encoded by *bcd*, *etfB*, and *etfA*, the *bcd* operon), a bifurcating enzyme consuming two moles of NADH and producing one mole of reduced ferredoxin, and the last two steps of butanol production were shown to be catalyzed by AdhE1 through its NADH-dependent aldehyde dehydrogenase activity and by BdhB, BdhC, and BdhA through their NADPH-dependent butanol dehydrogenase activity. These results had a strong impact on electron flux distribution, as it was demonstrated that both ferredoxin-NADP⁺ reductase and ferredoxin-NAD⁺ reductase activities were necessary to produce the NADPH and the extra NADH needed for butanol synthesis from acetyl-CoA (Fig. 2)¹⁹. Although the activities of these enzymes were previously detected in *C. acetobutylicum*⁶,²⁰, the encoding genes remained unknown¹⁹,²¹. Ferredoxin-NADP⁺ reductase enzymes (FNOR) (EC 1.18.1.3) are distributed over a variety of aerobic organisms from prokaryotes to eukaryotes, especially in plants²², but have never been

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

e-mail: soucaille@insa-toulouse.fr

<footnote n="1">¹TBI, Université de Toulouse, CNRS, INRAE, INSA, Toulouse, France. ²BBSRC/EPSRC Synthetic Biology Research Centre, School of Life Sciences, Centre for Biomolecular Sciences, University of Nottingham, Nottingham, UK. ³These authors contributed equally: Céline Foulquier, Antoine Rivière.</footnote>
<figure><img src="image_6.png" /></figure>

purified and characterized from any clostridial species. In contrast, ferredoxin-dependent NAD⁺ reduction coupled to proton export (Rnf) and ferredoxin-dependent transhydrogenases (Nfn) have been characterized from a molecular perspective in several clostridial species but no homologs have been found in *C. acetobutylicum*²³⁻²⁵. Moreover, *C. acetobutylicum* was shown to be unable to produce NADPH by the

oxidative pentose-phosphate pathway, as a gene encoding a glucose-6-P dehydrogenase, the first and key enzyme of this pathway, was missing¹⁹,²⁶.

In this work, we purify and identify the proteins responsible for the ferredoxin-NAD⁺ and ferredoxin-NADP⁺ reductase activities in *C. acetobutylicum* under solventogenic conditions. We further confirm
**Fig. 1** | **Central metabolism of Clostridium acetobutylicum.** The green box indicates the primary products under acidogenic conditions, whereas the red box indicates the primary products under solventogenic conditions. The letters in red and italics indicate the corresponding genes. *ack* acetate kinase, *adc* acetoacetate decarboxylase, *adhE1* aldehyde dehydrogenase, *adhE2* bifunctional aldehyde/alcohol dehydrogenase, *alsD* alpha-acetolactate decarboxylase, *alsS* acetolactate synthase, *bcd* butyryl-CoA dehydrogenase, *bdh* butanol dehydrogenase, *buk* butyrate kinase, *crt* crotonase, *ctfAB* CoA-transferase, *etf* electron transfer

flavoprotein, *hbd* 3-hydroxybutyryl-CoA dehydrogenase, *hyd* hydrogenase, *fnor* ferredoxin-NAD(P)\* oxidoreductases, *pdc* pyruvate decarboxylase, *pfor* pyruvate:ferredoxin oxidoreductase, *pta* phosphotransacetylase, *ptb* phosphotransbutyrylase, *thl* thiolase, *gapC* NADH-dependent glyceraldehyde-3-phosphate dehydrogenase, *gapN* nonphosphorylating NADPH-producing glyceraldehyde-3-phosphate dehydrogenase, Fd ox represents oxidized ferredoxin, whereas Fd red represents reduced ferredoxin.

<figure><img src="image_7.png" /><figcaption>**Fig. 2** | **Electron flux analysis of C. acetobutylicum in phosphate-limited chemostat cultures.** Under acidogenic (pH 6.3) (a) and solventogenic (pH 4.4) (b) conditions. All values (mmol/gDCW/h) are normalized to the flux of glucose</figcaption></figure>

consumption. Data were extracted from ref. 19. The black arrow represents the total flux of Fd<sub>red</sub> oxidation, the red arrow the hydrogenase flux, the blue arrow the Ferredoxin-NAD\* flux, and the green arrow the ferredoxin-NAD\* flux.

<table><caption>**Table 1** | **Purification of ferredoxin-NADP\+ reductases from C. acetobutylicum**</caption><thead><tr><td>Steps</td><td>Activity (units)</td><td>Protein (mg)</td><td>Specific activity (units/mg)</td><td>Recovery (%)</td><td>Purification fold</td></tr></thead><tbody><tr><td>Crude extract</td><td>2.7</td><td>150</td><td>0.02</td><td>100</td><td>1</td></tr><tr><td>Streptomycin sulfate supernatant</td><td>2.46</td><td>104</td><td>0.02</td><td>92</td><td>1.33</td></tr><tr><td>CAPTO-DEAE column</td><td>0.52</td><td>1.93</td><td>0.27</td><td>20</td><td>15</td></tr><tr><td>Resource Q column</td><td>0.066</td><td>0.22</td><td>0.3</td><td>2.44</td><td>12.5</td></tr><tr><td>Superose 12 column</td><td>0.02</td><td>0.012</td><td>1.73</td><td>0.74</td><td>96</td></tr></tbody></table>

their essential physiological role in butanol synthesis in *C. acetobutylicum* and demonstrate that butanol production is limited by the electron flux between reduced ferredoxin and both NADP\+ and NAD\+. By just overexpressing the ferredoxin-NADP\+ encoding gene, it is possible to produce butanol at yield that have only been obtained before by an extensive metabolic engineering approach¹⁸.

## Results

### Biochemical analysis of the glucose to *n*-butanol pathway in *C. acetobutylicum*

*C. acetobutylicum* can, in theory, convert one mole of glucose into one mole of *n*-butanol. The conversion of one mole of glucose to two moles of acetyl-CoA is associated with the production of two moles of NADH in the EMP (Emben–Meyerhof–Parnas) pathway and two moles of reduced ferredoxin during the decarboxylation of two moles of pyruvate to two acetyl-CoA using pyruvate ferredoxin oxidoreductase (PFOR) (Fig. 1). However, under solventogenic conditions, the conversion of two moles of acetyl-CoA to one mole of *n*-butanol produces one mole of reduced ferredoxin during the reduction of crotonyl-CoA to butyryl-CoA by butyryl-CoA dehydrogenase (BCD) and consumes four moles of NADH and one mole of NADPH to reduce butyraldehyde to *n*-butanol (Fig. 1) by the NADPH-dependent alcohol dehydrogenases BdhB, BdhC and BdhA¹⁹.

As the EMP pathway produces less NADH than the *n*-butanol pathway consumes, for each mole of *n*-butanol produced, two moles of reduced ferredoxin must be used to produce two moles of NADH using a ferredoxin-NAD\+ reductase. Furthermore, as *C. acetobutylicum* does not have an oxidative pentose-phosphate pathway²⁶ and as the nonphosphorylating NADPH-producing glyceraldehyde-3-P dehydrogenase encoded by *gapN* is expressed at a low level¹⁹, for each mole of *n*-butanol produced, one mole of reduced ferredoxin must also be used to produce one mole of NADPH using a ferredoxin-NADP\+ reductase.

From this analysis, it is clear that ferredoxin-NAD\+ and ferredoxin-NADP\+ reductases are key to providing electrons for the production of *n*-butanol (Fig. 2), but until now, the proteins involved remained totally unknown. We, therefore, decided to purify all the enzymes with ferredoxin-NAD\+ or ferredoxin-NADP\+ reductase activity.

### Purification and characterization of the ferredoxin-NADP\+ reductase of *C. acetobutylicum*

Proteins with ferredoxin-NADP\+ reductase activities were purified under strict anaerobic conditions from *C. acetobutylicum* ATCC 824 crude extract as described in the methods. Proteins were first captured using a Capto-DEAE matrix, and active eluted fractions were pooled and then purified using a Resource Q column. After concentration, active eluted fractions were finally loaded onto a Superose 12 column. The results of a traditional purification are presented in Table 1. Activities were measured at 340 nm by NADP\+ reduction using CA_C0303 reduced ferredoxin as the electron donor. During the purification process, the purified enzyme lost 60% of its activity after 48 h.

Active eluted fractions from gel filtration were then subjected to denaturing gel electrophoresis. As shown in Supplementary Fig. 1, the increase in ferredoxin-NADP\+ reductase activity seems to be linked to the concentration of a protein with an apparent molecular mass of 45 kDa on SDS-PAGE.

The region of the gel corresponding to this protein was eluted, digested by trypsin and analyzed by nano-LC-MS/MS as described in the methods. The encoding gene is *CA_C0764*, annotated as an NADPH-dependent glutamate synthase beta subunit. Identical treatments were applied to the 40 and 55 kDa proteins, which were identified as thiolase (*CA_C2873*) and phosphoribosylaminoimidazole carboxamide formyltransferase-IMP cyclohydrolase (*CA_C1395*). None of these proteins are oxidoreductases. The gene identification results were used to extract the quantitative transcriptomic and proteomic data performed¹⁹, and it was confirmed that the expression of *CA_C0764*
<figure><img src="image_8.png" /><figcaption>b</figcaption></figure>

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

**Fig. 3** | **Transcriptional regulation of CA_C0764.** *CA_C0764* chromosomal region of *C. acetobutylicum* (a) with putative -35 and -10 sequences of *CA_C0764* and *CA_C0765* analyzed with the BPROM tool. **b** Analysis of *CA_C0764* expression under acidogenesis and solventogenesis. Each error bar indicates the standard deviation

around the mean of three biological samples from chemostat cultures. Data were extracted from the supplementary materials data from ref. 19. Source data are provided as a Source Data file.

from a monocistronic operon was higher under solventogenic condition than under acidogenic condition (Fig. 3).

To validate the ferredoxin-NADP⁺ reductase activity of the CAC0764 protein, the *CA_C0764* gene was cloned into a replicative plasmid to be fused with a small tag (Strep-tag II) placed in the C-terminal position of CAC0764, as described in Supplementary Method 1. The recombinant protein was homologously overexpressed in MGCΔcacI502²⁷ (Supplementary Method 2) and then purified from the crude extract in a single step using affinity chromatography on a Strep-Tactin column. Recombinant proteins were eluted with deshiobiotin, and the purity of the eluted fraction was checked by SDS-PAGE with Coomassie blue staining (Supplementary Fig. 2).

As expected, denaturing gel electrophoresis showed a single band in the eluted fraction corresponding to an apparent molecular mass of 45 kDa, demonstrating that CAC0764 was pure. In vitro ferredoxin-NAD(P)⁺ reductase and NAD(P)H-ferredoxin reductase activities were evaluated in the recovered pure fraction using NADP⁺ reduction with reduced ferredoxin as an electron donor or methyl viologen reduction by NADPH. According to Fig. 4a, purified CAC0764-Strep-tag II exhibited both NADPH-ferredoxin reductase and ferredoxin-NADP⁺ reductase activities, and neither activity was observed in the presence of NADH and NAD⁺. These results confirmed that CAC0764 is an FNR enzyme that is strictly NADPH/NADP⁺ dependent. CAC0764 shares sequence homology (~40% identity) with one subunit, NfnB, of the NADH-dependent ferredoxin-NADP⁺ oxidoreductase (NfnAB) from *Clostridium kluyveri*.

NfnAB is an electron bifurcating enzyme that catalyzes the endergonic reduction of NADP⁺ with NADH coupled to the exergonic reduction of NADP⁺ with reduced ferredoxin. In order to investigate if CAC0764 could catalyze the Nfn bifurcating reaction, an NAD⁺ reduction assay with purified CAC0764 in the presence of both NADPH and Fd<sub>ox</sub> generating systems (see methods section) was performed but no activity could be detected demonstrating the CAC0764 is a true ferredoxin-NADP⁺ oxidoreductase.

Homologs of CAC0764 were found in all acetone-butanol-ethanol (ABE) producing Clostridia. For all the clostridia belonging to clade I<sup>27</sup>, *C. acetobutylicum*, *C. roseum/C. aurantibutyricum/C. felsineum*, and *a* *C. pasteurianum*, the CAC0764 homologs were expressed as a monocistronic operon while for the Clostridia of clade II, *C. saccharobutylicum*, *C. beijerinckii*, *C. diolis* DSM 15410, *C. pasteurianum* NRRL B-598, and *C. saccharoperbutylacetonicum*<sup>27</sup> they were expressed as a bicistronic operon with a gene homologous to *nfnA* from *C. kluyveri*.

## Purification and characterization of the ferredoxin-NAD⁺ reductase of *C. acetobutylicum*

Proteins with ferredoxin-NAD⁺ reductase activities were purified under strict anaerobic conditions from *C. acetobutylicum* ATCC 824 crude extract as described in the methods. Proteins were first captured using a Capto-DEAE matrix and eluted with a linear gradient of NaCl from 0.1 to 0.25 M, yielding two peaks (one minor and one major) of ferredoxin-NAD⁺ reductase activities. The peak activity eluted fractions were
<figure><img src="image_10.png" /><figcaption>**Fig. 4 | Enzyme activities of the purified ferredoxin-NAD(P)⁺ reductases.**
a CAC0764-Strep-tag and b Bcd-EtfB-Strep-tag-EtfA complex overexpressed and purified from *C. acetobutylicum* cells. Enzymatic activities were determined using oxidized methyl viologen (red) or reduced ferredoxin (blue) according to the Methods section. Source data are provided as a Source Data file.</figcaption></figure>

pooled separately and then concentrated before being finally loaded onto a Superose 12 or a Resource Q column. The results are presented in Table 2.

Active eluted fractions from gel filtration and Resource O were then loaded onto denaturing gel electrophoresis. As shown in Supplementary Fig. 3, ferredoxin-NAD⁺ reductase activity from gel filtration was associated with the presence of 3 proteins of 41, 37, and 34 kDa, and ferredoxin-NAD⁺ reductase activity from Resource Q was linked to the presence of two proteins of 167 and 53 kDa. All proteins were eluted, digested by trypsin and analyzed by nano-LC-MS/MS as described in Methods. The results indicated that the three proteins eluted from gel filtration were the three subunits of butyryl-CoA dehydrogenase encoded by *bcd*, *etfB* and *etfA*, and the two proteins eluted from Resource Q were the two subunits of NADH-dependent glutamate synthase encoded by *gltA* and *gltB*. The BCD enzyme complex was previously shown to have NADH-ferredoxin reductase activity in the presence of crotonyl-CoA¹⁹. This study shows clearly that BCD also has ferredoxin-NAD⁺ reductase activity in the absence of crotonyl-CoA or butyryl-CoA.

To determine if ferredoxin-NAD⁺ reductase activity could be obtained in the absence of Bcd, the *etfB*-etfA genes were cloned on a replicative plasmid with a sequence encoding Strep-tag II placed in the 3' position of *etfB* with and without *bcd* as the first gene of the synthetic

operon¹⁹. The recombinant proteins were produced in *E. coli* from the two plasmids and then purified on a Strep-Tactin column. When all three genes were expressed, an active complex could be purified, while when only two genes were expressed, only *EtfB* could be purified, indicating that *EtfB* and *EtfA* cannot form a complex in the absence of Bcd.

The ferredoxin-NAD⁺ reductase, ferredoxin-NADP⁺ reductase, NADH-ferredoxin, and NADPH-ferredoxin reductase activities of the Bcd-EtfB-EtfA complex were evaluated in the purified active fraction recovered from *C. acetobutylicum*. According to Fig. 4b, the purified complex exhibited ferredoxin-NAD⁺ reductase activity, but no ferredoxin-NADP⁺ reductase activity, no NADH-ferredoxin, and no NADPH-ferredoxin reductase activity were detected. These results confirmed that in addition to its butyryl-CoA dehydrogenase activity, the Bcd-EtfB-EtfA complex can exhibit ferredoxin-NAD⁺ reductase activity in the absence of any CoA derivative.

Finally, the ferredoxin-NAD⁺ reductase, ferredoxin-NADP⁺ reductase, NADH-ferredoxin reductase, and NADPH-ferredoxin reductase activities of the GltAB complex were evaluated in the purified active fraction recovered from *C. acetobutylicum*. Both ferredoxin-NAD⁺ reductase and NADH-ferredoxin reductase activities were detected, while very low ferredoxin-NADP⁺ reductase or NADPH-ferredoxin reductase activity was measured, indicating that GltAB is mainly NADH/NAD⁺ dependent (Supplementary Fig. 4).

## Construction of knockout mutants of the genes encoding enzymes with ferredoxin-NAD⁺ reductase and ferredoxin-NADP⁺ reductase activities

To investigate the role of ferredoxin-NAD(P)⁺ reductases in the production of butanol in vivo, group II intron-based ClosTron technology²⁸ was used to inactivate the *CA_C0764*, *gltB* and *etfB* genes in the *MGCΔcac1502* a useful mutant for functional genomics analysis²⁹. This technology uses the insertion of a group II intron into a genomic target site coupled to a retrotransposition-activated marker (erythromycin resistance), allowing stable gene inactivation. The retargeted introns were first directed to be inserted at position 407/408 on the sense strand of *CA_C0764* and at position 181/182 on the sense strand of *gltB* (Supplementary Method 1). After mutagenesis, the insertion mutants were checked by combining PCR screening, sequencing and Southern hybridization (Fig. 5) (Supplementary Method 3). Both *MGCΔcac1502-CA_C0764-408s::CT* and *MGCΔcac1502-gltb181s::CT* were successfully constructed following the procedure described in Supplementary Method 2.

A similar approach using the ClosTron method was used to inactivate the *etfB* gene. Despite repeated attempts and the use of at least two different retargeted ClosTron plasmids, insertions into the *etfB* gene could not be obtained. The low number of erythromycin-resistant clones that did arise had apparently inserted elsewhere in the genome, suggesting that the ferredoxin-NAD⁺ reductase activity of BCD is essential for *C. acetobutylicum*.

<table><caption>**Table 2 | purification of ferredoxin-NAD⁺ reductases from C. acetobutylicum**</caption><thead><tr><th>Steps</th><th>Eluted peak</th><th>Activity (units)</th><th>Protein (mg)</th><th>Specific activity (units/mg)</th><th>Recovery (%)</th><th>Purification fold</th></tr></thead><tbody><tr><td>Crude extract</td><td></td><td>13.4</td><td>203</td><td>0.07</td><td>100</td><td>1</td></tr><tr><td>Streptomycin sulfate supernatant</td><td></td><td>9.02</td><td>102</td><td>0.09</td><td>67</td><td>1.28</td></tr><tr><td rowspan="2">Capto-DEAE column</td><td>Major peak 1</td><td>7.3</td><td>5.44</td><td>1.34</td><td>54</td><td>19.14</td></tr><tr><td>Minor peak 2</td><td>1.27</td><td>4.06</td><td>0.31</td><td>9.5</td><td>4.49</td></tr><tr><td>Peak 1 on Superose 12 column</td><td></td><td>3.43</td><td>1.5</td><td>2.29</td><td>25</td><td>32.7</td></tr><tr><td>Peak 2 on Resource Q column</td><td></td><td>0.56</td><td>1.5</td><td>0.40</td><td>4.2</td><td>5.71</td></tr></tbody></table>
<figure><img src="image_11.png" /></figure>

**Fig. 5** | **Construction of the MGCΔcac1502-CA_C0764-408s::CT and MGCΔcac1502-gltb18Is::CT** mutants. a Schematic representation of the *CA_C0764* gene with a group II intron inserted at position 408 on the sense strand of *CA_C0764*. b PCR screening for the identification of putative *CA_C0764-408s::CT* mutants using gene-specific primers flanking the intron insertion site (lanes 2–5) (expected size 2060 bp) and PCR control with wild-type DNA (lane 7) (expected size 270 bp), lanes 1 and 8 DNA ladder. c Southern hybridization to demonstrate the presence of a single intron insertion in the selected *MGCΔcac1502cac0764-408s::CT* mutant. The intron probe was DIG-labeled and hybridized to HindIII-HF digested

genomic DNA of the *MGCΔcac1502cac0764-408s::CT* mutant (lane 2) with an expected size of 1970 bp. The HindIII-HF digested genomic DNA of the *C. acetobutylicum* *MGCΔcac1502* strain (lane 1) was also tested as a negative control. Lane 3 is a 1 kb DNA ladder. d PCR screening for the identification of one putative *gltb18Is::CT* mutant (among three colonies that were tested) using gene-specific primers flanking the intron insertion site (lanes 1–3) (expected size 2200 bp) and PCR control with wild-type DNA (lane 4) (expected size 500 bp), lane 1 DNA ladder. Source data are provided as a Source Data file.

## Role of ferredoxin-NADP⁺ reductase in the central metabolism of C. acetobutylicum

To better understand the role of ferredoxin-NADP⁺ reductase in the central metabolism of *C. acetobutylicum*, the growth and product formation of a strain with an inactivated *cac0764* gene (*MGCΔcac1502CA_C0764-408s::CT*) and a strain overexpressing *cac0764* (*MGCΔcac1502 (pCLFCA_C0764*)*) (Supplementary Methods 1 and 2) were compared to the control *MGCΔcac1502* strain and the *MGCΔcac1502 (pCons2-1)* strain containing an empty control plasmid, respectively. The final product yields of all strains after 10 days of culture are shown in Fig. 6a, b.

Figure 6b shows that the inactivation of the gene encoding the ferredoxin-NADP⁺ reductase caused a marked decrease in the butanol yield (40% of the theoretical yield) in comparison to the parental *MGCΔcac1502* strain (54% of the theoretical yield). Conversely, the acetone + acetoin yield was almost doubled (34% of the theoretical yield in the parental *MGCΔcac1502* strain versus 50% of the theoretical yield when *CA_C0764* was inactivated). On the other hand, the over-expression of the ferredoxin-NADP⁺ reductase-encoding gene (Fig. 6a) favored butanol production and significantly increased the butanol yield (from 50% to 72% of the theoretical yield) at the expense of acetone + acetoin production, whose yield was strongly reduced from 33% to 11% of the theoretical yield. Finally, the strain with the ferredoxin-NADP⁺ reductase-encoding gene inactivated could be complemented by the introduction of the *pCLF0764* plasmid over-expressing *CA_C0764*, which restored high butanol production (68% of the theoretical yield) (Supplementary Fig. 5). All these results were confirmed by ferredoxin-NADP⁺ reductase activity measurements (Fig. 7), which showed negligible ferredoxin-NADP⁺ reductase activity (<0.001 U/mg) when *CA_C0764* was inactivated and high ferredoxin-NADP⁺ reductase specific activity when *CA_C0764* was overexpressed from the *pCLF0764* plasmid (1.54 U/mg +/− 0.07, versus 0.0326 U/mg +/− 0.007 for the control).

These results clearly demonstrated the involvement of the ferredoxin-NADP⁺ reductase enzyme in the butanol production pathway of *C. acetobutylicum* under solventogenic conditions. Moreover, it was also demonstrated that the yield of *n*-butanol was limited by the level of ferredoxin-NADP⁺ reductase activity.

## Role of each of the two enzymes with ferredoxin-NAD⁺ reductase activity in the central metabolism of C. acetobutylicum

To better understand the role of ferredoxin-NAD⁺ reductases in the central metabolism of *C. acetobutylicum*, a strain with an inactivated *gltB* gene (*MGCΔcac1502-gltB18Is::CT*) and a strain overexpressing the genes encoding the Bcd-EtfB-EtfA complex (*MGCΔcac1502 (pCLF bcd-etb-etfa*)*) (Supplementary Methods 1 and 2) were grown anaerobically in liquid flasks containing SM with 60 g/l glucose under the same conditions, and the growth and product formation were measured. Comparative phenotypic analysis was performed by measuring both glucose consumption and the concentration of fermentation products (Fig. 6). Inactivation of *gltB* had no effect on the product profile suggesting that the GltAB enzyme complex plays a minor role in the production of NADH from reduced ferredoxin. In contrast, over-expression of the *bcd-etb-etfa* operon increased the butanol yield (from 50 to 56% of the theoretical yield) at the expense of acetone and acetoin production, which was reduced from 33 to 28% of the theoretical yield. This plus the fact that a knockout *etfB* mutant is not viable, strongly suggests that the BCD complex is responsible for the ferredoxin-NAD⁺ reductase activity in *C. acetobutylicum*.

## NADPH production in the MGCΔ1502-CA_C0764-408s::CT mutant

In this study, we showed that under solventogenic conditions, *CAC0764* is the sole enzyme responsible for ferredoxin-NADP⁺ reductase activity in *C. acetobutylicum*. As it was previously demonstrated that the oxidative pentose-phosphate pathway is missing in *C. acetobutylicum*<sup>21,26</sup>, we addressed the following question: how can the *MGCΔcac1502-CA_C0764-408s::CT* mutant generate the NADPH needed for residual butanol production and anabolic reactions? Another NADPH-producing enzyme already identified in *C. acetobutylicum* is glyceraldehyde-3-phosphate dehydrogenase (GapN), encoded by *CA_C3637*, a nonphosphorylating enzyme that catalyzes the oxidation of glyceraldehyde-3-phosphate to 3-phosphoglycerate (Fig. 1). In solventogenic chemostat cultures, this enzyme was shown to catalyze less than 5% of the total flux of the EMP pathway<sup>19</sup>. To determine whether *gapN* could be upregulated in the *MGCΔcac1502-CA_C0764-408s::CT* mutant in comparison to the *MGCΔ1502* control strain, quantitative
<figure><img src="image_12.png" /><figcaption>**Fig. 6** | Comparative final product yields in mol/mol of glucose consumed for all *C. acetobutylicum* strains. a *MGCΔcac1502* (*pCons2-1*) (blue), *MGCΔcac1502* (*pCLFCA_C0764*) (red), *MGCΔcac1502* (*pCLFbcd-ctfb-etfa*) (brown). b *MGCΔcac1502* (blue), *MGCΔcac1502-gltB181s::CT* (red), *MGCΔcac1502-CA_C0764-408s::CT* (brown). Each error bar indicates the standard deviation around the mean of four independent cultures in serum bottles at 37 °C in SM medium⁹. Asterisk indicates that the results were significantly different from the control according to the one-sided Student's *t* test, with *P* value ≤ 0.05, **P** value ≤ 0.01, and ***P** value ≤ 0.001. Source data are provided as a Source Data file.</figcaption></figure>

<figure><img src="image_13.png" /><figcaption>**Fig. 7** | Specific ferredoxin-NADP⁺ reductase activity in cell-free extracts. *MGCΔcac1502* (blue) was used as a control strain and compared to both the *MGCΔ1502-CA_C0764-408s::CT* (green) and *MGCΔcac1502(pCLFCA_C0764)* (red) mutants. Bar graphs are averages from two biological duplicates. For the enzyme assays, see "Methods". Source data are provided as a Source Data file.</figcaption></figure>

reverse-transcriptase PCR (RT-qPCR) analysis was carried out in both strains to determine the relative expression of the *gapN* and *gapC* (encoding the NADH-dependent GAPDH) genes using the *fabZ* gene as the normalization reference gene³⁰. Experiments were performed as described in Supplementary Method 4, and the normalized fold

<figure><img src="image_14.png" /><figcaption>**Fig. 8** | Normalized fold expression of *gapC* and *gapN* genes. *MGCΔcac1502* (blue) was used as a control strain and compared to the *MGCΔcac1502-CA_C0764-408s::CT* (red) mutant. Bar graphs are averages from two biological duplicates. The expression of *gapC* and *gapN* was measured using RT-qPCR and the *fabZ* gene (*CA_C357I*) as an internal standard according to Supplementary Method 4. For each biological sample, the average of the technical triplicates is represented (black circle). Source data are provided as a Source Data file.</figcaption></figure>

expression of both *gapN* and *gapC* in both the *MGCΔcac1502* and *MGCΔcac1502-CA_C0764-408s::CT* strains is presented in Fig. 8.

According to Fig. 8, the relative expression level of *gapN* in the *MGCΔ1502-CA_C0764-408s::CT* mutant was 4.3-fold higher than the relative expression level of *gapN* in *MGCΔcac1502*, showing that this gene was upregulated in the mutant strain.

Based on the phenotypic analysis performed under solventogenic conditions, the ferredoxin-NADP⁺ reductase activity determination in the cell crude extract, and the GapN expression results, a redox analysis was performed to determine how reduced ferredoxin is used to achieve the redox balance in each strain (Supplementary Table 3).

When CAC0764 is inactivated, reduced ferredoxin is mainly used for hydrogen production by HydA hydrogenase, a lower amount is used for NADH production (needed for ethanol, butyrate, lactate, and butanol synthesis), and 28% of the EMP flux is catalyzed by GapN to produce the NADPH needed for *n*-butanol synthesis and anabolism. In contrast, when CAC0764 is overexpressed, reduced ferredoxin is mainly used for the NADH and NADPH formation needed for ethanol, butyrate, and butanol production and anabolic reactions.

## Discussion

Ferredoxin-NAD⁺ and ferredoxin-NADP⁺ reductase activities were measured in *C. acetobutylicum* more than 40 years ago²⁰. Their activities under different physiological conditions were further studied by other groups, but the proteins were never purified and characterized²,⁶. Their key role in butanol production was also suggested long ago from a stoichiometric model of the metabolism³¹ but never demonstrated using a reverse genetic approach. More recently, using an updated genome-scale model constrained by transcriptomics and proteomics data¹⁹, the study quantified and demonstrated the need for both ferredoxin-NAD⁺ and ferredoxin-NADP⁺ reductases to produce butanol under solventogenic conditions. As all the efforts to identify these proteins through blast searches were unsuccessful, a classical purification protocol was developed to isolate both the ferredoxin-NAD⁺ and ferredoxin-NADP⁺ reductases of *C. acetobutylicum*. The applied strategy enabled the identification of only one protein catalyzing ferredoxin-NADP⁺ reductase activity (encoded by the *CA_C0764* gene) and two enzyme complexes (Bcd-EtfB-EtfA and GltAB) having a ferredoxin-NAD⁺ reductase activity.
The ferredoxin-NADP⁺ reductase from *C. acetobutylicum* does not share any amino acid identity with previously described ferredoxin-NADP⁺ reductases from other bacteria and shares 30% identity with GltB (the β chain of NADH-dependent glutamate synthase) (Supplementary Fig. 6), explaining its incorrect annotation. However, it is expressed as a monocistronic operon, which does not fit with a genetic organization as a bicistronic operon of glutamate synthase-encoding genes. CAC0764 was demonstrated to be strictly NADPH/NADP⁺ dependent, and FAD is required to retain full enzyme activity, as is generally described for the ferredoxin-NADP⁺ reductase enzyme family²². Homologs of CAC0764 were found in all ABE (acetone-butanol ethanol) producing Clostridia. For all the clostridia belonging to clade I²⁷, *C. acetobutylicum*, *C. roseum/C. aurantibutyricum/C. felsineum*, and a *C. pasteurianum* the CAC0764 homologs were expressed as a monocistronic operon while for the Clostridia of clade II, *C. saccharobutylicum*, *C. beijerinckii*, *C. diolis* DSM 15410, *C. pasteurianum* NRRL B-598, and *C. saccharoperbutylacetonicum*²⁷ they were expressed as a bicistronic operon with a gene homologous to *nfnA* from *C. kluyveri*. We can then conclude that for all the clade I solventogenic Clostridia, NADPH is produced by a ferredoxin-NADP⁺ oxidoreductase while for all the clade II Clostridia, NADPH is produced by a NADH-dependent ferredoxin-NADP⁺ oxidoreductase.

The two enzyme complexes with NADH-ferredoxin reductase activity BCD and GltAB were previously shown to be a butyryl-CoA dehydrogenase¹⁹, and an NADH-dependent glutamate synthase¹⁹, respectively. As the inactivation of *gltB* had no effect on the product profile, the GltAB enzyme complex probably plays a minor role in vivo in electrons transfer from reduced ferredoxin to NAD⁺. In contrast, the fact that (i) overexpression of the *bcd-etfB-etfA* operon increased the butanol yield and (ii) a knockout *etfB* mutant is not viable, strongly suggests that the BCD complex is responsible for the ferredoxin-NADP⁺ reductase activity in *C. acetobutylicum*. BCD is a bifurcating enzyme reducing crotonyl-CoA to butyryl-CoA with the consumption of two NADHs and the production of one reduced ferredoxin. This allows the endergonic transfer of electrons from NADH to oxidize ferredoxin and eliminate the excess NADH associated with acetate production in acidogenic conditions⁶,¹⁹. What is demonstrated here is that in the absence of crotonyl-CoA or butyryl-CoA, this enzyme can also carry exergonic electrons transfer from reduced ferredoxin to NAD⁺.

From a physiological perspective, inactivation of the ferredoxin-NADP⁺-encoding gene significantly decreased butanol production and increased acetone production, while its overexpression had the opposite effect, with a very high *n*-butanol yield from glucose of 72% of the theoretical value (Supplementary Fig. 7). This suggests that butanol production, under solventogenic conditions, is potentially limited by the flux of NADPH production and demonstrates that the carbon fluxes can also be modulated by manipulating the electron fluxes. Furthermore, in the absence of ferredoxin-NADP⁺ reductase, *C. acetobutylicum* maintains a certain flux of NADPH production by expressing a higher level of *gapN*, a gene coding for nonphosphorylating NADP⁺-dependent glyceraldehyde-3-P dehydrogenase. In this mutant, up to 28% of the EMP flux is catalyzed by *GapN*, resulting in lower ATP production.

Attempts to knock out *etfB* to abolish ferredoxin-NAD⁺ reductase activity have been unsuccessful. This suggests that *C. acetobutylicum* cannot grow in the absence of ferredoxin-NAD⁺ reductase activity. Consistently, all the electron flux analyses performed on chemostat cultures of *C. acetobutylicum* in acidogenic, solventogenic or alcohologenic conditions¹⁹ show a high flux in the reaction catalyzed by ferredoxin-NAD⁺ reductase. Furthermore, we can demonstrate using the *iCac967* genome-scale model that an *etfB* mutant that would be unable to produce butyrate and butanol would have no ferredoxin-NAD⁺ reductase activity could survive only if it produced lactate as the only fermentation product. Such a drastic redirection of

the metabolic fluxes is probably not possible by regulation of the metabolic pathway and might explain why the *etfB* mutant is not viable.

The discovery of the ferredoxin-NADP⁺ and ferredoxin-NAD⁺ reductase-encoding genes and the demonstration of their key role in butanol production present the possibility of complementary metabolic engineering strategies to create a homobutanologenic *C. acetobutylicum* strain.

## Methods

### Bacterial strains, plasmids, and oligonucleotides

All bacterial strains, plasmids, and oligonucleotides used in or derived from this study are listed in Supplementary Tables 1 and 2. Procedures for plasmid construction are detailed in Supplementary Method 1.

### Culture media and growth conditions

*E. coli* strains were grown aerobically at 37 °C in Luria-Bertani (LB) medium supplemented, when necessary, with ampicillin (100 µg/mL) and/or chloramphenicol (30 µg/mL). Agar (15 g/L) was added prior to sterilization on LB agar plates.

MGCΔcacI502, MGCΔcacI502-CA_C0764-408s::CT, and MGCΔcacI502 gltB181s::CT strains were kept in spore form at -20 °C in synthetic medium (SM). MGCΔI502-CA_C0764-408s::CT pCLFCA_C0764, MGCΔI502 pCLFCA_C0764, MGCΔI502 pCLF bcd-etfB-etfA and MGCΔI502 pCons2-1 strains were kept on glucose SM plates with thiamphenicol (10 µg/mL) and directly used to inoculate liquid flask cultures containing glucose SM with thiamphenicol (50 µg/mL). The liquid flask cultures of all *C. acetobutylicum* strains were grown anaerobically at 37 °C in 30 mL of SM⁷ with 55 g/L glucose. Liquid cultures of the recombinant strains were run at least in quadruplicate.

### Analytical procedures

The cell concentration was measured turbidimetrically by monitoring the optical density (OD) at 620 nm; an experimentally derived correlation factor of 0.3 g cellular dry weight per OD₆₂₀ nm was used for the biomass concentration calculations. Glucose, pyruvate, lactate, acetate, butyrate, acetoin, glycerol, ethanol, acetone, and butanol concentrations were measured in the culture supernatants using high-performance liquid chromatography (HPLC) analysis (Agilent 1200 series, Massy, France)¹⁷. In all figures, yields were expressed in mol of each product forms per mol of glucose consumed. For strain comparison purposes, the butanol and acetone yield were also expressed in % of the theoretical yield with a theoretical yield of butanol or acetone formation of one mol of *n*-butanol or acetone per mol of glucose consumed.

### Enzyme assays

Ferredoxin-NAD(P)⁺ reductase activity and NAD(P)H-ferredoxin reductase activity assays were performed in an anaerobic workstation under a nitrogen atmosphere. All reagent solutions were prepared in assay buffer (previously boiled and degassed with nitrogen) and kept under a nitrogen atmosphere. Specific activities were determined in a range where linearity with protein concentration was established. Each enzyme assay was done at least in duplicate. One unit of enzyme activity is defined as the amount of enzyme that catalyzes the conversion of 1 µM of substrate per min. The concentrations of components in the reaction mixtures (1 mL of total volume) are given below.

In vitro ferredoxin-NAD(P)⁺ reductase activity was assayed by measuring the reduction of NAD⁺ or NADP⁺ using electrons from reduced ferredoxin (CA_C0303) with H₂ as the reductant of ferredoxin (CA_C0303)³² in the presence of Fe-Fe hydrogenase from *Clostridium acetobutylicum* (CA_C0028)³³. Ferredoxin (CA_C0303) was produced from anaerobic culture of the recombinant strain *E. coli* Δ*iscRpthl-FdLL-C-Tag* and purified in the in the form of Strep-tag II fused ferredoxin using a prepacked Strep-Trap HP Column (Cityva, Sweden) and an
elution buffer containing 7.5 mM desthiobiotin³³. The reaction was performed anaerobically at 37 °C in 100 mM Tris-HCl buffer (pH 7) with 2 mM DTT, 25 µM FAD, 13 µM ferredoxin or 150 µM methyl viologen, 1,6 mM NAD⁺ or NADP⁺, 6 U (or more) of purified hydrogenase HydA from *C. acetobutylicum* and crude extract (or purified protein), followed by monitoring the increase in A₃₄₀ nm as an indication of the appearance of NADH or NADPH using a spectrophotometer (Hewlett Packard 8453). After a gentle stream with hydrogen in the quartz cuvette cells, assays were initiated by the addition of ferredoxin and then, after the reduction of ferredoxin (~5 min), by the addition of NAD⁺ or NADP⁺. In all reactions, nonenzymatic rates were subtracted from the observed initial reaction rates.

In vitro NAD(P)H-ferredoxin reductase activity was assayed by monitoring the increase in A₅₆₀nm as an indication of the reduction of methyl viologen using a spectrophotometer (Hewlett Packard 8453). The reaction was carried out anaerobically at 37 °C in quartz cuvette cells in 100 mM Tris-HCl buffer (pH 7,6) with 2 mM DTT, 10 µM FAD, 250 µM NADPH, or NADH, ethanol 3% vol/vol, 45 U Adh (*S. cerevisiae*), 10 mM methyl viologen, and crude extract or purified protein. Assays were initiated by the addition of methyl viologen. In all reactions, nonenzymatic rates were subtracted from the observed initial reaction rates.

In vitro NADPH-dependent ferredoxin-NAD⁺ oxidoreductase activity was assayed³⁴ by monitoring the increase in A₃₈₀ nm as an indication of the reduction of NAD⁺ using a spectrophotometer (Hewlett Packard 8453). The reaction was carried out anaerobically at 37 °C in quartz cuvette cells in 100 mM Tris-HCl buffer (pH 7,6) with 0.5 mM NADP⁺, 40 mM glucose-6-phosphate, and 2 U of glucose-6-phosphate dehydrogenase (NADPH regeneration system), 10 mM NAD⁺, 10 µM ferredoxin, and 1 U of purified HydA hydrogenase. N₂ was used for the gas phase. The reaction was started with CAC0764 addition. (ε = 1.2 mM⁻¹ cm⁻¹). The formation of one µmol of NADH per min was defined as representing one unit.

The extinction coefficients of methyl viologen at 560 nm, of NADH at 380 nm and of NADH and NADPH at 340 nm were 7.71 mM⁻¹ cm⁻¹, 1.2 mM⁻¹ cm⁻¹, 6.22 mM⁻¹ cm⁻¹, and 6.29 mM⁻¹ cm⁻¹, respectively. The total protein concentration of the cell-free extract or purified fraction was determined using the Bradford method (Bio-Rad reagent)³⁵ with bovine serum albumin as the standard.

## Purification of the ferredoxin-NAD⁺ and ferredoxin-NADP⁺ reductases in *C. acetobutylicum* under solventogenic conditions

The *C. acetobutylicum* ATCC 824 strain was kept in spore form at -20 °C in SM. The flask cultures of *C. acetobutylicum* strains were grown anaerobically in SM, inoculated with a spore stock at 10% (v/v), and heat-shocked at 80 °C for 15 min. Cells were grown at 37 °C to an OD₆₂₀ nm of ~2.0, and the pH was maintained by buffering the culture medium with calcium carbonate prior to inoculation of the bioreactor at 10% (v/v). pH-controlled batch fermentations were performed in SM. A 2 L Biostat B bioreactor (Sartorius, Aubagne, France) was used with a working volume of 1.3 L³⁶. After sterilization, the medium was sparged with O₂-free nitrogen for 30 min. During the course of the experiment, the medium was maintained under a slight nitrogen overpressure to avoid O₂ entry into the reactor. All tubing was made of butyl rubber, and the reactor gas outlet was protected with a pyrogallol arrangement. Cultures were stirred at 300 rpm, the temperature was set at 35 °C, and the pH was maintained at 4.8 with the automatic addition of NH₄OH (3 N). The cell concentration was measured turbidimetrically by monitoring the optical density (OD) at 620 nm (Biochrom libra S11), and product formation was measured in duplicate using HPLC analysis (Agilent 1200 series, Massy, France)¹⁷. When the OD₆₂₀ nm reached ~16, after the switch from the acidogenic to solventogenic phase, cells were harvested under hydrogen pressure and transferred into an anaerobic chamber. The cells were washed and concentrated 20 times in 100 mM Tris-HCl 2 mM DTT 10% glycerol (pH 7.6) buffer and frozen at -80 °C.

All purification procedures were performed under anaerobic conditions. All purification buffers were degassed in advance, and 10 µM FAD and 2 mM DTT were added to prevent nonreversible activity losses.

Frozen cells from solventogenic batch cultures of *C. acetobutylicum* ATCC 824 were thawed and broken by sonication using an ultrasonic disintegrator (Vibracell 72434, Bioblock) at 4 °C in four cycles of 30 s at 2-min intervals. Debris was removed by centrifugation at 8600×g for 10 min at 4 °C (Sigma centrifuge 2-16 K). Nucleic acids were precipitated by the addition of streptomycin sulfate (200 µg/mL) to the supernatant and removed by centrifugation as described above. The recovered extract was then diluted five times in 100 mM Tris-HCl buffer (pH 8) before loading on a 5 mL HiTrap Capto-DEAE matrix (Cytiva, ref. 28-9165-40) connected to an AKTA purifier (Cytiva, Sweden). Active fractions were screened with the ferredoxin-NAD⁺ and ferredoxin-NADP⁺ reductase assay using reduced ferredoxin as an electron donor. The column was equilibrated in 100 mM Tris-HCl buffer (pH 8), and elution was performed with a 3-step gradient of 100 mM Tris-HCl + 1 M NaCl buffer (pH 8): 1 CV 0–4%, 20 CV 4–16% (target elution) and 5 CV 16–100%; 2 mL fractions were collected. For ferredoxin-NADP⁺ reductase activity, the most active fractions from the Capto-DEAE column were pooled before being loaded on a Resource Q column equilibrated in 100 mM Tris-HCl buffer (pH 8) for a second chromatographic step. The most active eluted fractions were then collected and concentrated on a Vivaspin 15/10,000 MW (Sartorius Stedim, ref. VS1502) to reduce the sample volume to 150 µL by centrifugation at 3000×g for 15 min. For the last purification step, a 150 µL sample volume was loaded on a Superose 12, 10/300 GL column (GE Healthcare, ref. 17-5173-01) previously equilibrated in 100 mM Tris-HCl + 150 mM NaCl buffer (pH 7.6), and 400 µL fractions were collected. For ferredoxin-NAD⁺ reductase activity, the active fractions from the Capto-DEAE column from each eluted peak were pooled before being loaded onto a Resource Q equilibrated in 100 mM Tris-HCl buffer (pH 8) for a second chromatographic step or onto a Superose 12 10/300 GL column (GE Healthcare, ref. 17-5173-01) previously equilibrated in 100 mM Tris-HCl + 150 mM NaCl buffer (pH 7.6), and 400 µL fractions were collected. Finally, the total protein concentration of the cell-free extract or purified fractions was determined using the Bradford method (Bio-Rad reagent)³⁵ with bovine serum albumin as the standard.

The yields and purification factor of each step were calculated. The purity factor of the separate active fractions was also evaluated using SDS electrophoresis in 40-mL polyacrylamide gels.

## Identification of the gene coding for ferredoxin-NAD⁺ and ferredoxin-NADP⁺ reductase activities

Active eluted fractions collected after Superose 12 or Resource Q chromatography were loaded onto denaturing gel electrophoresis, and proteins were silver stained. For the ferredoxin-NADP⁺ reductase active fraction, the region of the gel corresponding to the protein at 45 kDa was used, and for the ferredoxin-NAD⁺ reductase active fractions, the regions of the gel corresponding to the proteins at (i) 41, 37, and 34 kDa and (ii) at 167 and 53 kDa were used. The gel regions were cut out using a sterile pipette tip. The gel plugs were then used for the identification of proteins by mass spectrometry. Each sample was subjected to trypsin digestion and analyzed by *nano-LC-MS/MS* on a CapLC-Q-TOF2 (Waters) and by MALDI on a MALDI MX (Waters). The candidate proteins were identified with ProteinLynx Global Server (Waters) and Mascot (Matrix Science) software using the Protein Data Bank entry for *C. acetobutylicum*. In both analyses, for ferredoxin-NADP⁺ reductase activity, only one protein was identified with a significant score (77% sequence coverage). For ferredoxin-NAD⁺ reductase activity,
four proteins were identified with significant scores: (a) butyryl-CoA dehydrogenase, 58.6% sequence coverage; (b) EtfB, 78.4% sequence coverage; (c) GltA, 78.4% sequence coverage, and (d) GltB, 50.7% sequence coverage.

## Purification of CAC0764 and Bcd-EtfB-EtfA fused with a Strep-tag

CAC0764 protein and Bcd-EtfB-CST-EtfA complex protein were produced and purified under strict anaerobic conditions in the form of Strep-tag II fused proteins using a prepacked Strep-Trap HP Column (Cityva, Sweden). The desthiobiotin concentration in the elution buffer was 7.5 mM elution buffer, and 25 µM FAD was added to all buffers³³.

## Reporting summary

Further information on research design is available in the Nature Research Reporting Summary linked to this article.

## Data availability

The authors declare that all data supporting the findings of this study are available within the paper and its Supplementary Information file. Source data are provided with this paper.

## References

1. Jones, D. T. Applied acetone-butanol fermentation. in *Clostridia: Biotechnology and Medical Applications* (eds Bahl, H. & Dürre, P.) 125-168 (Wiley-VCH Verlag GmbH, 2001).

2. Jones, D. T. & Woods, D. R. Acetone-butanol fermentation revisited. *Microbiol. Rev.* **50**, 484-524 (1986).

3. Jones, D. T. et al. Solvent production and morphological changes in *Clostridium acetobutylicum*. *Appl. Environ. Microbiol.* **43**, 1434-1439 (1982).

4. Dürre, P. Biobutanol: an attractive biofuel. *Biotechnol. J.* **2**, 1525-1534 (2007).

5. Ni, Y. & Sun, Z. Recent progress on industrial fermentative production of acetone-butanol-ethanol by *Clostridium acetobutylicum* in China. *Appl. Microbiol. Biotechnol.* **83**, 415-423 (2009).

6. Vasconcelos, I., Girbal, L. & Soucaille, P. Regulation of carbon and electron flow in *Clostridium acetobutylicum* grown in chemostat culture at neutral pH on mixtures of glucose and glycerol. *J. Bacteriol.* **176**, 1443-1450 (1994).

7. Girbal, L. & Soucaille, P. Regulation of *Clostridium acetobutylicum* metabolism as revealed by mixed-substrate steady-state continuous cultures: role of NADH/NAD ratio and ATP pool. *J. Bacteriol.* **176**, 6433-6438 (1994).

8. Girbal, L., Vasconcelos, I., Saint-Amans, S. & Soucaille, P. How neutral red modified carbon and electron flow in *Clostridium acetobutylicum* grown in chemostat culture at neutral pH. *FEMS Microbiol. Rev.* **16**, 151-162 (1995).

9. Girbal, L., Croux, C., Vasconcelos, I. & Soucaille, P. Regulation of metabolic shifts in *Clostridium acetobutylicum* ATCC 824. *FEMS Microbiol. Rev.* **17**, 287-297 (1995).

10. Girbal, L. & Soucaille, P. Regulation of solvent production in *Clostridium acetobutylicum*. *Trends Biotechnol.* **16**, 11-16 (1998).

11. Wiesenborn, D. P., Rudolph, F. B. & Papoutsakis, E. T. Phosphotransbutyrylase from *Clostridium acetobutylicum* ATCC 824 and its role in acidogenesis. *Appl. Environ. Microbiol.* **55**, 317-322 (1989a).

12. Wiesenborn, D. P., Rudolph, F. B. & Papoutsakis, E. T. Coenzyme A transferase from *Clostridium acetobutylicum* ATCC 824 and its role in the uptake of acids. *Appl. Environ. Microbiol.* **55**, 323-329 (1989b).

13. Sauer, U. & Dürre, P. Differential induction of genes related to solvent formation during the shift from acidogenesis to solventogenesis in continuous culture of *Clostridium acetobutylicum*. *FEMS Microbiol. Lett.* **125**, 115-120 (1995).

14. Soni, B. K., Soucaille, P. & Goma, G. Continuous acetone butanol fermentation: a global approach for the improvement in solvent productivity. *Appl. Microbiol. Biotechnol.* **25**, 317-321 (1987).

15. Soni, B. K., Soucaille, P. & Goma, G. Continuous acetone-butanol fermentation: influence of vitamins on the metabolic activity of *Clostridium acetobutylicum*. *Appl. Microbiol. Biotechnol.* **27**, 1-5 (1987).

16. Jang, Y. S. et al. Enhanced butanol production obtained by reinforcing the direct butanol-forming route in *Clostridium acetobutylicum*. *mBio* **3**, e00314-12 (2012).

17. Dusséaux, S., Croux, C., Soucaille, P. & Meynial-Salles, I. Metabolic engineering of *Clostridium acetobutylicum* ATCC 824 for the high-yield production of a biofuel composed of an isopropanol/butanol/ethanol mixture. *Metab. Eng.* **18**, 1-8 (2013).

18. Nguyen, N. P. T., Raynaud, C., Meynial-Salles, I. & Soucaille, P. Reviving the Weizmann process for commercial *n*-butanol production. *Nat. Commun.* **9**, 3682 (2018).

19. Yoo, M. et al. A quantitative system-scale characterization of the metabolism of *Clostridium acetobutylicum*. *mBio* **6**, e01808-e01815 (2015).

20. Petitdemange, H., Cherrier, C., Raval, G. & Gay, R. Regulation of the NADH and NADPH ferredoxin oxidoreductase in Clostridia of the butyric group. *Biochim. Biophys. Acta* **421**, 334-347 (1976).

21. Lee, J., Yun, H., Feist, A., Palsson, B. & Lee, S. Y. Genome-scale reconstruction and in silico analysis of the *Clostridium acetobutylicum* ATCC824 metabolic network. *Appl. Microbiol. Biotechnol.* **80**, 849-862 (2008).

22. Ceccarelli, E. A., Arakaki, A. K., Cortez, N. & Carrillo, N. Functional plasticity and catalytic efficiency in plant and bacterial ferredoxin-NADP(H) reductases. *Biochim. Biophys. Acta* **1698**, 155-165 (2004).

23. Müller, V., Chowdhury, N. P. & Basen, M. Electron bifurcation: a long-hidden energy-coupling mechanism. *Annu. Rev. Microbiol.* **72**, 331-353 (2018).

24. Liang, J., Huang, H. & Wang, S. Distribution, evolution, catalytic mechanism and physiological functions of the flavin-based electron-bifurcating NADH-dependent reduced ferredoxin:NADP⁺ oxidoreductase. *Front. Microbiol.* **10**, 373 (2019).

25. Lo, J. et al. Deletion of *nfnAB* in *Thermoanaerobacterium saccharolyticum* and its effect on metabolism. *J. Bacteriol.* **197**, 2920-2929 (2015).

26. Crown, S. B. et al. Resolving the TCA cycle and pentose-phosphate pathway of *Clostridium acetobutylicum* ATCC 824: Isotopomer analysis, in vitro activities and expression analysis. *Biotechnol. J.* **6**, 300-305 (2011).

27. Poehlein, A. et al. Microbial solvent formation revisited by comparative genome analysis. *Biotechnol. Biofuels* **10**, 58 (2017).

28. Heap, J. T., Pennington, O. J., Cartman, S. T., Carter, G. P. & Minton, N. P. The ClosTron: a universal gene knock-out system for the genus *Clostridium*. *J. Microbiol. Methods* **70**, 452-464 (2007).

29. Croux, C. et al. Construction of a restriction-less, marker-less mutant useful for functional genomic and metabolic engineering of the biofuel producer *Clostridium acetobutylicum*. *Biotechnol. Biofuels* **9**, 23 (2016).

30. Jones, S. W. et al. The transcriptional program underlying the physiology of clostridial sporulation. *Genome Biol.* **9**, R114 (2008).

31. Papoutsakis, E. T. Equations and calculations for fermentations of butyric acid bacteria. *Biotechnol. Bioeng.* **26**, 174-187 (1983).

32. Demuez, M., Cournac, L., Guerrini, O., Soucaille, P. & Girbal, L. Complete activity profile of *Clostridium acetobutylicum* [FeFe]-hydrogenase and kinetic parameters for endogenous redox partners. *FEMS Microbiol. Lett.* **275**, 113-121 (2007).

33. Gauvelin, C. et al. Roles of the F-domain in [FeFe] hydrogenase. *Biochim. Biophys. Acta* **1859**, 69-77 (2018).
34. Wang, S., Huang, H., Moll, J. & Thauer, R. K. NADP+ reduction with reduced ferredoxin and NADP+ reduction with NADH are coupled via an electron-bifurcating enzyme complex in *Clostridium kluyveri*. *J. Bacteriol.* **192**, 15115–15123 (2010).

35. Bradford, M. A rapid and sensitive method for the quantification of microgram quantities of protein utilizing the principle of protein-dye binding. *Anal. Biochem.* **7**, 248–254 (1976).

36. Meynial-Salles, I. et al. Metabolic engineering of *Clostridium acetobutylicum* for the industrial production of 1, 3-propanediol from glycerol. *Metab. Eng.* **7**, 329–336 (2005).

## Acknowledgements

This work was financially supported by the Agence Nationale de la Recherche (grants ANR acetoH2 PNRB 2006 ANR Bio6 BioE-001, Bio-butafuel ANR 088 BioE-012-01, Cellutanol ANR 14 CE05-0019-02) and by the BIOCORE European Project (grant no. FP7-241566). The authors thank AJE for the English language editing of the document (certification code 3DBE-123F-E411-E361-B16P).

## Author contributions

C.F., A.R., M.H., M.Y., S.D.R., C.P., M.P., and S.D. designed and performed the experiments. L.G. participated in the conception of this study. P.S. and I.M.S. conceived this study, analyzed the data, discussed the results, and wrote the manuscript.

## Competing interests

A patent application has been filed in 2015 with the following references: Soucaille P., Meynial-Salles I., Foulquier C., and Rivière A filed a patent application (European Patent application EP15306225) in 2015 with the title "New polypeptide having ferredoxin-NADP+ reductase activity polynucleotide encoding the same and uses thereof". The remaining authors declare no competing interests.

## Additional information

**Supplementary information** The online version contains supplementary material available at https://doi.org/10.1038/s41467-022-32269-1.

**Correspondence** and requests for materials should be addressed to Philippe Soucaille.

**Peer review information** *Nature Communications* thanks Milagros Medina, Chen Yang and the other, anonymous, reviewer(s) for their contribution to the peer review of this work.

**Reprints and permission information** is available at http://www.nature.com/reprints

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

**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) 2022