DOI QR코드

DOI QR Code

Rapid determination and quantification of hair-growth compounds in adulterated products by ultra HPLC coupled to quadrupole-orbitrap MS

  • Lee, Ji Hyun (Division of Advanced Analysis, Toxicological Evaluation and Research Department, National Institute of Food and Drug Safety Evaluation, Ministry of Food and Drug Safety) ;
  • Park, Han Na (Division of Advanced Analysis, Toxicological Evaluation and Research Department, National Institute of Food and Drug Safety Evaluation, Ministry of Food and Drug Safety) ;
  • Kang, Gihaeng (Division of Advanced Analysis, Toxicological Evaluation and Research Department, National Institute of Food and Drug Safety Evaluation, Ministry of Food and Drug Safety) ;
  • Kim, Nam Sook (Division of Advanced Analysis, Toxicological Evaluation and Research Department, National Institute of Food and Drug Safety Evaluation, Ministry of Food and Drug Safety) ;
  • Park, Seongsoo (Division of Advanced Analysis, Toxicological Evaluation and Research Department, National Institute of Food and Drug Safety Evaluation, Ministry of Food and Drug Safety) ;
  • Lee, Jongkook (College of Pharmacy, Kangwon National University) ;
  • Kang, Hoil (Division of Advanced Analysis, Toxicological Evaluation and Research Department, National Institute of Food and Drug Safety Evaluation, Ministry of Food and Drug Safety)
  • Received : 2019.01.27
  • Accepted : 2019.03.29
  • Published : 2019.04.25

Abstract

Recently, a number of adulterated products, which are advertised as hair-growth enhancer have been emerged among those who suffer hair loss disease. For continuous control of illegal products, in this study, a rapid and sensitive method for simultaneous screening of 12 compounds that enhance hair-growth was established to protect public health by ultrahigh-performance liquid chromatography coupled to quadrupole-orbitrap mass spectrometry (UHPLC-Q-Orbitrap-MS). Fragmentation pathways of them were proposed based on $MS^2$ spectral data obtained using the established method. In this analysis, the LODs and LOQs ranged from 0.05 to 50 ng/mL and from 0.17 to 167 ng/mL, respectively. The square of the linear correlation coefficient ($R^2$) was determined as more than 0.995. The intra- and inter-assay accuracies were respective 88-112 % and 88-115 %. Their precision values were measured within 5 % (intra-day) and 10 % (inter-day). Mean recoveries of target compounds in adulterated products ranged from 84 to 115%. The relative standard deviation of stability was less than 12 % at $4^{\circ}C$ for 48 h. The method was employed to screen 14 dietary supplements advertised to be effective for the treatment of hair loss. Some of the products (~21 %) were proven to contain synthetic drugs that promote hair growth such as triaminodil, minoxidil, and finasteride.

Keywords

1. Introduction

Hair loss has been one of the serious problems inlife quality. The market of the products for the treatments of hair loss thus grows rapidly. Several anti-alopeciadrugs have been used orally or topically to alleviatehair loss and promote hair growth.1-2 Some of the drugs such as finatsteride and minoxidil have been approved by the United States Food and Drug Administration (US FDA).3-4 It has been believed that 5α-reductase, which converts testosterone into dihydrotestosterone (DHT), plays an important role in hair loss.5 Finasteride, a 5α-reductase inhibitor, impedes scalp hair follicle miniaturization by lowering the blood level of DHT. Minoxidil, a vasodilator, enhances hair growth by providing hair follicles with more oxygen and nutrients.6,7 Several other drugshave also been prescribed for the treatments of hairloss by approval of authorities in some countries and used off-label worldwide. Dutasteride reduces the formation of DHT by inhibiting 5α-reductase.8Spironolactone, cyproterone acetate and flutamidebind androgen receptors in competition with DHT.9-10Not only scalp hairs but facial hairs including beards, moustaches, sideburns, and eyebrows are also considered as important factors for good out wardappearances. Contrary to scalp hair loss, whiskers, moustaches, and sideburns are stimulated to grow by testosterone, and androgenic-anabolic steroids (AAS) has thus been used for their growth.11,12 In addition, the cosmetic formulation of bimatoprost, a synthetic prostagland in analogue (PGA), to lengthen eyelashes was approved by US FDA and in high demand for women.13-14

