*Communication* **Phytochemical Analysis of the Fruits of Sea Buckthorn (***Hippophae rhamnoides***): Identification of Organic Acid Derivatives**

**Yong Hoon Lee 1,†, Hee Joo Jang 1,†, Kun Hee Park <sup>2</sup> , Seon-Hee Kim <sup>3</sup> , Jung Kyu Kim <sup>4</sup> , Jin-Chul Kim <sup>5</sup> , Tae Su Jang 6,\* and Ki Hyun Kim 1,\***


**Abstract:** *Hippophae rhamnoides* L. (Elaeagnaceae), commonly known as "Sea buckthorn" and "Vitamin tree", is a spiny deciduous shrub whose fruit is known for its nutritional composition, such as vitamin C, and is consumed as a dietary supplement worldwide. As part of our ongoing efforts to identify structurally new and bioactive constituents from natural resources, the phytochemical investigation of the extract of *H. rhamnoides* fruits led to the isolation of one malate derivative (**1**), five citrate derivatives (**2–6**), and one quinate derivative (**7**). The structures of the isolated compounds were elucidated by analysis of 1D and 2D nuclear magnetic resonance (NMR) spectroscopic data and high-resolution electrospray ionization (HR-ESI) liquid chromatography–mass spectrometry (LC/MS) data. Three of the citrate derivatives were identified as new compounds: (*S*)-1-butyl-5 methyl citrate (**3**), (*S*)-1-butyl-10 -methyl citrate (**4**), and (*S*)-1-methyl-10 -butyl citrate (**6**), which turned out to be isolation artifacts. The absolute configurations of the new compounds were established by quantum chemical electronic circular dichroism (ECD) calculation, which is an informative tool for verifying the absolute configuration of organic acid derivatives. The isolated compounds **1–7** were evaluated for their stimulatory effects on osteogenesis. Compounds **1**, **3**, **4**, **6**, and **7** stimulated osteogenic differentiation up to 1.4 fold, compared to the negative control. These findings provide experimental evidence that active compounds **1**, **3**, **4**, **6**, and **7** induce the osteogenesis of mesenchymal stem cells and activate bone formation.

**Keywords:** *Hippophae rhamnoides*; Elaeagnaceae; citrate derivatives; nuclear magnetic resonance (NMR); electronic circular dichroism (ECD)

#### **1. Introduction**

*Hippophae rhamnoides* L. (Elaeagnaceae), known as "Sea buckthorn" and "Vitamin tree", is a spiny deciduous shrub. It has historically been used for medicinal purposes worldwide, especially in Central and Southeastern Asia. In China, its pharmacological use was revealed more than a thousand years ago through the records of Sibu Yidian from the Tang Dynasty and Jing Zhu Ben Cao from the Qing dynasty [1]. Sea buckthorn oil was listed in the Russian Pharmacopeia in the 1950s and used as an anti-inflammatory and wound healing agent in Russia [2]. Its bark and fruits have also been used as Tibetan

**Citation:** Lee, Y.H.; Jang, H.J.; Park, K.H.; Kim, S.-H.; Kim, J.K.; Kim, J.-C.; Jang, T.S.; Kim, K.H. Phytochemical Analysis of the Fruits of Sea Buckthorn (*Hippophae rhamnoides*): Identification of Organic Acid Derivatives. *Plants* **2021**, *10*, 860. https://doi.org/10.3390/ plants10050860

Academic Editors: Narendra Singh Yadav and Jong Seong Kang

Received: 5 April 2021 Accepted: 23 April 2021 Published: 24 April 2021

**Publisher's Note:** MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

