# Same but different: Comparison of two system-specific molecular chaperones for the maturation of formate dehydrogenases

Nadine Schwanhold, Chantal Iobbi-Nivol, Angelika Lehmann, Silke Leimkuhler

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

## To cite this version:

Nadine Schwanhold, Chantal Iobbi-Nivol, Angelika Lehmann, Silke Leimkuhler. Same but different: Comparison of two system-specific molecular chaperones for the maturation of formate dehydrogenases. PLoS ONE, 2018, 13 (11), e0201935. doi:10.1371/journal.pone.0201935

HAL Id: hal-01927954

https://amu.hal.science/hal-01927954

Submitted on 20 Nov 2018

HAL is a multi-disciplinary open access archive for the deposit and dissemination of scientific research documents, whether they are published 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>

Distributed under a Creative Commons Attribution 4.0 International License
RESEARCH ARTICLE

# Same but different: Comparison of two system-specific molecular chaperones for the maturation of formate dehydrogenases

Nadine Schwanhold¹, Chantal Iobbi-Nivol², Angelika Lehmann¹, Silke Leimkühler¹*

1 Institute of Biochemistry and Biology, Department of Molecular Enzymology, University of Potsdam, Potsdam, Germany, 2 Aix-Marseille Université, CNRS, BIP UMR7281, Marseille, France

* sleim@uni-potsdam.de

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

## Abstract

The maturation of bacterial molybdoenzymes is a complex process leading to the insertion of the bulky bis-molybdopterin guanine dinucleotide (bis-MGD) cofactor into the apoenzyme. Most molybdoenzymes were shown to contain a specific chaperone for the insertion of the bis-MGD cofactor. Formate dehydrogenases (FDH) together with their molecular chaperone partner seem to display an exception to this specificity rule, since the chaperone FdhD has been proven to be involved in the maturation of all three FDH enzymes present in *Escherichia coli*. Multiple roles have been suggested for FdhD-like chaperones in the past, including the involvement in a sulfur transfer reaction from the L-cysteine desulfurase IscS to bis-MGD by the action of two cysteine residues present in a conserved CXXC motif of the chaperones. However, in this study we show by phylogenetic analyses that the CXXC motif is not conserved among FdhD-like chaperones. We compared in detail the FdhD-like homologues from *Rhodobacter capsulatus* and *E. coli* and show that their roles in the maturation of FDH enzymes from different subgroups can be exchanged. We reveal that bis-MGD-binding is a common characteristic of FdhD-like proteins and that the cofactor is bound with a sulfido-ligand at the molybdenum atom to the chaperone. Generally, we reveal that the cysteine residues in the motif CXXC of the chaperone are not essential for the production of active FDH enzymes.

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

OPEN ACCESS

**Citation:** Schwanhold N, Iobbi-Nivol C, Lehmann A, Leimkühler S (2018) Same but different: Comparison of two system-specific molecular chaperones for the maturation of formate dehydrogenases. PLoS ONE 13(11): e0201935.
https://doi.org/10.1371/journal.pone.0201935

**Editor:** Joel H. Weiner, University of Alberta, CANADA

**Received:** May 10, 2018

**Accepted:** August 21, 2018

**Published:** November 16, 2018

