Peniexpansones A–F: polyketides from Penicillium expansum DWS880 capable of enhancing the activity of fluconazole against Candida albicans
Abstract
Six new polyketides, peniexpansones A–F (1–4, 6, 7), along with an additional novel compound, (2E,4E,6E)-8-methyldeca-2,4,6-trienoic acid (5), and a known naphthopyrone, were isolated from the ethyl acetate extract of rice fermented with the soil fungus Penicillium expansum DWS880. Structurally, peniexpansones A–D feature a highly oxygenated tetrahydronaphthalene moiety linked to an acyl chain. The structures of the new compounds were elucidated by extensive spectroscopic analysis (1D/2D NMR and HRESIMS) and further supported by quantum chemical calculations. Peniexpansone C showed weak cytotoxicity against HCT116 cells (IC50 22.81 ± 0.42 μM) and moderate antimicrobial activity against pathogens including Staphylococcus aureus and Candida albicans. Moreover, peniexpansones A and B, can significantly improve the anti-fungal activity of fluconazole against resistant Candida albicans. Our results provided new structures for the future development of efficiency enhancement agents for anti-fungal drugs.Graphical Abstract

Keywords
Penicillium expansum Polyketides ECD calculation Anti-fungal activity enhancement1 Introduction
Fungi are well-known producers of diverse small-molecule natural products, such as terpenes, polyketides, alkaloids, and peptides [1, 2]. Among these, polyketides have attracted significant interest owing to their structural complexity and broad biological activities, establishing them as privileged scaffolds in drug discovery [3–5]. To date, over 10,000 natural polyketides have been identified, many of which possess substantial commercial and therapeutic value [6]. Prominent examples include lovastatin [7]—a leading therapeutic that comprises over 80% of the global lipid-lowering market—as well as the antibiotic erythromycin [8], the antifungal griseofulvin [9], and the anti-angiogenic compound fumagillin [10].
Since the initial report in 2009 of the antioxidative polyketide JBIR-12 [11]—characterized by a highly oxygenated tetrahydronaphthalene core linked to an acyl chain—fewer than ten natural analogues of this distinctive scaffold have been documented. In 2023, the absolute configuration of the tetrahydronaphthalene core was tentatively proposed via ECD comparison [12]; nevertheless, a more rigorous and experimentally supported stereochemical determination remained to be established.
In our ongoing investigation of bioactive metabolites from soil-derived fungi, Penicillium expansum DWS880 was isolated from a pristine humic layer collected in Daweishan National Forest Park. Chemical investigation of this strain led to the identification of six new polyketides, designated peniexpansones A–F. Structurally, peniexpansones A–D share the same rare architectural motif as JBIR-12. The complete and unequivocal absolute configuration for this structural family was established via 1D/2D NOE analysis and TDDFT-ECD calculations. All compounds were evaluated for cytotoxic and antimicrobial activities, providing a foundation for further exploration of their biological potential (Fig. 1).
The structures of compounds 1–8
2 Results and discussion
2.1 Structural analysis
Six new polyketides, peniexpansones A–F (1–4, 6, 7), were isolated from the ethyl acetate extract of the rice fermentation of the fungal strain. Structurally, compounds 1–4 each feature a 3-O-(2′E,4′E,6′E)-8′-methyldeca-2′,4′,6′-trienoate side chain with a chiral center at C-8′. Analysis of the NMR data indicated that the signals for this acyl side chain are highly consistent across compounds 1–4. Furthermore, a new (2E,4E,6E)-8-methyldeca-2,4,6-trienoic acid (5) was isolated as a yellow oil. The absolute configuration of 5 was established as 8S by comparative analysis of its chiroptical data with those of the (8R)-enantiomer [12]. The experimental ECD spectrum exhibited a negative Cotton effect at λmax 272 nm (Δε –5.7) and end absorption at 200 nm (Δε –7.8) (Table S81, Supporting Information), and a positive specific rotation was recorded ([α]D + 3.8)—both opposite in sign to the reported R enantiomer ([α]D –60, Cotton effects at 280 nm (Δε + 15) and 200 nm (Δε + 30) were both positive in the ECD spectrum). DFT calculations of specific rotation at the PBE1PBE/aug-cc-pvdz (IEFPCM, methanol) level further supported this assignment, yielding values of [α]D –124.3 for 8R and + 124.3 for 8S, which align with the respective literature and experimental signs. Accordingly, the side chain in 1–4 was determined to be 3-O-(2′E,4′E,6′E)-8′S-methyldeca-2′,4′,6′-trienoate. The structural differences among 1–4 arise solely from the substitution patterns on their tetrahydronaphthalene moieties.
Compound 1 was obtained as a yellow oil, and the molecular formula was established as C25H30O8 by HR-ESI–MS analysis based on the sodium adduct ion at m/z 481.1856 [M + Na]+ (calcd for C25H30O8Na+, 481.1833), indicating 11 degrees of unsaturation. The 1H NMR spectrum (Table 1) exhibited signals for six olefinic protons (δH 7.39, 6.55, 6.27, 6.12, 5.99, 5.86), one aromatic proton (δH 7.33), one methine proton (δH 2.15), two methylenes (δH 5.28, 4.98, 1.36), four methyl groups (δH 1.71, 1.57, 1.02, 0.87), and one methoxy group (δH 3.52). The 13C NMR and DEPT spectra revealed 25 carbon resonances, including two ketone carbonyls (δC 206.1, 191.8), one ester carbonyl (δC 165.4), six olefinic carbons (δC 147.6, 147.2, 142.9, 128.2, 127.6, 117.3), six aromatic carbons (δC 150.7, 144.2, 138.1, 122.1, 119.1, 111.5), two oxygenated non-protonated carbons (δC 83.3, 75.7), one methine (δC 38.8), two methylene groups (one oxygenated at δC 73.0), and five methyl groups (one oxygenated at δC 59.2). Analysis of COSY and HMBC correlations of 1 indicated two distinct substructures (Fig. 2): a highly oxygenated tetrahydronaphthalene moiety and a branched alkyl chain. The tetrahydronaphthalene core was confirmed by key HMBC correlations from H-9 to C-1/C-5/C-7/C-8, from H2-13 to C-5/C-6/C-7, from H3-11 to C-1/C-2/C-10, and from H3-12 to C-2/C-3/C-4. The 1D NMR data of 1 (Table 1) were highly similar to those of JBIR-12 [11], differing only by the replacement of a hydroxy group at C-13 in JBIR-12 with a methoxy group in 1, a modification supported by the observed HMBC correlation from H2-13 to C-14. The branched alkyl chain moiety was identical to that of JBIR-12, as evidenced by COSY correlations (H-10′/H-9′/H-8′/H-11′ and H-8′/H-7′/H-6′/H-5′/H-4′/H-3′/H-2′) and supporting HMBC data. The large coupling constants observed for H-2′/H-3′ (J = 15.3 Hz), H-4′/H-5′ (J = 14.8 Hz), and H-6′/H-7′ (J = 15.2 Hz) are consistent with trans-configured double bonds. This, combined with a key HMBC correlation from H-3′ to the ester carbonyl (C-1′), established the side chain as a (2E,4E,6E)-8-methyldeca-2,4,6-trienoate substituent.
The 1H NMR data (600 MHz) and 13C NMR data (150 MHz) of compounds 1–4 in CDCl3 (δ in ppm)
Key HMBC (H → C) and COSY (bold lines) correlations of compounds 1–7
The relative configuration of 1 was established through analysis of 2D NOESY and selective 1D NOE experiments. In alignment with the structural assignment of penijanthinone A [13], 2D NOESY correlations between H3-11 (δH 1.71) and H3-12 (δH 1.57) were absent. Furthermore, selective irradiation of H3-11 or H3-12 yielded no observable NOE enhancement of the other methyl group (Fig. S7, Supporting Information), indicating that H3-11 and H3-12 are oriented on opposite faces of ring A in 1. To determine the absolute configuration, the ECD spectra of 1R,3R,8′S-1a and 1S,3S,8′S-1a′ were calculated using the B3LYP/TZVP//B3LYP-D3(BJ)/TZVP level with the IEFPCM solvation model in methanol. As shown in Fig. 3, the calculated ECD curve for 1R,3R,8′S-1 closely matched the experimental spectrum. Thus, the absolute configuration of 1 was assigned as 1R,3R,8′S, and the compound was designated peniexpansone A.
The experimental and calculated ECD curves of 1 and 4
Compound 2 was obtained as a yellow oil, and the molecular formula was established as C23H26O8 by HR-ESI–MS at m/z 431.1721 [M + H]+ (calcd for C23H27O8+, 431.1700), indicating 11 degrees of unsaturation. Comparative analysis of the NMR data revealed that 1 and 2 share an identical core scaffold, differing structurally only at the C-6 position. The distinctive structural motif of three contiguous phenolic hydroxy groups at C-6, C-7, and C-8 in compound 2 was established by a definitive downfield shift of the C-6 resonance (δC 151.4) relative to that in 1 (δC 122.1), together with the key HMBC correlation from OH‑8 to C‑7/C‑8/C‑9. Due to signal overlap between H3-11 (δH 1.71) and H3-12 (δH 1.69), definitive NOE correlations for 2 could not be obtained. Attempts to determine the relative configuration via GIAO 13C NMR chemical shift calculations (using the STS protocol) [14] for the four possible diastereomers were unsuccessful; the computed chemical shifts showed insufficient dispersion and yielded Prel values below the 95% confidence threshold required for a reliable stereochemical assignment (Table S67). Consequently, TDDFT-ECD calculations were performed for the four diastereomeric candidates, 2a (1R,3R,8′S), 2a′ (1S,3S,8′S), 2b (1R,3S,8′S), and 2b′ (1S,3R,8′S) at the B3LYP/TZVP//B3LYP-D3(BJ)/TZVP level. The experimental ECD spectrum of 2 exhibited a positive Cotton effect (CE) at 200–300 nm and a negative Cotton effect at 300–350 nm. Among the computed spectra, only that of 2a reproduced both of these features, showing close agreement with the experimental data (Fig. 4). The calculated ECD curves for 2a and 2a′ (as well as for 2b and 2b′) displayed a nearly enantiomeric relationship, indicating that the configuration at C-8′ has a negligible influence on the overall ECD profile. Moreover, the calculated ECD curve of 2b showed negative CE at the region around 224 nm, which is the most significant discrepancy with the experimental curve. We further analyzed its dominant molecular orbitals contributing to the 224 nm band, and found that its negative CE is mainly caused by the transition from the n orbitals of oxygen atoms close to position 3 to the π* orbital of the 2-CO group (Fig. 5). While this transition is absent in 2a at the same band. Consequently, the configuration at C-3 should be the same as that in 2a, namely R. Thus, the absolute configuration of 2 was assigned as 1R,3R,8′S, and the compound was designated peniexpansone B.
The experimental and calculated ECD curves of 2 and 3
Orbital analysis of the characteristic CEs for structures 2a (3R) and 2b (3S). A The lowest-energy conformers of 2a and 2b, and their population. B Representative electronic transition orbitals in 2a contributing to the positive CE at 224.2 nm, and their percentage contributions. C Representative electronic transition orbitals in 2b contributing to the negative CE at 224.4 nm, and their percentage contributions
Compound 3 was obtained as a yellow oil. HR-ESI–MS analysis established its molecular formula as C24H26O8, based on the observed [M + H]+ ion at m/z 443.1713 (calcd for C24H27O8+, 443.1700), indicating 12 degrees of unsaturation. Comparative analysis of the 1H and 13C NMR data of 3 with those of 2 (Table 1) revealed an identical core scaffold. The key structural difference was the conversion of the C-6 substituent from a hydroxy group in 2 to an aldehyde group in 3, a change unequivocally supported by HMBC correlations from the aldehyde proton (CHO-13) to C-6 and C-7. Determination of the relative configuration of 3 was impeded by challenges similar to those with 2: signal overlap between H3-11 (δH 1.75) and H3-12 (δH 1.67) precluded definitive NOE correlations, and GIAO 13C NMR chemical shift calculations (STS protocol) [14] for the four possible diastereomers again exhibited insufficient dispersion (Prel < 95%; Table S68). Consequently, the absolute configuration was assigned by comparing the experimental ECD spectrum of 3 with TDDFT-ECD calculated curves. Geometry optimizations and subsequent ECD calculations were conducted at the B3LYP/TZVP//B3LYP-D3(BJ)/TZVP level in methanol for all candidate stereoisomers: 3a (1R,3R,8′S), 3a′ (1S,3S,8′S), 3b (1R,3S,8′S), and 3b′ (1S,3R,8′S). The experimental ECD spectrum of 3 displayed four major Cotton effects—positive, negative, positive, and negative—between 200 and 400 nm. Among all computed spectra, only that of diastereomer 3a (1R,3R,8′S) closely reproduced both the sign and magnitude of these effects (Fig. 4). Structures 3a (3R) and 3b (3S) are epimers at postion 3, and show opposite CE at the band around 278 nm, indicating that the chiral center at position 3 plays major role in their different optical features. We further analyzed the orbital contributions within this absorption band, and found that the transitions are dominated by excitation from the n orbitals of the oxygen atoms on the core skeleton to the π* orbital of the side chain, with the chiral center at position 3 serving as a pivotal hub (Fig. 6). The negative CE of the structure 3a in this band matches the experimental value, confirming that the configuration at position 3 is R. Consequently, the absolute configuration of 3 was assigned as 1R,3R,8′S, and the compound was designated peniexpansone C.
Orbital analysis of the characteristic CEs for structures 3a (3R) and 3b (3S). A The lowest-energy conformers of 3a and 3b, and their population. B Representative electronic transition orbitals in 3a contributing to the negative CE at 277.9 nm, and their percentage contributions. C Representative electronic transition orbitals in 3b contributing to the positive CE at 277.7 nm, and their percentage contributions
Compound 4 was obtained as a yellow oil, and the molecular formula was established as C24H28O6 by HR-ESI–MS analysis based on the detection of an [M + H]+ ion at m/z 413.1976 (calcd for C24H29O6+, 413.1959), indicating 11 degrees of unsaturation. The NMR spectra of 4 closely resembled those of 2, indicating an identical core structure with two key modifications: the hydroxy group at C-6 was supplanted with a methyl group (H3-13, δH 2.52), and the hydroxy substituent at C-7 was replaced by an aromatic proton (H-7, δH 6.70). The former modification was confirmed by an HMBC correlation from H3-13 to C-4/C-5/C-6, and the latter by correlations from H-7 to C-5/C-8/C-13 (Fig. 2). The relative configuration of 4 was determined through selective 1D NOE experiments, applying the strategy previously used for 1. Irradiation of H3-11 did not yield an observable enhancement of H3-12 (Fig. S35), indicating that these methyl groups are oriented on opposite faces of ring A. The absolute configuration of compound 4 was assigned as 1R,3R,8′S by comparing the experimental ECD spectrum with the TDDFT-calculated ECD curves of the two candidate diastereomers, 4a (1R,3R,8′S) and 4a′ (1S,3S,8′S) (Fig. 3). Thus, the complete structure of 4 was elucidated and assigned the name peniexpansone D. The consistent assignment of the 1R,3R,8′S configuration across these compounds — based on ECD analysis —suggests a highly conserved biosynthetic pathway, indicative of a common genetic origin within this fungal source.
Compounds 6 and 7 were identified as a pair of epimers, both exhibiting the molecular formula C12H14O3 as established by HR-ESI–MS analysis, with sodium adduct ions observed at m/z 229.0835 [M + Na]+ for compound 6 and m/z 229.0809 [M + Na]+ for compound 7 (calcd for C12H14O3Na+, 229.0835), consistent with 6 degrees of unsaturation. The similarity in their 1H and 13C NMR spectra (Table 2) suggested that both compounds share the same carbon skeleton. The key feature in the 1H NMR spectra of both compounds was a set of aromatic protons indicative of a 1,2,3-trisubstituted benzene ring. Additional signals were consistent with an oxygenated methine, a methylene group, and two methyl groups. The 13C NMR data confirmed the presence of 12 carbon resonances, which were classified as one ketone carbonyl, one oxygenated non-protonated carbon, one oxygenated methine, one methylene group, two methyl groups, and six aromatic carbons. COSY correlations revealed two spin systems: H2-4/H-3 and H-6/H-7/H-8. The planar structures were further elucidated by key HMBC correlations, observed from H-8 to C-6 and C-10; from H-7 to C-5 and C-9; from H-6 to C-4 and C-10; from H3-11 to C-8 and C-10; from H3-12 to C-1 and C-3; from H-4 to C-2 and C-10; and from H-3 to C-1 and C-5 (Fig. 2). On the basis of these results, compounds 6 and 7 were deduced to be tetralone derivatives [15].
The 1H NMR data (600 MHz) and 13C NMR data (150 MHz) of compounds 6 and 7 in CDCl3 (δ in ppm)
The relative configuration of 6 was established by a key NOESY correlation between H-3 and H3-12, indicating their co-facial orientation on the cyclohexanone ring (Fig. 7). For 7, the large coupling constant J3,4 (11.3 Hz) between H-4b (δH 3.02) and H-3 indicated a trans-diaxial relationship, while a NOESY correlation between H-4b and H3-12 placed H-3 and the methyl group on opposite faces. The absolute configurations of the cyclohexanone moiety in 6 and 7 were assigned as 2S,3S and 2S,3R by comparison of experimental and calculated ECD spectra (Fig. 8); thus, the compounds were designated peniexpansones E and F, respectively.
Key NOESY correlations of compounds 6 and 7
The experimental and calculated ECD curves of 6 and 7
By comparing spectroscopic properties with the reported data [16], the structure of 8 was elucidated as ligustrone B, first isolated from Cercospora ligustrina in 1975 [17].
2.2 Cytotoxicity and anti-microbial assays
Among the tested compounds, only 3 exhibited notable antimicrobial activity. It showed potent effects against Staphylococcus aureus (MIC = 25.0 μM), comparable to streptomycin (MIC = 22.0 μM), and inhibited Acinetobacter baumannii, Saccharomyces cerevisiae, and Candida albicans (MIC = 50.0 μM each). Notably, 3 was significantly more active than terbinafine (MIC = 439.2 μM) and fluconazole (MIC > 4 mM) against resistant C. albicans. None of the other compounds (1, 2, 4–8) were active (MIC > 50 μM) against the microbial strains tested.
Strikingly, despite their lack of intrinsic antifungal activity, compounds 1 and 2 demonstrated a profound ability to potentiate fluconazole against the resistant C. albicans isolate. The potentiation effect was concentration-dependent and remarkably strong for compound 1, which at 40 and 20 μM reduced the fluconazole MIC from > 4 mM to 10 μM and 25 μM, respectively. Similarly, compound 2 at 40 μM lowered the fluconazole MIC to 100 μM. In contrast, compounds 3 and 4 showed only weak potentiation, reducing the fluconazole MIC to 500 μM and 800 μM at 40 μM, respectively (Fig. S85-90, Supporting Information).
In cytotoxicity evaluations against HCT116, 4T1, and MHCC97H cancer cell lines, only compound 3 exhibited weak activity against HCT116 cells, with an IC50 value of 22.81 ± 0.42 μM. The other compounds showed no significant cytotoxicity across all tested cancer lines.
The selectivity and safety profiles of these compounds were further assessed using the human gastric epithelial cell line GES-1. Notably, compounds 1 and 2 displayed no cytotoxicity toward GES-1 cells at a concentration of 50 μM, underscoring their favorable safety profile as potential fluconazole potentiators. In contrast, compound 3 demonstrated moderate cytotoxicity against GES-1 cells.
3 Experimental
3.1 General experiment procedures
UV and CD spectra were acquired on a Chirascan™-plus circular dichroism spectrometer. Optical rotations were measured using a Rudolph Research Analytical Autopol Ⅳ automatic polarimeter. IR spectra were recorded on a Shimadzu Fourier transform infrared spectrometer with KBr pellets. High-resolution electrospray ionization mass spectrometry (HR-ESI–MS) spectra were acquired with an Agilent 6500 series Q-TOF mass spectrometer (Agilent Technologies, Singapore). NMR spectra were recorded at 600 MHz for 1H and 150 MHz for 13C on Bruker spectrometers (Bruker BioSpin, Rheinstetten, Germany). Column chromatography was performed using macroporous adsorbent resin D101 (Tianjin Haoju Resin Technology Co., LTD., Tianjin, China), silica gel (200 − 300 mesh, Qingdao Marine Chemical, Qingdao, China), and Sephadex LH-20 (GE Healthcare, Uppsala, Sweden). HPLC separations were conducted using an Agilent 1260 system, which included an automatic sampler, a pump, and a wavelength absorbance detector, fitted with an Agilent semipreparative column (9.4 mm × 150 mm, C18, 5 μm) (Agilent Technologies). Additionally, an EClassical P3500 prep-HPLC system equipped with a Sinochrom ODS-BP column (Elite, 5 μm, 10 mm × 250 mm), Phenyl-Hexyl (Shimadzu, 5 μm, 10 mm × 250 mm) and a DAD detector (EClassical P3500, Dalian Elite Analytical Instruments Co., Ltd., Dalian, China) was used. Thin-layer chromatography (TLC) separations were performed on precoated silica gel GF254 plates (Qingdao Marine Chemical), with spots visualized under UV light (254 or 365 nm) or by spraying with 10% H2SO4 in ethanol followed by heating.
3.2 Fungal material
The fungal strain DWS880 was isolated from a soil sample collected from Daweishan Forest Park, Hunan Province, China, in 2022. This strain was identified as P. expansum based on ITS sequencing, with a percentage identity of 99.50% to be P. expansum for ITS (Table S84). The fungal specimen is deposited at the Xiangya School of Pharmaceutical Sciences, Central South University, Changsha, China.
3.3 Extraction and isolation
The fungal strain was initially cultivated on potato dextrose agar for 8 days at 25 ℃. The mycelia were then cut into pieces and cultured on a cooked rice medium, prepared by mixing 100 g of rice with 80 mL of water in 500 mL Erlenmeyer flasks. The total weight of the rice used was 5 kg. The rice cultures were maintained at 25 ℃ for 71 d. After incubation at 25 ℃ for 71 d, the rice was chopped and exhaustively extracted with EtOAc (15 L, five cycles, room temperature). The solvent was then evaporated under vacuum to yield 90 g of crude extract. The crude extract was adsorbed onto D101 macroporous resin and eluted stepwise with water, 20% methanol, 80% methanol, and 100% methanol. The 80% aquatic MeOH portion was collected and condensed to yield 60 g residue. The 80% MeOH portion was subjected to silica gel column flushing with gradient petroleum ether and ethyl acetate (50:1 → 1:1, v/v) to afford fractions A–G. Fraction D was separated by Sephadex LH-20 (MeOH) to afford three sub-fractions (Fr. D1–Fr. D3). Fr. D2 was purified by semipreparative HPLC (MeOH–H2O, 60:40 → 82:18, v/v, 27 min, 3.0 mL/min, Phenyl-Hexyl column, 10 × 250 mm, 5 μm) to yield six subfractions (Fr. D21–Fr. D26). Fr. D21 was further purified by semipreparative HPLC ((MeOH–H2O, 43:57, v/v, 30 min, 3.0 mL/min, Phenyl-Hexyl, 10 mm × 250 mm, 5 μm) to yield compound 6 (1.2 mg, tR = 21.89 min) and 7 (2.0 mg, tR = 24.84 min). Fr. F was separated by Sephadex LH-20 (MeOH) to afford four sub-fractions (Fr. F1 − Fr. F4). Fr. F2 was purified by semipreparative HPLC (MeCN–H2O, 35:65 → 75:25, v/v, 35 min, 3.0 mL/min, ZORBAX-C18, 9.4 mm × 150 mm, 5 μm) to yield compound 5 (7.0 mg, tR = 20.55 min), 2 (6.0 mg, tR = 23.13 min), 3 (7.0 mg, tR = 25.54 min), and 4 (3.0 mg, tR = 28.61 min). Fr. F3 was purified by semipreparative HPLC (MeCN–H2O, 50:50 → 75:25, v/v, 25 min, 3.0 mL/min, Sinochrom ODS-BP, 10 mm × 250 mm, 5 μm) to yield compound 8 (10.0 mg, tR = 13.46 min). Fr. G was separated by Sephadex LH-20 (MeOH) to afford four sub-fractions (Fr. G1 − Fr. G4). Fr. G2 was further purified by semipreparative HPLC (MeCN–H2O, 50:50 → 71:29, v/v, 25 min, 3.0 mL/min, Sinochrom ODS-BP, 10 mm × 250 mm, 5 μm) to yield compound 1 (6.0 mg, tR = 23.75 min).
3.4 Compound characterization
3.4.1 Peniexpansones A (1)
Yellow oil. [α]D25 −71.2 (c 0.13, MeOH). UV (MeOH) λmax (log ε) 300 (3.17) nm. IR (KBr): υmax 3484, 2959, 1741, 1684, 1596, 1449, 1371, 1266, 1084, 739 cm−1. 1H NMR (600 MHz in CDCl3) and 13C NMR (150 MHz in CDCl3) data, see Table 1. HRESIMS m/z 481.1856 [M + Na]+ (calcd for C25H30O8Na+, 481.1833).
3.4.2 Peniexpansones B (2)
Yellow oil. [α]D25 −10.47 (c 0.09, MeOH). UV (MeOH) λmax (log ε) 310 (3.51) nm. IR (KBr): υmax 3434, 1704, 1614, 1302, 1094, 894 cm−1. 1H NMR (600 MHz in CDCl3) and 13C NMR (150 MHz in CDCl3) data, see Table 1. HRESIMS m/z 431.1721 [M + H]+ (calcd for C23H27O8+, 431.1700).
3.4.3 Peniexpansones C (3)
Yellow oil. [α]D25 −57.61 (c 0.09, MeOH). UV (MeOH) λmax (log ε) 250 (3.05) nm, 315 (3.08) nm. IR (KBr): υmax 3436, 2964, 1697, 1461, 1291, 1082, 1008, 890, 740 cm−1. 1H NMR (600 MHz in CDCl3) and 13C NMR (150 MHz in CDCl3) data, see Table 1. HRESIMS m/z 443.1713 [M + H]+ (calcd for C24H27O8+, 443.1700).
3.4.4 Peniexpansones D (4)
Yellow oil. [α]D25 −79.97 (c 0.57, MeOH). UV (MeOH) λmax (log ε) 288 (3.34) nm. IR (KBr): υmax 3428, 2966, 1692, 1604, 1464, 1320, 1096, 867 cm−1. 1H NMR (600 MHz in CDCl3) and 13C NMR (150 MHz in CDCl3) data, see Table 1. HRESIMS m/z 413.1976 [M + H]+ (calcd for C24H29O6+, 413.1959).
3.4.5 (2E,4E,6E)-8-methyldeca-2,4,6-trienoic acid (5)
Yellow oil. [α]D25 +3.77 (c 0.35, MeOH). UV (MeOH) λmax (log ε) 258 (2.72) nm. IR (KBr): υmax 3434, 2960, 1636, 1384, 498 cm−1. 1H NMR (600 MHz in CDCl3) and 13C NMR (150 MHz in CDCl3) data. HRESIMS m/z 181.1223 [M + H]+ (calcd for C11H17O2+, 181.1223). 1H NMR (600 MHz, CDCl3) δ 7.39 (dd, 1H, H-3), 6.58 (dd, 1H, H-5), 6.26 (dd, 1H, H-4), 6.12 (dd, 1H, H-6) 5.87, (ov, 1H, H-7), 5.85 (ov, 1H, H-2), 2.14 (m, 1H, H-8), 1.36 (m, 2H, H-9), 1.02 (d, 3H, Me-11) 0.87 (t, 3H, Me-10); 13C NMR (150 MHz, CDCl3) δ 172.4 (C-1), 147.4 (C-3), 147.2 (C-7), 142.7 (C-5), 128.1 (C-6), 127.6 (C-4), 118.8 (C-2), 38.8 (C-8), 29.5 (C-9), 19.7 (C-11), 11.7 (C-10).
3.4.6 Peniexpansones E (6)
Yellow oil. [α]D25 +93.33 (c 0.06, MeOH). UV (MeOH) λmax (log ε) 210 (3.44) nm. 1H NMR (600 MHz in CDCl3) and 13C NMR (150 MHz in CDCl3) data, see Table 2. HRESIMS m/z 229.0835 [M + Na]+ (calcd for C12H14O3Na+, 229.0835).
3.4.7 Peniexpansones F (7)
Yellow oil. [α]D25 +25.37 (c 0.07, MeOH). UV (MeOH) λmax (log ε) 210 (3.17) nm. 1H NMR (600 MHz in CDCl3) and 13C NMR (150 MHz in CDCl3) data, see Table 2. HRESIMS m/z 229.0809 [M + Na]+ (calcd for C12H14O3Na+, 229.0835).
3.5 Quantum chemical calculation
Crest [19] was used to search the conformational space of candidate structures at the GFNFF [20] level of theory, followed by optimization at the GFN2-xTB [21] level with a 4 kcal mol⁻1 energy window to discard high-energy conformers. Each conformer was then optimized and subjected to frequency analysis at the B3LYP-D3(BJ)/TZVP(IEFPCM) level of theory. ECD calculations for compounds 1–4, 6 and 7 were carried out at the B3LYP/TZVP (IEFPCM, methanol) level. DFT-GIAO 13C NMR shielding tensors were computed at the ωB97X-D/6-31G*(IEFPCM, chloroform) level, consistent with the reported STS protocol [14]. The calculated shielding tensors and ECD curves of all conformers were Boltzmann-averaged based on Gibbs free energies. The ECD spectra of compounds 1, 2, 6, and 7 were simulated using SpecDis v1.71 [22] with a sigma/gamma value of 0.35 eV, while values of 0.28 eV and 0.45 eV were applied for compounds 3 and 4, respectively. All DFT calculations were performed using the Gaussian 16 [23] software package. The DFT optimized geometry data, relative energies, and conformational populations of all optimized structures are provided in the Supporting Information. The 3D structures for compounds 6 and 7 were rendered with VMD 1.9.3 [24]. Molecular orbitals were computed with Multiwfn 3.8 (dev) [25] and subsequently visualized using VMD 1.9.3 [24].
3.6 Assay of antimicrobial activity
The antimicrobial activities of all compounds were evaluated against a panel of bacterial and fungal strains, including Bacillus subtilis ATCC 6051, Escherichia coli DH5α, Staphylococcus aureus ATCC 25923, Acinetobacter baumannii ATCC 19606, Pseudomonas aeruginosa ATCC 27853, Saccharomyces cerevisiae ATCC 9763, and Candida albicans ATCC 10231. Each sample was dissolved in 80% aqueous ethanol to achieve an initial concentration of 2 mM, followed by dilution to a final concentration of 7.5 × 105 CFU/mL. Bacterial suspensions were dispensed into individual wells of a 96-well culture plate. The plate was incubated at 37 ℃ for 48 h, after which microbial growth inhibition was assessed visually based on the absence of turbidity. The experimental design included medium blanks, bacterial controls, streptomycin as a positive control for bacteria, and terbinafine and fluconazole as positive controls for fungi.
The ability of the compounds to potentiate fluconazole against the resistant C. albicans was assessed using a broth microdilution method. Fixed concentrations of fluconazole were combined with serial dilutions of the test compounds in a 96-well plate. The fungal suspension was prepared and added as described previously. After incubation at 37 ℃ for 48 h, the minimum inhibitory concentration (MIC) of the combinations was determined visually by the absence of turbidity. The effective combinations are reported in Table 3.
Compound-mediated potentiation of fluconazole against resistant C. albicans
3.7 Cytotoxicity assay
The cytotoxicity was evaluated against the human gastric epithelial cell line GES-1, human colorectal cancer cell line HCT116, mouse breast cancer cell line 4T1, and highly metastatic liver cancer cell line MHCC97H. The assay was performed using the CCK-8 method, with doxorubicin serving as a positive control. Cells were cultured in DMEM supplemented with 10% fetal bovine serum and 1% penicillin–streptomycin, and maintained in a humidified atmosphere at 37 ℃ with 5% CO2. For the assay, cells were seeded into 96-well plates at a density of 10,000 cells per well and incubated for 24 h. Subsequently, the cells were exposed to the test compounds at varying concentrations ≤ 50 μM, prepared in culture medium, for a duration of 48 h. Alongside the experimental groups, blank controls (wells with CCK-8 solution but no cells) and negative controls (wells with vehicle) were also included. After the treatment, 10 μL of CCK-8 reagent was added to each well and incubated for 1 h. The optical density was then measured at 450 nm using a microplate reader (Feyond-A300, allsheng, Hangzhou, China). The IC50 values, indicative of the concentration of the compounds required to inhibit cell growth by 50%, were calculated using the GraphPad Prism 10.6.0 software.
4 Conclusions
In conclusion, the chemical investigation of P. expansum DWS880 led to the identification of six new polyketides, designated peniexpansones A–F. Among them, peniexpansones A–D (1–4) possess a highly oxygenated tetrahydronaphthalene unit connected to an acyl chain bearing a stereogenic carbon at C-8′, representing a unique substitution pattern. In biological evaluations, compound 3 exhibited weak cytotoxicity against the HCT116 cell line. Notably, it also demonstrated antimicrobial activity against Staphylococcus aureus, Acinetobacter baumannii, Saccharomyces cerevisiae, and Candida albicans. More importantly, compounds 1 and 2 showed no intrinsic antifungal activity but functioned as strong and selective potentiators of fluconazole against resistant C. albicans, while exhibiting minimal cytotoxicity against human normal cell lines, marking them as exceptionally promising candidates for combination therapy. These findings enhance the structural and bioactive diversity known among fungal polyketides.
Notes
Acknowledgements
The authors thank the Institute for Advanced Study, Central South University for NMR measurement.
Author contributions
Wen-Yu Lu: Writing–original draft, Methodology, Investigation. Qing-Hui Xiao: Investigation. Ai-Lin Liang: Methodology. Peng-Ju Xu: Methodology. Jing Li: Writing–review & editing, Investigation, Funding acquisition. Wen-Xuan Wang: Writing–review & editing, Writing–original draft, Funding acquisition, Conceptualization.
Funding
This work was supported by National Natural Science Foundation of China (82173713), and Natural Science Foundation of Hunan Province (2024JJ8127). Project of Hunan Administration of Traditional Chinese Medicine (B2023059), Changsha Municipal Natural Science Foundation (kq2403041). Fundamental Research Funds for the Central Universities of Central South University (No. 2024ZZTS0990).
Data availability
All data generated or analyzed during this study are included in this published article and its supplementary information files.
Declarations
Competing interests
The authors declare no conflict of interest in this manuscript.
References
-
1.Lv F, Zeng Y. Novel bioactive natural products from marine-derived Penicillium fungi: a review (2021–2023). Mar Drugs. 2024;22(5): 191. CrossRef PubMed Google Scholar
-
2.Newman DJ, Cragg GM. Natural products as sources of new drugs over the nearly four decades from 01/1981 to 09/2019. J Nat Prod. 2020;83(3): 770-803. CrossRef PubMed Google Scholar
-
3.Guo LD, Wu Y, Xu X, Lin Z, Tong R. Bent π-conjugation within a macrocycle: asymmetric total syntheses of Spirohexenolides A and B. Angew Chem Int Ed Engl. 2024;63(1): e202316259. CrossRef PubMed Google Scholar
-
4.Hao X, Wang ZJ, Shi X, Liu Y, Yan P, Yang Q, et al. Genome mining and targeted discovery of Twistane-Containing Angucyclinones from Marine Streptomyces phaeochromogenes XSA2. Org Lett. 2025;27(12): 2987-91. CrossRef PubMed Google Scholar
-
5.Tan Q, Fan RZ, Yang W, Zou G, Chen T, Wu J, et al. (+)/(−)-mycosphatide A, a pair of highly oxidized polyketides with lipid-lowering activity from the mangrove endophytic fungus Mycosphaerella sp. SYSU-DZG01. Chin Chem Lett. 2024;35(9): 109390. CrossRef PubMed Google Scholar
-
6.Zhang J, Jin SX, Yun Q, Qu XD. Biosynthesis of the unnatural extender units with polyketides and their structural modifications for applications in medicines. Synth Biol J. 2024;5(3): 561-70. PubMed Google Scholar
-
7.Zeller FP, Uvodich KC. Lovastatin for hypercholesterolemia. Drug Intell Clin Pharm. 1988;22(7–8): 542-5. PubMed Google Scholar
-
8.Brittain DC. Erythromycin. Med Clin North Am. 1987;71(6): 1147-54. CrossRef PubMed Google Scholar
-
9.De Carli L, Larizza L. Griseofulvin. Mutat Res. 1988;195(2): 91-126. CrossRef PubMed Google Scholar
-
10.Lefkove B, Govindarajan B, Arbiser JL. Fumagillin: an anti-infective as a parent molecule for novel angiogenesis inhibitors. Expert Rev Anti Infect Ther. 2007;5(4): 573-9. CrossRef PubMed Google Scholar
-
11.Izumikawa M, Nagai A, Doi T, Takagi M, Shin-Ya K. JBIR-12, a novel antioxidative agent from Penicillium sp. NBRC 103941. J Antibiot. 2009;62(4): 177-80. CrossRef PubMed Google Scholar
-
12.Wang Z, Li Z, Niu C, Yang L, Zhai Y, Li D, et al. Aniline-induced production of aniline-containing polyketides and related bicyclic polyketides by the Yellow River wetland-derived fungus Talaromyces funiculosus. Front Microbiol. 2023;14: 1200680. CrossRef PubMed Google Scholar
-
13.Chen M, Zheng YY, Chen ZQ, Shen NX, Shen L, Zhang FM, et al. NaBr-induced production of brominated azaphilones and related tricyclic polyketides by the marine-derived fungus Penicillium janthinellum HK1-6. J Nat Prod. 2019;82(2): 368-74. CrossRef PubMed Google Scholar
-
14.Li J, Liu JK, Wang WX. GIAO 13C NMR calculation with sorted training sets improves accuracy and reliability for structural assignation. J Org Chem. 2020;85(17): 11350-8. CrossRef PubMed Google Scholar
-
15.Liao HX, Zheng CJ, Huang GL, Mei RQ, Nong XH, Shao TM, et al. Bioactive polyketide derivatives from the mangrove-derived fungus Daldinia eschscholtzii HJ004. J Nat Prod. 2019;82(8): 2211-9. CrossRef PubMed Google Scholar
-
16.Liu T, Cai T, Huo J, Liu H, Li A, Yin M, et al. Force-enhanced sensitive and specific detection of DNA-intercalative agents directly from microorganisms at single-molecule level. Nucleic Acids Res. 2024;52(18): e86. CrossRef PubMed Google Scholar
-
17.Arnone A, Camarda L, Merlini L, Nasini G. The structure of ligustrone A, B and C, new metabolites from Cercospora ligustrina Boerema. Gazz Chim Ital. 1975;105(105): 1093-103. PubMed Google Scholar
-
18.Wang P, Fan RQ, Bai CY. Research on synergism of Xianglian lotion to Fluconazole against drug-resistance Candida albicans based on RNA-seq technology. CJGMCM. 2022;37(09): 1581-4. PubMed Google Scholar
-
19.Pracht P, Bohle F, Grimme S. Automated exploration of the low-energy chemical space with fast quantum chemical methods. Phys Chem Chem Phys. 2020;22(14): 7169-92. CrossRef PubMed Google Scholar
-
20.Spicher S, Grimme S. Robust atomistic modeling of materials, organometallic, and biochemical systems. Angew Chem Int Ed Engl. 2020;59(36): 15665-73. CrossRef PubMed Google Scholar
-
21.Bannwarth C, Ehlert S, Grimme S. GFN2-xTB—an accurate and broadly parametrized self-consistent tight-binding quantum chemical method with multipole electrostatics and density-dependent dispersion contributions. J Chem Theory Comput. 2019;15(3): 1652-71. CrossRef PubMed Google Scholar
-
22.Bruhn T, Schaumlöffel A, Hemberger Y, Bringmann G. SpecDis: quantifying the comparison of calculated and experimental electronic circular dichroism spectra. Chirality. 2013;25(4): 243-9. CrossRef PubMed Google Scholar
-
23.Frisch MJ TG, Schlegel HB, Scuseria GE, Robb MA, Cheeseman JR, Scalmani G, et al. Gaussian 16, revision B.01. Wallingford: Gaussian, Inc. 2016g. PubMed Google Scholar
-
24.Humphrey WF, Dalke A, Schulten K. VMD: visual molecular dynamics. J Mol Graph. 1996;14(1): 33-8. CrossRef PubMed Google Scholar
-
25.Lu T, Chen FW. Multiwfn: a multifunctional wavefunction analyzer. J Comput Chem. 2012;33(5): 580-92. CrossRef PubMed Google Scholar
Copyright information
© The Author(s) 2026
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/.









