DOI QR코드

DOI QR Code

Liquid chromatography-tandem mass spectrometric analysis of oleracone D and its application to pharmacokinetic study in mice

  • Lim, Dong Yu (College of Pharmacy, Dankook University) ;
  • Lee, Tae Yeon (College of Pharmacy, Dankook University) ;
  • Lee, Jaehyeok (BK21 FOUR Community-Based Intelligent Novel Drug Discovery Education Unit, Vessel-Organ Interaction Research Center (VOICE), College of Pharmacy and Research Institute of Pharmaceutical Sciences, Kyungpook National University) ;
  • Song, Im-Sook (BK21 FOUR Community-Based Intelligent Novel Drug Discovery Education Unit, Vessel-Organ Interaction Research Center (VOICE), College of Pharmacy and Research Institute of Pharmaceutical Sciences, Kyungpook National University) ;
  • Han, Young Taek (College of Pharmacy, Dankook University) ;
  • Choi, Min-Koo (College of Pharmacy, Dankook University)
  • Received : 2021.05.24
  • Accepted : 2021.07.26
  • Published : 2021.10.25

Abstract

We have demonstrated a sensitive analytical method of measuring oleracone D in mouse plasma using a liquid chromatography-tandem mass spectrometry (LC-MS/MS). Oleracone D and oleracone F (internal standard) in mouse plasma samples were processed using a liquid-liquid extraction method with methyl tertbutyl ether, resulting in high and reproducible extraction recovery (80.19-82.49 %). No interfering peaks around the peak elution time of oleracone D and oleracone F were observed. The standard calibration curves for oleracone D ranged from 0.5 to 100 ng/mL and were linear with r2 of 0.992. The inter- and intra-day accuracy and precision and the stability fell within the acceptance criteria. The pharmacokinetics of oleracone D following intravenous and oral administration of oleracone D at doses of 5 mg/kg and 30 mg/kg, respectively, were investigated. When oleracone D was intravenously injected, it had first-order elimination kinetics with high clearance and volume of distribution values. The absolute oral bioavailability of this compound was calculated as 0.95 %, with multi-exponential kinetics. The low aqueous solubility and a high oral dose of oleracone D may explain the different elimination kinetics of oleracone D between intravenous and oral administration. Collectively, this newly developed sensitive LC-MS/MS method of oleracone D could be successfully utilized for investigating the pharmacokinetic properties of this compound and could be used in future studies for the lead optimization and biopharmaceutic investigation of oleracone D.

Keywords

1. Introduction

Portulaca oleracea L. is a well-known annual plant found distributed in temperate, tropical and subtropical regions.1 It has been used for some time in many countries for the alleviation of various diseases due to its pharmacological effects that include anti- inflammatory, 2-4 antioxidative, 5, 6 anti-tumor, 7 antiaging, and hypolipidemic effects.1, 5 These effects may be due to the flavonoid constituents of Portulaca oleracea L.8-11 Yang et al.9, 10 were able to isolate a novel series of flavonoids, named oleracone C, D, E, and F (Fig. 1), from Portulaca oleracea L. and investigated their activities. Comparing the antioxidant activity of these oleracones by evaluating their radical scavenging effect, oleracone D was the most effective among the tested oleracones.5, 9, 10 Oleracone D also had the highest anticholinesterase activity; however, oleracone C and E also had anticholinesterase activity with lower inhibitory efficacy.10 Oleracone E and F had a similar lifespan extension potential from the lifespan assay using age-synchronized nematodes.5 The results demonstrated that these oleracone derivatives have potent antioxidative activity and may exhibit a therapeutic potential on the aging process. For further assays for investigating its activity and toxicity, the development of an analytical method for monitoring oleracone D concentration in a biological matrix is a necessary progression.

BGHHBN_2021_v34n5_193_f0001.png 이미지

Fig. 1. Chemical structure of oleracone D and oleracone F (internal standard, IS).

Pharmacokinetics of oleracone C in rats and the analytical detection method for plasma oleracone C concentrations have previously been reported.2 The plasma concentration of oleracone C had an elimination half-life (T1/2) of 0.3 h and showed low oral bio availability of 8.3%. For the determination of plasma oleracone C, the liquid chromatography tandem mass spectrometry (LC-MS/MS) method has been developed in the standard range of 10-2000 ng/mL, with a volume of 100 μL of rat plasma samples processed using four volumes of methanol precipitation method.2 Although oleracone D had better in vitro antioxidative and neuroprotective effects than those of oleracone C, the analytical method and pharmacokinetic study of oleracone D have not been previously reported.

Therefore, this study aimed to develop and validate a sensitive and reproducible LC-MS/MS method for measuring concentrations of oleracone D and apply this analytical method to a pharmacokinetic study of oleracone D in mice.

To achieve a substantial extraction recovery and negligible matrix effect from a biological matrix such as plasma, the most widely used sample preparation methods are protein precipitation (PPT) and liquid liquid extraction (LLE).12 Since LLE has an advantage in lowering interferences from the sample matrix and increasing the sensitivity of the analyte, 12, 13 LLE method was, therefore, utilized in this study for the sample preparation method for the analysis of oleracone D in 30 μL of mouse plasma samples. Additionally, we aimed to validate the analytical method for its linearity, selectivity, accuracy, precision, stability, recovery, and matrix effects according to the U.S. Food and Drug Administration Guideline for Bioanalytical Method.14

2. Experimental

2.1. Chemicals and reagents

Oleracone D and oleracone F (Fig. 1) were synthesized with a purity of >99.0% using a previously described method by Yoon et al.5 The purity was confirmed using nuclear magnetic resonance spectroscopy and mass spectroscopy and oleracone F was used as the internal standard (IS). Acetonitrile, water, and methanol were purchased from Tedia (Fairfield, CT, USA). Formic acid and methyl tert-butyl ether (MTBE) were purchased from Sigma-Aldrich Chemical Co. (St. Louis, MA, USA). All solvents and chemicals were of HPLC and of reagent grade.

2.2. Preparation of the stock and working solutions

Stock solutions were prepared individually by dissolving oleracone D and oleracone F in methanol at concentrations of 1 mg/mL. The oleracone D working solutions were prepared by serial dilution of the stock solution with acetonitrile to final concentrations of 0.5, 1, 2, 5, 20, 50, 100 ng/mL for the calibration standards and 1.5, 10, and 75 ng/mL for the quality control (QC) samples. The oleracone F solution was prepared at a concentration of 5 ng/mL, with a dilution of the stock solution with water. The stock and working solution of oleracone D were stored at -80 oC during the analysis.

2.3. Preparation of calibration standards and quality control samples

A volume of 30 μL of the working solution for calibration standards and QC samples were evaporated and reconstituted with a 30 μL aliquot of blank mouse plasma. The final concentrations of calibration standards were 0.5, 1, 2, 5, 20, 50, 100 ng/mL and 1.5 (low QC), 10 (middle QC), 75 (high QC) ng/mL, respectively, for the QC samples.

2.4. Sample preparation

The calibration standards and QC samples (30 μL) were combined with 20 μL of oleracone F solution (5 ng/mL in water) and 400 μL of MTBE. The mixture was vigorously vortexed for 15 min and subsequently centrifuged at 16, 000 × g for 5 min. The supernatant was transferred to a clean tube and evaporated to dryness until dry using a gentle stream of nitrogen. The residue was reconstituted in 200 μL of mobile phase, and subsequently, a 10 μL aliquot of the solution was injected into the LC-MS/MS system.

2.5. Instrument conditions

The LC system was an Agilent Infinity 1260 Infinite II HPLC system (Agilent Technologies, Santa Clara, CA, USA), and chromatographic separation was carried out using a Luna C18 (150 × 2.0 mm, 5 μm; Phenomenex, Torrance, CA, USA). The mobile phase was pumped with isocratic elution of water that contained 0.1% formic acid : acetonitrile containing 0.1% formic acid (20:80, v/v).

The flow rate was set to 0.2 mL/min with a column temperature maintained at 30 ºC. The Agilent 6430 triple quadrupole MS was equipped with an electrospray ionization source and was used for mass spectrometric detection and quantitative analysis. The mass spectrometer was operated in the positive ion mode, with multiple reaction monitoring (MRM) transitions at m/z 299.1→205.0 for oleracone D and at m/z 313.0 →219.1 for oleracone F (IS) with optimized fragmentor of 105 or 140 V, and collision energy of 20 eV, respectively.

2.6. Method validations

2.6.1. Selectivity and linearity

To determine the selectivity was determined by a comparison of the chromatogram of the blank plasma from six different mice and corresponding mouse blank plasma containing a lower limit of quantification (LLOQ) samples (0.5 ng/mL of oleracone D) and IS was performed.

The calibration standard curve for oleracone D (0.5-100 ng/mL) was obtained from the ratio of the peak areas of oleracone D and IS and plotted to the concentrations of oleracone D. The results were fitted using a least square linear regression analysis.

2.6.2. Precision and accuracy

Two sets of QC samples (low, middle, and high QC) on five independent days were used to determine the inter-day precision and accuracy. Intra-day precision and accuracy were determined by analyzing six sets of QC samples (low, middle, and high QC) on the same days. The precision was expressed as the coefficient of variance (CV, %), whereas accuracy was expressed as the percentage of the measured QC concentration to the nominal QC concentration.

2.6.3. Matrix effect and recovery

The matrix effect of oleracone D was determined by comparing the peak area obtained from the post extraction blank plasma that was spiked with three QC concentrations of oleracone D and with the peak area from the corresponding concentration of the neat solution. The matrix effect of IS was determined in the same manner as the procedure of oleracone D using a 5 ng/mL concentration of IS.

The extraction recovery of oleracone D was determined by the comparison of the peak area obtained from the pre-extraction samples spiked with three QC concentrations of oleracone D, with that of the post-extraction blank plasma spiked with corresponding QC concentrations of oleracone D. The extraction recovery of IS was determined in the same manner using the procedure of oleracone D using the concentration of 5 ng/mL of IS.

2.6.4. Stability

The stability of oleracone D in the mouse plasma was determined from the three sets of QC samples (low and high QC) under various conditions. The bench top stability was determined by incubating the QC samples at 25 oC for 6 h. Thereafter, the concentrations of the QC samples were determined using freshly prepared calibration standards. For the freeze-thaw stability, QC samples were analyzed after three cycles of freeze–thawing. One cycle of the freeze–thaw process involved storing the QC samples at -80 °C for over 12 h and thawing at 25 °C for 6 h. After undergoing three cycles of the freeze-thaw process, the concentrations of the QC samples were determined using freshly prepared calibration standards. The autosampler stability was evaluated by placing the extracted QC samples in the autosampler at 6 °C for 24 h, and the subsequently determined concentrations were calculated using freshly prepared calibration standards. The long term stability of oleracone D stock solution was evaluated by comparison the oleracone D stock solution stored at -80 °C for 6 months with freshly prepared oleracone D stock solution. Both oleracone D stock solutions were diluted to 500 ng/ mL using mobile phase and compared the peak area of corresponding oleracone D stock solutions.

2.7. Pharmacokinetic study

All experiments involving animals and the utilized protocols were approved by the Animal Care and Use Committee of the Kyungpook National University (No. 2019-0126) and were conducted in accordance with the National Institutes of Health guidelines for the care and use of laboratory animals.15

Male ICR mice (7-8 weeks old, 30-35 g) were purchased from Samtako (Osan, Korea) and acclimated to the animal facility of Kyungpook National University for 1 week with ad libitum access to food and water and fasted for 12 h prior to the pharmacokinetic experiments. After the intravenous administration of oleracone D (5 mg/kg dissolved in 2 mL mixture of DMSO : saline = 20 : 80 [v/v]) via the tail vein, blood samples were collected via the Retro-Orbital plexus using a heparinized collection tube at 0, 0.25, 0.5, 1, 2, 4, 8, and 24 h while mice were under anesthesia with isoflurane for 5 min. After the oral administration of oleracone D (30 mg/kg suspended in 5 mL of 0.5% carboxymethyl cellulose suspension) using oral gavage, blood samples were collected at 0, 0.25, 0.5, 1, 2, 4, 8, and 24 h, as described above. All blood samples were centrifuged at 16, 000 × g for 1 min to separate plasma, and the plasma sample was stored at -80 oC until analysis. The plasma samples (30 μL) were added to 20 μL of IS solution (5 ng/mL oleracone F in water) and 400 μL of MTBE. The mixture was vigorously vortexed for 15 min and then centrifuged at 16, 000 × g for 5 min. The supernatant was transferred to a clean tube and evaporated under dry using a gentle stream of nitrogen. The residue was reconstituted in 200 μL of mobile phase, and subsequently, 10 μL aliquot of the solution was injected into the LC-MS/ MS system.

2.8. Plasma protein binding

Protein binding of oleracone D (75 ng/mL, high QC) in mouse plasma was determined using a rapid equilibrium dialysis kit (ThermoFisher Scientific Korea, Seoul, Korea) according to the manufacturer’s protocol.7 Briefly, 100 μL of mouse plasma containing 75 ng/mL oleracone D was added to the inner sample chamber of a semipermeable membrane (molecular weight cut-off 8, 000 Da). To the outer chamber, 300 μL of phosphate buffered saline (PBS) was added. The samples were then incubated for 4 h at 37 °C on a shaking incubator at 300 rpm, followed by 30 μL aliquot being removed from both the sample and buffer chambers. Samples were treated with equal volumes of fresh PBS or blank plasma to match the sample matrices. The matrix-matched samples (60 μL) were mixed with 40 μL of IS solution (5 ng/mL oleracone F in water) and 800 μL of MTBE. The mixture was vigorously vortexed for 15 min and subsequently centrifuged at 16, 000 × g for 5 min. The supernatant was transferred to a clean tube and evaporated under dry using a gentle stream of nitrogen. The residue was reconstituted in 200 μL of mobile phase, and subsequently, 10 μL aliquot of the solution was injected into the LC-MS/MS system.

2.9. Plasma stability

The determination of the plasma stability of oleracone D (75 ng/mL, high QC) in mouse plasma was conducted.6 Briefly, 100 μL of mouse plasma containing 75 ng/mL oleracone D was incubated for 8 h at 37 °C on a shaking incubator at 300 rpm. From this, 30 μL aliquot was taken and added to 20 μL of oleracone F solution (5 ng/mL in water) and 400 μL of MTBE. The mixture was vigorously vortexed for 15 min and then centrifuged at 16, 000 × g for 5 min. The supernatant was transferred to a clean tube and evaporated under dry using a gentle stream of nitrogen. The residue was reconstituted in 200 μL of mobile phase, and subsequently, 10 μL aliquot of the solution was injected into the LC-MS/MS system.

2.10. Data analysis and statistics

Calculation of the pharmacokinetic parameters was conducted using WinNonlin 5.1 software (Pharsight Co., Mountain View, CA, USA). All data are expressed as the mean ± standard deviation (SD).

3. Results and Discussion

3.1. LC-MS/MS analysis of oleracone D

For optimization of the electrospray ionization conditions of oleracone D and IS, each compound was injected directly into the mass spectrometer ionization source. The oleracone D and IS had optimal ionization in the positive mode. The MRM transition of oleracone D was selected from the precursor ion ([M+H]+, m/z 299.1) and the most frequent product ion (m/z 205.0). The MRM transition of IS was selected from the precursor ion ([M+H]+, m/z 313.0) and the most frequent product ion (m/z 219.1) (Fig. 2).

BGHHBN_2021_v34n5_193_f0002.png 이미지

Fig. 2. Product ion scan of oleracone D and IS (oleracone F).

3.2. Selectivity and linearity

Fig. 3 shows the typical chromatograms of double blank, zero blank, LLOQ sample (0.5 ng/mL), and plasma sample at 8 h after intravenous administration of oleracone D. The retention times for oleracone D and IS were 3.4 and 2.7 min, respectively. The signal- to-noise (S/N) ratio of oleracone D was > 10.0 in the LLOQ samples, with no significant matrix interference at the retention times of oleracone D and IS in the blank samples, compared to the LLOQ samples, when observing the chromatograms (Fig. 3).

BGHHBN_2021_v34n5_193_f0003.png 이미지

Fig. 3. Representative multiple reaction monitoring (MRM) chromatograms of oleracone D and IS (oleracone F) in mouse plasma of double blank, zero blank, lower limit of quantification (LLOQ) sample (0.5 ng/mL), and plasma sample at 8 h following intravenous injection of oleracone D.

A calibration standard curve was linear within the oleracone D concentration ranges of 0.5-100 ng/mL, and a coefficient of determination (r2 ) was 0.992 determined from linear regression analysis, with a weighting of 1/concentration2 (Table 1).

Table 1. Calculated concentrations of oleracone D in calibration standards

BGHHBN_2021_v34n5_193_t0001.png 이미지

Data represented as mean ± SD (n = 6).

3.3. Precision and accuracy

Table 2 summarizes the intra- and inter-day precision and accuracy for oleracone D from three levels of QC samples. The intra- and inter-day precision was from 4.24 to 8.43% for oleracone D and the intra-and inter-day accuracy was from 89.38 to 108.43%, which fell in the acceptable criteria (less than 15%).

Table 2. Intra- and inter-day precision and accuracy of oleracone D

BGHHBN_2021_v34n5_193_t0002.png 이미지

Data represented as mean ± SD.

3.4. Matrix effect and recovery

The results of the extraction recoveries and matrix effects are summarized in Table 3. The extraction recoveries for oleracone D were calculated at three levels of QC samples and were high and reproducible in the range of extraction recovery from 80.19 to 82.49% and CV from 3.91 to 9.46%. Therefore, the LLE method utilized in this study was capable of efficiently extracting oleracone D from the mouse plasma. The matrix effects ranged from 89.21 to 106.23%, suggesting that the coeluting substances did not interfere with the ionization of oleracone D. The low CV in the matrix effect indicated that there were no significant differences for the peak areas of oleracone D in the three concentrations of the three samples of each plasma matrix. Therefore, it was possible to exclude any matrix effect on ion suppression or enhancement. The extraction recovery and matrix effects of IS was also high and reproducible.

Table 3. Extraction recoveries and matrix effects of oleracone D and IS

BGHHBN_2021_v34n5_193_t0003.png 이미지

Data represented as mean ± SD (n=3).

3.5. Stability

In Table 4, the results of stability experiments are presented. The accuracy of oleracone D ranged from 92.46% to 107.73% and the precision ranged from 2.33% to 5.96% for the three different conditions of stability test. As there was no significant difference observed for each measured concentration compared with nominal concentration under various stability treatment conditions, these results provided evidence that the oleracone D in the mouse plasma was stable up to 6 h at 25 ºC for bench-top stability, and was stable for 24 h in the autosampler stability after sample treatment, and also remained stable over three freeze-thaw cycles. The peak area for oleracone D stock solution stored at -80 °C for 6 months was 5025.67±170.12 (n = 3) and that of freshly prepared oleracone D stock solution was 5447.22±96.92 (n = 3). The mean peak area ratio was calculated as 92.26%, suggesting that the oleracone D stock solution was stable for storing at -80 oC for more than 6 months. 16 However, we should note that the long term plasma stability of oleracone D need to be evaluated for further preclinical studies using oleracone D.

Table 4. Stability of oleracone D

BGHHBN_2021_v34n5_193_t0004.png 이미지

Data represented as mean ± SD (n=3).

3.6. Pharmacokinetics of oleracone D

Temporal profiles for the plasma concentration of oleracone D after intravenous and oral administration are shown in Fig. 4, and the relevant pharmacokinetic parameters are detailed in Table 5. The plasma concentration of oleracone D after intravenous injection demonstrated first-order elimination kinetics, with an elimination half-life of 1.9±0.3 h. However, oleracone D had high clearance (CL) and volume of distribution (Vd) in mice, resulting in low plasma exposure, even with intravenous administration of oleracone D. In addition, the results suggested that oleracone D may undergo substantial metabolism or distribution even though the underlying mechanism requires further investigation. As oleracone D was stable for 8 h in the plasma and had very high plasma protein binding (99.6 ± 0.2%) (Table 6), the plasma instability could be excluded from the contributing factor for high CL and Vd of oleracone D in mice.

BGHHBN_2021_v34n5_193_f0004.png 이미지

Fig. 4. Plasma concentration-time profile of oleracone D in mice following an intravenous injection of oleracone D at a dose of 5 mg/kg (●) and an oral administration of oleracone D at a dose of 30 mg/kg (○). Data are expressed as mean ± SD (n=3).

Table 5. Pharmacokinetic parameters of oleracone D in mice following intravenous injection or oral administration of oleracone D

BGHHBN_2021_v34n5_193_t0005.png 이미지

C0: initial plasma concentration; AUC8h or AUC: Area under plasma concentration-time curve from zero to 8h or infinity; T1/2: elimination half-life; MRT: mean residence time; CL: clearance (Dose/AUC); Vd: volume of distribution (MRT×CL) BA: Bioavailability (Dose normalized AUCPO/dose normalized AUCIV×100) Data are expressed as mean ± SD (n=3).

Table 6. Physicochemical properties of oleracone D

BGHHBN_2021_v34n5_193_t0006.png 이미지

aLogP value was obtained from Chemdraw data base

The plasma concentration of oleracone D following oral administration reached Cmax within the initial sampling times (i.e. Tmax of 0.25 h), providing evidence that the gastrointestinal absorption of oleracone D is rapid, but was eliminated by multiexponential kinetics. Considering its rapid intestinal absorption, it had better collect plasma samples at earlier time points to obtain the accurate kinetic parameters for Cmax and Tmax. However, very low levels of the plasma concentrations of oleracone D after oral administration maintained even with high doses of oleracone D within the mice (30 mg/kg), resulted in low oral bioavailability of 0.95%.

We next investigated the different elimination kinetics of intravenous and oral administration, by measuring the aqueous solubility of oleracone D (Table 6). Its solubility (84.5±0.4 mg/mL) was much lower than the oral suspension (6 mg/mL). A dissolved fraction of oleracone D could be absorbed quickly due to its lipophilicity (LogP of 3.4); however, undissolved oleracone D in the intestine could be absorbed steadily when dissolved in the gastric fluid, which possibly increased the half-life of oleracone D compared with intravenous administration.

In agreement with our findings, Yang et al.2 reported that plasma concentrations of oleracone C, a structural analog of oleracone D, eliminated quickly, with a half-life of 0.35 h, and also showed the low oral bioavailability (8.3%) in rats following intravenous and oral administration (3 mg/mL).

4. Conclusions

We have developed and validated a sensitive analytical method for oleracone D using an LC-MS/ MS system through inter-day and intra-day precision and accuracy, matrix effect, extraction recovery, and stability tests. We have successfully investigated the pharmacokinetic characteristics and oral bioavailability of oleracone D in mice following its intravenous and oral administration using our validated and sensitive analytical method. Although the biological activities of oleracone D on the antioxidative and antiaging effects, it showed low oral bioavailability, high CL, and Vd values, which was partly attributed to its low aqueous solubility.

Acknowledgements

This work has no external grant.

References

  1. Y. X. Zhou, H. L. Xin, K. Rahman, S. J. Wang, C. Peng and H. Zhang, BioMed Research International, 2015, 925631 (2015).
  2. X. Yang, Z. Ying, F. He, X. Ying and G. Yang, Fitoterapia, 131, 44-49 (2018). https://doi.org/10.1016/j.fitote.2018.10.005
  3. C. Y. Li, Y. H. Meng, Z. M. Ying, N. Xu, D. Hao, M. Z. Gao, W. J. Zhang, L. Xu, Y. C. Gao and X. X. Ying, Journal of Agricultural and Food Chemistry, 64(29), 5837-5844 (2016). https://doi.org/10.1021/acs.jafc.6b02673
  4. Y. Meng, Z. Ying, Z. Xiang, D. Hao, W. Zhang, Y. Zheng, Y. Gao and X. Ying, The Journal of Pharmacy and Pharmacology, 68(3), 397-405 (2016). https://doi.org/10.1111/jphp.12526
  5. J. A. Yoon, C. Lim, D. S. Cha and Y. T. Han, Applied Sciences, 9(19), Article 4014 (2019).
  6. H. Son, K. Noh, C. Kang, M. Na, S. Oh, I. S. Song and W. Kang, Journal of Pharmaceutical and Biomedical Analysis, 166, 291-294 (2019). https://doi.org/10.1016/j.jpba.2019.01.030
  7. M. K. Choi, S. J. Nam, H. Y. Ji, M. J. Park, J. S. Choi and I. S. Song, Pharmaceutics, 12(3), e268 (2020).
  8. E. K. Akkol, I. E. Orhan and E. Yesilada, Food Chemistry, 131(2), 626-631 (2012). https://doi.org/10.1016/j.foodchem.2011.09.041
  9. X. Yang, Z. Ying, H. Liu, X. Ying and G. Yang, Natural Product Research, 33(24), 3500-3506 (2019). https://doi.org/10.1080/14786419.2018.1484465
  10. X. Yang, W. Zhang, X. Ying and D. Stien, Fitoterapia, 127, 257-262 (2018). https://doi.org/10.1016/j.fitote.2018.02.032
  11. Y. Duan, Z. Ying, M. Zhang, X. Ying and G. Yang, Natural Product Research, DOI: 10.1080/14786419.2020.1815742. Online ahead of print (2020).
  12. J. H. Jeon, J. Lee, M. K. Choi and I. S. Song, Archives of Pharmacal Research, 43(12), 1335-1346 (2020). https://doi.org/10.1007/s12272-020-01289-0
  13. M. K. Choi, S. Jin, J. H. Jeon, W. Y. Kang, S. J. Seong, Y. R. Yoon, Y. H. Han and I. S. Song, Journal of Ginseng Research, 44(2), 229-237 (2020). https://doi.org/10.1016/j.jgr.2018.10.006
  14. D. Zimmer, Bioanalysis, 6(1), 13-19 (2014). https://doi.org/10.4155/bio.13.298
  15. I. S. Song, S. J. Nam, J. H. Jeon, S. J. Park and M. K. Choi, Pharmaceutics, 13(5), e628 (2021).
  16. N. van de Merbel, N. Savoie, M. Yadav, Y. Ohtsu, J. White, M.F. Riccio, K. Dong, R. de Vries and J. Diancin, The AAPS Journal, 16(3), 392-399 (2014). https://doi.org/10.1208/s12248-014-9573-z