Ke
Ma
a,
Jie
Liu
a,
Zequan
Huang
a,
Mengyue
Wu
a,
Dong
Liu
a,
Jinwei
Ren
b,
Aili
Fan
*a and
Wenhan
Lin
*ac
aState Key Laboratory of Natural and Biomimetic Drugs, School of Pharmaceutical Sciences, Peking University, Beijing 100191, China. E-mail: fanaili@bjmu.edu.cn; whlin@bjmu.edu.cn
bState Key Laboratory of Mycology, Institute of Microbiology, Chinese Academy of Sciences, Beijing 100101, China
cNingbo Institute of Marine Medicine, Peking University, Ningbo 315832, Zhejiang, China
First published on 11th October 2024
Enzymes from the nuclear transport factor 2-like (NTF2-like) superfamily represent a rare group of biocatalysts with diverse catalytic functions facilitating intriguing skeleton formations. However, most proteins of this family remain enigmatic and await further elucidation. In this study, a combination of protein structural alignment with clustering analysis uncovers a new aldolase, AtoB, belonging to the NTF2-like superfamily. AtoB catalyzes the key intramolecular aldol reaction in linear tetracyclic meroterpenoid biosynthesis. The X-ray crystal structures of AtoB and AtoB-ligand complex are established at 1.9 Å and 1.6 Å resolution, respectively, revealing the rotation of the α4 helix and key residues in the active site for substrate binding. Molecular docking and site-directed mutagenesis demonstrate an acid–base pair involved in the AtoB-catalyzed aldol reaction, of which Arg59 is responsible for stereocontrol of hydroxylated C-10a during condensation. These findings provide valuable information for understanding the catalytic mechanisms of the AtoB-catalyzed aldol reaction. Additionally, a branching biosynthetic pathway of aspertetranones is elucidated during the exploration of the natural substrate of AtoB.
Among these understudied enzyme families, the nuclear transport factor 2-like (NTF2-like) superfamily, featuring a natural small alpha–beta fold first observed in the structure of the rat NTF2 protein, showcases astonishing attractiveness.8 Recently, a few fungal NTF2-like enzymes with distinct catalytic functions have been characterized. For example, SdnG, BvnE, NsrQ and Trt14 catalyze the Diels–Alder reaction, semipinacol rearrangement, isomerization, and intramolecular methoxy rearrangement in the biosynthesis of sordarin, brevianamides, blennolides and terretonins, respectively (Fig. 1A).9–12 It is noteworthy that genes encoding NTF2-like proteins are widely distributed among fungal and bacterial genomes, but with low similarity at the amino acid level. These facts imply that the undiscovered NTF2-like enzymes may hold the potential to catalyze intricate reactions, necessitating the development of a universal strategy to uncover new members of the NTF2-like superfamily.
Herein, stemming from protein structural alignment and CLANS clustering method, we identified a novel NTF2-like aldolase, AtoB, involved in the biosynthesis of aspertetranone A from the marine-derived fungus Aspergillus ochraceus LZDX-32-15. Through target gene deletion in A. ochraceus and heterologous expression in A. nidulans, a branching biosynthetic pathway of aspertetranones was charted and AtoB was found to catalyze an intramolecular adol reaction to form the linear tetracycle of aspertetranones. Surprisingly, a recent online article reported an isozyme OchL,13 which is assumed to catalyze the same reaction to form ochraceopones based on the evidence from three-protein (OchGJL) bioconversion (Fig. 1B). However, the cryptic substrate and the mechanistic details of this new class of intramolecular NTF2-like aldolase remain unknown. To investigate the molecular basis for the AtoB-catalyzed reaction, we performed in vitro characterization of this enzyme using the isolated natural substrate and solved the X-ray crystal structures of AtoB and its complex with the substrate analogue. Structural analyses indicated that the rotatable α4 helix played a key role in substrate binding, and Arg59–Tyr114 based acid–base chemistry was employed for AtoB to catalyze the aldol reaction. Notably, Arg59 is also critical for exerting stereochemical control at C-10a via hydrogen bonds. Combining the structural analyses and mutational studies, a mechanism for this novel NTF2-like aldolase AtoB was proposed.
Specifically, the superfamily members SdnG, BvnE, Trt14, and PrhC were used as probes to search for relevant NTF2-like enzymes against our in-house marine fungal genomic database. Ten putative protein sequences (with 10–40% amino acid identity) were obtained. However, the catalytic functions of these candidate proteins were hardly distinguished due to the diverse catalytic and non-catalytic roles and low peptide sequence similarities within the NTF2-like superfamily. Therefore, the structures of the 10 candidate proteins were predicted by RoseTTAFold2 for subsequent Foldseek searching.16 A total of 1432 deduplicated protein sequences were then obtained and processed for the sequence clustering by CLANS.14 The NTF2-like superfamily was segregated into nine main clusters, each containing more than 25 members (Fig. 2A). They presented diverse biological functions such as nuclear transport factors (cluster 1), catalytic enzymes (clusters 2–5), kinases (cluster 6), nutrient uptake related proteins (cluster 7), and cell shape related proteins (cluster 8). The known fungal NTF2-like enzymes (BvnE, NsrQ, and Trt14) and one of the candidate proteins (named AtoB) clearly constituted a distinct sequence cluster (cluster 2), suggesting that AtoB might have a potential catalytic function. The phylogenetic analysis of sequences in cluster 2 revealed that these NTF2-like members were classified into three clades, with Trt14, BvnE and PrhC in clade 1 and NsrQ and Dcr3 in clade 2 (Fig. 2B and S1†).
Fascinatingly, AtoB significantly differed from the known fungal NTF2-like enzymes and was located in the distinct clade 3, further implying that AtoB could process a unique function.
Further analysis of the up-stream and down-stream genes of atoB in the genome of A. ochraceus LZDX-32-15 enabled postulation of a biosynthetic gene cluster (ato cluster, accession number PQ310357) containing a polyketide synthase (AtoC), a prenyltransferase (AtoE), a FAD-dependent monooxygenase (AtoH) and a terpene cyclase (AtoK) and eight more putative proteins (AtoA, AtoD, AtoF, AtoG, AtoI, AtoJ, AtoL, and AtoM; Table S1†). Significantly, the ato cluster showed low homology to the known fungal meroterpenoid BGCs (Fig. S2†).17–23
To experimentally validate the products of ato cluster, the genes from atoA to atoM were co-expressed in A. nidulans LO8030, resulting in the mutant strain An-atocluster. As shown in Fig. 3, a product (1) was detected in the cultured An-atocluster and identified as aspertetranone A, while the A. ochraceus produced 1 and aspertetranone D (2). To verify the functions of atoC, atoE, atoH, and atoK as classical biosynthetic genes of fungal meroterpenoids, these four genes were co-expressed in A. nidulans to obtain the An-atoCEHK strain. Consistent with the hypothesis, An-atoCEHK generated a meroterpenoid precursor 3. These results suggested that ato genes might share comparable functions and similarity with och cluster, whose sequence is not released now (Fig. S3†).13 Compared to och cluster, AtoB might catalyze the same aldol reaction as the assumed isozyme OchL. However, the function and catalytic mechanism of OchL remained elusive due to the lack of a natural substrate.
Therefore, thirteen single- or double-knockout strains were constructed to explore the natural substrate of AtoB including ΔatoD, ΔatoF, ΔatoG, ΔatoI, ΔatoJ, ΔatoL, ΔatoM, ΔatoGΔatoD, ΔatoGΔatoL, ΔatoDΔatoL, ΔatoIΔatoM, ΔatoFΔatoM, and ΔatoBΔatoM (Fig. 3, S3 and S4†). The metabolites of these mutant strains revealed a branching biosynthetic pathway of aspertetranones and five new derivatives (7, 8, 10, 11, and 14), based on the substrate promiscuity of the short-chain dehydrogenase AtoL and cytochrome P450 AtoM (Fig. S3†).
The mutant strain ΔatoM accumulated 4a in addition to an unreported minor 4b, indicating that AtoM catalyzed the hydroxylation and ketonization of the allylic position (C-12) of 4a (Fig. 3 and Table S7†). In the double-knockout strain ΔatoBΔatoM and heterologous expression mutant An-atoCEHKGDLIF, three main products were detected and characterized by LC-MS and NMR analysis, demonstrating that both mutants failed to produce the linear tetracyclic compound 4a or 4b. Instead, compound 16 and two shunt products, 6a and 6b, with a γ-lactone ring significantly accumulated as reported before (Fig. 3A and S4†).13 These results confirmed that AtoB could catalyze the same aldol reaction as the assumed isozyme OchL. Inspiringly, a trace compound 5 together with two metabolic intermediates 19 and 23 were detected, which could be the natural substrate of AtoB (Fig. 3A and S4†). After a scale-up fermentation of ΔatoBΔatoM strain, 5 was obtained with a molecular weight of 420, similar to 4a and 4b. Thus, we hypothesized that AtoB catalyzed the aldol reaction with 5 as a natural substrate to form the linear tetracyclic epimers 4a and 4b.
To verify the hypothesis, the His6-tagged AtoB was expressed in Escherichia coli BL21(DE3) strain and purified by Ni-NTA chromatography. The recombinant AtoB was incubated with the isolated 5, 19, and 23, respectively. The LC-MS analysis of reaction mixtures exhibited the conversion of 5 into epimers 4a and 4b, confirming that AtoB readily accepted 5 as the substrate (Fig. 3C). These results also illustrated that AtoB catalyzed the stereoselective aldol reaction of 5 to generate the 10aS-products with the chirality of C-6a retained (4a) or flipped (4b).
Notably, comparison of the structures of apo-AtoB and AtoB–19 complex structures showed that they entirely merged, except for the position of α4 helix (Fig. 4A). The α4 helix from the AtoB–19 complex rotated closer to the active cavity for an approximate 6.2°, featuring the side chains of Arg72 and His68 extending into the active cavity (Fig. 4B). This conformational change implied that the hydrogen bonding interactions between Arg72/His68 and the substrate might play important roles in the substrate binding process.
For the aldol reaction catalyzed by AtoB, Tyr114 and Arg59 facilitated the condensation with Tyr114 positioned to interact with the Cα of the C-7 carbonyl group and Arg59 positioned to interact with the C-10a carbonyl group (Fig. 5B). The previously reported fungal NTF2-like enzymes usually share a similar acid–base pair to initiate the reaction, even though they catalyzed different reactions by the deprotonation of a hydroxyl group or Cα of a carbonyl group.10,11,26 Thus, Arg59–Tyr114 was considered as the catalytic acid–base pair to catalyze the aldol reaction and to form a 6/6/6/6 linear tetracyclic skeleton in aspertetranones.
To verify the binding mode of the substrate observed in the docking simulation, we systematically conducted the site-directed mutagenesis of AtoB and compared their activities to that of the wild type. The hydrophilic residues Arg59, Tyr114, Arg72, Glu136, Glu143, and His68 were substituted with Ala or other residues. Compared to the wild-type AtoB, the R72A mutant exhibited 6.65% yield of 4a, while R72Q maintained 51.20% yield of 4a (Fig. 5C and Table S6†). Similarly, the E136A, E136D, E136Q, E143A, E143D, and E143Q mutants showed the productivity of 4a with the range of 0% to 27.55%. These results indicated that Arg72, Glu136, and Glu143 played important roles in the substrate binding, especially the residue Glu136 which might be a key unit to fix the flexible ten-membered ring of 5. Additionally, the H68A mutant maintained 76.53% of 4a production, proving relatively weak π–π stacking interaction between the imidazole and the pyrone ring.
The activity of the Y114F variant was complete regression indicating that Tyr114 might trigger the aldol reaction as its ionized form to serve as a Lewis base to promote the deprotonation of Cα of the C-7 carbonyl group (Fig. 5D), following the spontaneous isomerization into the enol form with nucleophilicity. During this process, the hydrogen bond networks, especially the interactions between the carboxyl group of Glu136 and C-7 carbonyl/C-6a hydroxyl fix the conformation of the substrate to allow the nucleophilic enol attack at the C-10a carbonyl. Comparatively, Arg59 substituted R59A variant resulted in 5.91% yield of 4a. These in vitro assays revealed that Arg59 served as a proton donor for the C-10a carbonyl group, to finish the aldol reaction rather than as an initiator. Arg59 of AtoB is an unusual proton donor, as the calculated pKa value of its sidechain in solution is about 12.93.27 However, we noted that Arg329 of citrate synthase also acts as the proton donor with high-level QM/MM modelling evidence.28 Surprisingly, the R59Q variant maintained 26.74% yield of 4a. The side chain of glutamine might assist in proton transfer by interacting with a water molecule, which could also act as a proton donor. Furthermore, the interaction between Arg59 and C-10a carbonyl group was nearly perpendicular to the ten-member ring, setting an obstacle for conformational changes of the ten-membered ring and preventing the carbonyl group at C-10a from flipping to the opposite side of the ring plane (Fig. 5B and E). As a result, the stereocontrol of hydroxylated C-10a in the aldol product 4a and 4b was determined by the position of Arg59, suggesting that only a little distortion is satisfied when the nucleophile attacks from the si-face during the C–C bond formation step.
Interestingly, the enzyme reactions showed the ratios of 4a and 4b to be 78:22 for wild-type AtoB, while the ratios changed to 35:65 and 55:45 for E143D and E143Q, respectively (Fig. 5D and Table S6†). A plausible explanation is that variations of the hydrogen bond donor allowed a more stretched conformation of the ten-membered ring, resulting in reduction of the spatial hindrance of C-6a hydroxyl in 4b. Meanwhile, variants R72A, R72Q, and H68A produced a similar ratio of 4a and 4b (95:5). Accordingly, the loss of interactions between Arg72/His68 and the pyrone ring of 5 changed the location of the substrate and enhanced the hydrogen bonds between Glu136 and C-6a hydroxyl group, which prompted the retained chirality of C-6a. These observations indicated that the substrate was fixed in the active site through hydrogen bond networks, which governed the C-6a S or R configuration of 4a and 4b.
Footnote |
† Electronic supplementary information (ESI) available. See DOI: https://doi.org/10.1039/d4sc05590j |
This journal is © The Royal Society of Chemistry 2024 |