Next Article in Journal
Lipidomics Profiling of Metformin-Induced Changes in Obesity and Type 2 Diabetes Mellitus: Insights and Biomarker Potential
Previous Article in Journal
Metformin: The Winding Path from Understanding Its Molecular Mechanisms to Proving Therapeutic Benefits in Neurodegenerative Disorders
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

A Systematic Study of Yiqi Qubai Standard Decoction for Treating Vitiligo Based on UPLC-Q-TOF/MS Combined with Chemometrics, Molecular Docking, and Cellular and Zebrafish Assays

1
School of Medicine, Tongji University, Shanghai 200331, China
2
School of Pharmacy, Naval Medical University, Shanghai 200433, China
3
Department of Clinical Pharmacy, Shanghai General Hospital, Shanghai Jiaotong University School of Medicine, Shanghai 200080, China
4
School of Pharmacy, Shanghai Jiaotong University, Shanghai 200240, China
5
Department of Pharmacy, Shanghai Fourth People’s Hospital Affiliated to Tongji University School of Medicine, Shanghai 200434, China
6
School of Pharmacy, Anhui University of Chinese Medicine, Hefei 230012, China
*
Authors to whom correspondence should be addressed.
Pharmaceuticals 2023, 16(12), 1716; https://doi.org/10.3390/ph16121716
Submission received: 15 October 2023 / Revised: 27 November 2023 / Accepted: 7 December 2023 / Published: 11 December 2023
(This article belongs to the Special Issue Chemometrics in Pharmaceutical Research)

Abstract

:
The Yiqi Qubai (YQ) formula is a hospital preparation for treating vitiligo in China that has had reliable efficacy for decades. The formula consists of four herbs; however, the extraction process to produce the formula is obsolete and the active ingredients and mechanisms remain unknown. Therefore, in this paper, fingerprints were combined with the chemometrics method to screen high-quality herbs for the preparation of the YQ standard decoction (YQD). Then, the YQD preparation procedure was optimized using response surface methodology. A total of 44 chemical constituents, as well as 36 absorption components (in rat plasma) of YQD, were identified via UPLC-Q-TOF/MS. Based on the ingredients, the quality control system of YQD was optimized by establishing the SPE-UPLC-Q-TOF/MS identification method and the HPLC quantification method. Network pharmacological analysis and molecular docking showed that carasinaurone, calycosin-7-O-β-d-glucoside, methylnissolin-3-O-glucoside, genkwanin, akebia saponin D, formononetin, akebia saponin B, and apigenin may be the key active components for treating vitiligo; the core targets associated with them were AKT1, MAPK1, and mTOR, whereas the related pathways were the PI3K-Akt, MAPK, and FoxO signaling pathways. Cellular assays showed that YQD could promote melanogenesis and tyrosinase activity, as well as the transcription and expression of tyrosinase-associated proteins (i.e., TRP-1) in B16F10 cells. In addition, YQD also increased extracellular tyrosinase activity. Further efficacy validation showed that YQD significantly promotes melanin production in zebrafish. These may be the mechanisms by which YQD improves the symptoms of vitiligo. This is the first systematic study of the YQ formula that has optimized the standard decoction preparation method and investigated the active ingredients, quality control, efficacy, and mechanisms of YQD. The results of this study lay the foundations for the clinical application and further development of the YQ formula.

1. Introduction

Vitiligo is the most common hypopigmentation disorder. It is an acquired disease characterized by progressive loss of melanocytes. Its global prevalence is 0.5–2%, with no racial or sex differences [1]. Vitiligo is a complex disease whose pathogenesis remains unknown but may be related to genetic variation, autoimmunity, oxidative stress, the autophagy of melanocytes, and/or neuroendocrine function [2]. The goal of vitiligo treatment is to suppress depigmentation and stimulate repigmentation in the patients. Currently, the available treatments for vitiligo mainly include glucocorticoids, immunosuppressants, phototherapy, and surgical methods [3,4,5,6]. Unfortunately, the available therapies are not effective for all patients, and some have wide-ranging side effects that suppress further treatment [7]. Therefore, increasing research has focused on finding new multi-target therapeutic medicine with fewer side effects to fit the complex treatment mechanisms of vitiligo and the need for long-term medication.
There is a long history of using traditional Chinese medicine (TCM) to treat vitiligo in China, with which there have been positive results and few adverse effects [8]. The theory of using TCM in the treatment of vitiligo is that it benefits “qi” by activating “blood” and dispelling pathologic “wind” [9]. “Qi” is a theory of Chinese medicine that represents the substance that constitutes the human body, maintains human life activities, and promotes Yang, the driving force of metabolism and physiological function of the human body [10]. The Yiqi Qubai (YQ) formula is a TCM that has been utilized in the dermatology department of Shanghai General Hospital for decades [11] and has had remarkable curative results [12,13]. It consists of four herbs, Astragali radix (Huangqi, HQ), Akebiae fructus (Yuzhizi, YZZ), Leonuri fructus (Chongweizi, CWZ), and Caragana sinica roots (Jinquegen, JQG), in a weight ratio of 1:1:0.75:1.5. Each of the four herbs can soothe the liver and regulate “qi”, activate “blood”, and dispel “wind” [14,15,16]. The upregulation of mushroom tyrosinase (TYR) activity and melanogenesis of YQ was demonstrated using serum pharmacological methods [17]. Although YQ has been used clinically for many years, there is a lack of systematic research on it. Firstly, the decoction process does not follow the standard decoction principles, which may affect the uniformity of the preparation quality. Secondly, the composition of the decoction is unclear, which is not conducive to an in-depth study of its material basis and mechanisms. In addition, the identification methods in the existing quality standards are simple and lack quantitative methods to ensure quality.
Essentially, modern TCM preparations rely on the decoction, which is the oldest Chinese medicinal compound that can reflect the material basis of a prescription; therefore, the quality of the decoction is crucial. The Chinese Pharmacopoeia Commission formally proposed the “standard decoction” in 2016, which is characterized by the authenticity of the herbs, the uniformity of the extraction process, and the rigor of quality control [18]. The preparation principle of the standard decoction is based on the theory of TCM, following the boiling principle of TCM decoction, adopting standardized technology, and decoction with water. The preparation method is recommended to be optimized with reference to “the management standard of traditional Chinese medicine decoctions in medical institutions” published by the State Drug Administration of China [19]. The standard decoction of TCMs has a uniform process and a product with consistent quality; this will assist in the preparation of granules and allow preparations to be widely used in clinical practice. In order to ensure the accuracy of the medication and the consistency of the administered dose, it is necessary to establish a method for the preparation of a YQ standard decoction (YQD) that could be the basis for the clinical use of the medication.
As the raw material of a standard decoction, the quality of the TCM herbs is affected by origin, harvesting season, storage time, processing method, etc., which largely determines the efficacy of TCM preparations. The TCM fingerprint is a comprehensive, quantifiable technology that can provide more detailed information regarding the type and contents of the TCM components [20,21]. As data mining methods have developed, chemometrics has gained more and more attention because it can simplify complex data and uncover hidden patterns. The combination of chemometric methods such as hierarchical cluster analysis (HCA) and principal component analysis (PCA) with fingerprints is often used to classify and evaluate the quality of TCM herbs [22,23].
In recent years, mouse B16F10 melanoma cells have been used as a sensitive and reliable model for the quantitative analysis of melanin [24]. In addition, zebrafish have a high degree of chromosome homology with humans and have a translucent body that allows direct visualization of the development of melanin streaks, making it a suitable model organism for studying melanogenesis in vivo [25]. Melanocyte-specific enzymes, such as TYR and the TYR-associated proteins TRP-1 and TRP-2, are involved in melanin synthesis. TYR is the main rate-limiting enzyme in melanin synthesis, and any factor that activates TYR activity can effectively promote the production of melanin [26].
In this study, HPLC fingerprints were combined with chemometrics to establish a quality control system of raw herbs to screen high-quality herbs for the preparation of YQD. The preparation process of YQD was then optimized via response surface methodology (RSM). A UPLC-Q-TOF/MS technique was used to analyze the chemical constituents of YQD and the absorption components (in rat plasma). Next, the SPE-UPLC-Q-TOF/MS identification method and HPLC quantification method were developed to improve the quality control system. Network pharmacology analysis and molecular docking based on the identified components were used to predict the active ingredients, targets, and mechanisms of vitiligo treatment. Then, we investigated the effects of YQD on melanogenesis, TYR activity, and TRP-1 expression in B16F10 cells, as well as extracellular TYR activity, to explore the possible mechanisms. In addition, the efficacy was investigated by studying its effect on melanogenesis in zebrafish. To the best of our knowledge, this is the first systematic study on the preparation, composition, quality control, efficacy, and mechanism of YQD. The results of this study lay the foundations for the clinical application and further development of YQD. The whole study process is summarized in Figure 1.

2. Results and Discussion

2.1. Multivariate Statistical Analysis of the HPLC Fingerprint in Raw Herbs

2.1.1. HPLC-DAD Analysis and Similarity Evaluation

Biological and abiotic factors affected the chemical quality of plants and showed wide geographical variation. In order to ensure the stability of the raw materials, it is necessary to study the origins before establishing a standard decoction preparation process. Using HPLC-DAD, 15 batches of HQ (H-1–H-15), YZZ (Y-1–Y-15), CWZ (C-1–C-15), and JQG (J-1–J-15) from different origins were detected and their chromatograms were recorded. The precision, repeatability, and stability of the analytical method were validated. The results indicated that the fingerprint analysis method was stable and reliable. The similarity of the fingerprint chromatograms was assessed. Table S1 shows the similarity matching results for the fingerprints. A total of 11, 10, 4, and 9 common peaks were identified on the fingerprint spectra for HQ, YZZ, CWZ, and JQG, respectively (Figure 2A). The intragroup similarities of HQ, YZZ, CWZ, and JQG were 0.919–0.998, 0.773–0.997, 0.932–0.999, and 0.883–0.999, respectively. As shown in the results, each of the four herbs, when obtained from different production areas, differed within the group, especially YZZ and JQG; these differences may be due to the different geographies, harvesting times, and sources of the herbs.

2.1.2. Hierarchical Cluster Analysis (HCA)

In order to evaluate the similarity of the samples based on their fingerprint data, HCA was used to categorize them [22]. As shown in Figure 2B, HQ samples were divided into four clusters, of which H-11 (Zhangye, Gansu) and H-15 (Shanxi) were classified as group I and group II, respectively, H-12 (Zhangye, Gansu) and H-14 (Neimeng) were classified as group III, and the remaining samples (Dingxi, Gansu) were classified as group IV. The YZZ samples were divided into three clusters, of which Y-11 (Yichang, Hubei) was a single group. CWZ samples were divided into three clusters, C-11 (Nanyang, Henan) and C-14 (Anhui) were classified as group I and group II, respectively, and the other samples (Nanyang, Henan; Chifeng, Neimeng; and Fuxin, Liaoning) were classified as group III. The JQG samples were divided into three clusters, J-15 (Sichuan) was alone in group I, J-10 (Yichang, Hubei) and J-11 (Yichang, Hubei) were in group II, and the remaining samples (Jingmen/Xiangyang/Yuan’an, Hubei) were classified as group III. The findings reveal that samples from similar origins were more likely to be placed into the same group, meaning the content of the ingredients was comparable. Nevertheless, some samples from the same origin were categorized independently, probably because of batch variations. There was not a clear regional clustering of YZZ, which may be related to the fact that it is collected from a variety of sources (YZZ is the dry fruit of Akebia quinata (Thunb.) Decne., A. trifoliata (Thunb.) Koidz. var. australis (Diels) Rehd., and A. trifoliata (Thunb.) Koidz.). Differences in the chemical makeup of the herbs can affect their quality; therefore, the quality of herbs from various locations needs to be assessed to determine the best source of raw materials.

2.1.3. Principal Component Analysis (PCA)

As a statistical method of dimension reduction, PCA recombines the original variables into independent comprehensive variables [23]. Calculating the composite score of principal components (PCs) and sorting could help to determine the relative quality of each sample. In this work, PCA was performed using SIMCA (ver. 14.1). The eigenvalue and variance contribution rate of the correlation coefficient are shown in Table S2. With an eigenvalue >0.99, the first three PCs in HQ had the cumulative variance contribution rate of 81.609%; the first two PCs in YZZ had the cumulative variance contribution rate of 85.467%; the first one PC in CWZ was 73.803%; and the first three PCs in JQG had the cumulative variance contribution rate of 88.242%. These PCs contained the most information on the variables and fully reflected the original data. The factor load matrix and loadings plot (Figure S1) revealed that peaks 5 (calycosin-7-O-D-glucoside), 7, 8, and 11 in HQ; peaks 5 (cryptochlorogenic acid), 4, 7, and 2 in YZZ; peaks 1, 3, and 2 (chlorogenic acid) in CWZ; and peaks 7, 1, and 9 in JQG all demonstrated a higher contribution to PC1 than the other peaks. These peaks might be the primary chemical characteristic peaks that can be utilized to differentiate the quality of the different varieties of the four herbs.
The PC scores and comprehensive score values of the 15 batches of HQ, YZZ, CWZ, and JQG were calculated using the initial eigenvalue and weights of the PCs. A higher score means better quality. As shown in Table 1, the comprehensive score of H-12 (Zhangye, Gansu), Y-11 (Yichang, Hubei), C-13 (Chifeng, Neimeng), and J-11 (Yichang, Hubei) were highest for the HQ, YZZ, CWZ, and JQG samples, respectively. Similar to the HCA, samples with the highest and lowest comprehensive scores were individually clustered in the PCA score chart (Figure 2C). The ecological settings and herb quality are closely connected, and standard decoctions should be made from suitable and stable herbs. Our results showed that H-12, Y-11, C-13, and J-11 were better quality and should be used for the subsequent preparation of YQD.

