Next Article in Journal
mRNA, a Revolution in Biomedicine
Next Article in Special Issue
Preparation of Liposomal Formulations for Ocular Delivery of Thymoquinone: In Vitro Evaluation in HCEC-2 e HConEC Cells
Previous Article in Journal
New Properties and Mitochondrial Targets of Polyphenol Agrimoniin as a Natural Anticancer and Preventive Agent
Previous Article in Special Issue
Cyclodextrin as Functional Carrier in Development of Mucoadhesive Tablets Containing Polygoni cuspidati Extract with Potential for Dental Applications
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Isolation, Structure Determination of Sesquiterpenes from Neurolaena lobata and Their Antiproliferative, Cell Cycle Arrest-Inducing and Anti-Invasive Properties against Human Cervical Tumor Cells

1
Department of Pharmacognosy, Interdisciplinary Excellence Centre, University of Szeged, Eötvös u. 6., H-6720 Szeged, Hungary
2
Department of Organic Chemistry, University of Debrecen, H-4032 Debrecen, Hungary
3
Department of Pharmacodynamics and Biopharmacy, Faculty of Pharmacy, University of Szeged, Eötvös u. 6, H-6720 Szeged, Hungary
4
Clinical Institute of Pathology, Medical University of Vienna, Waehringer Guertel 18-20, A-1090 Vienna, Austria
5
Institute for Ethnobiology, Playa Diana, San José GT-170, Petén, Guatemala
6
Interdisciplinary Centre of Natural Products, University of Szeged, H-6720 Szeged, Hungary
*
Authors to whom correspondence should be addressed.
Pharmaceutics 2021, 13(12), 2088; https://doi.org/10.3390/pharmaceutics13122088
Submission received: 15 October 2021 / Revised: 16 November 2021 / Accepted: 26 November 2021 / Published: 5 December 2021

Abstract

:
Seven new germacranolides (13, 58), among them a heterodimer (7), and known germacranolide (4), eudesmane (9) and isodaucane (10) sesquiterpenes were isolated from the aerial parts of Neurolaena lobata. Their structures were determined by using a combination of different spectroscopic methods, including HR-ESIMS and 1D and 2D NMR techniques supported by DFT-NMR calculations. The enantiomeric purity of the new compounds was investigated by chiral HPLC analysis, while their absolute configurations were determined by TDDFT-ECD and OR calculations. Due to the conformationally flexible macrocycles and difficulties in assigning the relative configuration, 13C and 1H NMR chemical shift and ECD and OR calculations were performed on several stereoisomers of two derivatives. The isolated compounds (110) were shown to have noteworthy antiproliferative activities against three human cervical tumor cell line with different HPV status (HeLa, SiHa and C33A). Additionally, lobatolide C (6) exhibited substantial antiproliferative properties, antimigratory effect, and it induced cell cycle disturbance in SiHa cells.

1. Introduction

Sesquiterpene lactones (SLs) are one of the most prevalent and biologically significant classes of secondary metabolites of plants, comprising over 5000 known compounds. They are common in several families (e.g., Apiaceae, Solanaceae, Cactaceae and Euphorbiaceae), but the majority of them have been obtained from Asteraceae [1,2]. These lactones are classified biogenetically, according to the carbocyclic skeleton into four main groups: germacranolides, eudesmanolides, guaianolides and pseudoguaianolides. Besides these main types, there are a varieties of other lactones, formed by further modification of the carbon skeleton during biosynthesis. Germacranolides are considered as progenitors for other skeletal classes [3]. In traditionally-used medicinal plants, SLs often represents the active ingredients (e.g., achillin in Achillea millefolium, helenalin in Arnica montana and cynaropicrin in Cynara scolymus). Several SLs have been evaluated currently in cancer clinical trials, among these compounds, artemisinin, thapsigargin, parthenolide and its synthetic analogues that show low side effects and high antitumor potency, are promising anticancer leads or prodrugs [4]. Moreover, dimethylamino adduct of arglabin, an Artemisia glabella metabolite is registered as antitumor substance in the Republic of Kazakhstan [5]. In case of SLs, α-methylene-γ-lactone group is the molecular part mostly responsible for biological effects such as antitumor, anti-inflammatory and antimicrobial activities. This is a consequence of its alkylation potency acting on enzymes and transcription factors in different organisms. This ring is responsible, and often considered essential for the cytotoxic effect of SLs [6].
Neurolaena lobata (L.) R.Br. ex Cass. (Asteraceae) is a perennial flowering plant distributed in Central and South America and the West Indies. The leaves of the plant have been used by local communities primarily in case of malaria, parasitic ailments, and pain of various origins, but also for the treatment of inflammatory skin disorders, ulcers, different types of cancer and diabetes [7,8]. N. lobata is a rich source of sesquiterpene lactones. Previous phytochemical studies of the plant have led to the isolation of sesquiterpenes of germacranolide (neurolenins A-F, lobatin A and 3-epi-desacetylisovaleroylheliangine), seco-germacranolide (neurolobatins A and B), furanoheliangolide (lobatins B, and C, 8β-isovalerianyloxy-9α-hydroxycalyculatolide, 8β-isovalerianyloxy-9α-acetoxycalyculatolide and 5β-hydroxy-8β-isovaleroyloxy-9α-hydroxycalyculatolide) and eudesmanolide (3β-acetoxy-8β-isovaleroyloxyreynosin) types [9,10,11,12].
Previous pharmacological investigations of the crude leaf extract proved its anti-inflammatory [13,14], antiulcerogenic [15], antinociceptive [16], and antimicrobial [17,18,19] activities. Dichloromethane extract of the plant was tested for its antineoplastic activity in HL-60 cells and in anaplastic large cell lymphoma cell lines of human or mouse origin and antiproliferative and proapoptotic effects were detected. Furthermore, the extract induced tumor suppressors, down-regulated the expression of oncogenes, inhibited cell proliferation and triggered the apoptosis of malignant cells [20]. Moreover, neurolenin B specifically decreased pro-carcinogenic NPM/ALK expression in ALK+ ALCL cells, and attenuated tumor intra/extravasation into the lymphatics [21]. Pharmacological investigation of lobatin B revealed that this compound also inhibited the expression of NPM/ALK, and in addition, the expression of JunB and PDGF-Rβ, and attenuated proliferation of ALCL cells by arresting them in late M phase [22]. In vitro assays for cytotoxic activity of sesquiterpene lactones of N. lobata revealed many compounds to be effective against GLC4 human small cell lung carcinoma and COLO 320 human colorectal cancer cells. Lobatin B possessed the highest cytotoxicity (IC50 = 0.6 and 1.1 µM, in case of GLC4 and COLO 320, respectively), followed by neurolenin B (IC50 = 1.1 and 1.2 µM) [23].
The aim of the present work was the isolation of further biologically active sesquiterpenes from N. lobata. This paper reports the isolation and structure elucidation of seven new SLs (13, 58), including a dimeric lactone (7), and three known sesquiterpenes (4, 9, and 10) from the CH2Cl2 extract of the frozen herb. The isolated compounds were evaluated for their antiproliferative activities against three human cervical cancer cell lines with different human papillomavirus (HPV) status including HPV-negative C33A, as well as HPV-16 and 18-positive SiHa, and HeLa, respectively. The most promising natural product (6) was additionally investigated in order to reveal its cancer selectivity, action on cancer cell migration and cell cycle distribution.

2. Materials and Methods

2.1. General Experimental Procedures

The high-resolution MS spectra were acquired on a Thermo Scientific Q-Exactive Plus orbitrap mass spectrometer equipped with ESI ion source in positive ionization mode. The samples were dissolved in MeOH. The data were acquired and processed with the MassLynx software. Data acquisition and analysis were accomplished with Xcalibur software version 2.0 (Thermo Fisher Scientific, Waltham, MA, USA). A Bruker Avance DRX 500 spectrometer [500 MHz (1H) and 125 MHz (13C)] was used for recording the NMR spectra. The signals of the deuterated solvent CDCl3 were taken as reference. 2D NMR data were acquired and processed with standard Bruker software. Gradient-enhanced versions were used in the 1H–1H COSY, HSQC and HMBC experiments. Optical rotations were determined in CHCl3 by using a Perkin-Elmer 341 polarimeter. ECD spectra were recorded on a JASCO J-810 spectropolarimeter.
Rotational planar chromatography (RPC) was performed by using a Chromatotron apparatus (Harrison Research, Palo Alto, CA, USA) on self-coated silica plates (Kieselgel 60 GF254, 15 μm, Merck, Darmstadt, Germany). For column chromatography (CC), polyamide (MP Polyamide, 50–160 μm, MP Biomedicals, Irvine, CA, USA) and silica gel (Kieselgel 60, 63–200 μm, Merck, Darmstadt, Germany) were used. Thin-layer chromatography (TLC) was carried out using silica gel (Kieselgel 60 F254, Merck) and RP-C18 (F254s, Merck) pre-coated plates. The preparative TLC (prep-TLC) was performed on pre-coated silica gel plates (20 × 20 cm2, Kieselgel 60 F254, Merck). The compounds were detected by spraying with cc. vanillin sulfuric acid, followed by heating (120 °C).

2.2. Plant Material

Neurolaena lobata (L.) R.Br. ex Cass. (Asteraceae) was collected by R. Diaz, and R. O. Frisch (Institute for Ethnobiology, Playa Diana, GT-170 San José/Petén, Guatemala), in the flowering period, in the area of the Chakmamantok-rock formation (16 59′16″ N, 89 53′45″ W) in San José, Guatemala. Voucher specimens were archived at the herbarium of the Institute for Ethnobiology, San Jose, Guatemala, and at the Herbarium of the Department of Pharmacognosy, University of Szeged, Szeged, Hungary (No. 813). The fresh plant material (aerial parts) was air-dried and stored deep-frozen until subsequent extraction.

2.3. Extraction and Isolation

The dried and ground aerial parts of the plant (3.00 kg) was percolated with MeOH (50 L) at room temperature. The extract was concentrated under reduced pressure and solvent–solvent partition was performed with 5 × 1 L petroleum ether (A), 5 × 1 L of CH2Cl2 (B), and finally with 5 × 1 L of EtOAc (C). The CH2Cl2 phase (95.4 g) was separated on a polyamide column (287 g) with mixtures of MeOH and H2O (1:4, 2:3, 3:2 and 4:1, 3 L of each) as eluents to afford seven fractions (BI–BVII). Fraction BII (33.6 g), obtained with MeOH–H2O (1:4), was subjected to silica gel VLC, using a gradient system of cyclohexane–EtOAc–EtOH (from 3:1:0 to 5:5:1). The fractions were combined after TLC monitoring into nine fractions (BII/1-BII/9). Fraction BII/4 (98.7 mg) was separated by RPC with a CH2Cl2–acetone gradient system (from 99:1 to 9:1, and finally MeOH). The subfraction (35.6 mg) eluted with CH2Cl2–acetone 97:3 was subjected to preparative TLC, using a system of CH2Cl2–acetone (95:5), to give compound 9 (4.6 mg). Fracion BII/7 was separated by VLC with the gradient system of CH2Cl2–acetone (from 99:1 to 4:1, and finally MeOH) to yield 12 subfraction. Subfraction 10 (767 mg) was further chromatographed by VLC with cyclohexane–EtOAc–EtOH gradient elution (from 60:20:0.5 to 6:6:1) and fraction eluted with cyclohexane–EtOAc–EtOH 60:30:0.5 (338 mg) was separated by VLC with CH2Cl2–MeOH (from 100:1 to 100:3.5). Fractions eluted with CH2Cl2–MeOH 100:1.5 was combined and finally purified by RP TLC with MeOH–H2O 7:3 to yield compound 10 (2.1 mg) Fraction BIII (7.6 g) obtained from polyamide column with MeOH–H2O (3:2), was subjected to silica gel VLC, using a gradient system of cyclohexane–EtOAc–EtOH (from 30:5:0 to 30:30:2) to yield 10 main fraction (BIII/1–10). Fraction BIII/5 (230 mg) was further purified by RP VLC with MeOH–H2O gradient system (from 2:3 to 9:1). Fractions eluted with MeOH–H2O (7:3) were combined (40 mg) and finally preparative TLC (CHCl3–acetone 97:3) was applied to isolate compounds 4 (1.8 mg), 3 (3.4 mg) and 1 (1.98 mg). Fraction BIII/6 (2.7 g) was firstly separated by VLC with gradient system of CHCl3–acetone (from 99:1 to 9:1). Fractions were combined into 8 subfractions according to their components. Subfraction 1 (32 mg) was further purified by preparative TLC with the mobile phase toluene–EtOAc 3:2, affording the isolation of compounds 5 (1.6 mg), 6 (17.4 mg) and 2 (1.5 mg). Subfraction 4 was purified first by VLC with a gradient system of cyclohexane–EtOAc–EtOH (from 3:1:0 to 30:30:3). Fractions obtained with cyclohexane–EtOAc 3:2 was combined and further purified by preparative TLC with cyclohexane–EtOAc–EtOH (30:30:1) and compound 8 (1.5 mg) was isolated. Fraction BIII/8 (0.6 g) was chromatographed on VLC with a gradient system of cyclohexane–EtOAc–EtOH (from 3:1:0 to 30:30:3) and 10 main fractions were yielded. Fraction BIII/9 (0.7 g) was separated by VLC with a gradient system of CHCl3–MeOH (from 99:1 to 9:1), and 9 main fractions were yielded. Fraction 5 was purified by preparative TLC with cyclohexane–EtOAc–EtOH (30:30:1) as eluent and compound 7 (2.1 mg) was isolated.

