Verrucosidin Derivatives from the Deep Sea Cold-Seep-Derived Fungus Penicillium polonicum CS-252

Six novel verrucosidin derivatives, namely, poloncosidins A–F (1–6), together with one known analogue (7), were isolated and identified from the deep-sea-derived fungus Penicillium polonicum CS-252, which was obtained from cold-seep sediments collected in the South China Sea at a depth of 1183 m. Their structures were mainly established on the basis of a detailed interpretation of NMR spectroscopic and mass spectrometric data. The relative and absolute configurations of compounds 1–6 were determined by ECD calculations and a DP4+ probability analysis. Compounds 1–5 represent the first examples of verrucosidins with a 2,5-dihydrofuran ring which is uncommon among the known analogues. These compounds exhibited inhibitory activities against several human and aquatic pathogens with MIC values ranging from 4 to 32 μg/mL.


Introduction
Verrucosidins belong to a family of highly reduced polyketides, generally sharing a methylated α-pyrone, a conjugated polyene linker, and an epoxidated tetrahydrofuran ring [1]. All verrucosidin derivatives previously isolated and characterized from the fungal genus Penicillium were found to possess bioactivities such as antibacterial and GRP78/BiP down-regulator activities [2,3]. Penicillium polonicum is a ubiquitous fungus found as saprophyte in diverse environments [4], from which bioactive molecules with diverse chemical structures have been reported including diketopiperazine, α-pyrone derivative, dimeric anthraquinone, and taxol [5][6][7][8][9][10]. The verrucosidin biosynthesis gene cluster was also reported in P. polonicum, which confirmed the ability of the fungus to produce related secondary metabolites [11,12].
In our ongoing research for the discovery of bioactive secondary metabolites from marine-derived fungi [13][14][15], the fungal strain P. polonicum CS-252, which was obtained from sediments collected at the deep-sea cold seep area in the South China Sea, attracted our attention due to its unique HPLC profiles which showed peaks with different UV absorptions from those in our HPLC database. A chemical investigation of the culture extract of P. polonicum CS-252 led to the characterizations of six new verrucosidin derivatives, designated as poloncosidins A-F (1-6, Figure 1), together with a known analogue, deoxyverrucosidin (7). Among them, compounds 1-5 represent the first examples of verrucosidins with a 2,5-dihydrofuran ring. Details of the isolation, purification, structure elucidation, and biological evaluation of compounds 1-7 are provided herein.

Results and Discussion
The fungus P. polonicum CS-252 was cultured on a rice solid medium for chemical investigation. The obtained culture was extracted thoroughly with EtOAc to produce an organic extract that was further fractionated and purified with a combination of column chromatography over Lobar LiChroprep RP-18, Si gel, and Sephadex LH-20 as well as by prep. TLC to yield compounds 1-7.

Results and Discussion
The fungus P. polonicum CS-252 was cultured on a rice solid medium for chemical investigation. The obtained culture was extracted thoroughly with EtOAc to produce an organic extract that was further fractionated and purified with a combination of column chromatography over Lobar LiChroprep RP-18, Si gel, and Sephadex LH-20 as well as by prep. TLC to yield compounds 1-7.
Poloncosidin A (1) was obtained as a colorless oil and its molecular formula was assigned as C 24 H 32 O 4 on the basis of HRESIMS data, corresponding to nine indices of hydrogen deficiency. The 1 H-NMR and HSQC spectra of 1 (Table 1 and Figure S5) revealed the signals of nine methyls with one doublet (δ H 1. 16 Figure S3) indicated all 24 carbon resonances including nine methyls, five methines (with four olefinic at δ C 127.2, 137.1, 137.2, and 139.1 as well as an oxygenated at δ C 81.6), and ten nonprotonated carbons (with eight olefinic, one carbonyl, and one oxygenated). The aforementioned spectroscopic data suggested the characteristics of a verrucosidin derivative, which are similar to those of deoxyverrucosidin (7) previously described in Penicillium sp. in 2005, with the difference that the epoxide moiety in 7 was replaced by a double bond in 1. This difference was supported by the NMR evidence that signals for an oxygenated methine (δ H 3.45; δ C 67.5, CH-13) and an oxygenated nonprotonated carbon (δ C 67.4, C-14) in the NMR spectra of deoxyverrucosidin (7) [3] were not present in those of 1, while resonances for an olefinic methine (δ H 5.61, H-13; δ C 127.2, C-13) and an olefinic nonprotonated carbon (δ C 138.2, C-14) were observed in the NMR spectra of 1 (Tables 1 and 2). These data indicated that the epoxidated tetrahydrofuran ring in deoxyverrucosidin was replaced by a 2,5-dihydrofuran ring in 1. This deduction was supported by the key HMBC from H-11 (δ H 5.53) and H 3 -21 (δ H 1.33) to C-13 (δ C 127.2) and from H 3 -22 (δ H 1.66) to C-13, C-14 (δ C 138.2), and C-15 (δ C 81.6) ( Figure 2). The planar structure of 1 was thus identified as the first example of verrucosidins possessing a 2,5-dihydrofuran moiety. The relative configuration of compound 1 was established by an analysis of its NOESY data (Figure 3). Key NOESY enhancement from H-15 to H3-21 indicated the cofacial orientation of these protons, which established the relative configuration of the 2,5-dihydrofuran moiety. The geometries of the double bonds at C-6, C-8, and C-10 were determined to be E configurations by the NOESY correlations from H-7 to H-9 and H3-17   The relative configuration of compound 1 was established by an analysis of its NOESY data ( Figure 3). Key NOESY enhancement from H-15 to H 3 -21 indicated the cofacial orientation of these protons, which established the relative configuration of the 2,5-dihydrofuran moiety. The geometries of the double bonds at C-6, C-8, and C-10 were determined to be E configurations by the NOESY correlations from H-7 to H-9 and H 3 -17 and from H-9 to H-11 ( Figure 3). To establish the absolute configuration of 1, the time-dependent density functional (TDDFT)-ECD calculations at the BH&HLYP/TZVP level were performed in Gaussian 09. The calculated ECD spectrum for the (12S,15R)-isomer of compound 1 matched well with the experimental curve, leading to the establishment of the absolute configuration of 1 as 12S, 15R ( Figure 4). 2,5-dihydrofuran moiety. The geometries of the double bonds at C-6, C-8, and C-10 were determined to be E configurations by the NOESY correlations from H-7 to H-9 and H3-17 and from H-9 to H-11 ( Figure 3). To establish the absolute configuration of 1, the timedependent density functional (TDDFT)-ECD calculations at the BH&HLYP/TZVP level were performed in Gaussian 09. The calculated ECD spectrum for the (12S,15R)-isomer of compound 1 matched well with the experimental curve, leading to the establishment of the absolute configuration of 1 as 12S, 15R ( Figure 4).   Poloncosidin B (2) was obtained as a colorless oil and the molecular formula was assigned as C23H30O4 based on an analysis of the HRESIMS data with the same hydrogen deficiency index as that of 1. A comprehensive analysis of the 1 H-and 13 C-NMR data of 2 indicated a similar structure to that of 1, only differing in the absence of a substituent. Signals for a methyl (δH 1.93; δC 9.9, CH3-16) and an olefinic nonprotonated carbon (δC 108.6, C-2) in the NMR spectra of 1 were found to be absent in those of 2 (Tables 1 and 2). In contrast, compound 2 exhibited an additional resonance at δH 5.61 (1H, s, H-2) in the 1 H-NMR spectrum that correlates with an olefinic methine carbon at δC 87.9 (C-2), accounting for the shielding effect of the carbonyl group at C-1. The key HMBC correlations from H-2 to C-3 (δC 170.7) and C-4 (δC 106.0) as well as other key HMBC and COSY correlations ( Figure 2) established the structure of compound 2 as shown.
The relative configuration of the 2,5-dihydrofuran ring and the geometries of the double bonds at C-6, C-8, and C-10 were determined to be the same as those of compound 1 with an analysis of the NOESY correlations ( Figure 3). The absolute configuration of 2 was assigned as 12S and 15R, which was evidenced by a comparison of the overlaid experimental and calculated ECD curves at the BH&HLYP/TZVP level ( Figure 5). Poloncosidin B (2) was obtained as a colorless oil and the molecular formula was assigned as C 23 H 30 O 4 based on an analysis of the HRESIMS data with the same hydrogen deficiency index as that of 1. A comprehensive analysis of the 1 H-and 13 C-NMR data of 2 indicated a similar structure to that of 1, only differing in the absence of a substituent. Signals for a methyl (δ H 1.93; δ C 9.9, CH 3 -16) and an olefinic nonprotonated carbon (δ C 108.6, C-2) in the NMR spectra of 1 were found to be absent in those of 2 (Tables 1 and 2). In contrast, compound 2 exhibited an additional resonance at δ H 5.61 (1H, s, H-2) in the 1 H-NMR spectrum that correlates with an olefinic methine carbon at δ C 87.9 (C-2), accounting for the shielding effect of the carbonyl group at C-1. The key HMBC correlations from H-2 to C-3 (δ C 170.7) and C-4 (δ C 106.0) as well as other key HMBC and COSY correlations ( Figure 2) established the structure of compound 2 as shown.
The relative configuration of the 2,5-dihydrofuran ring and the geometries of the double bonds at C-6, C-8, and C-10 were determined to be the same as those of compound 1 with an analysis of the NOESY correlations ( Figure 3). The absolute configuration of 2 was assigned as 12S and 15R, which was evidenced by a comparison of the overlaid experimental and calculated ECD curves at the BH&HLYP/TZVP level ( Figure 5). COSY correlations (Figure 2) established the structure of compound 2 as shown.
The relative configuration of the 2,5-dihydrofuran ring and the geometries of the double bonds at C-6, C-8, and C-10 were determined to be the same as those of compound 1 with an analysis of the NOESY correlations ( Figure 3). The absolute configuration of 2 was assigned as 12S and 15R, which was evidenced by a comparison of the overlaid experimental and calculated ECD curves at the BH&HLYP/TZVP level ( Figure 5).    13 C-NMR data revealed 3 to be structurally related to compounds 1 and 2, featuring a methylated α-pyrone moiety and a 2,5-dihydrofuran ring linked with a shorter polyene linker. This deduction was supported by the NMR evidence that resonances for one of the trisubstituted double bonds in the conjugated polyene linker, comprising a singlet methyl (δ H 1.95; δ C 17.9, CH 3 -19), a singlet olefinic methine (δ H 5.93; δ C 137.1, CH-9), and an olefinic nonprotonated carbon (δ C 131.3, C-8) in the NMR spectra of 1 disappeared in those of 3. Combined with the key HMBC correlations from H-7 (δ H 6.08) to C-5 (δ C 158.4) and C-9 (δ C 138.3) as well as from H 3 -18 (δ H 1.33) to C-9 (δ C 138.3), C-10 (δ C 87.3) and C-11 (δ C 127.1) (Figure 2), the planar structure of compound 3 was thus established.
The observed NOESY correlations between H-13 and H 3 -18 helped to establish the relative configuration of the 2,5-dihydrofuran moiety ( Figure 3). The E-geometry of the double bonds at C-6 and C-8 was determined by the observed NOEs from H-7 to H-9 and H 3 -15 ( Figure 3). The absolute configuration of 3 was assigned as 10S and 13R based on the results the of TDDFT-ECD calculations at the BH&HLYP/TZVP level as shown in Figure 6. The observed NOESY correlations between H-13 and H3-18 helped to establish the relative configuration of the 2,5-dihydrofuran moiety (Figure 3). The E-geometry of the double bonds at C-6 and C-8 was determined by the observed NOEs from H-7 to H-9 and H3-15 ( Figure 3). The absolute configuration of 3 was assigned as 10S and 13R based on the results the of TDDFT-ECD calculations at the BH&HLYP/TZVP level as shown in Figure 6. Poloncosidin D (4) was isolated as a colorless oil. The molecular formula was established as C21H30O4 from its HRESIMS data with seven indices of hydrogen deficiency. The 1 H-and 13 C-NMR spectra of 4 revealed a 2,5-dihydrofuran moiety with three sets of Poloncosidin D (4) was isolated as a colorless oil. The molecular formula was established as C 21 H 30 O 4 from its HRESIMS data with seven indices of hydrogen deficiency. The 1 H-and 13 C-NMR spectra of 4 revealed a 2,5-dihydrofuran moiety with three sets of conjugated double bonds, conforming to the characteristics of verrucosidin derivatives. However, additional resonances for two carbonyl carbons (δ C 198.5, C-3 and 171.4, C-1) were observed in the 13 C-NMR spectrum of 4 that are quite different from those of compounds 1-3. The COSY correlation from H-2 (δ H 4.52) to H 3 -14 (δ H 1.21) as well as HMBC correlations from H 3 -14 to C-1 and C-3, from H-5 (δ H 7.19) to C-3, C-7 (δ C 140.8), and C-16 (δ C 17.8), and from H 3 -21 (δ H 3.60) to C-1 established the structure of 4, which possesses the same polyketide skeleton as that of compound 3, but without the formation of the α-pyrone moiety.
The relative configuration of 4 was established by the observed NOESY correlations (Figure 3), except for the stereogenic center of C-2. To establish the stereochemistry, TDDFT-ECD calculations for four possible isomers of 4 [4a- (2R,10S,13R), 4b-(2S,10S,13R),  4c-(2R,10R,13S), and 4d-(2S,10R,13S)] (Figure 7) were performed, and then two isomers [4a-(2R,10S,13R) and 4b-(2S,10S,13R)] were recommended for 4 ( Figure 8), which were further clarified by the DP4+ probability analysis. Both the proton and carbon data of the two possible isomers were calculated based on the DP4+ protocol and the results were analyzed with the experimental values [16]. The statistical analysis indicated that the isomer [4a-(2R,10S,13R)] was the equivalent structure with a probability of 99.91% (Figures 8 and S42). The absolute configuration of 4 was thus assigned. further clarified by the DP4+ probability analysis. Both the proton and carbon data of the two possible isomers were calculated based on the DP4+ protocol and the results were analyzed with the experimental values [16]. The statistical analysis indicated that the isomer [4a-(2R,10S,13R)] was the equivalent structure with a probability of 99.91% (Figures 8 and S42). The absolute configuration of 4 was thus assigned. The molecular formula of poloncosidin E (5) was assigned as C24H32O4 based on the accurate mass data measured by HRESIMS, which was identical to that of 1. The 1 H-and 13 C-NMR data of 5 were virtually identical to those obtained for 1 (Tables 1 and 2), except for the shielded location of C-5 at δC 155.9 and de-shielded location of C-18 at δC 22.7 in the 13 C-NMR spectrum of 5, conforming to the γ-gauche effect [17]. Moreover, resonances for C-6 (δC 125.6) and C-7 (δC 133.0) also moved to the shielded field in the 13 C-NMR spectrum of 5. A detailed analysis of the COSY and HMBC correlations of 5 ( Figure 2) revealed that it possessed the same planar structure as 1, but with Z-geometry for the double bond at C-6. The structure of compound 5 was thus established with the above further clarified by the DP4+ probability analysis. Both the proton and carbon data of the two possible isomers were calculated based on the DP4+ protocol and the results were analyzed with the experimental values [16]. The statistical analysis indicated that the isomer [4a-(2R,10S,13R)] was the equivalent structure with a probability of 99.91% (Figures 8 and S42). The absolute configuration of 4 was thus assigned. The molecular formula of poloncosidin E (5) was assigned as C24H32O4 based on the accurate mass data measured by HRESIMS, which was identical to that of 1. The 1 H-and 13 C-NMR data of 5 were virtually identical to those obtained for 1 (Tables 1 and 2), except for the shielded location of C-5 at δC 155.9 and de-shielded location of C-18 at δC 22.7 in the 13 C-NMR spectrum of 5, conforming to the γ-gauche effect [17]. Moreover, resonances for C-6 (δC 125.6) and C-7 (δC 133.0) also moved to the shielded field in the 13 C-NMR spectrum of 5. A detailed analysis of the COSY and HMBC correlations of 5 ( Figure 2) revealed that it possessed the same planar structure as 1, but with Z-geometry for the double bond at C-6. The structure of compound 5 was thus established with the above NMR data. The molecular formula of poloncosidin E (5) was assigned as C 24 H 32 O 4 based on the accurate mass data measured by HRESIMS, which was identical to that of 1. The 1 H-and 13 C-NMR data of 5 were virtually identical to those obtained for 1 (Tables 1 and 2), except for the shielded location of C-5 at δ C 155.9 and de-shielded location of C-18 at δ C 22.7 in the 13 C-NMR spectrum of 5, conforming to the γ-gauche effect [17]. Moreover, resonances for C-6 (δ C 125.6) and C-7 (δ C 133.0) also moved to the shielded field in the 13 C-NMR spectrum of 5. A detailed analysis of the COSY and HMBC correlations of 5 ( Figure 2) revealed that it possessed the same planar structure as 1, but with Z-geometry for the double bond at C-6. The structure of compound 5 was thus established with the above NMR data.
Unfortunately, the NOESY spectrum of 5 was not obtained due to instability and low abundance, and compound 5 was actually isolated as a highly impure mixture. To further define the geometry of the double bond at C-6, a DP4+ probability analysis was performed. The steric configuration of the 2,5-dihydrofuran moiety was determined as the same as that of compounds 1-4 based on the NMR data comparisons and biosynthetic origins. Two possible isomers [5a- (6E,12S,15R) and 5b-(6Z,12S,15R)] of compound 5 were thus generated for DP4+ calculations (Figure 9). The 13 C-NMR data calculated for the isomer [5b-(6Z,12S,15R)] matched well with the experimental data of 5 (probability = 100%), which led to the assignment of the steric configuration for 5. Poloncosidin F (6) was obtained as a colorless oil with a molecular formula of C23H32O6 as determined by the HRESIMS data. Its 1 H-and 13 C-NMR data (Tables 1 and 2) revealed a close structural relationship with deoxyverrucosidin (7) [3]. One of the visible differences were that the signals for a methoxy group 3-OCH3 (δH 3.83; δC 60.2) in the NMR spectra of deoxyverrucosidin (7) were absent in those of 6. Moreover, a de-shielded resonance for C-3 (δC 177.0) and a shielded field shift for C-2 (δC 89.9) were clearly observed in the 13 C-NMR spectrum of 6 ( Table 2), compared to those of 7 [3]. These spectroscopic features as well as in considering the molecular formula revealed that the ester linkage of the α-pyrone ring in deoxyverrucosidin (7) was hydrolyzed in 6. Supported by the HMBC correlations from H3-16 (δH 1.57) to C-1 (δC 165.3), C-2, and C-3 as well as from H3-17 (δC 1.72) to C-3, C-4 (δC 111.3), and C-5 (δC 155.4) (Figure 2), the planar structure of 6 was established.
The E-geometry of five conjugated double bonds at C-2, C-4, C-6, C-8, and C-10 were assigned by the NOE correlations from H3-17 to H3-16 and H3-18, from H3-18 to H3-19, and from H-9 to H-7 and H-11. NOE correlations from H-13 to H3-22 and H3-23 indicated the cofacial orientation of these groups, while the cross-peak from H-15 to H3-21 assigned these groups at the opposite face. Similarly, the absolute configuration of 6 was determined by the TDDFT-ECD calculations in Gaussian 09. The ECD spectrum calculated at the BH&HLYP/TZVP level for isomer (12S,13S,14R,15R)-6 matched well with the experimental curve (Figure 10), leading to the establishment of the absolute configuration of 6. Figure 10. Experimental ECD spectrum of 6 (black line) compared with that calculated for Poloncosidin F (6) was obtained as a colorless oil with a molecular formula of C 23 H 32 O 6 as determined by the HRESIMS data. Its 1 H-and 13 C-NMR data (Tables 1 and 2) revealed a close structural relationship with deoxyverrucosidin (7) [3]. One of the visible differences were that the signals for a methoxy group 3-OCH 3 (δ H 3.83; δ C 60.2) in the NMR spectra of deoxyverrucosidin (7) were absent in those of 6. Moreover, a de-shielded resonance for C-3 (δ C 177.0) and a shielded field shift for C-2 (δ C 89.9) were clearly observed in the 13 C-NMR spectrum of 6 ( Table 2), compared to those of 7 [3]. These spectroscopic features as well as in considering the molecular formula revealed that the ester linkage of the α-pyrone ring in deoxyverrucosidin (7) was hydrolyzed in 6. Supported by the HMBC correlations from H 3 -16 (δ H 1.57) to C-1 (δ C 165.3), C-2, and C-3 as well as from H 3 -17 (δ C 1.72) to C-3, C-4 (δ C 111.3), and C-5 (δ C 155.4) (Figure 2), the planar structure of 6 was established.
The E-geometry of five conjugated double bonds at C-2, C-4, C-6, C-8, and C-10 were assigned by the NOE correlations from H 3 -17 to H 3 -16 and H 3 -18, from H 3 -18 to H 3 -19, and from H-9 to H-7 and H-11. NOE correlations from H-13 to H 3 -22 and H 3 -23 indicated the cofacial orientation of these groups, while the cross-peak from H-15 to H 3 -21 assigned these groups at the opposite face. Similarly, the absolute configuration of 6 was determined by the TDDFT-ECD calculations in Gaussian 09. The ECD spectrum calculated at the BH&HLYP/TZVP level for isomer (12S,13S,14R,15R)-6 matched well with the experimental curve ( Figure 10), leading to the establishment of the absolute configuration of 6.
In addition to the new compounds 1-6, the known analogue deoxyverrucosidin (7) was also isolated and identified by a detailed spectroscopic analysis and comparisons with the reported data [3].
The E-geometry of five conjugated double bonds at C-2, C-4, C-6, C-8, and C-10 were assigned by the NOE correlations from H3-17 to H3-16 and H3-18, from H3-18 to H3-19, and from H-9 to H-7 and H-11. NOE correlations from H-13 to H3-22 and H3-23 indicated the cofacial orientation of these groups, while the cross-peak from H-15 to H3-21 assigned these groups at the opposite face. Similarly, the absolute configuration of 6 was determined by the TDDFT-ECD calculations in Gaussian 09. The ECD spectrum calculated at the BH&HLYP/TZVP level for isomer (12S,13S,14R,15R)-6 matched well with the experimental curve (Figure 10), leading to the establishment of the absolute configuration of 6. In addition to the new compounds 1-6, the known analogue deoxyverrucosidin (7) was also isolated and identified by a detailed spectroscopic analysis and comparisons with the reported data [3].
The obtained compounds, except for compound 5, were assayed for antimicrobial activities against human-, aqua-, and plant-pathogenic microbes. These compounds