A great number of products have been advertised and sold for the treatments of hair loss in dietary supplements. Consumers believe that these products are safe and do not contain any synthetic drugs or their analogues. However, several dietary supplement containing synthetic hair loss remedies have been manufactured and distributed illegally via the internetand black market without any toxicological declaration on their labels.15-16 Abuse and misuse of the products could pose a significant risk to public health due todiverse side effects of their adulterants.3,17 Therefore, routine screening of these products should be performed to protect consumers from the risk of illegal adulterants.

In recent years, analytical methods for detection of drugs that reverse hair loss have been developed by using capillary electrophoresis (CE), high performanceliquid chromatography (HPLC), and liquid chroma-tography-mass spectrometry (LC-MS/MS).13,18-23 In our laboratory, LC-MS method for determination of hair-growth compound were developed and applied to screen a total of 76 products advertised as hair growthenhancer. 24 However, a method for simultaneous iden-tification of active substances that promote hair growth using high resolution mass spectrometry (HR-MS) has yet to be established with high specificity and accuracy.

The aim of present study was to develop a newUPLC-Q-Orbitrap-MS based method that can simul-taneously screen 12 hair growths to identify their specific fragmentation pathway and screen illegally manufactured products for hair growth. The Q-Orbitrap-MS, as one of HR-MS group, is frequently applied for screening (non)-targeted compound incomplex matrices using MS2 spectra. Recently, ne wanalogue of minoxidil were detected in adulterated dietary supplement by Q-Orbitrap-MS, which wasidentified as triaminodil.25 Since new analogue is synthesised through minor chemical modification of the parent structures in an attempt to avoid theirdetection by authorities, we established fragmentation pathways for 12 hair growth compounds using MS2 spectra of Q-Orbitrap-MS. This method was validated and applied to screen 14 dietary supplement advertised to be effective for the treatment of hair loss.

 

2. Experimental

 

2.1. Chemicals and reagents

Triaminodil was synthesized in our laboratories following a synthetic procedure for minoxidil in the literature.26 In the synthesis of triaminodil, pyrrolidine was used instead of piperidine. Six standard compounds (alfatradiol, dutasteride, finasteride, flutamide, minoxidil, spironolactone and testosterone propionate) were purchased from USP (Rockville, MD, USA). Cyproterone acetate and diphenylcyclopropenone were supplied by Sigma-Aldrich (St, Louis, MO, USA). Bimatoprost and methyltestosterone were obtained from BOC Sciences (New York, USA) and TCI (Tokyo, Japan), respectively. Standard stock solutions (1000 μg/mL) were prepared in methanoland stored at 4 °C. Deionised water was prepared by using a Milli-Q-water purification system (Millipore, Billerica, MA, USA) at 18.2 MΩ cm-1. Ammoniumacetate was supplied by Biopure (Cornwall, UK). Methanol (MeOH) and acetonitrile (ACN) for the HPLC grade were purchased from Merck KGaA(Darmstadt, Germany). All solvents in this study were HPLC grade and were filtered through apoly (vinylidenedifluoride) (PVDF) filter (0.2 μm).

 

2.2. Sample preparation

All 14 dietary supplements, advertised as hair growth enhancers, were obtained from online and of fline markets in the last two years. They were in the form of capsules, tablets, powders, liquids. About 1g of a homogenized sample was dissolved in 50 mL methanol and degassed in a sonication bath for 30 min. The extract was then centrifuged for 10 min at 3000 rpm to remove any extra matrix materials. Analiquot of the supernatant was filtered through a 0.2 mmpolytetrafluoroethylene (PTFE) filter and diluted to appropriate concentration with MeOH for subsequentinstrumental analyses.

 

2.3. UHPLC-Q-Orbitrap-MS analysis