**Copyright:** © 2018 Schwanhold et al. This is an open access article distributed under the terms of the [Creative Commons Attribution License](https://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.

**Data Availability Statement:** All relevant data are within the paper.

**Funding:** This research was supported by Deutsche Forschungsgemeinschaft (LE1171/6-2 to SL). S.L. was supported by a visiting professor position at Aix-Marseille Université. The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript.

**Competing interests:** The authors have declared that no competing interests exist.

## Introduction

Molybdoenzymes comprise a large group of redox enzymes present in all kingdom of life [1]. Generally, molybdoenzymes are classified into three different families according to the ligands present at the molybdenum atom [2], namely the xanthine oxidase (XO) family, the sulfite oxidase (SO) family and dimethyl sulfoxide (DMSO) reductase family. While enzymes of the XO and SO families are present in pro- and eukaryotes, enzymes of the DMSO reductase family are exclusively present in prokaryotes. The molybdenum cofactor (Moco) in these group of enzymes is present in the form of the bis-molybdopterin guanine dinucleotide (bis-MGD) cofactor [3], generally coordinating two molybdopterin guanine dinucleotide moieties in
addition to a sulfido-/oxo-/ or selenido ligand and an amino acid ligand from the protein which can be a serine, cysteine, selenocysteine or aspartate. Further, the enzymes can contain additional cofactors such as iron-sulfur (Fe-S) clusters, cytochromes or flavin nucleotide cofactors [4]. Well-characterized examples of the DMSO reductase family include subfamily enzymes such as DMSO reductases, trimethylamine N-oxide (TMAO) reductases, formate dehydrogenases (FDH) or nitrate reductases (NR) [1]. The crystal structures of these enzymes revealed that the bis-MGD cofactor is deeply buried inside the proteins, at the end of a funnel-shaped passage giving access only to the substrate [5]. The localization of the cofactor within the enzyme suggested that chaperones are required to facilitate the insertion of the complex bis-MGD cofactor into their active sites by the involvement of a final folding step of the target enzyme after bis-MGD insertion [6–13]. For enzymes of the DMSO reductase family, these chaperones are also referred to as redox enzyme maturation proteins (REMPs) [14]. In general, these chaperones are highly specific for their target enzyme [15]. The chaperones for the DMSO reductase family enzymes were divided into different subfamilies: the NarJ-like chaperones for nitrate reductases, the DmsD-like chaperones for periplasmic DMSO reductases, the TorD-like chaperones for TMAO reductases, the YcdY-like chaperones for maturation of the YcdX proteins, and the FdhD-like chaperones for the maturation of several distinct FDHs [13, 16–18].

For the FdhD-subfamily, well-characterized representatives comprise *Escherichia coli* FdhD and *Rhodobacter capsulatus* FdsC [19–21]. FdhD displays an exception of the general proposed rule that all molybdochaperones are specific for one target molybdoenzyme [14], since the FdhD is essential for the maturation of all three *E. coli* FDHs, namely FdhF, FdoGHI, and FdnGHI [20]. The three *E. coli* FDH enzymes have in common that their catalytic subunits coordinate the bis-MGD cofactor with a selenocysteine ligand and a sulfido-ligand, and as additional cofactor a [4Fe-4S] cluster that is present in vicinity of the bis-MGD moiety [22, 23]. Additional subunits and the localization of the enzymes within the cell vary. FdnGHI and FdoGHI show high sequence similarities, both enzymes are trimers facing the periplasm that are membrane-bound via the I-subunit [24, 25]. The electrons are transferred from the substrate converted in the catalytic subunit (FdnG or FdoG) via the five Fe-S cluster-containing subunit (FdH or FdoH) to the membrane anchoring cytochrome-containing subunit (FdN or FdoI) finally to the quinone-pool [26]. FdnGHI expression is induced in the presence of nitrate under anaerobic conditions. In contrast, FdoGHI is also present at low levels under aerobiosis, during fermentative conditions and during nitrate respiration [26, 27]. The third FDH in *E. coli*, FdhF is facing the cytoplasm and forms the formate hydrogen lyase (FHL) complex together with the hydrogenase-3 (encoded by *hycBCDEFG*) [28]. The FHL complex is membrane bound and produces H₂ and CO₂ under fermentative conditions [29]. The expression of the gene is repressed by the presence of nitrate [29]. Recently, the role of the specific chaperone FdhD was investigated in detail for FdhF maturation, and it was shown that FdhD is essential for the insertion of the terminal sulfido-ligand present at the bis-MGD cofactor of FdhF in addition to FdnG and FdoG [19, 20]. FdhF was inactive when FdhD was absent during expression, but after purification the enzyme could be activated by chemical sulfuration [20]. For sulfuration of bis-MGD, FdhD specifically interacts with the L-cysteine desulfurase IscS in *E. coli* [20]. It has been suggested that IscS transfers the sulfur from L-cysteine to FdhD in form of a persulfide. The cysteines 121 and 124 located in a conserved CXXC motif of FdhD were proposed in the studies by Thomé et al. [20] and Arnoux et al. [19] to be involved in the sulfur transfer process from IscS to bis-MGD. Cysteine to alanine variants in these residues were shown to produce an inactive FdhF enzyme. FdhD was co-crystallized in complex with GDP and direct binding of bis-MGD to FdhD was suggested [19], but it has not been proven so far that the cofactor is bound in an active form.
Previously, the binding of bis-MGD to the FdhD-like chaperone FdsC from *R. capsulatus* has been investigated in detail [21]. Using *R. capsulatus* FdsC as source of bis-MGD cofactor in a direct reconstitution assay, it was shown that an intact bis-MGD cofactor was inserted into *E. coli* TMAO reductase (TorA) that reconstituted TorA activity [21]. In contrast to *E. coli*, *R. capsulatus* harbors only one FDH enzyme that is encoded by the *fdsGBACD* operon [30]. The enzyme was shown to be structurally different from the *E. coli* FDH enzymes, since a soluble (FdsGBA)₂ heterodimer is formed that is located in the cytoplasm. The purified enzyme was shown to be oxygen tolerant using NAD⁺ as terminal electron acceptor. However, the physiological role of this enzyme has not been characterized to date. Previous studies showed that the (FdsGBA)₂ enzyme expressed in the absence of FdsC contained a bis-MGD cofactor that lacked the terminal sulfido group [31, 32]. However, in contrast to *E. coli* FdhD, the cysteines 104 and 107 in FdsC (corresponding to Cys121 and Cys107 in FdhD) were not essential for the production of an active *R. capsulatus* FDH [21]. In *R. capsulatus*, three L-cysteine desulfurases are present: NifS2, NifS3 and NifS4 [33]. While the specific role of all three enzymes has not been characterized in detail so far, NifS4 was shown to be involved in the mobilization and transfer of sulfur to the Moco present in *R. capsulatus* xanthine dehydrogenase [33].

In this study, we compared the roles of FdhD and FdsC in the maturation of FDH enzymes and investigated the roles of the conserved cysteines in the conserved CXXC motif in this process. Since FdhD and FdsC share an amino acid sequence similarity of 46.3% and an identity of 32.4%, both chaperones were suggested to be functional homologues. We show that FdhD and FdsC have a common role that is not only specific to the FDH enzymes present in the respective organism. FdsC and FdhD can be interchanged on structurally very distinct FDH enzymes. Further, we demonstrate that the cysteine residues in the motif CXXC are not essential for the production of active FDH. Solely, FdhD-C121 was revealed to have a particular role in stimulating the activity of *E. coli* IscS exclusively. No similar effect was identified for *R. capsulatus* FdsC. We suggest that FdhD-Cys121 is assisting the specific interaction with *E. coli* IscS and enhances its L-cysteine desulfurase activity. By this process a higher sulfuration and transfer efficiency of bis-MGD is obtained, which might be only required for the *E. coli* system, since *E. coli* IscS has numerous interaction partners.

# Results

## A phylogenetic view of the FdhD family of molecular chaperones for FDH enzymes

Analysis of bacterial and archaeal genome sequence data in 5 phyla and 33 taxonomic families revealed a wide distribution of homologues to the *E. coli* FdhD protein (referred to as FdhD family of chaperones) (Fig 1). In bacteria like *E. coli*, where three FDH enzymes are present, only one FdhD-homologue has been identified. We aligned 41 FdhD-like sequences, which present an overall level of identity of about 18–34%. From the amino acid sequence alignment, an unrooted phylogenetic tree has been constructed (Fig 1). Overall, the phylogenetic tree can be divided into three groups. In group I only α-, β- and γ-proteobacterial FdhDs are present, including the ones from *E. coli* and *R. capsulatus*. Group II contains mainly FdhDs from Actinobacteria, but also from Cyanobacteria, Archaea and γ- or δ-Proteobacteria. Group III contains a larger group of archaeal FdhDs in addition to FdhD proteins from Bacilli and Clostridia or γ-, δ-, and ε-Proteobacteria. Overall, Firmicutes, Proteobacteria and Actinobacteria are known to contain a high to medium number of respiratory molybdoenzymes [34] which is obvious by the distribution of FdhD homologues in these organisms. Interestingly, γ-Proteobacteria are represented in all three groups. Many pathogens belong to this class of
<figure><img src="image_5.png" /><figcaption>**Fig 1. Phylogenetic tree of FdhD-like proteins.** Protein phylogeny of FdsC/FdhD homologues based on a full length multiple sequence alignment by Muscle [50]. The tree was constructed using the Maximum Likelihood method based on the Dayhoff matrix based model [51] [57]. The tree is drawn to scale, with branch lengths in the same units as those of the evolutionary distances used to infer the phylogenetic tree. The evolutionary distances were computed using the Dayhoff matrix based method and are in the units of the number of amino acid substitutions per site [52]. The scale bar indicates 0.2 substitutions per site. Numbers near branches indicate the bootstrap proportion for 100 replicas using the same method. The analysis involved 41 amino acid sequences. All positions containing gaps and missing data were eliminated. There was a total of 174 positions in the final dataset. Evolutionary analyses were conducted in MEGA7 [53]. Domains and classes of prokaryotes are marked in blue: α-Proteobacteria, β-Proteobacteria, γ-Proteobacteria, δ-Proteobacteria, ε-Proteobacteria. There are three main branches: group I, II and III. FdsC/FdhD homologues in group I and II harbor the conserved CXXC-motif, which is not present in group III (◆ - CC-motif, □ - only second cysteine, ■ - CXXXCXC-motif, ⊕ - no cysteine).</figcaption></figure>
Existing FDHs are listed as described in the box on the left and based on genomic analysis using www.ncbi.nlm.nih.gov or * img.jgi.doe.gov. If already described in literature, the corresponding active site metal (Mo or W) is displayed.

https://doi.org/10.1371/journal.pone.0201935.g001

bacteria which are known to harbor a high number of molybdoenzymes representing a higher
metabolic flexibility [34]. Thus, a wide distribution of FdhD-like chaperones is found here. To
provide more details, we also tried to include the different target FDH enzymes in the phyloge-
netic tree, to analyze whether FdhD has evolved together with its protein partner [14]. Here,
we divided the FDH enzymes into membrane-bound FdhF-type, dimeric FdhAB-type, and
cytosolic NAD⁺-dependent FDHs. Unfortunately, a lot of FDHs are not characterized to date
so that our division is mainly based on amino acid sequence homologies of characterized FDH
enzymes. The overall distribution revealed that the FdhD homologues for NAD⁺-dependent
FDHs are mainly present in group I, while group II predominantly contains FdhDs required
for monomeric FdhF-type enzymes. Group III can be divided into two subgroups, while one
subgroup contains mainly FdhD-homologues for FdhF-type enzymes, the second subgroup
predominantly contains FdhD proteins for FdhAB-type FDHs. A previous study reporting on
the evolution of the chaperones for the DMSO reductase family enzymes suggested that the
molecular chaperones evolved together with its molybdenum partner protein [14]. However, a
division into FdhDs for Mo- or W-containing FDHs was not possible, and FdhD proteins act-
ing on either Mo- or W-containing (or both) FDHs are distributed in all three groups of the
phylogenetic tree without any obvious division.

It has been reported previously that FdhD proteins are characterized by a highly conserved CXXC motif that is present on a disordered loop of the protein [19]. In *E. coli* FdhD, the cysteines 121 and 124 of this motif have been reported to be functionally important for the catalytic sulfur transferring activity from IscS for the formation of a sulfurated bis-MGD cofactor [20]. However, as obvious from the data in the phylogenetic tree, a detailed analysis of the amino acid sequence of the FdhD-like proteins showed that the CXXC motif is not highly conserved. FdhD proteins containing the CXXC motif are mainly present in group I and group II of the phylogenetic tree. Group III, in contrast, contains mainly FdhD sequences in which either only the second cysteine of the motif is conserved, which have a CC-motif instead or have no cysteines. Thus, this puzzling point encouraged us to understand the importance of the cysteines for the role of FdhD-like proteins from *E. coli* and *R. capsulatus*.

Copurification of *E. coli* FdhD and *R. capsulatus* FdsC variants with bis-MGD

For *E. coli* FdhD, it has been reported that especially the first cysteine of the two cysteines in
the conserved ¹²¹CXXC¹²⁴ motif is essential for its role in transferring the sulfur from IscS to
the bis-MGD cofactor, while Cys124 is additionally required for the production of an active
FdhF enzyme [20]. In contrast, it has been shown for *R. capsulatus* FdsC that the two con-
served cysteine residues are not essential for producing an active sulfido-containing
(FdsGBA)₂ enzyme [21]. To test whether the role of the cysteines is specific for the *E. coli* FDH
enzymes, which are structurally and functionally different from the soluble *R. capsulatus*
(FdsGBA)₂ enzyme, we constructed single and double cysteine variants of *E. coli* FdhD
(C121A, C124A, C121A/C124A) and compared their characteristics to *R. capsulatus* FdsC. For
*R. capsulatus* FdsC the corresponding cysteine variants were already available (C104A, C107A,
C104A/C107A) [21]. To directly compare the role of the two chaperones, the three cysteine
variants of both FdsC and FdhD in addition to the wild-type protein were purified from an *E*.
<figure><img src="image_6.png" /><figcaption>**Fig 2.** Amino acid sequence alignment of *E. coli* FdhD and *R. capsulatus* FdsC. The amino acid sequence alignment was generated using ClustalW. Identical amino acids and similar amino acids are marked with a black and a grey box, respectively. The conserved cysteine residues are labelled with red arrows. The blue bar over the alignment indicates the N-terminal domain (blue), the linker region (yellow) and the C-terminal domain (green) based on the structure of FdhD (PDB: 4pde) [19].</figcaption></figure>

https://doi.org/10.1371/journal.pone.0201935.g002

*coli* Δ*fdhD* strain. An amino acid sequence alignment of both proteins is shown in Fig 2, highlighting the position of the cysteines within the protein sequences.

The proteins were purified by Ni-NTA chromatography and the Coomassie stained SDS-gels showed that all eight proteins were purified with a purity of more than 90%. His₆-tagged FdhD has a molecular mass of 32.7 kDa, while His₆-tagged FdsC has a molecular mass of 28.6 kDa (Fig 3A and 3B). Size exclusion chromatography confirmed that all FdsC and FdhD variants were purified as a dimer in solution (data not shown). For FdsC it has been shown before that the protein was copurified with bound bis-MGD cofactor under our expression conditions [21]. To analyze the bis-MGD content for all FdsC and FdhD variants, the bis-MGD cofactor was extracted from the respective protein, oxidized by I₂/KI overnight at room temperature, and the obtained fluorescent derivative FormA-GMP was quantified after separation on a reversed-phase C18 HPLC column. The results in Fig 3C show that all FdsC and FdhD variants were copurified with the bis-MGD cofactor, however, the bis-MGD content of the FdsC wild type and variants was in average twice as high in comparison to the FdhD proteins. The bis-MGD binding ability thereby was not influenced by the cysteine to alanine exchanges in the variants.

## Reconstitution of TMAO reductase activity with bis-MGD provided by FdsC or FdhD

An *in vitro* assay system has been established previously by which the direct insertion of the bis-MGD cofactor bound to FdsC into apo-TorA can be monitored [21]. Very recent results showed that TorA is only active when the bis-MGD cofactor contains a terminal sulfido group [35]. Thus, the *in vitro* reconstitution of TorA activity can be used to determine the sulfuration level of the bis-MGD cofactor bound to FdhD or FdsC.

For the *in vitro* reconstitution assay, 8 µM of Moco-deficient apo-TorA were incubated with 50 µM of the FdhD or FdsC-variants and incubated for 2 hours at 37°C under anaerobic conditions without any further supplementations. The results in Fig 4 show that active TorA was obtained for both the FdhD and FdsC variants. Unexpectedly, while the TorA activity varied around 60–80 U/mg when FdsC wild-type or the C104A, C107A or C104A/C107A variants were used as bis-MGD source, almost 2-fold higher TorA activities were obtained with FdhD and the variants C121A, C124A or C121A/C124A with values around 130–150 U/mg. No clear
<figure><img src="image_7.png" /><figcaption>**Fig 3. Characterization of cysteine variants of FdsC and FdhD.** FdsC and FdhD and their cysteine variants were expressed in *E. coli* Δ*fdhD* strain and purified by Ni-NTA chromatography as described for FdsC previously [21]. A) 15 µg of purified FdsC, FdsC-C104A, FdsC-C107A and FdsC-C104/C107A were separated by SDS-PAGE. B) 15 µg of purified FdhD, FdhD-C121A, FdhD-C124A and FdhD-C121/C124A were separated by SDS-PAGE. C) 200–300 µM of protein was treated with I₂/KI-HCl to isolate Form A-GMP as described previously [47]. Form A-GMP was separated by HPLC. The corresponding peak area was normalized for protein concentration. The results represent the mean values from three independent measurements (±S.D.).</figcaption></figure>