General Experimental Procedures
The instruments, organic solvents, and reagents used in this work were the same as those used in our previous reports [13][14][15].

Fungal Material
The fungus P. polonicum CS-252 was isolated from the cold seep sediments, which were collected from South China Sea (depth 1183 m) in September 2020. The fungal strain was identified as Penicillium polonicum based on the ITS region sequence, which is the same (100%) as that of P. polonicum FNG3 (with accession no. MZ301257.1). The sequence data were deposited in GenBank with the accession no. OK175489.1. The strain is being preserved at the Key Laboratory of Experimental Marine Biology, Institute of Oceanology, Chinese Academy of Sciences (IOCAS).

Fermentation, Extraction and Isolation
For chemical investigations, 1 L Erlenmeyer flasks with a solid medium containing rice were autoclaved at 120 • C for 20 min before inoculation. Each flask contained 70 g rice, 0.3 g peptone, 0.5 g yeast extract, 0.2 g corn steep liquor, 0.1 g monosodium glutamate, and 100 mL of naturally sourced and filtered seawater obtained from the Huiquan Gulf of the Yellow Sea near the campus of IOCAS. The fresh mycelia of P. polonicum CS-252 were grown on a PDA medium at 28 • C for seven days, inoculated into the autoclaved medium and then statically cultured for 30 days at room temperature. The fermented rice substrate was fragmented mechanically and extracted thoroughly with EtOAc. The above EtOAc solutions were combined and concentrated under reduced pressure to produce 70.0 g organic extract.
The extract was subjected to vacuum liquid chromatography (VLC) over Si gel and fractionated using solvent mixtures of increasing polarity consisting of petroleum ether (PE) and EtOAc (20:1 to 1:1) and finally with CH 2 Cl 2 -MeOH (50:1 to 1:1) to yield 10 fractions