2.4. Physical Characteristics of New Compounds

8β-Hydroxy-9α-isovaleroyloxycalyculatolide (1): an amorphous solid; [α]28D +12 (c 0.1, CHCl3). ECD {CH3CN, λmax [nm] (Δε)}, c 2.84 × 10 − 4 M: 314 (+0.36), 256 (+1.72), 215 (−2.76); 1H NMR and 13C NMR data, see Table 1 and Table 2; HR-ESIMS m/z 379.1755 [M + H]+ (calcd for C20H27O7, 379.1757).
Lobatolide A (2): an amorphous solid; [α]27D +47 (c 0.05, CHCl3); 1H NMR and 13C NMR data, see Table 1 and Table 2; HR-ESIMS m/z 411.2019 [M − H2O + H]+, (calcd for C21H31O8, 411.2024).
Lobatin E (3): an amorphous solid; [α]27D −25 (c 0.2, CHCl3); ECD {CH3CN, λmax [nm] (Δε)}, c 7.37 × 10 − 4 M: 330sh (+0.84), 319 (+1.11), 288 (−0.08), 279 (+0.13), 256 (−1.06), 233 (+1.61), 207 (−6.10); 1H NMR and 13C NMR data, see Table 1 and Table 2; HR-ESIMS m/z 377.1597 [M + H]+ (calcd for C20H25O7, 377.1600).
Lobatin C (4): an amorphous solid; [α]27D −4 (c 0.1, CHCl3). ECD {CH3CN, λmax [nm] (Δε)}, c 2.75 × 10 − 4 M: 330sh (+0.80), 321 (+1.05), 288 (−0.14), 277 (+0.03), 256 (−0.49), 227 (+1.49), 204 (−4.48); 1H NMR and 13C NMR data, see Table 1 and Table 2; APCI-MS positive m/z 419 [M + H]+, 436 [M+NH4]+.
Lobatolide B (5): an amorphous solid; [α]28D −69 (c 0.2, CHCl3); ECD {CH3CN, λmax [nm] (Δε)}, c 4.51 × 10 − 4 M: 302 (−0.75), 231 (−9.49), 199 (+2.19); 1H NMR and 13C NMR data, see Table 1 and Table 2; HR-ESIMS m/z 455.2274 [M + H]+ (calcd for C23H35O9 455.2281).
Lobatolide C (6): a colorless gum; [α]28D +34 (c 0.2, CHCl3); ECD {CH3CN, λmax [nm] (Δε)}, c 4.80 × 10 − 4 M: 259 (−1.15), 221sh (+2.91), 196 (+21.51); 1H NMR and 13C NMR data, see Table 1 and Table 2; HR-ESIMS m/z 407.2063 [M + H]+ (calcd for C22H31O7, 407.2070).
Lobatolide D (7): white amorphous powder; [α]27D −30 (c 0.1, CHCl3); ECD {CH3CN, λmax [nm] (Δε)}, c 1.38 × 10 − 4 M: 328sh (+1.17), 315 (+1.99), 255sh (+5.76), 246 (+7.25), 217 (−27.65); 1H NMR and 13C NMR data, see Table 3; HR-ESIMS positive m/z 799.3551 [M + H]+ (calcd for C42H55O15, 799.3541); 821.3352 [M+Na]+ (calcd for C42H54O15Na, 821.3360).
8β-Isovaleroyloxyreynosin (8): a colorless oil; [α]28D +8 (c 0.1, CHCl3); ECD {CH3OH, λmax [nm] (Δε)}, c 4.59 × 10 − 4 M: 248 (−0.32), 222 (+0.17); 1H NMR and 13C NMR data, see Table 3; HR-ESIMS m/z 349.2017 [M + H]+, (calcd for C20H29O5, 349.2015).
Volenol (9): ECD {CH3CN, λmax [nm] (Δε)}, c 3.52 × 10 − 4 M: 194 (−3.50).

2.5. Computational Section

Mixed torsional/low-mode conformational searches were carried out by means of the Macromodel 10.8.011 software using the MMFF with an implicit solvent model for CHCl3 applying a 21 kJ/mol energy window [24]. Geometry re-optimizations of the resultant conformers [CAM-B3LYP/TZVP PCM/MeCN, CAM-B3LYP/TZVP PCM/CHCl3, B3LYP/6-31+G(d,p) level in vacuo], TDDFT-ECD [B3LYP/TZVP PCM/MeCN, BH&HLYP/TZVP PCM/MeCN, CAM-B3LYP/TZVP PCM/MeCN and PBE0/TZVP PCM/MeCN], or [B3LYP/TZVP PCM/CHCl3, BH&HLYP/TZVP PCM/CHCl3, CAM-B3LYP/TZVP PCM/CHCl3 and PBE0/TZVP PCM/CHCl3] and DFT-NMR calculations [mPW1PW91/6-311+G(2d,p) and mPW1PW91/6-311+G(2d,p) SMD/CHCl3] were performed with the Gaussian 09 package [25]. ECD spectra were generated as sums of Gaussians with 2400–3600 cm−1 width at half-height, using dipole-velocity-computed rotational strength values [26]. Computed NMR shift data were corrected with I = 185.4855 and S = −1.0306 for the carbons and I = 31.8996 and S = −1.0734 for the hydrogens in the gas-phase and I = 186.5242 and S = −1.0533 for the carbons and I = 31.8018 and S = −1.0936 for the hydrogens in the SMD calculations [27,28]. Boltzmann distributions were estimated from the CAM-B3LYP and the B3LYP energies. Visualization of the results was performed by the MOLEKEL software package [29].

2.6. Antiproliferative MTT Assay

Antiproliferative effect of the isolated compounds were determined in vitro using SiHa (HPV 16+), HeLa (HPV 18+), C33A (HPV negative) human cervical cell lines, NIH-3T3 mouse embryonic and MRC-5 human fibroblast cells by means of MTT ([3-(4,5-dimethylthiazol-2-yl)-2,5-diphenyltetrazolium bromide]) assay. Briefly, a limited number of human cancer cells (5000/well for SiHa and HeLa cells, 10,000/well in case of C33A cells) were seeded onto a 96-well microplate and became attached to the bottom of the well overnight. On the second day of the procedure, the test substances were added in two concentrations (10.0 and 30.0 µM) in order to obtain preliminary data and then the compounds were applied in serial dilution (applied final concentrations were: 0.1, 0.3, 1.0, 3.0, 10.0, 30.0 µM). Further details of the experiment were described in ref. [30]. All in vitro experiments were carried out on two 96-well dishes with at least five parallel wells [31].

2.7. Cell Cycle Analysis by Flow Cytometry

Cell cycle analysis was performed as described in [32,33].

2.8. Wound Healing Assay

In order to assess antimetastatic activity of the tested compound, wound healing assay was performed. The assay was done with specific wound healing assay chambers (Ibidi GmbH, Gräfelfing, Germany). SiHa cells were collected and 35,000 cells were seeded into both chambers of the insert. Cells were let to be attached to the plate surface during an overnight incubation at 37 °C in 5% CO2 atmosphere and then the inserts were removed. Cell debris was removed by a washing step with PBS. Test compounds were added to the wells in increasing concentrations in 2% FBS containing medium for 24 and 48 h. Migration of the cells into the wound site was visualized by phase-contrast inverted microscope (Axiovert 40, Zeiss, Thornwood, NY, USA). Images were taken with CCD camera at definite intervals and the migration of the cells was calculated as the ratio of wound closure by ImageJ software [34].

2.9. Statistical Analysis

Statistical analysis of the obtained data was performed by analysis of variance (ANOVA) followed by Dunnett’s test. All analyses were performed with GraphPad Prism 5 (GraphPad Software; San Diego, CA, USA).

3. Results and Discussion

Phytochemical investigation of the CH2Cl2-soluble phase of the MeOH extract, obtained from the aerial parts of N. lobata, resulted in the isolation of ten sesquiterpenes (110) (Figure 1). For the separation of the compounds a combination of different chromatographic methods, including open column chromatography (OCC), vacuum liquid chromatography (VLC), rotation planar chromatography (RPC), and preparative TLC were used. The structure elucidation was carried out by extensive spectroscopic analysis, using HR-ESIMS, and 1D and 2D NMR (1H–1H COSY, HSQC, HMBC and NOESY) spectroscopy and DFT-NMR calculations.

3.1. Structure Elucidation of the Isolated Compounds

3.1.1. Lobatolide A (1)

Compound 1 was isolated as an amorphous solid with [α]28D +12 (c 0.1, CHCl3). Its HR-ESIMS displayed a quasimolecular ion peak at m/z 379.1755 [M + H]+ (calcd for 379.1757), indicating the molecular formula, C20H26O7. The 1D and 2D NMR spectra exhibited resonances for a carbonyl group (δC 203.9), an α-methylene-γ-lactone ring [δH 4.78 dd (H-6), 3.30 dd (H-7), 6.30 d (H-13a) and 5.35 d (H-13b), δC 75.3 (C-6), 47.9 (C-7), 141.3 (C-11), 169.2 (C-12) and 122.6 (C-13)], a quaternary carbon (δC 88.9), a trisubstituted olefin (δH 5.60 d, δC 104.0 and 192.8), three methines (δH 3.01 m, 4.13 t and 5.27 d, δC 31.5, 73.8 and 77.1), and one secondary (δH 1.40 d, δC 16.1) and one tertiary methyl group (δH 1.37 s, δC 18.9) (Table 1 and Table 2). Further 1H and 13C NMR signals at δH 2.38 dd, 2.33 m, 2.19 m, 2 × 1.02 d, and δC 171.5, 43.1, 25.7, 2 × 22.4, were attributed to an isovaleroyloxy group. The chemical shifts and coupling constants of 1 were closely related to those of the 1-keto-furanoheliangolide derivative 8-isovaleroyloxy-9α-hydroxy-calyculatolide [9], with the only difference changing of the positions of isovaleroyloxy and hydroxy groups. The isovaleroyl group was assigned to C-9 with regard to the HMBC correlations between H-9 and the carbonyl carbon signal of isovaleroyl group. The relative configuration of compound 1 was determined by means of a NOESY experiment. The coupling constant of H-6 and H-7 (J6,7 = 5.0 Hz) indicated the β orientation of H-6 and α orientation of H-7, found in all sesquiterpenes isolated previously from this plant. NOESY correlation between H-6 and H-15 confirmed the β position of the 15-methyl group, while NOE effects observed between H-7 and H-5a, H-5a and H-4, H-7 and H-13b, and H-13b and H-8 dictated the α orientation of these protons. All of the above evidence supported the structure of this compound as 8β-hydroxy-9α-isovaleroyloxycalyculatolide, and named as lobatolide A (1).
The absolute configuration of 1 was determined by the solution TDDFT-ECD method [35,36]. Merck Molecular Force Field (MMFF) conformational search of (4R,6R,7R,8S,9R,10R)-1 resulted in 133 initial conformers in a 21 kJ/mol energy window. These structures were re-optimized at the CAM-B3LYP/TZVP [37] PCM/MeCN levels yielding 12 low-energy conformers over 1% Boltzmann population. ECD spectra computed at various levels of theory for these conformers reproduced well the experimental ECD spectrum (Figure 2). Furthermore, the low-energy conformers, differing only in the orientation of the OiVal and the OH groups (Figure 3), gave similar computed ECD spectra allowing the solid elucidation of the absolute configuration as (4R,6R,7R,8S,9R,10R)-1. To test the applicability of the DFT-NMR method [28,38], the 133 MMFF conformers were also re-optimized at the B3LYP/6-31+G(d,p) level yielding seven low-energy conformers. The NMR shift values were computed for these conformers both at the mPW1PW91/6-311+G(2d,p) and the mPW1PW91/6-311+G(2d,p) SMD/CDCl3 levels [39]. Computed and experimental 13C NMR chemical shift values gave a good overall agreement with CMAE (corrected mean absolute error) values [40] of 1.92 and 1.77 ppm, respectively (Table S1). Although the 13C NMR DFT calculation showed good agreement for compound 1, the exomethylene group caused larger deviations in the computed NMR shift values [41,42] of the other derivatives of this work even by testing at different levels of theory.

3.1.2. Lobatolide B (2)

Compound 2 (lobatolide B) was isolated as amorphous solid with [α]27D +47 (c 0.05, CHCl3). The molecular ion peak at m/z 411.2019 [M + H]+ (calcd for 411.2024) in the positive-ion HR-ESIMS showed its molecular formula to be C21H30O8. The 1H NMR spectrum displayed the presence of signals of an isovaleroyl side-chain (δH 2.14 dd, 2.10 brd, 2.02 m, 0.93 d, 0.91 d), and a methoxy signal (δH 3.40 s) (Table 1). From the 1H–1H COSY spectrum, one structural fragment with correlated protons was assigned: −CH(CH3)-CH2-CHR-CHR-CHR-CHR– (δH 3.01, 1.34, 2.63, 2.03, 4.29, 3.16, 5.07 and 3.94) (C-4–C-9). Long-range heteronuclear NMR correlations between the quaternary carbon C-3 and H-2, H-4, H-5b, and H3-15, between C-10 and H-2, H-8, and H3-14 signals, and between C-1 and H-2 and H-14 were observed in the HMBC spectrum. The two- and three-bond correlations revealed that the structural fragment together with C-1, C-3, C-10, the C-2 methine, and the 14-methyl group forms a furanoheliangolide skeleton. Detailed analysis of the spectral data suggested that compound 2 is very similar to 8β-isovalerianoyloxy-9α-hydroxycalyculatolide. The difference between the two compounds was the modified α-methylene-γ-lactone moiety. The upfield-shifted chemical shifts of H2-13 protons (δH 3.74 and 3.67) indicated an oxymethylene part in the structure instead of the exomethylene unit. Moreover, the presence of sp3 methine (δH-11 2.89, δC-11 49.3) resonances were observed. These observations led to the conclusion, that the C-11/C-13 double bond was saturated. 1H–1H COSY correlations between H-7 and H-11, H-11-and H-13, and long-range correlation between OMe and H2-13 proved that a methoxy group is connected at C-13.
The relative configuration of 2 corresponded well with the data reported for furanoheliangolides in the literature, containing an exomethylene moiety at the C-11 position [12]. The NOE correlations between H-8 and H-11 suggested β orientation of the methoxymethylene group at C-11. NOE correlations in a conformationally flexible macrolide ring have to be treated with precaution, since the large flexibility of the ring can lead to wrong assignment as reported in the literature for other related derivatives [43,44]. Thus the relative configuration of 2 was investigated by DFT-13C NMR calculations of eight selected stereoisomers with different configuration at the C-8, C-9 and C-11 chirality centers. The (4R,6R,7S,8S,9R,10R,11S), (4R,6R,7S,8S,9R,10R,11R), (4R,6R,7S,8R,9R,10R,11S), (4R,6R,7S,8R,9R,10R,11R), (4R,6R,7S,8R,9S,10R,11S), (4R,6R,7S,8R,9S,10R,11R), (4R,6R,7S,8S,9S,10R,11S) and (4R,6R,7S,8S,9S,10R,11R) stereoisomers of 2 were investigated at the mPW1PW91/6-311+G(2d,p)//B3LYP/6-31+G(d,p) and mPW1PW91/6-311+G(2d,p) SMD/CDCl3//B3LYP/6-31+G(d,p) levels, which were already tested for compound 1. The NMR calculations showed that the H-8/H-11 NOE correlation is also feasible with trans relative configuration of H-8 and H-11 due to the flexibility of the macrolide ring. Moreover, this NOE effect is much more probable with the β orientation of H-11 hydrogen than with the α one. The H-8–H-11 interatomic distance in the lowest-energy conformer of (4R,6R,7S,8S,9R,10R,11S)-2 (with α orientation of H-11) is 3.50, while this value is 2.19 Å for the (4R,6R,7S,8S,9R,10R,11R)-2 stereoisomer (with β orientation of H-11) (Figure 4). Thus, the 13C NMR calculations revealed that the H-8/H-11 NOE correlation unexpectedly derives from the trans relative configuration of the H-8 and H-11 protons.
The computed 13C NMR data showed clear preference for the (4R,6R,7S,8S,9R,10R,11R)-2 stereoisomer with a CMAE value of 1.61 vs. 2.36, 2.66, 2.51, 2.11, 2.22, 2.44 and 2.25 in the gas-phase calculations and 1.58 vs. 2.36, 2.71, 2.32, 2.27, 2.30, 2.37 and 2.23 in the SMD calculations (Tables S2 and S3). The DP4+ statistical analysis gave a 99.98% confidence in the gas-phase and 100.00% with SMD solvent model for this isomer [45,46]. The data above corroborate the proposed molecular formula of lobatolide B (2). The ECD spectrum of 2 was recorded both in MeCN and MeOH with different settings but nearly baseline ECD spectra were recorded most likely due to solubility problems. Since +47 value was measured for the specific rotation, this was used instead of the ECD to determine the absolute configuration [47]. CAM-B3LYP/TZVP PCM/CHCl3 re-optimization of the 136 MMFF conformers resulted in 17 low-energy conformers. OR values were computed for all the conformers at four different levels of theory and they had positive value in all the combinations and the Boltzmann weights were in the range of +53–+71 showing good agreement with the +47 experimental value (Table S4). Thus the absolute configuration of 2 was assigned as (4R,6R,7S,8S,9R,10R,11R) in line with the biosynthetic considerations.

3.1.3. Lobatolide C (3)

Compound 3 was isolated as an amorphous solid with [α]27D −25 (c 0.2, CHCl3). HR-ESIMS and detailed 1D and 2D NMR spectroscopic studies led to a furanoheliangolide structure for compound 3, found previously in case of lobatin B [10], on the basis of the following arguments. HR-ESIMS suggested that the molecular composition of the compound is C20H24O7, according to the quasimolecular ion peak at 377.1597 [M + H]+ (calcd for 377.1600). Analysis of the 1H and 13C NMR spectra confirmed that lobatin B and compound 3 are structurally related. The 1H and 13C NMR spectra revealed the presence of an isovaleroyl group (δH 2.37 dd, 2.31 dd, 2.16 m, 2 × 1.01 d; δC 171.7, 43.1, 25.6 and 2 × 22.4) and two quaternary carbon-bonded methyl groups (δH 1.41 s, 2.06 s, δC.18.4, 19.6) in the molecule (Table 1 and Table 2). Informative signals at δH 5.71 m (H-6), 3.57 dd (H-7), 6.29 d (H-13a), and 5.35 d (H-13b) and δC 74.9 (C-6), 45.2 (C-7), 141.2 (C-11), 170.0 (C-12), and 123.2 (C-13) additionally indicated the presence of a trans-fused α-methylene-γ-lactone ring with H-6 in the β position. The isovaleroyl group was assigned to C-9 with regard to the HMBC correlations between H-9 and the carbonyl carbon signal (isovaleroyl CO). The only difference between compound 2 and lobatin B is the connection of isovaleroyl (at C-9 in 3, and at C-8 in lobatin B), and hydroxy (at C-8 in 3, and at C-9 in lobatin B) groups.
The relative configuration of compound 3 was studied by means of a NOESY experiment. Cross-peak between H-5 and H3-15, and the small coupling constant (J5 = 1.8 Hz) detected at H-5 proved the Z-configuration of the double bond (C-4–C-5). NOE effects observed between H-7 and H-13b, and H-13b and H-8 dictated the α orientation of these protons. All of the above evidence supported the structure of 3 as 8β-hydroxy-9α-isovaleroyloxyatripliciolide, which was named as lobatolide C.
CAM-B3LYP PCM/MeCN re-optimization of the initial 92 MMFF conformers yielded 13 low-energy conformers, which can be sorted into two groups based on the axial or equatorial arrangement of the C-9 OiVal substituent (Figure 5). The conformers of a group showed difference in the rotation of the C-9 OiVal group (Figure 6). The conformers of groups A and B had markedly different computed ECD spectra, while conformers in the same group had similar curves. Although group B conformers had low populations, they had more intense computed ECD spectra than those of group A, and they governed the feature of the sum ECD curve above 240 nm (Figure 7). The Boltzmann-weighted ECD spectra reproduced well the major transitions of the experimental ECD spectrum allowing elucidation of the absolute configuration as (6R,7R,8S,9R,10R), which is in accordance with the biosynthetic considerations.

3.1.4. Lobatin C (4)

Compound 4 was isolated as an amorphous solid with [α]27D −4 (c 0.1, CHCl3). Based on 1D and 2D NMR data it was determined as lobatin C, isolated previously from N. lobata and N. macrocephala by Passreiter et al. [9,48]. As the compound was identified only by GC-MS analysis [9], and no NMR data were published previously in the literature, detailed NMR studies were performed affording the 1H and 13C NMR assignments (Table 1 and Table 2). CAM-B3LYP/TZVP PCM/MeCN re-optimization of the initial 109 MMFF conformers of (6R,7S,8S,9R,10R)-4 yielded 10 low-energy conformers and similarly to 3, two conformer groups could be identified with different conformations of the macrolide ring. In group A representing 90.2% sum population, the C-8 OiVal substituent had axial orientation, while it was equatorial in group B conformers (5.6% sum population). The Boltzmann-weighted CAM-B3LYP/TZVP PCM/MeCN ECD spectra reproduced well all the major transitions of the experimental ECD spectrum and hence the absolute configuration was determined as (6R,7S,8S,9R,10R) (Figure 8, Figures S21 and S22).

3.1.5. Lobatolide D (5)

Compound 5 was isolated as an amorphous solid with [α]28D −69 (c 0.2, CHCl3). The molecular formula C23H34O9 was determined by HR-ESIMS with a molecular ion peak at m/z 455.2274 [M + H]+ (calcd for C23H35O9 455.2281). Detailed analysis of spectral data afforded the structure, which was in good agreement with neurolenin B. The missing exomethylene signals (δH 6.31 and 5.81, br s each), as well as additional sp3 methine at δH 3.04 m, δC 40.2 (C-11) and oxymethylene resonances detected at δH 3.34 t and 3.68 dd and δC 66.3 (C-13), indicated the saturation of the C-11/C-13 double bond (Table 1). The H-7/H-11/H-13 spin system was established by means of 1H–1H COSY correlations and HMBC cross peaks of H-7/C-11, H2-13/C-7 and H2-13/C-11. The 1H-spectrum contained an additional methoxy signal (δH 3.39 s), which exhibited HMBC correlation with δC 66.3 and, therefore, must be situated at C-13. The relative configuration of compound 5, named as lobatolide D, was found to be identical with that of neurolenin B. The presence of trans-fused lactone-ring was demonstrated by NOE-effects of H-5a/H-6, H-5b/H-7, H-5b/H-4. The α position of H-11 was suggested by the NOE correlations H-11/H-7 and H-11/H-4, which implies a different relative configuration of the C-8, C-9, C-11 centers from that of compound 2. In order to assist the assignment of the relative configuration, DFT-13C NMR calculations were performed on eight stereoisomers differing in the C-8, C-9, C-11 chirality centers. The (4R,6R,7S,8S,9R,10R,11S), (4R,6R,7S,8S,9R,10R,11R), (4R,6R,7S,8R,9R,10R,11S), (4R,6R,7S,8R,9R,10R,11R), (4R,6R,7S,8R,9S,10R,11S), (4R,6R,7S,8R,9S,10R,11R), (4R,6R,7S,8S,9S,10R,11S) and (4R,6R,7S,8S,9S,10R,11R) stereoisomers of 5 were computed at the mPW1PW91/6-311+G(2d,p)//B3LYP/6-31+G(d,p) and mPW1PW91/6-311+G(2d,p) SMD/CDCl3//B3LYP/6-31+G(d,p) levels. The lowest-energy conformers of (4R,6R,7S,8S,9R,10R,11S)-5 and (4R,6R,7S,8S,9R,10R,11R)-5 differing in the configuration of the C-11 center showed that the H-11/H-7 NOE correlation is feasible for both of them with 2.33 Å interatomic distance for the (11S) and 2.86 Å for the (11R) epimer. The H-11/H-4 distance is found shorter (3.47 Å) for the (11S) epimer [4.96 Å in the (11R) epimer]. The computed 13C NMR chemical shift data did not show a clear preference for any stereoisomers, and thus the 1H NMR data were taken also into account (Tables S5–S8). The DP4+ statistical analysis of the combined 13C and 1H NMR data gave 99.37% confidence with the SMD solvent model for the (4R,6R,7S,8S,9R,10R,11S)-5 stereoisomer. Interestingly, the (4R,6R,7S,8S,9R,10R,11R) stereoisomer of 5, which is homochiral with 2, gave 0.00% probability at both NMR levels. This result suggested that C-11 chirality center of 2 and 5 was produced by the reduction of the exomethylene moiety, which afforded different C-11 configurations in compounds 2 and 5.
In the knowledge of the relative configuration, the absolute configuration was determined by the TDDFT-ECD protocol performed on the (4R,6R,7S,8S,9R,10R,11S)-5 stereoisomer (Figure 9). CAM-B3LYP re-optimization of the initial 89 MMFF conformers resulted in 12 low-energy conformers, which gave quite similar ECD spectra (Figure 10). The Boltzmann-weighted ECD spectra gave a good agreement with the major transitions of the experimental ECD spectrum at all the applied levels of theory allowing elucidation of the absolute configuration as (4R,6R,7S,8S,9R,10R,11S)-5.

3.1.6. Lobatolide E (6)

Compound 6 was obtained as a colorless gum with [α]28D +34 (c 0.2, CHCl3). From the molecular ion peak at m/z 407.2063 [M + H]+ (calcd for 407.2070) in the positive-mode HR-ESIMS, its molecular formula was determined to be C22H30O7. The 1H NMR spectrum showed the presence of signals due to an isovaleroyl and an acetyl group (Table 1). Further, the 13C-JMOD spectrum suggested that the skeleton consists of fifteen carbons (Table 2). From the 1H–1H COSY spectrum two structural fragments were assigned on the basis of correlated protons: −CHR-CH2-CH– (A) (δH 2.93, 2.35, 1.64 and 5.50) (C-1–C-3) and =CH-CHR-CHR-CHR-CH2– (B) (δH 5.66, 5.22, 3.21, 5.74, 2.63 and 1.40) (C-5–C-9). These structural elements, tertiary methyls and quaternary carbons were connected by examination of the long-range C–H correlations detected in the HMBC spectrum. The two- and three-bond correlations between H-3, H-6, H-15, H2-2 and the quaternary carbon C-4, and between H2-2, H-8, H2-9, H-14 and the quaternary C-10 and signals revealed that structural fragment A together with C-10, C-4, and the 14- and 15-methyl groups forms a germacrane skeleton. The lactone ring connected to the macrocycle in position C-6, C-7 was evident from the HMBC cross-peaks between H2-13, and C-12, and H2-13, H-7 and C-11. The position of the ester groups was proved by the long-range correlation between the ester carbonyl of isovaleroyl group (δC 173.2) and H-8 (δH 5.74 d) and the acetyl carbonyl (δC 171.8) and H-3 (δH 5.50 d). The remaining epoxy group, which were elucidated from the molecular composition, was placed at C-10–C-1, with regard to the 13C NMR chemical shifts (δC-1 65.6 and δC-10 61.4) and literature data for similar epoxy germacranolides [49]. The relative configuration of the chiral centers was studied by NOESY measurements. Diagnostic NOE correlations were detected between H-6 and H3-14, demonstrating the β orientation of these protons. Furthermore, NOESY cross-peaks were observed between H-1/H-2a, H-1/H-3, H-2a/H-3, H-1/H-9b, H-9b/H-7, H-7/H-8, H-14/H-9a, H-14/H-2b, H-6/H-14, indicating the α-oriented H-1, H-2a, H-7, H-8 and H-9b and the β-orientation of H-6, H-9a, H-2b and H-14. All of the above evidence was used to propose the structure of this compound as depicted in structural formula 6. Compound 6 is the 8-epimer of 8α-isovaleryloxy-8-desacylviguestenin, reported earlier from Viguiera procumbens (Asteraceae) [49] Moreover, in contrast to 8α-isovaleryloxy-8-desacylviguestenin, the trans position of double bond between C-4–C-5 was proved by the NOE correlations between H-6 and H3-15 (above the plane of the macrocycle) and between H-3 and H-5 (below the plane of the ring). Therefore, lobatolide E was elucidated as depicted structural formula 6.
ECD calculations performed for the 10 low-energy CAM-B3LYP conformers of (1R,3S,6R,7R,8R,10R)-6 obtained from the re-optimization of the initial 74 MMFF conformers gave nice overall agreement with the experimental ECD spectrum (Figure 11). In all the 10 computed conformers, the C-8 OiVal substituent adopted axial orientation and the conformers gave very similar ECD spectra allowing elucidation of the absolute configuration of 6 as (1R,3S,6R,7R,8R,10R) (Figure 12).

3.1.7. Lobatolide F (7)

Compound 7 was obtained as a white amorphous powder with [α]27D −30 (c 0.1, CHCl3). The molecular formula was established as C42H54O15 from its HR-ESIMS, which showed the quasimolecular ion peak at [M + H]+ 799.3551 (calcd for C42H55O15, 799.3541), requiring 16 degrees of unsaturation. The JMOD experiment displayed 42 resonances, including seven carbonyl, eight olefinic, nine methyl, four sp3 methylene, twelve sp3 methine (six of which oxymethine) and two oxygen-bearing sp3 quaternary carbons (Table 3). The combination of 1H and 1H–1H COSY spectra, revealed the presence of one cis-double-bond [H-2 (δH 6.22 d, J = 11.6 Hz), H-3 (δH 6.04 t, J = 11.6 Hz)], an isolated [H-2′ (δH 5.64, s)] and two exocyclic olefinic protons [H-13a’ (δH 6.29 d, J = 3.0 Hz); H-13b’ (δH 5.74 d, J = 3.0 Hz)], which are characteristic parts of the α-methylene-γ-lactone moiety (Table 3). Three ester residues were also identified as an acetyl [(δH 2.15 s (3H)] and two isovaleroyl-groups [H-2a′ (δH 2.27 dd, J = 15.8, 6.8 Hz); H-2b′ (δH 2.15, m); H-3′ (δH 2.02 m); H-4′ (δH 0.94 d, (3H) J = 6.7 Hz); H-5′ (δH 0.93 d, (3H), J = 6.7 Hz)], and [H-2″ (δH 2.08 m, (2H)]; H-3″ (δH 1.96 m); H-4″ (δH 0.90 d, (3H), J = 6.7 Hz); H-5″ (δH 0.87 d, (3H), J = 6.7 Hz)]. Comprehensive analysis of the 1D and 2D spectra led to the assignments of two different germacranolide units, designated as A and B.
Comparison of the NMR data of part A with those of sesquiterpenes previously isolated from the plant showed great similarities of unit A with neurolenin B. However, the absence of exomethylene and H-7 signals, as well as the appearance of proton resonances at δH 3.05 (dd, J = 14.5, 8.8 Hz) and δH 2.48 (dd, J = 14.5, 6.8) attributed to a saturated methylene suggested the relocation of the double bond from C-11–C-13 to C-7–C-11. This conclusion was supported by the significantly downfield-shifted signal of H-8 (δH 6.38) (neurolenin B: δH 5.56), in addition with HMBC correlations observed between the quaternary olefinic C-7 (δC 154.2) and H-6, H-8 and H-13a. Careful interpretation of the spectral data allowed determining part B as a 8β-isovaleroyloxy-9α-hydroxycalyculatolide unit. Chemical shifts of H-5′ (δH 3.54 m, 1H) and C-5 (δC 45.6) clearly indicated that part B is connected to A in position C-5. In addition, with the 1H–1H COSY correlations between H-13a/H-5′ and H-13b/H-5′ unambiguously established the linkage of units A and B through C-13/C-5′. HMBC cross-peaks of H-13a/C-4′, H-13a/C-5′ and H-13b/C-5′ provided further evidence of the connectivity.
The relative configuration of subunits corresponded well with the data reported for the monomers in the literature. The NOE correlations of H-4′/H-5′ and H-5′/H-7′ clearly indicated the α-orientation of H-5′. Overhauser effects between H-13a/H-8, H-13a/H-5′, H-6′/H-15′ and H-13b/H-15′ are also in agreement with the proposed molecular formula of 7, named as lobatolide D.

3.1.8. Lobatolide G (8)

Compound 8 was isolated as a colorless oil with [α]28D +8 (c 0.1, CHCl3). Its HR-ESIMS spectrum showed quasimolecular ion at m/z 349.2017 [M + H]+, (calcd for 349.2015), indicating the molecular formula of C20H28O5. In the 1H and 13C NMR spectra of 8, signals of an isovaleroyl group were identified (δH 2.17 br s, 2.16 br s, 2.06 m, 0.93 d and 0.92 d; δC 172.3, 43.8, 25.6 and 2 × 22.6). In addition, the 13C and JMOD spectra suggested a sesquiterpene core consisting of one methyl, five methylenes, five methines and four quaternary carbons (Table 3). Two structural elements could be assigned based on the 1H–1H COSY spectrum: −CHR-CH2-CH2− (A δH 3.50, 1.81, 1.58, 2.35 and 2.13) (C-1−C-3) and −CHR-CHR-CHR-CHR-CH2− (B δH 2.23, 4.50, 2.79, 5.75, 2.31 and 1.56) (C-5−C-9). Structural parts A and B were connected based on the HMBC correlations observed between C-4 and H-3, and H-5; C-10 and H-2a, H-5, H-9b and H-14; C-12 and H-13a, and H-13b). The location of the ester group was also determined via 3JC,H coupling between the isovaleroyl CO (δC 172.3) and skeletal proton H-8 (δH 5.75). Comparison of 8 with sesquiterpenes previously isolated from the plant, good agreement was observed with 3β-acetoxy-8β-isovaleroyloxyreynosin with the exception of the missing acetyl substitution at C-3.
The relative configuration of compound 8 corresponded well with the data reported in the literature. The NOE correlations of H-1/H-5, H-2a and H-9b, H-5/H-7, H-7/H-8 and H-9b clearly indicated the α-orientation of H-1, H-5, H-7, H-8 and H-9b and the β-orientation of isovaleroyl group. The β-orientation of H-6 and H3-14 corroborated with diagnostic Overhauser-effects of the methyl protons with H-6. The coupling constant pattern of 8 is in good agreement with of that of structurally related compounds [50]. All of the above evidence proved the eudesmanolide structure 8β-isovaleroyloxyreynosin for this compound, and named as lobatolide G. ECD calculations were performed for (1R,5S,6R,7R,8R,10R)-8 similarly to the above compounds, but only partial agreements were found, and the experimental spectrum was weak and noisy even after measuring in various solvents and concentrations/settings. Therefore, the absolute configuration was verified by computing for the closely related derivative 9 (see below).
Besides new compounds, volenol (9) and (+)-aphanamol I (10) were also isolated from N. lobata. Their 1D and 2D NMR data were agreed with those published in the literature. The eudesmane-type volenol (9) was previously isolated from the aerial parts of Artemisia rubripes [51], while the isodaucane-type sesquiterpenoid (+)-aphanamol I (10) was first identified from the fruits of Aphanamixis grandifolia [52,53].
Since the absolute configuration of volenol (9) was not described in the literature, the above TDDFT-ECD protocol was applied to (1R,5S,6S,7S,10R)-9 (Figure 13). CAM-B3LYP re-optimization of the initial 24 MMFF conformers resulted in 12 low-energy conformers, the ECD spectra of which resembled well the experimental ECD spectrum (Figure 14). Therefore, the absolute configuration of 9 was verified in accordance with the biosynthetic considerations as (1R,5S,6S,7S,10R).
Naturally occurring germacranolides arise from all-trans farnesyl pyrophosphate by a series of biochemical transformations. Compounds belonging to this group (e.g., costunolide) primarily contain two endocyclic double bonds with all-trans configuration [e.g., costunolide derivatives]. However, in many cases cis configuration also occurs, e.g., in lobatolide C (3). The stereochemical skeletal types have distinct conformations, even allowing a certain amount of conformational flexibility [54]. In many cases, one of the endocyclic double bonds has undergone epoxidation as can be found in lobatolide E (6) or give a furan-type germacranolide occur in calyculatolide derivatives, like compound 7 (unit B), and 14.

3.2. Antiproliferative Activity of Isolated Compounds In Vitro

The antiproliferative activities of the isolated compounds 17, 9, 10 and neurolenin B (unit A of compound 7) were tested on three cervical cancer cell lines of different human papilloma virus (HPV) of different status (SiHa, HeLa, and C33A). Many of the compounds were additionally tested on non-cancerous cell lines NIH/3T3 (mouse embryonic fibroblast) and MRC-5 (human fibroblast) to obtain information concerning the cancer selectivity of the isolated natural products (Table 4). Compound 4 and neurolenin B exhibited higher potency than clinically utilized reference agent cisplatin while compounds 1, 3, 6 and 7 exerted similar antiproliferative properties to those of cisplatin. None of the presented natural products displayed substantial cancer selectivity exception for 6 which resulted in IC50 values in the range of 3.2–3.6 μM on C33A and SiHa cells but above 11 μM on the two fibroblast cell lines. Neurolenin B was proved to be the most active molecule with IC50 values against cancer cells lower than 1.4 μM on all cancer cell lines and NIH/3T3 while its calculated IC50 was slightly higher against human fibroblasts, indicating that the cancer selectivity is limited. Based on these antiproliferative activities 6 was selected for additional experiments. HeLa was shown to be the least sensitive cell line from the three cervical cancers. Since HeLa and SiHa cells result from HPV-mediated transformation, HPV-positivity does not explain the different responses of the two cell lines. The most invasive SiHa cell line was chosen further to investigate the possible mechanisms behind the antiproliferative effects. SiHa cells carry HPV-16 genes that contribute to the formation of metastases of this cancer type. HPV-16 oncoproteins promote cervical cancer invasiveness by upregulating specific matrix metalloproteinases (MMPs). HPV-positive cell lines express higher levels of MMP-2, membrane type 1-MMP and tissue inhibitor of metalloproteinase 2 [55].
Based on the antiproliferative assays, cell cycle analysis was performed on SiHa cells. Compound 6 caused a significant concentration-dependent accumulation of cells in the S phase to the detriment of G0/G1 phase after 48 or 72 h treatment (Figure 15). Cell cycle arrest in the S phase might refer to the cells having a higher tendency to enter S phase due to the dysfunction of cyclin dependent kinases, or it can be a consequence of inhibited DNA synthesis. The molecular mechanism of this cell cycle blockade requires further elucidation. Iso-seco-tanapartholide 3-O-methyl ether elicited similar S phase accumulation in HL-60 leukemia cells after 24 h of incubation at low concentration (5 μM). This phenomenon was accompanied by G2/M cell cycle phase arrest and subG1 accumulation, which could not be observed in the cell cycle profile of compound 6 [56]. G2/M phase accumulation appears to be a more characteristic effect in cervical cell lines treated by SLs. ROS generation, caspase-3 activation, inhibition of Bcl-2, and enhancement of Bax protein transcription seems a common background of in vitro anticancer activities [57,58]. A sesquiterpene lactone component of Elephantopus mollis (EM23) caused cell accumulation in the S phase parallel to the G2/M growth in K562 and HL-60 leukemia cells. ROS generation seems a grounded explanation for the pro-apoptotic effects since ROS quenching with N-acetyl-L-cysteine almost completely reversed these observations [59].
Compound 6 with high antiproliferative potency and modest selectivity to cancerous cells were chosen in order to assess its antimetastatic activity. The migration of cancer cells is an important property to form metastases. Therapeutical modalities that target the dissemination of tumor cells from the primary tumor could significantly improve the success of cancer therapy. Wound healing assay represents a simple tool to follow the antimigratory effects of SLs. Wound healing assay was performed on SiHa cell line with compound 6 at 1.5 µM and 3 µM concentrations. The inhibition of wound healing was evaluated by the measurement of the wound surface. 30% wound closure was detected after 24 h in the control wells, while 3 µM treatment with compound 6 resulted in the inhibition of wound healing by causing only 12% closure. After 48 h, the controls showed 70% wound closure; however, the 1.5 and 3 µM treated cells migrated only to 47% and 18% of the wound surface, respectively (Figure 16). These results follow the literature on the effects of sesquiterpenes. Parthenolide and other SLs inhibit cell migration by suppressing the FAK1 signaling pathway in breast cancer cultures [60]. Artesunate, an artemisinin derivative, promotes the downregulation of COX-2 expression associated with lymph node metastasis in cervical cancers. Conversely, COX-2 overexpression and consequent PGE2 formation can be responsible for the enhanced expression of metalloproteinases and vascular endothelial growth factors [61].

4. Conclusions

In this study, seven new sesquiterpenoids (13 and 58), together with three known ones (4, 9, and 10), were isolated from the dichloromethane phase of methanolic extract of the medicinal plant N. lobata. Four of the isolated compounds (14), and unit B of 7 are furanoheliangolide-type compounds, 5, 6 and unit A of 7 are germacranolides, 8 and 9 are eudesmanes, while 10 is an isodaucane. Compounds 1, 3, 4, and 68 contain an α-methylene-γ-lactone ring, while in case of 2 and 5, the exomethylene unit of the α-methylene-γ-lactone ring was modified to an oxymethylene part. Compounds are substituted with hydroxy, acetyl and isovaleroyl groups.
In the antiproliferative assay, lobatolides C (3), and E (6), neurolenin B (unit A of 7), and lobatolide F (7) were proved to be the most active compounds against HeLa, C33A, and SiHa cells. Although, neurolenin B showed the lowest IC50 values against cancer cells, it also inhibited normal cells at low concentration. Considering the selectivity of the compounds, lobatolide E (6) seems to be the most promising one. Moreover, this compound caused a significant, concentration-dependent accumulation of SiHa cells in the S phase during the cell cycle analysis investigation. The antimetastatic activity of compound 6 was also confirmed by wound healing assay.

Supplementary Materials

The following are available online at https://www.mdpi.com/article/10.3390/pharmaceutics13122088/s1, Figures S1–S55: 1D and 2D NMR, MS and ECD spectra of compounds 1–9. Tables S1–S3, S5, S6 and S8: Comparison of the experimental 13C NMR data. Table S4: Boltzmann populations and optical rotations of the low-energy conformers of (4R,6R,7S,8S,9R,10R,11R)-2 computed at various levels for the CAM-B3LYP/TZVP PCM/CHCl3 re-optimized MMFF conformers. Table S7: Comparison of the experimental 1H NMR data.

Author Contributions

Conceptualization, A.V., A.M., T.K. (Tibor Kurtán), G.K., R.F. and J.H.; methodology, A.V., A.M. and J.H.; software, A.M. and Z.S.; investigation, I.L., N.K., S.B.K., T.K. (Tibor Kovács), N.B., N.N., Z.S., A.M.; resources, T.K. (Tibor Kurtán), I.Z. and J.H.; data curation, A.V. and A.M.; writing—original draft preparation, A.V., A.M., Z.S. and J.H.; writing—review and editing, T.K. (Tibor Kurtán), I.Z. and J.H.; supervision, T.K. (Tibor Kurtán), and J.H.; project administration, J.H. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the Economic Development and Innovation Operative Programme (GINOP-2.3.2-15-2016-00012); the grant 20391-3/2018/FEKUSTRAT by the Ministry of Human Capacities, Hungary; and the National Research, Development and Innovation Office (NKFI FK 134653 and K 138672).

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Data sharing not applicable.

Acknowledgments

A.M. thanks for the János Bolyai Research Scholarship of the Hungarian Academy of Sciences. The authors thank to Imre Ocsovszki (Department of Biochemistry, Faculty of Medicine, University of Szeged, Dom ter 9., H-6720 Szeged, Hungary) for the flow cytometry investigations. The Governmental Information-Technology Development Agency (KIFÜ) is acknowledged for CPU time.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Chadwick, M.; Trewin, H.; Gawthrop, F.; Wagstaff, C. Sesquiterpenoids lactones: Benefits to plants and people. Int. J. Mol. Sci. 2013, 14, 12780–12805. [Google Scholar] [CrossRef] [Green Version]
  2. Rodriguez, E.; Towers, G.N.H.; Mitchell, J.C. Biological activities of sesquiterpene lactones. Phytochemistry 1976, 15, 1573–1580. [Google Scholar] [CrossRef]
  3. Seaman, F.C. Sesquiterpene lactones as taxonomic characters in the Asteraceae. Bot. Rev. 1982, 48, 121–595. [Google Scholar] [CrossRef]
  4. Ren, Y.; Yu, J.; Kinghorn, A.D. Development of anticancer agents from plant-derived sesquiterpene lactones. Curr. Med. Chem. 2016, 23, 2397–2420. [Google Scholar] [CrossRef] [PubMed]
  5. Zhumakayeva, A.; Rakhimov, K.; Sirota, V.; Arystan, L.; Madiyarov, A.; Adekenov, S. Long-term results of combination therapy for locally advanced breast cancer. Georgian Med. News 2018, 282, 30–35. [Google Scholar]
  6. Babaei, G.; Aliarab, A.; Abroon, S.; Rasmi, Y.; Gholizadeh, S.; Aziz, G. Application of sesquiterpene lactone: A new promising way for cancer therapy based on anticancer activity. Biomed Pharm. 2018, 106, 239–246. [Google Scholar] [CrossRef]
  7. Giron, L.M.; Freire, V.; Alonz, A.; Caceres, A. Ethnobotanical survey of the medicinal flora used by the Caribs of Guatemala. J. Ethnopharmacol. 1991, 34, 173–187. [Google Scholar] [CrossRef]
  8. Hartwell, J.L. Plants used against cancer. A survey. J. Nat. Prod. 1968, 31, 71–170. [Google Scholar]
  9. Passreiter, C.M.; Wendisch, D.; Gondol, D. Sesquiterpene lactones from Neurolaena Lobata. Phytochemistry 1995, 39, 133–137. [Google Scholar] [CrossRef]
  10. Borges-del-Castillo, J.; Manresa-Ferrero, M.T.; Rodríguez-Luis, F.; Vázquez-Bueno, P. Panama flora II. New sesquiterpene lactones from Neurolaena lobata. J. Nat. Prod. 1982, 45, 762–765. [Google Scholar] [CrossRef]
  11. Manchand, P.S.; Blount, J.F. Stereostructures of neurolenins A and B, novel germacranolide sesquiterpenes from Neurolaena lobata (L.) R.Br. J. Org. Chem. 1978, 43, 4352–4354. [Google Scholar] [CrossRef]
  12. Lajter, I.; Vasas, A.; Béni, Z.; Forgo, P.; Binder, M.; Bochkov, V.; Zupkó, I.; Krupitza, G.; Frisch, R.; Kopp, B.; et al. Sesquiterpenes from Neurolaena lobata and their antiproliferative and anti-inflammatory activities. J. Nat. Prod. 2014, 77, 576–582. [Google Scholar] [CrossRef]
  13. De Las Heras, B.; Slowing, K.; Benedi, J.; Carretero, E.; Ortega, T.; Toledo, C.; Bermejo, P.; Iglesias, I.; Abad, M.J.; Gomez-Serranillos, P.; et al. Antiinflammatory and antioxidant activity of plants used in traditional medicine in Ecuador. J. Ethnopharmacol. 1998, 61, 161–166. [Google Scholar] [CrossRef]
  14. Walshe-Roussel, B.; Choueiri, C.; Saleem, A.; Asim, M.; Caal, F.; Cal, V.; Rojas, M.O.; Pesek, T.; Durst, T.; Arnason, J.T. Potent anti-inflammatory activity of sesquiterpene lactones from Neurolaena lobata (L.) R. Br. ex Cass., a Q’eqchi’ Maya traditional medicine. Phytochemistry 2013, 92, 122–127. [Google Scholar] [CrossRef] [PubMed]
  15. Gracioso, J.S.; Hiruma-Lima, C.A.; Souza Brito, A.R.M. Antiulcerogenic effect of a hydroalcoholic extract and its organic fractions of Neurolaena lobata (L.) R.Br. Phytomedicine 2000, 7, 283–289. [Google Scholar] [CrossRef]
  16. Gracioso, J.S.; Paulo, M.Q.; Hiruma-Lima, C.A.; Souza Brito, A.R. Antinociceptive effect in mice of a hydroalcoholic extract of Neurolaena lobata (L.) R. Br. and its organic fractions. J. Pharm. Pharmacol. 1998, 50, 1425–1429. [Google Scholar] [CrossRef] [PubMed]
  17. Fujimaki, Y.; Kamachi, T.; Yanagi, T.; Caceres, A.; Maki, J.; Aoki, Y. Macrofilaricidal and microfilaricidal effects of Neurolaena lobata, a Guatemalan medicinal plant, on Brugia pahangi. J. Helminthol. 2005, 79, 23–28. [Google Scholar] [CrossRef]
  18. Chinchilla, M.; Valerio, I.; Sanchez, R.; Mora, V.; Bagnarello, V.; Martinez, L.; Gonzalez, A.; Vanegas, J.C.; Apestegui, A. In vitro antimalarial activity of extracts of some plants from a biological reserve in Costa Rica. Rev. Biol. Trop. 2012, 60, 881–891. [Google Scholar] [CrossRef] [Green Version]
  19. Bedoya, L.M.; Alvarez, A.; Bermejo, M.; Gonzalez, N.; Beltran, M.; Sanchez-Palomino, S.; Cruz, S.M.; Gaitan, I.; del Olmo, E.; Escarcena, R.; et al. Guatemalan plants extracts as virucides against HIV-1 infection. Phytomedicine 2008, 15, 520–524. [Google Scholar] [CrossRef] [PubMed]
  20. Unger, C.; Popescu, R.; Giessrigl, B.; Laimer, D.; Heider, S.; Seelinger, M.; Diaz, R.; Wallnofer, B.; Egger, G.; Hassler, M.; et al. The dichloromethane extract of the ethnomedicinal plant Neurolaena lobata inhibits NPM/ALK expression which is causal for anaplastic large cell lymphomagenesis. Int. J. Oncol. 2013, 42, 338–348. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  21. Unger, C.; Kiss, I.; Vasas, A.; Lajter, I.; Kramer, N.; Atanasov, A.G.; Nguyen, C.H.; Chatuphonprasert, W.; Brenner, S.; Krieger, S.; et al. The germacranolide sesquiterpene lactone neurolenin B of the medicinal plant Neurolaena lobata (L.) R.Br. ex Cass inhibits NPM/ALK-driven cell expansion and NF-κB-driven tumour intravasation. Phytomedicine 2015, 22, 862–874. [Google Scholar] [CrossRef] [PubMed]
  22. Kiss, I.; Unger, C.; Huu, C.N.; Atanasov, A.G.; Kramer, N.; Chatruphonprasert, W.; Brenner, S.; McKinnon, R.; Peschel, A.; Vasas, A.; et al. Lobatin B inhibits NPM/ALK and NF-κB attenuating anaplastic-large-cell-lymphomagenesis and lymphendothelial tumour intravasation. Cancer Lett. 2015, 356, 994–1006. [Google Scholar] [CrossRef]
  23. Francois, G.; Passreiter, C.M.; Woerdenbag, H.J.; van Looveren, M. Antiplasmodial activities and cytotoxic effects of aqueous extracts and sesquiterpene lactones from Neurolaena lobata. Planta Med. 1996, 62, 126–129. [Google Scholar] [CrossRef] [PubMed]
  24. MacroModel. Schrödinger, LLC. 2015. Available online: http://www.schrodinger.com/MacroModel (accessed on 30 September 2021).
  25. Frisch, M.J.; Trucks, G.W.; Schlegel, H.B.; Scuseria, G.E.; Robb, M.A.; Cheeseman, J.R.; Scalmani, G.; Barone, V.; Mennucci, B.; Petersson, G.A.; et al. Gaussian 09, Revision E. 01; Gaussian: Wallingford, CT, USA, 2013. [Google Scholar]
  26. Stephens, P.J.; Harada, N. ECD cotton effect approximated by the Gaussian curve and other methods. Chirality 2010, 22, 229–233. [Google Scholar] [CrossRef] [PubMed]
  27. CHESHIRE CCAT, the Chemical Shift Repository for Computed NMR Scaling Factors. Available online: http://cheshirenmr.info/index.htm (accessed on 30 September 2021).
  28. Lodewyk, M.W.; Siebert, M.R.; Tantillo, D.J. Computational prediction of 1H and 13C chemical shifts: A useful tool for natural product, mechanistic, and synthetic organic chemistry. Chem. Rev. 2012, 112, 1839–1862. [Google Scholar] [CrossRef] [PubMed]
  29. Varetto, U. MOLEKEL 5.4; Swiss National Supercomputing Centre: Manno, Switzerland, 2009. [Google Scholar]
  30. Kulmány, Á.E.; Frank, É.; Papp, D.; Szekeres, A.; Szebeni, G.J.; Zupkó, I. Biological evaluation of antiproliferative and anti-invasive properties of an androstadiene derivative on human cervical cancer cell lines. J. Steroid Biochem. Mol. Biol. 2021, 214, 105990. [Google Scholar] [CrossRef]
  31. Mosmann, T. Rapid colorimetric assay for cellular growth and survival: Application to proliferation and cytotoxicity assays. J. Immunol. Methods 1983, 65, 55–63. [Google Scholar] [CrossRef]
  32. Vermes, I.; Haanen, C.; Reutelingsperger, C. Flow cytometry of apoptotic cell death. J. Immunol. Methods 2000, 243, 167–190. [Google Scholar] [CrossRef]
  33. Latif, D.A.; Gonda, T.; Vágvölgyi, M.; Kúsz, N.; Kulmány, Á.; Ocsovszki, I.; Zomborszki, P.Z.; Zupkó, I.; Hunyadi, A. Synthesis and in vitro antitumor activity of naringenin oxime and oxime ether derivatives. Int. J. Mol. Sci. 2019, 20, 2184. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  34. Behrens, J.; Kameritsch, P.; Wallner, S.; Pohl, U.; Pogoda, K. The carboxyl tail of Cx43 augments p38 mediated cell migration in a gap junction-independent manner. Eur. J. Cell Biol. 2010, 89, 828–838. [Google Scholar] [CrossRef] [PubMed]
  35. Superchi, S.; Scafato, P.; Gorecki, M.; Pescitelli, G. Absolute configuration determination by quantum mechanical calculation of chiroptical spectra: Basics and applications to fungal metabolites. Curr. Med. Chem. 2018, 25, 287–320. [Google Scholar] [CrossRef] [PubMed]
  36. Maándi, A.; Kurtaán, T. Applications of OR/ECD/VCD to the structure elucidation of natural products. Nat. Prod. Rep. 2019, 36, 889–918. [Google Scholar] [CrossRef]
  37. Yanai, T.; Tew, D.P.; Handy, N.C. A new hybrid exchange-correlation functional using the Coulomb-attenuating method (CAM-B3LYP). Chem. Phys. Lett. 2004, 393, 51–57. [Google Scholar] [CrossRef] [Green Version]
  38. Li, W.S.; Yan, R.J.; Yu, Y.; Shi, Z.; Mándi, A.; Shen, L.; Kurtán, T.; Wu, J. Determination of the absolute configuration of super-carbon-chain compounds by a combined chemical, spectroscopic, and computational approach: Gibbosols A and B. Angew. Chem. Int. Ed. 2020, 59, 13028–13036. [Google Scholar] [CrossRef] [PubMed]
  39. Adamo, C.; Barone, V. Exchange functionals with improved long-range behavior and adiabatic connection methods without adjustable parameters: The mPW and mPW1PW models. J. Chem. Phys. 1998, 108, 664–675. [Google Scholar] [CrossRef]
  40. Qiu, S.; de Gussem, E.; Tehrani, K.A.; Sergeyev, S.; Bultinck, P.; Herrebout, W. Stereochemistry of the tadalafil diastereoisomers: A critical assessment of Vibrational Circular Dichroism, Electronic Circular Dichroism, and Optical Rotatory Dispersion. J. Med. Chem. 2013, 56, 8903–8914. [Google Scholar] [CrossRef]
  41. Dračínský, M.; Buděšínský, M.; Warżajtis, B.; Rychlewska, U. Solution and solid-state effects on NMR chemical shifts in sesquiterpene lactones: NMR, X-ray, and theoretical methods. J. Phys. Chem. A 2012, 116, 680–688. [Google Scholar] [CrossRef] [PubMed]
  42. Li, W.S.; Mándi, A.; Liu, J.J.; Shen, L.; Kurtán, T.; Wu, J. Xylomolones A–D from the Thai mangrove Xylocarpus moluccensis: Assignment of absolute stereostructures and unveiling a convergent strategy for limonoid biosynthesis. J. Org. Chem. 2019, 84, 2596–2606. [Google Scholar] [CrossRef] [PubMed]
  43. Liang, L.F.; Lan, L.F.; Taglialatela-Scafati, O.; Guo, Y.W. Sartrolides A–G and bissartrolide, new cembranolides from the South China Sea soft coral Sarcophyton trocheliophorum Marenzeller. Tetrahedron 2013, 69, 7381–7386. [Google Scholar] [CrossRef]
  44. Liang, L.F.; Kurtán, T.; Mándi, A.; Yao, L.G.; Li, J.; Lan, L.F.; Guo, Y.W. Structural, stereochemical, and bioactive studies of cembranoids from Chinese soft coral Sarcophyton trocheliophorum. Tetrahedron 2018, 74, 1933–1941. [Google Scholar] [CrossRef] [Green Version]
  45. Smith, S.G.; Goodman, J.M. Assigning stereochemistry to single diastereoisomers by GIAO NMR calculation: The DP4 probability. J. Am. Chem. Soc. 2010, 132, 12946–12959. [Google Scholar] [CrossRef] [PubMed]
  46. Grimblat, N.; Zanardi, M.M.; Sarotti, A.M. Beyond DP4: An improved probability for the stereochemical assignment of isomeric compounds using quantum chemical calculations of NMR shifts. J. Org. Chem. 2015, 80, 12526–12534. [Google Scholar] [CrossRef] [PubMed]
  47. Sun, P.; Cai, F.Y.; Lauro, G.; Tang, H.; Su, L.; Wang, H.L.; Li, H.H.; Mándi, A.; Kurtaán, T.; Riccio, R.; et al. Immunomodulatory biscembranoids and assignment of their relative and absolute configurations: Data set modulation in the Density Functional Theory/Nuclear Magnetic Resonance approach. J. Nat. Prod. 2019, 82, 1264–1273. [Google Scholar] [CrossRef] [PubMed]
  48. Passreiter, C.M.; Stöber, S.; Ortega, A. Furanoheliangolides from leaves of Neurolaena macrocephala. Z. Naturforsch. 2000, 55c, 1026–1029. [Google Scholar] [CrossRef] [PubMed]
  49. Bohlmann, F.; Jakupovic, J.; Ahmed, M.; Grenz, M.; Suding, H.; Robinson, H.; Kino, R.M. Germacranolides and diterpenes from Viguiera species. Phytochemistry 1981, 20, 113–116. [Google Scholar] [CrossRef]
  50. Scotta, R.; Zdero, C.; Bohlmann, F. Germacranolides, guaianolides and eudesmanolides from Greenmaniella resinosa. Phytochemistry 1987, 26, 1999–2006. [Google Scholar] [CrossRef]
  51. Lee, K.H.; Min, Y.D.; Choi, S.Z.; Kwon, H.C.; Cho, O.R.; Lee, K.C.; Lee, K.R. A new sesquiterpene lactone from Artemisia rubripes Nakai. Arch. Pharm. Res. 2004, 27, 1016–1019. [Google Scholar] [CrossRef] [PubMed]
  52. Nishizawa, M.; Inoue, A.; Hayashi, Y.; Sastrapradja, S.; Kosela, S.; Iwashita, T. Structure of aphanamol I and II. J. Org. Chem. 1984, 49, 3660–3662. [Google Scholar] [CrossRef]
  53. Hansson, T.; Wickberg, B. A short enantiospecific route to isodaucane sesquiterpenes from limonene. On the absolute configuration of (+)-aphanamol I and II. J. Org. Chem. 1992, 57, 5370–5376. [Google Scholar] [CrossRef]
  54. Krishnaswamy, N.R. Chemistry of Natural Products: A Unified Approach; Universities Press: Hyderabad, India, 1999; pp. 66–67. [Google Scholar]
  55. Schröpfer, A.; Kammerer, U.; Kapp, M.; Dietl, J.; Feix, S.; Anacker, J. Expression pattern of matrix metalloproteinases in human gynecological cancer cell lines. BMC Cancer 2010, 10, 553. [Google Scholar] [CrossRef] [Green Version]
  56. Molnár, J.; Szebeni, G.J.; Csupor-Löffler, B.; Hajdú, Z.; Szekeres, T.; Saiko, P.; Ocsovszki, I.; Puskás, L.G.; Hohmann, J.; Zupkó, I. Investigation of the antiproliferative properties of natural sesquiterpenes from Artemisia asiatica and Onopordum acanthium on HL-60 cells in vitro. Int. J. Mol. Sci. 2016, 17, 83. [Google Scholar] [CrossRef] [Green Version]
  57. Jeyamohan, S.; Moorthy, R.K.; Kannan, M.K.; Velanganni Arockiam, A.J. Parthenolide induces apoptosis and autophagy through the suppression of PI3K/Akt signaling pathway in cervical cancer. Biotech. Lett. 2016, 38, 1251–1260. [Google Scholar] [CrossRef]
  58. Abu-Izneid, T.; Abdur Rauf, A.; Shariati, M.A.; Khalil, A.A.; Imran, M.; Rebezov, M.; Uddin, S.; Mahomoodally, M.F.; Rengasamy, K.R.R. Sesquiterpenes and their derivatives-natural anticancer compounds: An update. Pharmacol. Res. 2020, 161, 105165. [Google Scholar] [CrossRef]
  59. Li, H.; Li, M.; Wang, G.; Shao, F.; Chen, W.; Xia, C.; Wang, S.; Li, Y.; Zhou, G.; Liu, Z. EM23, a natural sesquiterpene lactone from Elephantopus mollis, induces apoptosis in human myeloid leukemia cells through thioredoxin- and reactive oxygen species–mediated signalling pathways. Front. Pharmacol. 2016, 7, 77. [Google Scholar] [CrossRef] [Green Version]
  60. Berdan, A.C.; Ho, R.; Lehtola, H.S.; To, M.; Hu, X.; Huffman, T.R.; Petri, Y.; Altobelli, C.R.; Demeulenaere, S.G.; Olzmann, J.A.; et al. Parthenolide covalently targets and inhibits focal adhesion kinase in breast cancer cells. Cell. Chem. Biol. 2019, 26, 1027–1035. [Google Scholar] [CrossRef]
  61. Zhang, L.; Qian, H.; Sha, M.; Luan, Z.; Lin, M.; Yuan, D.; Li, X.; Huang, J.; Ye, L. Downregulation of HOTAIR ex-pression mediated antimetastatic effect of artesunate on cervical cancer by inhibiting COX-2 expression. PLoS ONE 2016, 11, e0164838. [Google Scholar] [CrossRef] [Green Version]
Figure 1. Structures of compounds 110 isolated from N. lobata.
Figure 1. Structures of compounds 110 isolated from N. lobata.
Pharmaceutics 13 02088 g001
Figure 2. Comparison of the experimental ECD spectrum of 1 measured in MeCN with the CAM-B3LYP/TZVP PCM/MeCN spectrum of (4R,6R,7R,8S,9R,10R)-1 (level of optimization: CAM-B3LYP/TZVP PCM/MeCN). The bars represent the rotational strength values of the lowest energy conformer.
Figure 2. Comparison of the experimental ECD spectrum of 1 measured in MeCN with the CAM-B3LYP/TZVP PCM/MeCN spectrum of (4R,6R,7R,8S,9R,10R)-1 (level of optimization: CAM-B3LYP/TZVP PCM/MeCN). The bars represent the rotational strength values of the lowest energy conformer.
Pharmaceutics 13 02088 g002
Figure 3. Overlapped geometries of the twelve low-energy CAM-B3LYP/TZVP PCM/MeCN conformers of (4R,6R,7R,8S,9R,10R)-1.
Figure 3. Overlapped geometries of the twelve low-energy CAM-B3LYP/TZVP PCM/MeCN conformers of (4R,6R,7R,8S,9R,10R)-1.
Pharmaceutics 13 02088 g003
Figure 4. The characteristic H-8–H-11 NOE correlation indicated on the lowest-energy conformers of the (a) (4R,6R,7S,8S,9R,10R,11S)-2 (3.50 Å) and (b) (4R,6R,7S,8S,9R,10R,11R)-2 (2.19 Å) stereoisomers. Level of optimization: B3LYP/6-31+G(d,p).
Figure 4. The characteristic H-8–H-11 NOE correlation indicated on the lowest-energy conformers of the (a) (4R,6R,7S,8S,9R,10R,11S)-2 (3.50 Å) and (b) (4R,6R,7S,8S,9R,10R,11R)-2 (2.19 Å) stereoisomers. Level of optimization: B3LYP/6-31+G(d,p).
Pharmaceutics 13 02088 g004
Figure 5. Comparison of the experimental ECD spectrum of 3 measured in MeCN with the CAM-B3LYP/TZVP PCM/MeCN spectrum of (6R,7R,8S,9R,10R)-3 (level of optimization: CAM-B3LYP/TZVP PCM/MeCN). The bars represent the rotational strength values of the lowest-energy conformer.
Figure 5. Comparison of the experimental ECD spectrum of 3 measured in MeCN with the CAM-B3LYP/TZVP PCM/MeCN spectrum of (6R,7R,8S,9R,10R)-3 (level of optimization: CAM-B3LYP/TZVP PCM/MeCN). The bars represent the rotational strength values of the lowest-energy conformer.
Pharmaceutics 13 02088 g005
Figure 6. Overlapped geometries of the conformers of the two conformer groups of (6R,7R,8S,9R,10R)-3. Group A: Confs. A-F, I, J; group B: Confs. G, H, K-M. Level of optimization: CAM-B3LYP/TZVP PCM/MeCN.
Figure 6. Overlapped geometries of the conformers of the two conformer groups of (6R,7R,8S,9R,10R)-3. Group A: Confs. A-F, I, J; group B: Confs. G, H, K-M. Level of optimization: CAM-B3LYP/TZVP PCM/MeCN.
Pharmaceutics 13 02088 g006
Figure 7. Comparison of the experimental ECD spectrum of 3 measured in MeCN with the CAM-B3LYP/TZVP PCM/MeCN spectra of conformers A and G of (6R,7R,8S,9R,10R)-3, as the lowest-energy representatives of groups A and B. Level of optimization: CAM-B3LYP/TZVP PCM/MeCN.
Figure 7. Comparison of the experimental ECD spectrum of 3 measured in MeCN with the CAM-B3LYP/TZVP PCM/MeCN spectra of conformers A and G of (6R,7R,8S,9R,10R)-3, as the lowest-energy representatives of groups A and B. Level of optimization: CAM-B3LYP/TZVP PCM/MeCN.
Pharmaceutics 13 02088 g007
Figure 8. Comparison of the experimental ECD spectrum of 4 measured in MeCN with the CAM-B3LYP/TZVP PCM/MeCN spectrum of (6R,7S,8S,9R,10R)-4 (level of optimization: CAM-B3LYP/TZVP PCM/MeCN). The bars represent the rotational strength values of the lowest energy conformer.
Figure 8. Comparison of the experimental ECD spectrum of 4 measured in MeCN with the CAM-B3LYP/TZVP PCM/MeCN spectrum of (6R,7S,8S,9R,10R)-4 (level of optimization: CAM-B3LYP/TZVP PCM/MeCN). The bars represent the rotational strength values of the lowest energy conformer.
Pharmaceutics 13 02088 g008
Figure 9. Comparison of the experimental ECD spectrum of 5 measured in MeCN with the BH&HLYP/TZVP PCM/MeCN spectrum of (4R,6R,7S,8S,9R,10R,11S)-5 (level of optimization: CAM-B3LYP/TZVP PCM/MeCN). The bars represent the rotational strength values of the lowest energy conformer.
Figure 9. Comparison of the experimental ECD spectrum of 5 measured in MeCN with the BH&HLYP/TZVP PCM/MeCN spectrum of (4R,6R,7S,8S,9R,10R,11S)-5 (level of optimization: CAM-B3LYP/TZVP PCM/MeCN). The bars represent the rotational strength values of the lowest energy conformer.
Pharmaceutics 13 02088 g009
Figure 10. Overlapped geometries of the twelve low-energy CAM-B3LYP/TZVP PCM/MeCN conformers of (4R,6R,7S,8S,9R,10R,11S)-5.
Figure 10. Overlapped geometries of the twelve low-energy CAM-B3LYP/TZVP PCM/MeCN conformers of (4R,6R,7S,8S,9R,10R,11S)-5.
Pharmaceutics 13 02088 g010
Figure 11. Comparison of the experimental ECD spectrum of 6 measured in MeCN with the PBE0/TZVP PCM/MeCN spectrum of (1R,3S,6R,7R,8R,10R)-6 (level of optimization: CAM-B3LYP/TZVP PCM/MeCN). The bars represent the rotational strength values of the lowest energy conformer.
Figure 11. Comparison of the experimental ECD spectrum of 6 measured in MeCN with the PBE0/TZVP PCM/MeCN spectrum of (1R,3S,6R,7R,8R,10R)-6 (level of optimization: CAM-B3LYP/TZVP PCM/MeCN). The bars represent the rotational strength values of the lowest energy conformer.
Pharmaceutics 13 02088 g011
Figure 12. Overlapped geometries of the ten low-energy CAM-B3LYP/TZVP PCM/MeCN conformers of (1R,3S,6R,7R,8R,10R)-6.
Figure 12. Overlapped geometries of the ten low-energy CAM-B3LYP/TZVP PCM/MeCN conformers of (1R,3S,6R,7R,8R,10R)-6.
Pharmaceutics 13 02088 g012
Figure 13. Comparison of the experimental ECD spectrum of 9 measured in MeCN with the BH&HLYP/TZVP PCM/MeCN spectrum of (1R,5S,6S,7S,10R)-9 (level of optimization: CAM-B3LYP/TZVP PCM/MeCN). The bars represent the rotational strength values of the lowest energy conformer.
Figure 13. Comparison of the experimental ECD spectrum of 9 measured in MeCN with the BH&HLYP/TZVP PCM/MeCN spectrum of (1R,5S,6S,7S,10R)-9 (level of optimization: CAM-B3LYP/TZVP PCM/MeCN). The bars represent the rotational strength values of the lowest energy conformer.
Pharmaceutics 13 02088 g013
Figure 14. Overlapped geometries of the twelve low-energy CAM-B3LYP/TZVP PCM/MeCN conformers of (1R,5S,6S,7S,10R)-9.
Figure 14. Overlapped geometries of the twelve low-energy CAM-B3LYP/TZVP PCM/MeCN conformers of (1R,5S,6S,7S,10R)-9.
Pharmaceutics 13 02088 g014
Figure 15. Cell cycle distributions of SiHa cells after treatment with compound 6 for 48 (upper pane) or 72 h (lower panel). Distribution of cell populations in different cell cycle phases. *, ** and *** indicate p < 0.05, p < 0.01 and p < 0.001, respectively, by means of one-way ANOVA followed by Dunnett’s post-hoc test.
Figure 15. Cell cycle distributions of SiHa cells after treatment with compound 6 for 48 (upper pane) or 72 h (lower panel). Distribution of cell populations in different cell cycle phases. *, ** and *** indicate p < 0.05, p < 0.01 and p < 0.001, respectively, by means of one-way ANOVA followed by Dunnett’s post-hoc test.
Pharmaceutics 13 02088 g015
Figure 16. Effect of compound 6 on the migration of SiHa cancer cells after 24 and 48 h of incubation (left and right upper panels, respectively). ** and *** indicate p < 0.01 and p < 0.001, respectively, by means of one-way ANOVA followed by Dunnett’s post-hoc test. Wound healing assay. Lower panel: representative images of reduced wound healing at 0, 24 and 48 h post-treatment.
Figure 16. Effect of compound 6 on the migration of SiHa cancer cells after 24 and 48 h of incubation (left and right upper panels, respectively). ** and *** indicate p < 0.01 and p < 0.001, respectively, by means of one-way ANOVA followed by Dunnett’s post-hoc test. Wound healing assay. Lower panel: representative images of reduced wound healing at 0, 24 and 48 h post-treatment.
Pharmaceutics 13 02088 g016
Table 1. 1H NMR data of compounds 16 in CDCl3 at 500 MHz (δ in ppm, mult. J in Hz).
Table 1. 1H NMR data of compounds 16 in CDCl3 at 500 MHz (δ in ppm, mult. J in Hz).
Position123456
1 2.93 dd (11.3, 2.3)
25.60 s5.52 s5.64 s5.65 s6.58 d (11.7)2.35 m, 1.64 m
3 5.97 t (11.7)5.50 dd (11.4, 5.7)
43.01 m3.01 m 3.09 m
52.59 ddd (14.2, 9.3, 4.8) (α)
2.05 d (14.2) (β)
2.63 ddd (14.2, 9.6, 7.0)
2.03 m
6.00 dd (4.2, 1.8)5.99 d (3.7, 1.8)1.79 dt (5.1, 12.3)
1.63 dt (5.1, 12.1)
5.66 d (9.6)
64.78 dd (9.3, 4.7)4.29 dd (9.6, 6.9)5.71 m5.24 m4.56 dd (12.1, 5.3)5.22 t (9.4)
73.30 m3.16 dd (8.9, 6.9)3.57 dd (7.6, 4.0)3.86 dd (4.6, 1.6)3.04 ddd (10.3, 8.1, 5.3)3.21 br d (9.4)
84.13 t (4.7)5.07 d (5.4)4.04 m5.01 dd (5.2, 1.6)5.57 d (9.6)5.74 d (5.7)
95.27 d (4.7)3.94 d (5.4)5.27 d (3.6)5.31 d (5.2)5.57 d (9.6)2.63 dd (15.2, 5.7)
1.40 dd (15.2, 1.7)
11 2.89 m 2.23 br d (8.1)
13a6.30 d (3.0)3.74 dd (9.4, 4.1)6.29 d (3.0)6.35 d (3.0)3.68 dd (10.2, 5.3)6.23 d (3.5)
13b5.35 d (3.0)3.67 m5.35 d (3.0)5.48 d (3.0)3.34 t (10.2)5.65 s
141.37 s1.45 s1.41 s1.43 s1.31 s1.21 s
151.40 d (7.0)1.34 d (7.3)2.06 s2.08 s1.15 d (6.3)1.87 d (1.3)
OMe 3.40 s 3.39 s
iVal CO 1′
2a′2.38 dd (15.1, 7.2)2.14 dd (15.3, 7.2)2.37 dd (15.2, 7.3)2.15 m2.00 m2.27 m
2b′2.33 dd (15.1, 7.2)2.10 m2.31 dd (15.2, 7.2)2.12 m2.00 m2.22 m
3′2.19 m2.02 m2.16 m2.01 m2.07 m2.07 m
4′1.02 d (6.7)0.93 d (7.0)1.01 d (6.6)0.92 d (7.1)0.91 d (6.7)0.95 d (6.4)
5′1.02 d (6.7)0.91 d (6.9)1.01 d (6.6)0.90 d (7.1)0.89 d (6.7)0.94 d (6.4)
8-OH2.82 d (4.5) 4.05 s
3-OAc 2.10 s
9-OAc 2.23 s2.11 s
Table 2. 13C NMR data of compounds 16 in CDCl3 at 125 MHz (δC in ppm).
Table 2. 13C NMR data of compounds 16 in CDCl3 at 125 MHz (δC in ppm).
Position123456
1203.9204.2203.9202.3204.765.6
2104.0103.8104.6104.3125.631.4
3192.8192.4185.8185.4147.677.1
431.531.4130.9134.628.7145.3
540.741.7134.3134.438.1123.9
675.375.474.975.577.475.4
747.946.145.244.937.553.3
873.875.175.975.768.968.5
977.174.377.573.674.843.3
1088.990.788.388.579.361.4
11141.349.3141.2139.040.2138.2
12169.2174.2170.0168.7174.2171.2
13122.672.1123.2123.866.3121.7
1418.918.418.418.623.620.2
1516.116.019.621.219.813.0
OMe 59.6 59.2
iVal CO 1′171.5171.7171.7169.3171.2173.2
2′43.143.043.142.642.944.1
3′25.725.625.625.324.726.7
4′22.422.522.422.422.422.7
5′22.422.622.422.422.522.7
Ac CO 168.8170.4171.8
Ac Me 21.320.620.8
Table 3. 1H (500 MHz) and 13C (125 MHz) NMR data of compounds 7 and 8 in CDCl3 (δ in ppm).
Table 3. 1H (500 MHz) and 13C (125 MHz) NMR data of compounds 7 and 8 in CDCl3 (δ in ppm).
Position7 (Unit A)7 (Unit B)8
δH, Mult. (J in Hz)δCδH, Mult. (J in Hz)δCδH, Mult. (J in Hz)δC
1-205.0-203.83.50 dd (11.4, 4.7)78.8
26.22 d (11.6)123.85.64 s105.41.81 m, 1.58 m31.1
36.04 t (11.6)147.3-193.02.35 m, 2.13 m33.5
42.14 m25.02.81 m36.8-142.0
52.13 m, 1.96 m38.93.54 m44.82.23 d (11.0)53.7
65.07 d (5.0)79.64.44 t (4.9)75.04.50 t (11.1)75.2
7-154.23.68 m45.62.79 dd (11.1, 2.7)52.2
86.38 d (8.9)66.95.20 br d (5.1)76.15.75 d (2.7)65.9
95.69 d (8.9)73.14.05 br s74.112.31 dd (15.2, 2.2), 1.56 m40.6
10-80.7-90.7-42.8
11-135.7-138.0-134.8
12-172.2-168.9-170.0
13a3.05 dd (14.5, 8.8)23.06.29 d (3.0),124.16.15 d (3.2)119.6
13b2.48 dd (14.5, 6.8) 5.74 d (3.0) 5.44 d (3.0)
141.32 s23.91.48 s19.25.01 s, 4.94 s111.1
151.14 d (6.4)22.51.40 d (6.8)10.00.96 s13.8
iVal-CO 1′s172.0-171.5-172.3
2′2.27 dd (15.8, 6.8), 2.15 m (2H)43.12.08 m (2H)42.82.17 br s, 2.16 br s43.8
3′2.02 m25.51.96 m25.32.06 m25.6
4′0.94 d (6.7)22.50.90 d (6.7)22.40.93 d (4.2)22.6
5′0.93 d (6.7)22.50.87 d (6.7)22.40.92 d (4.2)22.6
9-OAc 170.3
2.15 s20.7
Table 4. Antiproliferative activity of the isolated compounds determined by MTT assay. Growth inhibition and calculated IC50 values of the tested compounds and their 95% confidence intervals (CI).
Table 4. Antiproliferative activity of the isolated compounds determined by MTT assay. Growth inhibition and calculated IC50 values of the tested compounds and their 95% confidence intervals (CI).
Comp.Conc. (μM)Growth Inhibition Values (%) ± SEM at 10 and 30 μM
Calculated IC50 Values (μM) [95% CI]
HeLaC33ASiHaNIH/3T3MRC-5
11020.39 ± 0.6290.00 ± 2.5477.93 ± 0.56n.t. an.t.
3081.78 ± 1.7091.79 ± 0.2292.10 ± 0.73
IC5016.78 [15.57–18.09]5.89 [2.32–5.52]4.86 [4.45–5.31]
210<20 b<20<20n.t.n.t.
3021.55 ± 2.4236.60 ± 2.95<20
31091.42 ± 0.3479.15 ± 3.3360.65 ± 1.5892.96 ± 0.5688.13 ± 1.22
3096.58 ± 1.1794.46 ± 0.2178.77 ± 1.2398.15 ± 0.1397.58 ± 0.10
IC504.76 [4.27–5.32]6.82 [5.97–7.79]7.81 [6.81–8.96]4.12 [3.68–4.61]6.67 [5.35–8.30]
41083.94 ± 1.8097.27 ± 0.1879.77 ± 1.47n.t.n.t.
3090.10 ± 1.0397.19 ± 0.4194.62 ± 0.27
IC505.10 [5.10–6.51]2.05 [2.91–2.18]2.22 [1.92–2.56]
510<2050.50 ± 3.49<20n.t.n.t.
3063.86 ± 3.2389.95 ± 1.465.52 ± 2.52
610<2096.55 ± 0.2284.19 ± 2.3240.67 ± 2.92<20
3091.88 ± 1.2396.72 ± 0.3384.76 ± 1.8996.80 ± 0.2573.05 ± 1.64
IC5015.69 [12.51–19.69]3.53 [2.83–4.39]3.24 [2.95–3.55]11.17 [10.64–11.74]11.94 [10.62–13.43]
71091.87 ± 0.7294.66 ± 0.1977.17 ± 1.0490.65 ± 1.20<20
3096.44 ± 0.1095.18 ± 0.3193.50 ± 0.4298.06 ± 0.1997.63 ± 0.27
IC506.00 [5.44–6.62]6.62 [4.38–10.00]6.02 [5.53–6.54]5.68 [5.36–6.02]15.42 [12.71–18.71]
910<20<20<20n.t.n.t.
3037.27 ± 2.6526.69 ± 2.7335.69 ± 1.29
1010<20<20<20n.t.n.t.
3036.49 ± 1.55<20<20
Neu B1094.00 ± 0.5998.29 ± 0.0388.73 ± 0.3698.11 ± 0.0994.26 ± 0.23
3094.81 ± 0.3797.91 ± 0.1593.42 ± 0.3398.02 ± 0.0597.46 ± 0.21
IC501.24 [1.07–1.45]0.47 [4.59–4.93]1.33 [1.14–1.54]1.31 [1.22–1.41]5.18 [4.69–5.73]
Cispl.1032.23 ± 1.1677.17 ± 1.2160.98 ± 0.9273.88 ± 0.4671.24.± 2.89
3093.70 ± 0.8297.50 ± 0.1388.95 ± 0.5397.10 ± 0.1570.65 ± 1.34
IC5012.14 [10.18–14.46]5.85 [5.37–6.38]4.29 [3.72–4.95]5.50 [4.46–6.35]5.77 [4.30–7.74]
a Not tested. b Cancer cell growth inhibition values less than 20% were considered negligible and are not given numerically.
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Vasas, A.; Lajter, I.; Kúsz, N.; Király, S.B.; Kovács, T.; Kurtán, T.; Bózsity, N.; Nagy, N.; Schelz, Z.; Zupkó, I.; et al. Isolation, Structure Determination of Sesquiterpenes from Neurolaena lobata and Their Antiproliferative, Cell Cycle Arrest-Inducing and Anti-Invasive Properties against Human Cervical Tumor Cells. Pharmaceutics 2021, 13, 2088. https://doi.org/10.3390/pharmaceutics13122088

AMA Style

Vasas A, Lajter I, Kúsz N, Király SB, Kovács T, Kurtán T, Bózsity N, Nagy N, Schelz Z, Zupkó I, et al. Isolation, Structure Determination of Sesquiterpenes from Neurolaena lobata and Their Antiproliferative, Cell Cycle Arrest-Inducing and Anti-Invasive Properties against Human Cervical Tumor Cells. Pharmaceutics. 2021; 13(12):2088. https://doi.org/10.3390/pharmaceutics13122088

Chicago/Turabian Style

Vasas, Andrea, Ildikó Lajter, Norbert Kúsz, Sándor Balázs Király, Tibor Kovács, Tibor Kurtán, Noémi Bózsity, Nikolett Nagy, Zsuzsanna Schelz, István Zupkó, and et al. 2021. "Isolation, Structure Determination of Sesquiterpenes from Neurolaena lobata and Their Antiproliferative, Cell Cycle Arrest-Inducing and Anti-Invasive Properties against Human Cervical Tumor Cells" Pharmaceutics 13, no. 12: 2088. https://doi.org/10.3390/pharmaceutics13122088

APA Style

Vasas, A., Lajter, I., Kúsz, N., Király, S. B., Kovács, T., Kurtán, T., Bózsity, N., Nagy, N., Schelz, Z., Zupkó, I., Krupitza, G., Frisch, R., Mándi, A., & Hohmann, J. (2021). Isolation, Structure Determination of Sesquiterpenes from Neurolaena lobata and Their Antiproliferative, Cell Cycle Arrest-Inducing and Anti-Invasive Properties against Human Cervical Tumor Cells. Pharmaceutics, 13(12), 2088. https://doi.org/10.3390/pharmaceutics13122088

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