**Copyright:** © 2021 by the authors. Licensee MDPI, Basel, Switzerland. This article is an open access article distributed under the terms and conditions of the Creative Commons Attribution (CC BY) license (https:// creativecommons.org/licenses/by/ 4.0/).

folk medicine to treat pulmonary disorders, cough, fever, and tumors [3]. For the past 50 years, several medicinal preparations based on *H. rhamnoides* have been clinically used to treat radiation damage, oral inflammation, burns, and gastric ulcers, and more than 300 preparations have been reported in the literature [4].

*H. rhamnoides* bears yellowish or orange berries, which are commonly found as a powdered drink mix in markets as a functional food and are used to make fruit sauce and wine [5]. Vitamin C is a major nutritional component of *H. rhamnoides* fruits [6], the level of which exceeds that of lemons and oranges. Previous pharmacological studies of *H. rhamnoides* have reported that its extracts exhibited therapeutic properties, including anti-platelet effects via the inhibitory mechanism of thrombin-activated platelets and antimicrobial effects against *S. aureus* and *C. albicans* through the inhibition of adhesion and biofilm formation [7,8]. In addition, previous phytochemical investigations on *H. rhamnoides* have demonstrated the presence of diverse types of bioactive substances, including vitamins (A, E, K, riboflavin, and folic acid), carotenoids, phytosterols, minerals, organic acids, polyunsaturated fatty acids, and amino acids [9,10]. Because of the nutritional importance of *H. rhamnoides* fruits, many research groups have attempted to investigate the phytochemical constituents of these fruits, in which flavonoids, proanthocyanidins, alkaloids, and phenolic acids have been identified [11,12]. Nevertheless, previous findings suggested that only few studies have been carried out on the chemical constituents of *H. rhamnoides* fruits, despite their health benefits.

In this context, as a part of continuous research programs to discover structurally new bioactive compounds from diverse natural sources [13–19], our group has focused on the potential bioactive composition of *H. rhamnoides* fruits. In our recent study of *H. rhamnoides* fruits, we identified six chemical compounds, including one citric acid derivative, two flavonoids, one phenolic compound, and two megastigmane derivatives, which were evaluated for their anti-inflammatory effects by determining LPS-induced NO production in RAW 264.7 cells [20]. As a result, we found that a citric acid derivative, 1,5-dimethyl citrate, effectively inhibited LPS-induced NO production and simultaneously showed antiinflammatory effects by inhibiting the expression of pro-inflammatory mediators, iNOS and COX-2, and the activity of pro-inflammatory cytokines, IL-6 and TNF-α [20]. Our recent findings provide experimental evidence that the organic acid derivative present in *H. rhamnoides* fruits could function as a promising bioactive substance. In the present study, further phytochemical investigation of the extract of *H. rhamnoides* fruits was performed to identify potential bioactive organic acid constituents based on recent preliminary findings. Phytochemical analysis of the BuOH-soluble fraction derived from the extract led to the isolation of seven organic acid derivatives: one malate derivative (**1**), five citrate derivatives (**2–6**), three of which were identified as new compounds, and one quinate derivative (**7**). The structures of the isolated organic acid derivatives (**1–7**) were determined by analysis of 1D and 2D nuclear magnetic resonance (NMR) spectroscopic data and high-resolution electrospray ionization (HR-ESI) liquid chromatography–mass spectrometry (LC/MS). In addition, the absolute configurations of the new compounds were established using quantum chemical electronic circular dichroism (ECD) calculations. Herein, we describe the isolation and structural elucidation of the organic acid derivatives **1–7** and the evaluation of their biological activities for stimulatory effects on osteogenesis.

#### **2. Results and Discussion**

#### *2.1. Isolation of Compounds*

The aqueous extract powder of *H. rhamnoides* fruits was suspended for solvent partitioning in water and then fractionated with four organic solvents, namely, hexane, dichloromethane, ethyl acetate, and *n*-butanol. LC/MS-based analysis combined with our in-house built UV library, and TLC analysis of the solvent-partitioned fractions suggested that the BuOH-soluble fraction contained the majority of the organic acid derivatives. Phytochemical analysis of the BuOH-soluble fraction was performed using repeated column chromatography with silica gel 60, RP-C<sup>18</sup> silica gel, Diaion HP-20, Sephadex LH-20,

and high-performance liquid chromatography (HPLC). The final semi-preparative HPLC separation afforded seven organic acid derivatives (**1–7**) (Figure 1). umn chromatography with silica gel 60, RP-C<sup>18</sup> silica gel, Diaion HP-20, Sephadex LH-20, and high-performance liquid chromatography (HPLC). The final semi-preparative HPLC separation afforded seven organic acid derivatives (**1–7**) (Figure 1).

that the BuOH-soluble fraction contained the majority of the organic acid derivatives. Phytochemical analysis of the BuOH-soluble fraction was performed using repeated col-

*Plants* **2021**, *10*, x FOR PEER REVIEW 3 of 11

**Figure 1.** Chemical structure of compounds **1–7**. **Figure 1.** Chemical structure of compounds **1–7**.

#### *2.2. Elucidation of Compound Structures 2.2. Elucidation of Compound Structures*

Compound **3** was isolated as a yellowish gum. The molecular formula was determined to be C11H18O<sup>7</sup> from the protonated molecular ion peak [M + H]<sup>+</sup> at *m*/*z* 263.1122 (calculated for C11H19O7, 263.1131) in the positive-ion mode of HRESIMS (Figure S1). The <sup>1</sup>H NMR spectrum (Figure S2) of **3** (Table 1) showed the presence of characteristic signals of an oxygenated methylene group at <sup>H</sup> 4.08 (2H, t, *J* = 6.5 Hz), a methoxy group at <sup>H</sup> 3.66 (3H, s), two pairs of relatively deshielded methylene groups overlapped at <sup>H</sup> 2.84 and 2.92, another two pairs of methylene groups at <sup>H</sup> 1.61 (2H, m) and 1.40 (2H, m), and a terminal methyl group at <sup>H</sup> 0.94 (3H, t, *J* = 7.5 Hz). The 13C NMR data of **3** (Table 1) assigned with the assistance of the HSQC experiment (Figure S4) confirmed the 11 carbon signals, which corresponded to three carboxyl carbons (*δ*<sup>C</sup> 170.3, 170.5, and 175.3), two oxygenated carbons (*δ*<sup>C</sup> 64.2 and 72.9), four sets of methylene groups (*δ*<sup>C</sup> 18.6, 30.2, 42.7, and 42.7), one methoxy-oriented carbon (*δ*<sup>C</sup> 50.7), and one methyl carbon (*δ*<sup>C</sup> 12.5). Based on this evidence, it was predicted that compound **3** was a citric acid derivative, since the characteristic NMR data of **3** was similar to that of the citric acid derivatives [21,22]. The gross structure of **3** was determined by 2D NMR experiments (1H-<sup>1</sup>H COSY; Figure S3 and HMBC; Figure S5). The HMBC correlations of H2-2/C-1, H2-2/C-1', H2-2/C-4, H2-4/C-1', H2- 4/C-2, and H2-4/C-5 verified the citric acid moiety. Furthermore, the HMBC correlation of 5-OCH<sup>3</sup> with C-5 confirmed that the methoxy group was attached at C-5 of the citric acid moiety, and that of H2-1'' with C-1, along with 1H-<sup>1</sup>H cross-peaks between H2-1''/ H2-2''/ H2-3''/ H3-4'', revealed that a butyl group was attached at the other end of C-1 (Figure 2). Thus, the planar structure of **3** was established, as shown in Figure 2. Finally, the absolute configuration at C-3 was assigned using quantum chemical ECD calculations [23–25]. Two possibilities with absolute configurations of 3*R* and 3*S* were calculated. The weighted calculated ECD spectrum (green line) of **3***S* (3*S*) agreed well (although not identical) with the experimental curve of **3**, compared to that (red line) of **3***R* (3*R*) (Figure 3), thus determining the absolute configuration of C-3 as *S*. The chemical structure of **3**, including its absolute Compound **3** was isolated as a yellowish gum. The molecular formula was determined to be C11H18O<sup>7</sup> from the protonated molecular ion peak [M + H]<sup>+</sup> at *m*/*z* 263.1122 (calculated for C11H19O7, 263.1131) in the positive-ion mode of HRESIMS (Figure S1). The <sup>1</sup>H NMR spectrum (Figure S2) of **3** (Table 1) showed the presence of characteristic signals of an oxygenated methylene group at *δ*<sup>H</sup> 4.08 (2H, t, *J* = 6.5 Hz), a methoxy group at *δ*<sup>H</sup> 3.66 (3H, s), two pairs of relatively deshielded methylene groups overlapped at *δ*<sup>H</sup> 2.84 and 2.92, another two pairs of methylene groups at *δ*<sup>H</sup> 1.61 (2H, m) and 1.40 (2H, m), and a terminal methyl group at *δ*<sup>H</sup> 0.94 (3H, t, *J* = 7.5 Hz). The <sup>13</sup>C NMR data of **3** (Table 1) assigned with the assistance of the HSQC experiment (Figure S4) confirmed the 11 carbon signals, which corresponded to three carboxyl carbons (*δ*<sup>C</sup> 170.3, 170.5, and 175.3), two oxygenated carbons (*δ*<sup>C</sup> 64.2 and 72.9), four sets of methylene groups (*δ*<sup>C</sup> 18.6, 30.2, 42.7, and 42.7), one methoxy-oriented carbon (*δ*<sup>C</sup> 50.7), and one methyl carbon (*δ*<sup>C</sup> 12.5). Based on this evidence, it was predicted that compound **3** was a citric acid derivative, since the characteristic NMR data of **3** was similar to that of the citric acid derivatives [21,22]. The gross structure of **3** was determined by 2D NMR experiments (1H-1H COSY; Figure S3 and HMBC; Figure S5). The HMBC correlations of H2-2/C-1, H2-2/C-10 , H2-2/C-4, H2-4/C-10 , H2-4/C-2, and H2-4/C-5 verified the citric acid moiety. Furthermore, the HMBC correlation of 5-OCH<sup>3</sup> with C-5 confirmed that the methoxy group was attached at C-5 of the citric acid moiety, and that of H2-1" with C-1, along with <sup>1</sup>H-1H cross-peaks between H2-1"/ H2-2"/ H2-3"/ H3-4", revealed that a butyl group was attached at the other end of C-1 (Figure 2). Thus, the planar structure of **3** was established, as shown in Figure 2. Finally, the absolute configuration at C-3 was assigned using quantum chemical ECD calculations [23–25]. Two possibilities with absolute configurations of 3*R* and 3*S* were calculated. The weighted calculated ECD spectrum (green line) of **3***S* (3*S*) agreed well (although not identical) with the experimental curve of **3**, compared to that (red line) of **3***R* (3*R*) (Figure 3), thus determining the absolute configuration of C-3 as *S*. The chemical structure of **3**, including its absolute configuration, was elucidated as (*S*)-1-butyl-5-methyl citrate.

configuration, was elucidated as (*S*)-1-butyl-5-methyl citrate. **Figure 2.** Key 1H-<sup>1</sup>H COSY ( ) and HMBC ( ) correlations of compounds **3**, **4**, and **6**. Compound **4** was obtained as a white amorphous powder. The molecular formula was determined to be the same as that of compound **3**, C11H18O7, from HRESIMS data (Figure S6), which showed a protonated molecular ion peak [M + H]<sup>+</sup> at *m*/*z* 263.1119 (calculated for C11H19O7, 263.1131). The <sup>1</sup>H (Figure S7) and <sup>13</sup>C NMR data (Figure S9) of **4** (Table 1) were also almost identical to those of compound **3**, except for apparent differences in the chemical shifts responsible for the citric acid moiety—the <sup>1</sup>H chemical shift of the methoxy group (*δ*<sup>H</sup> 3.66 in **3**, *δ*<sup>H</sup> 3.75 in **4**). Detailed analysis of the 2D NMR spectra (1H-1H COSY; Figure S8 and HMBC; Figure S10) assigned the gross planar structure of **4** (Figure 2), where the key HMBC correlation of 10 -OCH<sup>3</sup> with C-1<sup>0</sup> was observed, indicating that the methoxy group was attached at C-10 in compound **4** (Figure 2). Finally, to confirm the absolute configuration of C-3 in **4**, the computationally calculated ECD data of the

two enantiomers, 3*R-* and 3*S*-isomers, were compared with the experimental ECD data of **4** (Figure 4). The ECD data of **4** showed a good fit with the calculated ECD spectrum (green line) of **4***S* (3*S*) (Figure 4). This finding allowed us to assign the structure of **4** as (*S*)-1-butyl-10 -methyl citrate. *2.2. Elucidation of Compound Structures* Compound **3** was isolated as a yellowish gum. The molecular formula was determined to be C11H18O<sup>7</sup> from the protonated molecular ion peak [M + H]<sup>+</sup> at *m*/*z* 263.1122 (calculated for C11H19O7, 263.1131) in the positive-ion mode of HRESIMS (Figure S1). The *2.2. Elucidation of Compound Structures* Compound **3** was isolated as a yellowish gum. The molecular formula was determined to be C11H18O<sup>7</sup> from the protonated molecular ion peak [M + H]<sup>+</sup> at *m*/*z* 263.1122 (calculated for C11H19O7, 263.1131) in the positive-ion mode of HRESIMS (Figure S1). The *2.2. Elucidation of Compound Structures* Compound **3** was isolated as a yellowish gum. The molecular formula was determined to be C11H18O<sup>7</sup> from the protonated molecular ion peak [M + H]<sup>+</sup> at *m*/*z* 263.1122 (calculated for C11H19O7, 263.1131) in the positive-ion mode of HRESIMS (Figure S1). The

<sup>1</sup>H NMR spectrum (Figure S2) of **3** (Table 1) showed the presence of characteristic signals

<sup>1</sup>H NMR spectrum (Figure S2) of **3** (Table 1) showed the presence of characteristic signals

<sup>1</sup>H NMR spectrum (Figure S2) of **3** (Table 1) showed the presence of characteristic signals

<sup>1</sup>H COSY, HSQC,

<sup>1</sup>H

that the BuOH-soluble fraction contained the majority of the organic acid derivatives. Phytochemical analysis of the BuOH-soluble fraction was performed using repeated column chromatography with silica gel 60, RP-C<sup>18</sup> silica gel, Diaion HP-20, Sephadex LH-20, and high-performance liquid chromatography (HPLC). The final semi-preparative HPLC

that the BuOH-soluble fraction contained the majority of the organic acid derivatives. Phytochemical analysis of the BuOH-soluble fraction was performed using repeated column chromatography with silica gel 60, RP-C<sup>18</sup> silica gel, Diaion HP-20, Sephadex LH-20, and high-performance liquid chromatography (HPLC). The final semi-preparative HPLC

that the BuOH-soluble fraction contained the majority of the organic acid derivatives. Phytochemical analysis of the BuOH-soluble fraction was performed using repeated column chromatography with silica gel 60, RP-C<sup>18</sup> silica gel, Diaion HP-20, Sephadex LH-20, and high-performance liquid chromatography (HPLC). The final semi-preparative HPLC


**Table 1.** <sup>1</sup>H (850 MHz) and <sup>13</sup>C (212.5 MHz) NMR data of compounds **3**, **4**, and **6** in CD3OD <sup>a</sup> . of an oxygenated methylene group at <sup>H</sup> 4.08 (2H, t, *J* = 6.5 Hz), a methoxy group at <sup>H</sup> 3.66 of an oxygenated methylene group at <sup>H</sup> 4.08 (2H, t, *J* = 6.5 Hz), a methoxy group at <sup>H</sup> 3.66 of an oxygenated methylene group at <sup>H</sup> 4.08 (2H, t, *J* = 6.5 Hz), a methoxy group at <sup>H</sup> 3.66

*Plants* **2021**, *10*, x FOR PEER REVIEW 3 of 11

*Plants* **2021**, *10*, x FOR PEER REVIEW 3 of 11

*Plants* **2021**, *10*, x FOR PEER REVIEW 3 of 11

separation afforded seven organic acid derivatives (**1–7**) (Figure 1).

separation afforded seven organic acid derivatives (**1–7**) (Figure 1).

separation afforded seven organic acid derivatives (**1–7**) (Figure 1).

**Figure 1.** Chemical structure of compounds **1–7**.

**Figure 1.** Chemical structure of compounds **1–7**.

**Figure 1.** Chemical structure of compounds **1–7**.

<sup>a</sup> Coupling constants (in Hz) are given in parentheses. Assignments based on <sup>1</sup>H-1H COSY, HSQC, and HMBC. <sup>b</sup> Overlapped. the absolute configuration of C-3 as *S*. The chemical structure of **3**, including its absolute configuration, was elucidated as (*S*)-1-butyl-5-methyl citrate. the absolute configuration of C-3 as *S*. The chemical structure of **3**, including its absolute configuration, was elucidated as (*S*)-1-butyl-5-methyl citrate. the absolute configuration of C-3 as *S*. The chemical structure of **3**, including its absolute configuration, was elucidated as (*S*)-1-butyl-5-methyl citrate. 1' - 175.3 - 174.1 - 173.9 <sup>b</sup> 1'-OCH<sup>3</sup> - - 3.75 s 51.5 - - 1-OCH<sup>3</sup> - - - - 3.67 s 50.2

**Figure 2.** Key <sup>1</sup>H- **Figure 2.** Key <sup>1</sup>H COSY ( ) and HMBC ( ) correlations of compounds **3**, **4**, and **6**. <sup>1</sup>H-1H COSY ( **Figure 2.** Key 1H-<sup>1</sup>H COSY ( ) ) and HMBC ( and HMBC ( ) correlations of compounds **3**, **4**, and **6**. **Figure 2.** Key 1H-<sup>1</sup>H COSY ( ) and HMBC ( ) correlations of compounds ) correlations of compounds **3 3**,, **4 4** , and , and **6 6**.. and HMBC. <sup>b</sup>Overlapped.

Compound **4** was obtained as a white amorphous powder. The molecular formula

<sup>H</sup> 3.75 in **4**). Detailed analysis of the 2D NMR spectra (1H-

(Figure S6), which showed a protonated molecular ion peak [M + H]<sup>+</sup> at *m*/*z* 263.1119 (calculated for C11H19O7, 263.1131). The 1H (Figure S7) and 13C NMR data (Figure S9) of **4** (Table 1) were also almost identical to those of compound **3**, except for apparent differences in the chemical shifts responsible for the citric acid moiety—the 1H chemical shift of the

COSY; Figure S8 and HMBC; Figure S10) assigned the gross planar structure of **4** (Figure 2), where the key HMBC correlation of 1′-OCH<sup>3</sup> with C-1′ was observed, indicating that the methoxy group was attached at C-1′ in compound **4** (Figure 2). Finally, to confirm the absolute configuration of C-3 in **4**, the computationally calculated ECD data of the two enantiomers, 3*R-* and 3*S*-isomers, were compared with the experimental ECD data of **4**  (Figure 4). The ECD data of **4** showed a good fit with the calculated ECD spectrum (green line) of **4***S* (3*S*) (Figure 4). This finding allowed us to assign the structure of **4** as (*S*)-1-

**Figure 3.** Experimental and calculated ECD spectra of **3**. **Figure 3.** Experimental and calculated ECD spectra of **3**.

<sup>H</sup> 3.66 in **3**,

methoxy group (

butyl-1'-methyl citrate.

**Figure 4.** Experimental and calculated ECD spectra of **4**. **Figure 4.** Experimental and calculated ECD spectra of **4**. **Figure 4.** Experimental and calculated ECD spectra of **4**.

Compound **6**, obtained as a white amorphous powder, was determined to have the molecular formula C11H18O7, similar to those of compounds **3** and **4**, from the protonated molecular ion peak [M + H]<sup>+</sup> at *m*/*z* 263.1145 (calculated for C11H19O7, 263.1131) in the positive-ion mode of HRESIMS (Figure S11). The 1H (Figure S12) and 13C NMR data (Figure S14) (Table 1) of **6** were also quite similar to those of compounds **3** and **4**, with the only noticeable differences being that the 13C chemical shift of C-1 was deshielded to *δ*<sup>C</sup> 173.9, compared to these two compounds (*δ*<sup>C</sup> 170.2 in **3**, *δ*<sup>C</sup> 170.0 in **4**), and the 1H chemical shift of H-1'' was deshielded to *δ*<sup>H</sup> 4.18, compared to **3** and **4** (*δ*<sup>H</sup> 4.08 in **3**, *δ*<sup>H</sup> 4.07 in **4**). This finding suggested the possibility that the butyl group in **6** was attached to the different carboxyl groups of compounds **3** and **4**. The gross planar structure of **6** was verified by 2D NMR analysis (1H-<sup>1</sup>H COSY; Figure S13 and HMBC; Figure S15) (Figure 2). The key HMBC correlation between 1-OCH<sup>3</sup> and C-1 supported the presence of the methoxy group at C-1, and the location of the butyl group was determined to be C-1' by the key HMBC correlation of H-1'' and C-1' (Figure 2). Finally, the absolute configuration of **6** was assigned using quantum chemical ECD calculations. As shown in Figure 5, the experimental ECD curve of **6** is consistent with the calculated ECD spectrum (green line) of **6***S* (3*S*). Therefore, the chemical structure of **6** was characterized as (*S*)-1-methyl-1'-butyl citrate. Compound **6**, obtained as a white amorphous powder, was determined to have the molecular formula C11H18O7, similar to those of compounds **3** and **4**, from the protonated molecular ion peak [M + H]<sup>+</sup> at *m*/*z* 263.1145 (calculated for C11H19O7, 263.1131) in the positive-ion mode of HRESIMS (Figure S11). The <sup>1</sup>H (Figure S12) and <sup>13</sup>C NMR data (Figure S14) (Table 1) of **6** were also quite similar to those of compounds **3** and **4**, with the only noticeable differences being that the <sup>13</sup>C chemical shift of C-1 was deshielded to *δ*<sup>C</sup> 173.9, compared to these two compounds (*δ*<sup>C</sup> 170.2 in **3**, *δ*<sup>C</sup> 170.0 in **4**), and the <sup>1</sup>H chemical shift of H-1" was deshielded to *δ*<sup>H</sup> 4.18, compared to **3** and **4** (*δ*<sup>H</sup> 4.08 in **3**, *δ*<sup>H</sup> 4.07 in **4**). This finding suggested the possibility that the butyl group in **6** was attached to the different carboxyl groups of compounds **3** and **4**. The gross planar structure of **6** was verified by 2D NMR analysis (1H-1H COSY; Figure S13 and HMBC; Figure S15) (Figure 2). The key HMBC correlation between 1-OCH<sup>3</sup> and C-1 supported the presence of the methoxy group at C-1, and the location of the butyl group was determined to be C-10 by the key HMBC correlation of H-1" and C-10 (Figure 2). Finally, the absolute configuration of **6** was assigned using quantum chemical ECD calculations. As shown in Figure 5, the experimental ECD curve of **6** is consistent with the calculated ECD spectrum (green line) of **6***S* (3*S*). Therefore, the chemical structure of **6** was characterized as (*S*)-1-methyl-10 -butyl citrate. Compound **6**, obtained as a white amorphous powder, was determined to have the molecular formula C11H18O7, similar to those of compounds **3** and **4**, from the protonated molecular ion peak [M + H]<sup>+</sup> at *m*/*z* 263.1145 (calculated for C11H19O7, 263.1131) in the positive-ion mode of HRESIMS (Figure S11). The 1H (Figure S12) and 13C NMR data (Figure S14) (Table 1) of **6** were also quite similar to those of compounds **3** and **4**, with the only noticeable differences being that the 13C chemical shift of C-1 was deshielded to *δ*<sup>C</sup> 173.9, compared to these two compounds (*δ*<sup>C</sup> 170.2 in **3**, *δ*<sup>C</sup> 170.0 in **4**), and the 1H chemical shift of H-1'' was deshielded to *δ*<sup>H</sup> 4.18, compared to **3** and **4** (*δ*<sup>H</sup> 4.08 in **3**, *δ*<sup>H</sup> 4.07 in **4**). This finding suggested the possibility that the butyl group in **6** was attached to the different carboxyl groups of compounds **3** and **4**. The gross planar structure of **6** was verified by 2D NMR analysis (1H-<sup>1</sup>H COSY; Figure S13 and HMBC; Figure S15) (Figure 2). The key HMBC correlation between 1-OCH<sup>3</sup> and C-1 supported the presence of the methoxy group at C-1, and the location of the butyl group was determined to be C-1' by the key HMBC correlation of H-1'' and C-1' (Figure 2). Finally, the absolute configuration of **6** was assigned using quantum chemical ECD calculations. As shown in Figure 5, the experimental ECD curve of **6** is consistent with the calculated ECD spectrum (green line) of **6***S* (3*S*). Therefore, the chemical structure of **6** was characterized as (*S*)-1-methyl-1'-butyl citrate.

**Figure 5.** Experimental and calculated ECD spectra of **6**. **Figure 5.** Experimental and calculated ECD spectra of **6**. **Figure 5.** Experimental and calculated ECD spectra of **6**.

The structures of the known compounds (Figure 1) were determined to be (*S*) -dimethyl malate (**1**) [26], (*S*)-5,1′-dimethyl citrate (**2**) [21], 1′-butyl-1,5-dimethyl citrate (**5**) [22], The structures of the known compounds (Figure 1) were determined to be (*S*) -dimethyl malate (**1**) [26], (*S*)-5,1′-dimethyl citrate (**2**) [21], 1′-butyl-1,5-dimethyl citrate (**5**) [22], The structures of the known compounds (Figure 1) were determined to be (*S*)-dimethyl malate (**1**) [26], (*S*)-5,10 -dimethyl citrate (**2**) [21], 10 -butyl-1,5-dimethyl citrate (**5**) [22], and butyl quinate (**7**) [27] by spectroscopic methods, including <sup>1</sup>H and <sup>13</sup>C NMR spectra, by

comparing their spectroscopic data with those previously reported in the literature and MS data obtained from LC/MS analysis. To the best of our knowledge, compounds **2**, **5**, and **7** are reported from *H. rhamnoides* for the first time in this study. by comparing their spectroscopic data with those previously reported in the literature and MS data obtained from LC/MS analysis. To the best of our knowledge, compounds **2**, **5**, and **7** are reported from *H. rhamnoides* for the first time in this study. In the meantime, it appeared likely that compounds **3–7**, including the new com-

and butyl quinate (**7**) [27] by spectroscopic methods, including 1H and 13C NMR spectra,

In the meantime, it appeared likely that compounds **3–7**, including the new compounds, arose from the addition of a butyl group to the corresponding carboxylic acid during the experimental procedure, because *n*-BuOH was used as the organic solvent for solvent partitioning. In order to verify whether compounds **3–7** were genuine natural compounds or artifacts, the aqueous extract of *H. rhamnoides* fruits and fractions derived from the extract were subjected to LC/MS analysis, both alone and co-injected with each compound. As a result, there was no peak with molecular ions corresponding to compounds **3–7** in the crude aqueous extract, whereas all the compounds were detected in the BuOH-soluble fraction, suggesting that the compounds with the butyl group may be artifacts produced by the addition of a butyl group during solvent partition with *n*-BuOH. pounds, arose from the addition of a butyl group to the corresponding carboxylic acid during the experimental procedure, because *n*-BuOH was used as the organic solvent for solvent partitioning. In order to verify whether compounds **3–7** were genuine natural compounds or artifacts, the aqueous extract of *H. rhamnoides* fruits and fractions derived from the extract were subjected to LC/MS analysis, both alone and co-injected with each compound. As a result, there was no peak with molecular ions corresponding to compounds **3–7** in the crude aqueous extract, whereas all the compounds were detected in the BuOH-soluble fraction, suggesting that the compounds with the butyl group may be artifacts produced by the addition of a butyl group during solvent partition with *n*-BuOH.

#### *2.3. Evaluation of Biological Activities of Compounds 1–7 2.3. Evaluation of Biological Activities of Compounds 1–7*

*Plants* **2021**, *10*, x FOR PEER REVIEW 6 of 11

Mesenchymal stem cells (MSCs) in the bone marrow are pluripotent cells that are known to differentiate into osteocytes as well as adipocytes. Since microenvironmental changes cause alterations in the regulation of gene expression in MSC differentiation, where alterations in the expression of the related genes might disturb the balance between osteoprogenitor and adipocyte progenitor cells in osteoporosis patients [28–30], a therapy that is able to adjust the gene expression in MSCs would be promising in the management of postmenopausal osteoporosis. To determine the therapeutic effects of the isolated compounds **1–7** in promoting osteogenesis, all the compounds were examined for their effects on the differentiation of murine MSCs into osteoblasts. The murine mesenchymal stem cell line C3H10T1/2 was treated with 20 µM of the compounds during osteogenesis, and the differentiated cells were stained to indicate alkaline phosphatase (ALP) production, which is considered a distinctive marker of osteoblast differentiation [31]. As a result, compounds **1–7** slightly stimulated differentiation of mesenchymal stem cells into osteocytes (Figure 6). Although all the compounds did not show the superior reaction of ALP staining as in the oryzativol A-treated one, the positive control group with compounds **1–7** showed more than 1.2 times the bone differentiation-promoting effects compared to the negative control group. In particular, compounds **1**, **3**, **4**, **6**, and **7** clearly stimulated the differentiation of osteogenic premature cells up to 1.4 times compared to the negative control group (Figure 6). Mesenchymal stem cells (MSCs) in the bone marrow are pluripotent cells that are known to differentiate into osteocytes as well as adipocytes. Since microenvironmental changes cause alterations in the regulation of gene expression in MSC differentiation, where alterations in the expression of the related genes might disturb the balance between osteoprogenitor and adipocyte progenitor cells in osteoporosis patients [28–30], a therapy that is able to adjust the gene expression in MSCs would be promising in the management of postmenopausal osteoporosis. To determine the therapeutic effects of the isolated compounds **1–7** in promoting osteogenesis, all the compounds were examined for their effects on the differentiation of murine MSCs into osteoblasts. The murine mesenchymal stem cell line C3H10T1/2 was treated with 20 µM of the compounds during osteogenesis, and the differentiated cells were stained to indicate alkaline phosphatase (ALP) production, which is considered a distinctive marker of osteoblast differentiation [31]. As a result, compounds **1–7** slightly stimulated differentiation of mesenchymal stem cells into osteocytes (Figure 6). Although all the compounds did not show the superior reaction of ALP staining as in the oryzativol A-treated one, the positive control group with compounds **1– 7** showed more than 1.2 times the bone differentiation-promoting effects compared to the negative control group. In particular, compounds **1**, **3**, **4**, **6**, and **7** clearly stimulated the differentiation of osteogenic premature cells up to 1.4 times compared to the negative control group (Figure 6).

**Figure 6.** The effects of compounds **1–7** on the differentiation of mesenchymal stem cells (MSCs) toward osteoblasts. (**A**) Stimulatory effect of compounds **1–7** on osteogenic differentiation of MSC. Fully-differentiated C3H10T1/2 cells were stained with alkaline phosphatase (ALP) at 9 d after osteogenic differentiation with 20 µM of compounds **1–7**. (**B**) ALP enzyme activity was measured in osteogenically differentiated C3H10T1/2 cells treated with compounds **1–7**. The values were calculated relatively by setting the untreated negative control to 1. Oryzativol A (5 μM) was added to the experimental set as a positive control and marked as OryA. \* denotes *p* < 0.05 and \*\*\* denotes *p* < 0.001. **Figure 6.** The effects of compounds **1–7** on the differentiation of mesenchymal stem cells (MSCs) toward osteoblasts. (**A**) Stimulatory effect of compounds **1–7** on osteogenic differentiation of MSC. Fully-differentiated C3H10T1/2 cells were stained with alkaline phosphatase (ALP) at 9 d after osteogenic differentiation with 20 µM of compounds **1–7**. (**B**) ALP enzyme activity was measured in osteogenically differentiated C3H10T1/2 cells treated with compounds **1–7**. The values were calculated relatively by setting the untreated negative control to 1. Oryzativol A (5 µM) was added to the experimental set as a positive control and marked as OryA. \* denotes *p* < 0.05 and \*\*\* denotes *p* < 0.001.

#### **3. Materials and Methods**

#### *3.1. General Experimental Procedures*

Optical rotations were measured using a JASCO P-2000 polarimeter (JASCO, Easton, MD, USA). Ultraviolet (UV) spectra were acquired on an Agilent 8453 UV-visible spectrophotometer (Agilent Technologies, Santa Clara, CA, USA). Electronic circular dichroism (ECD) spectra were measured using a Jasco J-1500 spectropolarimeter (Jasco). Nuclear magnetic resonance (NMR) spectra were recorded using a Bruker AVANCE III (Bruker). HRESIMS spectra were recorded on an Agilent 1290 Infinity II series with a 6545 LC/Q-TOF mass spectrometer (Agilent Technologies). Medium-pressure liquid chromatography (MPLC) was performed on a Smart Flash AKROS (Yamazen, Osaka, Japan) using an analytical Universal ODS-SM 120 Å column (3.0 × 20.0 cm, 50 µm) (Yamazen) and Universal Premium Silica gel 60 Å column (2.3 × 12.3 cm, 30 µm) (Yamazen). Preparative and semi-preparative HPLC were performed on a Waters 1525 binary HPLC pump with a Waters 996 photodiode array detector (Waters Corporation, Milford, CT, USA) using an analytical Agilent Eclipse XDB-C18 column (250 × 21.2 mm, 7 µm) and a Phenomenex Luna Phenyl-hexyl 100 Å column (250 × 10 mm, 10 µm), respectively. LC/MS analysis was performed on an Agilent 1200 series HPLC system with a diode array detector and a 6130 Series ESI mass spectrometer using an analytical Kinetex C<sup>18</sup> 100 Å column (100 mm × 2.1 mm i.d., 5 µm; flow rate: 0.3 mL/min) (Phenomenex). All HRESIMS data were obtained using an Agilent G6545B quadrupole time-of-flight (Q-TOF) mass spectrometer (Agilent Technologies) with an Agilent EclipsePlus C<sup>18</sup> column (2.1 mm × 50 mm i.d., 1.8 µm; flow rate: 0.3 mL/min) maintained at 20 ◦C. Silica gel 60 (230–400 mesh; Merck, Darmstadt, Germany), RP-C<sup>18</sup> silica gel (Merck, 230–400 mesh), silica Sep-Pak Vac 6 cc cartridges (Waters), and Diaion HP-20 (Mitsubishi Chemical Co. Ltd., Tokyo, Japan) were used for column chromatography. Sephadex LH-20 (Pharmacia, Uppsala, Sweden) was used as the packing material for molecular sieve column chromatography. Thin-layer chromatography (TLC) was performed using precoated silica gel F<sup>254</sup> plates and RP-C<sup>18</sup> F254s plates (Merck), and spots were detected under UV light or by heating after spraying with anisaldehyde-sulfuric acid.

#### *3.2. Plant Material*

Aqueous extract powder of *H. rhamnoides* fruits was purchased from Korea Beauty and Healthcare Co., Ltd. in October 2018. The material was authenticated by one of the authors (K. H. K.). A voucher specimen of the material (VT-2018) was deposited in the herbarium of the School of Pharmacy, Sungkyunkwan University, Suwon, Korea.

#### *3.3. Extraction and Isolation*

Aqueous extract powder (270 g) of *H. rhamnoides* fruits was suspended in distilled water (700 mL) and then sequentially partitioned with hexane (HX), dichloromethane (MC), ethyl acetate (EA), and *n*-butanol (BuOH). Four layers with different polarities were obtained: HX-soluble (2.7 g), MC-soluble (3.4 g), EA-soluble (7.8 g), and BuOHsoluble (25.2 g) fractions. The BuOH-soluble fraction (25.2 g) was subjected to opencolumn chromatography on a Diaion HP-20 column and eluted stepwise with distilled water (3 L) and 100% methanol (MeOH). The respective fractions were evaporated to dryness. The resulting MeOH fraction (9.7 g) was chromatographed on silica gel with a gradient solvent system of MC/MeOH (10:1 → 1:1), yielding five subfractions (A–E). Subfraction B (1.2 g) was subjected to an RP-C<sup>18</sup> column with a gradient solvent system of MeOH/H2O (15:85 → 100:0), which afforded six subfractions (B1–B6). Subfraction B2 (393.3 mg) was separated by MPLC on a Yamazen Universal Premium silica gel column with MC/MeOH (99:1 → 95:5) to obtain two subfractions (B21 and B22). Subfraction B21 (219.0 mg) was purified by preparative HPLC (MeOH/H2O, 30:70 → 80:20 → 100:0) to obtain four subfractions (B211–B214). Subfraction B212 (58.5 mg) was purified using semi-preparative HPLC (MeCN/H2O, 5:95) to furnish compounds **1** (*t*<sup>R</sup> 29.4 min, 2.2 mg) and **2** (*t*<sup>R</sup> = 43.2 min, 7.4 mg). Subfraction B4 (356.1 mg) was fractionated by preparative

HPLC (MeOH/H2O, 70:30 → 100:0) to yield six subfractions (B41–B46). Subfraction B45 (39.3 mg) was purified using semi-preparative HPLC (MeOH/H2O, 48:52), and compound **6** (*t*<sup>R</sup> = 43.1 min, 2.8 mg) was obtained. Subfraction B46 (58.0 mg) was purified using semipreparative HPLC (MeCN/H2O, 25:75), and compounds **4** (*t*<sup>R</sup> = 35.8 min, 3.1 mg) and **5** (*t*<sup>R</sup> 71.5 min, 1.8 mg) were obtained. Subfraction C (1.2 g) was fractionated by MPLC using an ODS-SM column with MeOH/ H2O (30:70 → 100:0) to yield five subfractions (C1–C5). Subfraction C2 (338.8 mg) was chromatographed using silica Sep-Pak with an isocratic solvent system of MC/MeOH (5:1) to give four subfractions (C21–C24). Subfraction C21 (21.9 mg) was further purified using semi-preparative HPLC (MeOH/H2O, 35:65) to yield compound **7** (*t*<sup>R</sup> = 36.9 min, 7.5 mg). Subfraction D (1.4 g) was separated by preparative HPLC (MeOH/H2O, 25:75 → 90:10 → 100:0), and two subfractions were obtained (D1 and D2). Subfraction D2 (446.2 mg) was fractionated by preparative HPLC again with a gradient solvent system of MeOH/H2O (70:30 → 85:15), which afforded three subfractions (D21–D23). Subfraction D23 (140.1 mg) was purified using semi-preparative HPLC (MeOH/H2O, 40:60 → 50:50) to yield compound **3** (*t*<sup>R</sup> = 50.8 min, 5.7 mg).

#### 3.3.1. (*S*)-1-Butyl-5-Methyl Citrate (**3**)

Yellowish gum; [α] 25 D +6.0 (*c* 0.50, MeOH); UV (MeOH) λmax (log ε) 200 (3.9), 220 (1.5) nm; ECD (MeOH, 0.89% *w*/*v*) λmax (∆ε) 217 (+1.58) nm; <sup>1</sup>H (850 MHz) and <sup>13</sup>C (212.5 MHz) NMR data, see Table 1; ESIMS (positive-ion mode) *m*/*z* 263.1 [M + H]<sup>+</sup> ; HRESIMS (positive-ion mode) *m*/*z* 263.1122 [M + H]<sup>+</sup> (calcd for C11H19O7, 263.1131).

#### 3.3.2. (*S*)-1-Butyl-10 -Methyl Citrate (**4**)

White amorphous powder; [α] 25 <sup>D</sup> −5.4 (*c* 1.65, MeOH); UV (MeOH) λmax (log ε) 201 (3.8), 220 (1.4) nm; ECD (MeOH, 0.94% *w*/*v*) λmax (∆ε) 222 (+0.72) nm; <sup>1</sup>H (850 MHz) and <sup>13</sup>C (212.5 MHz) NMR data, see Table 1; ESIMS (positive-ion mode) *m*/*z* 263.1 [M + H]<sup>+</sup> ; HRESIMS (positive-ion mode) *m*/*z* 263.1119 [M + H]<sup>+</sup> (calcd for C11H19O7, 263.1131).

#### 3.3.3. (*S*)-1-Methyl-10 -Butyl Citrate (**6**)

White amorphous powder; [α] 25 <sup>D</sup> −24.3 (*c* 1.41, MeOH); UV (MeOH) λmax (log ε) 198 (3.9), 218 (1.5) nm; ECD (MeOH, 0.92% *w*/*v*) λmax (∆ε) 217 (+1.71) nm; <sup>1</sup>H (850 MHz) and <sup>13</sup>C (212.5 MHz) NMR data, see Table 1; ESIMS (positive-ion mode) *m*/*z* 263.1 [M + H]<sup>+</sup> ; HRESIMS (positive-ion mode) *m*/*z* 263.1145 [M + H]<sup>+</sup> (calcd for C11H19O7, 263.1131).

#### *3.4. Computational Analysis*

All conformers proposed in this study were acquired through the MacroModel (version 2019-3, Schrödinger LLC) module with mixed torsional/low-mode sampling, implemented with the MMFF94 force field. All searches were initially set in the gas phase, with a 10 kJ/mol energy window limit and a maximum of 10,000 steps to thoroughly explore all potential conformers. The Polak–Ribiere conjugate gradient protocol was established with 10,000 maximum iterations and a 0.001 kJ (mol Å)−<sup>1</sup> convergence threshold on the rms gradient to minimize conformers [23–25]. Conformers proposed in this study (within 5 kJ/mol found in the MMFF force field) were selected for geometry optimization by Tmolex 4.3.1, with the DFT settings of B3-LYP/6-31+G(d,p). ECD calculations for **3***R***/3***S***, 4***R***/4***S,* and **6***R***/6***S* conformers were performed at identical theoretical levels and basis sets. The calculated ECD spectra were simulated by superimposing each transition, where σ is the width of the band at a height of 1/e. ∆*E*<sup>i</sup> and *R*<sup>i</sup> are the excitation energies and rotatory strengths for transition *i*, respectively. In this study, the value of *σ* was 0.10 eV. The excitation energies and rotational strengths for ECD spectra were calculated based on the Boltzmann populations of the conformers, and ECD visualization was carried out using SigmaPlot 14.0.

$$\Delta\epsilon(E) = \frac{1}{2.297 \times 10^{-39}} \frac{1}{\sqrt{2\pi\sigma}} \sum\_{A}^{i} \Delta E\_i R\_i e^{\left[-\left(E - \Delta E\_i\right)^2/(2\sigma)^2\right]}$$

#### *3.5. Cell Culture and Differentiation*

Mouse mesenchymal stem cell line C3H10T1/2 cells were cultured in Dulbecco's modified Eagle's medium (DMEM), 10% fetal bovine serum (FBS), and 1% penicillin-streptomycin (P/S) at 37 ◦C in a 5% CO<sup>2</sup> incubator. For osteogenic differentiation, C3H10T1/2 cells were exposed to DMEM (5% FBS, 1% P/S) media with 10 mM β-glycerophosphate and 50 µg/mL ascorbic acid for 9 d. During osteogenic differentiation, 20 µM of compounds **1–7** were added to the cells, and 5 µM of oryzativol A was used as the positive control.

To investigate the degree of differentiation into osteocytes, cells were stained for alkaline phosphatase (ALP) after 9 days of osteogenic differentiation. After fixation with 4% formaldehyde solution, cells were exposed to 0.4 mg/mL of nitroblue tetrazolium and 0.2 mg/mL of 5-bromo-4chloro-30 -indolyphosphate in alkaline phosphatase (AP) buffer (100 mM Tris-HCl, pH 9.5; 100 mM NaCl, 5 mM MgCl2) in the dark for 4 h. Next, the reaction was stopped with 5 mM EDTA. To measure APL activity, cell lysates were collected from the differentiated osteogenic cells and treated with *p*-nitrophenyl phosphate (*p*-NPP) solution (Alkaline Phosphatase Assay Kit (ab83369; Abcam, Cambridge, MA, USA). Color changes of the samples were measured at 405 nm with a SpectraMax M2 Microplate Spectrophotometer.

#### *3.6. Statistical Analysis*

Each sample was tested in triplicate, and the test was repeated three times. Data are presented as the mean ± standard deviation (S.D.). Differences between the control and experimental groups were analyzed using a two-tailed unpaired Student's *t*-test*;* statistical significance was defined as *p* < 0.05.

#### **4. Conclusions**

In this study, phytochemical examination of the extracts of *H. rhamnoides* fruits resulted in the isolation of seven organic acid derivatives: one malate derivative (**1**), five citrate derivatives (**2–6**), and one quinate derivative (**7**). The structure of the isolates was elucidated by analysis of 1D and 2D NMR data and HR-ESIMS data as well as quantum chemical ECD calculations. Three of the citrate derivatives were identified as new compounds: (*S*)-1-butyl-5-methyl citrate (**3**), (*S*)-1-butyl-10 -methyl citrate (**4**), and (*S*)-1 methyl-10 -butyl citrate (**6**), which turned out to be isolation artifacts. To the best of our knowledge, compounds **2**, **5**, and **7** are reported for the first time from *H. rhamnoides* in this study. An osteogenesis-promoting activity bioassay in C3H10T1/2 cells demonstrated that compounds **1–7** showed little more than 1.2 times the bone differentiation promoting effect, compared to the negative control group. In particular, compounds **1**, **3**, **4**, **6**, and **7** clearly stimulated the differentiation of osteogenic premature cells up to 1.4 times, compared to the negative control group. Although the stimulatory effect of the active compounds on osteogenic differentiation was far behind that of the positive control, oryzativol A, compounds **1**, **3**, **4**, **6**, and **7** apparently helped to differentiate the mesenchymal stem cells into osteocytes.

**Supplementary Materials:** 1D and 2D NMR, and HR-ESIMS of compounds **3**, **4**, and **6** are available online at https://www.mdpi.com/article/10.3390/plants10050860/s1, Figure S1: HR-ESIMS data of **3**, Figure S2: <sup>1</sup>H NMR spectrum of **3** (CD3OD, 700 MHz), Figure S3: <sup>1</sup>H-1H COSY spectrum of **3** (CD3OD), Figure S4: HSQC spectrum of **3** (CD3OD), Figure S5: HMBC spectrum of **3** (CD3OD), Figure S6: HR-ESIMS data of **4**, Figure S7: <sup>1</sup>H NMR spectrum of **4** (CD3OD, 700 MHz), Figure S8: <sup>1</sup>H-1H COSY spectrum of **4** (CD3OD), Figure S9: HSQC spectrum of **4** (CD3OD), Figure S10: HMBC spectrum of **4** (CD3OD), Figure S11: HR-ESIMS data of **6**, Figure S12: <sup>1</sup>H NMR spectrum of **6** (CD3OD, 700 MHz), Figure S13: <sup>1</sup>H-1H COSY spectrum of **6** (CD3OD), Figure S14: HSQC spectrum of **6** (CD3OD), Figure S15: HMBC spectrum of **6** (CD3OD), Table S1: Gibbs free energies and Boltzmann distribution of conformers **3***S* (3*S*), Table S2: Gibbs free energies and Boltzmann distribution of conformers **4***S* (3*S*), Table S3: Gibbs free energies and Boltzmann distribution of conformers **6***S* (3*S*). **Author Contributions:** Conceptualization, T.S.J. and K.H.K.; formal analysis, Y.H.L., H.J.J., K.H.P., S.-H.K., J.K.K., and J.-C.K.; investigation, Y.H.L., H.J.J., and K.H.P.; writing—original draft preparation, Y.H.L., H.J.J., S.-H.K., and K.H.K.; writing—review and editing, J.K.K. and K.H.K.; visualization, H.J.J. and S.-H.K.; supervision, T.S.J. and K.H.K.; project administration, T.S.J. and K.H.K.; funding acquisition, T.S.J. and K.H.K. All authors have read and agreed to the published version of the manuscript.

**Funding:** This work was supported by a grant from the National Research Foundation of Korea (NRF), funded by the Korean government (MSIT) (grant number: 2019R1A5A2027340) and a Korea Institute of Science and Technology intramural research grant (2E30641). This research was also supported by the Ministry of Trade, Industry and Energy (MOTIE), Korea Institute for Advancement of Technology (KIAT) through the Program for Smart Specialization Infrastructure Construction.

**Conflicts of Interest:** The authors declare no conflict of interest.

#### **References**