2.2. Optimization of YQD Preparation Process

2.2.1. Parameter Optimization through RSM

RSM is a collection of statistical and mathematical methodologies for developing and optimizing procedures in which many variables impact a desired response; it has the aim of optimizing that response [27]. The Box–Behnken design (BBD), a type of RSM, reduces the number of tests and investigates the interactions between the influencing factors to obtain the optimal process parameters quickly and effectively [28]. The experimental parameters were optimized using BBD based on the results of the single-factor studies of the solid–liquid ratio, decoction time, and concentration process. Three independent variables, namely, the first solid–liquid ratio (X1), the first extraction time (X2), and the second extraction time (X3), were selected to be included in the response surface modeling, and a three-factor, three-level (−1, 0, +1) BBD was carried out. Synthesizing all indicators (calycosin-7-O-β-d-glucoside and chlorogenic acid contents of the extract and its total solid content), the overall desirability (OD) value was used to express the overall impact results. According to the Hassan method [29], the OD value was calculated as follows:
di = (Yi − Ymin)/(Ymax − Ymin)
The di value of each indicator was used to calculate the total OD as follows [30]:
OD = (d1 × d2…dx) 1/n
A total of 17 experiments were conducted (n = 17) and the BBD test design and response values are shown in Table 2. To determine the predicted response values, the following equations were fitted using a quadratic mode:
OD = 0.8612 + 0.3228X1 + 0.1158X2 + 0.1248X3 − 0.1440X1X2 − 0.1435X1X3 − 0.0205X2X3 − 0.2179X12 − 0.0159X22 − 0.0058X32
where X1 is the first liquid–solid ratio, X2 is the first extraction time, and X3 is the second extraction time.
A further ANOVA was performed on the regression model (Table S3). The table showed that the model was statistically significant (p < 0.05) and that the linear terms (X1, X2, and X3) had a significant (p < 0.05) effect on extraction. The interaction surfaces and contour plots of X1, X2, and X3 are shown in Figure 3A–D. The plots showed that when the first extraction time (X2) or the second extraction time (X3) was fixed at a certain level, the OD value increased rapidly alongside the increase in the first liquid–solid ratio (X1) and then decreased slightly. This indicated that the increase in the first liquid–solid ratio was not always conducive to the response and that saturation may have occurred. By adjusting the liquid–solid ratio appropriately, two compounds could be yielded more efficiently. The interaction surfaces and contour plots between X2 and X3 are shown in Figure S2. In cases where the second extraction time was fixed at a certain level, the OD score increased with the first extraction time; however, this increase was small. Perhaps the increased extraction time resulted in a more efficient destruction of the plant cells, resulting in a higher extraction rate; however, this interaction was not the primary cause.
The optimal extraction conditions for YQD were calculated by solving the model’s extreme values and analyzing the contours. The first liquid–solid ratio was 8.086, the first extraction time was 59.998 min, and the second extraction time was 39.998 min. The optimal extraction process resulted in an OD value of 1.06. Due to actual manufacturing considerations, the optimum conditions of YQD were as follows: a liquid–solid ratio of 8, the first extraction time was 60 min, and the second extraction time was 40 min.

2.2.2. Experimental Validation of the Model

The optimized YQD extraction process was as follows: 100 g of the prescription amounts of the HQ, YZZ, CWZ, and JQG pieces were blended together and extracted with eight times the amount of water. After being soaked for 30 min, fire-boiled, and decocted on a small fire for 60 min, the mixture was immediately filtered through a 200-mesh screen. The filter residue was then extracted with seven times the amount of water, fire-boiled, decocted on a small fire for 40 min, and immediately filtered again. The twice-combined extraction decoction was concentrated to 0.2 g/mL under reducing pressure at 60 °C.
Six validation experiments were conducted to verify the above optimal solution. The OD value was 0.96 and the average contents of calycosin-7-O-β-d-glucoside and chlorogenic acid were 0.107 mg/g and 0.338 mg/g, respectively. Table S4 contains the detailed extraction rates.

2.3. Identification of the Constituents of YQD and Optimization of the Quality Control Methods

2.3.1. Identification of Chemical Constituents

UPLC-Q-TOF/MS was used for qualitative analysis of YQD. First, the MS data were matched with the Natural Products HR-MS/MS Spectral Library 1.0 database. Compounds not included in the database were identified in the literature and mass spectrometry pyrolysis laws. Finally, 44 chemical compositions (C1–C44) were identified in YQD (Figure 4A,B and Table 3). Among them, C12 from HQ and C13, C24, C30, and C31 from YZZ had higher response values. Typical compounds of flavonoids and triterpenoid saponins were used to analyze the cleavage pattern as follows:
C12 (C22H22O10, 30.86 min) showed a quasi-molecular ion at m/z 447.1269. The neutral loss of glucose (C6H10O5, −162 Da) in the MS2 spectrum formed a fragment ion m/z 285.0750, which further lost a methyl (•CH3, −15 Da), and the methoxy group (CO, −28 Da) formed the fragment ions m/z 270.0517 and m/z 242.0541. Fragment ion m/z 285.0750 may also lose methanol (CH3OH, −32 Da) and further lose the methoxy group (CO, −28 Da) to obtain the fragment ion m/z 225.0522. Fragment m/z 137.0217 may be obtained via the RDA reaction of the compound or the deglycosylated glycoside element. C12 was identified as calycosin-7-O-β-d-glucoside by comparing its MS2 cleavage behavior and retention time with those of the reference compounds.
C24 (eluted at 42.79 min) showed a quasi-molecular ion at m/z 1091.5812 [M + H]+, with a molecular formula of C53H86O23 and fragment ions at m/z 959.5186, 813.4651, 651.4084, and 489.3562, which were generated by the quasi-molecular ion peaks losing one molecule of xylose, one molecule of deoxymannose, and one molecule of glucose in turn; the fragment ions at m/z 489.3562 was arjunolic acid sapogenin. Using data from reference [31,32], C24 was identified as abekia saponin. The MS2 spectrum and fragmentation pathway of C12 and C24 are shown in Figure 5. All of the above results are beneficial for the further discovery of the absorbed ingredients after oral administration.

2.3.2. Identification of Absorbed Components in Rat Plasma

A total of 36 components of YQD were identified from rat plasma using background reduction, baseline correction, and molecular characterization extraction methods, of which, 4 were prototype components (P1–P4) and 32 were metabolites (M1–M32). Results are shown in Figure 4C and Table 4. The prototype components were stachydrine (P1), carasinaurone (P2), calycosin-7-O-β-d-glucoside (P3), and akebia saponin B (P4). Among them, P1 and P3 are the quality control indicators of HQ and CWZ, respectively, in the Chinese Pharmacopoeia (2022) [33]. P2 and P4 are the signature active ingredients of JQG and YZZ, respectively [34,35]. Due to low oral bioavailability and a short half-life, the components with high responses in vitro (C13, C24, C30, and C31) were not found as prototypes or metabolites in vivo [36,37,38]. By studying the metabolic pathways in detail, we found that the metabolites were mainly from glucuronidation, sulfation, hydroxylation, and demethylation reactions of the flavonoids in YQD.
Glucuronidation: M2 showed a quasi-molecular ion at m/z 621.1487 in the negative ion mode that was 176 Da higher than calycosin-7-O-β-d-glucoside at m/z 447.1286, suggesting that calycosin-7-O-β-d-glucoside had undergone a glucuronidation reaction. The fragment ions of M2 at m/z 283.0620 and 268.0344 were like calycosin-7-O-β-d-glucoside; therefore, it was speculated that M2 was a glucuronidation of calycosin-7-O-β-d-glucoside. Similarly, the metabolites M10, M23, M6, and M27 were speculated to be glucuronide products of apigenin, formononetin, genkwanin, and methylnissolin, respectively.
Sulfation: M32 showed a quasi-molecular ion at m/z 347.0230 in the negative ion mode that was 80 Da higher than formononetin at m/z 267.0663, indicating that a sulfation reaction had occurred in formononetin. Fragment ions of M32 at m/z 267.0650 and 252.0431 were similar to formononetin; therefore, it was speculated that M32 was a sulfated metabolite of formononetin.
Glucuronidation + Sulfation: in the negative ion mode, M5 showed a quasi-molecular ion at m/z 525.0359, a fragment ion at m/z 269.0465, and was speculated to be the sulfated and glucuronidated product of apigenin because the quasi-molecular ion was 176 and 80 Da higher than that of apigenin at m/z 269.0455, and the fragment ion was similar to apigenin.
Demethylation: M24 was tentatively identified as a demethylated product by comparing the quasi-molecular ion between M24 and formononetin; the gap was 14 Da in negative ion mode.
In addition, M30 showed a quasi-molecular ion at 303.1212 and fragment ions at m/z 303.1024, 167.0699, 133.0677, and 123.0434 and so was hypothesized to be isoflavan transformed by methylnissolin [39]. The [M − H] m/z of M9 was 541.0301, which was 176, 80, and 17 Da higher than that of apigenin. It was speculated that M9 was the hydroxylation, sulfation, and glucuronidation product of apigenin. The possible metabolic pathways are shown in Figure 6.

2.3.3. Establishment of the SPE-UPLC-Q-TOF/MS Identification Method

Identification is an important part of the quality standard of TCM. Traditional TLC or HPLC methods often have low sensitivity, many interfering factors, and complicated pretreatment methods that cannot meet the identification needs of TCM compound preparations. In this study, the SPE-UPLC-Q-TOF/MS technique was used to identify the characteristic peaks of YQD; the separation of the components was better after SPE extraction. Asperulosidic acid (C3) and carasinaurone (C11) in the H2O eluting phase; neochlorogenic acid (C5), chlorogenic acid (C7), and cryptochlorogenic acid (C8) in the 30% methanol elution phase; and calycosin-7-O-β-d-glucoside (C12) and calceolarioside B (C13) in the 60% methanol elution phase were selected as the characteristic components because they were from four separate herbs with high response values, complete separation, and good peak shapes.
The results of sample identification showed that, in the H2O elution section, the retention times of C3 and C11 were 21.81 and 36.12 min, respectively; in the 30% methanol elution section, the retention times of C5, C7, and C8 were 23.46, 27.17, and 27.91 min, respectively; and in the 60% methanol elution section, the retention times of C1 and C13 were 36.41 and 38.38 min, respectively (Figure 7I).

2.3.4. Establishment of HPLC Quantitative Method

An HPLC method was established for the determination of calycosin-7-O-β-d-glucoside (CAG, C12) and chlorogenic acid (CA, C7). CAG is a representative component of HQ [33]. CA may promote melanogenesis in melanoma cells [40] and had a higher response value in the UPLC/MS chromatogram. Therefore, these two substances were selected as quantitative components to compensate for the lack of quantitative analysis in the existing quality standards.
As can be seen from Figure 7II, the response value of the target peak was high, the separation degree was good, and the other peaks did not interfere with the measurement. The standard curves were Y = 0.292X–0.3234, r = 0.9998 (CAG) and Y = 0.099X–0.0567, r = 1.000 (CA). The linear ranges for CAG and CA were 5.21–208.4 μg/mL and 10.32–412.8 μg/mL, respectively. The CAG and CA control solutions were taken separately, and the peak areas were measured via successive injections six times according to the established HPLC conditions. The results showed that the RSDs of CAG and CA were 0.25% and 1.63%, respectively (n = 6), which indicated that the precision of the instrument was good. For the evaluation of repeatability, the RSDs of the peak areas from six samples of the same batch were calculated and were 0.77% and 3.31%, respectively (n = 6), indicating that the repeatability of the method was good. The stability of the samples was determined by re-analyzing the RSDs of the peak area from one prepared sample at 0, 2, 4, 6, 8, 12, and 24 h. The RSDs of CAG and CA were 0.54% and 3.97%, showing that the sample solutions were stable within 24 h at room temperature. The recoveries were calculated using three concentrations (equivalent to 0.5, 1.0, and 1.5 times the original amounts in the samples), and the recovery rate of CAG was 92.25–98.63% and RSD was 0.87–2.30%; the recovery rate of CA was 97.78–105.30% and the RSD was 2.97–3.46%, which indicated that the method had a high recovery rate and good accuracy.
Using the newly established method, ten batches of YQD samples were evaluated. The results (Table S5) show that the content of CAG and CA were 0.333–0.484 mg/g and 0.275–0.383 mg/g, respectively. The RSD was less than 15%, indicating that the method could conduct the quantitative analysis of the two components. Under the analytical conditions, the response value and resolution of CAG and CA met the quantitative requirements well.

2.4. Network Pharmacology Analysis and Molecular Docking Studies

2.4.1. Bioactive Ingredients and Potential Targets of YQD

