Discovery And Identification Of Potential Anti-melanogenic Active Constituents Of Bletilla Striata

ali.ma@wecistanche.com

Yiyuan Luo1, Juan Wang1, Shuo Li1, Yue Wu1, Zhirui Wang1, Shaojun Chen1 and Hongjiang Chen1,2*

Abstract

Background: Bletilla striata is the main medicine of many skin whitening classic formulas in traditional Chinese medicine (TCM) and is widely used in cosmetic industry recently. However, its active ingredients are still unclear and its fibrous roots are not used effectively. The aim of the present study is to discover and identify its potential anti-melanogenic active constituents by zebrafish model and molecular docking.

Methods: The antioxidant activities were evaluated by 2,2-diphenyl-1-picrylhydrazyl (DPPH) radical scavenging activity, 2,2′-azino-bis-(3-ethylbenthiazoline-6-sulphonic acid) (ABTS) radical scavenging activity and ferric reducingantioxidant power (FRAP) assay. The anti-melanogenic activity was assessed by tyrosinase inhibitory activity in vitro and melanin inhibitory in zebrafish. The chemical profiles were performed by ultra-high-performance liquid chromatography combined with quadrupole time-of-flight tandem mass spectrometry (UPLC-Q-TOF-MS/MS). Meanwhile, the potential anti-melanogenic active constituents were temporary identified by molecular docking.

Results: The 95% ethanol extract of B. striata fibrous roots (EFB) possessed the strongest DPPH, ABTS, FRAP and tyrosinase inhibitory activities, with IC50 5.94 mg/L, 11.69 mg/L, 6.92 mmol FeSO4/g, and 58.92 mg/L, respectively. In addition, EFB and 95% ethanol extract of B. striata tuber (ETB) significantly reduced the melanin synthesis of zebrafish embryos in a dose-dependent manner. 39 chemical compositions, including 24 stilbenoids were tentatively identified from EFB and ETB. Molecular docking indicated that there were 83 (including 60 stilbenoids) and 85 (including 70 stilbenoids) compounds exhibited stronger binding affinities toward tyrosinase and adenylate cyclase.

Conclusion: The present findings supported the rationale for the use of EFB and ETB as natural skin-whitening agents in pharmaceutical and cosmetic industries.

Keywords: Bletilla striata, Antioxidant, Anti-melanogenic activity, UPLC-Q-TOF-MS/MS, Zebrafish, Molecular docking

Click to cistanche side effects and benefits for antioxidant

Background

Bletilla striata (Thunb.) Reichb. f. is a herbaceous perennial plant widely distributed in Asia, such as China, Korea, and Japan [1]. The dried tuber of B. striata, also known as Baiji, firstly recorded in Shennong’s Classic of Materia Medica, has been widely used as a traditional Chinese medicine (TCM) for thousands of years in China. Chinese pharmacopeias states that it possesses the capability of astringency upon hemostasis and analgesis, therefore, it was widely used for the treatment of hematemesis, hemoptysis, traumatic bleeding, ulcers, swelling and chapped skin [2, 3]. Pharmacological studies demonstrated that B. striata possessed a wide spectrum of biological activities, such as wound healing [4, 5], antiulcer [6, 7], hemostatic [8], anti-inflammation [9], antioxidant [10], antibacterial [11], anti-influenza viral [12], and antiaging [13]. At the same time, B. striata contains various classes of chemical compositions, including polysaccharides [14], bibenzyls, phenanthrenes, anthraquinones,flavonoids, and 2-isobutylmalates [3, 15], etc.

B. striata is the main medicine of many skin whitening classic formulas in TCM [16] and it is widely used in cosmetic industry recently [17]. However, its active ingredients are still unclear and even some research results were contrary. For instance, the research results given by Chen et al. indicated that 95% ethanol extract of B. striata possessed higher tyrosinase inhibitory activity than those of its water extract with inhibition rate of 68.36% in vitro [18], while the result of Huang et al. showed that the inhibitory activity of B. striata water extract on tyrosinase was stronger than that of 95% ethanol extract with inhibition rate of 62% in vitro [19]. The research results byLuetal.[20]andLinghuetal.[21]demonstratedthat both water and 95% ethanol extract of B. striata, especially its chloroform fractionation, could inhibit B16 cell growth and induce its apoptosis in a concentration dependent manner.

