Next Article in Journal
The In Vitro Anti-Tumor Activity of Phycocyanin against Non-Small Cell Lung Cancer Cells
Next Article in Special Issue
New Cembranoids and a Biscembranoid Peroxide from the Soft Coral Sarcophyton cherbonnieri
Previous Article in Journal
Structural Elucidation of Irish Organic Farmed Salmon (Salmo salar) Polar Lipids with Antithrombotic Activities
Previous Article in Special Issue
Anticancer Activity of Euplotin C, Isolated from the Marine Ciliate Euplotes crassus, Against Human Melanoma Cells
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Communication

Asperversins A and B, Two Novel Meroterpenoids with an Unusual 5/6/6/6 Ring from the Marine-Derived Fungus Aspergillus versicolor

Hubei Key Laboratory of Natural Medicinal Chemistry and Resource Evaluation, Tongji Medical College, Huazhong University of Science and Technology, Wuhan 430030, China
*
Authors to whom correspondence should be addressed.
Mar. Drugs 2018, 16(6), 177; https://doi.org/10.3390/md16060177
Submission received: 17 April 2018 / Revised: 14 May 2018 / Accepted: 14 May 2018 / Published: 23 May 2018
(This article belongs to the Special Issue Terpenoids from Marine Organisms)

Abstract

:
Asperversins A (1) and B (2), two novel meroterpenoids featuring an uncommon 5/6/6/6 ring system, along with five new analogues (37) and a known compound asperdemin (8), were obtained from the marine-derived fungus Aspergillus versicolor. Their structures and absolute configurations were confirmed by extensive spectroscopic analyses, single-crystal X-ray diffraction studies, and electronic circular dichroism (ECD) calculation. All new compounds were tested for their acetylcholinesterase enzyme (AChE) inhibitory activities and cytotoxic activities, of which compound 7 displayed moderate inhibitory activity against the AChE with an IC50 value of 13.6 μM.

1. Introduction

Meroterpenoids are hybrid secondary metabolites that partially derive from the terpenoid pathways [1]. Naturally occurring meroterpenoids have been isolated from a variety of sources including animals, plants, bacteria, and fungi, [2] and are exemplified by ubiquinone-10 (coenzyme Q10) [3], α-tocopherol (vitamin E) [4], vinblastine [5], merochlorin A [6,7], and teleocidin B-4 [8].
Biosynthetically, the complex structures of fungal meroterpenoids are mostly derived from simple precursors like a linear isoprenoid or the C-2 carbon unit acetyl-CoA, via a series of chemical transformations catalyzed by two enzyme families, terpene cyclases and polyketide synthases (PKSs) [9,10]. In addition to the enormous structural diversity, fungal meroterpenoids have attracted wide interest from the scientific community due to their broad spectrum of pharmacological activities [11,12,13,14].
During our investigation of bioactive metabolites from fungi [15,16,17,18,19], many novel meroterpenoids were isolated from various fungi. Asperterpene A, obtained from Aspergillus terreus, represents a potential lead compound and a versatile scaffold for the development of drugs for the treatment of Alzheimer’s disease [15]. Apiroaspertrione A, bearing a unique spiro[bicyclo[3.2.2]nonane-2,1′-cyclohexane] carbocyclic skeleton, shows an effective activity against methicillin-resistant Staphylococcus aureus (MRSA) [16]. The metabolites of the fungus Aspergillus versicolor, that were obtained from mud from the South China Sea, are reported here. This endeavor resulted in the isolation of seven undescribed meroterpenoids, called as asperversins A–G (17), and a structurally related known compound, asperdemin (8) [20]. Notably, asperversins A (1) and B (2) have a unique 5/6/6/6 ring system that features a tetrahydrofuran ring. Herein, the isolation, structure elucidation, plausible biosynthetic pathway, and the biological evaluations of compounds 17 are presented.

2. Results and Discussion