Network pharmacology emphasizes the multi-level and multi-angle discovery of potential targets and signaling pathways. This includes identifying active ingredients, their targets, disease targets, and screening the core genes and then constructing networks within the compound–target–disease–signaling pathway and evaluating the effects of compounds on these networks [41,42]. The ingredients comprising TCM preparations are numerous; however, the components present at the highest amounts or those that can be readily absorbed by the body may be the active components. Therefore, the compounds found in rat plasma were subjected to further network pharmacology analysis. As the metabolites of daudzein and equol (M1, M3, M4, M7, M8, M14, M17, M19, M21, and M22) are derived from diverse sources [43], they were not specific metabolic components of YQD and were therefore excluded. Additionally, apigenin (C42), a prototype component that produces many metabolites, and four components (C13, C24, C30, and C31) with high response values in the UPLC/MS chromatogram were added to the pharmacological prediction model. The 31 compounds used for network pharmacology analysis are shown in Table 5. In total, 466 targets related to the compounds were found (Table S6), including AKT1, TNF, ALB, mTOR, MAPK1, IL2, PIK3CA, IGF1, and others.

2.4.2. PPI and CTD Network Analysis

The target genes related to vitiligo were searched, and 78 overlapping targets were identified by matching the 466 YQD targets with the 1059 vitiligo genes (Figure 8I(A), Table S7). Based on these 78 targets, a PPI network was successfully constructed and visualized (Figure 8I(B), Table S8). The top 20 targets with degree values were selected as potential core targets (Figure 8I(C)) and included AKT1, mTOR, and MAPK1. Phosphorylated AKT1 has been reported to increase the transcription of MITF, a major regulator of melanin synthesis, which further upregulates the expression of TYR, TRP-1, and TRP-2 and ultimately stimulates melanin synthesis [44]. mTOR activation was beneficial for melanocyte protection from dendritic loss under oxidative stress, thereby maintaining melanin levels in skin pigmentation [45]. Early expression of phosphorylated MAPK1 can upregulate the transcriptional activity of MITF, thereby promoting melanin synthesis [46].
To investigate the underlying relationship mechanism of YQD on the therapy of vitiligo, a compound–target–disease (CTD) network was constructed, as shown in Figure 8I(D). There were 111 nodes in the network, with 31 components, 78 targets, 1 illness (vitiligo), and 1 drug (YQD). The average degree value of the compounds was 6.79, with 16 compounds having a degree value greater than 6 (Table S9). Among them, M30, M24, C42, M32, M10, M13, M15, M6, P3, and C31 ranked in the top ten; the seven metabolites were derived from C42, C16, C32, C12, and C21, respectively. Research has shown that apigenin (C42) may target TRP-1 and TRP-2 through the MAPK pathway, thereby promoting melanin production [47]. Akebia saponin D (C31), formononetin (C32), and genkwanin (C21) have demonstrated good anti-inflammatory, antioxidant, and immunomodulatory activities [48,49,50,51].

2.4.3. GO and KEGG Enrichment Analysis

According to GO classification enrichment analysis, the vitiligo targets were categorized as being involved in the stress response, activation of cytokines, biological regulation, and cell communication (Figure 8II(A)). The top 20 vitiligo-related pathways included the TNF, FoxO, PI3K-Akt, and mTOR signaling pathways according to KEGG enrichment (Figure 8II(B) and Table S10).
To validate that the ingredients could regulate the targets through the aforementioned pathways to produce the therapeutic effects, we constructed a compound–target–pathway (CTP) network using 41 core target genes (in the top 20 pathways) (Figure 8II(C)). The network had 92 nodes, including 31 components, 20 pathways, and 41 targets. M30, C42, M10, M15, M24, M13, M6, and P3 had higher than average degree values in the network, whereas the targets with higher degree values were AKT1, MAPK1, PIK3CA, and mTOR. The results were consistent with the PPI and CTD network analyses. The PI3K-Akt, TNF, MAPK, and FoxO signaling pathways were enriched for more than 10 targets (Table S11) that have been reported in the literature as potentially critical pathways for vitiligo treatment [52,53,54].
The results of network pharmacology showed that YQD is a multi-component, multi-target, and multi-pathway treatment for vitiligo. Calycosin-7-O-β-d-glucoside (C12), methylnissolin-3-O-glucoside (C16), genkwanin (C21), akebia saponin D (C31), formononetin (C32), and apigenin (C42) may be the key components that regulate multiple targets to perform therapeutic roles.

2.4.4. Molecular Docking Study

Molecular docking is an alternative technique for predicting the binding affinity of molecules to the active site of a particular receptor without the use of animals; it enables high-throughput screening of components [55,56]. It is generally accepted that the lower the energy when the conformation of the ligand binding to the receptor is stable, the great-er the probability of action. Target prediction in network pharmacology is heavily reliant on the predictive capacity of network prediction tools, and the data are primarily generated from experiments [57]. However, because of the restricted number of experimentally confirmed drug targets, certain less-studied TCMs projected fewer targets, causing their constituents to be excluded from the active ingredients.
In this study, the predicted active compounds with high (M30, M24, and C42), medium (M32, M10, M15, M6, and C31), and low (P4, C30, and P2) degree values in network pharmacology were docked with potential core targets (AKT1, MAPK1, and mTOR) to verify the predicted results and explore the possibility of compounds with a low degree value participating in the binding of targets. As shown in Figure 9A, the binding energies of the 11 ligands to the core targets ranged between −4.7 and −11.8 kcal/mol, which implies that all of them had good binding activities to the proteins except for P2/MAPK1 (−4.7 kcal/mol). We found that the strongest binding ability to AKT1 was with M6 and M30, which formed hydrogen bonds with residues ALA-58 (Figure 9B) and TYR-272 (Figure 9C) that had binding energies of −10.2 and −9.0 kcal/mol, respectively. The component with the strongest binding capacity to MAPK1 was M15, and the hydrogen bond binding sites were at ARG-359, VAL-21, GLN-355, ARG-91, and ASP-100 (Figure 9D). M10/MAPK1 ranked second and formed hydrogen bonds with residues ASP-167, GLY-169, and TYR-36 (Figure 9E). The binding energy of the 11 ligands to mTOR ranged between −8.2 and −11.8 kcal/mol, with C31/mTOR having the lowest binding energy and forming hydrogen bonds with five residues (Figure 9F). In addition, P4 is predicted to bind strongly to mTOR, forming hydrogen bonds with GLU-2032, ASP-2102, and ARG-73 (Figure 9G). The results suggest that the components with lower predicted activity determined using network pharmacology may also have a strong binding capacity to core proteins and thus exert therapeutic effects.
In conclusion, combining the above computational results and considering the metabolic pathways of the metabolites, we hypothesized that C11 (carasinaurone, JQG), C12 (calycosin-7-O-β-d-glucoside, HQ), C16 (methylnissolin-3-O-glucoside, HQ), C21 (genkwanin, CWZ), C31 (akebia saponin D, YZZ), C32 (formononetin, HQ), C40 (akebia saponin B, YZZ), and C42 (apigenin, HQ) may be the pharmacodynamic markers of YQD that deserve further in-depth study.

2.5. Effects of YQD on Melanogenesis In Vitro and In Vivo

2.5.1. Effects of YQD on Melanogenesis in B16F10 Cells

The viability of the murine B16F10 melanoma cell line was examined in the presence of different concentrations of YQD using the CCK-8 assay. The results showed that B16F10 cells treated with YQD at concentrations of 1–1000 µg/mL for 48 h did not display any changes in cell viability when compared with the control (Figure 10A).
To determine the effects of YQD on pigmentation, we examined the melanin content of B16F10 cells treated with different concentrations (1–1000 µg/mL) of YQD and a positive control drug (8-MOP) for 48 h. As shown in Figure 10B, YQD increased the melanin content in B16F10 cells in a dose-dependent manner.
To investigate its effect on the melanogenic pathway in B16F10 cells, we explored whether YQD affected TYR activity and the expression of the melanin biosynthetic gene (TRP-1). As shown in Figure 10C, the activity of cellular TYR increased significantly in a concentration-dependent manner following incubation with YQD. The results of the qRT-PCR experiments showed that the mRNA level of TRP-1 was elevated after YQD treatment (Figure 10D); moreover, Western blotting showed that the protein level was also increased significantly in a dose-dependent manner (Figure 10E). These results indicate that YQD induced pigmentation by enhancing TYR activity and upregulating TYR family genes.

2.5.2. Effects of YQD on Extracellular TYR Activity

An in vitro TYR activation assay showed that YQD could activate mushroom TYR and that the activation rate was concentration dependent (Table S12). The TYR activation rate for 8 mg/mL of YQD reached 79.66%, which was equivalent to the activation rate for 0.5 mg/mL of 8-MOP.

2.5.3. Effects of YQD on Melanogenesis in Zebrafish

After observing the effects of YQD in mammalian cells, we wanted to determine whether a model organism had a similar response. Zebrafish embryos were used to investigate the role of YQD in promoting melanogenesis in vivo. Firstly, the acute toxicity and dose were examined by observing pericardial edema, delayed yolk absorption, and mortality in zebrafish embryos within 96 h postfertilization (hpf) in the presence of YQD. The statistical results are shown in Table S13, and the condition of zebrafish 72 hpf in the presence of YQD is depicted in Figure S3A. The results show that concentrations of YQD below 1000 μg/mL had no influence on zebrafish embryos. However, a considerable number of juvenile fish died in both groups at 2000 μg/mL (the mortality rates were 50% and 100%). All embryos had died 72 hpf when the concentration was increased to 4000 μg/mL. The calculated LC50 was 1997 μg/mL (Figure S3B).
Next, we explored the effects of YQD on melanogenesis. Zebrafish embryos were separated into eight groups and constantly treated with blank culture medium, varying doses of YQD (200–1200 μg/mL), or a positive control drug (8-MOP) from 24 hpf to 84 hpf. The density and distribution of the melanin granules were observed and recorded every day. We found that the melanin density in the head and tail gradually increased as the embryos grew. The difference between groups peaked at 84 hpf; melanin production in the head, dorsal spine, and abdomen of zebrafish treated with YQD was more noticeable than in other parts of the body, whereas the survival and development of juvenile fish were not affected in the high-concentration treatment group (Figure 11A). We measured the melanin density in the tails of zebrafish and found that the melanin density in the YQD-treated groups (400–1200 μg/mL) was significantly higher than that in the control group (Figure 11B).

3. Materials and Methods

3.1. Materials and Regents

HQ, YZZ, CWZ, and JQG were collected from different parts of China. The origins and batch numbers are shown in Table S14. Methanol and acetonitrile (both liquid chromatography grade) were purchased from Merck (Darmstadt, Germany) and formic acid was acquired from ANPEL Lab Tech. (Shanghai, China). 8-Methoxypsoralen (8-MOP) was purchased from Macklin Biochemical Co., Ltd. (Shanghai, China). The solid-phase extraction (SPE) cartridges (Oasis HLB, 3 mL) were purchased from Waters (Milford, MA, USA). Pure water (18.2 MΩ/cm) was deionized using a Milli-Q system (Millipore, Bedford, MA, USA).

3.2. Statistical Analysis of the HPLC Fingerprints

Chromatographic data from 15 different batches of each herb were mined and analyzed. Similarity analysis was performed using the “Chinese Materia Medica chromatographic fingerprint similarity evaluation system” (ver. 2012) [20]. The time window width was set to 0.1 min, and the median method was selected for full-spectrum peak matching. The simulated average fingerprint R was generated, the similarity between the sample and the average fingerprint R was calculated, and then the common peaks that existed in each herb were observed [21]. A chemometric analysis including HCA and PCA was used to examine the similarities and differences between the samples [22]. SIMCA software (ver. 14.1) was used to categorize the samples, and dendrograms were drawn to characterize the classification results. Significant PCs were extracted under the condition that the corresponding eigenvalues were greater than 0.99. A score plot was used to assess sample variation and PC comprehensive scores were calculated to evaluate the sample quality.

3.3. Preparation of the HQ, YZZ, CWZ, and JQG Decoction

An amount of 100 g of the HQ, YZZ, CWZ, or JQG pieces was extracted using 10 times the amount of water. After being soaked for 30 min and decocted for 60 min, the mixture was immediately filtered through a 200-mesh screen. The filtrate residue was then extracted with 8 times the amount of water, decocted for 40 min, and immediately filtered again. The twice-combined extraction decoction was concentrated to 0.2 g/mL through reducing pressure [19].

3.4. Optimization and Validation of the YQD Preparation Conditions

A three-factor with three levels (−1, 0, +1) BBD was applied to optimize the critical parameters, such as the first solid–liquid ratio, first extraction time, and second extraction time [28]. The three responses used to evaluate the influence of the extraction parameters were the calycosin-7-O-β-d-glucoside and chlorogenic acid contents of the extract and its total solid content, and the OD value was used to express the overall impact results. Multiple regression analysis was used to develop the best-fit polynomial quadratic model [22]. The p-value, lack of fit value, and coefficient of determination (R2) from the ANOVA were used to assess the model adequacy. Experimental design, statistics, response surfaces, and contour plots were obtained using Design-Expert software (ver. 12.0).
After optimizing the preparation conditions of YQD, six validation experiments were conducted for the optimal solution, with a first liquid–solid ratio of 8, a second liquid–solid ratio of 7, a first extraction time of 60 min, and a second extraction time of 40 min. The experimental values were compared with the predicted ones to validate the model.