The UHPLC-Q-Orbitrap-MS experiments wereconducted by using a Q Exactive Orbitrap mass spectrometer equipped with a Thermo Dionex UltiMate 3000 LC (Thermo Scientific, San Jose, CA). The eluents were 1 mM ammonium acetate in distilled water (A) and acetonitrile (B), and the column was a Poroshell 120 EC-C18 (100 × 2.1 mm, i.d. 2.7 µm) maintained at 30 °C. The injected sample volume was 1 mL, and the flow rate was 0.4 mL·min−1. Aftereluent B was initially maintained at 10 % for 2 min, the gradient profile was as follows: 2.0-6.0 min (A: 90-0 %, B: 10-100 %), 6.0-8.0 min (A: 0 %, B:100 %), 8.0-8.1 min (A: 0-90 %, B: 100-10 %), and 8.1-10.0min (A: 90 %, B:10 %). Full MS/ddMS2 (data-dependent MS2) experiments were conducted. The mass spectrometer parameters were as follows: ionsource, heated ESI (HESI); ion mode, positive10except for flutamide and alfatradiol; spray voltage, 3.5 kV(+) and 3.0 kV (−); capillary temperature, 320 ℃ C; sheath gas, 42 arbitrary units; auxiliary gas, 10 arbitrary units; probe heater temperature, 350 °C (+) and 300 °C (−); S-lens RF level, 50; resolution, 70,000(full scan), 17,500 (MS/MS); automatic gain control (AGC) target, 3e6 (full scan), 1e5 (MS/MS); scan range, 50 to 1000 m/z; maximum infusion time (IT), 100 ms (full scan) and 50 ms (MS/MS); microscans, 1; loop count, 5; MSX count, 1; Top N, 5; isolation window, 4 m/z; underfill ratio, 1.0 %; intensity threshold, 2e4; exclude isotopes, on; and dynamic exclusion, 10.0 s. The mass spectrometer was calibrated according to the manufacturer’s instructions. Data were analysed by using Xcalibur 3.0 software (Thermo Scientific, San Jose, CA).

 

2.4. Method validation

Several analytical parameters of target compoundsincluding specificity, limit of detection (LOD), limitof quantification (LOQ), linearity, precision and accuracy, recovery, and stability were evaluated by running three replicates. The stock solution of each component was added to a 10 mL volumetric flask inan amount of 500 μL to prepare a 50 μg/mL working solution. The spiking process was performed by dissolving the three-type (solid, liquid, cream) matrix-blank sample in MeOH during the sample preparation and then adding the working solution appropriately to the concentration required for the analysis. Thespecificity was confirmed by comparing the blank with the spiked standards. LODs and LOQs of targetcompounds were determined using signal-to-noise (S/N) ratio of the lowest detectable concentration of standard compounds spiked in matrices, and defined as the analyte amount that yields respective S/N ratios 3 and 10. Linearity was evaluated using the square of the linear correlation coefficient (R2) obtained by plotting the peak areas of six different concentrations. Serial dilutions for 10, 50, 100, 250, 500, and 1000 ng/mL were performed to obtain linearity except foral fatradiol. The linearity of alfatradiol was determined at 6 points of 200, 250, 500, 600, 750, and 1000 ng/mL. The inter- and intra-day parameters were also evaluated by performing three replicate experiments on three different days and in a day, respectively. Their accuracy was determined by comparing the average concentration calculated from the linearequation with the theoretical concentration. The precision was expressed by using the relative standard deviation (RSD). Percent recovery was calculated by comparing the peak areas of an analyte and its spiked standard. Stability of the target compounds wasevaluated for 24 h and 48 h at 4 °C.

 

3. Results and Discussion

 

3.1. Optimization of instrument conditions

To search the optimum compositions of mobile phase for chromatographic separation, several aqueoussolutions were examined with ACN including 5 m Mammonium acetate, 1 mM ammonium acetate, 5 m Mammonium formate, and 0.1 % formic acid. We could not observe any significant differences insensitivities of the target compounds in positive ion mode except for alfatradiol when 5 mM ammoniumacetate, 5 mM ammonium formate, and 0.1 % formicacid solutions. Alfatradiol was resoluted better when 5 mM ammonium acetate solution was used compared to 5 mM ammonium formate and 0.1 % formic acidsolutions. Peak tailings and sensitivities of all the compounds improved when 1 mM ammonium acetatesolution was used. The results made us to choose 1 mM ammonium acetate in distilled water (A) and ACN (B) as the compositions of mobile phase forrapid and efficient chromatographic separation of the target compounds.

Several mass parameters were investigated toattain suitable selectivity and sensitivity for the targetcompounds. Considering that mass selectivity and sensitivity move in opposite directions according t ofull width half maximum (FWHM) values, the full MS resolution was optimized at 70,000 FWHM for the majority of the analytes to obtain suitable selectivity and sensitivity. Likewise, the dd-MS2 resolution wasset at 17,500 FWHM to get suitable selectivity and sensitivity. The mass tolerance window was within 5 ppm considering the detection capability, signal intensity, and matrix interferences. All mass errors for the protonated molecules ([M + H]+) were ranged from -2.3 to 3.5 ppm, demonstrating that their massaccuracy was highly reliable in the UHPLC-Q-Orbitrap-MS (Table 1).

 

Table 1. Quasimolecular ions and MS2 fragment ions of hair growth remedies

BGHHBN_2019_v32n2_56_t0001.png 이미지

 

3.2. Fragmentation of hair growth

Fragmentations of 12 hair growth compounds were obtained from the full-MS/ddMS2 of Orbitrap-MS which is composed of a full MS scan followed by 5 data-dependent scans with normalized collision energy of 23-75 NCE. In Table 1, the characteristic fragmentions of significant abundance are donated for hair growth compounds. In order to better understand the observed fragmentation in the MS2 spectrum, a fragmentation pathway is proposed for each compoundas shown in Figs. 1 and 2. The fragmentation pathways were carefully envisaged using the following several factors in previous literature.27 First, the carbon atom(s) with branches or strain conducted toinductive cleavage. Second, stability depends oncarbocation order (tertiary > secondary > primary) and positive charge stabilized by resonance or inductive effects. Eight of 12 compounds were classified aschemical structure. Minoxidil and triaminodil had asimilar fragment ion pattern (1-8) with a difference of 14 Da depending on the chemical of R, which is described by the sequential cleavage in the R structure. Fragmentation pathways (9-14) of finasteride and dutasteride were distinguished by their amide group. Methyltestosterone and testosterone propionate werefound to produce common fragment patterns (15-22), eliminating R1, R2, ketone, and methyl group. Cyproterone acetate and spironolactone had fragment patterns depend on R1 and R2. Cyproterone acetategenerated specific ions 25 to 27 because of sequential breaks of acetyl, methyl and ketone, and chloride, whereas spironolactone showed 28 to 30 from ions m/z 341 [M-SCOCH3]+ which is dominated with higherintensity and stability. Fragmentation process of al fa-tradiol, bimatoprost, diphencyclopropenone, flutamide were proposed in Fig. S1, respectively. These fragment peaks can be possibly used as a transition in the MRM analysis.

 

BGHHBN_2019_v32n2_56_f0001.png 이미지

Fig. 1. Proposed fragmentation pathways for (a) minoxidil and triaminodil, (b) dutasteride and finasteride.

 

BGHHBN_2019_v32n2_56_f0002.png 이미지

Fig. 2. Proposed fragmentation pathways for (a) methyltesto-sterone and testosterone propionate, (b) cyproterone acetate and spironolactone.

 

3.3. Method validation

Specificity was investigated by comparing several types of matrix blanks and matrix spikes for the 12 target compounds. As shown in Fig. 3, no interferencepeaks were observed in the chromatograms. We thusconcluded that our UHPLC-Q-Orbitrap-MS experiments could provide a resolution high enough to distinguish the analytes from their isobaric ions. Several analytical parameters for the target compounds are indicated in Table S1 including LODs, LOQs, and linearity. LODs and LOQs for the target compounds in three types of matrices were determined as the lowest concentrations yielding S/N ratios of 3 and 10, respectively. LODs were in range of 0.05-50 ng/mL for three types of matrices, and the range of LOQs was from 0.17 to 167 ng/mL. The linearity wasevaluated by plotting the peak area corresponding tosix serial concentrations ranging from 10 to 1000 ng/mL. All the R2 values for the target compounds were higher than 0.995. The value indicates highly good linearity of the method. Accuracy and precision wereassessed at three concentrations (low, medium, and high) for intra- and inter-day comparisons. Accuracy was determined by the percent recovery while precision was evaluated with the inter-day repeatability and reproducibility using the relative standard deviation (RSD). The intra- and inter- day accuracy ranged from 88 to 112 % and 88 to 115 %, respectively (Table S2). The precision was within 5 % (intra-day) and 12 % (inter-day).

 

BGHHBN_2019_v32n2_56_f0003.png 이미지

Fig. 3. The extracted parent ion chromatograms of target compounds spiked in solid type matrix-blank sample by using UHPLC-Q-Orbirap-MS.