Fermentation, Extraction and Isolation
Conformational searches were carried out using molecular mechanics with the MM+ method and HyperChem 8.0 software. Afterward, the geometries were optimized at the gas-phase B3LYP/6-31G(d) level in Gaussian09 software to identify the energy-minimized conformers. Then, the optimized conformers were subjected to the calculations of ECD spectra using the TD-DFT at the BH&HLYP/TZVP level. Simultaneously, the solvent effects of the MeOH solution were evaluated at the same DFT level using the SCRF/PCM method [18].

ECD Calculations
Conformational searches were carried out via molecular mechanics with the MM+ method in HyperChem 8.0 software. Afterward, the geometries were optimized at the gas-phase B3LYP/6-31G(d) level in Gaussian09 software to identify the energy-minimized conformers. Then, the optimized conformers were subjected to calculations of the ECD spectra using the TDDFT at the BH&HLYP/TZVP level. Simultaneously, solvent effects of the MeOH solution were evaluated at the same DFT level using the SCRF/PCM method. [18] 3.6. Computational NMR Chemical Shift Calculations and DP4+ Analyses All the theoretical calculations were performed using the Gaussian 09 program package. Conformational searches for possible isomers were conducted with molecular mechanics using the MMFF method with Macromodel software (Schrödinger, LLC) and the corresponding stable conformer, from which distributions higher than 2% were collected. Then, the B3LYP/6-31G(d) PCM level in DMSO was used to optimize the conformers. After that, the NMR shielding tensors of all optimized conformers were calculated using the DFT method at the mPW1PW91\6-31+G(d) PCM level on DMSO and then averaged based on Boltzmann's distribution theory [16]. GIAO (gauge-independent atomic orbital) NMR chemical calculations were performed using an equation described previously [19]. Finally, the NMR chemical shifts and shielding tensors ( 1 H and 13 C) were analyzed and compared with the experimental chemical shifts using DP4+ probability [16].

Conclusions
In summary, six new verrucosidin derivatives (1-6) were identified from the deep sea cold-seep-derived Penicillium polonicum CS-252. Compounds 1-5 represent the first examples of verrucosidins with a 2,5-dihydrofuran ring. The structures of these compounds were determined using a combination of methods (NMR, MS, ECD calculations and DP4+ probability analysis).