Asperversin A (1), a white amorphous powder, was established to have a molecular formula of C24H32O8 by HRESIMS ion at m/z 449.2173 ([M + H]+, calcd. for C24H33O8, 449.2175), indicating 9 degrees of unsaturation. The strong IR absorptions at 1741 and 1703 cm−1 suggested the presence of ester and unsaturated ester groups. The 1H NMR data (in methanol-d4, Table 1) displayed an olefinic proton at δH 5.91 (d, J = 1.0 Hz, H-15), a methoxy group at δH 3.70 (s, OMe-24) and six methyl groups at δH 2.20 (d, J = 1.0 Hz, H3-17), 2.07 (s, H3-23), 1.43 (H3-18), 1.31 (H3-19), 1.26 (H3-20), and 1.18 (H3-21). The 13C and DEPT spectra revealed the presence of three carbonyl carbons (δC 173.6, 171.6, and 167.1) and two enol systems (δC 165.1, 162.2, 102.2, and 98.5), three quaternary sp3 carbons (two oxygenated), four methines (two oxygenated), three methylene carbons, and seven methyl groups.
Comparison of the 1H and 13C NMR data (Table 1 and Table 2) of 1 with the known compound asperdemin (8) [20] revealed their close similarities, except for the presence of an additional methoxy group (δH 3.70, δC 52.3) and an acetoxyl unit. Further analyses of the 1H–1H COSY and HMBC spectra of 1 indicated that 1 shared the same B/C/D ring system as that of 8 (Figure 1). The key HMBC correlation from H-6 to C-22 (δC 171.6) indicated the acetoxyl group at C-6. Furthermore, the 1H–1H COSY cross-peaks of H-1/H2-2 and the HMBC connections from H-1 to C-3, from H2-2 to C-3, and from OMe to C-3 confirmed the existence of a methyl ester group in 1, rather than a seven-membered lactone. Although there was no obvious HMBC correlation to connect C-1 and C-4, the absence of hydroxyl at C-1 confirmed by the 1D and 2D NMR data in DMSO-d6 and the chemical shifts of C-1 (δC 81.9 for 1; 68.0 for 8) and C-4 (δC 77.7 for 1; 84.8 for 8) indicated that C-1 and C-4 were connected with an ether bond to establish a tetrahydrofuran ring, which was consistent with the unsaturation degrees deduced by HRESIMS. Thus, the unusual tetrahydrofuran A ring replaced the seven-membered lactone ring, which distinguished compound 1 from other metabolites reported in this compound class.
The relative configuration of 1 was partially elucidated by the examination of its Nuclear Overhauser Effect SpectrocopY (NOESY) spectrum (Figure 2). The configuration of Me-19 was arbitrarily assigned as same as that of 8 in a β-orientation. Furthermore, the NOESY correlations of H3-18/H-11b and H-11b/H3-19 suggested their cofacial and β-orientation. In contrast, H-9 should be α-oriented. The NOESY correlations of H-1/H-9, H-9/H-5, and H-1/H-5 indicated their α-orientation. The spin-coupling constant (J = 2.3 Hz) of H-5 and H-6 suggested H-6 was equatorial. To confirmed the absolute configuration, the ECD calculation of 1 [1S,5S,6R,8R,9R,10R-1 (1a) and 1R,5R,6S,8S,9S,10S-1 (1b)] was conducted in MeOH using the time-dependent density functional theory (TD-DFT) method at the B3LYP/6-311+g (d, p) level. The experimental ECD of 1 matched well with the calculated ECD curve of 1a, indicating the absolute configuration of 1 as 1S,5S,6R,8R,9R,10R (Figure 3).
The HRESIMS spectrum suggested a molecular formula of C23H30O8 for compound 2, indicating a molecular weight of 14 mass units lower than 1. The 1H and 13C NMR data of 2 (Table 1 and Table 2) closely resembled those of compound 1, except for the difference of chemical shift of carbonyl (C-3: δC 176.7 for 2; 173.6 for 1), and the absence of the methoxy group. Compound 2 was speculated to be the carboxylic acid analogue of 1. The structure of 2 was confirmed by extensive analyses of its 1H–1H COSY and HMBC spectra. Compound 2 was confirmed to have the same relative and absolute configurations as 1 by the NOESY spectra and their similar ECD curves (Figure 1, Figure 3 and Figure 4).
Asperversin C (3), a white amorphous powder, has a molecular formula of C26H34O9, which was determined on the basis of the HRESIMS ion peak at m/z 491.2269 ([M + H]+, calcd. for C26H35O9, 491.2281). The similar IR, UV, 1D, and 2D NMR spectra of 3 and 1 suggested both compounds shared the same B/C/D rings. A comparison of their 1D and 2D NMR data indicated the presence of a double bond (δH 5.05 and 5.08, δC 117.7, CH2-20; δC 144.4, C-4) and the absence of a methyl group in compound 3. The location of the double bond between C-4 and C-20 was confirmed by the HMBC correlations from H2-20 to C-4, C-5, and C-21, and from H3-21 to C-4 and C-20. Moreover, the HMBC correlations from H-1 to C-24 and H3-25 to C-24 defined the location of an additional acetoxyl unit at C-1. Finally, the planar structure of 3 with the A ring opened was confirmed as shown (Figure 1). The NOE correlations from H-1/H3-19, H3-18/H3-19 suggested H-1, Me-18, and Me-19 were β-oriented. In contrary, H-5, H-6 and H-9 were α-oriented based on the correlations of H-5/H-6 and H-5/H-9. In consideration of the similar specific rotation and ECD spectra (Figure 4), the absolute configuration of 3 was assigned as 1R,5S,6R,8R,9R,10R.
Asperversin D (4) possesses a molecular formula of C21H28O6 with 16 units less than that of compound 8 on the basis of HRESIMS data. Comparison of the 1H and 13C NMR data of 4 with those of 8 revealed the presence of an additional methylene group and the absence of an oxygenated methine group in 4. Inspection of the 2D NMR of 4 indicated this methylene group (δC 23.3; δH 1.67 and 1.92) should be assigned to C-6 by the 1H–1H COSY correlations of H-5/H2-6/H-7. The gross structure of 4 and its relative configuration were then confirmed by comprehensive analyses of the 1H–1H COSY, HMBC, and NOESY spectra. Finally, the absolute configuration of 4 was ascertained as 1R,5S,8R,9R,10R by a similar ECD spectrum to that of 1 (Figure 4).
Asperversin E (5), colorless crystals, exhibits a molecular formula of C23H30O8 with nine degrees of unsaturation from a HRESIMS ion peak at m/z 457.1830 ([M + Na]+, calcd. for 457.1838). Close similarities were observed in the 13C NMR data (Table 2) between compounds 5 and 4. Inspection of the 1D NMR spectra, HSQC, and HMBC spectra revealed that 5 contained an additional acetoxyl group (δH 2.06, δC 21.4, and 170.2). Detailed HMBC correlations from H3-23 (δH 2.06) to C-22 (δC 170.2) and from H-6 (δH 5.52) to C-22 confirmed the location. The X-ray crystallographic diffraction (Figure 5) revealed the configuration of 5 was 1R,5S,6R,8R,9R,10R.
A protonated molecular ion peak at m/z 477.2098 ([M + H]+, calcd. 477.2125) was obtained from the HRESIMS experiment and indicated asperversin F (6) has the molecular formula C25H32O9. The 1H and 13C NMR spectra of 6 showed there were two more acetoxyl groups (δH 2.14, δC 21.3, and 171.2; δH 1.99, δC 20.6, and 171.1) than 4. The locations of two acetoxyls were assigned to C-1 and C-6 by the HMBC correlations from H-1 (δH 4.96) to C-24 and from H-6 (δH 5.64) to C-22, respectively. The relative configuration and absolute configuration of 6 were determined by analyzing the NOESY spectrum and by comparing its ECD spectrum with that of 5 (Figure 4).
Asperversin G (7), colorless cubic crystals, has a molecular formula of C23H28O6 based on the HRESIMS spectrum. The 1H and 13C NMR data (Table 1 and Table 2) of 7 resembled those of 5, with the main difference being the presence of an α,β-unsaturated ketone group (δH 7.08, d, J = 10.2 Hz and 5.87, d, J = 10.2 Hz; δC 203.2, 155.1, and 125.1). The α,β-unsaturated ketone further constructed an 2-cyclohexen-1-one (ring A), as ascertained by the HMBC correlations from H3-19 to C-1, C-5, C-9, and C-10, from H-1 to C-3, from H-2 to C-3 and C-4, and from H3-20 and H3-21 to C-3, C-4, and C-5. The remainder of 7 was identical to 5 as elucidated by 2D NMR spectra. The absolute configuration of 7 was established unambiguously as 5S,6R,8R,9R,10R by a single-crystal X-ray diffraction (Figure 5).
Asperversins A (1) and B (2) were found to be the first examples of pyrone meroterpenoids featuring an exclusive 5/6/6/6 ring system compared to all known meroterpenoids [2]. Therefore, a plausible biosynthetic pathway is proposed to illustrate the generation of 1 and 2 (Scheme 1). Asperversins A (1) and B (2) are probably biosynthesized via a polyketide and mevalonate hybrid biogenetic pathway [2,21,22]. After a series of oxidations, the A ring of I was opened, and then, followed by the key electrophilic addition of OH-4 with the α,β-unsaturated carboxylic acid bond, the intermediate III, featuring a tetrahydrofuran ring, was derived from II. Compounds 2 and 1 were produced subsequently by an additional acetylation and methylation.
All compounds (18) were inactive for cytotoxicity against A549, MCF-7, HepG2, and HL-60 cells, and did not have any antimicrobial effects against methicillin-resistant Staphylococcus aureus (MRSA), methicillin-sensitive Staphylococcus aureus (MSSA), Escherichia coli, or Pseudomonas aeruginosa. The acetylcholinesterase (AChE) inhibition effects of these compounds (18) were assessed by Ellman’s spectrophotometric method [23] using human recombinant AChE with galanthamine as control compound (Table 3). Compound 7 exhibited the strongest inhibition to AChE with an IC50 value of 13.6 μM, and compounds 16 and 8 showed no activities up to a concentration of 40 μM.
The molecular docking studies were conducted in order to get an insight into the binding pattern and extent of binding of compound 7 with the target enzyme (Figure 6). This was mainly attributed to the basic skeleton of compound 7 which provided better binding prospects with the formation of hydrogen bonding interactions to the amino acid residues Glu 291 and Tyr 341 of the protein. So we speculated the α,β-unsaturated ketone group was the key pharmacophore for the acetylcholinesterase (AChE) inhibition effect.

