Terminal Cyclohexane-Type Meroterpenoids from the Fruiting Bodies of Ganoderma cochlear

Eleven new cyclohexane-type meroterpenoids (1, 3–5, 7, 8, 11–15) and four known similar meroterpenoids (2, 6, 9, and 10) were isolated from Ganoderma cochlear. Their structures and absolute configurations at stereogenic centers were elucidated by using HRESIMS, NMR spectroscopy and computational methods. In addition, the structure of the known meroterpenoid, cochlearol G (2), was revised, and the absolute configurations at the stereogenic centers of known meroterpenoids 9 and 10 were determined. All the isolated meroterpenoids were evaluated for their activities against renal fibrosis and triple negative breast cancer, and their insulin resistance. The results of the renal fibrosis study showed that meroterpenoid 11 inhibits over-expression of fibronectin, collagen I and α-SMA. Results of the wound healing study revealed that 4, 6 and 8 significantly inhibit migration of BT549 cells. Observations made in Western blotting experiments showed that 6 decreases the levels of TWIST1 and ZEB1, and increases the level of E-cadherin. Finally, meroterpenoids 7, 9, 11, and 15 significantly up-regulate p-AMPK protein expression in normal L6 myotubes cells.
Introduction<!>General<!>Fungal Material<!>Extraction and Isolation<!>Compound Characterization<!><!>Compound Characterization<!><!>Compound Characterization<!>Renal Fibrosis Activity Assay<!>Cell Culture<!>Cell Viability Assay<!>Wound Healing Assays<!>Western Blot Analysis<!>Statistical Analysis<!>Cell Culture<!>Western Blotting Analysis<!>Results and Discussion<!><!>Results and Discussion<!><!>Results and Discussion<!><!>Results and Discussion<!><!>Results and Discussion<!><!>Results and Discussion<!><!>Conclusion<!>Data Availability Statement<!>Author Contributions<!>Funding<!>Conflict of Interest<!>Publisher’s Note<!>Supplementary Material<!>
<p>Ganoderma is not only a famous Chinese medicine, it is also used globally as a food, in the form of tea, coffee and other beverages, and in syrups and dietary supplements (Wang et al., 2020; Kumar, 2021). Polysaccharides and triterpenoids are representative of the important biologically active components of Ganoderma (Wang et al., 2020). Recent ongoing research studies exploring Ganoderma demonstrated that it contains meroterpenoid components that possess extensive biological activities, such as renal protection and neuroprotection, anti-inflammation, -tumor and -oxidation properties, and analgesic effects (Jiang et al., 2021). These studies led to a deeper understanding of the components of Ganoderma and insight into active ingredients responsible for its traditional medical properties.</p><p>Ganoderma meroterpenoids comprise a class of substances with great potential, not only because they contain a variety of structural subtypes, but also because many members possess a host of biological activities. For instance, in 2009 ganomycin I was found to inhibit HIV-1 protease, and then in 2014 it was shown to have inhibitory effects on the production of monocyte chemotactic protein 1 (MCP-1) and fibronectin. Following these discoveries, ganomycin I was observed to inhibit NSC proliferation in 2015, and then in 2017 it was discovered to display hypoglycemic, hypolipidemic and insulin-sensitizing effects (Dine et al., 2009; Luo et al., 2015; Yan et al., 2015; Wang et al., 2017). As a result, we have been engaged in a program to isolate and identify new Ganoderma meroterpenoids and to assess their unique biological activities. In a previous effort, we found that these substances have inhibitory effects on renal fibrosis (Luo et al., 2015; Meng et al., 2021). In the current investigation, we isolated fifteen terminal cyclohexane-type meroterpenoids (1–15) from Ganoderma cochlear and evaluated their biological activities against renal fibrosis and insulin resistance. Breast cancer is one of the malignant cancer, and the morbidity and mortality is highest in women, with triple negative breast cancer (TNBC) being. Based on the fact that Ganoderma has been used to treat cancer, we also investigated the activities of the meroterpenoids against cells of triple negative breast cancer, which is difficult to treat cancers and has an extremely high mortality rate (Collignon et al., 2016). The results of this study are described below.</p><!><p>Optical rotations were determined using an Anton Paar MCP-100 digital polarimeter. UV and CD spectra were recorded on a Chirascan instrument. NMR spectra were obtained by using a Bruker Avance Ⅲ 600 MHz or 500 MHz spectrometer, with TMS as an internal standard. All NMR chemical shifts are given in ppm. HRESIMS were recorded using a Shimazu LC-20AD AB SCIEX triple TOF 6600+ MS spectrometer (Shimadzu Corporation, Tokyo, Japan). MCI gel CHP 20P (75–150 μm, Mitsubishi Chemical Industries, Tokyo, Japan), C-18 silica gel (40–60 μm; Daiso Co., Japan), and Sephadex LH-20 (Amersham Pharmacia, Uppsala, Sweden) were used for column chromatography. Silica gel (Qingdao Marine Chemical Inc., Qingdao, China) was used for vacuum column chromatography (VLC). Preparative HPLC was carried out using a Saipuruisi chromatograph equipped with a Thermo Hypersil GOLD-C18 column (250 × 21.2 mm, i. d., 5 μm). Semi-preparative HPLC was carried out using a Saipuruisi chromatograph with a YMC-Pack ODS-A column (250 × 10 mm, i. d., 5 μm). Chiral HPLC analysis was run using an Agilent 1260 chromatograph with a Daicel Chiralpak column (IC, 250 mm × 10 mm, i. d., 5 μm) or a Daicel Chiralpak column (IC, 250 mm × 4.6 mm, i. d., 5 μm).</p><!><p>The dried fruiting bodies of G. cochlear were purchased from Tongkang Pharmaceutical Co. Ltd. Guangdong province, China, in July 2014. This fungus was authenticated by Prof. Zhu-Liang Yang at Kunming Institute of Botany, Chinese Academy of Sciences, China, and a voucher specimen (CHYX-0589) is deposited at Institute for Inheritance-Based Innovation of Chinese Medicine, Shenzhen University Health Science Center, China.</p><!><p>Powders of G. cochlear (200 kg) fruiting bodies were extracted with refluxing 80% EtOH (3 × 120 L, 4, 3, 3 h) and the extract was concentrated under reduced pressure to afford a crude residue. An aliquot (8 kg of the residue corresponding to 95 kg fungal material) was suspended in water and extracted three times with EtOAc, followed by concentration of the combined extracts to afford an EtOAc soluble residue (4 kg). The residue was subjected to silica gel column using an eluant comprised of increasing amounts of acetone in petroleum ether to provide four parts (Fr.1–Fr.4). Fr.2 (860 g) was further divided into six parts (Fr.2.1–Fr.2.6) by MCI gel CHP 20P (MeOH/H2O, 60–100%). Fr.2.2 (120.0 g) was divided into five parts (Fr.2.2.1−Fr.2.2.5) by using RP-18 column chromatography (eluting solvent: MeOH/H2O, 40–100%). Fr.2.2.3 (21.0 g) was subjected to a MCI gel CHP 20P column chromatography (MeOH/H2O, 40–100%) to obtain five parts (Fr.2.2.3.1−Fr.2.2.3.5). Fr.2.2.3.3 (3.7 g) was gel filtered through Sephadex LH-20 (MeOH) to produce three parts (Fr.2.2.3.3.1−Fr.2.2.3.3.3.). Fr.2.2.3.3.1 (1.9 g) was subjected to preparative HPLC (MeOH/H2O, 65–100%) to produce three parts (Fr.2.2.3.3.1.1−Fr.2.2.3.3.1.3). Semi-preparative HPLC (eluting solvent: MeOH/H2O containing 0.05% TFA, 60–100%, flow rate: 3 mL/min) of Fr.2.2.3.3.1.2 (0.9 g) gave three subfractions (Fr.2.2.3.3.1.2.1−Fr.2.2.3.3.1.2.3). Semi-preparative HPLC (eluting solvent: MeOH/H2O containing 0.05% TFA, 70%, flow rate: 3 mL/min) of subfraction Fr.2.2.3.3.1.2.3 (300 mg) yielded 12 (5.6 mg, tR = 20.5 min). Semi-preparative HPLC (eluting solvent: MeOH/H2O containing 0.05% TFA, 65%, flow rate: 3 mL/min) of Fr.2.2.3.3.1.3 (200 mg) produced subfraction Fr.2.2.3.3.1.3.3 (13.1 mg, tR = 20.7 min), which by using chiral Daicel Chialpak IC column (n-hexane/ethanol containing 0.05% TFA, 92:8, flow rate: 1 mL/min) yielded 14 (4.0 mg, tR = 13.9 min) and 15 (5.9 mg, tR = 20.0 min).</p><p>Fr.3 (780 g) was fractionated into eight parts (Fr.3.1–Fr.3.8) by using a MCI gel CHP 20P column chromatography (MeOH/H2O, 40–100%). Fr.3.4 (120 g) was divided into three portions (Fr.3.4.1−Fr.3.4.3) by using RP-18 column chromatography (MeOH/H2O, 20%–100%). Submission of Fr.3.4.1 (115 g) to a MCI gel CHP 20P column chromatography (MeOH/H2O, 40–100%) provided six portions (Fr.3.4.1.1–Fr.3.4.1.6). Fr.3.4.1.3 (21.6 g) was separated by using Sephadex LH-20 (MeOH) into four parts (Fr.3.4.1.3.1–Fr.3.4.1.3.4). Fr.3.4.1.3.3 (6.5 g) was subjected to RP-18 (50 µm) column chromatography to obtain five parts (Fr.3.4.1.3.3.1–Fr.3.4.1.3.3.5). Fr.3.4.1.3.3.1 (1.1 g) was gel filtered through Sephadex LH-20 (MeOH) to produce two subfractions (Fr.3.4.1.3.3.1.1 and Fr.3.4.1.3.3.1.2). Fr.3.4.1.3.3.1.1 (0.7 g) was separated into three portions by using preparative TLC (CH2Cl2:MeOH = 7:1). Submission of Fr.3.4.1.3.3.1.1.2 (Rf = 0.6, 100 mg) to semi-preparative HPLC (acetonitrile:MeOH:H2O containing 0.05% TFA, 25.5:25.5:49.0, flow rate: 3 mL/min) yielded 2 (34.8 mg, tR = 18.5 min) and 1 (11.4 mg, tR = 20.0 min). Subjection of Fr.3.4.1.3.3.1.1.3 (Rf = 0.8, 200 mg) to chiral Daicel Chialpak IC column chromatography (n-hexane/ethanol containing 0.05% TFA, 91:9, flow rate: 3 mL/min) yielded 13 (76.9 mg, tR = 16.0 min). Fr.3.4.1.3.3.3 (1.4 g) was filtered through Sephadex LH-20 (MeOH) to produce four fractions (Fr.3.4.1.3.3.3.1–Fr.3.4.1.3.3.3.4). Fr.3.4.1.3.3.3.1 (79 mg) was subjected to semi-preparative HPLC (MeOH:H2O containing 0.05% TFA, 62%, flow rate: 3 mL/min) to yield 11 (7.2 mg, tR = 37.4 min). Fr.3.4.1.3.3.3.3 (0.8 g) was subjected to preparative HPLC (acetonitrile:H2O containing 0.05% TFA, 42%, flow rate: 8 mL/min) to generate six fractions (Fr.3.4.1.3.3.3.3.1–Fr.3.4.1.3.3.3.3.6). Fr.3.4.1.3.3.3.3.2 (178.3 mg) was subjected to chiral Daicel Chiralpak IC column chromatography (n-hexane/isopropanol containing 0.05% TFA, 91:9, flow rate: 3 mL/min) to yield 7 (18.3 mg, tR = 34.6 min) and 8 (18.0 mg, tR = 47.2 min). Fr.3.4.1.3.3.3.3.3 (0.4 g) subjected to chiral Daicel Chiralpak IC column chromatography (n-hexane/ethanol containing 0.05% TFA, 90:10) to produce 9 (8.82 mg, tR = 32.5 min) and 10 (6.15 mg, tR = 24.0 min).</p><p>Fr.3.4.1.3.3.4 (1.0 g) was filtered through Sephadex LH-20 (MeOH) to produce four fractions (Fr.3.4.1.3.3.4.1–Fr.3.4.1.3.3.4.4). Fr.3.4.1.3.3.4.4 (110 mg) was subjected to semi-preparative HPLC (MeOH:H2O containing 0.05% TFA, 60%, flow rate: 3 mL/min) to form four fractions, the fourth fraction of which was subjected to preparative HPLC (acetonitrile:MeOH:H2O containing 0.05% TFA, 28:28:44, flow rate: 12 mL/min) to yield 6 (37.4 mg, tR = 24.2 min).</p><p>Fr.3.4.1.1 (39.4 g) was gel filtered though Sephadex LH-20 (MeOH) and then subjected to C-18 silica gel column chromatography (MeOH/H2O, 20–100%) to produce eight fractions (Fr.3.4.1.1.1–Fr.3.4.1.1.8). Fr.3.4.1.1.5 (10 g) was separated into seven parts (Fr.3.4.1.1.5.1–Fr.3.4.1.1.5.7) by using silica gel column chromatography using an eluant comprised of increasing amounts of methanol in dichloromethane. Fr.3.4.1.1.5.5 (2.0 g) was gel filtered through Sephadex LH-20 (MeOH) to obtain five fractions (Fr.3.4.1.1.5.5.1-Fr.3.4.1.1.5.5.5). Fr.3.4.1.1.5.5.5 (0.6 g) was subjected to semi-preparative HPLC (MeOH:H2O containing 0.05% TFA, 55%, flow rate: 3 mL/min) to produce four subfractions (Fr.3.4.1.1.5.5.5.1–Fr.3.4.1.1.5.5.5.4). Fr.3.4.1.1.5.5.5.2 (116 mg) was subjected to semi-preparative HPLC (acetonitrile:MeOH:H2O containing 0.05% TFA, 19:19:62, flow rate: 3 mL/min) to yield 3 (10.6 mg, tR = 24.2 min). Fr.3.4.1.1.5.5.4 (0.3 g) was divided into two subfractions by using semi-preparative HPLC (MeOH:H2O containing 0.05% TFA, 55%, flow rate: 3 mL/min), the first subfraction was subjected to chiral Daicel Chiralpak IC column chromatography (n-hexane/ethanol containing 0.05% TFA, 90:10) to yield 4 (18.2 mg, tR = 19.7 min) and 5 (17.6 mg, tR = 24.1 min).</p><p>Meroterpenoid 6 (2.5 mg) was also isolated from 40 kg Ganoderma lucidum cultivated in Yongsheng County of Yunnan Province, China. This material was also authenticated by Prof. Zhu-Liang Yang at Kunming Institute of Botany, Chinese Academy of Sciences, China. A voucher specimen (CHYX-0609) is deposited at Institute for Inheritance-Based Innovation of Chinese Medicine, Health Science Center, Shenzhen University, China.</p><!><p>Ganodercin G (
1
), yellow gum; [α]D
20 –14.5 (c 0.06, MeOH); UV (MeOH) λ
max (log ε) 362 (3.40), 201 (4.29) nm; CD (MeOH) Δ
ε
202–1.83, HRESIMS m/z 375.1812 [M + H]+, (calcd for C21H27O6, 375.1802), 1H and 13C NMR data, see Tables 1 and 2.</p><!><p>1H NMR (500 MHz) data of 1, 3–5, 7, and 8 (δ in ppm, J in Hz).</p><p>Record in methanol-d4.</p><p>Record in acetone-d6.</p><p>13C NMR data of 1, 3–5, 7–15 (δ in ppm).</p><p>Record in methanol-d4 at 125 MHz.</p><p>Record in acetone-d6 at 125 MHz.</p><p>Record in methanol-d4 at 150 MHz.</p><!><p>Ganodercin H (

3

), yellow gum; [α]D
20 +10.6 (c 0.09, MeOH); CD (MeOH) Δ
ε
361–1.52, Δ
ε
211 +1.18; UV (MeOH) Î»
max (log ε) 361 (3.52), 261 (3.74), 216 (4.26) nm, HRESIMS m/z 375.1804 [M + H]+, (calcd for C21H27O6, 375.1802), 1H and 13C NMR data, see Tables 1 and 2.</p><p>Ganodercin I (
4
), yellow gum; [α]D
20 +37.6 (c 0.11, MeOH); UV (MeOH) Î»
max (log ε) 297 (3.66), 223 (4.19), 202 (4.33) nm, CD (MeOH) Δ
ε
299 +0.83, Δ
ε
260–0.31, Δ
ε
216 +4.43, HRESIMS m/z 361.2012 [M + H]+, (calcd for C21H29O5, 361.2010), 1H and 13C NMR data, Tables 1 and 2.</p><p>Ganodercin J (
5
), yellow gum; [α]D
20 –5.9 (c 0.10, MeOH); UV (MeOH) λ
max (log ε) 296 (3.54), 222 (3.83), 201 (4.31) nm, CD (MeOH) Δ
ε
204 +2.89, HRESIMS m/z 361.2008 [M + H]+, (calcd for C21H29O5, 361.2010), 1H and 13C NMR data, see Tables 1 and 2.</p><p>Ganodercin K (
7
), yellow gum; [α]D
20 –35.0 (c 0.08, MeOH); CD (MeOH) Δ
ε
257–0.40, Δ
ε
231–1.38, Δ
ε
215 +1.62, Δ
ε
200–4.86; UV (MeOH) λ
max (log ε) 364 (3.55), 256 (3.80), 225 (4.14), 201 (4.21) nm, HRESIMS m/z 377.1964 [M + H]+, (calcd for C21H29O6, 377.1959), 1H and 13C NMR data, see Tables 1 and 2.</p><p>3-epi-Ganodercin K (
8
), yellow gum; [α]D
20 –5.0 (c 0.02, MeOH); CD (MeOH) Δ
ε
259 + 0.32, Δ
ε
227 +1.38, Δ
ε
216 +2.18, Δ
ε
201–3.21; UV (MeOH) λ
max (log ε) 364 (3.61), 256 (3.87), 225 (4.21), 201 (4.30) nm, HRESIMS m/z 377.1964 [M + H]+, (calcd for C21H29O6, 377.1959), 1H and 13C NMR data, see Tables 1 and 2.</p><p>Ganodercin L (
9
), yellow gum; [α]D
20 +12.0 (c 0.03, MeOH); CD (MeOH) Δ
ε
255 +0.38, Δ
ε
228 +1.48, Δ
ε
205 +0.75, Δ
ε
196–1.02; UV (MeOH) λ
max (log ε) 363 (3.48), 256 (3.72), 226 (4.00), 202 (4.01) nm, HRESIMS m/z 377.1961 [M + H]+, (calcd for C21H29O6, 377.1959), 1H and 13C NMR data, see Tables 2 and 3.</p><!><p>1H NMR data of 9–15 in methanol-d
4 (δ in ppm, J in Hz).</p><p>Record in acetone-d6 at 500 MHz.</p><p>Record in methanol-d4 at 500 MHz.</p><p>Record in methanol-d4 at 600 MHz.</p><!><p>3-epi-Ganodercin L (
10
), yellow gum; [α]D
20 –17.3 (c 0.08, MeOH); CD (MeOH) Δ
ε
256–0.13, Δ
ε
228–0.49, Δ
ε
201 +0.31, Δ
ε
192–1.74; UV (MeOH) λ
max (log ε) 363 (3.35), 256 (3.59), 227 (3.88), 201 (3.89) nm, HRESIMS m/z 377.1964 [M + H]+, (calcd for C21H29O6, 377.1959), 1H and 13C NMR data, see Tables 2 and 3.</p><p>Ganodercin M (
11
), yellow gum; [α]D
20 –16.7 (c 0.03, MeOH); CD (MeOH) Δ
ε
256–0.18, Δ
ε
230–0.70, Δ
ε
202 +0.66, Δ
ε
192–1.75; UV (MeOH) λ
max (log ε) 363 (3.28), 257 (3.58), 227 (3.88), 201 (3.89) nm, HRESIMS m/z 391.2120 [M + H]+, (calcd for C22H31O6, 391.2115), 1H and 13C NMR data, see Tables 2 and 3.</p><p>Ganodercin N (
12
), yellow gum; [α]D
20 +22.5 (c 0.04, MeOH); UV (MeOH) Î»
max (log ε) 386 (3.38), 267 (3.91), 227 (3.89), 200 (4.09) nm, HRESIMS m/z 375.1812 [M + H]+, (calcd for C21H27O6, 375.1802), 1H and 13C NMR data, see Tables 2 and 3.</p><p>Ganodercin O (
13
), yellow gum; [α]D
20 +10.0 (c 0.06, MeOH); CD (MeOH) Δ
ε
228–1.19, Δ
ε
216–2.12, Δ
ε
202 +1.65; UV (MeOH) Î»
max (log ε) 364 (3.49), 219 (4.34) nm, HRESIMS m/z 375.1805 [M + H]+, (calcd for C21H27O6, 375.1802), 1H and 13C NMR data, see Tables 2 and 3.</p><p>Ganodercin P (
14
), yellow gum; [α]D
20 +24.0 (c 0.03, MeOH); CD (MeOH) Δ
ε
254 +0.28, Δ
ε
229 +0.83; UV (MeOH) Î»
max (log ε) 363 (3.45), 256 (3.71), 226 (4.00) nm, HRESIMS m/z 377.1965 [M + H]+, (calcd for C21H29O6, 377.1959), 1H and 13C NMR data, see Tables 2 and 3.</p><p>3-epi-Ganodercin P (
15
), yellow gum; [α]D
20 –4.0 (c 0.05, MeOH); CD (MeOH) Δ
ε
256–0.33, Δ
ε
226–1.24; UV (MeOH) λ
max (log ε) 363 (3.38), 256 (3.63), 227 (3.92) nm, HRESIMS m/z 377.1968 [M + H]+, (calcd for C21H29O6, 377.1959), 1H and 13C NMR data, see Tables 2 and 3.</p><p>Crystal Data for C21H28O6 (M = 376.43 g/mol): monoclinic, space group P21 (no. 4), a = 6.62350(10) Å, b = 11.59270(10) Å, c = 13.40440(10) Å, β = 100.4300(10)°, V = 1012.240(19) Å3, Z = 2, T = 99.9(9) K, μ(CuKα) = 0.737 mm−1, Dcalc = 1.235 g/cm3, 9814 reflections measured (6.704° ≤ 2Θ ≤ 148.606°), 3973 unique (R
int = 0.0305, R
sigma = 0.0308) which were used in all calculations. The final R
1 was 0.0324 (I > 2σ(I)) and wR
2 was 0.0865 (all data). Crystallographic data of 3-epi-ganodercin P (15) have been deposited at the Cambridge Crystallographic Data Centre (CCDC deposition no. 2105988).</p><!><p>TGF-β1 induced rat renal proximal tubular cells (NRK-52E) were used to assess expression of the target gene. The cell culture method, and cell viability and western blotting assays were conducted following our previously reported protocols (Meng et al., 2021).</p><!><p>MDA-MB-231 (ATCC, Cat No. HTB-26, United States), BT549 (ATCC, Cat No. HTB-122, United States) and HCC1806 (ATCC, Cat No. CRL-2335, United States) cells were cultured in Dulbecco's Modified Eagle Medium (HyClone, Cat No.SH30243.01, United States) or RPMI Medium Modified (HyClone, Cat No.SH30809.01, United States) supplemented with 10% fetal bovine serum (HyClone, Cat No.SV30160.03, United States) and 1% penicillin/streptomycin (Hyclone, Cat No.SV30010, United States) in CO2 (5%) incubator at 37°C.</p><!><p>Cell viability was evaluated by using a CCK-8 assay kit (meilunbio, Cat No. MA0218, China) according to the manufacturer's instructions. MDA-MB-231, BT549 and HCC1806 cells were seeded into 96-well plates with 3 × 103 cells per well. After 72 h exposure in the medium containing desired compounds (20 μM), 100 μL fresh medium including 10 μL CCK-8 reagent was added, and then the plates were incubated at 37°C for 2 h. The absorbances of CCK-8 in each well was measured at 450 and 600 nm by using a Cytation5 (BioTek, United States). DMSO was used as a control. Each sample was plated in triplicate.</p><!><p>Confluent MDA-MB-231, BT549, and HCC1806 cells were wounded by scratching using Wound Making Tool-Auto Scratch (BioTek, United States). After exposure for 24 h in the medium containing desired compounds (20 μM), the scratched gap area of each cell monolayer was photographed by using Cytation5 (BioTek, United States) and quantified by using Image-Pro Plus 6.0. (http://rsb.info.nih.gov/ij/download.html). The migration efficiency of each cell was then determined as average percentage of closure of the scratch area. Each sample was plated in triplicate. All product recommended protocols were followed.</p><!><p>The total proteins of the whole cell lysates were separated by using SDS-PAGE and transferred to a PVDF membrane (Millipore, Cat No. IPVH00010, United States). After being probed with primary antibodies overnight at 4°C, the proteins of interest were then detected using HRP-conjugated IgG (CST, Cat No.7074P2, United States) and visualized by using ECL substrate (4ABiotech, Cat No.4AW011, China) based imaging with a Minichemi™ chemiluminescence imaging system (SAGECREATION, China).</p><!><p>Analysis of statistical data, obtained from triplicate measurements, was performed by using the Student's t-test for two groups or by one-way ANOVA for multiple groups. #p < 0.05 was considered to be significant.</p><!><p>L6 myotubes cells were maintained in α-MEM culture medium supplemented with 10% (v/v) fetal bovine serum (FBS), 100 U/mL penicillin-streptomycin, and incubated at 37°C in an atmosphere of 5% CO2.</p><!><p>Cells were washed with pre-cold PBS, lysed using RIPA buffer supplemented with Proteinase Inhibitors Cocktail (MCE), and PMSF (MCE). AMPK, p-AMPK, AKT, p-AKT and GAPDH antibody were incubated overnight and the secondary antibody was incubated for 1 h at room temperature. Signals were detected using a Western blotting Imaging System according to the manufacturer's specifications. The following primary antibodies were used for blotting: p-AKT (4060S, CST), AKT (9272S, CST), p-AMPK (2535S, CST), AMPK (2532, CST) and GAPDH (5174S, CST).</p><!><p>Powders of the dried fruiting bodies of the fungus G. cochlear (200 kg) were first extracted with refluxing 80% EtOH and partitioned in water to obtain an ethyl acetate soluble fractions, an qliquote was then subjected to multiple chromatographic separation steps. This procedure led to isolation of fifteen the cyclohexane-type meroterpenoids 1–15 (Figure 1).</p><!><p>Structures of 1–15 isolated from G. cochlear.</p><!><p>Ganodercin G (1), obtained as a yellow gum, has a molecular formula of C21H26O6 (nine degrees of unsaturation) established by using HRESIMS (m/z 375.1812 [M + H]+, calcd for 375.1802), and 13C NMR and DEPT spectroscopy. The presence of a 1,2,4-trisubstituted benzene ring in 1 was assigned by analysis of the 1H NMR spectrum of 1 (Table 1), which contain resonances in the aromatic region (δ
H 7.37, d, J = 2.2 Hz, H-3′; δ
H 7.02, dd, J = 8.9, 2.2 Hz, H-5′; δ
H 6.81, d, J = 8.9 Hz, H-6′). The 13C NMR and DEPT data of 1 (Table 2) revealed the presence of three methyls, four sp3 methylenes, five methines (four sp2 and one sp3) and nine nonprotonated carbons (one ketone, one carbonyl, six aromatic including two oxygenated and one aliphatic). The NMR spectra of 1 are similar to those of (–)-ganotheaecoloid F (Luo et al., 2018), which was isolated from the fruiting bodies of G. theaecolum, with the main difference being the presence of a Δ
3(4) double bond in 1 instead of the Δ
2(3) double bond in (–)-ganotheaecoloid F. This difference is supported by the 1H–1H COSY correlation of H-4 (δ
H 6.90)/H2-5 (Î´
H 2.97 and 2.91) and the HMBC (Figure 2) correlations of H2-2 (Î´
H 4.12)/C-3 (Î´
C 126.6), C-15 (Î´
C 170.5), H-4/C-15, and H2-5/C-3. The ROESY correlation (Figure 3) of H-2/H-5 suggested that the Δ
3(4) double bond in 1 has E configuration. Chiral HPLC analysis revealed that compound 1 is enantiomerically pure. Time-dependent density functional theory (TDDFT) ECD calculations were used to determine the absolute configuration at the stereogenic center in 1. It was found that the calculated ECD curve of (10S)-1 obtained at the CAM-B3LYP/def2SVP level matches well with the experimental CD spectrum (Figure 4).</p><!><p>Key COSY and HMBC correlations of 1, 3–5, 7–15.</p><p>Key ROESY correlations of 1, 3–5, 9–15.</p><p>Comparison of the calculated ECD and experimental CD spectra in MeOH. (A): The calculated ECD spectrum of (10S)-1 at CAM-B3LYP/def2SVP level, σ = 0.20 eV; shift = –13 nm. (B): The calculated ECD spectrum of (6R,10S)-4 at CAM-B3LYP/def2SVP level, σ = 0.28 eV; shift = +8 nm. (C): The calculated ECD spectrum of (6R,10S)-5 at B3LYP/6-31G (d,p) level, σ = 0.25 eV; shift = +20 nm. (D): The calculated ECD spectrum of (3S,10S)-7a at B3LYP/6-31G (d,p) level, σ = 0.25 eV; shift = +20 nm and calculated ECD spectrum of (3R,10S)-7b at B3LYP/6-31G (d,p) level, σ = 0.23 eV; shift = +20 nm. (E): The calculated ECD spectrum of (3S,6R,10S)-9a at CAM-B3LYP/def2SVP level, σ = 0.28 eV; shift = +13 nm and calculated ECD spectrum of (3R,6R,10S)-9b at B3LYP/6-31G (d,p) level, σ = 0.20 eV; shift = +16 nm. (F): The calculated ECD spectrum of (6S,6R,10S)-13 at CAM-B3LYP/def2SVP level, σ = 0.20 eV; shift = –2 nm.</p><!><p>Cochlearol G (2) has been previously described as a component of fruiting bodies of G. cochlear by Wang and co-workers (Wang et al., 2019a). Although the ROESY experiment performed by Wang et al. led to assignment of the relative configurations of the two stereogenic centers in the terminal six-member ring of 2, the configuration of the Δ
3(4) double bond remained unresolved. As a result, we carefully analyzed the ROESY spectrum and found the existence of a correlation H2-2 (Î´
H 4.16)/H2-5 (Î´
H 2.65 and 2.42). This finding demonstrated that the Δ
3(4) double bond in 2 has E configuration, rather than the Z configuration depicted by Wang et al. (Wang et al., 2019a).</p><p>Ganodercin H (3) was isolated as a yellow gum. The HRESIMS of 3 gave a molecular ion m/z 375.1804 [M + H]+ consistent with the molecular formula C21H26O6, indicating nine degrees of unsaturation. Detailed analysis of the 1D and 2D NMR data of 3 and ganotheaecoloid D shows that the E Δ
3(4)-double bond is present in 3 rather than Z-configuration found in ganotheaecoloid D (8). This conclusion is supported by the ROESY correlation in 3 of H-2 (Î´
H 4.09)/H-5 (Î´
H 2.45). The relative configurations of their stereogenic centers in the terminal cyclohexane ring in 3 were assigned as 6R*,10S* by using ROESY correlation of H-6 (Î´
H 1.94)/H-10 (Î´
H 3.38). Moreover, the absolute configurations at these centers were assigned as 6R,10S by comparing its CD spectrum to that of ganotheaecoloid D.</p><p>Ganodercin I (4), isolated as a yellow gum, has the molecular formula C21H28O5 (eight degrees of unsaturation) based on HRESIMS analysis (m/z 361.2012 [M + H]+; calcd for C21H29O5
+, 361.2010). 1D and 2D NMR analysis shows that 4 has a structure that is similar to that of 3. The presence of resonances associated with a terminal double bond (Î´
H 4.84, Ha-14, Î´
H 4.61, Hb-14) and two singlet methyls (Î´
H 0.97, H3-12, Î´
H 0.64, H3-13) in 1H NMR spectrum reveal that 4 contains the same terminal cyclohexane part that is present in 3. The differences between 3 and 4 are that a methylene appears at C-1 in the latter instead of a ketone in the former and the double bond Δ
3(4) of 3 is at Δ
2(3) in 4. The 1H-1H COSY correlation of H2-1 (Î´
H 3.73 and 3.54)/H-2 (Î´
H 5.89) and the HMBC correlations of H-3′ (δ
H 6.58)/C-1 (Î´
C 31.7), H2-1/C-3 (Î´
C 133.6), H-2/C-15 (Î´
C 172.3) and H2-4/C-2 (Î´
C 140.5), C-3, C-15 support the above conclusions. The relative configurations of the two stereogenic centers in 4 were assigned as 6R*,10S* based on ROESY correlation of H-6 (Î´
H 1.66)/H-10 (Î´
H 3.29). Further examination of ROESY correlation of H-2/H-4 indicate that the Δ
2(3) double bond is Z-configuration. Finally, the experimental CD spectrum matched the calculated (CAM-B3LYP/def2SVP level) spectrum of (6R,10S)-4.</p><p>Ganodercin J (5) has the same molecular formula as 4 (HRESIMS). Comariasons of 1H, 13C NMR, and DEPT spectra of 5 and 4 suggest they have similar structures. The difference between 5 and 4, inferred from 1D NMR data, is that the double bond on the terminal ring in 5 is at the Δ
7(8), whereas in compound 4 it is at Δ
7(14). This conclusion is supported by the observation of HMBC correlations of H3-14 (Î´
H 1.70)/C-6 (Î´
C 50.2), C-7 (Î´
C 137.8), C-8 (Î´
C 120.1) and H-8 (Î´
H 5.22)/C-6. The configuration of the Δ
2(3) double bond was determined as Z by observing a ROESY correlation of H-2/H2-4. Moreover, the observation of ROESY correlation of H-6 (Î´
H 1.66)/H-10 (Î´
H 3.36) indicated that 5 is the 6R*,10S* diastereomer. The absolute configurations at the stereogenic centers in 5 were assigned as 6R 10S by comparing experimental CD spectrum with calculated ECD curve (Figure 4).</p><p>Two isolated meroterpenoids 7 (ganodercin K, yellow gum) and 8 (3-epi-ganodercin K, yellow gum) were found to have the same molecular formula (C21H28O6) and nearly identical NMR data, suggesting that they have the same planar structure with both lacking unsaturation in chains connecting the aryl and cyclohexane moieties. Indeed, a detailed comparison of their 1D NMR data with those of 1 shows that 7 and 8 lack the Δ
3(4) double bond in 1, a conclusion supported by 1H–1H COSY correlations of H-2/H-3/H-4/H-5, and the HMBC correlations of H-2/C-1, C-3, C-4, C-15 and H-3/C-4, C-15. It was likely that 7 and 8 are diastereomers having 3S*,10S* and 3R*,10S* relative configurations at their stereogenic centers. Analysis of the CD spectra of 7 and 8 enabled the use of computational methods to elucidate the absolute configurations at their stereogenic centers. As seen by viewing Figure 4, the experimental CD spectrum of 7 well matches the calculated ECD curve of 3R,10S. In addition, the calculated ECD spectrum of 3S,10S matches the experimental CD spectrum of 8.</p><p>Compounds 9 and 10 were previously isolated from G. theaecolum by Luo et al. (Luo et al., 2018) and characterized the enantiomers of ganotheaecoloid C. However, the absolute configurations at the three stereogenic centers in these substances were not determined in the earlier effort. As a result, we used TDDFT-ECD calculations to determine their absolute configurations. It can be seen from viewing Figure 4, that the experimental CD spectrum of 9 agrees well with the calculated ECD curve of the 3S,6R,10S stereoisomer. Importantly, the experimental CD spectrum of 10 is in close agreement with the calculated ECD spectrum of the 3R,6R,10S stereoisomer, showing that 9 and 10 are actually epimers rather than enantiomers. By carefully analyzing the 1D NMR data reported by Luo, it was found that pairs of signals are present in the 13C NMR spectrum, supporting our conclusion. The revised structures of 9 and 10 were renamed as ganodercin L for 9 and 3-epi-ganodercin L for 10.</p><p>Ganodercin M (11) was obtained as a yellow gum and has the molecular formula C22H30O6 (HRESIMS ion observed at m/z 391.2120 [M + H]+, calcd for 391.2115). The NMR data of 11 are similar to those of 9, except that the free C-15 carboxylic acid group in 9 is a methyl ester in 11. This conclusion is confirmed by observations of the HMBC correlation of OCH
3/C-15. The ROESY correlation of H-6/H-10 showed that the relative configurations at the stereogenic centers in the terminal ring in 11 are 6R*,10S*. The computational methods applied to 9 and 10 were used to determine the absolute configurations at the stereogenic centers in 11. Matching experimental and calculated CD spectra showed that 11 is the 3S,6R,10S stereoisomer.</p><p>The NMR spectroscopic data of ganodercin N (12) (yellow powder, C21H26O6) are similar to those of (–)-ganotheaecoloid F (Luo et al., 2018). One difference is that a methine (δ
C 57.1) and a nonprotonated carbon (Î´
C 88.6) are present in 12 instead of the double bond (Î´
C 136.7 and Î´
C 128.6) in (–)-ganotheaecoloid F. This proposal is supported by the HMBC correlations of H3-14 (δ
H 1.28)/C-6 (Î´
C 57.1), C-7 (Î´
C 88.6), C-8 (Î´
C 39.7). Also, the downfield chemical shift of C-10 at 87.6 ppm, indicated that C-10 is the oxygenated. Moreover, the HMBC correlation of H-10/C-7 indicated that C-10 is linked to C-7 via an oxygen bridge. The relative configurations at the stereogenic centers in the terminal bicyclic ring in 12 were determined by using ROESY data. ROESY correlation of H-6 (δ
H 1.26)/H3-12 (Î´
H 1.02) suggested that both H-6 and CH3-12 have a β-orientation. Accordingly, the ROESY correlation of H-10 (δ
H 3.70)/H3-13 (Î´
H 0.97) revealed that H-10 and CH3-13 have a α-orientation. The C-14 methyl was determined to have a α-orientation, when consideration is given to the existence of the oxygen bridge. This proposal was further supported by the ROESY correlation of H2-4 (δ
H 2.58 and 2.52)/H3-14. Thus, the relative configurations at the stereogenic centers in 12 were assigned as 6S*,7R*,10S*. Since no strong correlation of H-2/H-4 was observed, we hypothesized that the Δ
2(3) double bond in 12 is E. Density functional theory (DFT, B3LYP/6-311G(d,p) level) was used to verify this proposal. The experimental NMR data of 12 were compared with those calculated for 12a (E-isomer) and 12b (Z-isomer). It was found that the calculated 13C NMR data for 12a possesses the highest R
2 value and a 100% probability in DP4+ analysis. Because no observable Cotton effects are present in the experimental CD spectrum of 12, it was not possible to use ECD calculations to assign stereochemistry. However, the absolute configuration at the stereogenic centers in 12 were assigned as 6S,7R,10S by comparing the calculated specific optical rotation +22.2 for (6S,7R,10S) with the experimental one of 12 ([α]D
20+22.5).</p><p>Ganodercin O (13) has the same molecular formula as 12 and similar NMR data, differing only by resonances associated with the location of the double bond. The position of double bond in 13 was demonstrated to be Δ
3(4) by the observation of 1H–1H COSY correlations of H-4/H-5/H-6, and the HMBC correlations of H-2/C-1, C-3, C-4, C-15 and H-4/C-15. The double bond Δ
7(8) was deduced to have the E-configuration by existence of the ROESY correlation of H-2/H-4. Similar to that in 12, H-6 and CH3-12 were determined to have a β-orientation by presence of the ROESY correlation of H-6/H3-12. Moreover, H-10 and CH3-14 were assigned to have a α-orientation by the ROESY correlations between H-10/H3-13, H2-5/H-13, H-14. As a result, the relative configurations of the stereogenic centers in 13 were determined to be 6S*,7R*,10S*. Finally, a match between the computed [B3LYP/6-31 g(d,p)] ECD spectrum of 6S,7R,10S-
13 and the experimental one provided the absolute configurations at the stereogenic centers in 13 (Figure 1).</p><p>Ganodercin P (14) and 3-epi-ganodercin P (15) were found to have the same molecular formula of C21H28O6 by using HRESIMS. Careful analysis revealed that their NMR data were similar, suggesting that they have the same planar structures. The observations of 1H–1H COSY correlations of H-2/H-3/H-4, and the HMBC correlations of H-2/C-1, C-15 and H-3/C-15 clearly revealed that one methylene (C-2) and one methine (C-3) in both 14 and 15 replace the E-double bond in 12. Because 14 and 15 contain the same terminal bicyclic ring that is present in 12 and 13, the relative configurations at the stereogenic centers were determined to be 6S*,7R*,10S* by the presence of ROESY correlations H-6/H3-12, H-10/H3-13, H2-5/H3-13, H3-14. The relatively configuration of C-3 is difficult to determined owing to the flexible nature of the chain. As a result, both 14 and 15 can be the two possible diastereomers 3R*,6S*,7R*,10S* and 3S*,6S*,7R*,10S*. Computational methods were used to determine the absolute configurations at the stereogenic centers in 14 and 15. The calculated curves of these substances showed that the terminal oxygen bridge doe not contribute to the Cotton effects (Figure 5). As a result, the absolute configuration of C-3 in 14 is R and that of compound 15 is S.</p><!><p>Comparison of the calculated ECD spectrum of (3R,6S,7R,10S)-14a at CAM-B3LYP/def2SVP level, σ = 0.20 eV, shift = +10 nm; and the calculated ECD spectrum of (3S,6S,7R,10S)-14b at B3LYP/6-31G (d,p) level, σ = 0.20 eV, shift = +10 nm; with the experimental CD spectra of 14 and 15 in MeOH (A). Ortep plot of the X-ray crystallographic data of 15. Anisotropic displacement ellipsoids display 50% probability levels (B).</p><!><p>X-ray crustallographic analysis of 15 was carried out, which was generated by crystallization from Kappa single diffractometer in CuKα with a Flack parameter of –0.08 (7). Analysis of the crystal data leads to assignment of the absolute configurations at the stereogenic centers in 15 to be 3S,6S,7R,10S. Since 14 and 15 are C-3 epimers, the absolute configurations at the centers in 14 are 3R,6S,7R,10S.</p><p>The final substance isolated was (−)-ganotheaecoloid F (6), which was identified by comparison its spectroscopic data with those reported in the literature (Luo et al., 2018).</p><p>In this study, we have determined the structures and absolute configurations at the stereogenic centers in 15 meroterpenoids isolated from dried fruiting bodies of G. cochlear. It is worthy noting that all of these substances contain a terminal cyclohexane ring and all are enantiomerically pure. Meroterpenoids of this type isolated from Ganoderma have been also found to exist in enantiomerically pure (Luo et al., 2018; Wang et al., 2019a; Wang et al., 2019b). Among them, the absolute configuration of four structures was successfully determined. The absolute configurations at the centers in these substances as well as the fact that C-10 in all have the S-configuration might be related to their biosynthetic origin, and might aid subsequent structural identification of analogues.</p><p>The protective activity of the meroterpenoids isolated in this effort were assessed by observing the expression of renal fibroblast biomarkers in TGF-β1–induced NRK-52E, which plays an important role in stimulation of renal fibroblasts (Zhang et al., 2012; Meng et al., 2015). The results showed that 11 displays selective inhibitory activity in that it significantly inhibits over-expression of fibronectin, collagen I and α-SMA at the protein level at a concentration of 40 μM (Figures 6A,B). In addition, western blot analysis, carried out in a dose-concentration dependent manner, demonstrated that the optimum activity of 11 is 40 μM (Figure 6D). To explore the mechanism underlying the antifibrotic effect of, Smad2/3 phosphorylation was investigated. We found that 11 has no effect on phosphorylation of Smad2 or Smad3 in TGF-β1-induced NRK-52E cells at concentrations of 10 μM, 20 and 40 μM (Figure 6C), suggesting that its effects could play potential roles in renal fibrosis through a non-Smad pathway.</p><!><p>Compounds 1–15 inhibit renal fibrosis in vitro. (A,B): Protein expression of collagen I, fibronectin, and α-SMA was examined by Western blot analysis. (C): Compound 11 at 0 μM, 10 μM, 20 μM, and 40 μM affected TGF-β1-induced Smad phosphorylation. GAPDH was used as an internal control. Three independent experiments were performed. (D): Compound 11 dose-dependently inhibits fibronectin and collagen I expression. GAPDH was used as an internal control. GW was used as a positive control. GW: GW788388. Three independent experiments were performed.</p><!><p>We investigated the cellular phenotype promoted by the isolated substances in breast cancer cells by using a cell viability and a wound healing assay in three TNBC cell lines including MDA-MB-231, BT549 and HCC1806. Inspection of the plots given in Figure 7A shows that treatment with 7 and 9 results in significant decreases in cell viability of 29.2 and 26.3% in MDA-MB-231 and HCC1806 cells, respectively. The other meroterpenoids have negligible inhibitory effects on cell viability even at concentrations as high as 20 μM (Figure 7A). Interestingly, although all the fifteen substances have rather low cytotoxicities toward these three breast cancer cell lines, three of the isolates including 4, 6 and 8 significant inhibit the migration ability of BT549 cells. All of the other meroterpenoids display no significant effects on MDA-MB-231 and HCC1806 cells (Figures 7B,C). Moreover, we observed that treatment with 4 and 8 markedly decreases the protein level of TWIST1 and ZEB1 without noticeably affecting the E-cadherin level in BT549 cell lysates. Meroterpenoid 6 not only decreases the protein level of TWIST1 and ZEB1, it also increases the protein level of E-cadherin (Figure 7D). ZEB1, E-cadherin and TWIST1 are generally acknowledged to be transcriptional factors driving the Epithelial–Mesenchymal Transition (EMT), one of the most important pathogenic events occurring in the initiation of cancer metastasis. Thus, the data indicate that 4, 6 and 8 suppress the metastatic potential of TNBC cells through down regulation of EMT, and consequently they are promising lead compounds for the development of the anti-cancer drugs against metastasis of TNBC.</p><!><p>(A) Cell viabilities of MDA-MB-231, BT549 and HCC1806 were analyzed by using CCK8. All compounds concentrations were 20 μM; (B) Quantification of the migratory ability of MDA-MB-231, BT549 and HCC1806 cells. All compounds concentrations were 20 μM; (C) Representative images of the wound healing assay in BT549 cells; (D) Compounds 4, 6 and 8 (40 μM) treated BT549 celllysates were subjected to immunobloting with the indicated antibodies respectively. The results were represented as mean ± S.D. of biological triplicates. #p < 0.05.</p><!><p>All the isolated meroterpenoids at a concentration of 20 μM were exposed to normal L6 myotubes cells for 24 h. The control group was treated with insulin (INS, 100 nM, 15 min) and Berberine (BBR, 10 μM, 24 h). Notably, INS and BBR treatment led to significant increases in the phosphorylation AKT and AMPK, respectively. Compared with the control group, the meroterpenoids have no significant impact on the AKT pathway. Meanwhile, 7, 9, 11 and 15 significantly up-regulate p-AMPK protein expression (Figure 8). This finding suggests that these substances might enhance insulin sensitivity by activating AMP-activated protein kinase (AMPK) in normal L6 myotubes cells.</p><!><p>The effects of compounds on L6 myotubes. (A–F): L6 myotubes were treated with vehicle or 20 μM compounds or 10 μM BBR for 24 h (n = 3). Data are for three replicate experiments. (n = 3) *p < 0.05, **p < 0.01 ***p < 0.001, one-way ANOVA.</p><!><p>In conclusion, the study described above resulted in the isolation of eleven new meroterpenoids, and four known meroterpenoids. The structure of cochlearol G (2) was revised, and the absolute configurations at the stereogenic centers in 9 and 10 were determined by using ECD calculations. Biological studies related to renal fibrosis showed that (1) 11 inhibits over-expression of fibronectin, collagen I and α-SMA, (2) 4, 6 and 8 significantly inhibitof the migration ability of BT549 cells, (3) 6 decreases the protein level of TWIST1 and ZEB1 and increases the protein level of E-cadherin, and (4) 7, 9, 11 and 15 significantly up-regulate p-AMPK protein expression in normal L6 myotubes cells.</p><!><p>The original contributions presented in the study are included in the article/Supplementary Material, further inquiries can be directed to the corresponding authors.</p><!><p>YC designed the research. FQ conducted chemical experiments. TX, YL, and DC conducted biological experiments in vitro. FQ, TX, HZ, LL, and YC analyzed data. FQ and YC wrote and revised the manuscript. All authors discussed the results and commented on the manuscript at all stages.</p><!><p>We are indebted to the NSFC-Joint Foundation of Yunnan Province (U1702287), National Natural Science Foundation of China (82030115), Shenzhen Fundamental Research Program (JCYJ20200109113803838), National Science Fund for Distinguished Young Scholars (81525026), Guangdong Key Laboratory for Functional Substances in Medicinal Edible Resources and Healthcare Products (2021B1212040015), and National Natural Science Foundation of China (82104036).</p><!><p>The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.</p><!><p>All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.</p><!><p>The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fchem.2021.783705/full#supplementary-material</p><!><p>Click here for additional data file.</p><p>Click here for additional data file.</p><p>Click here for additional data file.</p><p>Click here for additional data file.</p>