3.5. Preparation of the Sample and Standard Solution

HPLC and UPLC-TOF/MS samples: an amount of 1 mL of each decoction sample (0.2 g/mL) was accurately extracted, centrifuged at 13,000 rpm for 10 min, and then the supernatant was filtered through a 0.22 μm membrane before use.
SPE-UPLC-TOF/MS samples: an amount of 500 μL of YQD was added to 500 μL of deionized water, mixed thoroughly, and centrifuged at 12,000 rpm for 15 min. Then, 1 mL of the supernatant was subjected to the activated SPE column and eluted successively with H2O, 30%, 60%, and 100% methanol (1 mL for each). The flow rate was controlled uniformly (about 10 drops/min). The eluent was collected and dried, redissolved in 100 μL of 50% methanol, vortexed for 3 min, sonicated for 3 min, and centrifuged at 12,000 rpm for 15 min, and the supernatant was then taken for analysis.
Standard solutions: stock solutions of calycosin-7-O-β-d-glucoside and chlorogenic acid were separately prepared in methanol. The stock solutions of each analyte were diluted in 60% methanol to prepare calibration solutions and filtered through 0.22 μm membranes before use.

3.6. Animals and Drug Administration

Sprague–Dawley rats (180–200 g, 6–7 weeks old) were provided by Shanghai Silaike Experimental Animal Co., Ltd. (Shanghai, China, approval number: 2021-0009). The rats were kept in a room (24 ± 2 °C with 60% ± 5% humidity) with a 12 h light/dark cycle for one week and had free access to water and conventional laboratory food to acclimate. The rats were randomly assigned to the control group or YQD group (n = 6 each). The rats in the YQD group were orally administered with YQD (dissolved in purified water) once a day at a dose of 15 g/kg (decoction pieces weight/body weight) for 5 consecutive days. The same volume of saline solution was given to the rats in the control group. Blood was collected 1 h after the last dosing.
All experiments were conducted in accordance with the Guiding Principles for the Care and Use of Laboratory Animals and were approved by the Institutional Animal Care and Use Committee (IACUC) of Shanghai Jiao Tong University (approval number: A2021145, approval date: 23 December 2021). All work was performed following the recommendation of the ARRIVE guidelines.

3.7. Plasma Collection and Preparation

The blood samples were collected and centrifuged at 4000 rpm for 10 min at 4 °C. The plasma samples were then transferred to Eppendorf tubes and kept at −80 °C for further analysis. To precipitate protein, 200 μL of plasma was treated with 800 μL of methanol, vortex mixed for 1 min, and centrifuged at 12,000 rpm for 15 min at 4 °C. The supernatant was transferred to a new tube and evaporated to dryness. For further analysis, the residue was redissolved in 80 μL of 50% methanol and centrifuged at 12,000 rpm for 15 min.

3.8. UPLC-Q-TOF/MS and HPLC Analysis

For qualitative analysis, the samples were analyzed on an Agilent 1290 Infinity UPLC system (Milford, MA, USA) equipped with a Waters UPLC column (ACQUITY HSS T3, 2.1 × 100 mm, 1.8 μm). The mobile phase consisted of solvent A (0.1% formic acid aqueous solution) and solvent B (0.1% formic acid acetonitrile), and the gradient elution program was set as follows: 0–10 min, 100% A; 10–27 min, 100–85% A; 27–62 min, 85%–60% A; 62–70 min, 60%–5% A; 70–75 min, 5% A; and 75.01–78 min, 100% A. The flow rate was 0.3 and the column temperature was 30 °C. The injection volume was set at 5 μL.
Mass spectrometric analyses were performed using an AB 4600 Q-TOF-MS/MS system (AB SCIEX, Framingham, MA, USA) equipped with an ESI source. The parameters for the mass spectrometry were as follows: scan range, 50–1800 m/z (MS) and 50–1250 m/z (MS/MS); ion spray voltage floating, −4.5 kV (negative mode) and 5.0 kV (positive mode); declustering potential, 100 V; source temperature, 500 °C; nebulizer gas, 50 psi; auxiliary heater gas, 50 psi; curtain gas, 35 psi; and collision energy, 10 eV (MS) and ±40 eV (MS/MS). To obtain the fragment information, N2 was used as the collision gas.
Quantitative analysis was performed on an UltiMate 3000 (Thermo, Waltham, MA, USA) liquid chromatograph system using an YMC C18 column (250 mm × 4.6 mm, 5 μm). The mobile phase was composed of solvent A (0.1% aqueous formic acid) and solvent B (acetonitrile), with the gradient elution program: 0–15 min, 95–85% A; 15–20 min, 85–80% A; 20–40 min, 80–65% A; 40–45 min, 65–95% A; 45–60 min, 95% A. The column temperature was maintained at 25 °C and flow rate was set at 1 mL/min. An amount of 10 μL of the solutions was injected and the wavelength was set to 260 nm.

3.9. Network Pharmacological Analysis

