Citation

Luo Pan, Lv Jian-Ming, Zhen Hong-Ting, Zhao Ying-Qi, Liu Jing-Yuan, Hong Jin-Yu, Li Shao-Yang, Wang Gao-Qian, Chen Guo-Dong, Zhang Shui-Xing, Hu Dan, Gao Hao. Identification of two novel sesterterpene skeletons offers the first experimental evidence for the cyclization mechanism of mangicdiene synthase[J]. Chin. Chem. Lett., 2026, 37(1): 111042  

Pan Luo^a,b,c,1, Jian-Ming Lv^b,1, Hong-Ting Zhen^b,1, Ying-Qi Zhao^b, Jing-Yuan Liu^b, Jin-Yu Hong^b, Shao-Yang Li^b, Gao-Qian Wang^b, Guo-Dong Chen^b,c, Shui-Xing Zhang^a,*, Dan Hu^b,*, Hao Gao^a,b,*     

Terpenoids are a remarkable class of natural products derived from linear polyprenyl diphosphate precursors. Impressively, they exhibit diverse structural complexity and a wide range of biological functions [1–3]. These fascinating molecules are formed through terpene synthases catalyzed cyclization cascade reactions, followed by tailoring enzymes mediated post-cyclization modifications [4–11], where terpene synthases play a pivotal role in constructing complicated polycyclic terpene skeletons. Therefore, there has been considerable interests in study on the terpene cyclization mechanism [12–14], which can set the stage for generation of novel terpene skeletons through protein engineering.

Mangicol-type sesterterpenoids are a small group of terpenoids, which are characterized with the 5–5–6–5 tetracyclic skeleton and exclusively found in fungi [15–17]. Notably, the tetracyclic skeleton of mangicols closely resembles that of the diterpene allokutznerene, except for the position of the double bond and an additional isoprene unit [18]. Thus far, 13 mangicol-type sesterterpenoids have been reported (Fig. S1 in Supporting information) [15–17], some of which possess anti-cancer and anti-inflammatory activities [15,17]. Among them, mangicol J showed comparable inhibitory activity to the pharmaceutical agent indomethacin in terms of treating phorbol myristate acetate-induced oedema in the mouse ear edema assay [17]. Therefore, biosynthesis of mangicols draws the extensive attention. Previously, the genetic basis for formation of the mangicol-type skeleton was unveiled. The chimeric sesterterpene synthase FgMS, which is composed of the prenyltransferase domain and terpene cyclase domain as other known fungal sesterterpene synthases [19], firstly elongates short prenyl chains into geranylfarnesyl diphosphate (GFPP), and then catalyzes cyclization of GFPP to mangicdiene [16]. In addition, several cytochrome P450 enzymes have also been identified to be responsible for biosynthesis of mangicol-type sesterterpenoids [17].

Regarding the formation of mangicdiene, it is well accepted that the cyclization begins with the departure of the pyrophosphate group of GFPP, followed by a concerted C1-C11/C10-C14 bicyclization to generate carbocation A⁺ (Fig. 1), which has been proposed to undergo a sequential 1,2-alkyl rearrangement to yield the C10 carbocation C⁺ via B⁺ [16]. However, the density functional theory (DFT) calculation revealed that A⁺ is likely to undergo two asynchronous events, including the 1,2-alkyl shift and C—C bond cleavage, to give D⁺ (Fig. 1), as the intermediate B⁺ is not a minimum on the potential energy surface [20]. This is in line with isolation of variecoltetraene, which is likely derived from E⁺ via B⁺. Subsequently, D⁺ is converted to the nonclassical C10 carbocation intermediate F⁺, wherein the single bond between C14 and C15, as well as the double bond between C10 and C14, is partially broken. After C10 carbocation formation, both C6-C10 cyclization (path a) and C2-C10 cyclization (path b) are able to lead to formation of mangicdiene (Fig. 1). The path a involves C6-C10 cyclization to G⁺, 1,2-hydride shift followed by C2-C6 cyclization to H⁺, 1,2-hydride shift to I⁺, and the final deprotonation of I⁺ (Fig. 1) [16]. Previously, the diterpene allokutznerene was also hypothesized to be formed in a similar manner [18]. In contrast, the path b seems a little prolix yet energetically reasonable, which requires C2-C10 cyclization to J⁺, 1,2-hydride shift followed by C2-C6 cyclization to K⁺, 1,2-hydride and alkyl shift to I⁺, and the final deprotonation of I⁺ (Fig. 1) [20]. Nevertheless, due to the lack of experimental evidence, it remains inconclusive whether mangicdiene is formed via the concise path a or the energetically favorable path b.