3. Experimental Section

3.1. General Experimental Procedures

The NMR spectra were recorded on Bruker AM-400 and 600 spectrometers (Bruker, Karlsruhe, Germany). The 1H and 13C NMR chemical shifts were referenced to the solvent or solvent impurity peaks for methanol-d4 (δH 3.31 and δC 49.0), and acetone-d6 (δH 2.05 and δC 206.3 and 29.8). The FT-IR spectra were determined with a Bruker Vertex 70 instrument (Bruker, Karlsruhe, Germany). The UV and ECD spectra were measured using a Perkin Elmer Lambda 35 UV spectrophotometer (PerkinElmer, Inc., Fremont, CA, USA) and a JASCO-810 ECD spectrometer (JASCO Co., Ltd., Tokyo, Japan), respectively. Optical rotations were measured by an AUTOPOL IV-T Automatic polarimeter (Rudolph Research Analytical, Hackettstown, NJ, USA). HRESIMS data were obtained in the positive ion mode on a Thermo Fisher LTQ XL spectrometer (Thermo Fisher, Palo Alto, CA, USA). Crystal X-ray diffraction data were measured on a Rigaku XtaLAB PRO MM007HF (Rigaku, Tokyo, Japan). Semipreparative HPLC was carried out using an Agilent 1200 quaternary system with a UV detector using a reversed-phase C18 column (5 μm, 10 × 250 mm, Welch Ultimate XB-C18). Column chromatography (CC) was performed using silica gel (100–200 and 200–300 mesh; Qingdao Marine Chemical Inc., Qingdao, China), ODS (50 μm, YMC, Kyoto, Japan), and Sephadex LH-20 (GE Healthcare Bio-Sciences AB, Uppsala, Sweden). Thin-layer chromatography (TLC) was performed on silica gel 60 F254 (Yantai Chemical Industry Research Institute, Yantai, China).

3.2. Fungi Material

The fungus Aspergillus versicolor was isolated from mud of the South China Sea. The sequence data for this strain have been submitted to DDBJ/EMBL/GenBank under accession no. 2081031. A voucher sample has been deposited in the culture collection of Tongji Medical College, Huazhong University of Science and Technology, P. R. China.

3.3. Fermentation and Extraction

The Aspergillus versicolor strain was cultured on potato dextrose agar (PDA) at 28 °C for 7 days to prepare the seed culture and inoculated into 100 sterilized Erlenmeyer flasks (1 L), which contained 200.0 g of rice and 160.0 mL of distilled water, then incubated at 28 °C for 28 days. Following incubation, the rice was soaked with 95% EtOH many times until the solvent was near-colorless at room temperature. Finally, solvent removal afforded a brown extract that was suspended in water (2 L) and extracted with EtOAc (1:1, v/v) three times. After removal of EtOAc, a dark brown gum (170 g) was obtained.

3.4. Isolation

The extracts (170.0 g) were subjected to a silica gel chromatography column (CC) eluting with PE/ EtOAc (20:1–0:1) progressively to obtain six fractions (Fr. 1–Fr. 6). Fr.3 was further separated with silica gel CC to yield four subfractions (Fr. 3.1–Fr. 3.4). The subfraction Fr. 3.1 was subjected to a Sephadex LH-20 (MeOH) to afford three parts (Fr. 3.1a–Fr. 3.1c). After crude purification of Fr. 3.1a by ODS column (MeOH-H2O, 20:80–70:30, v/v), the second part of Fr. 3.1a (Fr. 3.1a-2, MeOH–H2O, 60:40, v/v) was purified by repeated semi-preparative HPLC (MeCN–H2O, 55:45, 2 mL/min) to yield 1 (9.3 mg, tR 25.2 min), 2 (6.2 mg, tR 28.5 min), 7 (23.6 mg, tR 30.1 min).
Three subfractions (Frs. 3.2a–3.2c) were obtained from Fr. 3.2 by chromatography on an ODS column (MeOH–H2O, 20:80–70:30, v/v). Fr.3.2b was chromatographed on ODS (MeOH–H2O, 40:60–70:30) to afford a mixture (4, 5, 6, 8) in 40% yield. The mixture was purified by semi-preparative HPLC (MeCN–H2O, 40:60–60:40, 2 mL/min) to yield 8 (16.3 mg, tR 25.5 min, 45:55), 5 (200.0 mg, tR 27.2 min, 45:55), 4 (3.7 mg, tR 21.0 min, 40:60), and 6 (10.2 mg, tR 21.0 min, 40:60). Fr. 3.2c was chromatographed on ODS (MeOH–H2O, 30:70–80:20, v/v) to yield four parts. Fr. 3.2c-2 (the second part of Fr. 3.2c, MeOH–H2O, 50:50–60:40, v/v) was purified by repeated semi-preparative HPLC (MeCN–H2O, 46:54, 2 mL/min) to yield 3 (2.78 mg, tR 15.1 min).