As the in vitro tests of tyrosinase inhibition and melanoma cell lines, didn’t involve the complex physiological in vivo conditions, and the absorption, metabolism, distribution and excretion of the test sample, many of the samples showed significant inhibitory activity against tyrosinase and melanoma cell lines in vitro, nevertheless, exhibited weakly effective or even ineffective in vivo [22,23]. The glabridin showed significant tyrosinase inhibitory activity with IC50 0.43μmol/L, which was 176 times stronger than that of kojic acid. Nevertheless, it was no inhibitory effects on the pigmentation of zebrafish in vivo [24].

Meanwhile, B. striata fibrous roots (FB) were the byproducts generated during the processing of B. striatatuber (TB). Modern researches indicated that the FB contains similar compounds with TB and with higher content of phenolic [25]. In addition, the antibacterial [26], antioxidant and anti-tyrosinase activities of FB extract were stronger than those of TB extract [25]. However, the resources of FB were not effectively used and were abandoned in the farmland, which led to the waste of FB resources and environmental pollution [27].

Zebrafish (Danio rerio), a small tropical freshwater fish, is an emerging animal model for pharmacology

and toxicology research in vivo, with many advantages including low cost, short life cycle, high transparency, easy to maintain, high fertility, and required fewer test samples [28, 29]. In addition, zebrafish has melanin pigments on the surface, allowing simple observation of the pigmentation process, and was widely used in the anti-melanogenic study [28, 30].

In the present study, we compared the antioxidant and tyrosinase inhibitory activities of crude polysaccha- ride and 95% ethanol extract of TB and FB in vitro, and anti-melanogenic activities in zebrafish model. Mean- while, the potential anti-melanogenic active constituents were temporary identified by molecular docking. We discovered that 95% ethanol extract of TB (ETB) and FB (EFB) possesses the significantly antioxidant and anti-tyrosinase activities, and can reduce melanin synthesis of zebrafish embryos in a dose-dependent manner. Thus, ETB and EFB can be used as natural skin-whitening agents in pharmaceutical and cosmetic industries.

Methods

Chemicals and reagents

Tyrosinase, 3-(3,4-dihydroxyphenyl)-L-alanine (L-DOPA) and arbutin were purchased from Aladdin reagent company (Shanghai, China). 2,2-diphenyl-1-picrylhydrazyl (DPPH), 2,2′-azino- bis-(3-ethylbenthiazoline-6-sulphonic acid) (ABTS), 6-hydroxy-2,5,7,8-tetramethylchroman-2-carboxylic acid (Trolox), and 2,4,6-tris (2-pyridyl)-s-triazine (TPTZ) were purchased from Sigma-Aldrich (St. Louis, MO, USA). Acetonitrile and methanol (HPLC grade) were purchased from Tedia (Fairfield, OH, USA). Deionized water was prepared by the Milli-Q water system (18.2 MΩ, Millipore, USA).

Plant collection and extraction procedure

B. striata was collected from Quzhou YiNianTang Agriculture and Forestry Technology Co., Ltd. (Quzhou, Zhejiang, China). Voucher specimens were deposited at Herbarium of Zhejiang Pharmaceutical College (Accession no. 190615). The whole plant was divided into tubers and fibrous roots. Subsequently, they were cut into small pieces and dried by vacuum freeze-drying, respectively.

The dried and powdered samples (TB and FB) were reflux extracted with 95% ethanol for three times (1.5h for each time). The extract was filtered using a What man filter paper (No.1). The filtrate was concentrated to dryness in a rotary evaporator (Hei-VAP, Heidolph, Germany) at 50°C and was subseqently dried by a vacuum freeze dryer. The ethanol extraction of TB (ETB) and FB (EFB) yields were 5.21 and 6.53%, respectively. The dregs were used to extract the crude polysaccharide by dispersed in 80°C water for 4h and precipitation with ethanol [31]. The polysaccharide yields of TB (PTB) and FB (PFB) were 14.75 and 6.45%, respectively.

cistanche extract

Chemical analysis

The chemical analysis was performed on an ultra-high-performance liquid chromatography combined with quadrupole time-of-flight tandem mass spectrometry(UPLC-Q-TOF-MS/MS). The chromatographic separation was performed on a Waters Acquity UPLC tem (Waters Corp., Milford, MA, US) with a Waters BEH Shield RP C18 column (100×2.1mm, 1.7μm) at 30°C. The mobile phase was consisted of 0.1% formic acid in water (A) and acetonitrile (B), with a linear gradient: 0–3min, 5–16% B; 3–8min, 16–30% B; 8–10min, 30–35% B; 10–15min, 35–55% B; 15–18min, 55–80% B; 18–19 min, 80–5% B; 19–20 min, 5% B. The flow rate was 0.3 mL/min, and the injection volume was 2 μL.

The TOF-MS/MS experiments were performed by using an AB SCEIX Triple TOF 5600 mass spectrometry (AB SCIEX, Foster City, CA, USA) equipped with an electrospray ionization (ESI). The MS detection was con- ducted in negative ionization method. The parameters were set as follows: source and desolvation temperature were 100°C and 450°C respectively; desolvation gas flow rate was 900L/h; capillary voltage was 2kV; cone voltage was 40 V; collision energy was 22 eV; and the full scan spectra was from 100 to 2000 Da.

Antioxidant assays DPPH radical scavenging activity

The DPPH radical scavenging activity was determined in accordance with Ali et al. with slightly modifications [32]. Briefly, 50 μL test sample solution was mixed with 150 μL DPPH ethanol solution (0.2mM) in 96-well plates and was kept in the dark for 30min. The absorbance of the mixture at 517nm (A1) was assessed by microplate spectrophotometer (Thermo Fisher Scientific, America). The blank solution without test sample mixed with DPPH ethanol solution was set as the positive group (A0). The blank solution without DPPH mixed with test sample solution was set as the blank group (A2). All experiments were executed in triplicate. The DPPH radical scavenging rate was calculated using the following formula:

DPPH radical scavenging rate (%) = [A0-(A1-A2)]/ A0 × 100%.

ABTS radical scavenging activity

The ABTS radical scavenging capacity was measured using the method described by Ali et al. [32]. Briefly, 7mM ABST aqueous solution was mixed with 2.45mM potassium persulfate at a ratio of 1:1 (v/v), and incubating at room temperature in the dark for 16 h to obtain the ABST+ stock solution. The stock solution was diluted with phosphate buffered saline (PBS) to adjust the absorbance of(0.72±0.2) at 734nm. 20μL test samples at different concentrations were mixed thoroughly with 180μL ABTS+ solution. The reactive mixtures were kept in the dark for 5min and subsequently measured the absorbance (B1) at 734nm. The absorbance of the mixtures without test sample and ABST+ were set as B and B , respectively. All experiments were executed in triplicate. The ABTS radical scavenging rate was calculated according to the following formula:

ABTS radical scavenging rate (%) = [B0-(B1-B2)]/B0 × 100%.

Ferric reducing antioxidant power assay (FRAP)

The ferric iron reducing activity was determined according the procedures by Kosakowska et al. [33]. 300mM acetate buffer, 10mM TPTZ (2,4,6-tris (2-pyridyl)-s-triazine) and 20mM FeCl3 were mixed at a (v/v/v) ratio of 10:1:1 to prepared the working reagent. 100μL of each test sample solution was mixed with 100μL TPTZ working reagent, incubated at 37 °C for 20 min and the absorbance was determined at 593 nm. Meanwhile, a series of FeSO4 standard solutions in the concentration ranges of 0–1000μg/ mL were used to prepare the calibration curve. The ferric reducing capacity was calculated in the formation of an intense Fe2+-TPTZ blue complex. The results were expressed as Fe2+ antioxidant capacity (mmol FeSO4/g of extract).

Tyrosinase inhibitory activity

Tyrosinase inhibitory activities were determined by spectrophotometric method using L-DOPA as substrate [34]. Briefly, 50μL sample solutions at different concentrations were mixed with 50 μL tyrosinase (200 units/mL, dissolved in pH 6.8 phosphate buffer) in 96-well plates and incubated for 15 min at 25 °C. The reaction was then initiated with the addition of L-DOPA (50 μL). After incubating of 30 min at 25 °C, the absorbance at 490 nm (Aa) was determined using a multiskan sky microplate spectrophotometer (Thermo Fisher Scientific, America). Similarly, the absorbance of the sample wells without tyrosinase (Ab) and the control wells with enzyme but without sample (Ab) were detected at the same time. The tyrosinase inhibitory activity was calculated by the equation as below:

Tyrosinase inhibitory rate (%) = [1− (Aa − Ab)/Ac] × 100%.

The IC50 was calculated using the GraphPad Prism equation as below:

Y = min + (max − min)/1 + 10(x−logIC50)×Hill slope

X represented the inhibitor concentration; Y represented the inhibition data (%) along with their minimal (min) and maximal (max) values; Hill slope is the slope factor.

Cistanche inhibits tyrosinase activity.

Melanin inhibitory in zebrafish

Wild-type AB line zebrafish (provided by Hunter Bio-technology, Inc., Hangzhou, Zhejiang Province) were housed in fish water (0.2% instant ocean salt in deionized water, pH 6.9–7.2, conductivity 480–510 mS/cm, and hardness 53.7–71.6mg/L CaCO3) in a 14/10h light/ dark photoperiod at a constant temperature 28±0.5°C, and fed live brine shrimp twice daily. The zebrafish embryos were obtained from natural spawning and collected within 30min. The zebrafish assay was accredited by the Association for Assessment and Accreditation of Laboratory Animal Care (AAALAC). The present study was approved by the IACUC (Institutional Animal Care and Use Committee) at Hunter Biotechnology, Inc. and the IACUC approval number was 001458.

Embryos at the 6h post fertilization (hpf) stage were treated each extract at six concentrations (10, 30, 62.5, 125, 250 and 500 mg/L) for 72 h to evaluating the maximum non-lethal concentration (MNLC). As a result, the MNLCs of all the extract were more than 62.50mg/L. 10 embryos were placed in 6-well and were exposed to tested samples at concentrations of 10 and 30 mg/L from 6 hpf to 54 hpf (48h exposure). DMSO (0.05%, v/v) and arbutin (10 and 30mg/L) were used as the normal and positive control, respectively.

Synchronized embryos were collected and observed under a stereomicroscope (SZX7, Olympus, Japan) equipped with a digital camera (VertA1, China). Melanin accumulation, directly correlated to the zebrafish pigmentation areas, was calculated using the GNU Image Manipulation Program (GIMP version 2.10.2) and expressed as percentage with respect to the negative control (melanin accumulation 100%) [28, 30].

Molecular docking study

158 compounds isolated from B. striata in the literature [2], including 19 glycosides, 28 bibenzyls, 19 phenan- threnes, 18 biphenanthrenes, 23 dihydrophenanthrenes, 5 anthocyanins, 11 steroids, 8 triterpenoids, 12 phe- nolic acids, 5 quinones, and 10 other compounds, were used as the ligands. The 3D structures were drawn by ChemBio3D and were optimized using MM2 and Autodock Tools. The 3D crystal structure of tyrosinase (PDB ID:5M8N) and adenylate cyclase (PDB ID:5IV3) were retrieved from RCSB Protein Data Bank (www.rcsb. org/pdb/home/home.do). All water molecules and hetero atoms were removed from the crystal structures by using Chimera and MGL Tools. Finally, Autodock

Vina was used to dock the receptor protein with the small molecule ligands. Parameters of the receptor pro- tein docking site were set as −36.700, 7.034, −19.104, and − 17.467, − 22.807, 2.786, respectively, according to the original ligand mimosine and LRE1 [35]. In order to achieve higher computational accuracy, 20 exhaustiveness parameters for each receptor were generated, and the conformation with the highest affinity was selected as the final docking conformation and visualized in Pymol 2.3. Meanwhile, the original ligands mimosine and LRE1 were taken as the positive control [36].

In silico ADMET prediction

The top 3 hits from B. striata possessed better binding affinity with tyrosinase and adenylate cyclase, were subjected to ADMET analysis. QikProp in Schrodinger suite was used to calculate their physical properties and drug-related characteristics, including molecular weight (MW), predicted octanol/water partition coefficient (QPlogPo/w), predicted aqueous solubility (QPlogS), predicted apparent Caco-2 cell permeability for gut- blood barrier (QPPCaco), predicted brain/blood partition coefficient (QPlogBB), predicted skin permeability (QPlogKp), predicted IC50 for blockage of HERG K+ channels (QPlog HERG) and predicted human oral absorption. The properties were assessed based on Lipin- ski’s rule of five and literature [37, 38].

Molecular dynamics simulation

In order to examine the stability and dynamic fluctuations of the ligand-protein complex under a simulated biological environment, and further validate the docking results, the compounds blestrin D (75) was chosen as the ligand for the molecular dynamics (MD) simulation study, based on the molecular docking result. The complex of blestrin D with tyrosinase and adenylate cyclase were performed MD simulations and run for 100ns, using TIP3P water mode. The root mean square deviation (RMSD) values of protein backbone atoms relative to the initial structure were calculated to examine the protein stability over the course of simulation period [39].

Results

Chemical composition

The typical base peak chromatogram chromatograms of ETB and EFB were shown in Fig. 1. The PeakView software was used to identify the constituents. The molecular formula was accurately assigned within mass error of 10ppm. The exact molecular weight and the fragment ions were used to identify the components according the literatures and the free chemical structure database, such as ChemSpider and Massbank. As a result, 39 chemical compositions, including 24 stilbenoids (bibenzyls, phenanthrenes and their derivatives), 6 glycosides, 4 phenolic acids, 3 quinones, 1 steroid, and 1 other compound were tentatively identified from EFB and ETB. At the same time, their relative semi-quantification was performed by measuring peak areas of each compound in MS mode using the extracted ion chromatograms [40]. The detailed information of the identification and their relative contents (heat map highlights, the darker the color, the higher the concentration) were summarized in Table 1.

Antioxidant capacity

It is well known that Ultraviolet A (UVA)-irradiation can induce reactive oxygen species (ROS) production and mediate excessive melanogenesis in skin cells. Natural polyphenol was the inhibitors of ROS generation and could be responsible for the anti-melanogenic activity of plant extracts [41, 42]. Thus, in the present study, the antioxidant capacity was assessed through DPPH, ABTS radical scavenging activity and FRAP assay. The results (Table 2 and Fig. 2) showed that the EFB (IC50=5.94mg/L) possesses the strongest DPPH radical scavenging activity in vitro compared to the PTB (IC50 = 548.24 mg/L), PFB (IC50 = 285.81 mg/L), and ETB (IC50=65.25mg/L). A similar trend was observed in ABTS and FRAP assay. The 95% ethanolic extracts of TB and PB had stronger antioxidation activity in vitro than their crude polysaccharides. In addition, EFB exhibit stronger antioxidation activity than ETB, which was consistent with the previous research results [25].

Tyrosinase inhibitory activity

Tyrosinase was the key rate-limiting enzyme in melanin biosynthesis pathway and was widely used to identify the potential of natural products with anti-melanogenesis effect [43]. The results (Fig. 3) showed that ETB and EPB had stronger tyrosinase inhibitory activity than PTB and PPB in vitro, which was consistent with previous research results [25]. EFB exhibited stronger tyrosinase inhibition activity in a dose dependent manner (between 20 and 200mg/L) with IC50=58.92mg/L than ETB (IC50=75.44mg/L) and the positive compound arbutin (IC50 = 66.02 mg/L).

Melanin inhibitory in zebrafish

Zebrafish are recognized as a highly advantageous vertebrate model system to evaluate anti-melanogenesis activity, as it possesses similar organ systems and gene sequences to human beings [28]. Therefore, zebrafish was used to evaluate the anti-melanogenesis activity. As shown in Fig. 4, a large number of melanin (black spots) deposited in the zebrafish embryo in the control and 0.05% DMSO group. The microscopy images showed there was no significant difference in total melanin content between two groups (p>0.05), indicating that the vehicle 0.05% DMSO didn’t affect the melanin production in zebrafish embryos. However, after exposure to the tested samples and arbutin for 48h, zebrafish embryos exhibited varying degrees of melanin deposition (Figs. 3B,4), excepting PFB 10 mg/L group. It is noticeable that the relative melanin contents in ETB 10mg/L (78.38%) and 30 mg/L (64.72%) groups were significant lower than that in PTB 10mg/L (93.06%) and 30mg/L (82.02%) groups (p<0.01). It demonstrated that the anti-melanogenesis activity of ETB was superior to that of PTB, which was consistent with the tyrosinase inhibitory activity. Whereafter, it is noteworthy that the anti-melanogenesis activities of ETB (81.65, 75.61%) and EFB (78.38, 64.72%) were stronger than that of arbutin (93.57, 81.48%) at the same dose of 10 mg/L and 30 mg/L (p < 0.01 or 0.05).

Molecular docking

Tyrosinase is the rate-limiting enzyme in melanin biosynthesis: hydroxylation of tyrosine to 3,4-dihydroxy- phenylalanine (DOPA), and oxidation of DOPA to dopaquinone. Therefore, tyrosinase has been considered as a critical target for the development of melanogenesis inhibitors [44]. Additionally, adenylate cyclase is key enzyme of cAMP-induced melanogenesis. The binding of melanotropin alpha polypeptide (α-MSH) to MC1R receptors induced the activation of adenylate cyclase, increase in the cAMP level, up-regulated expression of tyrosinase gene, and subsequently increase the melanin synthesis [45]. Therefore, we performed a molecular docking study of the 158 compounds isolated from B. striata previously [2], with tyrosinase and adenylate cyclase. The binding energies of the studied ligands with tyrosinase and adenylate cyclase were summaried in Table S1. There were 83 (including 60 stilbenoids) and 85 (including 70 stilbenoids) compounds exhibit stronger binding affinities toward tyrosinase and adenylate cyclase, comparing to those of the original ligands mimosine (− 6.4 kcal/mol) and LRE1 (− 8.5 kcal/mol). 1,8-bis(p-hydroxybenzyl)- 4-methoxyphenanthrene-2,7-diol (56), blestrin D (75), and 2,7-dihydroxy-1,6-bis(p-hydroxybenzyl)-4-meth- oxy-9,10-dihydrophenanthrene (93) were the top three ligands toward tyrosinase, with binding energy −10.2, 10.0, and 9.7kcal/mol, respectivey. While blestrin D (75), blestrin B (73), and 3,3′-dihydroxy-5-methoxy- 2,5′,6-tris(p-hydroxybenzyl) bibenzyl (42) were the top three ligands toward adenylate cyclase, with binding energy − 12.1, − 11.9, and − 11.6 kcal/mol, respectivey. Their theoretical binding affinities and ligand-amino acid interactions were summaried in Table 3.

It was worth noting that blestrin D (75), a biphenanthrene compound, entered into the binding pocket of tyrosinase and adenylate cyclase (Fig. 5), and exhibited significant binding affinities toward both tyrosinase and adenylate cyclase. The Van der Walls interactions contributes to overall energy interaction value, meanwhile hydroxyl group produces hydrogen bonds with amino acid residues Glu 451, Arg114 of tyrosinase, and Val 167 of adenylate cyclase, respectivey (Table 3 and Fig. S1). Compound 93 also exhibited significant binding affinities toward tyrosinase through formation of 6 hydrogen bonds with amino acid residues Glu232, Glu451, Ser106, Cys113, Lys223, and Gln236.

In silico ADMET prediction

The predicted values of several important parameters along with their acceptable range were summarized in Table 4. Apart from the QPlogS of 3,3′,5-trimethoxybibenzyl, blestrin B and blestrin D, and QPlogBB of 3,3′,5-trimethoxybibenzyl, all the other calculated ADMET properties of the top three hits were within the expected ranges, meanwhile, the QPlog HERG of all the top three hits were less than −5, indicating they were drugable potential for human use. The QPlog Kp of all the top three hits for both tyrosinase and adenylate cyclase, excepting 3,3′,5-trimethoxybibenzyl, were in the favorable range, indicating these compounds can easily penetrate into skin.

Molecular dynamic simulation

The RMSD plot showed that the RMSD values of the protein-ligand complex did not fluctuate significantly throughout the entire simulation period. The RMSD values (Fig. 6A) of blestrin D with tyrosinase and ade- nylate cyclase complex were determined ranging from 1.5 to 2.8Å, and from 2.0 to 3.6Å, respectively. The RMSF value reflecting the flexibility of each residue during simulations. The residues with higher RMSF values revealed they were more flexible (Fig. 6B). The RMSF values of the amino acid residues in loop regions were found to highly fluctuate. While, the active site residues kept a small RMSF value (less than 1.0Å), such as Arg114, Val167, Tyr226, Leu229, and Arg230 in tyrosinase, and Val167, Leu102, Phe45, Lys95 and Leu166 in adenylate cyclase, indicating the binding pockets were kept stable. The molecular interactions of blestrin D with tyrosinase and adenylate cyclase were shown in Fig. 7.

The binding free energies were calculated by molecular mechanics-generalized born surface area (MM-GBSA) method [39]. The predicted binding free energy of blestrin D with tyrosinase and adenylate cyclase were−96.94kcal/mol and−139.69kcal/mol, respectively, implying the binding interactions were spontaneous (Table 5). Furthermore, the contributions favoring ligand binding were nonpolar solvation interaction, van der Waals and electrostatic interaction.

It was worth noting that the types and contents of the chemical components in EFB were more than that in ETB. The peak areas of all the identified compounds in EFB were approximately two-fold of that in ETB. Mili- tarine was the most abundant compounds in both EFB and ETB, which possessed significant antioxidant, anti-inflammatory and neuroprotection activities, and was chosen as the chemical markers for the quality control of B. striata in Chinese Pharmacopoeia 2020 edition [46]. The relative peak area of militarine in EFB (38.39×106) was slightly higher than that in ETB (30.48×106). Gastrodin is a well-known compound in many Chinese herbal medicines, such as Gastrodia elata, which exhibits significant tyrosinase inhibitory and radi- cal scavenging effects [30]. The relative peak area of gastrodin in EFB (1.43×106) was slightly lower than that in ETB (2.17 × 106).

The 95% ethanolic extracts of TB and PB had stronger antioxidation activity in vitro than their crude polysac- charides. In addition, EFB exhibit stronger antioxidation activity than ETB, which was consistent with the previous research results [25]. This phenomenon may be attributed to that the EFB contains higher level of antioxidant compounds, such as militarine and stilbenoids (natural plant polyphenols). For example, the relative peak area of 4,8,4′,8′-tetramethoxy-[1,1′-biphenanthrene]-2,7,2′,7′- tetrol, a biphenanthrene with four phenolic hydroxyl groups, was 9.15×106 in FB, while it was 0.01×106 in TB.

In the present study, PPB and PTB showed slight tyrosinase inhibition activity in a dose dependent manner. However, it demonstrated that no activity was detected in the aqueous extract of PB (contains PPB) [25]. This phenomenon may be due to the different sample preparation. The aqueous extract of PB was used to the evaluate tyrosinase inhibitory activity in Jiang’s experiment [25], while in the present study, aqueous extract was precipitation with ethanol to obtain the crude polysaccharide.

Recently, the zebrafish model widely used in evaluation of anti-melanogenesis effect in vivo [47, 48]. Arbu- tin, a hydroquinone glucoside compound existing in various plants, which has been commercially used in the cosmetic industry, was set as positive control. The anti-melanogenesis activities of ETB and EFB were stronger than that of arbutin, indicating that ETB and EFB can be applied as skin-lightening agents in cosmetic industry.

The relative peak area of blestrin D in EFB(1.91 × 106), is approximately four times of that in ETB (0.46 × 10 ). The relative peak area of compound 56 (8.73 × 106) and 93 (0.64 × 106) in EFB is also significant higher than that in ETB (less than 0.01×106 and 0.10 × 106). This may be one of the reasons why EFB exhibits stronger anti-melanogenic effects than ETB. It is interesting to note that top three ligands toward tyrosinase and adenylate cyclase all belong to stilbenoids (bibenzyls, phenanthrenes and its derivative). Stilbenoids are the natural plant polyphenols, which has attracted great interest in the last years because of their remarkable bioactivities such as anti-inflammatory, antimicrobial and antioxidant activity [49]. Resveratrol is the most common stilbenoid. The previous research indicated that resveratrol and its derivatives can significantly inhibit the catalytic activity, gene expression, and posttranslational modifications of tyrosinase. Thus, they showed potentially useful as skin lightening and antiaging agents in cosmetics [41, 50].

Cistanche showed potentially useful as skin lightening and antiaging agents in cosmetics.

For more information, please click the picture.

Conclusions

In this study, the antioxidant and anti-melanogenic activity of the crude polysaccharide from the B. striatatuber (PTB) and fibrous roots (FPB), and 95% ethanol extract of B. striata tuber (ETB) and fibrous roots (EFB) were systematically investigated. The results showed that EFB possessed the strongest DPPH, ABTS radical scavenging activity and ferric iron reducing activity. In addition, ETB and EFB can significantly reduce tyrosinase activity in vitro and melanin synthesis of zebrafish embryos in a dose-dependent manner. Molec- ular docking indicated that there were a large number of compounds, mostly belonged to stilbenoids, exhib- iting stronger binding affinities toward tyrosinase and adenylate cyclase, compared to the binding affinity of the original ligand. The present findings supported the rationale for the use of ETB and EFB as natural skin- whitening agents in pharmaceutical and cosmetic industries.