As shown in Table S3, mean recoveries of the target compounds from the three types of matrices were in ranges of 85-104 % (solid), 86-112 % (liquid), and 84-115 % (cream). In addition, their precisions were less than 6 %. These results demonstrate that this newly developed method is highly efficient to analyze the target compounds with relatively low matrix effects. Stability was examined at differentstorage times at 4 °C. Each RSD was within 8 % (24h) and 12 % (48 h), which indicates that all targetcompounds were stable during UHPLC-Q-Orbitrap-MS experiments (Table S4).

 

BGHHBN_2019_v32n2_56_f0004.png 이미지

Fig. 4. UHPLC-Q-Orbirap-MS/MS product ion spectra (high resolution) of triaminodil, monoxidil, and finasteride in (a) the standard solutions and (b) the adulterated products.

 

3.4. Sample application

Several types of 14 dietary supplements, advertised as hair loss remedies, were rapidly screened by UHPLC-Q-Orbitrap-MS. The target compounds were detected within 10 ppm mass tolerance and furtherconfirmed by performing MS/MS experiments As summarized in Table 2, approximately 21 % of samples (3/14) were adulterated with triaminodil, finasteride, and minoxidil in solid-type. The retentiontime of detected compounds were ranged from 1.9 to 5.6 min and mass accuracies were < 4.5 ppm, whichindicated positive results. Diagnostic fragment ions of samples corresponded to those of standard solutions. The amounts of the drugs or a drug analogue in the adulterated products ranged from 5.9 to 16.4 mg/g, which indicated the possibility of considerable risk to the health of public.

 

Table 2. Detected numbers of dietary supplements (n=14)

BGHHBN_2019_v32n2_56_t0002.png 이미지

 

4. Conclusions

A rapid and sensitive method for simultaneousidentification of 12 drugs or their analogues for the treatments of hair growth in dietary supplement wasestablished by UHPLC-Q-Orbitrap-MS/MS. Specific fragmentation pathways were proposed by interpreting the MS2 spectra of protonated ions. We have demonstrated the accuracy and practicality of the method by identifying 3 adulterated products from 14 products. Several drugs and a drug analogue were detected in adulterated products including minoxidil, triaminodil, and finasteride in amounts that ranged from 5.9 to 16.4 mg/g. These adulterated products contained large amounts of drugs and a drug analogueenough to cause serious side effects. This newly developed method will thus be highly useful forauthorities to screen dietary supplements advertised ashair loss remedies. Also, the results of this fragmentationstudy may be useful monitoring for rapid identification of new substances, which should contribute to efforts to safe guard food safety and public health.

 

Conflict of Interest

The authors declare that they have no conflict of interest.

 

Acknowledgements

This research was supported by Research Grants (15181MFDS521& MFDSAAT2017) from the Ministry of Food and Drug Safety (MFDS) in Korea. This work was supported by the Ministry of Food and Drug Safety and Kangwon National University. The authors thank the Central Laboratory of Kangwon National University for providing them with technicalassistance for the spectroscopic experiments.

References

  1. U. Blume-Peytavi, C. Kunte, A. Krisp, N. G. Bartels, U. Ellwanger, and R. Hoffmann, J. Dtsch. Dermatol Ges., 5, 391-395 (2007). https://doi.org/10.1111/j.1610-0387.2007.06295.x
  2. B. Patel, M. A. M. Velasco, F. T. Gutierrez, and D. Khesin, J. Bioequiv. Bioavail., 3, 1-9 (2017).
  3. S. Sheikh, A. Ahmad, S. M. Ali, M. U. Ahmad, M. Paithankar, D. Saptarshi, P. Kale, K. Maheshwari, H. V. Barkate, H. L. Pathak, M. Mushtaq, and I. Ahmad, J Clin. Exp. Dermatol., 6, 253 (2015).
  4. A. Adil and M. Godwin, J. Am. Acad. Dermatol., 77, 136-141 (2017). https://doi.org/10.1016/j.jaad.2017.02.054
  5. S. Inui and S. Itami, Exp Dermatol, 22, 168-171 (2013). https://doi.org/10.1111/exd.12024
  6. M. N. Pereira, H. L. Schulte, N. Duarte, E. M. Lima, L. L. Sa-Barreto, T. Gratieri, G. M. Gelfuso, and M. S. Cunha-Filho, Eur. J. Pharm. Sci., 96, 411-419 (2017). https://doi.org/10.1016/j.ejps.2016.10.016
  7. H. Wolff, T. W. Fischer, and U. B. Peytavi, Dtsch. Arztebl. Int., 113, 377-386 (2016).
  8. Y. Tsunemi, R. Irisawa, H. Yoshiie, B. Brotherton, H. Ito, R. Tsuboi, M. Kawashima, and M. Manyak, J. Dermatol., 43, 1051-1058 (2016). https://doi.org/10.1111/1346-8138.13310
  9. R. Sinclai, M. Wewerinke, and D. Jolley, Br. J. Dermatol., 152, 466-473 (2005). https://doi.org/10.1111/j.1365-2133.2005.06218.x
  10. A. Yazdabadi and R. Sinclair, Aust J Dermatol., 52, 132-134 (2011). https://doi.org/10.1111/j.1440-0960.2010.00735.x
  11. B. J. W. Dixson and M. J. Rantala, Arch. Sex. Behav., 45, 877-889 (2016). https://doi.org/10.1007/s10508-015-0588-z
  12. D. G. Barceloux and R. B. Palmer, Dis. Mon., 59, 226-248 (2013). https://doi.org/10.1016/j.disamonth.2013.03.010
  13. E. Marchei, D. D. Orsi, C. Guarino, M. C. Rotolo, S. Graziano, and S. Pichini, Cosmetics, 12, 1-8 (2016).
  14. P. M. Krishna, B. T. Rao, M. B. Raju, C. N. Rao, R. K. Kumar, and P. Venkateswarlu, J. Pharm. Res., 4, 2381-2383 (2011).
  15. R. Yuan and Y. Lin, Pharmacol Ther., 86, 191-198 (2000). https://doi.org/10.1016/S0163-7258(00)00039-5
  16. J. Wang, B. Chen, and S. Yao, Food Addit. Contam A, 25, 822-830 (2008). https://doi.org/10.1080/02652030801946553
  17. R. D. Sinclair and R. D. Dawber, Clin. Dermatol., 19, 167-178 (2001). https://doi.org/10.1016/S0738-081X(00)00128-0
  18. C. Shi, Z. Chen, Y. Chen, and R. Yu, Anal. Methods, 7, 6840-6849 (2015).
  19. C. Wu, Y. Jin, J. Zhang, Y. Ren, and Z. Jia, Chin. Chem. Lett., 24, 509-511 (2013). https://doi.org/10.1016/j.cclet.2013.03.037
  20. D. D. Orsi, M. Pellegrini, S. Pichini, D. Mattioli, E. Marchei, and L. Gagliardi, J. Pharm. Biomed Anal., 48, 641-648 (2008). https://doi.org/10.1016/j.jpba.2008.06.008
  21. P. Ptacek, J. Macek, and J. Klima, J. Chromatogr B, 738, 305-310 (2000). https://doi.org/10.1016/S0378-4347(99)00543-5
  22. P. B. Phapale, H. W. Lee, M. S. Lim, E. H. Kim, S. D. Kim, J. Park, M. Lee, S. K. Hwang, and Y. R. Yoon, J. Chromatogr B, 878, 1718-1723 (2010). https://doi.org/10.1016/j.jchromb.2010.04.029
  23. S. C. Patterson and T. Ramstad, K. A. Mills, Il Farmaco., 60, 547-554 (2005). https://doi.org/10.1016/j.farmac.2005.04.002
  24. H. N. Park, J. H. Lee, S. Park, S. Y. Beak, Y. Lee, and H. Kang, J. Forensic Sci. Int., 284, 129-135 (2018). https://doi.org/10.1016/j.forsciint.2017.12.042
  25. J. H. Lee, H. N. Park, H. J. Park, N. S. Kim, S. Park, J. Lee, and S. Y. Beak, Food Addictives & Contaminants: Part A, 35, 2-9 (2018). https://doi.org/10.1080/19440049.2017.1397295
  26. M. Maltes, J. Org. Chem., 60, 2436-3442 (1995). https://doi.org/10.1021/jo00113a025
  27. J. Do, E. Noh, S. Yoon, J. H. Lee, S. Park, S. Mandava, S. Y. Baek, and J. Lee, Arch. Pharm. Res., 40, 736-745 (2017). https://doi.org/10.1007/s12272-017-0921-2