3.4.1. Asperversin A (1)

C24H32O8; a white amorphous powder; [α]25D +46.7 (c 0.34, MeOH); UV (MeOH) λmax (log ε) 205 (4.51) and 286 (3.84) nm; IR (KBr) νmax 3445, 2958, 1741, 1710, 1654, 1588, 1444, 1403, 1384, 1237, 1194, 996 cm−1; ECD (MeOH) λmaxε) = 205 (+25.02) and 286 (+1.86) nm; For 1H NMR (400 MHz) and 13C NMR (100 MHz) data, see Table 1 and Table 2; HRESIMS 449.2173 ([M + H]+, calcd. for C24H33O8, 449.2175).

3.4.2. Asperversin B (2)

C23H30O8; a white amorphous powder; [α]25D +100.3 (c 0.29, MeOH); UV (MeOH) λmax (log ε) 202 (4.78) and 285 (4.14) nm; IR (KBr) νmax 3433, 2963, 2929, 1740, 1710, 1652, 1584, 1446, 1407, 1235, 1194, 1028 cm−1; ECD (MeOH) λmaxε) = 207 (+21.37) and 286 (+0.11) nm; For 1H NMR (400 MHz) and 13C NMR (100 MHz) data, see Table 1 and Table 2; HRESIMS m/z 457.1846 ([M + Na]+, calcd. for C23H30NaO8, 457.1838).

3.4.3. Asperversin C (3)

C25H32O9; a white amorphous powder; [α]25D +100.9 (c 0.23, MeOH); UV (MeOH) λmax (log ε) 205 (4.51) and 284 (3.87) nm; IR (KBr) νmax 3425, 2922, 1743, 1711, 1655, 1591, 1233, 1192, 1032 cm−1; ECD (MeOH) λmaxε) = 203 (+19.88) and 286 (+2.34) nm; For 1H NMR (400 MHz) and 13C NMR (100 MHz) data, see Table 1 and Table 2; HRESIMS m/z 491.2269 ([M + H]+, calcd. for C26H35O9, 491.2281).

3.4.4. Asperversin D (4)

C21H28O6; a white amorphous powder; [α]25D +98.4 (c 0.25, MeOH); UV (MeOH) λmax (log ε) 205 (4.77) and 286 (3.93) nm; IR (KBr) νmax 3429, 2928, 1702, 1653, 1586, 1449, 1408, 1386, 1313, 1274, 1154, 1113 cm−1; ECD (MeOH) λmaxε) = 205 (+23.13) and 285 (+2.24) nm; For 1H NMR (600 MHz) and 13C NMR (150 MHz) data, see Table 1 and Table 2; HRESIMS 377.1919 ([M + H]+, calcd. for C22H29O8, 377.1964).

3.4.5. Asperversin E (5)

C23H30O8; colorless crystals; [α]25D +75.5 (c 0.20, MeOH); UV (MeOH) λmax (log ε) 205 (4.57) and 286 (3.89) nm; IR (KBr) νmax 3431, 2988, 2931, 1714, 1653, 1588, 1447, 1402, 1233, 1128, 1086 cm−1; ECD (MeOH) λmaxε) = 208 (+28.96) and 281 (+1.86) nm; For 1H NMR (600 MHz) and 13C NMR (150 MHz) data, see Table 1 and Table 2; HRESIMS 457.1830 ([M + Na]+, calcd. for C23H30O8Na, 457.1838).

3.4.6. Asperversin F (6)

C25H32O9; a white amorphous powder; [α]25D +100.6 (c 0.48, MeOH); UV (MeOH) λmax (log ε) 205 (4.61) and 286 (3.77) nm; IR (KBr) νmax 3462, 2992, 2934, 1736, 1655, 1590, 1448, 1409, 1374, 1230, 1202, 1131, 1024 cm−1; ECD (MeOH) λmaxε) = 215 (+23.68) and 282 (+0.34) nm; For 1H NMR (400 MHz) and 13C NMR (100 MHz) data, see Table 1 and Table 2; HRESIMS 477.2143 ([M + H]+, calcd. for C25H33O9, 477.2098).

3.4.7. Asperversin G (7)

C23H28O6; colorless cubic crystals; [α]25D +110.2 (c 0.21, MeOH); UV (MeOH) λmax (log ε) 205 (4.60), 228 (4.18) and 285 (3.88) nm; IR (KBr) νmax 3434, 2975, 2927, 2862, 1793, 1706, 1674, 1587, 1448, 1408, 1257, 1232, 1141 cm−1; ECD (MeOH) λmaxε) = 215 (+16.60), 291 (+1.11) and 342 (–1.01) nm; For 1H NMR (400 MHz) and 13C NMR (100 MHz) data, see Table 1 and Table 2; HRESIMS [M + H]+ m/z 401.1977 (calcd. for C23H29O6, 401.1964).

3.5. X-ray Crystal Structure Analysis

The crystallographic data for 5 (deposition no. CCDC 1816743), and 7 (deposition no. CCDC 1813401) have been deposited in the Cambridge Crystallographic Data Centre. Copies of the data can be obtained free of charge from the Cambridge Crystallographic Data Centre.

3.5.1. Crystal Data for Asperversin E (5)

C23H30O8∙H2O (M = 452.48 g/mol), orthorhombic, space group P212121 (no. 19), a = 9.95215(4) Å, b = 12.89469(4) Å, c = 17.07108(6) Å, V = 2190.729(13) Å3, Z = 4, T = 100.01(10) K, μ (CuKα) = 0.881 mm−1, Dcalc = 1.372 g/cm3, 23,143 reflections measured (8.594 ≤ 2Θ ≤ 147.672), 4396 unique (Rint = 0.0181, Rsigma = 0.0097) which were used in all calculations. The final R1 was 0.0274 (I > 2σ(I)) and wR2 was 0.0772. The Flack parameter was 0.039(16).

3.5.2. Crystal Data for Asperversin G (7)

C23H28O6 (M = 400.45 g/mol), monoclinic, space group P21, a = 8.77354(16) Å, b = 13.09891(17) Å, c = 9.64200(17) Å, α = 90.0°, β = 116.686(2)°, γ = 90.0°, V = 990.06(3) Å3, Z = 2, T = 100.01(10) K, μ (CuKα) = 0.790 mm−1, Dcalc = 1.541 g/cm3, 23,143 reflections measured (10.268 ≤ 2Θ ≤ 147.12), 3462 unique (Rint = 0.0195, Rsigma = 0.0162) which were used in all calculations. The final R1 was 0.0276 (I > 2σ(I)) and wR2 was 0.0738. The Flack parameter was 0.11(7).

3.6. Computational Methods

Monte Carlo conformational searches were carried out by the Spartan’s 10 software using the Merck Molecular Force Field (MMFF). The conformers with a Boltzmann-population of over 5% were chosen for ECD calculations, and then the conformers were initially optimized at B3LYP/6-31+g (d, p) level in MeOH using the CPCM polarizable conductor calculation model. The theoretical calculation of ECD was conducted in MeOH using Time-dependent Density functional theory (TD-DFT) at the B3LYP/6-311+g (d, p) level for all conformers of compounds 1a and 1b. Rotatory strengths for a total of 50 excited states were calculated. ECD spectra were generated using the program SpecDis 1.6 (University of Würzburg, Würzburg, Germany) and GraphPad Prism 5 (University of California San Diego, USA) from dipole-length rotational strengths by applying Gaussian band shapes with sigma = 0.3 eV [24].

3.7. Acetylcholinesterase Inhibitory Activities Assay

The inhibitory activities of compounds 18 against AChE were evaluated by Ellman’s method with slight modification [25]. Briefly, 20 µL AChE (0.05 U/mL) and 20 µL different concentrations of compounds (0.1 M phosphate buffer, pH 8.0) containing 0.1% Triton X-100) were pre-incubated at 37 °C for 30 min, then 160 µL of the assay solution which consisted of a 0.1 M phosphate buffer pH 8.0, with the addition of 600 µM 5,5’-dithio-bis (2-nitrobenzoic acid) and 550 µM acetylthiocholine iodide (ATCh) were added in the 96-well plates. The absorbance value at 425 nm was recorded for 15 min and enzyme activity was calculated from the slope of the obtained linear trend. The reaction rates were compared and the % inhibition due to the presence of tested inhibitors was calculated. Each concentration was analyzed in triplicate, and IC50 values were determined graphically from log concentration inhibition curves (Graphpad Prism 5, GraphPad Software Inc., San Diego, CA, USA). Galanthamine was used as a standard inhibitor.

3.8. Docking Studies