Figure 1

	Download: JPG larger image
Fig. 1. Previously proposed and herein experimentally supported biosynthetic mechanisms of mangicdiene (1).

Due to being able to unequivocally determine the origin of carbons and hydrogens in the terpene skeleton, isotopic labelling experiments have been widely employed to unveil the terpene cyclization mechanism [21–28]. Alternatively, as carbocation intermediates are likely to be quenched via deprotonation or hydration to yield shunt products at any stage of the cascade cyclization reaction, identification of shunt products generated by the terpene cyclase or its variants can be used to confirm the presence of key carbocation intermediates, and then facilitate the elucidation of the cyclization mechanism [25,26,29–32]. For the case of biosynthesis of mangicdiene, as all the carbons and hydrogens in mangicdiene from the path a or path b share the same origin, isotopic labelling experiments cannot discern the specific pathway through which mangicdiene is formed. Therefore, the viable approach is to endeavor in the identification of shunt products related to mangicdiene. Previously, when FgMS was expressed in Escherichia coli (E. coli), the bicyclic shunt product variecoltetraene was concomitantly identified (Fig. 1) [16]. In addition, several trace shunt products were also detected, but none of them were isolated as a result of their rather low yields [16]. Considering the potential variability of the terpene cyclase product profile in different heterologous expression hosts, such as the notable difference observed when the amorpha-4,11-diene synthase (ADS) was expressed in Aspergillus nidulans and E. coli [33], the product profile of the mangicdiene synthase might also be altered in different hosts, thus resulting in isolation and identification of new shunt products.

During our continuous mining of terpene cyclases [6,21,34–38], we located a candidate sesterterpene cyclase ManD from Fusarium sp. JNU-XJ070152–01, sharing the 93% amino acid identity with FgMS, also termed as MgcD [16,17]. Moreover, the tailoring enzymes encoded by genes in the flanking region of manD are also highly homologous (>90% identity) to those encoded by the mgc cluster (Fig. S2 in Supporting information), which are dedicated to biosynthesis of mangicols [17]. Based on these, we proposed that the man cluster (Genbank accession number: PQ623309) is functionally the same as the mgc cluster, and thus decided to utilize the mangicdiene synthase gene manD in hand to probe its product profile in Saccharomyces cerevisiae and A. oryzae, both of which have been widely used for heterologous expression of fungal sesterterpene synthases [17,25,27,29,31,35,39–45]. To this end, the intronless manD gene was prepared (Supplementary methods and Table S1 in Supporting information), and then incorporated into the engineered S. cerevisiae strain [34]. Upon gas chromatography-mass spectrometry (GC–MS) analysis, we found that the resulting transformant can only produce one sesterterpene hydrocarbon, the MS fragment pattern of which is completely consistent with that of mangicdiene (Fig. S3 in Supporting information) [16]. The result indicated that heterologous expression of manD in the engineered S. cerevisiae strain cannot accumulate any other related shunt products. We then turned to express manD in A. oryzae NSAR1 [46]. Intriguingly, high-performance liquid chromatography (HPLC) analysis showed that the A. oryzae NSAR1 transformant possessing manD could yield one major product 1 along with three minor products 2−4 (Fig. 2). Among them, 1 and 2 were individually identified as mangicdiene and variecoltetraene by comparing the NMR data and specific optical rotation with those reported values (Supplementary methods and Figs. S4−S7 in Supporting information) [16], while 3 and 4 were elucidated as two novel sesterterpene skeletons by exhaustive nuclear magnetic resonance (NMR) analysis, and named fusadiene and fusatriene, respectively (Tables S2−S4 and Figs. S8−S22 in Supporting information).

Figure 2

	Download: JPG larger image
Fig. 2. HPLC analysis of mycelial metabolites from the A. oryzae NSAR1 transformant harboring manD and the parent strain.

Fusadiene (3) was obtained as colorless oil, and its molecular formula was determined to be C₂₅H₄₀ (six degrees of unsaturation) based on the high resolution electron ionization mass spectrometry (HR-EI-MS) ion at m/z 340.31253 [M]^•+ (calcd. for C₂₅H₄₀, 340.31245). In the ¹H and ¹³C NMR spectra, one olefinic proton signal [δ_H 5.12 (1H, br t, J = 7.2 Hz)] and four olefinic carbon signals (δ_C 140.8, 130.8, 125.4, 123.1) were observed, implying that 3 possesses the tetracyclic skeleton with two double bonds. The ¹H−¹H homonuclear correlation spectroscopy (COSY) spectrum showed the presence of four spin coupling systems C-24–C-3–C-4–C-5, C-8–C-9–C-10–C-14, C-12–C-13, and C-16–C-17–C-18. And in the heteronuclear multiple-bond correlation (HMBC) spectrum, the correlations from H₃–20 to C-18/C-19/C-21, from H₃–21 to C-18/C-19/C-20, from H₃–22 to C-13/C-14/C-15/C-16, from H₃–23 to C-6/C-7/C-8, from H₃–24 to C-2/C-3/C-4, from H₃–25 to C-1/C-11/C-12/C-14, from Hb-1 to C-2/C-10, from Ha-5/Hb-5 to C-6, and from H-10 to C-2/C-3/C-6 were observed. The assignments of C-20 and C-21 were determined on the basis of the nuclear Overhauser effect (NOE) correlation between H-18 and H₃–20. Accordingly, the planar structure of 3 was established (Fig. 3A). The NOE correlations between Ha-8 and H-14/H₃–25 indicated that the rings B and C are fused in a cis configuration, and that H-10/C-3 and H-14/C-25 are located on the different sides of the ring C. The NOE correlations between H-10 and H₃–22, and between H-14 and H₂–17 revealed that in comparison with H-14 and C-25, C-22 is oriented in the opposite direction. However, the relative configuration of C-3 cannot be determined by NOE analysis. We then performed quantum chemical calculation of NMR chemical shifts. By comparison of the experimental ¹³C NMR data of 3 with those calculated values of two putative stereoisomers (Fig. 4), the relative configuration of 3 was determined to be 2S*,3S*,10S*,11R*,14R*,15S*. As both 1 and 3 are produced by ManD, and they should share the early-stage cyclization process, we thus assigned the absolute configuration of 3 as 2S,3S,10S,11R,14R,15S.

Figure 3

	Download: JPG larger image
Fig. 3. Key 1H−1H COSY, HMBC, and NOE correlations of 3 (A) and 4 (B).

Figure 4

	Download: JPG larger image
Fig. 4. Comparison of the experimental 13C chemical shifts of 3 with calculated values of (2S*,3S*,10S*,11R*,14R*,15S*)−3 and (2S*,3R*,10S*,11R*,14R*,15S*)−3.

Fusatriene (4) was obtained as colorless oil, and the HR-EI-MS ion at m/z 340.31253 [M]^•+ (calcd. for C₂₅H₄₀, 340.31245) revealed its molecular formula as C₂₅H₄₀ (six degrees of unsaturation). The presence of six olefinic carbon signals (δ_C 149.2, 136.3, 130.9, 125.3, 123.2, 112.0), as well as four olefinic proton signals [δ_H 5.30 (1H, br dd, J = 12.0, 3.2 Hz), 5.12 (1H, br t, J = 6.6 Hz), 4.83 (1H, br s), 4.77 (1H, br s)], in the 1D NMR spectra demonstrated that 4 is featured with the tricyclic skeleton and three double bonds. The ¹H−¹H COSY spectrum showed that there are four spin coupling systems C-1–C-2–C-10(–C-14)–C-9–C-8, C-4–C-5–C-6, C-12–C-13, and C-16–C-17–C-18. In combination with HMBC correlations from H₃–20 to C-18/C-19/C-21, from H₃–21 to C-18/C-19/C-20, from H₃–22 to C-13/C-14/C-15/C-16, from H₃–23 to C-6/C-7/C-8, from H₃–25 to C-1/C-11/C-12/C-14, and from Ha-24/Hb-24 to C-2/C-3/C-4, and the NOE correlation between H-18 and H₃–20, the planar structure of 4 was established (Fig. 3B). The NOE correlations between H-2 and Hb-13/H₃–22 revealed that the rings B and C are cis fused, with H-2 and H-14/C-25 being positioned on the opposite sides of the ring B, and C-22 and H-14/C-25 being situated on the opposite sides of the ring C. And based on the NOE correlation between H-10 and H₃–22, it can be determined that H-10 is oriented in a direction distinct from that of H-14. Taken together, the relative configuration of 4 was defined as 2R*,10R*,11R*,14R*,15S*. In addition, the geometrical configuration of the Δ⁶-double bond was assigned as E, due to the NOE correlation between Ha-5 and H₃–23. Likewise, according to the biosynthetic relationship between 1 and 4, the absolute configuration of 4 was assigned as 2R,10R,11R,14R,15S.

Structural characterization of 3 and 4 proves the occurrence of J⁺ and K⁺ in the cyclization process from GFPP to mangicdiene, experimentally supporting the involvement of the path b for construction of mangicdiene (Fig. 1). Generally, isotopic labelling experiments serve as a powerful approach to unambiguously elucidate the cyclization mechanism of terpene synthases regardless of identification of shunt products [43,47,48]. Taking biosynthesis of the diterpene cattleyene as example, several carbons and hydrogens of cattleyene from two plausible pathways did not originate from the same source, and thus the cyclization mechanism for cattleyene was elucidated through extensive isotopic labelling experiments [47]. However, for some exceptional cases, where isotopic labelling experiments are not applicable, identification of shunt products is required to uncover the mystery [49,50]. For example, two plausible routes were presumed to be responsible for biosynthesis of the sesquiterpene β-terrecyclene, neither of which could be ruled out by isotopic labelling experiments. Therefore, to probe the mechanism involved in the formation of β-terrecyclene, several potential key amino acids in the active pocket of the β-terrecyclene synthase TerA was mutated, which led to detection of two additional products. Consequently, in combination with the isotopically sensitive branching experiments, the mechanism for β-terrecyclene biosynthesis was demonstrated [50]. Our work represents another rare example showing the significance of identifying shunt products, instead of isotopic labelling experiments, in elucidating the terpene cyclization mechanism. The other intriguing point is that heterologous expression of the mangicdiene synthase gene in different hosts can lead to alteration of its product profile, which might be attributed to the variation of intracellular conditions (such as pH, temperature, and metal ion species) in different hosts [38,51–53]. However, the similar phenomenon was few reported [33]. Notably, as ManD are not completely identical to FgMS, we could not rule out the possibility that alteration of the product profile was induced by the variance of their amino acid sequences.

In conclusion, we identified a candidate sesterterpene synthase ManD from Fusarium sp. JNU-XJ070152–01, which is highly homologous to the mangicdiene synthase FgMS. Through heterologous expression, we found that ManD is functionally the same as FgMS, but the product distribution of ManD in A. oryzae NSAR1 slightly differs from the previously reported product profile of FgMS in E. coli. As a result, we isolated and identified two novel sesterterpene skeletons fusadiene and fusatriene from the A. oryzae NSAR1 transformant expressing manD. The results, for the first time, experimentally confirmed that during the cyclization process from GFPP to mangicdiene, the C10 cation intermediate preferred to undergo C2-C10 cyclization (path b) rather than C6-C10 cyclization (path a).

The authors declare that they have no known competing financial interests or personal relationships that could have appeared to influence the work reported in this paper.

Pan Luo: Writing – original draft, Validation, Investigation. Jian-Ming Lv: Writing – original draft, Validation. Hong-Ting Zhen: Investigation. Ying-Qi Zhao: Investigation. Jing-Yuan Liu: Validation. Jin-Yu Hong: Formal analysis. Shao-Yang Li: Validation. Gao-Qian Wang: Validation. Guo-Dong Chen: Formal analysis. Shui-Xing Zhang: Writing – review & editing, Supervision, Conceptualization. Dan Hu: Writing – review & editing, Validation, Supervision, Conceptualization. Hao Gao: Writing – review & editing, Supervision, Conceptualization.

We thank Prof. K. Gomi (Tohoku University) and Prof. K. Kitamoto (The University of Tokyo) for providing the A. oryzae NSAR1 heterologous expression system. This work was financially supported by grants from the National Key Research and Development Program of China (No. 2024YFE0102000), the National Natural Science Foundation of China (Nos. 81925037, 82321004, U22A20371, U24A20782, 32170060, 22177037, 22207039, 22307045), the Guangdong Major Project of Basic and Applied Basic Research (No. 2023B0303000026), the Guangdong Natural Science Funds for Distinguished Young Scholars (No. 2022B1515020028, China), the Guangdong International Science and Technology Cooperation Base (No. 2021A0505020015, China), the Guangdong Basic and Applied Basic Research Foundation (Nos. 2023B1515040016, 2023A1515110388), the Innovative and Research Teams Project of Guangdong Higher Education Institution (No. 2021KCXTD001, China), the Guangzhou Science and Technology Project (Nos. 202206010020, 2024A04J6241, 2023A04J0080, China), the Fundamental Research Funds for the Central Universities (Nos. 21623105, 21624210).

Supplementary material associated with this article can be found, in the online version, at doi:10.1016/j.cclet.2025.111042.