Vicky Tzu-Jung
Juan†
a,
Patrick
Bajan†
b,
Chatchakorn
Eurtivong
c,
Tianyang
Liu
b,
Christopher J.
Squire
*de,
Evelyn Yu-Wen
Huang
*b and
Ivanhoe K. H.
Leung
*abe
aSchool of Chemical Sciences, The University of Auckland, Auckland, 1142, New Zealand
bSchool of Chemistry and Bio21 Molecular Science & Biotechnology Institute, The University of Melbourne, Parkville, Victoria 3052, Australia. E-mail: evelyn.huang@unimelb.edu.au; ivanhoe.leung@unimelb.edu.au
cDepartment of Pharmaceutical Chemistry, Faculty of Pharmacy, Mahidol University, Bangkok, 10400, Thailand
dSchool of Biological Sciences, The University of Auckland, Auckland, 1142, New Zealand. E-mail: c.squire@auckland.ac.nz
eMaurice Wilkins Centre for Molecular Biodiscovery, The University of Auckland, Auckland, 1142, New Zealand
First published on 24th October 2025
Rv3406 evolved from the ubiquitous taurine-catabolising enzyme TauD and functions as a sulfur-scavenging protein in Mycobacterium tuberculosis. Structural and biochemical analyses reveal specific changes that shape its chemical environment for ligand interaction and explain its broad substrate range. These findings show how amino acid substitutions redefine protein function and drive adaptation to the unique metabolic context of Mycobacteria.
In Mtb, Rv3406 functions as an alkyl sulfatase that hydroxylates alkyl sulfate esters, liberating inorganic sulfate for metabolic use (Fig. 1a).6 Rv3406 is the only alkyl sulfatase characterised in Mtb and is a member of the Fe(II)- and 2-oxoglutarate (2OG)-dependent oxygenase superfamily.7 Sequence alignment shows homology to other bacterial sulfur-scavenging oxygenases, the closest matches being to AtsK, an alkyl sulfatase from Pseudomonas putida (54% sequence identity),8 and TauD, a taurine dioxygenase from Escherichia coli (36% sequence identity).9 These homologues are functionally equivalent to Rv3406, and also allow P. putida and E. coli to scavenge sulfur from alternative, non-inorganic sulfate, sources.8,9
![]() | ||
| Fig. 1 Rv3406 is an Fe(II) and 2OG-dependent oxygenase in Mtb. (a) Rv3406 catalyses the hydroxylation of alkyl sulfate esters to release inorganic sulfate and an aldehyde. (b) NMR monitoring of Rv3406 activity in real time: signals for the substrates, 2OG and (R)-2EHS, and products, succinate and 2-ethylhexanal, can be observed directly. Spectra are at 8 min intervals, and the first spectrum was recorded 4 min after sample mixing. (c) 1H NMR time course allows the concurrent measurement of 2OG oxidation and substrate hydroxylation. The reaction mixture contained 2 μM Rv3406, 50 μM Fe(II), 1 mM ascorbate, 500 μM 2OG, 500 μM (R)-2EHS, 50 mM Tris-D11 (pH 8.0) in 10% D2O/90% H2O. We were unable to accurately integrate the resonances for 2-ethylhexanal due to signal overlap. | ||
Despite their homology and related biological functions, Rv3406, AtsK, and TauD, have different substrate specificity. While TauD accepts taurine, a small, non-proteinogenic amino acid,9 Rv3406 and AtsK prefer alkyl sulfates that contain large, hydrophobic sidechains, e.g., (R)-2-ethylhexyl sulfate ((R)-2EHS).6,8 The substrate preferences of AtsK and TauD likely reflect the sulfur sources available to their respective hosts.6,8,9 Understanding the structural determinants of substrate specificity in these closely related enzymes is key to clarifying the role of Rv3406 in Mtb sulfur metabolism and the nutrient conditions within macrophages. Here, we provide evidence that structural divergence of Rv3406 from its homologue TauD underlies its broader substrate specificity and hints at its refined role as a sulfur-scavenging enzyme in Mtb.
Previous biochemical studies of Rv3406 suggested broad substrate specificity but were limited to testing only medium-chain alkyl sulfates.6 In the standard 2OG oxygenase mechanism, oxidative decarboxylation of 2OG produces a reactive Fe(IV)
O intermediate that oxidizes the substrate.10,11 The coupling ratio and relative rates of 2OG and substrate oxidation are key indicators of substrate preference, as a weak substrate may not bind close enough to the metal centre or remain bound long enough for oxidation. To further explore the substrate preference of Rv3406, we produced recombinant Rv3406 and tested a wider range of alkyl sulfates (methyl, hexyl, octyl and dodecyl sulfate) (Table 1 and Fig. S1, SI). The reaction was characterised using a nuclear magnetic resonance (NMR)-based assay, which allows concurrent monitoring of substrate consumption and product formation (Fig. 1b and Fig. S2–S5, SI). This method is superior to other spectrophotometry-based assays as it allows the coupling ratio between substrate and co-substrate oxidation to be measured.
| Substrate | Activity | 2OG-to-substrate coupling | Relative rate to (R)-2EHS (%) |
|---|---|---|---|
| (R)-2EHS | Yes | 1 : 1 |
100 |
| Methyl sulfate | No | — | — |
| Hexyl sulfate | Yes | 90 : 1 |
0.3 |
| Octyl sulfate | Yes | 3.5 : 1 |
15.0 |
| Dodecyl sulfate | Yes | 1 : 1 |
31.7 |
| Methanesulfonate | No | — | — |
| Propane-1-sulfonate | No | — | — |
| Heptane-1-sulfonate | Yes | 8 : 1 |
3.8 |
| Octane-1-sulfonate | Yes | 5 : 1 |
4.5 |
| Dodecane-1-sulfonate | Yes | 1 : 1 |
39.5 |
| 4-Dodecylbenzenesulfonate | No | — | — |
| Taurine | No | — | — |
| Isethionate | No | — | — |
| 2-Mercaptoethane-1-sulfonate | No | — | — |
Mtb Rv3406 did not catalyse the turnover of methyl sulfate but displayed activity towards longer alkyl sulfate analogues. Among these, the highest activity was observed with (R)-2EHS (Table 1), for which 2OG oxidation was tightly coupled to hydroxylation of the substrate (Fig. 1c). Interestingly, although the alkyl moiety in hexyl sulfate has the same length as the principal chain in (R)-2EHS, without the additional ethyl branch, the relative turnover rate of hexyl sulfate was dramatically decreased, from 100% to ∼0.3% (based on 2OG turnover) (Table 1). The 2OG
:
substrate coupling ratio switched from 1
:
1 to 90
:
1 (Table 1), indicating that, for every 90 rounds of 2OG oxidation, only a single molecule of hexyl sulfate is hydroxylated. The trend in the sulfate series of substrates is for the relative turnover rate to increase and the 2OG
:
substrate coupling ratio to decrease, as the alkyl chain length increases (Table 1). Alkyl chain length and branching are clearly critical factors governing substrate turnover in Rv3406. Our results are also consistent with those reported by Sogi et al., who found that Rv3406 could catalyse the turnover of (R)-2EHS, hexyl sulfate, heptyl sulfate, and, to a smaller extent, pentyl sulfate.6
Dodecyl sulfate, the longest analogue tested, showed a coupling ratio of 1
:
1, the same as for (R)-2EHS, but the relative rate does not approach the 100% control, instead only reaching 31.7%. Interestingly, no aldehyde peak was detected in the reaction (Fig. S7, SI). This indicates that hydroxylation occurs at a different position, possibly away from the carbon directly bonded to the sulfate group. We propose that the longer alkyl chain of dodecyl sulfate forces the molecule to bend within the enzyme's binding pocket and shifts the site of hydroxylation. Frustratingly, we were also unable to observe a clear set of signals that correspond to the hydroxylated product in the NMR spectra. This could be due to the limited solubility and low yield of the product. Nonetheless, our findings clearly show that hydroxylation does not necessarily happen immediately adjacent to the ester oxygen on the substrate.
We then evaluated alkyl sulfonates which lack the bridging oxygen atom present in alkyl sulfates. Sulfonate substrates followed a similar trend to their sulfate analogues (Table 1). Sulfonates with less than seven carbons in their alkyl chain remain unmodified (Table 1). Heptane-1-sulfonate has slow turnover and an 8
:
1 coupling ratio; maximal turnover of the twelve carbon molecule affords only 39.5% conversion but a 1
:
1 coupling ratio (Table 1). The 4-dodecylbenzenesulfonate substrate shows that size is not the only determining factor. This largest sulfonate analogue, is not turned over and suggests the enzyme only recognises linear or simple branched chains, and not aromatic groups as substrate (we note that this chain is equivalent in length to a 15 carbon alkyl chain).
To define the structural determinants governing substrate binding to Rv3406 and turnover, we carried out extensive crystallography experiments. While crystals complexed with (R)-2EHS could not be obtained, we succeeded in crystallising Rv3406 using Ni(II) as a substitute for Fe(II) together with N-oxalylglycine (NOG) as an unreactive analogue/substitute of 2OG (Table S1 and S2, SI). The enzyme structure (Fig. 2a) displays the canonical structure displayed by earlier apo and metal bound structures (PDB: 4FFA and 4CVY), including tetrameric assembly and the conserved metal binding active site architecture (Fig. S8, SI). The Ni(II) is coordinated by the triad HXD/E…H motif (His97, Asp99, and His251) (Fig. 2b), conserved in all 2OG-dependent oxygenases,10 while NOG coordinates the Ni(II) in a bidentate manner (Fig. 2b) with the final metal coordination site opposite His251, occupied by a water molecule as a proxy for the substrate hydroxylation site (Fig. 2b). Comparison with the crystal structures of E. coli TauD (PDB: 1GY9)12 and P. putida AtsK (PDB: 1OIK),13 both protein–Fe(II)–2OG–substrate complexes, shows that the Ni(II) in our structure occupies the same position and has similar ligands as the Fe(II) in TauD and AtsK. NOG binds in a manner analogous to 2OG in these enzymes (Fig. S7 and S8, SI).
![]() | ||
| Fig. 2 Structure of Rv3406 with co-substrate mimic NOG. (a) Rv3406-Ni(II)-NOG structure (PDB: 8EVN). The molecular structure of NOG is depicted in the inset. The overall backbone and metal binding residues of our Rv3406-Ni(II)-NOG complex are the same as those of the published Rv3406-Fe(II) structure.6 (b) NOG binds to the active site metal ion in a bidentate manner. The final coordination site on Ni(II) opposite His251, is occupied by a water molecule as a proxy for the hydroxylation site for substrate oxidation. (c) The active site of TauD-Fe(II)-2OG with taurine binding (PDB: 1GY9) defines Asn95TauD and Tyr73TauD as critical residues stabilising the ligand. These two amino acids are substituted by alanine residues in Rv3406 and explain the inability of RV3406 to turnover taurine as substrate. Docking of (R)-2EHS to Rv3406-Ni(II)-NOG shows a binding model similar to that seen in the AtsK-Fe(II)-2OG-(R)-2EHS structure (PDB: 1OIK). | ||
Using our Rv3406–Ni(II)–NOG structure, we docked (R)-2EHS into the active site and found it positioned in the same location and orientation as in the AtsK–Fe(II)–2OG–(R)-2EHS structure (PDB: 1OIK). Our model places the sulfate of (R)-2EHS close to the water-coordination site of Ni(II) and held in place by two hydrogen bond interactions between sulfate oxygens and Val100 and Arg266 (Fig. 2c). The non-polar alkyl chain putatively extends into a non-polar pocket of the enzyme although a loop in this section of the X-ray crystal structure (residues 82–102) is disordered and not modelled. The location of disordered loops is similar to that observed in the AtsK-Fe(II)-2OG-(R)-2EHS structure (PDB: 1OIK) (S9, SI).13 In TauD (PDB: 1GY9), the equivalent sequence is an ordered β-strand structure (Fig. S10, SI).12 The flexibility within the disordered loop(s) may allow Rv3406 to accommodate alkyl chains of different length in the active site.
The amine group of taurine in the TauD-Fe(II)-2OG-taurine structure (PDB: 1GY9),12 hydrogen bonds with Tyr73TauD and Asn95TauD (Fig. 2c), residues that are substituted with alanine in Rv3406 (Fig S11, SI). These substitutions likely increase the size and hydrophobicity of the substrate binding pocket and could explain the observed preference for alkyl sulfates over taurine by Rv3406. To further understand the roles of Tyr73TauD and Asn95TauD in governing the substrate preference between TauD and Rv3406, we produced TauD Y73A, TauD N95A, and TauD Y73A/N95A mutants and compared them with WT TauD and Rv3406.
All three TauD mutants were highly active and were able to catalyse uncoupled 2OG turnover. Y73A/N95A TauD shows ∼2.5 times more 2OG turnover compared to the other mutants and both WT enzymes (Fig. 3a). The reason for this is not known. A slower turnover rate for uncoupled 2OG turnover may be due to 2OG binding in multiple orientations at the active-site metal, as there is no substrate to anchor 2OG in the “productive” orientation. A reduction in polar amino acids at the active site could prevent binding in non-productive orientations, since fewer hydrogen bond donors or acceptors would be available to stabilise 2OG in these positions. We then tested the mutants and WT enzymes with taurine as substrate (Fig. 3b). As expected, WT TauD turns over the most taurine, with equimolar amounts of succinate and aminoacetaldehyde produced at the reaction end point. WT Rv3406 is unable to utilise taurine as a substrate (Fig. 3b). All three TauD mutants accept taurine as a substrate; N95A TauD shows the highest activity, with an almost 1
:
1 2OG
:
substrate coupling ratio, whereas Y73A and Y73A/N95A TauD exhibit reduced activity compared to N95A (Fig. 3b). These results suggest that the tyrosine residue is important for taurine catalysis and likely acts synergistically with asparagine, as indicated by the cumulative loss of activity in the Y73A/N95A double mutant. In contrast, it is striking that with a more hydrophobic and spacious environment, TauD mutants were able to utilise (R)-2EHS as a substrate (Fig. 3c). This demonstrates the influence of the alanine substitutions in diversifying the sulfatase behaviour of the Rv3406 protein from that of its predecessor in taurine catabolism.
![]() | ||
Fig. 3 Mutagenesis reveals residues that govern substrate specificity. (a) Uncoupled turnover of 2OG for WT and mutant proteins expressed as a percentage of starting substrate concentration. The reaction mixture contained 2 μM enzyme, 50 μM Fe(II), 1 mM ascorbate, 500 μM 2OG, 5 μM catalase, 50 mM Tris-D11 (pH 8.0) in 10% D2O : 90% H2O. Reaction time was 1 hour. (b) Catalytic activity of TauD WT, Rv3406 WT and mutant TauD towards taurine. Coupling ratios are shown above the product formation plots. The reaction conditions are the same as above in the 2OG only reaction in (a), but with the addition of 500 μM taurine. (c) Catalytic activity of TauD WT, Rv3406 WT and mutant TauD protein towards (R)-2EHS. Coupling ratios are shown above the product formation plots. The reaction conditions are the same as above for 2OG but with the addition of 500 μM (R)-2EHS. | ||
We examined Rv3406's substrate preference and expanded our experimental range to include medium- and long-chain alkyl sulfates and sulfonates. Since many microorganisms use these compounds as alternative sulfur sources,14 our findings suggest that Mtb may similarly exploit them through the extended substrate profile of Rv3406. Through cross-species mutagenesis, we probed the structural basis for this promiscuous substrate preference and identified specific amino acid residues driving sulfatase binding and catalysis in Rv3406. Sequence similarity and structural analysis suggest that TauD is ancestral to Rv3406 (Fig. S12, SI). Comparing the two is intriguing because, while they perform similar biological roles – utilising alternative sulfur sources – they rely on entirely different organic substrates. The intracellular niche of Mtb within macrophages contrasts sharply with the commensal existence of E. coli in the gastrointestinal tract, exposing each organism to distinct nutritional landscapes, particularly in relation to sulfur availability. These environmental differences likely drive selection for enzymes capable of exploiting diverse substrates. Although the availability of different sulfate and sulfonate substrates inside host macrophage is unclear, it is possible that the broad substrate specificity of Rv3406 may represent an adaptive strategy to access varied sulfur sources within the macrophage. This warrants further investigation as little is known about the abundance, composition, and physiological relevance of these compounds. In addition to its role in sulfur scavenging, Rv3406 has been implicated in the oxidation and inactivation of several antitubercular agents.15–17 This work therefore provides valuable structure–activity relationship insights to guide drug development, either by avoiding susceptible functional groups or by designing inhibitors as potential adjuvants.
We acknowledge use of the Nuclear Magnetic Resonance Centre and Mass Spectrometry Hub at the University of Auckland, and the Melbourne Magnetic Resonance Platform, Melbourne Protein Characterisation Platform, and Mass Spectrometry and Proteomics Facility at the University of Melbourne for this work. We thank Dr Nicholas Demarais and Martin Middleditch of the University of Auckland for their technical input for mass spectrometry analyses. We thank Jingsong Zhou for his help with NMR experiments. This research was undertaken in part using the MX2 beamline at the Australian Synchrotron, part of ANSTO, and made use of the Australian Cancer Research Foundation (ACRF) detector. EYWH was supported by the Melbourne Research Scholarship, Rowden White Scholarship, Norma Hilda Scholarship, Dame Margaret Blackwood Soroptimist Scholarship, and Dr Albert Shimmins Postgraduate Writing-Up Award. IKHL acknowledge the University of Auckland and the University of Melbourne for financial support of this project.
Crystallographic data for Rv3406-Ni(II)-NOG has been deposited at the Protein Data Bank (PBD) under identification code 8EVN and can be obtained from https://doi.org/10.2210/pdb8EVN/pdb.
Footnote |
| † Vicky Tzu-Jung Juan and Patrick Bajan contributed equally to this study. |
| This journal is © The Royal Society of Chemistry 2025 |