Cocrystal structure of human AChE with compound 1 (PDB code: 4EY7) was obtained from the Protein Data Bank for docking tests (http://www.rcsb.org) [26]. The virtual screening was implemented in the Surflex-Dock module of the Sybyl software. Molecules were built with Chemdraw and optimized at molecular mechanical and semiempirical level by using Open Babel GUI. The crystallographic ligand was extracted from the active site and the designed ligands were modelled. All the hydrogen atoms were added to define the correct ionization and tautomeric states, and the carboxylate, phosphonate and sulfonate groups were considered in their charged form. In the docking calculation, the default FlexX scoring function was used for exhaustive searching, solid body optimization, and interaction scoring. Finally, the ligands with the lowest-energy and the most favorable orientation were selected [27].

4. Conclusions

In conclusion, eight meroterpenoids were isolated from the marine derived fungus Aspergillus versicolor, including seven new ones, asperversins A–E (17), and a known one, asperdemin (8). Notably, asperversins A (1) and B (2) possess a unique 5/6/6/6 ring system with a constructed tetrahydrofuran ring. The absolute structures of all new compounds were confirmed by comprehensive spectroscopic analyses: compound 1 was elucidated by experimental and calculated ECD, and compounds 5 and 7 were ascertained by the X-ray crystallographic diffractions. All isolated meroterpenoids were tested for their acetylcholinesterase enzyme (AChE) inhibitory activities, and compound 7 displayed moderate inhibitory activity against the AChE with an IC50 value of 13.6 μM.

Supplementary Materials

The following are available online at https://www.mdpi.com/1660-3397/16/6/177/s1, HRESIMS, 1D and 2D NMR, IR, and UV spectra of 17.

Author Contributions

H.L. conducted the main experiments, data analyzes, and wrote the manuscript; M.D. assisted the isolated experiments; W.S. performed the bioactivity analyzes; C.C. and H.Z. were responsible for the analysis of single-crystal X-ray crystallography. J.W. conducted the fermentation of fungi; C.Q. and Z.L. analyzed the spectroscopic data; J.W., Y.X., and Y.Z. designed the experiments and commented the manuscript. All authors reviewed the manuscript.

Funding

This work was financially supported by the Program for the Changjiang Scholars of Ministry of Education of the People’s Republic of China (No. T2016088); the National Science Fund for Distinguished Young Scholars (No. 8172500151); the National Natural Science Foundation of China (Nos. 81573316, 31770379, and 31670354); Innovative Research Groups of the National Natural Science Foundation of China (No. 81721005).

Acknowledgments

We thank the Analytical and Testing Center at Huazhong University of Science and Technology for assistance in testing of ECD, UV, and IR spectra.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Geris, R.; Simpson, T.J. Meroterpenoids produced by fungi. Nat. Prod. Rep. 2009, 26, 1063–1094. [Google Scholar] [CrossRef] [PubMed]
  2. Matsuda, Y.; Abe, I. Biosynthesis of fungal meroterpenoids. Nat. Prod. Rep. 2016, 33, 26–53. [Google Scholar] [CrossRef] [PubMed]
  3. Morton, R.A. Ubiquinone. Nature 1958, 182, 1764. [Google Scholar] [CrossRef] [PubMed]
  4. Fernholz, E. On the constitution of α-tocopherol. J. Am. Chem. Soc. 1938, 60, 700–705. [Google Scholar] [CrossRef]
  5. Moncrief, J.W.; Lipscomb, W.N. Structures of leurocristine (vincristine) and vincaleukoblastine. X-ray analysis of leurocristine methiodide. J. Am. Chem. Soc. 1965, 87, 4963–4964. [Google Scholar] [CrossRef] [PubMed]
  6. Kaysser, L.; Bernhardt, P.; Nam, S.J.; Loesgen, S.; Ruby, J.G.; Skewes-Cox, P.; Jensen, P.R.; Fenical, W.; Moore, B.S. Merochlorins A–D, Cyclic Meroterpenoid Antibiotics Biosynthesized in Divergent Pathways with Vanadium-Dependent Chloroperoxidases. J. Am. Chem. Soc. 2012, 134, 11988–11991. [Google Scholar] [CrossRef] [PubMed]
  7. Kaysser, L.; Bernhardt, P.; Nam, S.J.; Loesgen, S.; Ruby, J.G.; Skewes-Cox, P.; Jensen, P.R.; Fenical, W.; Moore, B.S. Correction to “Merochlorins A–D, Cyclic Meroterpenoid Antibiotics Biosynthesized in Divergent Pathways with Vanadium-Dependent Chloroperoxidases”. J. Am. Chem. Soc. 2014, 136, 14626. [Google Scholar] [CrossRef]
  8. Awakawa, T.; Zhang, L.H.; Wakimoto, T.; Hoshino, S.; Mori, T.; Ito, T.; Ishikawa, J.; Tanner, M.E.; Abe, I. A Methyltransferase Initiates Terpene Cyclization in Teleocidin B Biosynthesis. J. Am. Chem. Soc. 2014, 136, 9910–9913. [Google Scholar] [CrossRef] [PubMed]
  9. Fischbach, M.A.; Walsh, C.T. Assembly-line enzymology for polyketide and nonribosomal peptide antibiotics: Logic, machinery, and mechanisms. Chem. Rev. 2006, 106, 3468–3496. [Google Scholar] [CrossRef] [PubMed]
  10. Itoh, T.; Tokunaga, K.; Matsuda, Y.; Fujii, I.; Abe, I.; Ebizuka, Y.; Kushiro, T. Reconstitution of a fungal meroterpenoid biosynthesis reveals the involvement of a novel family of terpene cyclases. Nat. Chem. 2010, 2, 858–864. [Google Scholar] [CrossRef] [PubMed]
  11. Macias, F.A.; Varela, R.M.; Simonet, A.M.; Cutler, H.G.; Cutler, S.J.; Dugan, F.M.; Hill, R.A. Novel bioactive breviane spiroditerpenoids from Penicillium brevicompactum Dierckx. J. Org. Chem. 2000, 65, 9039–9046. [Google Scholar] [CrossRef] [PubMed]
  12. Zhang, Y.C.; Li, C.; Swenson, D.C.; Gloer, J.B.; Wicklow, D.T.; Dowd, P.F. Novel antiinsectan oxalicine alkaloids from two undescribed fungicolous Penicillium spp. Org. Lett. 2003, 5, 773–776. [Google Scholar] [CrossRef] [PubMed]
  13. Guo, C.J.; Knox, B.P.; Chiang, Y.M.; Lo, H.C.; Sanchez, J.F.; Lee, K.H.; Oakley, B.R.; Bruno, K.S.; Wang, C.C.C. Molecular Genetic Characterization of a Cluster in A. terreus for Biosynthesis of the Meroterpenoid Terretonin. Org. Lett. 2012, 14, 5684–5687. [Google Scholar] [CrossRef] [PubMed]
  14. Nakamura, H.; Matsuda, Y.; Abe, I. Unique chemistry of non-heme iron enzymes in fungal biosynthetic pathways. Nat. Prod. Rep. 2018. [Google Scholar] [CrossRef] [PubMed]
  15. Qi, C.; Bao, J.; Wang, J.; Zhu, H.; Xue, Y.; Wang, X.; Li, H.; Sun, W.; Gao, W.; Lai, Y.; et al. Asperterpenes A and B, two unprecedented meroterpenoids from Aspergillus terreus with BACE1 inhibitory activities. Chem. Sci. 2016, 7, 6563–6572. [Google Scholar] [CrossRef] [PubMed]
  16. He, Y.; Hu, Z.; Sun, W.; Li, Q.; Li, X.N.; Zhu, H.; Huang, J.; Liu, J.; Wang, J.; Xue, Y.; et al. Spiroaspertrione A, a Bridged Spirocyclic Meroterpenoid, as a Potent Potentiator of Oxacillin against Methicillin-Resistant Staphylococcus aureus from Aspergillus sp. TJ23. J. Org. Chem. 2017, 82, 3125–3131. [Google Scholar] [CrossRef] [PubMed]
  17. He, Y.; Hu, Z.X.; Li, Q.; Huang, J.F.; Li, X.N.; Zhu, H.C.; Liu, J.J.; Wang, J.P.; Xue, Y.B.; Zhang, Y.H. Bioassay-Guided Isolation of Antibacterial Metabolites from Emericella sp. TJ29. J. Nat. Prod. 2017, 80, 2399–2405. [Google Scholar] [CrossRef] [PubMed]
  18. Qiao, Y.B.; Zhang, X.T.; He, Y.; Sun, W.G.; Feng, W.Y.; Liu, J.J.; Hu, Z.X.; Xu, Q.Q.; Zhu, H.C.; Zhang, J.W.; et al. Aspermerodione, a novel fungal metabolite with an unusual 2,6-dioxabicyclo[2.2.1] heptane skeleton, as an inhibitor of penicillin-binding protein 2a. Sci. Rep. 2018, 8, 5454. [Google Scholar] [CrossRef] [PubMed]
  19. Zhou, P.; Wu, Z.; Tan, D.; Yang, J.; Zhou, Q.; Zeng, F.; Zhang, M.; Bie, Q.; Chen, C.; Xue, Y.; et al. Atrichodermones A–C, three new secondary metabolites from the solid culture of an endophytic fungal strain, Trichoderma atroviride. Fitoterapia 2017, 123, 18–22. [Google Scholar] [CrossRef] [PubMed]
  20. Yurchenko, A.N.; Smetanina, O.F.; Kalinovsky, A.I.; Pivkin, M.V.; Dmitrenok, P.S.; Kuznetsova, T.A. A new meroterpenoid from the marine fungus Aspergillus versicolor (Vuill.) Tirab. Russ. Chem. Bull. 2010, 59, 852–856. [Google Scholar] [CrossRef]
  21. Wang, J.F.; Wei, X.Y.; Qin, X.C.; Tian, X.P.; Liao, L.; Li, K.M.; Zhou, X.F.; Yang, X.W.; Wang, F.Z.; Zhang, T.Y.; et al. Antiviral Merosesquiterpenoids Produced by the Antarctic Fungus Aspergillus ochraceopetaliformis SCSIO 05702. J. Nat. Prod. 2016, 79, 59–65. [Google Scholar] [CrossRef] [PubMed]
  22. Kaur, A.; Raja, H.A.; Swenson, D.C.; Agarwal, R.; Deep, G.; Falkinham, J.O.; Oberlies, N.H. Talarolutins A–D: Meroterpenoids from an endophytic fungal isolate of Talaromyces minioluteus. Phytochemistry 2016, 126, 4–10. [Google Scholar] [CrossRef] [PubMed]
  23. Zhan, G.Q.; Liu, J.J.; Zhou, J.F.; Sun, B.; Aisa, H.A.; Yao, G.M. Amaryllidaceae alkaloids with new framework types from Zephyranthes candida as potent acetylcholinesterase inhibitors. Eur. J. Med. Chem. 2017, 127, 771–780. [Google Scholar] [CrossRef] [PubMed]
  24. Chen, S.H.; Chen, D.N.; Cai, R.L.; Cui, H.; Long, Y.H.; Lu, Y.J.; Li, C.Y.; She, Z.G. Cytotoxic and Antibacterial Preussomerins from the Mangrove Endophytic Fungus Lasiodiplodia theobromae ZJ-HQ1. J. Nat. Prod. 2016, 79, 2397–2402. [Google Scholar] [CrossRef] [PubMed]
  25. Ellman, G.L.; Courtney, K.D.; Andres, V., Jr.; Feather-Stone, R.M. A new and rapid colorimetric determination of acetylcholinesterase activity. Biochem. Pharmacol. 1961, 7, 88–95. [Google Scholar] [CrossRef]
  26. Cheung, J.; Rudolph, M.J.; Burshteyn, F.; Cassidy, M.S.; Gary, E.N.; Love, J.; Franklin, M.C.; Height, J.J. Structures of Human Acetylcholinesterase in Complex with Pharmacologically Important Ligands. J. Med. Chem. 2012, 55, 10282–10286. [Google Scholar] [CrossRef] [PubMed]
  27. Feng, X.; Ambia, J.; Chen, K.M.; Young, M.; Barth, P. Computational design of ligand-binding membrane receptors with high selectivity. Nat. Chem. Biol. 2017, 13, 715–723. [Google Scholar] [CrossRef] [PubMed]
Figure 1. The structures of compounds 18.
Figure 1. The structures of compounds 18.
Marinedrugs 16 00177 g001
Figure 2. Key 2D correlations of compounds 1, 3, 4 and 7.
Figure 2. Key 2D correlations of compounds 1, 3, 4 and 7.
Marinedrugs 16 00177 g002
Figure 3. The experimental and calculated ECD of 1.
Figure 3. The experimental and calculated ECD of 1.
Marinedrugs 16 00177 g003
Figure 4. The experimental ECD of 28.
Figure 4. The experimental ECD of 28.
Marinedrugs 16 00177 g004
Figure 5. ORTEP drawings of compounds 5 and 7.
Figure 5. ORTEP drawings of compounds 5 and 7.
Marinedrugs 16 00177 g005
Scheme 1. Postulated biosynthetic pathway of compounds 1 and 2.
Scheme 1. Postulated biosynthetic pathway of compounds 1 and 2.
Marinedrugs 16 00177 sch001
Figure 6. Low-energy binding conformations of compounds bound to AChE generated by virtual ligand docking. (A) Compound 7 was observed to occupy the active site with significant scores and adopted a conformation similar to that of galanthamine. (B) Compound 7 had the ability to form key hydrogen bonding interaction with residues Ser125, Tyr133, Glu202, and His447. (C,D) Compounds 2 (C) and 5 (D) were observed to occupy the active site with pretty low scores comparing with that of galanthamine.
Figure 6. Low-energy binding conformations of compounds bound to AChE generated by virtual ligand docking. (A) Compound 7 was observed to occupy the active site with significant scores and adopted a conformation similar to that of galanthamine. (B) Compound 7 had the ability to form key hydrogen bonding interaction with residues Ser125, Tyr133, Glu202, and His447. (C,D) Compounds 2 (C) and 5 (D) were observed to occupy the active site with pretty low scores comparing with that of galanthamine.
Marinedrugs 16 00177 g006
Table 1. 1H NMR Data for Compounds 17 (δ in ppm, J in Hz).
Table 1. 1H NMR Data for Compounds 17 (δ in ppm, J in Hz).
No.1 a2 a3 a4 b5 b6 a7 a
13.98, dd (9.6, 2.3)3.98, dd (9.7, 2.0)5.42, dd (11.0, 1.8)3.87, t (7.0)3.79, dd (7.0, 5.2)4.96, d (7.3)7.08, d (10.2)
22.46, dd (15.7, 9.6)
2.68, m
2.39, m
2.63, dd (15.7, 2.0)
2.43, dd (14.8, 3.9)
2.88, dd (14.8, 1.8)
2.84, m
3.33, m
2.79, dd (16.0, 7.0)
3.32, d (16.0)
3.10, dd (16.2, 7.3)
3.55, d (16.2)
5.87, d (10.2)
52.13, d (2.3)2.12, d (2.5)2.49, d (3.0)2.21, dd (12.1, 2.4)2.38, d (2.3)2.49, d (2.6)1.98, d (1.8)
65.48, dd (4.5, 2.3)5.48, m5.27, m1.67, m
1.92, m
5.52, dd (6.3, 3.0)5.64, m5.60, m
72.11, m
2.35, m
2.11, m
2.35, m
2.04, m
2.30, dd (14.4, 2.9)
1.72, m
2.06, m
1.92, dd, (14.6, 6.3)
2.26, dd, (14.6, 3.0)
2.05, m
2.40, dd (14.7, 3.0)
1.91, dd (14.7, 3.0)
2.42, m
91.97, dd (12.6, 5.2)1.98, dd (11.3, 6.4)1.90, dd (12.7, 4.6)2.32, dd (12.5, 4.6)2.34, m2.15, m1.81, dd (12.9, 4.3)
112.33, m
2.41, m
2.40, m
2.41, m
2.40, m
2.80, dd (17.4, 4.6)
2.08, m
2.43, dd (16.2, 4.6)
2.09, m
2.35, m
2.29, m
2.30, m
2.42, dd (16.5, 12.9)
2.76, dd (16.5, 4.3)
155.91, d (1.0)5.92, d (1.1)5.91, d (1.1)5.81, d (1.1)5.72, d (0.6)5.91, d (1.1)5.70, d (1.1)
172.20, d (1.0)2.21, d (1.1)2.21, d (1.1)2.15, d (1.1)2.05, d (0.6)2.20, d (1.1)2.18, d (1.1)
181.43, s1.43, s1.42, s1.30, s1.32, s1.41, s1.40, s
191.31, s1.32, s1.33, s1.14, s1.35, s1.53, s1.46, s
201.26, s1.26, s5.05, m
5.08, m
1.41, s1.31, s1.49, s1.23, s
211.18, s1.18, s1.94, s1.50, s1.57, s1.73, s1.19, s
232.07, s2.07, s2.07, s 2.06, s2.14, s2.09, s
243.70, s 2.08, s
25 1.99, s
a Recorded at 400 MHz in methanol-d4; b recorded at 600 MHz in acetone-d6.
Table 2. 13C and DEPT NMR Data of Compounds 17 (δ in ppm, J in Hz).
Table 2. 13C and DEPT NMR Data of Compounds 17 (δ in ppm, J in Hz).
No.1 a2 a3 a4 b5 b6 a7 a
184.0, CH84.5, CH77.2, CH68.8, CH69.3, CH72.6, CH155.1, CH
237.3, CH237.9, CH236.2, CH239.8, CH239.9, CH236.3, CH2125.1, CH
3173.6, C176.7, C173.2, C171.2, C171.0, C172.9, C203.2, C
479.7, C79.6, C144.4, C84.8, C84.1, C86.7, C45.1, C
561.4, CH61.5, CH51.3, CH50.2, CH49.5, CH51.3, CH53.0, CH
668.3, CH68.3, CH73.4, CH23.3, CH272.6, CH72.7, CH68.8, CH
745.3, CH245.3, CH243.3, CH240.6, CH243.3, CH243.3, CH243.4, CH2
881.5, C81.6, C81.2, C80.9, C79.9, C80.4, C79.3, C
948.1, CH48.1, CH44.2, CH44.8, CH44.4, CH45.1, CH46.9, CH
1046.5, C46.4, C45.6, C45.3, C45.4, C45.1, C39.4, C
1119.4, CH219.4, CH219.7, CH217.5, CH217.2, CH217.6, CH216.9, CH2
1298.5, C98.6, C99.0, C98.2, C98.3, C98.3, C97.3, C
13167.1, C167.2, C167.1, C164.5, C164.5, C164.6, C164.8, C
14162.2, C162.2, C162.2, C160.8, C161.0, C162.3, C160.5, C
15102.2, CH102.0, CH101.8, CH100.6, CH100.6, CH101.7, CH100.3, CH
16165.1, C165.2, C164.8, C163.3, C162.8, C166.8, C162.7, C
1719.5, CH319.5, CH319.5, CH319.5, CH319.5, CH319.5, CH319.7, CH3
1823.0, CH323.0, CH322.1, CH320.9, CH321.9, CH322.0, CH322.0, CH3
1912.4, CH312.5, CH315.3, CH315.2, CH316.2, CH316.0, CH318.1, CH3
2030.5, CH330.6, CH3117.7, CH234.9, CH333.9, CH333.9, CH326.6, CH3
2125.5, CH325.6, CH326.5, CH324.0, CH326.4, CH326.2, CH322.8, CH3
22171.6, C171.7, C171.5, C 170.2, C171.2, C169.7, C
2321.3, CH321.3, CH320.9, CH3 21.4, CH321.3, CH321.5, CH3
2452.3, OCH3 172.2, C 171.1, C
25 21.4, CH3 20.6, CH3
a Recorded at 100 MHz in methanol-d4; b recorded at 150 MHz in acetone-d6.
Table 3. The inhibitory value and docking score of compounds 18 with AChE.
Table 3. The inhibitory value and docking score of compounds 18 with AChE.
CompoundsInhibitory Ratio at 100 μMIC50 Value (μM) aDocking Score b
135.7%>403.25
236.2%>403.89
327.6%>402.76
450.2%>405.32
548.7%>404.98
662.9%>405.87
792.4%13.67.85
825.3%>402.58
Galanthamine 3.57
a IC50 value of compounds against AChE (μM); b Docking score/interaction potential of compounds with target protein.

Share and Cite

MDPI and ACS Style

Li, H.; Sun, W.; Deng, M.; Qi, C.; Chen, C.; Zhu, H.; Luo, Z.; Wang, J.; Xue, Y.; Zhang, Y. Asperversins A and B, Two Novel Meroterpenoids with an Unusual 5/6/6/6 Ring from the Marine-Derived Fungus Aspergillus versicolor. Mar. Drugs 2018, 16, 177. https://doi.org/10.3390/md16060177

AMA Style

Li H, Sun W, Deng M, Qi C, Chen C, Zhu H, Luo Z, Wang J, Xue Y, Zhang Y. Asperversins A and B, Two Novel Meroterpenoids with an Unusual 5/6/6/6 Ring from the Marine-Derived Fungus Aspergillus versicolor. Marine Drugs. 2018; 16(6):177. https://doi.org/10.3390/md16060177

Chicago/Turabian Style

Li, Huaqiang, Weiguang Sun, Mengyi Deng, Changxing Qi, Chunmei Chen, Hucheng Zhu, Zengwei Luo, Jianping Wang, Yongbo Xue, and Yonghui Zhang. 2018. "Asperversins A and B, Two Novel Meroterpenoids with an Unusual 5/6/6/6 Ring from the Marine-Derived Fungus Aspergillus versicolor" Marine Drugs 16, no. 6: 177. https://doi.org/10.3390/md16060177

APA Style

Li, H., Sun, W., Deng, M., Qi, C., Chen, C., Zhu, H., Luo, Z., Wang, J., Xue, Y., & Zhang, Y. (2018). Asperversins A and B, Two Novel Meroterpenoids with an Unusual 5/6/6/6 Ring from the Marine-Derived Fungus Aspergillus versicolor. Marine Drugs, 16(6), 177. https://doi.org/10.3390/md16060177

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop