Acute myocardial infarction (AMI) is a leading cause of death and disability worldwide, especially in developing countries, where its rates are increasing [
1]. This is largely due to lifestyle changes, poor diets, and lack of medical resources [
1]. Given the concerning rise in AMI incidence, it is crucial to explore drugs specifically aimed at improving the condition of patients with heart disease. Recent research indicates that the natural product formononetin can help protect heart cells by reducing oxidative stress and promoting cell survival, which can lessen damage from ischemia-reperfusion injury [
2,
3]. Natural products have been proven to be valuable source of bioactive molecules for drug development [
4], continued research in this area could provide more therapeutic options for patients suffering from cardiovascular issues.
Daphniphyllum alkaloids (DAs) are a diverse group of natural products from the
Daphniphyllum genus, characterized by unique azapolycyclic structures with multiple continuous stereogenic centers [
5]. These alkaloids exhibit a variety of scaffold types, with this diversity arising from their distinctive biosynthesis, which involves the repeated fission of C–C and/or C–N bonds, followed by rearrangements, decarboxylation, and other reactions [
5]. Approximately 360 natural DAs have been reported thus far, some of which exhibit significant biological effects, such as cytotoxic and anti-human immunodeficiency virus (HIV) activities [
5–
9]. Challenging and caged polycyclic architectures and the promising biological profiles make DAs intriguing targets for natural products and synthetic chemistry [
10–
16].
Our group has a longstanding interest in exploring medicinal plants to identify natural products with bioactive properties [
7,
30–
32]. The plants of
D. calycinum have a long history of use in traditional Chinese medicine for treating fever, asthma, inflammation, and influenza [
7,
9]. Given this medicinal significance and inspired by previous studies revealing the significant structural diversity of DAs in
D. calycinum [
7,
17–
18], we set out to investigate this species and discovered four novel alkaloids (daphcalycines A–D;
Fig. 1), as well as their biogenetic precursor (caldaphnidine C). Compound
1 possesses a unique 13-oxa-17-aza-pentacyclo[7.6.4.1
12,15.0
4,8.0
9,15] eicosane core, with the linkage between C-3 and C-14 and the disruption of the C-2–C-18 bond being unprecedented in previously reported DAs. Compounds
2–
4 have a rearranged 6/5 bicyclic (rings A and B), while
2 also possesses a novel C
23N skeleton generated through the ring formation of C-23/C-15 as well as eleven continuous stereogenic centers. Remarkably, compound
2 markedly enhanced the survival of H9c2 cardiomyocytes under oxygen glucose deprivation and reoxygenation (OGD/R) conditions. Preliminary mechanistic study revealed that
2 displayed cardioprotective activity through affecting the kelch-like ECH-associated protein 1-nuclear factor erythroid 2-related factor 2/heme oxygenase-1 (Keap1-Nrf2/HO-1) antioxidant pathway, which could provide a candidate molecule for the further development of treatments for AMI. Herein, the structure elucidation, proposed biosynthetic pathway, and bioactivity evaluation of compounds
1–
4 are presented.
The molecular formula of compound 1, C21H29NO4, was established from high resolution electrospray ionization mass spectroscopy (HRESIMS) (m/z 360.2175 [M + H]+, calcd. 360.2169), corresponding to eight indices of hydrogen deficiency (IHD). The 1H nuclear magnetic resonance (NMR) spectrum showed one methyl group [δH 2.09 (3H, s)] and one olefinic proton (δH 5.82). Combined with its distortionless enhancement by polarization transfer (DEPT), heteronuclear single quantum correlation (HSQC), and heteronuclear multiple-bond correlation (HMBC) spectra, the 13C NMR data of 1 (Table S1 in Supporting information) revealed 21 carbon signals, including two carbonyls (δC 174.5 and 204.4), two olefinic carbons (δC 153.0 and 130.0). Apart from the aforementioned four carbons, the remaining seventeen carbon signals were assigned to three non-protonated carbons (one hemiacetal carbon at δC 106.1), two methines, eleven methylenes (one oxygenated carbon at δC 76.2, two nitrogen-bearing carbons at δC 60.5 and 54.3), and one methyl carbon (δC 27.4). Considering the presence of two carbonyl groups and one olefinic group, the occurrence of five unassigned IHD indicated that compound 1 possessed a pentacyclic framework.
Compound
1 was determined to feature an unprecedented 13-oxa-17-aza-pentacyclo[7.6.4.1
12,15.0
4,8.0
9,15] eicosane core by detailed analysis of its 2D NMR spectrum. Two independent spin coupling systems could be deduced from the
1H–
1H correlation spectroscopy (COSY) spectrum of
1 (
Fig. 2). The HMBC correlations from H
2–1 to C-2, C-5, C-8 and C-9; from H
2–7 to C-2 and C-5; from H-6 and H
2–12 to C-5; from H-15 to C-8 and C-10, coupled with the comparison of chemical shifts of the congeners with similar fragment [
19–
21], sufficiently facilitated the construction of a 7/7/5 (rings C–E) tricyclic systems. Furthermore, the HMBC correlations from H
2–13 to C-3, C-5, and C-8; from H
2–14 to C-3, C-4, and C-8; and from H
2–4 to C-3, C-5, C-8, and C-14 revealed the presence of cyclohexane structure (ring A). In addition, a 2-oxopropyl side chain was inferred to be attached to the nitrogen atom based on the chemical shift of C-19 (
δC 60.5) and the HMBC correlations from H
3–20 to C-18 and C-19 together with from H
2–19 to C-2 and C-7. Lastly, the HMBC correlations from H
2–21 to C-3, C-4, C-5, and C-6, combined with the remaining one IHD and the chemical shifts of C-3 (
δC 106.1) and C-21 (
δC 76.2) disclosed the presence of a furan ring (ring B). The planar structure of compound
1 was thereby established.
The relative configuration of compound
1 was assigned through analysis of the nuclear Overhauser effect spectroscopy (NOESY) spectrum (
Fig. 2), where H
2–21 at
δH 4.40 and 3.97 correlated with H
2–13a and H-6, respectively, confirming the relative configurations of C-5, C-6, and C-8, with H-6 designated as
β-oriented. Due to overlapping signals between H-10 and H
2–4a at
δH 2.84 and between H
2–1a and H
2–19b at
δH 3.66 in pyridine-
d5, we employed methanol-
d4 as the solvent for NMR testing of compound
1, which provided distinct signals for these protons (Table S1). The NOESY correlation of H
2–21a and H-10 confirmed that H-10 was also
β-oriented. Additionally, the presence of a methylene bridge between C-3 and C-5 rendered the 6-oxabicyclo[3.2.1] octane (rings A and B) in compound
1 a rigid structure, naturally fixing the orientation of the C-3–C-4 bond in the
α-direction. The inference was supported by the NOESY correlations of H
2–4a /H
2–1a and H
2–7a, as well as H
2–4b /H
2–7a and H-6. Finally, the absolute configurations of compound
1 were established using electronic circular dichroism (ECD) calculation. As illustrated in Fig. S1 (Supporting information), the calculated ECD curves for compound
1 closely aligned with the experimental one, allowing for the assignment of its absolute configuration, as depicted in
Fig. 1.
The HRESIMS data of compound 2 revealed a molecular formula of C23H33NO6, corresponding to eight IHD. The 13C NMR spectrum resolved 23 carbon resonances corresponding to five quaternary carbons (two oxygenated carbons at δC 101.8 and 79.7, and one carbonyl carbon at δC 173.6), seven methines (one oxygenated carbon at δC 73.3, one nitrogen-bearing carbon at δC 84.0, and one hemiacetal carbon at δC 104.9), ten methylenes (one oxygenated carbon at δC 75.2, two nitrogen-bearing carbons at δC 59.9 and 54.3), and one methyl, which were distinguished with the aid of HSQC and DEPT experiments. The above-mentioned information suggested compound 2 to be a heptacyclic alkaloid.
The 1H–1H COSY spectrum revealed the presence of three spin systems in compound 2 (Fig. S2 in Supporting information). In the HMBC spectrum, the correlations from H-1 to C-2, C-3, C-4, C-8, C-9, C-19, and C-23; from H2–4 to C-1, C-2, C-3, C-5, C-6, and C-8; from H-6 to C-5 and C-8; from H-10 to C-9 and C-15; and from H2–23 to C-9 and C-19 led to the establishment of a partial structure of rings A–E. Furthermore, the HMBC correlations from H2–21 to C-5, C-6, C-7, and C-8, coupled with the chemical shifts of C-7 (δC 104.9) and C-21 (δC 75.2) supported the presence of a furan ring (ring F). Lastly, the HMBC correlations from H2–13 to C-8, C-9, and C-22, along with the remaining one IHD and the chemical shifts of C-9 (δC 101.8) and C-22 (δC 173.6) revealed the existence of a δ-lactone (ring G). Hence, compound 2 was successfully constructed as a novel 6/5/7/5/6/5/6 seven-ring skeleton.
The network of NOESY correlations including H
2–13/H
2–21b, H-1, and H-15, H
2–21b/H-6, H
2–19b/H-1 and H
3–20 revealed that these protons were on the same side of the molecular plane and were arbitrarily assigned in
β-orientation. Therefore, the NOESY correlations of H
2–4a/H-2 and H-10, as well as H
2–4b/H-7 showed that they were
α-oriented. The relative configuration of OH-3 could not be accurately determined due to insufficient NOESY signals. Fortunately, after numerous attempts, crystals of compound
2 suitable for X-ray diffraction experiments were successfully obtained from an optimized binary solvent system (MeOH/CH
2Cl
2/H
2O, 9:1:0.1). The X-ray structure (
Fig. 3) not only confirmed the planar structure inferred from the NMR data but also definitively established its relative configuration. As shown in Fig. S1, the ECD calculation of (1
R, 2
S, 3
R, 5
R, 6
S, 7
S, 8
R, 9
R, 10
S, 15
S, 18
S)-
2 was consistent with the experimental result, which assigned its absolute configuration.
The molecular formula of C
22H
29NO
5 for compound
3 was determined by the NMR data and (+)-HRESIMS ion at
m/
z 388.2123 [M + H]
+ (calcd. 388.2118), indicative of nine IHD. The 1D NMR data of
3 (Table S2 in Supporting information) closely resembled that of the known compound daphnimacropodine A [
21], with the primary difference being that the 5/6 bicyclic system (rings A and B) in daphnimacropodine A was substituted with a 6/5 bicyclic system (rings A and B) in
3. The HMBC correlations from H-1 to C-2, C-3, C-4, C-8, and C-19; from H
2–4 to C-2, C-3, C-5, and C-8; and from H
3–20 to C-2, along with the
1H–
1H COSY correlations of H
3–20/H-18/H
2–19 enabled the establishment of the rearranged rings A and B. The partial relative configuration of compound
3 was elucidated by its NOESY data (Fig. S3 in Supporting information). The single crystal was successfully obtained from methanol solvent at room temperature, enabling the single-crystal X-ray diffraction analysis (Cu K
α) to confirm the planar structure and the absolute configuration as depicted in
Fig. 1, with a Flack’s parameter of 0.07(7).
The structure of compound 4 was established by the comprehensive spectroscopic data analyses. Detailed structural elucidation was described in the Supporting information.
A plausible biosynthetic pathway for compounds
1–
4 was proposed in
Scheme 1. Initially, compound
5 (caldaphnidine C) could be converted to
1iii through a series of hydroxylation, dehydration, and oxidation reactions. Following this,
1iv was formed through an Mannich-like reaction and simultaneous proton migration [
13,
22]. Next, the C-2–C-3 bond broke through further dehydration and oxidation reactions. The cyclohexane (ring A) structure was produced through subsequent transformation steps, including aldol and decarboxylation reactions [
23–
26]. Finally, through oxidation and intramolecular hemiacetal reactions [
21], compound
1 was generated (path A). Compounds
2–
4 were proposed to be derived from compound
5 and its C-10 isomer (path B). Starting from the hydroxylation at positions C-2, C-3, C-7, and C-21 [
21], the key intermediate
2i could be yielded. The rearranged rings A and B were then formed through a pinacol rearrangement reaction [
27]. Subsequently, cleavage of the C-7–N bond and intramolecular hemiacetal reaction could lead to the generation of
2iv [
21]. From this critical intermediate
2iv, two distinct pathways diverge: on the one hand, intermediate
2iv underwent a cascade of
N-methylation, oxidation, and Polonovski reactions to generate
2vi [
24,
25], which underwent aza-Prins, esterification, and reduction reactions to form compound
2 [
28,
29]. On the other hand, intermediate
2iv underwent formation of the C-22–N bond and oxidation reactions to generate compound
3. Following this, compound
3 was subjected to double bond isomerization [
33], oxidation, and dehydration reactions to yield
2ix. Ultimately, through keto-enol tautomerization [
34,
35], oxidation, reduction, and ethylation reactions, compound
4 was generated.
We evaluated the protective effects of the isolated compounds on H9c2 cardiomyocytes under OGD/R conditions. We first assessed cell viability using the CCK-8 assay, and the results showed that compound
2 significantly enhanced cell survival (
Fig. 4A). Further dose-dependent experiments revealed that the protective effect of compound
2 increased with concentration, demonstrating a clear dose-response relationship (
Fig. 4B), suggesting its potential as a cardioprotective agent. To further investigate whether compound
2 could inhibit OGD/R-induced reactive oxygen species (ROS) production, we measured intracellular ROS levels using flow cytometry and the fluorescent probe DCFH-DA. The results indicated that compound
2 significantly reduced ROS generation in H9c2 cells following OGD/R treatment (
Figs. 4C and
D), indicating its ability to mitigate oxidative stress and protect cardiomyocytes from hypoxic injury. Additionally, we examined the impact of compound
2 on apoptosis using flow cytometry. The data showed that compound
2 markedly decreased the proportion of early and late apoptotic cells (
Figs. 4E and
F), demonstrating its capacity to inhibit apoptosis and further supporting its protective role. To elucidate the molecular mechanisms underlying the protective effects of compound
2, we performed Western blot analysis, focusing on key proteins in the Keap1/Nrf2/HO-1 signaling pathway. The results revealed that compound
2 dose-dependently inhibited Keap1 protein expression, stabilized Nrf2 protein levels, and promoted HO-1 expression (
Figs. 4G–
J). These findings suggested that compound
2 exerted its protective effects by activating the Nrf2/HO-1 antioxidant pathway, thereby enhancing cellular antioxidant capacity and alleviating oxidative stress induced by hypoxia.
In conclusion, we have described four novel alkaloids from D. calycinum with intricate and diverse polycyclic skeletons, substantially enriching the structural diversity of DAs. The biosynthetic pathways for the four unusual alkaloids have been also proposed. Additionally, we initially discovered that DAs possessed the potential for cardioprotective activity, and related mechanism studies have been conducted. The discovery of daphcalycines A–D (1–4) not only enriches the diversity of DAs skeletons but also opens up new perspectives for development on DAs as cardioprotective agents.