https://doi.org/10.1371/journal.pone.0201935.g003

differences in activity were obtained between FdsC and the cysteine substitution variants or FdhD and the cysteine substitution variants.

Since the results in Fig 3C showed that overall, FdhD and the C121A, C124A or C121A/C124A variants were loaded with lower levels of bis-MGD, we compared the bis-MGD content of the reconstituted TorA proteins using FdhD or FdsC proteins. We analyzed only the bis-MGD loading of TorA after reconstitution using FdsC and FdhD wild-type proteins, since the
<figure><img src="image_8.png" /><figcaption>**Fig 4.** Reconstitution of *E. coli* apo-TMAO reductase using FdsC or FdhD as bis-MGD source. 50 µM of each purified chaperone variant was incubated with apo-TMAO reductase (8 µM) to reconstitute enzyme activity. After 2 hours TMAO reductase activity was measured under anaerobic conditions following the oxidation of reduced benzyl viologen at 600 nm in the presence of 0.1 µM TMAO. Kinetic Data are mean values from three independent measurements (±S.D.).</figcaption></figure>

https://doi.org/10.1371/journal.pone.0201935.g004

analysis requires a lot of source TorA protein for bis-MGD quantification and differences were not expected between the wild-type and the cysteine variants according to Fig 3. FdsC-reconstituted TorA contained about two-times the amount of bis-MGD in comparison to FdhD-reconstituted TorA (Fig 5A). However, when the obtained specific TMAO reductase activities from Fig 4 of FdsC and FdhD wild-type were related to their bis-MGD levels (detected as For-mmA-GMP), it became obvious that about 5-fold higher relative TorA activity in relation to bis-MGD bound to the protein was obtained with FdhD in comparison to FdsC. Since both proteins bind the bis-MGD cofactor, the differences in activity can only be explained by a different saturation level of the sulfido-ligand at the molybdenum atom. In a recent report by Kaufmann et al. [35], it was shown that the presence of a terminal sulfido-ligand at the active site of TorA contributed to a 20-fold increased TorA activity. In contrast, enzyme preparations containing an oxo-ligand instead were mainly inactive. Therefore, it has to be concluded that the bis-MGD cofactor on FdhD is more saturated with the sulfido-ligand than the cofactor present on FdsC, which of course can only be qualitatively speculated from this assay.

## Influence of FdhD and FdsC variants on the maturation of FDH enzymes from different classes

So far, the results showed that the cysteine substitutions in either FdsC or FdhD did not influence the bis-MGD binding ability or transfer of the cofactor to TorA. For *E. coli* FdhD, however, it has been reported before that Cys121 and Cys124 are essential for the production of an active FdhF enzyme [20]. In contrast, for the corresponding C104A and C107A substitution variants of *R. capsulatus* FdsC, no effect on (FdsGBA)₂ activity was reported, using a heterologous expression system in *E. coli* [21]. To differentiate whether the observed differences in the role of the cysteine residue is based on differences in the chaperone or the type of FDH enzyme used, we analyzed the effect of *E. coli* FdhD on *R. capsulatus* (FdsGBA)₂ activity and of *R. capsulatus* FdsC on *E. coli* (FdoGHI)₃ activity.
<figure><img src="image_9.png" /><figcaption>**Fig 5. Characterization of *E. coli* apo-TMAO reductase reconstituted with FdsC or FdhD.** For the reconstitution of enzyme activity, 50 µM of each FdhD, FdsC or their variants was incubated with apo-TMAO reductase (8 µM) in a total volume of 8 ml. After 7 hours, the samples were concentrated to 500 µl before size exclusion chromatography. The results represent the mean values from three independent measurements (±S.D.). A) Relative Form A-GMP content (MGD in LU\*s per mg TMAO reductase) of 3 µM TMAO reductase was analyzed by HPLC. B) Activity of the fraction containing TMAO reductase was measured under anaerobic conditions following the oxidation of reduced benzyl viologen at 600 nm in the presence of 0.1 µM TMAO. TMAO reductase activity was normalized for the MGD content.</figcaption></figure>

https://doi.org/10.1371/journal.pone.0201935.g005

**Fig 6** shows the results on *R. capsulatus* (*FdsGBA*)₂ activity after coexpression with *E. coli* FdhD wild-type or the cysteine variants in an *E. coli* Δ*fdhD* deletion strain. The results show that active *R. capsulatus* (*FdsGBA*)₂ enzymes were obtained, showing that *E. coli* FdhD can replace FdsC in its function. Surprisingly, the FDH activity was even 2.6-fold higher with FdhD wild-type in comparison to FdsC wild-type. When the cysteine variants C121A, C124A or the double variant C121A/C124A of FdhD were coexpressed instead, a 1.4 times lower activity was obtained in comparison to FdhD wild-type, showing that the cysteine substitutions indeed have an effect on the role of FdhD for the maturation of (*FdsGBA*)₂. However, the results show that FdhD-C121 or FdhD-C124 are not essential for the maturation of the *R. capsulatus* enzyme, since FDH activity was obtained. Rather, the presence of both cysteines was required to obtain the highest FDH activity.

For comparison, the roles of the cysteines in FdhD were additionally analyzed for *E. coli* (FdoGHI)₃ activity. Here, we made use of an in-gel activity assay that directly detects the activity of (FdoGHI)₃ after separation in native gels. To suppress the activity of (FdnGHI)₃, cells were grown with 0.2% sodium formate in the absence of nitrate. The gel was incubated at 37°C under aerobic conditions with a staining solution containing sodium formate, NBT and PMS.
<figure><img src="image_10.png" /><figcaption>**Fig 6. Influence of FdsC and FdhD on FDH activity.** A) *R. capsulatus* (FdsGBA)₂ was expressed in *A*fdhD strain in the presence of FdhD WT (pNB14), FdhD-C121A (pNB15), FdhD-C124A (pNB16), FdhD-C121A/C124A (pNB17) or FdsC WT (pTHfds03) and purified by Ni-NTA affinity chromatography as described previously [21]. *R. capsulatus* (FdsGBA)₂ activity was detected photometrically by the increase in NADH recorded at 340 nm and 1 U/mg is defined as the reduction of 1 µmol NAD⁺/min/mg of enzyme at room temperature. Kinetic Data are mean values from three independent measurements (±S.D.). The enzyme expressed in the absence of FdsC or FdhD was not active (data not shown), as reported previously [21]. B) 55 ml cultures of each strain were grown anaerobically at 37°C for 16 hours in LB media in the presence of 20 µM IPTG, 10 µM sodium molybdate and antibiotic as needed. Equivalent amounts (40–60 µg) of Triton-X100 treated crude extracts were applied to each lane and separated by non-denaturing PAGE. The gels were stained with 1 mM NBT, 0.5 mM PMS and 50 mM formate in 50 mM potassium phosphate, pH 6.8. Lanes from left to right: BW25113, ΔfdoG, ΔfdhD, +pFdhD: *A*fdhD deficient strain complemented with plasmid pNB14, +pFdhD-C121A: *A*fdhD deficient strain complemented with plasmid pNB15, +pFdhD-C124A: *A*fdhD deficient strain complemented with plasmid pNB16, +pFdhD-C121A/C124A: *A*fdhD deficient strain complemented with plasmid pNB15, +pFdsC: *A*fdhD deficient strain complemented with plasmid pTHfds14, +pFdsC-C104A: *A*fdhD deficient strain complemented with plasmid pNBfds04, +pFdsC-C107A: *A*fdhD deficient strain complemented with plasmid pNBfds05, +pFdsC-C104A/C107A: *A*fdhD deficient strain complemented with plasmid pNBfds06.</figcaption></figure>

https://doi.org/10.1371/journal.pone.0201935.g006

The (FdoGHI)₃ activity obtained after expression of FdhD wild-type or the C121A, C124A, C121A/C124A variants in an *E. coli* ΔfdhD strain showed that the highest activity was obtained
with the FdhD wild-type protein (Fig 6A). In contradiction to previous reports, we obtained an (FdoGHI)₃ activity when the FdhD variants C121A, C124A or the double variant C121A/C124A were used [19, 20]. While the in-gel activities of the FdhD-C124A and the FdhD-C121A/C124A variants were comparable, the FdhD-C121A variant showed a slightly reduced activity (Fig 5A). However, since the FdhD-C121A/C124A double variant was active, the results show that Cys121 is not essential for FdoGHI activity.

In contrast to FdhD, the cysteine variants of *R. capsulatus* FdsC, C104A, C107A, and the double variant C104A/C107A showed no difference on FdoGHI activity in comparison to FdsC wild-type (Fig 6B). Additionally, the in-gel activities were overall comparable to the FdhD cysteine variants. These results reveal that differences in the maturation of FdoGHI exist between FdsC and FdhD wild-type proteins. However, when the cysteines in FdhD were substituted, the maturation efficiency of the FdhD variants was comparable to FdsC. This implies that the cysteines in FdhD have a role that is specific to the FdhD protein.

## Influence of FdhD and FdsC on the L-cysteine desulfurase activity of *E. coli* IscS

The results shown above indicated that the enhanced (FdoGHI)₃ activity by FdhD wild-type in comparison to its cysteine variants might be based on an enhancing effect of FdhD on the L-cysteine desulfurase activity of IscS. IscS was shown before to provide the sulfur for the formation of the sulfido ligand on the bis-MGD cofactor. It has been reported previously that FdhD stimulates the activity of IscS, while this stimulating effect was lost in the FdhD-C121A single and FdhD-C121A/C124A double variants [19, 20]. We were able to reproduce these results (Fig 7A), however, when we tested the effect of *R. capsulatus* FdsC on IscS activity, we did not observe a stimulating effect (Fig 7A). Further, we tested the effect of both FdhD and FdsC on the L-cysteine desulfurase proteins present in *R. capsulatus*. *R. capsulatus* contains the three L-cysteine desulfurases NifS2, NifS3 and NifS4 as house-keeping L-cysteine desulfurases, while a fourth one, NifS, was shown to be specific for nitrogenase [33]. The results in Fig 7B show that neither FdhD nor FdsC could stimulate the L-cysteine desulfurase activity of any of the three *R. capsulatus* L-cysteine desulfurases NifS2, NifS3 or NifS4. A truncated version of Nfs2 was used, since the full-length protein was shown to be unstable [33]. As positive control, we used the *R. capsulatus* SufE protein, the physiological partner of NifS4 that was able to stimulate the activity of NifS4 but not of the other L-cysteine desulfurases. The results imply that FdhD-C121 is specifically required for stimulating the activity of *E. coli* IscS.

## Interaction of FdhD and FdsC with L-cysteine desulfurases

It has been reported previously that FdhD forms a complex with IscS that can be copurified [19, 20]. To analyze whether the FdhD-C121A/C124A variant effects the interaction with IscS, we performed analytical size exclusion chromatography. For comparison, we also analyzed the complex formation of IscS with FdsC and the FdsC-C104A/C107A variant.

The chaperones FdhD or FdsC were incubated with IscS and complex formation was analyzed after separation on a Superdex 200 column. When a complex between IscS and one of the chaperones is formed, the complex elutes more rapidly from the size exclusion column due to its increase in molecular mass. A 1:1 complex between the IscS dimer and FdhD is expected, as reported previously [19, 20]. The results in Fig 8A show that FdhD and IscS readily formed a complex that co-eluted at 11 ml after separation on a Superdex 200 column. The peak fractions eluting 10-12 ml were separated by SDS-PAGE and showed the presence of both proteins, confirming the formation of a protein-complex. The complex formation was not
<figure><img src="image_11.png" /><figcaption>**Fig 7.** Influence of FdsC and FdhD on L-cysteine desulfurase activity. L-cysteine desulfurase activity was measured by determination of total sulfide produced [48]. A) 1 µM IscS from *E. coli* was mixed with either 2 µM FdhD/FdhD variants or FdsC/FdsC variants and incubated for 10 min in the presence of 1 mM L-cysteine at 30°C. The activity of IscS alone was set to 1. B) L-cysteine desulfurase IscS (1 µM) from *E. coli* or NifS4, NifS3 and NifS2-Δ1-188 from *R. capsulatus* (2 µM) were mixed in a 1:2 ratio with FdsC, FdhD or SufE, respectively and incubated for 10 min in the presence of 1 mM L-cysteine at 30°C. The fold induction of IscS, NifS4, NifS3 or NifS2-Δ1-188 activity by incubation with FdsC, FdhD or SufE is relative to the activity of the L-cysteine desulfurase alone, respectively. Kinetic Data are mean values from three independent measurements (±S.D.).</figcaption></figure>

https://doi.org/10.1371/journal.pone.0201935.g007
<figure><img src="image_12.png" /></figure>

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

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

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

**Fig 8.** Interaction of FdhD with IscS. Complex formation of (A) 40 µM FdhD + 20 µM IscS, (B) 40 µM FdhD-C121A/C124A + 20 µM IscS, (C) 40 µM FdsC + 20 µM IscS and (D) 40 µM FdsC-C104A/C107A + 20 µM IscS were analyzed by size exclusion chromatography on a Superdex Increase 200 column (GE Healthcare) equilibrated in 100 mM potassium phosphate buffer, pH 8.0, 100 mM NaCl, 10 mM β-mercaptoethanol. The elution of proteins was followed at 280 nm. Indicated fractions in a range of 10.5–15.5 ml (0.5 ml fractions) were analyzed for their protein content by 15%SDS–PAGE.

https://doi.org/10.1371/journal.pone.0201935.g008

impaired when using the FdhD-C121A/C124A variant showing that the two cysteines are not essential for complex formation, since a comparable elution volume of 11 ml was obtained (Fig 8B). In contrast, no complex formation was determined between IscS and FdsC or the FdsC-C104A/C107A variant (Fig 8C and 8D), as obvious from the elution profile and the corresponding fractions separated by SDS-PAGE, since no peak with a different elution time was obtained after incubation of the proteins. The results therefore reveal that FdhD forms complex with IscS that can be copurified, in contrast to FdsC, where no stable complex with IscS was formed.
Discussion

FdhD-like chaperones are found in a vast number of prokaryotes. They are often encoded in the same operon together with the structural genes for FDH enzymes or in vicinity to genes involved in Moco biosynthesis [22]. So far, the best characterized chaperones are FdhD from *E. coli* and FdsC from *R. capsulatus* [19–21].

Here, we addressed the questions whether FdsC from *R. capsulatus* and FdhD from *E. coli* are functional homologues and have exchangeable roles. Both chaperones are essential for their corresponding FDH enzyme(s), but so far, their exchangeable roles have not been investigated in detail. The crystal structure of dimeric *E. coli* FdhD in complex with GDP showed two symmetrical binding sites for nucleotides coming from the arm in the N-terminal domain of one monomer and from the C-terminal domain of the adjacent monomer (Fig 2) [19]. Based on the GDP binding sites the bis-MGD binding site was modeled on FdhD [19].

In this study, we compared the copurification of the chaperones with bound bis-MGD directly. We concluded that the bis-MGD cofactor bound to FdhD is present in a form with a higher saturation of the terminal sulfido-ligand. This was confirmed by a higher activity of TorA obtained after reconstitution with the cofactor bound to FdhD in comparison to FdsC. Recently it was shown that TorA contains a terminal sulfido ligand at the molybdenum atom that is essential for obtaining a fully active enzyme [35]. Thus, our conclusions seem feasible, since the higher activity was not based on a better transfer of the bis-MGD cofactor from FdhD to TorA.

We further investigated the roles of the cysteines in the conserved CXXC motif present in groups I and II of FdhD-like chaperones. This motif is not present or is modified in enzymes of group III (Fig 1). We showed that the cysteines did not impact the copurification of FdhD or FdsC with bis-MGD (Fig 3). However, both cysteines of FdhD had an impact on the activity of either *E. coli* (FdoGHI)₃ or *R. capsulatus* (FdsGBA)₂ (Fig 6). For *R. capsulatus* FdsC, in contrast, no difference in FDH activity was observed when the cysteine variants were used for the maturation of either *E. coli* (FdoGHI)₃ or *R. capsulatus* (FdsGBA)₂ (Fig 6B and reference [21]). Since no differences in bis-MGD transfer to a recipient protein was obtained between the FdhD-like chaperone and its cysteine variants (Fig 4), we investigated the interaction and sulfur transferase activity with IscS, the interaction partner protein from which the sulfur atom for bis-MGD sulfuration is obtained.

While FdhD wild-type was able to enhance the L-cysteine desulfurase activity of IscS, this effect was absent in the cysteine substitution variants (Fig 7) [20]. Further, this stimulating effect is an exclusive characteristic for FdhD with its partner protein IscS, since FdhD was unable to stimulate the L-cysteine desulfurase activity of the three L-cysteine desulfurases present in *R. capsulatus*. Moreover, the *R. capsulatus* FdsC chaperone does not provide this stimulating effect neither with the L-cysteine desulfurases from its own organism nor with the IscS protein from *E. coli*. Studies on the interactions between FdhD and FdsC and the different L-cysteine desulfurases revealed that FdhD and IscS form tight complex that can be copurified (Fig 8). Cys121 in *E. coli* FdhD seems to have an additional exceptional role in enhancing the L-cysteine desulfurase activity exclusively of IscS. This role might be especially required for the *E. coli* IscS protein, for which numerous interaction partners have been identified. IscS was additionally shown to interact with a number of acceptor proteins for delivery of sulfur including the involvement of (i) IscU, CyaY, Fdx and IscX for Fe-S cluster formation, (ii) TusA for either the (c)mnm⁵s²U34 modifications of tRNA or the biosynthesis of Moco, and (iii) ThiI for the synthesis of thiamine or the s⁴U8 modification of tRNA [36, 37]. Different binding sites for some of these molecules were mapped [38], ensuring either the simultaneous binding or a competitive binding on overlapping binding sites. While the binding site of FdhD on IscS has
not been mapped so far, the stimulating effect on the activity of IscS might ensure that the sul-
fur is shuttled in the direction of bis-MGD sulfuration after binding of FdhD. It is possible that
the stimulation of the L-cysteine desulfurase activity of IscS is triggered by a specific interaction
with the disordered loop at the active site of IscS [39]. This has to be investigated in detail in
future studies, i.e. by solving the co-crystal structure of IscS with FdhD. The reason why other
FdhD-like chaperones do not comprise the enhancing effect on the L-cysteine desulfurase
activity might be based on differences in the interaction site. As shown for R. capsulatus FdsC
in this study, the interaction between both enzymes is weaker, since protein complexes were
not copurified (Fog 8). However, since the bis-MGD cofactor on FdsC was identified to be
present in its sulfureted form, E. coli IscS is able to deliver the sulfur to FDsC. Since in R. capsu-
latus three L-cysteine desulfurases are present which in total might have a lower number of
interacting proteins [33], a more transient interaction might be sufficient to shuttle the sulfur
directly in the direction of bis-MGD sulfuration. Since the specific roles of the three L-cysteine
desulfurases have not been characterized in R. capsulatus so far, this assumption is of course
speculative. In the organisms in which the CXXC motif is completely absent in the FdhD-like
chaperone, the sulfur transfer pathway and the involvement of the corresponding L-cysteine
desulfurase might differ. In organisms like Bacillus subtilis, numerous L-cysteine desulfurases
were identified with a specific role for the synthesis of one sulfur-containing molecule
[40]. Since these proteins are consequently not involved for sulfur transfer to several accep-
tor proteins and are specific for one particular pathway [41], the mode of interaction and
the involvement of an additional cysteine on FdhD might serve a different role. Thus, for
each organism the interaction of the FdhD-like chaperone with its partner L-cysteine desul-
furase needs to be investigated in detail. Unfortunately, so far in many organisms the
FdhD-like enzyme, the partner FDH enzymes and the corresponding L-cysteine desul-
furases are largely uncharacterized.

In total, we conclude that the chaperones FdsC and FdhD have exchangeable roles.
Both FdsC and FdhD can substitute each other in the maturation of the FDH partner
enzyme(s) (Fig 9). This was not expected, since it has been suggested before that the sys-
tem-specific chaperone coevolved with its corresponding partner molybdoenzyme. E. coli
FdhD was shown before to be the chaperone for all three E. coli FDH enzymes, FdhF,
(FdnGHI)₃ and (FdoGHI)₃ [20]. The crystal structures of FdhF and (FdnGHI)₃ were
solved and revealed to be highly identical [28, 42]. The crystal structure of R. capsulatus
(FdsGBA)₂ is not available so far. From amino acid sequence alignments, however, it
becomes clear that the bis-MGD binding domain of FdsA is highly conserved and might
provide a similar fold as the E. coli enzymes. Due to the high conservation of the bis-MGD
containing domain of all FDH enzymes both FdhD and FdsC are expected to bind to a
common motif present in all FDH enzymes, which then provides the basis of the
exchangeable roles of both chaperones. The specific binding site of the FdhD-like chaper-
one on its FDH target enzyme in addition to the mechanism of bis-MGD insertion needs
to be clarified in more detail in future studies.

**Materials and methods**

**Bacterial strains, plasmids, media and growth conditions**

Bacterial strains, plasmids and primers are listed in Table 1. *E. coli* BW25113 Δ*fdhD* cells were obtained from the Keio collection [43]. Cells were grown at 30°C under anaerobic conditions in LB medium containing 1 mM molybdate, 20 µM isopropyl β-D-1-thiogalactopyranoside (IPTG), 150 µg/mL ampicillin for 24 h. For site-directed mutagenesis and construction of the FdhD variants C121A, C124A and C121A/C124A, the expression vector pNB10 [21] was used
<figure><img src="image_16.png" /><figcaption>Fig 9. Model for bis-MGD sulfuration and insertion into the target enzymes by involvement of FdsC and FdhD. FdsC and FdhC bind the bis-MGD cofactor synthesized by the Moco biosynthesis machinery and provide the sulfurated bis-MGD for the FDHs of *R. capsulatus* and *E. coli*, respectively. Both chaperones are involved in the sulfuration of bis-MGD. FdhD interacts with the L-cysteine desulfurase IscS. The corresponding L-cysteine desulfurase for *R. capsulatus* has not been identified to date. FdsC and FdhD are functional homologues since the can replace each other in the maturation of the respective FDH.</figcaption></figure>

https://doi.org/10.1371/journal.pone.0201935.g009

as a template and base-pair exchanges were introduced by polymerase chain reaction mutagenesis. In general, to guarantee comparable growth conditions for the *E. coli* and *R. capsulatus* used, and protein expression was generally carried out at 30°C.

## Expression and purification of recombinant proteins

FdsC, its variants (from plasmids pTHfds02, pNBfds07, pNBfds08, pNBfds09, Table 1), FdhD and its variants (from Plasmids pNB10, pNB11, pNB12, pNB13, Table 1) were purified as described before for FdsC [21]. For co-expression of *R. capsulatus* FDH (FdsGBA + FdsD) with FdsC, FdhD or its cysteine variants, *E. coli* BW25113 Δ*fdhD* cells were transformed with pTHfds15 and either pTHfds14 pNB14, pNB15 pNBfds16 or pNBfds17. *R. capsulatus* (FdsGBA)₂ was expressed and purified as described previously [30]. Apo-TorA was purified from soluble extract of RK5200 strain transformed with pTorA grown at 30°C 130 rpm by Ni-
<table><caption>Table 1. Bacterial strains and plasmids used in this study.</caption><thead><tr><th></th><th>genotype/characterization</th><th>refs</th></tr></thead><tbody><tr><td><strong>plasmids</strong></td><td></td><td></td></tr><tr><td>pNBfdsC04</td><td>fdsC<sup>C104A</sup> gene cloned into NcoI/SacI site of pTrcHis, Amp<sup>R</sup></td><td>[21]</td></tr><tr><td>pNBfdsC05</td><td>fdsC<sup>C107A</sup> gene cloned into NcoI/SacI site of pTrcHis, Amp<sup>R</sup></td><td>[21]</td></tr><tr><td>pNBfdsC06</td><td>fdsC<sup>C104A/C107A</sup> gene cloned into NcoI/SacI site of pTrcHis, Amp<sup>R</sup></td><td>[21]</td></tr><tr><td>pNBfdsC07</td><td>fdsC<sup>C104A</sup> gene cloned into NdeI/SacI site of pTrcHis, expressing N-terminally tagged His<sub>6</sub>.FdsC<sup>C104A</sup>, Amp<sup>R</sup></td><td>this work</td></tr><tr><td>pNBfdsC08</td><td>fdsC<sup>C107A</sup> gene cloned into NdeI/SacI site of pTrcHis, expressing N-terminally tagged His<sub>6</sub>.FdsC<sup>C107A</sup>, Amp<sup>R</sup></td><td>this work</td></tr><tr><td>pNBfdsC09</td><td>fdsC<sup>C104A/C107A</sup> gene cloned into NdeI/SacI site of pTrcHis, expressing N-terminally tagged His<sub>6</sub>.FdsC<sup>C104A/C107A</sup>, Amp<sup>R</sup></td><td>this work</td></tr><tr><td>pNB10</td><td>fdhD gene cloned into NdeI/SalI site of pTrcHis, expressing N-terminally tagged His<sub>6</sub>-FdhD, Amp<sup>R</sup></td><td>this work</td></tr><tr><td>pNB11</td><td>fdhD<sup>C121A</sup> gene cloned into NdeI/SalI site of pTrcHis, expressing N-terminally tagged His<sub>6</sub>-FdhD<sup>C121A</sup>, Amp<sup>R</sup></td><td>this work</td></tr><tr><td>pNB12</td><td>fdhD<sup>C124A</sup> gene cloned into NdeI/SalI site of pTrcHis, expressing N-terminally tagged His<sub>6</sub>-FdhD<sup>C124A</sup>, Amp<sup>R</sup></td><td>this work</td></tr><tr><td>pNB13</td><td>fdhD<sup>C121A/C124A</sup> gene cloned into NdeI/SalI site of pTrcHis, expressing N-terminally tagged His<sub>6</sub>-FdhD<sup>C121A/C124A</sup>, Amp<sup>R</sup></td><td>this work</td></tr><tr><td>pNB14</td><td>fdhD gene cloned into NcoI/SalI site of pTrcHis, Amp<sup>R</sup></td><td>this work</td></tr><tr><td>pNB15</td><td>fdhD<sup>C121A</sup> gene cloned into NcoI/SalI site of pTrcHis, Amp<sup>R</sup></td><td>this work</td></tr><tr><td>pNB16</td><td>fdhD<sup>C124A</sup> gene cloned into NcoI/SalI site of pTrcHis, Amp<sup>R</sup></td><td>this work</td></tr><tr><td>pNB17</td><td>fdhD<sup>C121A/C124A</sup> gene cloned into NcoI/SalI site of pTrcHis, Amp<sup>R</sup></td><td>this work</td></tr><tr><td>pMN20</td><td>nifS4 gene cloned into NdeI/XhoI site of pET28a, Km<sup>R</sup></td><td>[33]</td></tr><tr><td>pMN51</td><td>nifS3 gene cloned into NheI/XhoI site of pET28a, Km<sup>R</sup></td><td>[33]</td></tr><tr><td>pMN54</td><td>NifS2Δ1-188 gene cloned into NdeI/SalI site of pET28a, Km<sup>R</sup></td><td>[33]</td></tr><tr><td>pSL209</td><td>iscS gene cloned into NcoI/BamHI sites of pET15b, Amp<sup>R</sup></td><td>[54]</td></tr><tr><td>pTHfds02</td><td>fdsC gene cloned into NheI/SacI site of pTrcHis, expressing N-terminally tagged His<sub>6</sub>.FdsC, Amp<sup>R</sup></td><td>[30]</td></tr><tr><td>pTHfds14</td><td>fdsC gene cloned into NcoI/SacI site of pTrcHis, Amp<sup>R</sup></td><td>this work</td></tr><tr><td>pTHfds15</td><td>fdsGBA genes cloned into SacI-/SalI site (MSCI) and fdsD gene cloned into NdeI/XhoI (MSCII) site of pACYC-duet1, co-expressing N-terminally tagged His<sub>6</sub>-FdsGBA and FdsD, Cm<sup>R</sup></td><td>[21]</td></tr><tr><td>pTorA</td><td>torA gene cloned into XbaI/HindIII site of JF119EH, Amp<sup>R</sup></td><td>[55]</td></tr><tr><td><strong>strains</strong></td><td></td><td></td></tr><tr><td>ΔfdhD (DE3)</td><td>JW3866-3, Δ(araD-araB)567, ΔlacZ4787(::rrnB-3), λ-, rph-1, ΔfdhD758::kan, Δ(rhaD-rhaB)568, hsdR514 λ(DE3)</td><td>[43]</td></tr><tr><td>ΔfdoG</td><td>JW3865-2, F-, Δ(araD-araB)567, ΔlacZ4787(::rrnB-3), λ-, rph-1, ΔfdoG757::kan, Δ(rhaD-rhaB)568, hsdR514</td><td>[43]</td></tr><tr><td>BL21 (DE3)</td><td>ompT gal dcm lon hsdSB(rB- mB-) λ(DE3)</td><td>Novagen</td></tr><tr><td>RK5200</td><td>F- araD139Δ(argF lac)U169 deoC1 flbB5201 gyrA219 relA1 rpsL150 non-9 ptsF25 chlA200::Mucts ΔmoaA::Km</td><td>[56]</td></tr></tbody></table>

https://doi.org/10.1371/journal.pone.0201935.t001

NTA chromatography as described previously [44, 45]. Expression and purification of IscS (pSL209), NifS2Δ1-188 (pNM54), NifS3 (pNM51), and NifS4 (pNM20) were performed as described previously [33, 46].
## Interaction study by size exclusion chromatography

20 µM IscS was mixed with 40 µM FdhD or FdsC and incubated in 100 mM potassium phosphate buffer, pH 8.0 for 30 min at room temperature. The proteins were loaded onto a Superdex 200 column equilibrated in the same buffer. Proteins were separated at a flow rate of 1 mL/min, and the elution profile was recorded at 280 nm. The proteins in the elution fractions were separated by 15% SDS-PAGE.

## Cofactor analysis

200 µM of purified chaperone (FdsC, FdsC-C104A, FdsC-C107A, FdsC-C104A/C107A, FdhD, FdhD-C121A, FdhD-C124A, FdhD-C121A/C124A) or 2 µM reconstituted TMAO reductase were incubated in acidic $I_2$/KI-HCl at room temperature to release and oxidize bis-MGD to its fluorescent degradation product FormA-GMP. Form A-GMP was separated by a C18 reversed-phase high pressure liquid chromatography column (4.6 x 250-mm Hypersil ODS, 5-µm particle size) after the method described previously [47]. Metal analysis was performed using PerkinElmer Life Sciences Optima 2100DV inductively coupled plasma optical emission spectrometer as described earlier [33].

## In vitro reconstitution of apo-TorA

8 µM apo-TorA was incubated with 50 µM FdsC, FdsC-C104A, FdsC-C107A, FdsC-C104A/C107A, FdhD, FdhD-C121A, FdhD-C124A, or FdhD-C121A/C124A in 100 mM potassium phosphate buffer for 2 hours at room temperature under anaerobic conditions. TMAO reductase activity was measured under anaerobic condition following the oxidation of reduced benzyl viologen in the presence of 0.1 µM TMAO ($ε_{NADH} = 7,400 M^{-1}cm^{-1}$) as described by Kaufmann et al. [35]. For detailed analysis of reconstituted TMAO-reductase with FdsC and FdhD, TorA was incubated with FdsC or FdhD for 7 hours at 37°C and purified from the mixtures by gel filtration on a Superdex 200 10/300GL column (GE Healthcare). The containing TorA were combined and used for cofactor analysis and activity measurements.

## Enzyme assays

*R. capsulatus* ($FdsGBA_2$) activity was measured with 120–300 nM enzyme in the assay by recording the production of NADH at 340 nm ($ε_{NADH} = 6,220 M^{-1}cm^{-1}$) as described previously [30].

L-cysteine desulfurase activity of IscS, NifS2Δ1-188, NifS3, and NifS4 was measured by determination of the rate of sulfide production as described previously [33, 48]. 1 µM L-cysteine desulfurase was mixed in a 1:2 ratio with FdsC, FdhD, and SufE in a total volume of 480 µl containing 100 mM potassium phosphate, 200 mM NaCl, 10 µM pyridoxal phosphate, and 1 mM dithiothreitol, pH 8.0. After incubation for 10 min at 30°C, the reactions were stopped by the addition of 60 µl of 20 mM N,N-dimethyl-p-phenylenediamine in 7.2 M HCl and 60 µl of 30 mM $FeCl_3$ in 1.2 M HCl. After additional incubation for 20 min at 30°C, precipitated protein was removed by centrifugation, and methylene blue was measured at 670 nm. A standard curve was generated using sodium sulfide in a range of 0–25 µM.

## In-gel staining of FdhO activity

For detection of FdhO activity, BW25113 cells, *ΔfdhO* cells, *ΔfdhD* cells and *ΔfdhD* cells complemented with plasmids coding for FdhD and FdsC as well as for their cysteine variants (pNB14, pNB15, pNB16, pNB17, pTHfds14, pNBfds04, pNBfds05, pNBfds06) were grown anaerobically in 55 ml LB medium supplemented with 10 µM sodium molybdate, 0.2% sodium
formate and 20 µM IPTG for 16 hours (starting with OD$_{600}$ of 0.05). After harvesting cells were resuspended in 1.5 ml 50 mM potassium phosphate buffer, pH 6.8 containing 1% (v/v) Triton X-100 and were disrupted by sonication. Cell debris were removed by centrifugation and the protein concentration of the supernatant was measured by the BCA method as described by the manufacturer (Micro BCA Protein Assay Reagent Kit, Pierce). 40–60 µg of total proteins were separated by non-denaturing PAGE using Tris-Glycin with 0.1% (v/v) as running buffer. Detection of FdhO activity was performed as described using 5 mM sodium formate as substrate, 1 mM nitroblue tetrazolium as electron acceptor and 0.5 mM phenazine methosulfate as mediator in 50 mM potassium phosphate, pH 6.8 [49]. The gels were incubated in staining solution for 30 min at 37°C under aerobic conditions.

## Phylogenetic data analysis

Amino acid sequences of 41 FdhD-like sequences from various taxa (bacteria: actinobacteria, firmicutes, cyanobacteria and proteobacteria, and archaea, Table 2) were extracted from ENSEMBLE and amino acid sequence alignments were created with Muscle [50]. The FdhD/

<table><caption>**Table 2. Protein accession numbers used for the phylogenetic tree.** NCBI protein accession numbers for each FdhD/fdsC amino acid sequence in alphabetical order of the organisms.</caption><thead><tr><td>Accession Number</td><td>Organism</td></tr></thead><tbody><tr><td>ABR74626</td><td>Actinobacillus succinogenes</td></tr><tr><td>KAJ33787</td><td>Agrobacterium tumefaciens</td></tr><tr><td>ELT44815</td><td>Arthrobacter nitrophenolicus</td></tr><tr><td>ACO76640</td><td>Azotobacter vinelandii</td></tr><tr><td>CAB15688</td><td>Bacillus subtilis</td></tr><tr><td>KHS53770</td><td>Brevibacterium linens</td></tr><tr><td>AFS79899</td><td>Clostridium acidurici 9a</td></tr><tr><td>CAF19235</td><td>Corynebacterium glutamicum</td></tr><tr><td>SCV00736</td><td>Cupravidus necator</td></tr><tr><td>AGW12659</td><td>Desulfovibrio gigas</td></tr><tr><td>YP_009800</td><td>Desulfovibrio vulgaris str. Hildenborough</td></tr><tr><td>AIZ53701</td><td>Escherichia coli K13</td></tr><tr><td>BAD74746</td><td>Geobacillus kaustophilus</td></tr><tr><td>ACY12706</td><td>Haliangium ochraceum</td></tr><tr><td>ADJ24264</td><td>Hyphomicrobium denitrificans</td></tr><tr><td>AEB94994</td><td>Metallosphaera cuprina</td></tr><tr><td>AIS30994</td><td>Methanobacterium formicicum</td></tr><tr><td>ABQ86500</td><td>Methanobrevibacter smithii</td></tr><tr><td>ABQ87597</td><td>Methanobrevibacter smithii</td></tr><tr><td>WP_048064104</td><td>Methanococcus maripaludis</td></tr><tr><td>ABR54454</td><td>Methanococcus vannielii SB</td></tr><tr><td>ACK84408</td><td>Methylobacterium extorquens CM5</td></tr><tr><td>ACT50427</td><td>Methylovorus glucosetrophus</td></tr><tr><td>CDO44190</td><td>Mycobacterium bovis</td></tr><tr><td>ACV78244</td><td>Nakamurella multipartita</td></tr><tr><td>ACC79387</td><td>Nostoc punctiforme</td></tr><tr><td>BAF60791</td><td>Pelotomaculum thermopropionicum</td></tr><tr><td>ABG41355</td><td>Pseudoalteromonas atlantica</td></tr><tr><td>ABR85479</td><td>Pseudomonas aeruginosa</td></tr></tbody></table>

(Continued)
<table><caption>**Table 2.** (Continued)</caption><thead><tr><td>Accession Number</td><td>Organism</td></tr></thead><tbody><tr><td>ETD87413</td><td>Rhodobacter capsulatus</td></tr><tr><td>OPF94773</td><td>Rhodopseudomonas palustris</td></tr><tr><td>ABI70605</td><td>Shewanella frigidimarina</td></tr><tr><td>AAP18794</td><td>Shigella flexneri</td></tr><tr><td>ADD40018</td><td>Stackebrandtia nassauensis</td></tr><tr><td>CAI81841</td><td>Staphylococcus aureus</td></tr><tr><td>BAG17950</td><td>Streptomyces griseus</td></tr><tr><td>BAB65040</td><td>Sulfolobus tokodaii</td></tr><tr><td>ABK15737</td><td>Syntrophobacter fumaroxidans MPOB</td></tr><tr><td>YP_845333.1</td><td>Syntrophobacter fumaroxidans MPOB</td></tr><tr><td>ACP05777</td><td>Vibrio cholerae</td></tr><tr><td>WP_011137999</td><td>Wolinella succinogenes</td></tr></tbody></table>

https://doi.org/10.1371/journal.pone.0201935.t002

FdsC amino acid sequences were selected from different organisms including characterized and uncharacterized FDH proteins to have a broad selection of sequences which was supposed to enable the analysis whether a splitting in the different groups of FDHs occurs [22]. The tree was constructed using the Maximum Likelihood method based on the Dayhoff matrix based model [51] The tree is drawn to scale, with branch lengths in the same units as those of the evolutionary distances used to infer the phylogenetic tree. The evolutionary distances were computed using the Dayhoff matrix based method and are in the units of the number of amino acid substitutions per site [52]. The scale bar indicates 0.2 substitutions per site. Numbers near branches indicate the bootstrap proportion for 100 replica using the same method. All positions containing gaps and missing data were eliminated. There was a total of 174 positions in the final dataset. Evolutionary analyses were conducted with the software MEGA7 [53].

## Acknowledgments

Kristin Gennermann is thanked for help with the phylogenetic data analysis.

## Author Contributions

**Conceptualization:** Chantal Iobbi-Nivol, Silke Leimkühler.

**Data curation:** Nadine Schwanhold, Angelika Lehmann.

**Formal analysis:** Nadine Schwanhold, Silke Leimkühler.

**Funding acquisition:** Silke Leimkühler.

**Investigation:** Nadine Schwanhold, Angelika Lehmann.

**Methodology:** Chantal Iobbi-Nivol, Angelika Lehmann.

**Project administration:** Silke Leimkühler.

**Resources:** Chantal Iobbi-Nivol, Silke Leimkühler.

**Supervision:** Chantal Iobbi-Nivol, Silke Leimkühler.

**Validation:** Nadine Schwanhold, Silke Leimkühler.

**Visualization:** Angelika Lehmann.
**Writing – original draft:** Nadine Schwanhold, Silke Leimkühler.

**Writing – review & editing:** Chantal Iobbi-Nivol, Silke Leimkühler.

## References

1. Leimkühler S, Iobbi-Nivol C. Bacterial molybdoenzymes: old enzymes for new purposes. FEMS micro-biology reviews. 2016; 40(1):1–18. Epub 2015/10/16. https://doi.org/10.1093/femsre/fuv043 PMID: 26468212.

2. Hille R. The mononuclear molybdenum enzymes. Chemical Rev. 1996; 96:2757–816.

3. Palmer T, Vasishta A, Whitty PW, Boxer DH. Isolation of protein FA, a product of the *mob* locus required for molybdenum cofactor biosynthesis in *Escherichia coli*. Eur J Biochem. 1994; 222:687–92. PMID: 8020507

4. Yokoyama K, Leimkühler S. The role of FeS clusters for molybdenum cofactor biosynthesis and molybdoenzymes in bacteria. Biochim Biophys Acta. 2014. Epub 2014/10/01. https://doi.org/10.1016/j.bbamcr.2014.09.021 PMID: 25268953.

5. Kisker C, Schindelin H, Rees DC. Molybdenum-cofactor-containing enzymes: structure and mechanism. Ann Rev Biochem. 1997; 66:233–67. https://doi.org/10.1146/annurev.biochem.66.1.233 PMID: 9242907

6. Blasco F, Dos Santos JP, Magalon A, Frixon C, Guigliarelli B, Santini CL, et al. NarJ is a specific chaperone required for molybdenum cofactor assembly in nitrate reductase A of *Escherichia coli*. Mol Micro-biol. 1998; 28:435–47. PMID: 9632249

7. Genest O, Mejean V, Iobbi-Nivol C. Multiple roles of TorD-like chaperones in the biogenesis of molybdoenzymes. FEMS microbiology letters. 2009; 297(1):1–9. Epub 2009/06/13. https://doi.org/10.1111/j.1574-6968.2009.01660.x PMID: 19519768.

8. Genest O, Neumann M, Seduk F, Stocklein W, Mejean V, Leimkühler S, et al. Dedicated metallochaperone connects apoenzyme and molybdenum cofactor biosynthesis components. J Biol Chem. 2008; 283 (31):21433–40. Epub 2008/06/05. https://doi.org/10.1074/jbc.M802954200 PMID: 18522945.

9. Chan CS, Chang L, Rommens KL, Turner RJ. Differential Interactions between Tat-specific redox enzyme peptides and their chaperones. J Bacteriol. 2009; 191(7):2091–101. Epub 2009/01/20. https://doi.org/10.1128/JB.00949-08 PMID: 19151138; PubMed Central PMCID: PMC2655534.

10. Chan CS, Chang L, Winstone TM, Turner RJ. Comparing system-specific chaperone interactions with their Tat dependent redox enzyme substrates. FEBS Lett. 2010; 584(22):4553–8. Epub 2010/10/27. https://doi.org/10.1016/j.febslet.2010.10.043 PMID: 20974141; PubMed Central PMCID: PMC3285697.

11. Lorenzi M, Sylvi L, Gerbaud G, Mileo E, Halgand F, Walburger A, et al. Conformational selection underlies recognition of a molybdoenzyme by its dedicated chaperone. PloS one. 2012; 7(11):e49523. Epub 2012/11/28. https://doi.org/10.1371/journal.pone.0049523 PMID: 23185350; PubMed Central PMCID: PMC3501500.

12. Neumann M, Leimkühler S. The role of system-specific molecular chaperones in the maturation of molybdoenzymes in bacteria. Biochemistry research international. 2011; 2011:850924. Epub 2010/12/15. https://doi.org/10.1155/2011/850924 PMID: 21151514; PubMed Central PMCID: PMC2997495.

13. Turner RJ, Papish AL, Sargent F. Sequence analysis of bacterial redox enzyme maturation proteins (REMPs). Canadian journal of microbiology. 2004; 50(4):225–38. Epub 2004/06/24. https://doi.org/10.1139/w03-117 PMID: 15213747.

14. Bay DC, Chan CS, Turner RJ. NarJ subfamily system specific chaperone diversity and evolution is directed by respiratory enzyme associations. BMC Evol Biol. 2015; 15:110. Epub 2015/06/13. https://doi.org/10.1186/s12862-015-0412-3 PMID: 26067063; PubMed Central PMCID: PMCPMC4464133.

15. Lemaire ON, Bouillet S, Mejean V, Iobbi-Nivol C, Genest O. Chaperones in maturation of molybdoenzymes: Why specific is better than general? Bioengineered. 2017; 8(2):133–6. Epub 2016/09/01. https://doi.org/10.1080/21655979.2016.1218579 PMID: 27580420; PubMed Central PMCID: PMCPMC5398582.

16. Ilbert M, Mejean V, Iobbi-Nivol C. Functional and structural analysis of members of the TorD family, a large chaperone family dedicated to molybdoproteins. Microbiology. 2004; 150(Pt 4):935–43. https://doi.org/10.1099/mic.0.26909-0 PMID: 15073303.

17. Chan CS, Bay DC, Leach TG, Winstone TM, Kuzniatsova L, Tran VA, et al. 'Come into the fold': A comparative analysis of bacterial redox enzyme maturation protein members of the NarJ subfamily. Biochim Biophys Acta. 2014; 1838(12):2971–84. Epub 2014/08/27. https://doi.org/10.1016/j.bbamem.2014.08.020 PMID: 25157671.
18. Cherak SJ, Turner RJ. Assembly pathway of a bacterial complex iron sulfur molybdoenzyme. Biomol Concepts. 2017; 8(3-4):155-67. Epub 2017/07/09. https://doi.org/10.1515/bmc-2017-0011 PMID: 28688222.

19. Arnoux P, Ruppelt C, Oudouhou F, Lavergne J, Siponen MI, Toci R, et al. Sulphur shuttling across a chaperone during molybdenum cofactor maturation. Nature communications. 2015; 6:6148. Epub 2015/02/05. https://doi.org/10.1038/ncomms7148 PMID: 25649206.

20. Thome R, Gust A, Toci R, Mendel R, Bittner F, Magalon A, et al. A sulfurtransferase is essential for activity of formate dehydrogenases in Escherichia coli. J Biol Chem. 2012; 287(7):4671-8. Epub 2011/12/24. https://doi.org/10.1074/jbc.M111.327122 PMID: 22194618; PubMed Central PMCID: PMC3281601.

21. Böhmer N, Hartmann T, Leimkühler S. The chaperone FdsC for Rhodobacter capsulatus formate dehydrogenase binds the bis-molybdopterin guanine dinucleotide cofactor. FEBS Lett. 2014; 588(4):531-7. Epub 2014/01/22. https://doi.org/10.1016/j.febslet.2013.12.033 PMID: 24444607.

22. Hartmann T, Schwanhold N, Leimkühler S. Assembly and catalysis of molybdenum or tungsten-containing formate dehydrogenases from bacteria. Biochim Biophys Acta. 2014. Epub 2014/12/17. https://doi.org/10.1016/j.bkapap.2014.12.006 PMID: 25514355.

23. Hille R, Hall J, Basu P. The mononuclear molybdenum enzymes. Chemical reviews. 2014; 114 (7):3963-4038. Epub 2014/01/29. https://doi.org/10.1021/cr400443z PMID: 24467397; PubMed Central PMCID: PMC4080432.

24. Benoit S, Abaibou H, Mandrand-Berthelot MA. Topological analysis of the aerobic membrane-bound formate dehydrogenase of Escherichia coli. J Bacteriol. 1998; 180(24):6625-34. Epub 1998/12/16. PMID: 9852007; PubMed Central PMCID: PMC107766.

25. Jormakka M, Tornroth S, Byrne B, Iwata S. Molecular basis of proton motive force generation: structure of formate dehydrogenase-N. Science. 2002; 295(5561):1863-8. Epub 2002/03/09. https://doi.org/10.1126/science.1068186 PMID: 11884747.

26. Sawers G. The hydrogenases and formate dehydrogenases of Escherichia coli. Antonie Van Leeuwenhoek. 1994; 66(1-3):57-88. Epub 1994/01/01. PMID: 7747941.

27. Jormakka M, Byrne B, Iwata S. Formate dehydrogenase—a versatile enzyme in changing environments. Curr Opin Struct Biol. 2003; 13(4):418-23. PMID: 12948771.

28. Boyington JC, Gladyshev VN, Khangulov SV, Stadtman TC, Sun PD. Crystal structure of formate dehydrogenase H: catalysis involving Mo, molybdopterin, selenocysteine, and an Fe₄S₄ cluster. Science. 1997; 275:1305-8. PMID: 9036855

29. Sawers G, Suppmann B. Anaerobic induction of pyruvate formate-lyase gene expression is mediated by the ArcA and FNR proteins. J Bacteriol. 1992; 174(11):3474-8. Epub 1992/06/01. PMID: 1592804; PubMed Central PMCID: PMCPMC206030.

30. Hartmann T, Leimkühler S. The oxygen-tolerant and NAD(+) -dependent formate dehydrogenase from Rhodobacter capsulatus is able to catalyze the reduction of CO2 to formate. The FEBS journal. 2013; 280(23):6083-96. Epub 2013/09/17. https://doi.org/10.1111/febs.12528 PMID: 24034888.

31. Hartmann T, Schrapers P, Utesch T, Nimtz M, Rippers Y, Dau H, et al. The Molybdenum Active Site of Formate Dehydrogenase Is Capable of Catalyzing C-H Bond Cleavage and Oxygen Atom Transfer Reactions. Biochemistry. 2016; 55(16):2381-9. Epub 2016/04/08. https://doi.org/10.1021/acs.biochem.6b00002 PMID: 27054466.

32. Schrapers P, Hartmann T, Kositzki R, Dau H, Reschke S, Schulzke C, et al. Sulfido and cysteine ligation changes at the molybdenum cofactor during substrate conversion by formate dehydrogenase (FDH) from Rhodobacter capsulatus. Inorg Chem. 2015; 54(7):3260-71. Epub 2015/03/25. https://doi.org/10.1021/ic502880y PMID: 25803130.

33. Neumann M, Stöcklein W, Walburger A, Magalon A, Leimkühler S. Identification of a Rhodobacter capsulatus L-cysteine desulfurase that sulfurates the molybdenum cofactor when bound to XdhC and before its insertion into xanthine dehydrogenase. Biochemistry. 2007; 46(33):9586-95. https://doi.org/10.1021/bi700630p PMID: 17649978.

34. Harel A, Häggblom MM, Falkowski PG, Yee N. Evolution of prokaryotic respiratory molybdoenzymes and the frequency of their genomic co-occurrence. FEMS Microbiol Ecol. 2016; 92:187.

35. Kaufmann P, Duffus BR, Mitrova B, Iobbi-Nivol C, Teutloff C, Nimtz M, et al. Modulating the Molybdenum Coordination Sphere of Escherichia coli Trimethylamine N-Oxide Reductase. Biochemistry. 2018; 57(7):1130-43. Epub 2018/01/16. https://doi.org/10.1021/acs.biochem.7b01108 PMID: 29334455.

36. Hidese R, Mihara H, Esaki N. Bacterial cysteine desulfurases: versatile key players in biosynthetic pathways of sulfur-containing biofactors. Applied microbiology and biotechnology. 2011; 91(1):47-61. Epub 2011/05/24. https://doi.org/10.1007/s00253-011-3336-x PMID: 21603932.
37. Leimkühler S, Bühning M, Beilschmidt L. Shared Sulfur Mobilization Routes for tRNA Thiolation and Molybdenum Cofactor Biosynthesis in Prokaryotes and Eukaryotes. Biomolecules. 2017; 7(1). Epub 2017/01/19. https://doi.org/10.3390/biom7010005 PMID: 28098827; PubMed Central PMCID: PMC5372717.

38. Shi R, Proteau A, Villarroya M, Moukadiri I, Zhang L, Trempe JF, et al. Structural basis for Fe-S cluster assembly and tRNA thiolation mediated by IscS protein-protein interactions. PLoS biology. 2010; 8(4): e1000354. Epub 2010/04/21. https://doi.org/10.1371/journal.pbio.1000354 PMID: 20404999; PubMed Central PMCID: PMC2854127.

39. Cupp-Vickery JR, Urbina H, Vickery LE. Crystal structure of IscS, a cysteine desulfurase from Escherichia coli. Journal of molecular biology. 2003; 330(5):1049-59. Epub 2003/07/16. PMID: 12860127.

40. Black KA, Dos Santos PC. Shared-intermediates in the biosynthesis of thio-cofactors: Mechanism and functions of cysteine desulfurases and sulfur acceptors. Biochim Biophys Acta. 2015; 1853(6):1470-80. Epub 2014/12/03. https://doi.org/10.1016/j.bbamcr.2014.10.018 PMID: 25447671.

41. Black KA, Dos Santos PC. Abbreviated Pathway for Biosynthesis of 2-Thiouridine in Bacillus subtilis. J Bacteriol. 2015; 197(11):1952-62. Epub 2015/04/01. https://doi.org/10.1128/JB.02625-14 PMID: 25825430; PubMed Central PMCID: PMCPMC4420905.

42. Jormakka M, Tornroth S, Abramson J, Byrne B, Iwata S. Purification and crystallization of the respiratory complex formate dehydrogenase-N from Escherichia coli. Acta crystallographica Section D, Biological crystallography. 2002; 58(Pt 1):160-2. Epub 2001/12/26. PMID: 11752799.

43. Baba T, Ara T, Hasegawa M, Takai Y, Okumura Y, Baba M, et al. Construction of Escherichia coli K-12 in-frame, single-gene knockout mutants: the Keio collection. Molecular systems biology. 2006; 2:2006 0008. Epub 2006/06/02. https://doi.org/10.1038/msb4100050 PMID: 16738554; PubMed Central PMCID: PMC1681482.

44. Pommier J, Mejean V, Giordano G, lobbi-Nivol C. TorD, a cytoplasmic chaperone that interacts with the unfolded trimethylamine N-oxide reductase enzyme (TorA) in Escherichia coli. J Biol Chem. 1998; 273 (26):16615-20. Epub 1998/06/20. PMID: 9632735.

45. Reschke S, Sigfridsson KG, Kaufmann P, Leidel N, Horn S, Gast K, et al. Identification of a Bis-molybdopterin Intermediate in Molybdenum Cofactor Biosynthesis in Escherichia coli. J Biol Chem. 2013; 288 (41):29736-45. Epub 2013/09/05. https://doi.org/10.1074/jbc.M113.497453 PMID: 24003231; PubMed Central PMCID: PMC3795271.

46. Leimkühler S, Rajagopalan KV. An *Escherichia coli* NifS-like sulfurtransferase is required for the transfer of cysteine sulfur in the *in vitro* synthesis of molybdopterin from precursor Z. J Biol Chem. 2001; 276:22024-31. https://doi.org/10.1074/jbc.M102072200 PMID: 11290749

47. Johnson JL, Hainline BE, Rajagopalan KV, Arison BH. The pterin component of the molybdenum cofactor. Structural characterization of two fluorescent derivatives. J Biol Chem. 1984; 259:5414-22. PMID: 6546929

48. Urbina HD, Silberg JJ, Hoff KG, Vickery LE. Transfer of sulfur from IscS to IscU during Fe/S cluster assembly. J Biol Chem. 2001; 276(48):44521-6. https://doi.org/10.1074/jbc.M106907200 PMID: 11577100.

49. Enoch HG, Lester RL. The purification and properties of formate dehydrogenase and nitrate reductase from Escherichia coli. J Biol Chem. 1975; 250(17):6693-705. Epub 1975/09/10. PMID: 1099093.

50. Edgar RC. MUSCLE: a multiple sequence alignment method with reduced time and space complexity. BMC Bioinformatics. 2004; 5:113. Epub 2004/08/21. https://doi.org/10.1186/1471-2105-5-113 PMID: 15318951; PubMed Central PMCID: PMCPMC517706.

51. Schwarz R, Dayhoff MO. Matrices for detecting distant relationships. Atlas of protein sequences. 1979:353-58.

52. Dayhoff MO, Orcutt BC. Methods for identifying proteins by using partial sequences. Proc Natl Acad Sci U S A. 1979; 76(5):2170-4. Epub 1979/05/01. PMID: 287054; PubMed Central PMCID: PMCPMC383559.

53. Kumar S, Stecher G, Tamura K. MEGA7: Molecular Evolutionary Genetics Analysis Version 7.0 for Bigger Datasets. Mol Biol Evol. 2016; 33(7):1870-4. Epub 2016/03/24. https://doi.org/10.1093/molbev/msw054 PMID: 27004904.

54. Leimkühler S, Rajagopalan KV. A sulfurtransferase is required in the transfer of cysteine sulfur in the *in vitro* synthesis of molybdopterin from precursor Z in Escherichia coli. J Biol Chem. 2001; 276 (25):22024-31. Epub 2001/04/06. https://doi.org/10.1074/jbc.M102072200 PMID: 11290749.

55. Ilbert M, Mejean V, Giudici-Orticoni MT, Samama JP, Lobbi-Nivol C. Involvement of a mate chaperone (TorD) in the maturation pathway of molybdoenzyme TorA. J Biol Chem. 2003; 278:28787-92. https://doi.org/10.1074/jbc.M302730200 PMID: 12766163
56. Stewart V, MacGregor CH. Nitrate reductase in *Escherichia coli* K-12: involvement of *chlC*, *chlE*, and *chlG* loci. J Bacteriol. 1982; 151:788-99. PMID: 7047497

57. Saitou N, Nei M. The neighbor-joining method: a new method for reconstructing phylogenetic trees. Mol Biol Evol. 1987; 4(4):406-25. Epub 1987/07/01. https://doi.org/10.1093/oxfordjournals.molbev. a040454 PMID: 3447015.