The components discovered via UPLC-Q-TOF/MS were used to establish the chemical information database of YQD for network pharmacology research. We searched for the component-related genes with a possibility > 0.1 or z-score > 1.5 using the SwissTargetPrediction (https://www.swisstargetprediction.ch) and PharmMapper (https://www.lilab-ecust.cn/pharmmapper (both accessed on 2 April 2022)) online tools. Vitiligo-related genes were prepared from the GeneCards (https://www.genecards.org (accessed on 2 April 2022)), OMIM (https://www.omim.org (accessed on 24 March 2022)), DisGeNET (https://www.disgenet.org (accessed on 2 April 2022)), TTD (https://www.bidd.nus.edu.sg/group/cjttd (accessed on 24 March 2022)), and DrugBank (https://www.drugbank.ca (accessed on 24 March 2022)) databases. These target genes were standardized using the Uniprot database (https://www.uniprot.org (accessed on 2 April 2022)). Then, component-related targets and vitiligo-related targets were crossed using Venny (ver. 2.1.0) to obtain the potential targets of YQD in vitiligo treatment.
The protein–protein interaction (PPI) network was identified using the String Database (https://www.string-db.org (accessed on 28 May 2022)) and then imported into and visualized using Cytoscape (ver. 3.7.2). The Gene Ontology (GO) function and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichment were performed via the Metascape database (https://metascape.org (accessed on 28 May 2022)). The compound–target–disease (CTD) and component–target–pathway (CTP) networks were constructed and analyzed using Cytoscape software (ver. 3.7.2).

3.10. Molecular Docking

The 3D structures of the ingredients and the crystal structures of targets were obtained from the Pubchem (https://www.ncbi.nlm.nih.gov) and PDB (https://www.rcsb.org (both accessed on 30 May 2022)) databases, respectively. The Autodock Tool (ver. 1.5.6) was employed for ingredient–protein molecular docking and calculating the binding affinity. The docking results were visualized using PyMOL (ver. 2.5.2).

3.11. In Vitro Cells Assay

3.11.1. Cell Cultures and Viability Assay

Mouse B16F10 melanoma cells (Chinese Academy of Sciences, Shanghai, China) were cultured in RPMI 1640 (Gibco, MA, USA) medium, supplemented with 10% FBS and a penicillin–streptomycin antibiotic mix in 5% CO2 at 37 °C.
The cell viability was measured using CCK-8 solution (Macklin, Shanghai, China) [24]. Briefly, B16F10 cells were seeded in 96-well plates (1 × 104 cells/well) for 24 h. The cells were then treated with 8-MOP (100 μmol/L) as a positive control and different concentrations of YQD for 48 h. Then, the culture medium was removed and 10 μL of CCK-8 solution was added to each well. After incubation for 1 h at 37 °C, the absorbance was measured at 450 nm using a Varioskan Flash spectrometer (Thermo Scientific, Waltham, MA, USA). All the assays were performed in triplicate and the cell viability was calculated as (Asample/Acontrol) × 100%.

3.11.2. Melanin Measurement and TYR Activity Assay

Melanin release from B16F10 cells was measured using a previously described method with slight modifications [24]. B16F10 cells were seeded at a density of 1 × 104 cells/well in a 6-well plate overnight. Then, the cells were treated with different concentrations of YQD and 8-MOP (100 μmol/L) for 48 h. Following the treatment, the cells were washed with ice-cold PBS and dissolved in 200 μL of 1 mol/L NaOH (with 10% DMSO) at 70 °C for 1 h. The absorbance at 490 nm was measured via a multiplate reader. The protein concentration of each sample was determined using the BCA Protein Assay Kit (Biomed, Beijing, China). The melanin content of the YQD-treated cells was calculated relative to the untreated cells.
The assay for TYR activity was carried out as previously described [26], with a slight modification. B16F10 cells were incubated with test samples for 48 h and lysed with 1% Triton X-100 solution containing 1% sodium deoxycholate for 15 min at 0 °C. Each lysate was centrifuged at 12,000 rpm for 10 min to obtain the supernatant. After protein quantification and adjustment, 100 µL of the cell lysate was mixed with 100 µL of 0.1% L-DOPA (Yuanye Bio-Technology, Shanghai, China) and incubated at 37 °C for 30 min. The absorbance was measured at 475 nm using a multiplate reader, and the TYR activity in YQD-treated cells was presented relative to the activity in untreated cells.

3.11.3. Quantitative Real-Time PCR

Total cellular RNA was prepared from B16F10 cells treated with YQD using TRIzol reagent (solarbio, Beijing, China). cDNA was synthesized using Hifair III (Yeasen, Shanghai, China) according to the manufacturer’s instructions. qRT-PCR was performed at 95 °C for 5 min, followed by 40 cycles at 95 °C for 10 s and 60 °C for 30 s. The PCR primers used were as follows: 5′-GCCATCTTTGTCACCAGGTT-3′ (forward) and 5′-GCTCGAGCAAACTTCCATTC-3′ (reverse) for TRP1; 5′-CGCAGCCACTGTCGAGTC-3′ (forward) and 5′-GTCATCCATGGCGAACTGGT-3′ (reverse) for β-actin. The relative expression was determined by normalizing the data to the β-actin mRNA level.

3.11.4. Western Blot Analysis

B16F10 cells were treated with different concentrations of YQD for 48 h. Cells were then lysed in cold RIPA buffer for 30 min on ice. The lysates were centrifuged at 12,000 rpm for 15 min at 4 °C and the supernatants were collected. The proteins were separated by 10% SDS-PAGE and then transferred onto PVDF membranes. The membranes were blocked with 5% skim milk, incubated with appropriate antibodies (TRP-1, ABclonal, Wuhan, China) at 4 °C overnight, and then incubated with the corresponding secondary antibodies for 2 h at room temperature. The targeted proteins were detected using the ECL detection system (GE Healthcare, Beijing, China) and visualized with the ChemiDoc XRS+ Imaging System (Bio-Rad, Hercules, CA, USA).

3.12. Extracellular TYR Activity

The Tyrosine Assay Kit (Yuanye Bio-Technology, Shanghai, China) was used to measure the extracellular TYR activity. L-DOPA, TYR, PBS, and different concentrations of YQD were added to the tubes according to the manufacturer’s instructions. The mixture was incubated for 10 min at 37 °C in a water bath; then, TYR was added and incubated for 2 min at 37 °C. The absorbance was determined at 475 nm and the activation rate was calculated.

3.13. In Vivo Melanogenic Effect in Zebrafish

Wild-type zebrafish embryos (AB strain) were collected from Shanghai Jiao Tong University. The embryos were incubated with E3 solution as a control and different concentrations of YQD from 24 h postfertilization (hpf) to 96 hpf at 28 °C. Two replications were set up for all groups (n = 30 in each group). The pericardium edema rate, yolk absorption delay rate, and mortality rate of zebrafish within 96 hpf were observed and recorded to investigate the toxicity of YQD. The LC50 was calculated using GraphPad Prism (ver. 8.0.2).
In the melanin measurement assay [24], the embryos were incubated with 8-MOP (100 µmol/L) as positive control or variable concentrations of YQD from 24 hpf to 84 hpf. Three replications were set up for all groups (n = 10 in each group). A stereomicroscope (Olympus Instruments, Tokyo, Japan) was used to observe and record the melanin in the heads and tails of the zebrafish daily. The melanin density in the tail was measured using Image-Pro Plus (ver. 6.0).

3.14. Statistical Analysis

All data were reported as the mean and standard deviation. The results were analyzed and illustrated using GraphPad Prism (ver. 8.0.2) software. The level of statistical significance is described in the respective figure legends (* p < 0.05, ** p < 0.01, and *** p < 0.001, as indicated by a one-way ANOVA).

4. Conclusions

In this study, a combination of fingerprints and chemometrics was used to select the optimal herbs for the preparation of YQD. The preparation process for YQD was optimized, and the chemical constituents and the absorbed constituents were identified. The SPE-UPLC-Q-TOF/MS identification method and HPLC quantification method were established to improve the quality standard of YQD. It was predicted that several signature components and metabolites from the four herbs might be important active components, reflecting the rationality of the YQ formula. In addition, B16F10 cell and zebrafish assays confirmed the role of YQD in promoting melanogenesis. The possible mechanism is that YQD can increase the expression and transcription of TRP-1 and enhance the activity of TYR, which is meaningful for the treatment of vitiligo. The results of this study fill many gaps in the research of YQ formulations and lay a good foundation for the clinical application and further development of YQD.

Supplementary Materials

The following supporting information can be downloaded at: https://www.mdpi.com/article/10.3390/ph16121716/s1. Tables S1: Similarity analysis of the fingerprint chromatograms; Table S2: Eigenvalue and variance contribution rate; Table S3: ANOVA table for the OD (Box-Behnken design); Table S4: Experimental validation results; Table S5: The content of CAG and CA in 10 batches of YQD; Table S6: Compounds and target used for network pharmacology analysis; Table S7: Overlapping of drug-related targets and disease-related targets; Table S8: Data of PPI network diagrams; Table S9: The degree value in the CTD network; Table S10: Data of KEGG analysis; Table S11: The degree value in the CTP network; Table S12: Activation rate of YQD to mushroom tyrosinase. Table S13: The toxicity of YQD to zebrafish embryos; Table S14: Origin and batch information of HQ, YZZ, CWZ, and JQG. Figure S1: Loading plot of HQ (A); YZZ (B); CWZ (C); and JQG (D); Figure S2: The interaction surface (A) and the contour diagram (B) of the first extraction time (X2) and the second extraction time (X3); Figure S3: (A) The condition of zebrafish 72 hpf in the presence of YQD; (B) The calculated LC50 of YQD in zebrafish.

Author Contributions

L.C.: conceptualization, software, investigation, resources, and writing—original draft; C.M.: investigation and resources; W.S.: investigation; C.Y.: validation; J.W.: investigation and resources; Z.W.: project administration and supervision; Y.L.: project administration and supervision; G.F.: funding acquisition, project administration, supervision, conceptualization, and writing—review and editing. All authors have read and agreed to the published version of the manuscript.

Funding

This work is supported by the National Natural Science Foundation of China (81793289); National Key Clinical Specialty Discipline Construction Program of China (02. DY11.03.19.03); special fund for clinical research of Wu Jieping Medical Foundation (320.6750.2021-2-75).

Institutional Review Board Statement

Ethics Approval: Rats were maintained, handled, and bred according to standard protocols from the Institutional Animal Care and Use Committee of Shanghai Jiao Tong University (IACUC A2018075). Adult and larvae zebrafish were maintained, handled, and bred according to standard protocols from the Institutional Animal Care and Use Committee of Shanghai Jiao Tong University (IACUC 202101310). Animal studies were carried out in accordance with the Guiding Principles for Care and Use of Laboratory Animals and followed the guidelines for the ethical review of laboratory animal welfare People’s Republic of China National Standard GB/T 35892-2018.

Informed Consent Statement

Not applicable.

Data Availability Statement

The data used to support the findings of this study are included in the article.

Conflicts of Interest

The authors declare no conflict of interest.

Abbreviations

YQDYiqi Qubai standard decoction
TCMTraditional Chinese medicine
HCAHierarchical cluster analysis
PCAPrincipal component analysis
RSMResponse surface methodology
BBDBox–Behnken design
UPLC-Q-TOF/MSUltrahigh performance liquid chromatography coupled with quadrupole-time of flight mass spectrometry
HPLCHigh performance liquid chromatography
SPESolid-phase extraction
CTDCompound–target–disease
PPIProtein–protein interaction
CTPComponent–target–pathway
KEGGKyoto Encyclopedia of Genes and Genomes
GOGene Ontology
mTORMammalian target of rapamycin
MAPK1Mitogen-activated protein kinase 1
PI3KPhosphoinositide 3-kinase
AKTSerine-threonine kinase
FoxOForkhead box, class O
TNFTumor necrosis factor
MITFMicrophthalmia-associated transcription factor
TYRTyrosinase
TRP-1Tyrosinase-related proteins 1
TRP-2Tyrosinase-related proteins 2

References

  1. Krüger, C.; Schallreuter, K.U. A review of the worldwide prevalence of vitiligo in children/adolescents and adults. Int. J. Dermatol. 2012, 51, 1206–1212. [Google Scholar] [CrossRef] [PubMed]
  2. Bergqvist, C.; Ezzedine, K. Vitiligo: A Review. Dermatology 2020, 236, 571–592. [Google Scholar] [CrossRef] [PubMed]
  3. Zubair, R.; Hamzavi, I.H. Phototherapy for Vitiligo. Dermatol. Clin. 2020, 38, 55–62. [Google Scholar] [CrossRef] [PubMed]
  4. Razmi, T.M.; Afra, T.P.; Parsad, D. Vitiligo surgery: A journey from tissues via cells to the stems. Exp. Dermatol. 2019, 28, 690–694. [Google Scholar] [CrossRef] [PubMed]
  5. Seneschal, J.; Duplaine, A.; Maillard, H.; Passeron, T.; Andreu, N.; Lassalle, R.; Favary, C.; Droitcourt, C.; Taïeb, A.; Ezzedine, K. Efficacy and Safety of Tacrolimus 0.1% for the Treatment of Facial Vitiligo: A Multicenter Randomized, Double-Blinded, Vehicle-Controlled Study. J. Investig. Dermatol. 2021, 141, 1728–1734. [Google Scholar] [CrossRef] [PubMed]
  6. Wada-Irimada, M.; Tsuchiyama, K.; Sasaki, R.; Hatchome, N.; Watabe, A.; Kimura, Y.; Yamasaki, K.; Aiba, S. Efficacy and safety of i.v. methylprednisolone pulse therapy for vitiligo: A retrospective study of 58 therapy experiences for 33 vitiligo patients. J. Dermatol. 2021, 48, 1090–1093. [Google Scholar] [CrossRef] [PubMed]
  7. Picardo, M.; Dell’Anna, M.L.; Ezzedine, K.; Hamzavi, I.; Harris, J.E.; Parsad, D.; Taieb, A. Vitiligo. Nat. Rev. Dis. Primers 2015, 1, 15011. [Google Scholar] [CrossRef]
  8. Chen, L.L.; Chen, S.G.; Sun, K.K.; Sun, P.; Zou, Y.X.; Zhan, Z.S.; Wang, J.F. Research Progress on Treatment of Vitiligo with Traditional Chinese Medicine. J. Basic Chin. Med. 2022, 28, 822–826. [Google Scholar] [CrossRef]
  9. Shi, N.; Chen, Y.J.; Wang, J.; Ni, H. Clinical observation on the effect of Zengse Pill in treating patients with vitiligo of qistagnancy and blood-stasis syndrome type. Chin. J. Integr. Med. 2008, 14, 303–306. [Google Scholar] [CrossRef]
  10. Chen, Z.; Liu, L.; Gao, C.; Chen, W.; Vong, C.T.; Yao, P.; Yang, Y.; Li, X.; Tang, X.; Wang, S.; et al. Astragali Radix (Huangqi): A promising edible immunomodulatory herbal medicine. J. Ethnopharmacol. 2020, 258, 112895. [Google Scholar] [CrossRef]
  11. Zhu, G.D. Treatment of vitiligo based on the principle of pacifying liver by resolving stasis and activating blood circulation, plus exorcising “wind”. An observation on therapeutic effects in 100 cases. J. Tradit Chin. Med. 1982, 2, 71–75. [Google Scholar] [PubMed]
  12. Zhang, C.; Zhou, L.; Huang, J.; Shi, W. A combination of Yiqiqubai granule and 308-nm excimer laser in treatment of segmental vitiligo: A prospective study of 233 patients. J. Dermatol. Treat. 2017, 28, 668–671. [Google Scholar] [CrossRef] [PubMed]
  13. Yang, K.Y. The Effect of 308 nm Excimer Light Combined with Yiqiqubai Granule on the Treatment of Vitiligo. Master’s Thesis, Shanghai Jiaotong University, Shanghai, China, 2016. [Google Scholar]
  14. Meng, Q.; Niu, Y.; Niu, X.; Roubin, R.H.; Hanrahan, J.R. Ethnobotany, phytochemistry and pharmacology of the genus Caragana used in traditional Chinese medicine. J. Ethnopharmacol. 2009, 124, 350–368. [Google Scholar] [CrossRef]
  15. Ma, B.; Yang, S.; Tan, T.; Li, J.; Zhang, X.; Ouyang, H.; He, M.; Feng, Y. An integrated study of metabolomics and transcriptomics to reveal the anti-primary dysmenorrhea mechanism of Akebiae Fructus. J. Ethnopharmacol. 2021, 270, 113763. [Google Scholar] [CrossRef] [PubMed]
  16. Huang, Y.X.; Wang, J.; Wang, Y.; Kuang, W.; Xie, M.J.; Zhang, M. Pharmacological mechanism and clinical study of Qiming granules in treating diabetic retinopathy based on network pharmacology and literature review. J. Ethnopharmacol. 2023, 302, 115861. [Google Scholar] [CrossRef] [PubMed]
  17. Wu, J.Q.; Sun, Y.; Han, X.D.; Zhou, F.J.; Zhu, G.D.; Liu, G.L. Effects of extract of Yiqiqubai Compound on mushroom tyrosinase activity and melanogenesis. Chin. J. Tradit. Chin. Med. Pharm. 2013, 28, 2792–2794. [Google Scholar]
  18. Commission, C.P. Technical Requirements for Quality Control and Standard Formulation of Chinese Herbal Formula Granules (Draft for Comment). 2016. Available online: http://www.chp.org.cn/view/ff808081559d1ea301565ac24962560d?a=XWJX (accessed on 5 August 2016).
  19. Chen, S.-L.; Liu, A.; Li, Q.; Toru, S.; Zhu, G.-W.; Sun, Y.; Dai, Y.-T.; Zhang, J.; Zhang, T.-J.; Takehisa, T.; et al. Research strategies in standard decoction of medicinal slices. Chin. J. Chin. Mater. Med. 2016, 41, 1367–1375. [Google Scholar] [CrossRef]
  20. Liang, W.; Chen, W.; Wu, L.; Li, S.; Qi, Q.; Cui, Y.; Liang, L.; Ye, T.; Zhang, L. Quality Evaluation and Chemical Markers Screening of Salvia miltiorrhiza Bge. (Danshen) Based on HPLC Fingerprints and HPLC-MS(n) Coupled with Chemometrics. Molecules 2017, 22, 478. [Google Scholar] [CrossRef]
  21. Li, B.; Lima, M.; Nie, Y.; Xu, L.; Liu, X.; Yuan, H.; Chen, C.; Dias, A.C.; Zhang, X. HPLC-DAD Fingerprints Combined with Multivariate Analysis of Epimedii Folium from Major Producing Areas in Eastern Asia: Effect of Geographical Origin and Species. Front. Pharmacol. 2021, 12, 761551. [Google Scholar] [CrossRef]
  22. Shi, W.; Wu, Z.; Wu, J.; Jia, M.; Yang, C.; Feng, J.; Lou, Y.; Fan, G. A comprehensive quality control evaluation for standard decoction of Smilax glabra Roxb based on HPLC-MS-UV/CAD methods combined with chemometrics analysis and network pharmacology. Food Chem. 2023, 410, 135371. [Google Scholar] [CrossRef]
  23. Fan, X.; Hong, T.; Yang, Q.; Wang, D.; Peng, J.; Xiao, W.; Yang, X.; Hu, X.; Yu, C.; Du, S.; et al. Quality assessment of fried licorice based on fingerprints and chemometrics. Food Chem. 2022, 378, 132121. [Google Scholar] [CrossRef]
  24. Molagoda, I.M.N.; Karunarathne, W.A.H.M.; Park, S.R.; Choi, Y.H.; Park, E.K.; Jin, C.-Y.; Yu, H.; Jo, W.S.; Lee, K.T.; Kim, G.-Y. GSK-3β-Targeting Fisetin Promotes Melanogenesis in B16F10 Melanoma Cells and Zebrafish Larvae through β-Catenin Activation. Int. J. Mol. Sci. 2020, 21, 312. [Google Scholar] [CrossRef]
  25. Singh, A.P.; Nüsslein-Volhard, C. Zebrafish stripes as a model for vertebrate colour pattern formation. Curr. Biol. 2015, 25, R81–R92. [Google Scholar] [CrossRef] [PubMed]
  26. Sun, L.J.; Guo, C.Y.; Yan, L.T.; Li, H.J.; Huo, X.P.; Xie, X.; Hu, J. Syntenin regulates melanogenesis via the p38 MAPK pathway. Mol. Med. Rep. 2020, 22, 733–738. [Google Scholar] [CrossRef] [PubMed]
  27. Alshammari, F.; Alam, M.B.; Naznin, M.; Javed, A.; Kim, S.; Lee, S.H. Profiling of Secondary Metabolites of Optimized Ripe Ajwa Date Pulp (Phoenix dactylifera L.) Using Response Surface Methodology and Artificial Neural Network. Pharmaceuticals 2023, 16, 319. [Google Scholar] [CrossRef] [PubMed]
  28. Ji, P.; Li, C.; Wei, Y.; Wu, F.; Liu, S.; Hua, Y.; Yao, W.; Zhang, X.; Yuan, Z.; Wen, Y. Screening study of blood-supplementing active components in water decoction of Angelica sinensis processed with yellow rice wine based on response surface methodology. Pharm. Biol. 2020, 58, 1167–1176. [Google Scholar] [CrossRef]
  29. Hassan, E.E.; Parish, R.C.; Gallo, J.M. Optimized formulation of magnetic chitosan microspheres containing the anticancer agent, oxantrazole. Pharm. Res. 1992, 9, 390–397. [Google Scholar] [CrossRef]
  30. Wei, W.U.; Cui, G.; Bin, L.U. Optimization of multiple evariables: Application of central composite design and overall desirability. Chin. Pharm. J. 2020, 35, 530–533. [Google Scholar]
  31. Oda, K.; Matsuda, H.; Murakami, T.; Katayama, S.; Ohgitani, T.; Yoshikawa, M. Adjuvant and haemolytic activities of 47 saponins derived from medicinal and food plants. Biol. Chem. 2000, 381, 67–74. [Google Scholar] [CrossRef]
  32. Jiang, D.; Gao, Q.P.; Shi, S.P.; Tu, P.F. Triterpenoid Saponins from the Fruits of Akebiae quinata. Chem. Pharm. Bull. 2006, 54, 595–597. [Google Scholar] [CrossRef]
  33. Commission, C.P. Pharmacopoeia of the People’s Republic of China; China Medical Science Press: Beijing, China, 2020; pp. 252–315. [Google Scholar]
  34. Jin, Q.; Han, X.H.; Hong, S.S.; Lee, C.; Choe, S.; Lee, D.; Kim, Y.; Hong, J.T.; Lee, M.K.; Hwang, B.Y. Antioxidative oligostilbenes from Caragana sinica. Bioorg. Med. Chem. Lett. 2012, 22, 973–976. [Google Scholar] [CrossRef] [PubMed]
  35. Ma, B.L.; Yang, S.L.; Li, J.M.; Ouyang, H.; He, M.Z.; Feng, Y.L.; Tan, T. A four-step filtering strategy based on ultra-high-performance liquid chromatography coupled to quadrupole-time-of-flight tandem mass spectrometry for comprehensive profiling the major chemical constituents of Akebiae Fructus. Rapid Commun. Mass Spectrom. 2019, 33, 1464–1474. [Google Scholar] [CrossRef]
  36. Xue, Z.; Yang, B. Phenylethanoid Glycosides: Research Advances in Their Phytochemistry, Pharmacological Activity and Pharmacokinetics. Molecules 2016, 21, 991. [Google Scholar] [CrossRef] [PubMed]
  37. Chen, W.; Zheng, Y.; Yan, L.; Yuan, M.; Ouyang, L.; Li, J.; Zhou, G.; Zhong, R. Validated UHPLC-MS/MS method for simultaneous determination of four triterpene saponins from Akebia trifoliata extract in rat plasma and its application to a pharmacokinetic study. Biomed. Chromatogr. 2019, 33, e4585. [Google Scholar] [CrossRef]
  38. Li, P.; Peng, J.; Li, Y.; Gong, L.; Lv, Y.; Liu, H.; Zhang, T.; Yang, S.; Liu, H.; Li, J.; et al. Pharmacokinetics, Bioavailability, Excretion and Metabolism Studies of Akebia Saponin D in Rats: Causes of the Ultra-Low Oral Bioavailability and Metabolic Pathway. Front. Pharmacol. 2021, 12, 621003. [Google Scholar] [CrossRef] [PubMed]
  39. Liu, A.X.; Wang, J.J.; Zhang, G.J.; Liu, Y.; Lu, L.N.; Wang, Y.; Zhang, J.X.; Li, X.Y.; Ma, X.Y. Selection of quality control components from Astragali Radix based on multi-component metabolism in rats. Drug Eval. Res. 2018, 41, 216–222. [Google Scholar] [CrossRef]
  40. Li, H.R.; Habasi, M.; Xie, L.Z.; Aisa, H.A. Effect of chlorogenic acid on melanogenesis of B16 melanoma cells. Molecules 2014, 19, 12940–12948. [Google Scholar] [CrossRef] [PubMed]
  41. Jin, N.; Liu, Y.; Xiong, P.Y.; Zhang, Y.Y.; Mo, J.W.; Huang, X.S.; Zhou, Y. Exploring the Underlying Mechanism of Ren-Shen-Bai-Du Powder for Treating Inflammatory Bowel Disease Based on Network Pharmacology and Molecular Docking. Pharmaceuticals 2022, 15, 1038. [Google Scholar] [CrossRef]
  42. Tiwari, P.; Ali, S.A.; Puri, B.; Kumar, A.; Datusalia, A.K. Tinospora cordifolia Miers enhances the immune response in mice immunized with JEV-vaccine: A network pharmacology and experimental approach. Phytomedicine 2023, 119, 154976. [Google Scholar] [CrossRef]
  43. Chen, J.X.; Chen, R.Y.; Lian, Y.; Yao, X.S.; Huang, Y.X.; Gao, Y.X.Y. Progress in Microbial Conversion and Functional Activity of Soy Isoflavones. Shipin Yanjiu Yu Kaifa 2021, 42, 176–182. [Google Scholar] [CrossRef]
  44. Zang, D.; Niu, C.; Aisa, H.A. Amine derivatives of furocoumarin induce melanogenesis by activating Akt/GSK-3β/β-catenin signal pathway. Drug Des. Dev. Ther. 2019, 13, 623–632. [Google Scholar] [CrossRef]
  45. Wan, J.; Lin, F.; Zhang, W.; Xu, A.; DeGiorgis, J.; Lu, H.; Wan, Y. Novel approaches to vitiligo treatment via modulation of mTOR and NF-κB pathways in human skin melanocytes. Int. J. Biol. Sci. 2017, 13, 391–400. [Google Scholar] [CrossRef]
  46. Shin, S.; Ko, J.; Kim, M.; Song, N.; Park, K. Morin Induces Melanogenesis via Activation of MAPK Signaling Pathways in B16F10 Mouse Melanoma Cells. Molecules 2021, 26, 2150. [Google Scholar] [CrossRef]
  47. Liu-Smith, F.; Meyskens, F.L. Molecular mechanisms of flavonoids in melanin synthesis and the potential for the prevention and treatment of melanoma. Mol. Nutr. Food Res. 2016, 60, 1264–1274. [Google Scholar] [CrossRef]
  48. Gu, M.; Jin, J.; Ren, C.; Chen, X.; Gao, W.; Wang, X.; Wu, Y.; Tian, N.; Pan, Z.; Wu, A.; et al. Akebia Saponin D suppresses inflammation in chondrocytes via the NRF2/HO-1/NF-κB axis and ameliorates osteoarthritis in mice. Food Funct. 2020, 11, 10852–10863. [Google Scholar] [CrossRef]
  49. Huang, J.; Chen, X.; Xie, A. Formononetin ameliorates IL-13-induced inflammation and mucus formation in human nasal epithelial cells by activating the SIRT1/Nrf2 signaling pathway. Mol. Med. Rep. 2021, 24, 832. [Google Scholar] [CrossRef]
  50. Sugimoto, M.; Ko, R.; Goshima, H.; Koike, A.; Shibano, M.; Fujimori, K. Formononetin attenuates H2O2-induced cell death through decreasing ROS level by PI3K/Akt-Nrf2-activated antioxidant gene expression and suppressing MAPK-regulated apoptosis in neuronal SH-SY5Y cells. Neurotoxicology 2021, 85, 186–200. [Google Scholar] [CrossRef]
  51. Wang, X.; Song, Z.J.; He, X.; Zhang, R.Q.; Zhang, C.F.; Li, F.; Wang, C.Z.; Yuan, C.S. Antitumor and immunomodulatory activity of genkwanin on colorectal cancer in the APC(Min/+) mice. Int. Immunopharmacol. 2015, 29, 701–707. [Google Scholar] [CrossRef]
  52. Niu, C.; Yin, L.; Aisa, H.A. Novel Furocoumarin Derivatives Stimulate Melanogenesis in B16 Melanoma Cells by Up-Regulation of MITF and TYR Family via Akt/GSK3β/β-Catenin Signaling Pathways. Int. J. Mol. Sci. 2018, 19, 746. [Google Scholar] [CrossRef]
  53. Lei, Z.; Yu, S.; Ding, Y.; Liang, J.; Halifu, Y.; Xiang, F.; Zhang, D.; Wang, H.; Hu, W.; Li, T.; et al. Identification of key genes and pathways involved in vitiligo development based on integrated analysis. Medicine 2020, 99, e21297. [Google Scholar] [CrossRef]
  54. Singh, M.; Mansuri, M.S.; Kadam, A.; Palit, S.P.; Dwivedi, M.; Laddha, N.C.; Begum, R. Tumor Necrosis Factor-alpha affects melanocyte survival and melanin synthesis via multiple pathways in vitiligo. Cytokine 2021, 140, 155432. [Google Scholar] [CrossRef] [PubMed]
  55. Gupta, M.; Sharma, R.; Kumar, A. Docking Techniques in Toxicology: An Overview. Curr. Bioinform. 2020, 15, 600–610. [Google Scholar] [CrossRef]
  56. Saikia, S.; Bordoloi, M. Molecular Docking: Challenges, Advances and its Use in Drug Discovery Perspective. Curr. Drug Targets 2019, 20, 501–521. [Google Scholar] [CrossRef] [PubMed]
  57. Meng, F.C.; Tang, L.D. Challenges and prospect in research of Chinese materia medica network pharmacology. Chin. Tradit. Herb. Drugs 2020, 51, 2232–2237. [Google Scholar] [CrossRef]
Figure 1. Schematic diagram of the study process.
Figure 1. Schematic diagram of the study process.
Pharmaceuticals 16 01716 g001
Figure 2. Statistical analysis based on the fingerprints of 15 batches of raw herbs. (A) Similarity matching of HPLC fingerprints (numbers represent the common peaks); (B) results of the hierarchical cluster analysis (HCA); and (C) principal component analysis (PCA) score scatters (HQ: Huangqi; YZZ: Yuzhizi; CWZ: Chongweizi; JQG: Jinquegen).
Figure 2. Statistical analysis based on the fingerprints of 15 batches of raw herbs. (A) Similarity matching of HPLC fingerprints (numbers represent the common peaks); (B) results of the hierarchical cluster analysis (HCA); and (C) principal component analysis (PCA) score scatters (HQ: Huangqi; YZZ: Yuzhizi; CWZ: Chongweizi; JQG: Jinquegen).
Pharmaceuticals 16 01716 g002aPharmaceuticals 16 01716 g002b
Figure 3. The interaction surface (A) and the contour diagram (B) of the first liquid–solid ratio (X1) and the first extraction time (X2); the interaction surface (C) and the contour diagram (D) of the first liquid–solid ratio (X1) and the second extraction time (X3).
Figure 3. The interaction surface (A) and the contour diagram (B) of the first liquid–solid ratio (X1) and the first extraction time (X2); the interaction surface (C) and the contour diagram (D) of the first liquid–solid ratio (X1) and the second extraction time (X3).
Pharmaceuticals 16 01716 g003
Figure 4. (A) Base peak chromatogram (BPC) of the YQD sample in negative ion mode; (B) base peak chromatogram (BPC) of the YQD sample in positive ion mode; and (C) extracted ion chromatograms (EIC) of absorbed components and metabolites after oral administration of YQD at 1 h.
Figure 4. (A) Base peak chromatogram (BPC) of the YQD sample in negative ion mode; (B) base peak chromatogram (BPC) of the YQD sample in positive ion mode; and (C) extracted ion chromatograms (EIC) of absorbed components and metabolites after oral administration of YQD at 1 h.
Pharmaceuticals 16 01716 g004
Figure 5. MS2 and the fragmentation pathways of (A) calycosin-7-O-β-d-glucoside and (B) abekia saponin.
Figure 5. MS2 and the fragmentation pathways of (A) calycosin-7-O-β-d-glucoside and (B) abekia saponin.
Pharmaceuticals 16 01716 g005
Figure 6. Possible metabolic pathways of (A) calycosin-7-O-β-d-glucoside; (B) methylnissolin-3-O-glucoside; (C) apigenin; (D) genkwanin; and (E) astraisoflavan-7-O-β-d-glucoside.
Figure 6. Possible metabolic pathways of (A) calycosin-7-O-β-d-glucoside; (B) methylnissolin-3-O-glucoside; (C) apigenin; (D) genkwanin; and (E) astraisoflavan-7-O-β-d-glucoside.
Pharmaceuticals 16 01716 g006
Figure 7. (I): EIC diagram of YQD detected after SPE segmentation. (II): HPLC chromatograms of reference (A,B) and the YQD sample (C). (1. chlorogenic acid; 2. calycosin-7-O-β-d-glucoside).
Figure 7. (I): EIC diagram of YQD detected after SPE segmentation. (II): HPLC chromatograms of reference (A,B) and the YQD sample (C). (1. chlorogenic acid; 2. calycosin-7-O-β-d-glucoside).
Pharmaceuticals 16 01716 g007
Figure 8. I: (A) Venn diagram of targets related to YQD and vitiligo. (B) Visualization analysis of the PPI diagram—the size of the circle represents the target degree. (C) Adjacent nodes of the overlapping targets between YQD and vitiligo—the X-axis represents the number of adjacent proteins of the target and the Y-axis represents different targets. (D) Compound–target–disease (CTD) network—31 components (blue) were connected to 78 targets (green) and YQD (red), and 78 target genes were connected to 31 components and vitiligo (red). II: (A) GO analysis. (B) Top 20 enriched KEGG pathways. (C) Compound–target–pathway (CTP) network—31 components (blue) were connected to 41 targets (orange), and the 41 targets were connected to 20 pathways (yellow) and 31 components.
Figure 8. I: (A) Venn diagram of targets related to YQD and vitiligo. (B) Visualization analysis of the PPI diagram—the size of the circle represents the target degree. (C) Adjacent nodes of the overlapping targets between YQD and vitiligo—the X-axis represents the number of adjacent proteins of the target and the Y-axis represents different targets. (D) Compound–target–disease (CTD) network—31 components (blue) were connected to 78 targets (green) and YQD (red), and 78 target genes were connected to 31 components and vitiligo (red). II: (A) GO analysis. (B) Top 20 enriched KEGG pathways. (C) Compound–target–pathway (CTP) network—31 components (blue) were connected to 41 targets (orange), and the 41 targets were connected to 20 pathways (yellow) and 31 components.
Pharmaceuticals 16 01716 g008aPharmaceuticals 16 01716 g008b
Figure 9. (A) Heatmap of the predicted molecular docking of the predicted active compounds with three core targets; (B) M6 and AKT1; (C) M30 and AKT1; (D) M15 and MAPK1; (E) M10 and MAPK1; (F) C31 and mTOR; and (G) P4 and mTOR.
Figure 9. (A) Heatmap of the predicted molecular docking of the predicted active compounds with three core targets; (B) M6 and AKT1; (C) M30 and AKT1; (D) M15 and MAPK1; (E) M10 and MAPK1; (F) C31 and mTOR; and (G) P4 and mTOR.
Pharmaceuticals 16 01716 g009
Figure 10. YQD promotes melanogenesis in B16F10 cells. The B16F10 cells were treated with YQD at different concentrations (1–1000 µg/mL), medium alone (control), or 8-MOP (positive control, 100 μmol/L) for 48 h: (A) effect on cell viability examined using the CCK-8 assay; (B) melanin content of cells measured by NaOH assay; (C) tyrosinase activity measured by tyrosinase activity assay; (D) expression of the TRP-1 gene analyzed via qRT-PCR; and (E) levels of TRP-1 proteins evaluated using Western blotting. (ns: no significance, * p < 0.05, ** p < 0.01, and *** p < 0.001).
Figure 10. YQD promotes melanogenesis in B16F10 cells. The B16F10 cells were treated with YQD at different concentrations (1–1000 µg/mL), medium alone (control), or 8-MOP (positive control, 100 μmol/L) for 48 h: (A) effect on cell viability examined using the CCK-8 assay; (B) melanin content of cells measured by NaOH assay; (C) tyrosinase activity measured by tyrosinase activity assay; (D) expression of the TRP-1 gene analyzed via qRT-PCR; and (E) levels of TRP-1 proteins evaluated using Western blotting. (ns: no significance, * p < 0.05, ** p < 0.01, and *** p < 0.001).
Pharmaceuticals 16 01716 g010
Figure 11. YQD promotes melanogenesis in zebrafish. The zebrafish embryos were exposed to E3 buffer (control), varying doses of YQD (200–1200 μg/mL), or 8-MOP (positive control, 100 μmol/L) from 24 hpf (hours postfertilization) to 84 hpf. The melanin densities (IOD) in the tails were measured using Image-Pro Plus: (A) melanin in the tails of zebrafish 84 hpf; (B) melanin densities in zebrafish. (ns: no significance, * p < 0.05, *** p < 0.001).
Figure 11. YQD promotes melanogenesis in zebrafish. The zebrafish embryos were exposed to E3 buffer (control), varying doses of YQD (200–1200 μg/mL), or 8-MOP (positive control, 100 μmol/L) from 24 hpf (hours postfertilization) to 84 hpf. The melanin densities (IOD) in the tails were measured using Image-Pro Plus: (A) melanin in the tails of zebrafish 84 hpf; (B) melanin densities in zebrafish. (ns: no significance, * p < 0.05, *** p < 0.001).
Pharmaceuticals 16 01716 g011
Table 1. Rank of principal component scores and comprehensive scores for HQ, YZZ, CWZ, and JQG.
Table 1. Rank of principal component scores and comprehensive scores for HQ, YZZ, CWZ, and JQG.
Sample
(HQ/CWZ)
PC1
Scores
PC2
Scores
PC3
Scores
Comprehensive ScoresRankSample
(YZZ/JQG)
PC1
Scores
PC2
Scores
PC3
Scores
Comprehensive ScoresRank
H-1−0.960.34−1.04−0.5812Y-12.87−0.28-1.922
H-2−0.900.390.81−0.297Y-2−0.281.91-0.138
H-30.041.480.920.396Y-32.34−0.59-1.505
H-4−0.87−0.09−0.37−0.519Y-4−2.66−0.25-−1.8713
H-5−1.550.351.45−0.5410Y-5−2.11−1.06-−1.6312
H-60.331.171.380.555Y-61.822.11-1.604
H-7−1.06−0.710.17−0.6413Y-7−1.94−0.52-−1.4211
H-8−0.84−0.69−0.17−0.5711Y-8−1.61−1.53-−1.3610
H-9−1.100.310.31−0.478Y-9−1.66−0.51-−1.229
H-101.56−0.622.130.994Y-100.232.18-0.537
H-114.63−3.23−0.231.832Y-115.09−0.56-3.401
H-123.871.04−2.001.901Y-121.200.50-0.916
H-13−2.161.14−2.09−1.2114Y-13−2.95−1.10-−2.2115
H-143.201.40−0.391.832Y-14−3.511.34-−2.1814
H-15−4.22−2.29−0.88−2.6715Y-153.16−1.64-1.893
C-11.81--1.332J-1−1.361.620.81−0.536
C-20.64--0.477J-2−1.57−0.23−0.87−1.1215
C-3−0.83--−0.6113J-3−1.242.071.09−0.365
C-40.17--0.1310J-4−1.170.410.20−0.668
C-51.00--0.743J-5−1.200.00−0.50−0.8111
C-60.94--0.704J-6−1.420.84−1.01−0.8813
C-7−0.34--−0.2511J-7−1.15−0.660.12−0.819
C-80.41--0.308J-8−1.18−1.160.79−0.8212
C-90.81--0.606J-9−0.71−1.36−0.60−0.7110
C-100.28--0.219J-104.94−0.39−0.213.032
C-11−5.01--−3.7015J-115.54−0.291.243.581
C-120.86--0.645J-12−0.55−1.78−0.21−0.637
C-131.87--1.381J-13−1.17−0.69−0.62−0.9014
C-14−2.08--−1.5414J-14−0.62−0.041.76−0.204
C-15−0.55--−0.4012J-152.851.66−1.981.813
Table 2. Three-factor Box–Behnken design (BBD) with factorial levels shown in actual values for the YQD.
Table 2. Three-factor Box–Behnken design (BBD) with factorial levels shown in actual values for the YQD.
RunIndependent VariablesResponsesOD(Y)
X1X2X3CAGCCACTSR
18(0)60(1)40(1)0.1070.3850.38500.968
28(0)45(0)30(0)0.0870.3410.36340.887
39(1)45(0)40(1)0.1100.4140.39020.988
49(1)30(−1)30(0)0.0980.3320.39300.937
58(0)45(0)30(0)0.0910.3220.38160.907
67(−1)30(−1)30(0)0.0650.2870.33420.000
77(−1)45(0)40(1)0.0760.2450.33300.633
88(0)60(1)20(−1)0.0810.3010.37120.856
99(1)45(0)20(−1)0.0980.3370.38140.929
108(0)45(0)30(0)0.0860.2640.37220.825
118(0)30(−1)20(−1)0.0710.2650.33340.670
129(1)60(1)30(0)0.0990.3710.40420.967
138(0)45(0)30(0)0.0810.3190.34180.809
148(0)30(−1)40(1)0.0780.3430.36540.864
157(−1)60(1)30(0)0.0690.2540.33120.606
168(0)45(0)30(0)0.0860.3170.36880.878
177(−1)45(0)20(−1)0.0870.2390.32920.000
X1: first liquid–solid ratio (mL/g); X2: first extraction time (min); X3: second extraction time (min). CAGC: calycosin-7-O-β-d-glucoside content (mg/g); CAC: chlorogenic acid content (mg/g); TSR: total solid ratio (g); OD: overall desirability.
Table 3. Characterization of chemical constituents of YQD by UPLC-Q-TOF/MS.
Table 3. Characterization of chemical constituents of YQD by UPLC-Q-TOF/MS.
IDRT (Min)AdductionMeasured MassCalculated MassppmFormulaNameMS/MS (m/z)SourceType
C10.99[M + H]+118.0858118.0863−4.23C5H11NO2betaine118.0862; 58.0644HQalkaloid
C21.16[M + H]+144.1014144.1019−3.47C7H13NO2stachydrine144.1014; 84.0801; 58.0645CWZalkaloid
C318.55[M − H]431.1208431.11953.02C18H24O12asperulosidic acid431.1207; 137.0244; 93.0354HQiridoid glycosides
C418.96[M − H]313.0923313.0929−1.92C14H18O8glucovanillin151.0409; 123.0446HQglycosides
C519.49[M − H]353.0890353.08783.40C16H18O9neochlorogenic acid353.0900; 191.0574; 179.0356; 135.0459CWZorganic acids
C620.56[M − H]593.2094593.20871.18C25H38O16leonuriside B593.2054; 461.1655; 135.0437CWZphenylethanoid glycosides
C723.06[M − H]353.0884353.08781.70C16H18O9chlorogenic acid191.0562; 85.0292CWZorganic acids
C823.74[M − H]353.0881353.08780.85C16H18O9cryptochlorogenic acid353.0915; 191.0567; 179.0364; 173.0469; 135.0455CWZorganic acids
C925.53[M − H]463.1254463.12461.73C22H24O11hesperetin 7-O-glucoside463.1243; 301.0800; 257.0800CWZflavonoids-glycosides
C1028.62[M + H]+312.1540312.1554−4.48C14H21N3O5leonurine312.1588; 181.0482; 114.0997CWZalkaloid
C1130.41[M − H]301.0715301.0718−1.00C16H14O6carasinaurone301.0723; 257.0778; 139.00404; 124.0175JQGflavonoids
C1230.86[M + H]+447.1269447.1286−3.80C22H22O10calycosin-7-O-β-d-glucoside285.0750; 270.0517; 253.0507; 242.0541; 225.0522; 137.0217HQflavonoids-glycosides
C1333.14[M − H]477.1419477.14023.56C23H26O11calceolarioside B477.1409; 315.1073; 161.0253; 133.0296YZZphenylethanoid glycosides
C1434.08[M − H]515.1211515.11953.11C25H24O12isochlorogenic acid B515.1205; 353.0871; 191.0561; 173.0452; 135.0451CWZorganic acids
C1534.81[M − H]515.1212515.11953.30C25H24O12isochlorogenic acid A515.1207; 353.0871; 191.0565; 179.0363; 135.0456CWZorganic acids
C1635.70[M − H]461.1461461.14531.73C23H26O10methylnissolin-3-O-glucoside461.1452; 179.0371; 161.0245; 133.0303HQflavonoids-glycosides
C1736.68[M − H]515.1202515.11951.36C25H24O12isochlorogenic acid C515.1195; 353.0877; 191.0558; 173.0460; 135.0456CWZorganic acids
C1838.03[M + FA − H]475.1247475.12460.21C22H22O9ononin267.0691; 252.0458CWZ/HQflavonoids-glycosides
C1939.17[M + H]+1075.53051075.532−1.21C52H82O23mutongsaponin B943.4894; 797.4332; 635.3815;YZZtriterpenoid saponins
C2040.41[M + FA − H]987.4837987.48063.14C47H74O19saponin Ph987.4972; 941.4768; 469.1586YZZtriterpenoid saponins
C2142.03[M − H]283.0617283.06121.77C16H12O5genkwanin283.0595; 268.0382; 211.0388; 195.0456CWZflavonoids
C2242.08[M − H]463.1608463.161−0.43C23H28O10astraisoflavan-7-O-β-d-glucoside463.1615; 301.1061; 179.0690HQflavonoids-glycosides
C2342.61[M + FA − H]1003.51161003.5119−0.30C48H78O19scheffoleoside A1003.5117; 957.4984; 487.3411; 469.1534YZZtriterpenoid saponins
C2442.79[M + H]+1091.56301091.5633−0.27C53H86O23akebia saponin959.5186; 813.4651; 651.4084; 489.3562YZZtriterpenoid saponins
C2543.60[M − H]431.1001431.09843.94C21H20O10apigenin 7-O-glucoside431.0965; 269.0419HQflavonoids-glycosides
C2643.96[M + FA − H]1003.51091003.5119−1.00C48H78O19asiaticoside1003.5229; 957.5041; 487.3487; 469.1570YZZtriterpenoid saponins
C2744.10[M + H]+1091.55791091.5633−4.95C53H86O2323-hydroxyakemisaponin A959.5174; 813.4715; 651.3946; 489.3537YZZtriterpenoid saponins
C2844.73[M + FA − H]1003.51541003.51193.49C48H78O19(2α,3β,6β)-2,3,6-trihydroxy-Olean-12-en-28-oic acid-O-6-deoxy-α-L-mannopyranosyl-(1→4)-O-β-d-glucopyranosyl-(1→6)-β-d-glucopyranosyl ester1003.5123; 957.5028; 487.3442; 469.1551YZZtriterpenoid saponins
C2945.53[M + FA − H]969.2802969.27643.92C56H44O13kobophenol A923.2608; 801.2261JQGdiphenylethenes
C3049.65[M + FA − H]1105.54491105.54361.18C52H84O22akebia saponin E1105.5482; 1059.5363; 735.4314YZZtriterpenoid saponins
C3150.55[M + FA − H]973.5012973.5014−0.21C47H76O18akebia saponin D973.4997; 927.4915; 603.3873; 323.0960; 179.0557YZZtriterpenoid saponins
C3251.88[M − H]267.0656267.0663−2.62C16H12O4formononetin267.0668; 252.0436; 195.0444; 135.0122CWZ/HQflavonoids
C3355.19[M + Cl]933.4640933.46202.14C46H74O17(3β,4α)-3-[(O-α-L-arabinopyranosyl-(1→2)-O-[β-d-glucopyranosyl-(1→4)]-α-L-arabinopyranosyl)oxy]-23-hydroxy-Olean-12-en-28-oic acid933.4631; 735.4207YZZtriterpenoid saponins
C3456.36[M + FA − H]1089.55231089.54873.30C52H84O21yuzhizioside IV1089.5504; 1043.5381; 719.4348; 323.0966YZZtriterpenoid saponins
C3559.25[M + FA − H]871.4718871.46972.41 C43H70O15astragaloside II871.4730; 825.4697; 765.4370HQtriterpenoid saponins
C3661.63[M + FA − H]871.4689871.4697−0.92 C43H70O15astrasieversianin VII871.4678; 825.4468HQtriterpenoid saponins
C3763.61[M + FA − H]871.4731871.46973.90 C43H70O15cyclosiversioside D871.4665; 825.4716; 765.4327HQtriterpenoid saponins
C3863.95[M + FA − H]957.5079957.5123−4.60 C47H76O17pulsatilla Saponin D911.4976; 749.4245YZZtriterpenoid saponins
C3964.27[M − H]487.3428487.3429−0.21 C30H48O5arjunolic acid487.3422; 363.0096; 322.9947YZZtriterpenoids
C4064.84[M + FA − H]781.4380781.43800.00 C40H64O12akebia saponin B781.4406; 735.4318; 603.3912YZZtriterpenoid saponins
C4165.12[M + FA − H]913.4847913.48024.93 C45H72O16astragaloside I913.4814; 867.4702; 765.4267HQtriterpenoid saponins
C4265.53[M − H]269.0457269.04550.74 C15H10O5apigenin269.0430; 225.0556HQflavonoids
C4365.71[M + FA − H]913.4814913.48021.31 C45H72O16isoastragaloside I913.4861; 867.4913HQtriterpenoid saponins
C4466.41[M + FA − H]913.4811913.48020.99 C45H72O16neoastragaloside I/HQtriterpenoid saponins
Table 4. Identification of absorbed prototype components and metabolites in rat plasma by UPLC-Q-TOF/MS (1 h after the last administration).
Table 4. Identification of absorbed prototype components and metabolites in rat plasma by UPLC-Q-TOF/MS (1 h after the last administration).
IDRT (Min)AdductionMeasured MassCalculated MassppmFormulaNameMS/MS (m/z)
P1(C2)1.4[M + H]+144.1015144.1019−2.78C7H13NO2stachydrine144.1019; 84.0803; 58.0649
P2(C11)30.8[M − H]301.0727301.07182.99C16H14O6carasinaurone273.0734; 241.0522; 193.09778; 139.0410; 124.0178
P3(C12)31.58[M − H]447.1281447.1286−1.12C22H22O10calycosin-7-O-β-d-glucoside447.1179; 285.0760; 270.0554; 253.0519
P4(C40)65.48[M + FA − H]781.4405781.43803.20C40H64O12akebia saponin B781.4385; 735.4349; 603.4002; 471.3462
M125.78[M − H]509.0412509.03953.34C21H18O13Sdaidzein + sulfation + glucuronidation509.0401; 333.0079; 253.0511
M227.92[M − H]621.1487621.14614.19C28H30O16calycosin-7-O-β-d-glucoside + glucuronidation621.1492; 459.0964; 283.0620; 268.0344
M328.14[M − H]509.0415509.03953.93C21H18O13Sdaidzein + sulfation + glucuronidation509.0419; 429.0845; 333.0045; 253.0527
M428.81[M − H]429.0836429.08272.10C21H18O10daidzein + glucuronidation253.0516; 175.0232
M528.81[M − H]525.0359525.03452.67C21H18O14Sapigenin + sulfation + glucuronidation525.0370; 349.0028; 269.0465
M629.71[M − H]459.0946459.09332.83C22H20O11genkwanin + glucuronidation459.0964; 283.0616; 268.0389; 240.0447
M730.75[M − H]497.0775497.07593.22C21H22O12Sequol + sulfation + glucuronidation497.0759; 417.1175; 321.0430; 241.0857; 175.0232
M831.56[M − H]429.0843429.08273.73C21H18O10daidzein + glucuronidation429.0854; 253.0514; 129.0216; 113.0250
M931.75[M − H]541.0301541.02941.29C21H18O15Sapigenin + hydroxylation + glucuronidation + sulfation541.0318; 461.0719; 285.0413
M1032.96[M − H]445.0794445.07764.04C21H18O11apigenin + glucuronidation445.0820; 269.0482; 224.0483; 175.0253; 113.0251
M1133.69[M − H]433.1158433.1144.16C21H22O102′,4′,7-trihydroxyisoflavan + glucuronidation433.12775; 296.9257; 257.0741; 175.0232; 135.0435
M1233.65[M − H]491.1575491.15593.26C24H28O117-hydroxy-2′,3′,4′-trimethoxyisoflavan + glucuronidation491.1598; 315.1253; 271.1348
M1334.78[M − H]445.0798445.07764.94C21H18O11apigenin + glucuronidation445.2026; 269.0479; 180.0724
M1434.84[M − H]333.0089333.00744.50C15H10O7Sdaidzein + sulfation333.0095; 253.0519; 224.0475; 117.0347
M1534.88[M − H]447.0940447.09331.57C21H20O117,3′,4′-trihydroxyisoflavanone + glucuronidation447.0942; 271.0625; 198.9095; 165.0206
M1635.55[M − H]459.0950459.09333.70C22H20O11calycosin + glucuronidation459.0962; 283.0612; 268.0391; 175.0251
M1735.55[M − H]417.1207417.11913.84C21H22O9equol + glucuronidation417.1223; 399.1180; 241.0886; 175.0258
M1835.91[M − H]447.0943447.09332.24C21H20O117,3′,4′-trihydroxyisoflavanone + glucuronidation447.1322; 271.0639; 151.0070; 113.0255
M1936.18[M − H]417.1210417.11914.56C21H22O9equol + glucuronidation417.1217; 399.1116; 241.0871; 175.0269
M2036.11[M − H]639.1957639.19314.07C29H36O16astraisoflavan-7-O-β-d-glucoside + glucuronidation639.1943; 463.1845; 301.1089; 286.0901
M2138.58[M − H]417.1207417.11913.84C21H22O9equol + glucuronidation417.1214; 241.0878; 175.0247; 121.0303
M2238.63[M − H]321.0449321.04383.43C15H14O6Sequol + sulfation321.0489; 242.0913; 121.0305
M2338.92[M − H]443.0999443.09843.39C22H20O10formononetin + glucuronidation443.1701; 267.0670; 252.0441; 223.0398; 179.0250
M2440.04[M − H]253.0511253.05061.98C15H10O4formononetin + demethylation253.0518; 225.0558; 135.0094
M2540.86[M − H]433.1153433.1143.00C21H22O102′,4′,7-trihydroxyisoflavan + glucuronidation433.2095; 257.0832; 175.0253; 113.0252
M2641.42[M − H]433.1157433.1143.93C21H22O102′,4′,7-trihydroxyisoflavan + glucuronidation433.1166; 257.0826; 136.0182
M2741.99[M − H]475.1267475.12464.42C23H24O11methylnissolin + glucuronidation475.1302; 301.1060; 284.0711; 269.0862
M2842.76[M − H]459.0951459.09333.92C22H20O11calycosin + glucuronidation459.0979; 283.0662; 268.0368; 255.0686
M2943.13[M − H]477.1421477.14023.98C23H26O11astraisoflavan + glucuronidation477.1445; 301.1100; 286.0859; 271.0616
M3043.28[M + H]+303.1212303.1227−4.95C17H18O5astraisoflavan303.1024; 167.0699; 133.0677; 123.0434
M3144.51[M − H]445.1154445.1143.15C22H22O10formononetin + hydrogenation + glucuronidation445.1165; 269.0850; 254.0608; 217.0108; 113.0268
M3245.68[M − H]347.0230347.0231−0.29C16H12O7Sformononetin + sulfation347.0279; 267.0650; 252.0431; 223.0376
Table 5. Information table on active components of YQD.
Table 5. Information table on active components of YQD.
IDNameFormulaMW
C13Calceolarioside BC23H26O11478.15
C24Akebia saponin PJ1C53H86O231090.56
C30Akebia saponin EC52H84O221060.54
C31Akebia saponin DC47H76O18928.50
C42ApigeninC15H10O5270.05
P1StachydrineC7H13NO2143.09
P2CarasinauroneC16H14O6302.08
P3Calycosin-7-O-β-d-glucosideC22H22O10446.12
P4Akebia saponin BC40H64O12736.44
M2Calycosin-7-O-β-d-glucoside+glucuronidationC28H30O16622.15
M5Apigenin+sulfation+glucuronidationC21H18O14S526.04
M6Genkwanin+glucuronidationC22H20O11460.10
M9Apigenin+hydroxylation+glucuronidation+sulfationC21H18O15S542.04
M10Apigenin+glucuronidationC21H18O11446.08
M112′,4′,7-trihydroxyisoflavan+glucuronidationC21H22O10434.12
M127-hydroxy-2′,3′,4′-trimethoxyisoflavan+glucuronidationC24H28O11492.16
M13Apigenin+glucuronidationC21H18O11446.08
M157,3′,4′-trihydroxyisoflavanone+glucuronidationC21H20O11448.10
M16Calycosin+glucuronidationC22H20O11460.10
M187,3′,4′-trihydroxyisoflavanone+glucuronidationC21H20O11448.10
M20Astraisoflavan-7-O-β-d-glucoside+glucuronidationC29H36O16640.20
M23Formononetin+glucuronidationC22H20O10444.11
M24Formononetin+demethylationC15H10O4254.06
M252′,4′,7-trihydroxyisoflavan+glucuronidationC21H22O10434.12
M262′,4′,7-trihydroxyisoflavan+glucuronidationC21H22O10434.12
M27Methylnissolin+glucuronidationC23H24O11476.13
M28Calycosin+glucuronidationC22H20O11460.10
M29Astraisoflavan+glucuronidationC23H26O11478.15
M30AstraisoflavanC17H18O5302.12
M31Formononetin+hydrogenation+glucuronidationC22H22O10446.12
M32Formononetin+sulfationC16H12O7S348.03
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Cui, L.; Ma, C.; Shi, W.; Yang, C.; Wu, J.; Wu, Z.; Lou, Y.; Fan, G. A Systematic Study of Yiqi Qubai Standard Decoction for Treating Vitiligo Based on UPLC-Q-TOF/MS Combined with Chemometrics, Molecular Docking, and Cellular and Zebrafish Assays. Pharmaceuticals 2023, 16, 1716. https://doi.org/10.3390/ph16121716

AMA Style

Cui L, Ma C, Shi W, Yang C, Wu J, Wu Z, Lou Y, Fan G. A Systematic Study of Yiqi Qubai Standard Decoction for Treating Vitiligo Based on UPLC-Q-TOF/MS Combined with Chemometrics, Molecular Docking, and Cellular and Zebrafish Assays. Pharmaceuticals. 2023; 16(12):1716. https://doi.org/10.3390/ph16121716

Chicago/Turabian Style

Cui, Lijun, Cui Ma, Wenqing Shi, Chen Yang, Jiangping Wu, Zhenghua Wu, Yuefen Lou, and Guorong Fan. 2023. "A Systematic Study of Yiqi Qubai Standard Decoction for Treating Vitiligo Based on UPLC-Q-TOF/MS Combined with Chemometrics, Molecular Docking, and Cellular and Zebrafish Assays" Pharmaceuticals 16, no. 12: 1716. https://doi.org/10.3390/ph16121716

APA Style

Cui, L., Ma, C., Shi, W., Yang, C., Wu, J., Wu, Z., Lou, Y., & Fan, G. (2023). A Systematic Study of Yiqi Qubai Standard Decoction for Treating Vitiligo Based on UPLC-Q-TOF/MS Combined with Chemometrics, Molecular Docking, and Cellular and Zebrafish Assays. Pharmaceuticals, 16(12), 1716. https://doi.org/10.3390/ph16121716

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop