Thank you for visiting nature.com. You are using a browser version with limited support for CSS. To obtain the best experience, we recommend you use a more up to date browser (or turn off compatibility mode in Internet Explorer). In the meantime, to ensure continued support, we are displaying the site without styles and JavaScript.
Scientific Reports volume 12, Article number: 14737 (2022 ) Cite this article 3 Aminopropylmethylamine
Deoxynivalenol (DON) is a mycotoxin, produced by filamentous fungi such as Fusarium graminearum, that causes significant yield losses of cereal grain crops worldwide. One of the most promising methods to detoxify this mycotoxin involves its enzymatic epimerization to 3-epi-DON. DepB plays a critical role in this process by reducing 3-keto-DON, an intermediate in the epimerization process, to 3-epi-DON. DepBRleg from Rhizobium leguminosarum is a member of the new aldo–keto reductase family, AKR18, and it has the unusual ability to utilize both NADH and NADPH as coenzymes, albeit with a 40-fold higher catalytic efficiency with NADPH compared to NADH. Structural analysis of DepBRleg revealed the putative roles of Lys-217, Arg-290, and Gln-294 in NADPH specificity. Replacement of these residues by site-specific mutagenesis to negatively charged amino acids compromised NADPH binding with minimal effects on NADH binding. The substrate-binding site of DepBRleg is larger than its closest structural homolog, AKR6A2, likely contributing to its ability to utilize a wide range of aldehydes and ketones, including the mycotoxin, patulin, as substrates. The structure of DepBRleg also suggests that 3-keto-DON can adopt two binding modes to facilitate 4-pro-R hydride transfer to either the re- or si-face of the C3 ketone providing a possible explanation for the enzyme’s ability to convert 3-keto-DON to 3-epi-DON and DON in diastereomeric ratios of 67.2% and 32.8% respectively.
Mycotoxins are toxic secondary metabolites produced by filamentous fungi that pose a serious health safety risk for humans and livestock1. Deoxynivalenol (DON), belongs to a class of Fusarium mycotoxins termed trichothecenes which share a tetracyclic sesquiterpene carbon skeleton2. Commercial cereal grain crops infected by toxigenic Fusarium species show diminished quality, and yield, and often contain high levels of DON. The FDA estimates crop losses due to mycotoxins in the US alone, to be around $932 million annually, fueling the need for more effective management strategies3.
DON is a potent protein synthesis inhibitor that blocks the A-site of ribosomes, hindering aminoacyl tRNA binding and preventing peptidyl transferase activity4. The ribosome interaction is mediated through hydrogen bond contacts from the C12, 13 epoxide ring, C3 hydroxyl, and van der Waal contacts between the C9, C10 double bond5. Together, these substituent groups contribute to DON’s overall toxicity. Certain bacteria, such as Devosia mutans 17-2-E-8 enzymatically detoxify DON by stereochemically inverting the C3 OH from the (S)-configuration to the (R)-configuration producing the diastereomer, 3-epi-DON6. This impairs binding to the ribosome and overall, reduces its toxicity; as evidenced by MTT bioassays, 3-epi-DON possesses an IC50 357 times higher than DON7. This process, termed ‘DON epimerization’ (Dep) occurs in two stages8. First DepA, a PQQ-dependent alcohol dehydrogenase oxidizes the C3 OH group of DON to form 3-keto-DON9. A second enzyme, DepB, then reduces 3-keto-DON at the re-face of this stereogenic center to produce 3-epi-DON8,10 (Fig. 1). The intermediate, 3-keto-DON is only moderately less toxic than DON with an IC50 3-times higher than DON7. Therefore, the critical step of DON detoxification is catalyzed by DepB.
The Dep pathway. The C3 OH group of DON (shown in red) is stereochemically inverted to produce the diastereomer, 3-epi-DON. In the expanded image of this stereogenic center, R1 and R2 represent the priority groups attached to the sp2 hybridized C3. Hydride attack on the re-face at the prochiral C3 center generates 3-epi-DON.
DepB belongs to the aldo–keto reductase (AKR) family of enzymes. These ubiquitous oxidoreductases occur in all domains of life, and they catalyze the reduction of diverse carbonyl substrates. They play vital roles in steroid hormone metabolism11, oxidative stress responses12,13, xenobiotic detoxification14, and are implicated in certain cancers15 and metabolic disorders such as diabetes16,17. Most AKRs are monomeric proteins with molecular weights ranging from 34 to 37 kDa and predominantly utilize NADP(H) as their coenzyme.
Here we present the crystal structure of DepB from Rhizobium leguminosarum (herein designated as DepBRleg). We determined that DepBRleg belongs to a recently discovered family of AKRs and although it can use either NADH or NADPH as a coenzyme, the structural basis for its strong preference for NADPH is unknown. Although the stereochemistry of hydride transfer is well established for AKRs, the stereoselectivity of the reduction process, which is key to detoxification, is a less explored avenue. We, therefore, probed the determinants of DepBRleg’s unusual coenzyme specificity and propose features of the substrate-binding site which dictate the stereochemical outcome of 3-keto-DON reduction through molecular modeling studies. In addition to its role in the Dep pathway, DepBRleg is also active towards a wide range of aldehyde and ketone substrates including the mycotoxin patulin, making it a versatile biocatalyst for the decontamination of various xenobiotics.
AKRs are traditionally classified into families based on sequence identity and phylogenetic analysis. The nomenclature system consists of the root symbol AKR, followed by an arabic numeral designating the family and an alphabet representing the subfamily18. The number of AKRs in the protein database continues to expand and to date 18 families have been identified. Protein sequences in each family share at least 40% sequence identity while sub-families share 60% sequence identity19. We used a protein sequence similarity network (SSN)20 as an alternate, less computationally demanding method to determine sequence relationships of DepBRleg with other AKR family members. For ease of computation, sequences were iteratively clustered to produce a 40% representative node network and filtered based on sequence length. The total number of sequences following the reduction of the initial dataset is 878. In this network, nodes represent protein sequences sharing 40% or more sequence identity and edges represent the pairwise alignments between protein sequences. Edges were drawn between nodes if they exceeded the prescribed stringency threshold or BLAST E-value of e−57.
Previous phylogenetic analysis showed that AKRs can be divided into two large groups with AKR1-AKR5 forming one branch of the phylogenetic tree and the other AKRs forming a separate branch21. This is generally reflected in the SSN topology where AKR1 to AKR5 family members cluster together (Fig. 2A). At this stringency threshold, there was good delineation of AKR7-11 and AKR 13-15 into distinct, isofunctional clusters. DepBRleg and its homologs from D. mutans 17-2-E-810 and AKR18A1 from Sphingomonas S3-422 clustered with AKR6, AKR12, and AKR14 families. At higher stringency thresholds (e−67), DepBRleg and its closest homologs resolve into a single isofunctional cluster (Fig. 2B). Curiously, in vitro assays previously revealed that AKR18A1 reduces 3-keto-DON to DON rather than 3-epi-DON in the presence of NADH22. The enzyme also reduces the C7 ketone group of the estrogenic mycotoxin, zearalenone (ZEN) to produce α-zearalenol (α-ZOL) and β-zearalenol (β-ZOL). Aside from this, AKR18A1’s substrate specificity towards other endogenous aldehydes and ketones has not been examined, nor has its crystal structure been solved.
A 40% representative sequence similarity network (SSN) at a threshold of e−57. (A) Nodes within each cluster contain a representative protein sequence of a collection of sequences that share 40% or more sequence identity. Clusters have been annotated and color-coded based on curated protein sequences from the AKR Database21. Grey nodes represent AKRs that have no functional annotation. Metabolic pathways associated with the function of AKR enzymes have been color-coded to cross-reference with each cluster and numerically labeled to designate each AKR family number. The nodes of AKRs possessing biochemical data on coenzyme specificity and crystal structures have also been accordingly labeled as outlined in the legend. (B) SSN generated at a threshold of e−67 depicting complete segregation of the AKR18 family from AKR6, AKR12 and AKR14 members. (C) Genome Neighborhood Diagram of DepBRleg. The GND was generated using the EFI-GNT server which depicts the coding sequence region of depBRleg along with putative functions of upstream and downstream genes.
Besides AKR18A1, the next closest homologs of DepBRleg identified from this SSN analysis were AKR12 family members based on the first stringency cut-off threshold (e−57). As per AKR classification rules, DepBRleg is not an AKR12 family member as it shares less than 40% sequence identity with enzymes from this family. In addition, these AKR12 members are involved in the biosynthetic pathways of polyketide macrolide antibiotics such as tylosin, erythromycin, and avermectin23–25. A genome neighborhood diagram (GND) of DepBRleg revealed the presence of a putative monoamine oxidase (MAO) upstream of the DepBRleg gene along with several putative spermidine/putrescine transporter genes downstream of DepBRleg (Fig. 2C). The neighboring genes reported for the depB homolog in D. mutans 17-2-E-8 differ from depBRleg but are also not polyketide synthesis genes10.
Recombinant N-terminal His-tagged DepBRleg was purified to homogeneity by Ni-NTA chromatography. A band corresponding to a molecular weight of approximately 39 kDa was observed on a 10% SDS-PAGE gel (See Supplementary Fig. S1). The substrate specificity of the purified enzyme towards 3-keto-DON, aliphatic and aromatic carbonyl compounds, as well as nicotinamide cofactors were assessed by steady-state kinetics (Table 1). DepBRleg possesses the highest specificity constant (kcat/Km) with the diketone 9,10-phenanthrenequinone (9,10-PQ) on the order of 27 times higher relative to 3-keto-DON and 35 times higher relative to the smaller diketone, isatin. DepBRleg is also active towards endogenous toxic aldehydes derived from oxidative stress responses such as lipid peroxidation26,27. These aldehydes include acrolein, methylglyoxal, 4-oxo-2-nonenal, 4-hydroxy-2-nonenal, hexanal, and butanal. Among the aldehydes, DepBRleg displayed high specificity constants for the smaller unsaturated aldehydes such as acrolein and methylglyoxal, although the Km for methylglyoxal was too high to be determined reliably. Among the αβ unsaturated aliphatic aldehydes, 60-fold higher specificity constant was observed for 4-oxo-2-nonenal compared to 4-hydroxy-2-nonenal. The apparent Km for NADPH is about 100-fold lower than NADH (15.2 ± 1.41 µM vs. 1560 ± 399 µM) when evaluated using DL-glyceraldehyde as a substrate at a fixed concentration of 10 mM, while the apparent kcat with NADPH was only 13-fold higher than with NADH (0.337 ± 0.00757 s−1 vs. 0.0242 ± 0.00344 s−1). DepBRleg was also determined to reduce the mycotoxin patulin but not citrinin. Both these mycotoxins are produced by Penicillium expansum which is the causative agent of blue mold rot in apples28. LC–MS/MS analysis indicated that patulin was indeed transformed to E-ascladiol, a by-product previously reported to be significantly less cytotoxic compared to patulin29,30. The mechanistic basis for this transformation has been proposed to involve the spontaneous opening of the hemiacetal ring of patulin followed by the reduction of the aldehyde to the alcohol31 (See Fig. S2 for LC–MS/MS results).
DepBRleg was successfully crystallized, and its structure was solved by molecular replacement using an AKR from Polaromonas sp. JS666 as a search model (PDB ID: 4XK2). Structural and refinement statistics are summarized in Table 2.
Four protomers of DepBRleg are present in the asymmetric unit designated as chains A, B, C, and D. Overall, the three-dimensional structure of each protomer adopts the classic triose-phosphate isomerase (TIM) barrel fold (α/β)8, except that the presence of Pro-204 in the region that would otherwise form strand β7 disrupts the hydrogen bonding pattern leaving this region as an extended loop. Two auxiliary helices designated as H1 and H2 are found at the periphery of the barrel, while the N-terminal end of the barrel is capped by an N-terminal β-hairpin loop. The TIM barrel scaffold contains βα loops that connect β-strands to α-helices and the αβ loops that connect α-helices to β-strands. βα loops are longer than αβ loops and impart structural and hence functional diversity to TIM barrel proteins such as AKRs32. The β3α3 loop (designated as Loop A) and the C-terminal loop (designated as Loop C) are proposed to govern substrate binding in AKRs while the primary role of the longer β7H1 loop (designated as Loop B) is for coenzyme binding33 (Fig. 3A). Due to a lack of clear electron density, the following residues within loop B could not be modeled: Chain A (214–241), Chain B (217–244), Chain C (219–245), and Chain D (218–245).
Crystal structure of DepBRleg protomer and octameric assembly. (A) Cartoon representation of the protomer with a view looking down the top of the TIM barrel rendered in black and white with the β7 annotated in red. Loop A (orange) caps the top of the substrate-binding cleft, loop B (dashes) could not be modeled but its position has been highlighted in this figure, and finally, loop C (green) which extends into the active site of the opposing protomer in the dimer. A sequence logo highlighting the conservation of the catalytic residues has also been provided. (B) Surface and cartoon rendition of the DepBRleg octamer. This assembly was predicted using the PISA server. View rotated at a 90° angle depicting 3D domain swapping of loop C at the dimer interface.
A pair of protomers (A and C, B and D) associate together forming a dimer with an average interface area of 2362.9 Å2, stabilized by 24 hydrogen bonds and 13 salt bridges. Dimerization involves 3D domain swapping with the exchange of the α-helix of loop C which interacts with loop A on the associating protomer. These dimers interact around the crystallographic four-fold symmetry axis to form a D4 symmetric octamer (Fig. 3B), with key interactions mediated by the N-terminal β-hairpin, α3, α4, and the β4α4 loop. PISA34 estimates that these interactions bury 748.6 Å2 and result in a solvation-free energy gain of -809.8 kcal/mol. The absolute molecular weight of DepBRleg was determined experimentally by SEC-MALS to be 325.6 ± 12.1 kDa which is consistent with the 317.5 kDa predicted for a DepBRleg octamer. (See Supplementary Fig. S3).
The closest structural homologs of DepBRleg, are the functionally uncharacterized AKRs from the structural genomics initiative, Polaromonas sp. AKR (PDB ID: 4XK2, 44% sequence identity, 1.6 Å RMSD), an Escherichia coli AKR, Tas (PDB ID: 1LQA, 35% sequence identity, 2.0 Å RMSD), that has broad substrate specificity35, the beta subunit of the voltage-gated potassium channel from Rattus norvegicus (AKR6A2; PDB ID: 3EAU, 32% sequence identity, 2.1 Å RMSD), E. coli AKR14A1 (PDB ID: 4AUB, 35% sequence identity, 2.4 Å RMSD) that is specific towards methylglyoxal and diketones such as isatin36 and mithramycin side chain reductase from Streptomyces argillaceus (PDB ID: 6OW0, 33% sequence identity, 2.3 Å RMSD). Mithramycin side chain reductase is a ketoreductase that reduces the 4’ ketone side chain of mithramycin DK, an intermediate in the biosynthesis of the tricyclic antitumor polyketide, mithramycin37.
Diversity in certain loop regions was observed among these structural homologs. Differences in the conformation of loop B and the auxiliary helices, H1 and H2 could be due to conformational changes upon binding of NADP(H) and/or substrates. Loop B in the apo-structures of DepBRleg and Polaromonas sp. AKR is disordered, but clear electron densities for this region were observed for the NADP+ binary complexes of the other structural homologs. This is in agreement with the finding that loop B becomes ordered upon coenzyme binding38. In addition to these conformational differences, variations in loop length are also evident, particularly for the E. coli Tas protein which possesses a long α4β4 loop (24 residues) compared with DepBRleg (8 residues). The length and structure of DepBRleg loop C (26 residues) also differed with AKR6A2 (6 residues), Tas (5 residues), AKR14A1 (10 residues), and mithramycin side chain reductase (6 residues). In the hydroxysteroid reductase subfamily of AKRs truncation of this loop led to reduced specific activity towards steroid substrates, therefore the diversity in loop lengths among these AKRs likely contributes to corresponding differences in substrate specificity39.
DepBRleg reduced 3-keto-DON using NADH and NADPH as coenzymes, although the catalytic efficiency with NADH was 40-times lower than with NADPH. The dissociation constant (Kd) for NADH, as determined by tryptophan fluorescence quenching experiments, was about sixfold higher than NADPH. In contrast, most AKRs are NADPH specific with NADPH dissociation constants on average, 1000-fold lower than NADH40. DepBRleg, therefore, shares the rare dual coenzyme specificity with a small number of enzymes from different AKR families, including the hydroxysteroid dehydrogenases (AKR1C)41,42, xylose reductases (AKR2B)43,44,45, beta subunit of voltage-gated potassium ion channels (AKR6A)46, aflatoxin reductases (AKR7A)47, and pyridoxal dehydrogenases (AKR15A)48. The structural basis for NADH specificity was studied for Candida tenuis xylose reductase (AKR2B5), revealing a key Glu-227 residue on loop B which forms bidentate hydrogen bonds with the 2’ and 3’ hydroxyl groups of the adenosine moiety of NAD+44. In DepBRleg, this residue is also conserved and corresponds to Glu-222. In the NADP+ bound complex of AKR2B5, the 2’ monophosphate of NADP+ interacts with a positively charged Arg-280, the peptide backbone of Asn-276, and Ser-275. These residues are not strictly conserved in DepBRleg and are instead replaced with Gln-294, Arg-290, and Ala-288, underscoring the general variability in NADP+ interactions among the AKRs.
Although NADP+ was present at 400 µM during the crystallization of DepBRleg, there was no electron density corresponding to this coenzyme. The putative coenzyme binding pocket was suitably defined and superimposed well with coenzyme binding sites of AKR binary complexes. AKRs bind NADPH in an extended anti-conformation to achieve the stereospecificity of hydride transfer33. An aromatic residue typically forms π-stacking interactions with the nicotinamide ring while the carboxamide moiety is oriented by highly conserved glutamine, serine, and asparagine residues. In DepBRleg, these residues correspond to Trp-206, Gln-178, Ser-152, and Asn-153 respectively. In AKR6A2, the closest structural homolog of DepBRleg, the 2’monophosphate of NADP+ interacts with the side chains of Gln-62, Lys-254, Ser-325, and Gln-329, and (Fig. 4A). In DepBRleg, these residues correspond to Glu-26, Lys-217, Arg-290, and Gln-294 respectively. To examine the conservation of these residues, a multiple sequence alignment was constructed with AKR family representatives. Lys-217 which is present on loop B is not strictly conserved and is often substituted for a basic residue or a smaller hydrophobic residue. Arg-290 is partially conserved across the AKR families and this position is often frequented by a basic residue or a polar hydrophilic residue. Gln-294 is strictly conserved across AKR6A2, AKR7A1, AKR9A1, AKR11A1, AKR14A1, AKR15A1, and AKR18A1 but for the other families, it is replaced with either a basic or aromatic residue (Fig. 4B). Lastly, Glu-26 was the least conserved and was therefore not selected for further analysis.
Model of the coenzyme binding site highlighting the conservation of putative coenzyme binding residues. (A) NADP+ was modeled by superimposing the ternary complex of AKR6A2 with apo-DepBRleg. White sticks represent the residues from DepBRleg complex, while cyan sticks are residues corresponding to AKR6A2. Trp-206 has been re-oriented in the apo-DepBRleg structure to prevent a steric clash with the nicotinamide ring of NADP+. π-stacking interactions between the nicotinamide ring and Trp-206 are shaded in blue. In apo-DepBRleg, Lys-217 points away from the entrance of the coenzyme binding pocket. However, by analogy to Lys-254 in AKR6A2, this residue may interact with at least one of the oxygen atoms of 2’monophosphate. (B) MSA depicting conservation of residues near the 2’monophosphate of the modeled NADP+ (See Supplementary Table S1 for accession codes).
All recombinant N-terminal His-tagged DepBRleg coenzyme variants were purified by Ni-NTA chromatography with expected molecular weights of 39 kDa (See Supplementary Fig. S1). Replacements of Lys-217, Arg-290, and Gln-294 for negatively charged glutamate significantly altered the Kd for NADPH relative to wild type DepBRleg, with no significant effect on the Kd of NADH. Among these variants, the R290E variant displayed an 11-fold increase, followed closely by Q294E with a ninefold increase and finally K217E with only a sixfold change (Table 3). Substitutions for neutral residues as in the case for R290N displayed a moderate change in the Kd by fivefold, while for K217M, the Kd change relative to the wild type was more pronounced with a nearly sevenfold change. Overall, for all the mutants, there was no appreciable change in Kd for NADH.
The catalytic mechanism of AKRs involves a stereospecific 4-pro-R hydride transfer from the nicotinamide ring of the coenzyme to the substrate carbonyl carbon, followed by protonation of the carbonyl oxygen by a tyrosyl residue49. This is facilitated by a proton relay with histidine or a lysine-aspartate pair40. In DepBRleg, these catalytic residues are conserved: Asp-48, Tyr-53, Lys-81, and His-122. A water molecule occupies the position equivalent to the carbonyl oxygen group of the substrate and forms hydrogen bonds (2.6–2.9 Å) to the side chains of both His-122 and Tyr-5349. DepBRleg’s substrate-binding pocket is lined with the following residues: Met-21, Asp-48, Val-52, Tyr-53, Lys-81, Arg-83, Phe-84, His-122, Ala-123, Ser-152, Asn-153 and Gln-180. Overall, DepBRleg’s substrate-binding pocket is slightly larger and more hydrophobic in comparison with the NADP+- cortisone bound complex of AKR6A2. The physiological substrate of AKR6A2 is unknown, however, much like DepBRleg it also displays broad substrate specificity with aldehyde and ketone substrates50.
The stereospecificity of hydride transfer is strictly conserved in AKRs. In DepBRleg, the residues Ser-152, Asn-153, and Gln-178 are critical for maintaining the anti-conformation of the nicotinamide ring to achieve this stereospecificity51. However, hydride attack may occur either on the re-face or si-face of the prochiral carbonyl group. We determined by HPLC analysis that DepBRleg reduced 3-keto-DON to produce a diastereomeric ratio of 67.2% for 3-epi-DON and 32.8% for DON (See Supplementary Fig. S4 for HPLC chromatograms). To examine the molecular basis for this diastereoselectivity, we superimposed apo-DepBRleg with the recently solved Debaryomyces nepalensis xylose reductase (PDB ID: 5ZCM) complexed with a DTT-NADP+ adduct52. In this complex, the geometry of the C4N of NADP+ strongly resembles the puckered ring conformation of reduced NADPH. 3-keto-DON was then modeled into this complex in an orientation poised for hydride transfer to the re-face to produce 3-epi-DON. The carbonyl oxygen of 3-keto-DON was placed within hydrogen-bonding distance of Tyr-53 and His-122, while Arg-343 from Loop C was re-positioned to provide hydrogen bond contacts with the C8 ketone of 3-keto-DON. Phe-84, which is present on Loop A could potentially stack with the cyclohexene ring of 3-keto-DON (Fig. 5A). This residue corresponds to Trp-86 in Thermotoga maritima AKR, an enzyme that catalyzes the reduction of ethyl 2-oxo-4-phenylbutyrate to the R- and S-enantiomer of ethyl-2-hydroxy-4-phenylbutyrate53. In that study, the replacement of Trp-86 with smaller amino acids alleviates space constraints and enabled the enzyme to increase the production of the R-alcohol53. To produce DON, 3-keto-DON would require a 180° flip in the substrate-binding pocket (Fig. 5B). An observation gleaned from sequence analysis with other AKRs, is the presence of a bulky residue following the catalytic histidine (Fig. 5C). However, in DepBRleg this position is occupied by Ala-123 which eliminates any space constraints allowing for the accommodation of the C15 primary alcohol of 3-keto-DON. We note that in AKR18A1 which reduces 3-keto-DON to DON, a glycine residue occupies the positions corresponding to Phe-84 and Ala-123 providing a great deal of flexibility in its respective active site.
Model of the ternary complex of DepBRleg. (A) Binding orientation for 3-keto-DON for hydride attack at the re-face. π-stacking interactions between Phe-84 and the cyclohexene ring are shaded yellow. (B) Binding orientation for 3-keto-DON for hydride attack at the si-face. (C) Multiple sequence alignments depicting the sequence variability at positions corresponding to Phe-84 and Ala-123 in DepBRleg with that of various AKR family representatives.
A comprehensive biochemical and structural analysis was conducted for DepBRleg, a close homolog of DepB from D. mutans 17-2-E-8. This AKR catalyzes the critical step during DON detoxification to afford the significantly less toxic diastereomer, 3-epi-DON. SSN cluster analysis revealed that DepBRleg is a member of the recently discovered AKR18 family. Within this family, only preliminary biochemical characterization was conducted for AKR18A1 from Sphingomonas S3-4 suggesting its involvement in a DON detoxification pathway22. At a cutoff threshold of e−57, DepBRleg localizes to a functionally divergent cluster in the SSN comprised of AKR6, AKR12, and AKR14 families. AKR6 members comprise the cytosolic subunit of voltage-gated potassium ion channels and are proposed to modulate these channels via coenzyme binding46. AKR12 members are reductases involved in macrolide antibiotic synthesis23,24,25, while the AKR14A1 member identified in the SSN is proposed to be involved in methylglyoxal detoxification54. The gene organization as depicted in the GND indicates that depBRleg is not located in a similar biosynthetic operon. Instead, it is flanked by a putative MAO and several putative putrescine/spermidine transporter protein genes.
DepBRleg’s substrate specificity was examined across a diverse range of model carbonyl substrates with its preferred coenzyme, NADPH. Amongst the ketone substrates, DepBRleg possesses the highest specificity constant for 9,10-PQ while comparable specificity constants were observed for the other diketone, isatin and 3-keto-DON. 9,10-PQ is a xenobiotic air pollutant that triggers the generation of reactive oxygen species (ROS) in vivo12 while isatin is an endogenous indole derived from tryptophan or phenylalanine metabolism by gut bacteria and has cytotoxic effects on microbes and various human cancer cell lines55. Overall, higher specificity constants were observed for ketones, versus aldehyde substrates. Among the lipid peroxidation substrates, DepBRleg displayed the highest specificity constant for 4-oxo-2-nonenal, followed by methylglyoxal and acrolein. DepBRleg had low specificity constants for 4-hydroxy-nonenal and DL-glyceraldehyde. The latter is typically a model substrate used for aldose reductases and aldehyde reductases from AKR1 families17. Surprisingly, DepBRleg can also reduce the mycotoxin patulin to E-ascladiol, an activity that was previously described only for the short-chain dehydrogenase family of oxidoreductases56,57. In summary, the broad substrate specificity of this AKR, therefore, lends credence to its secondary role in general detoxification processes.
A consensus among most AKRs is their preference for the coenzyme, NADPH over NADH49. Investigations of the coenzyme binding affinities of DepBRleg revealed that it can utilize both NADPH and NADH, albeit with preference for NADPH over NADH. Given that NADPH is more expensive than NADH, the ability of DepBRleg to utilize NADH as a coenzyme will be beneficial for its use as a biocatalyst to detoxify DON. Based on comparative structural analysis, the residues Lys-217, Arg-290, and Gln-294 are proposed as key contributing elements for DepBRleg’s preference for NADPH. These structural elements likely contribute to a hydrogen bond network encompassing the 2’monophosphate of NADPH. Precedents for coenzyme specificity are well established for short-chain dehydrogenases (SDRs). In NADP+ dependent SDRs, a basic residue in the glycine-rich motif of the Rossman fold lies governs NADP+ specificity while for NAD+ dependent SDRs, typically an aspartate or glutamate is present58. Similarly, in DepBRleg, the basic residues, Arg-290 and Lys-217 contribute the necessary salt-bridge interactions with the 2’monophosphate, while Gln-294 also provides hydrogen bond interactions and provides planar stacking interactions between the amide group of its side chain and the adenosine ring. Not surprisingly, substitutions for negatively charged residues generated electrostatic repulsion with the 2’monophosphate group resulting in an increased NADPH Kd value for R290E, Q294E, and K217E. While the change in Kd for the R290N variant was not as dramatic, we rationalize that replacement of the arginine for a polar, shorter residue like asparagine likely altered the hydrogen bond network coordinating the 2’monophosphate which consequently affected the binding mode of NADPH. K217M and K217E displayed a similar increase in Kd for NADPH, however, the effect was not as pronounced as in the Arg-290 variants. As Lys-217 is present on Loop B, which displays a great deal of sequence variance, it is plausible that a proximal, basic residue may have compensated for the absence of the lysine residue in the K217M and K217E variants.
Reduction of 3-keto-DON by DepBRleg produced 3-epi-DON and DON in diastereomeric ratios of 67.2% and 32.8% respectively. Modeling suggests that the formation of the two diastereomers is due to alternate 3-keto-DON binding modes that enable the re or si-face to approach the 4-pro-R hydrogen of the nicotinamide cofactor. Although DepBRleg does not show strict diastereoselectivity towards 3-epi-DON formation, coupling the enzyme with DepA which has an apparent Km of 32 ± 4 µM for DON9, would drive the equilibrium towards 3-epi-DON formation.
In conclusion, this comprehensive analysis of DepBRleg provides extensive insights into the AKR18 family. As evidenced by this structure–function study, DepBRleg presents as a versatile biocatalyst owing to its broad substrate specificity and a large active site which can accommodate a range of toxic endogenous aldehydes and ketones. The availability of the crystal structure of the enzyme should also facilitate future protein engineering work in targeting residues in the vicinity of the proposed substrate-binding site to improve the diastereoselectivity of the enzyme. While the feasibility of employing DepB in DON detoxification applications remains a contentious issue due to its NADPH dependence, this study provides a translational framework to modify its coenzyme specificity.
3-keto-DON, DON, and E-ascladiol were purchased from TripleBond (ON, Canada). Patulin, citrinin, 4-oxo-2-nonenal, and 4-hydroxynonenal were purchased from Cayman Chemicals (MI, USA). 3-epi-DON was previously purified and the identity was confirmed via LC-MS6. All other chemicals unless otherwise stated were from Sigma-Aldrich and Thermofisher.
His-tagged DepBRleg and variants were recombinantly expressed in E.coli BL21 LOBSTR (low background strain). This engineered strain minimizes background contamination by histidine-rich E. coli proteins such as SlyD and ArnA during Ni-NTA chromatography59. An overnight starter culture was used to inoculate 4L of LB media and cultures were grown at 37 °C with shaking speed at 200 RPM. Recombinant protein expression was induced with 1 mM isopropyl β-D-1-thiogalactopyranoside (IPTG) and incubated at 15 °C overnight with shaking. Cells were harvested by centrifugation and washed with 20 mM HEPES pH 8.0. The pellet was resuspended in 20 mM HEPES pH 8.0 buffer containing up to 1 mg/mL DNase I and lysed by 7–8 passages through a French press at 15,000 lb/in2. Cell debris was removed through centrifugation at 4 °C and the clarified lysate filtered through a 0.45 µm filter before incubation for 1 h at 4 °C with Ni2+-NTA resin in buffer (50 mM sodium phosphate buffer pH 8.0, 300 mM NaCl) with binding buffer (20 mM imidazole pH 8, 150 mM NaCl). The mixture was loaded onto a gravity column and washed with the same binding buffer. His-tagged proteins were eluted with 150 mM imidazole pH 8. Buffer exchange with 20 mM HEPES pH 7.5 with 10% glycerol and 150 mM NaCl was conducted in a stirred cell equipped with a YM10 filter (Amicon). Protein concentration was determined using a Bradford Assay60 with Bovine Serum Albumin used as a standard. The purity of the recombinant enzyme and the molecular weight of recombinant DepBRleg and coenzyme variants was estimated using Coomassie blue-stained SDS-PAGE.
The absolute molecular weight of DepBRleg was determined using Size Exclusion Chromatography coupled with Multi-Angle Light Scattering (SEC-MALS). The column utilized for size exclusion was a P3000 single-pore GPC/SEC column (Malvern Panalytical). The system was equilibrated with the following buffer system: 20 mM HEPES, 150 mM NaCl, 0.02% NaN3 at pH 7. The concentration of the protein samples was 0.5 mg/mL and samples were filtered through 0.2 µm filters prior to injection.
Crystallization conditions for DepBRleg were screened using the JCSG-plus™ kit (Molecular Dimensions) and an initial hit was obtained in condition B1 (0.1 M citrate pH 4, 0.8 M ammonium sulfate). Condition B1 was optimized and a crystal was obtained in the following optimized condition: 0.8 M ammonium sulfate, 20 mM ammonium acetate, and 0.1 M TRIS citrate. The protein was crystallized at 11 mg/mL in the presence of 400 µM NADP+ via sitting drop vapor diffusion at 288 K. X-ray diffraction data were collected at the Canadian Macromolecular Crystallography Facility (CMCF-BM), processed using XDS, and scaled using XSCALE61. Crystals diffracted to 2.5 Å, and belonged to the space group P4212, with unit-cell parameters a = b = 159.45, c = 163.9 Å. The crystal structure of DepBRleg was solved using the automated molecular replacement pipeline, MrBUMP which generates an exhaustive list of search models. First, multiple sequence alignments are conducted between the query and homologous protein structures, second, the top hits are subjected to model preparation using a suite of CCP4 programs and finally, the models are submitted to Phaser for molecular replacement62. The best search model prepared from this process is the AKR from Polaromonas sp. JS666. AutoBuild and subsequent refinements were conducted in PHENIX63 and final model building was conducted in COOT64.
Structural alignments were conducted via the DALI server65. The ternary complex of the AKR6A2 member, Rattus norvegicus voltage-dependent K+ channel beta subunit (Kvβ) (PDB ID: 3EAU) was superimposed on to the apo-DepBRleg structure. Structures were visualized in PyMOL (Schrödinger, Inc.) and torsion angles of residues were adjusted in the DepBRleg structure to prevent steric clash with the 2’ monophosphate and the nicotinamide ring of NADP+. Putative residues involved in NADPH specificity were selected within a 4–5 Å distance from the 2’ monophosphate attached to the adenosine ribose. Multiple sequence alignments were constructed in UGENE66 using representative sequences from each AKR family and aligned using the MUSCLE algorithm67 Based on this analysis, the following residues were selected as candidates for mutagenesis: Lys-217, Arg-290 and Gln-294.
Site-specific mutagenesis and mutagenic primer design were conducted as per the Single-Primer Reactions IN Parallel (SPRINP) protocol previously described68. The wild type DepBRleg plasmid construct (pET28aDepBRleg) was used as a template with single point mutation oligonucleotide primers (Integrated DNA Technologies) and Q5 High-Fidelity DNA Polymerase (New England Biolabs). All variant constructs were confirmed by a double restriction digest (NdeI/HindIII) followed by sequence analysis conducted at (AFL) Laboratory Services Division at the University of Guelph. All primers used for mutagenesis are listed in Supplementary Table S2.
The AKR SSN was constructed as previously described69. The master list of AKR sequences was downloaded from UniProt using the PFAM identifier (PF00248). In total, 318,574 sequences were retrieved, and the dataset was reduced by iteratively clustering the sequences using the H-CD-HIT server70. Sequences were then filtered based on an amino acid length cutoff between 320–350 to remove fragments and AKRs part of large multi-domain complexes. Finally, curated sequences from the AKR Database (https://hosting.med.upenn.edu/akr/) and the PDB were also added to the list of sequences, and any duplicates were removed by submission to CD-HIT with the sequence identity cut-off parameter set at 100%. An in-house BLAST search against this finalized database was then conducted at varying E-values and the resulting network was visualized using Cytoscape71. An optimal stringency threshold was selected to minimize the extensive fragmentation of the network.
The thermal stability of wild-type DepBRleg and variants were assessed using SYPRO Orange while fluorescence was monitored with a StepOnePlus Real-Time PCR system (Applied Biosystems, Foster City, U.S.A.). Experiments were carried out in triplicate in 100 mM HEPES pH 7.5 with a final protein concentration of 0.3 mg/mL and 1X stock of 5000X Sypro Orange as previously described72 (See Supplementary Table 3).
Dissociation constants (Kd) of DepBRleg and each coenzyme was determined using a tryptophan fluorescence-quenching assay with the PTI Fluorimeter (with FelixGX software). Excitation and emission wavelengths of 295 nm and 335 nm respectively were used while slit widths were fixed at 4 nm. Titrations were performed in 20 mM HEPES pH 7.5 with 120 μg of protein. The absorbance of only protein and buffer was measured before the assay to minimize any inner filter effects. All assays were conducted in triplicate. The change in fluorescence (ΔF) relative to the control (protein and buffer) was determined and then plotted against the concentration of each respective coenzyme. The data were fitted to the following equation73 by nonlinear regression using GraphPad Prism Version 8 software to determine the Kd:
F is the protein fluorescence intensity at varying concentrations of either NADPH or NADH represented by [L], and F0 is the fluorescence intensity of the protein in the absence of coenzyme.
For substrate specificity experiments, reactions were carried out in 1 mL reaction volumes consisting of 100 mM HEPES pH 7.5 buffer with the NADPH concentration fixed at 1 mM. The kinetic parameters for NADPH and NADH were determined with DL-glyceraldehyde concentration fixed at 10 mM. Substrate utilization was monitored by the decrease in absorbance of NADPH or NADH at 380 nm (ɛ = 1132 M−1 cm−1) using the Varian Cary 100 spectrophotometer at 298 K. The rates of NADPH oxidation in the absence of enzyme were subtracted from the measured rates. NADPH was maintained under saturating conditions, for this reason, reactions were monitored at 380 nm due to the high absorbance of NADPH at 340 nm. Reactions were conducted in triplicate and initiated upon the addition of the enzyme ranging between 6 and 787 µg. One unit of enzyme is defined as the amount of enzyme required to oxidize 1 µmol of NADPH per minute. Each assay was monitored for a time period of 10 min. Stock solutions of substrates were prepared in dimethyl sulfoxide (DMSO) and the amount of substrate added to the cuvette was adjusted to ensure that DMSO concentrations were about 10% of the reaction volume to minimize inhibition of enzyme activity (no more than 9% reduction in enzyme activity).
Samples and standards were analyzed using an HPLC system (Agilent Technology 1200 Series, Palo Alto, CA, USA) equipped with a quaternary pump, an inline degasser, and a diode array detector set at 218 nm. A Phenomenex® 4 µ Jupiter Proteo 90A (250 × 4.6 mm) with a C18 guard column (Torrance, CA, USA) was used for the separation. DON, 3-keto-DON, and 3-epi-DON were eluted using a binary mobile phase set at a flow rate of 1.0 mL/min. The composition of the mobile phase was acetonitrile: water (10:90) for DON/3-epi-DON and acetonitrile:water (20:80) for 3-keto-DON. The injection volume was 10 μL.
The LC–MS/MS analysis was conducted using a Thermo® Scientific Q-Exactive™ Orbitrap mass spectrometer equipped with a Vanquish™ Flex Binary UPLC System (Waltham, MA, USA). Kinetex F5 100A column (150 × 4.6 mm, 2.6 μm) (Phenomenex) was used for separation. The binary mobile phase consisted of solvent A (99.9% H2O/0.1% formic acid, and solvent B (99.9% acetonitrile/0.1% formic acid). The chromatographic elution conditions were as follows: 0–16 min—5% B; 16–17 min—5–100% B; 17–21 min—100% B; 21–22 min—100–5% B; 22–28 min—5% B. The column compartment was kept at 22 °C; the flow rate was set at 0.3 mL/min; the injection volume was 10 μL; UV = 276 nm. The heated electrospray ionization (HESI) source was used in positive mode for the ionization of the target compounds; DDMS (top 10) mode was used, with NCE set at 30.
DepBRleg at a final concentration of 9 µg/mL was incubated with 170 µM of 3-keto-DON, 400 µM NADPH and 50 mM TRIS pH 8, 150 mM NaCl. The reaction was incubated overnight and quenched with acetonitrile. Tubes were vortexed and left to stand for 10 min before filtration with a 0.45 µm filter. Samples were stored at − 20 °C until ready for HPLC analysis.
Coordinates for the DepBRleg structure have been deposited in the Worldwide Protein Data Bank (wwPDB) with the PDB accession code: 7UTF. Data from this study can be requested from the corresponding author Stephen YK Seah.
Keller, N. P. Fungal secondary metabolism: Regulation, function and drug discovery. Nat. Rev. Microbiol. 17, 167–180 (2019).
Article CAS PubMed PubMed Central Google Scholar
Desjardins, A. E., Hohn, T. M. & McCormick, S. P. Trichothecene biosynthesis in Fusarium species: Chemistry, genetics, and significance. Microbiol. Rev. 57, 595–604 (1993).
Article CAS PubMed PubMed Central Google Scholar
Desjardins, A. et al. Mycotoxins: Risks in plant, animal, and human systems. Mycopathologia 65, 2 (2003).
Pierron, A. et al. Microbial biotransformation of DON: Molecular basis for reduced toxicity. Sci. Rep. 6, 1–13 (2016).
Article ADS CAS Google Scholar
Wang, W. et al. The ribosome-binding mode of trichothecene mycotoxins rationalizes their structure: Activity relationships. Int. J. Mol. Sci. 22, 1–17 (2021).
He, J. W. et al. An epimer of deoxynivalenol: Purification and structure identification of 3-epi-deoxynivalenol. Food Addit. Contam. Part A Chem. Anal. Control. Expo. Risk Assess. 32, 1523–1530 (2015).
Article CAS PubMed Google Scholar
He, J. W. et al. Toxicology of 3-epi-deoxynivalenol, a deoxynivalenol-transformation product by Devosia mutans 17-2-E-8. Food Chem. Toxicol. 84, 250–259 (2015).
Article CAS PubMed Google Scholar
Hassan, Y. I. et al. The enzymatic epimerization of deoxynivalenol by Devosia mutans proceeds through the formation of 3-keto-DON intermediate. Sci. Rep. 7, 1–11 (2017).
Carere, J., Hassan, Y. I., Lepp, D. & Zhou, T. The enzymatic detoxification of the mycotoxin deoxynivalenol: Identification of DepA from the DON epimerization pathway. Microb. Biotechnol. 11, 1106–1111 (2018).
Article CAS PubMed Google Scholar
Carere, J., Hassan, Y. I., Lepp, D. & Zhou, T. The identification of DepB: An enzyme responsible for the final detoxification step in the deoxynivalenol epimerization pathway in Devosia mutans 17-2-E-8. Front. Microbiol. 9, 1–9 (2018).
Hoog, S. S., Pawlowski, J. E., Alzari, P. M., Penning, T. M. & Lewis, M. Three-dimensional structure of rat liver 3α-hydroxysteroid/dihydrodiol dehydrogenase: A member of the aldo–keto reductase superfamily. Proc. Natl. Acad. Sci. U. S. A. 91, 2517–2521 (1994).
Article ADS CAS PubMed PubMed Central Google Scholar
Matsunaga, T. et al. Involvement of an aldo–keto reductase (AKR1C3) in redox cycling of 9,10-phenanthrenequinone leading to apoptosis in human endothelial cells. Chem. Biol. Interact. 181, 52–60 (2009).
Article CAS PubMed Google Scholar
Li, D. & Ellis, E. M. Aldo–keto reductase 7A5 (AKR7A5) attenuates oxidative stress and reactive aldehyde toxicity in V79–4 cells. Toxicol. Vitro 28, 707–714 (2014).
Guengerich, F. P. et al. Reduction of aflatoxin B1 dialdehyde by rat and human aldo–keto reductases. Chem. Res. Toxicol. 14, 727–737 (2001).
Article CAS PubMed Google Scholar
Laffin, B. & Petrash, J. M. Expression of the aldo–keto reductases AKR1B1 and AKR1B10 in human cancers. Front. Pharmacol. 3 JUN, 1–7 (2012).
Tang , WH , Martin , KA & Hwa , J. Aldose reductase , oxidative stress , and diabetic mellitus .Front.Pharmacol.3(MAY), 1–8 (2012).
Vander Jagt, D. L., Robinson, B., Taylor, K. K. & Hunsaker, L. A. Reduction of trioses by NADPH-dependent aldo–keto reductases. Aldose reductase, methylglyoxal, and diabetic complications. J. Biol. Chem. 267, 4364–4369 (1992).
Article CAS PubMed Google Scholar
Hyndman, D., Bauman, D. R., Heredia, V. V. & Penning, T. M. The aldo–keto reductase superfamily homepage. Chem. Biol. Interact. 143–144, 621–631 (2003).
Article PubMed CAS Google Scholar
Jez, J. M., Flynn, T. G. & Penning, T. M. A nomenclature system for the aldo–keto reductase superfamily. Adv. Exp. Med. Biol. 414, 579–589 (1997).
Article CAS PubMed Google Scholar
Gerlt, J. A. et al. Enzyme function initiative-enzyme similarity tool (EFI-EST): A web tool for generating protein sequence similarity networks. Biochim. Biophys. Acta Proteins Proteom. 1854, 1019–1037 (2015).
Jez, J. M. & Penning, T. M. The aldo–keto reductase (AKR) superfamily: An update. Chem. Biol. Interact. 130–132, 499–525 (2001).
He, W. J. et al. An aldo–keto reductase is responsible for Fusarium toxindegrading activity in a soil Sphingomonas strain. Sci. Rep. 7, 1–13 (2017).
Bate, N., Butler, A. R., Smith, I. P. & Cundliffe, E. The mycarose-biosynthetic genes of Streptomyces fradiae, producer of tylosin. Microbiology 146, 139–146 (2000).
Article CAS PubMed Google Scholar
Summers, R. G. et al. Sequencing and mutagenesis of genes from the erythromycin biosynthetic gene cluster of Saccharapolyspora erythraea that are involved in l -mycarose and d -desosamine production. Microbiology 143, 3251–3262 (1997).
Article CAS PubMed Google Scholar
Ikeda, H., Nonomiya, T., Usami, M., Ohta, T. & Omura, S. Organization of the biosynthetic gene cluster for the polyketide anthelmintic macrolide avermectin in Streptomyces avermitilis. Proc. Natl. Acad. Sci. U. S. A. 96, 9509–9514 (1999).
Article ADS CAS PubMed PubMed Central Google Scholar
Wonisch, W., Schaur, R. J., Bilinski, T. & Esterbauer, H. Assessment of growth inhibition by aldehydic lipid peroxidation products and related aldehydes by Saccharomyces cerevisiae. Cell Biochem. Funct. 13, 91–98 (1995).
Article CAS PubMed Google Scholar
Semchyshyn, H. M. Reactive carbonyl species in vivo: Generation and dual biological effects. Sci. World J. 2014, 27–31 (2014).
Ciegler, A., Vesonder, R. F. & Jackson, L. K. Production and biological activity of patulin and citrinin from Penicillium expansum. Appl. Environ. Microbiol. 33, 1004–1006 (1977).
Article ADS CAS PubMed PubMed Central Google Scholar
Tannous, J. et al. Patulin transformation products and last intermediates in its biosynthetic pathway, E- and Z-ascladiol, are not toxic to human cells. Arch. Toxicol. 91, 2455–2467 (2017).
Article CAS PubMed Google Scholar
Tannous, J. et al. Sequencing, physical organization and kinetic expression of the patulin biosynthetic gene cluster from Penicillium expansum. Int. J. Food Microbiol. 189, 51–60 (2014).
Article CAS PubMed Google Scholar
Ianiri, G. et al. Searching for genes responsible for patulin degradation in a biocontrol yeast provides insight into the basis for resistance to this mycotoxin. Appl. Environ. Microbiol. 79, 3101–3115 (2013).
Article ADS CAS PubMed PubMed Central Google Scholar
Kadumuri, R. V. & Vadrevu, R. Diversity in αβ and βα Loop connections in TIM barrel proteins: Implications for stability and design of the fold. Interdiscip. Sci. Comput. Life Sci. 10, 805–812 (2018).
Borhani, D. W., Harter, T. M. & Petrash, J. M. The crystal structure of the aldose reductase·NADPH binary complex. J. Biol. Chem. 267, 24841–24847 (1992).
Article CAS PubMed Google Scholar
Krissinel, E. & Henrick, K. Detection of protein assemblies in crystals.In 163–174.https://doi.org/10.1007/11560500_15 (2005).
Obmolova, G. et al. Crystal structure of the Escherichia coli Tas protein, an NADP(H)-dependent aldo–keto reductase. Proteins Struct. Funct. Genet. 53, 323–325 (2003).
Article CAS PubMed Google Scholar
Lapthorn, A. J., Zhu, X. & Ellis, E. M. The diversity of microbial aldo/keto reductases from Escherichia coli K12. Chem. Biol. Interact. 202, 168–177 (2013).
Article CAS PubMed Google Scholar
Wheeler, R. et al. Discovery of a cryptic intermediate in late steps of Mithramycin biosynthesis. Angew. Chemie Int. Ed. 59, 826–832 (2020).
Sanli, G., Dudley, J. I. & Blaber, M. Structural biology of the aldo–keto reductase. Cell Biochem. Biophys. 38, 79–101 (2003).
Article CAS PubMed Google Scholar
Kwon, S. G. et al. Important role of the C-terminal region of pig aldo–keto reductase family 1 member C1 in the NADPH-dependent reduction of steroid hormones. Indian J. Biochem. Biophys. 50, 237–241 (2013).
Penning, T. M. The aldo–keto reductases (AKRs): Overview. Chem. Biol. Interact. 234, 236–246 (2015).
Article CAS PubMed Google Scholar
Ishikura, S. et al. Characterization of two isoforms of mouse 3(17)α-hydroxysteroid dehydrogenases of the aldo–keto reductase family. Biol. Pharm. Bull. 27, 1939–1945 (2004).
Article CAS PubMed Google Scholar
Zhang, Y., Dufort, I. & Rheault, P. Characterization of a human 20-hydroxysteroid dehydrogenase. J. Mol. Endocrinol. 25, 221–228 (2000).
Handumrongkul, C., Ma, D.-P. & Silva, J. L. Cloning and expression of Candida guilliermondii xylose reductase gene (xyl1) in Pichia pastoris. Appl. Microbiol. Biotechnol. 49, 399–404 (1998).
Article CAS PubMed Google Scholar
Kavanagh, K. L., Klimacek, M., Nidetzky, B. & Wilson, D. K. Structure of xylose reductase bound to NAD+ and the basis for single and dual co-substrate specificity in family 2 aldo–keto reductases. Biochem. J. 373, 319–326 (2003).
Article CAS PubMed PubMed Central Google Scholar
Amore , R. , Kötter , P. , Küster , C. , Ciriacy , M. & Hollenberg , CP Cloning and expression of the NAD(P)H-dependent xylose reductase-encoding gene (XYL1) in Saccharomyces cerevisiae . assimilating yeast Pichia stipitis.Gene 109, 89–97 (1991).
Article CAS PubMed Google Scholar
Liu, S. Q., Jin, H., Zacarias, A., Srivastava, S. & Bhatnagar, A. Binding of pyridine nucleotide coenzymes to the β-subunit of the voltage-sensitive K+ channel. J. Biol. Chem. 276, 11812–11820 (2001).
Article CAS PubMed Google Scholar
Ellis, E. M. & Hayes, J. D. Substrate specificity of an aflatoxin-metabolizing aldehyde reductase. Biochem. J. 312, 535–541 (1995).
Article CAS PubMed PubMed Central Google Scholar
Yokochi, N., Yoshikane, Y., Trongpanich, Y., Ohnishi, K. & Yagi, T. Molecular cloning, expression, and properties of an unusual aldo–keto reductase family enzyme, pyridoxal 4-dehydrogenase, that catalyzes irreversible oxidation of pyridoxal. J. Biol. Chem. 279, 37377–37384 (2004).
Article CAS PubMed Google Scholar
Jez, J. M., Bennett, M. J., Schlegel, B. P., Lewis, M. & Penning, T. M. Comparative anatomy of the aldo–keto reductase superfamily. Biochem. J. 326, 625–636 (1997).
Article CAS PubMed PubMed Central Google Scholar
Tipparaju, S. M., Barski, O. A., Srivastava, S. & Bhatnagar, A. Catalytic mechanism and substrate specificity of the β-subunit of the voltage-gated potassium channel. Biochemistry 47, 8840–8854 (2008).
Article CAS PubMed Google Scholar
Ma, H., Ratnam, K. & Penning, T. M. Mutation of nicotinamide pocket residues in rat liver 3α-hydroxysteroid dehydrogenase reveals different modes of cofactor binding. Biochemistry 39, 102–109 (2000).
Article CAS PubMed Google Scholar
Paidimuddala, B., Mohapatra, S. B., Gummadi, S. N. & Manoj, N. Crystal structure of yeast xylose reductase in complex with a novel NADP-DTT adduct provides insights into substrate recognition and catalysis. FEBS J. 285, 4445–4464 (2018).
Article CAS PubMed Google Scholar
Wang, Z. et al. Semi-rational engineering of a thermostable aldo–keto reductase from Thermotoga maritima for synthesis of enantiopure ethyl-2-hydroxy-4-phenylbutyrate (EHPB). Sci. Rep. 7, 1–11 (2017).
ADS PubMed PubMed Central CAS Google Scholar
Grant, A. W., Steel, G., Waugh, H. & Ellis, E. M. A novel aldo–keto reductase from Escherichia coli can increase resistance to methylglyoxal toxicity. FEMS Microbiol. Lett. 218, 93–99 (2003).
Article CAS PubMed Google Scholar
Streeter, C. C., Lin, Q. & Firestine, S. M. Isatins inhibit N 5 -CAIR synthetase by a substrate depletion mechanism. Biochemistry 58, 2260–2268 (2019).
Article CAS PubMed Google Scholar
Xing, M., Chen, Y., Li, B. & Tian, S. Characterization of a short-chain dehydrogenase/reductase and its function in patulin biodegradation in apple juice. Food Chem. 348, 129046 (2021).
Article CAS PubMed Google Scholar
Chan, E. T. S., Zhu, Y., Li, X.-Z., Zhou, T. & Seah, S. Y. K. Characterization of two dehydrogenases from Gluconobacter oxydans involved in the transformation of patulin to ascladiol. Toxins (Basel). 14, 423 (2022).
Baker, P. J., Britton, K. L., Rice, D. W., Rob, A. & Stillman, T. J. Structural consequences of sequence patterns in the fingerprint region of the nucleotide binding fold. Implications for nucleotide specificity. J. Mol. Biol. 228, 662–671 (1992).
Article CAS PubMed Google Scholar
Andersen, K. R., Leska, N. C., Schwartz, T. U. & Optimized, E. coli expression strain LOBSTR eliminates common contaminants from His-tag purification. Proteins 81, 1857–1861 (2013).
Article CAS PubMed PubMed Central Google Scholar
Bradford, M. M. A rapid and sensitive method for the quantitation of microgram quantities of protein utilizing the principle of protein-dye binding. Anal. Biochem. 72, 248–254 (1976).
Article CAS PubMed Google Scholar
Kabsch, W. XDS. Acta Crystallogr. Sect. D Biol. Crystallogr. 66, 125–132 (2010).
Keegan, R. M. & Winn, M. D. MrBUMP: An automated pipeline for molecular replacement. Acta Crystallogr. Sect. D Biol. Crystallogr. 64, 119–124 (2007).
Liebschner, D. et al. Macromolecular structure determination using X-rays, neutrons and electrons: Recent developments in Phenix. Acta Crystallogr. Sect. D Struct. Biol. 75, 861–877 (2019).
Emsley, P., Lohkamp, B., Scott, W. G. & Cowtan, K. Features and development of Coot. Acta Crystallogr. Sect. D Biol. Crystallogr. 66, 486–501 (2010).
Holm, L., Kääriäinen, S., Wilton, C. & Plewczynski, D. Using Dali for structural comparison of proteins. In Current Protocols in Bioinformatics (eds Baxevanis, A. D. et al.) (Wiley, 2006). https://doi.org/10.1002/0471250953.bi0505s14.
Okonechnikov, K. et al. Unipro UGENE: A unified bioinformatics toolkit. Bioinformatics 28, 1166–1167 (2012).
Article CAS PubMed Google Scholar
Edgar, R. C. MUSCLE: A multiple sequence alignment method with reduced time and space complexity. BMC Bioinform. 5, 1–19 (2004).
Edelheit, O., Hanukoglu, A. & Hanukoglu, I. Simple and efficient site-directed mutagenesis using two single-primer reactions in parallel to generate mutants for protein structure-function studies. BMC Biotechnol. 9, 1–8 (2009).
Copp, J. N., Akiva, E., Babbitt, P. C. & Tokuriki, N. Revealing unexplored sequence-function space using sequence similarity networks. Biochemistry 57, 4651–4662 (2018).
Article CAS PubMed Google Scholar
Li, W., Jaroszewski, L. & Godzik, A. Clustering of highly homologous sequences to reduce the size of large protein databases. Bioinformatics 17, 282–283 (2001).
Article CAS PubMed Google Scholar
Shannon, P. et al. Cytoscape: A software environment for integrated models. Genome Res. 13, 426 (1971).
Turner, M., Heney, K. A. & Rod Merrill, A. The N-terminus of Paenibacillus larvae C3larvinA modulates catalytic efficiency. Biosci. Rep. 41, 1–16 (2021).
Wang, W. & Seah, S. Y. K. Purification and biochemical characterization of a pyruvate-specific class II aldolase, HpaI. Biochemistry 44, 9447–9455 (2005).
Article CAS PubMed Google Scholar
The authors would like to thank Xiu-Zhen Li and Lili Mats (Agriculture and AgriFood Canada) as well as Dr. Rod Merrill and Chau Tran (University of Guelph). The financial support from Agriculture and Agri-Food Canada and Grain Farmers of Ontario through the Industry-led Research and Development Project (J-002433) is acknowledged and much appreciated.
Department of Molecular and Cellular Biology, University of Guelph, Guelph, Canada
Nadine Abraham, Kurt L. Schroeter, Jonathan Chan, Natasha Evans, Matthew S. Kimber & Stephen YK Seah
Guelph Research and Development Centre, Agriculture and Agri-Food Canada, Guelph, ON, Canada
Nadine Abraham, Yan Zhu, Jonathan Chan, Natasha Evans, Jason Carere & Ting Zhou
You can also search for this author in PubMed Google Scholar
You can also search for this author in PubMed Google Scholar
You can also search for this author in PubMed Google Scholar
You can also search for this author in PubMed Google Scholar
You can also search for this author in PubMed Google Scholar
You can also search for this author in PubMed Google Scholar
You can also search for this author in PubMed Google Scholar
You can also search for this author in PubMed Google Scholar
You can also search for this author in PubMed Google Scholar
N.A. and S.Y.K.S. designed the study. N.A. and S.Y.K.S. wrote the manuscript and prepared the figures. N.A., K.S., Y.Z. J.C. and N.E. conducted the reported experiments. N.A. and M.S.K. performed molecular modeling. N.A. and S.Y.K.S. interpreted the results. N.A. and S.Y.K.S. integrated all research. J.C. reviewed and edited the manuscript. T.Z. obtained research funding and is involved in the conception and supervision of the research project. All authors discussed the results and commented on the manuscript.
Correspondence to Stephen Y. K. Seah.
The authors declare no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/.
Abraham, N., Schroeter, K.L., Zhu, Y. et al. Structure–function characterization of an aldo–keto reductase involved in detoxification of the mycotoxin, deoxynivalenol. Sci Rep 12, 14737 (2022). https://doi.org/10.1038/s41598-022-19040-8
DOI: https://doi.org/10.1038/s41598-022-19040-8
Anyone you share the following link with will be able to read this content:
Sorry, a shareable link is not currently available for this article.
Provided by the Springer Nature SharedIt content-sharing initiative
By submitting a comment you agree to abide by our Terms and Community Guidelines. If you find something abusive or that does not comply with our terms or guidelines please flag it as inappropriate.
Scientific Reports (Sci Rep) ISSN 2045-2322 (online)
6as)-4-Formylhexahydro-2-Oxo-2h-Cyclopenta[B]Furan-5-Yl-Ester Sign up for the Nature Briefing: Translational Research newsletter — top stories in biotechnology, drug discovery and pharma.