Design, Synthesis, In Vitro, and In Silico Insights of 5-(Substituted benzylidene)-2-phenylthiazol-4(5H)-one Derivatives: A Novel Class of Anti-Melanogenic Compounds

(Z)-5-Benzylidene-2-phenylthiazol-4(5H)-one ((Z)-BPT) derivatives were designed by combining the structural characteristics of two tyrosinase inhibitors. The double-bond geometry of trisubstituted alkenes, (Z)-BPTs 1–14, was determined based on the 3JC,Hβ coupling constant of 1H-coupled 13C NMR spectra. Three (Z)-BPT derivatives (1–3) showed stronger tyrosinase inhibitory activities than kojic acid; in particular, 2 was to be 189-fold more potent than kojic acid. Kinetic analysis using mushroom tyrosinase indicated that 1 and 2 were competitive inhibitors, whereas 3 was a mixed-type inhibitor. The in silico results revealed that 1–3 could strongly bind to the active sites of mushroom and human tyrosinases, supporting the kinetic results. Derivatives 1 and 2 decreased the intracellular melanin contents in a concentration-dependent manner in B16F10 cells, and their anti-melanogenic efficacy exceeded that of kojic acid. The anti-tyrosinase activity of 1 and 2 in B16F10 cells was similar to their anti-melanogenic effects, suggesting that their anti-melanogenic effects were primarily owing to their anti-tyrosinase activity. Western blotting of B16F10 cells revealed that the derivatives 1 and 2 inhibited tyrosinase expression, which partially contributes to their anti-melanogenic ability. Several derivatives, including 2 and 3, exhibited potent antioxidant activities against ABTS cation radicals, DPPH radicals, ROS, and peroxynitrite. These results suggest that (Z)-BPT derivatives 1 and 2 have promising potential as novel anti-melanogenic agents.


Introduction
Human skin has three layers: the epidermis, dermis, and subcutaneous layer [1]. The epidermis is the top layer of the skin and acts as the first line of defense against external factors, such as pathogens and UV irradiation. The epidermis primarily consists of keratinocytes and melanocytes in the basal layer [1]. Melanocytes contain special organelles called melanosomes, which produce melanin. Melanin produced in the melanosomes of melanocytes is transferred to keratinocytes [2,3]. In mammals, there are

Synthesis of (Z)-BPT Derivatives 1-14
As described in Scheme 1, the key intermediate 16 [43] for synthesizing (Z)-BPT derivatives 1-14 was prepared via a two-step reaction: i) cyclization and ii) HBr salt removal. The reflux of thiobenzamide and bromoacetic acid in ethyl acetate produced 2-phenylthiazol-4(5H)-one hydrobromide 15 with a 69% yield, and treatment with pyridine produced 2-phenylthiazol-4(5H)-one 16 with a 92% yield. The condensation of 16 with diverse substituted benzaldehydes (a-n) in the presence of piperidine (0.3 equiv.) in refluxed ethyl alcohol generated the desired (Z)-BPT derivatives 1-14 as single products with yields of 41-81%. The structures of the final compounds were confirmed using NMR ( 1 H and 13 C) spectroscopy and mass spectrometry (Supplementary Materials). The double-bond stereochemistry of the trisubstituted alkenes 1-14, generated from the condensation between benzaldehydes and 16, was determined based on the 3 J (between vicinal 1 H and 13 C) coupling constants of 1 H-coupled 13 C NMR spectra. Nair et al. previously reported that the double-bond geometries in trisubstituted enamides could be differentiated based on the observed vicinal 1 H and 13 C-coupling constant ( 3 JC,Hβ) of the amide carbonyl carbon [44]. According to the report, vicinal coupling constants ranging from 3.6 to 7.0 Hz indicate that Hβ and the carbonyl carbon of the amide are on the same side. On the other hand, the vicinal coupling constants (≥7.5 Hz), which are usually greater than 10 Hz, indicate that Hβ and the carbonyl carbon are on the opposite sides. To confirm the double-bond geometry of the (Z)-BPT derivatives, the 13 C NMR of 2 was measured in the 1 H-coupled NMR mode. The amide carbonyl carbon had a vicinal 1 H and 13 C-coupling constant of 5.3 Hz

Results and Discussion
2.1. Synthesis of (Z)-BPT Derivatives 1-14 As described in Scheme 1, the key intermediate 16 [43] for synthesizing (Z)-BPT derivatives 1-14 was prepared via a two-step reaction: i) cyclization and ii) HBr salt removal. The reflux of thiobenzamide and bromoacetic acid in ethyl acetate produced 2-phenylthiazol-4(5H)-one hydrobromide 15 with a 69% yield, and treatment with pyridine produced 2-phenylthiazol-4(5H)-one 16 with a 92% yield. The condensation of 16 with diverse substituted benzaldehydes (a-n) in the presence of piperidine (0.3 equiv.) in refluxed ethyl alcohol generated the desired (Z)-BPT derivatives 1-14 as single products with yields of 41-81%. The structures of the final compounds were confirmed using NMR ( 1 H and 13 C) spectroscopy and mass spectrometry (Supplementary Materials). The doublebond stereochemistry of the trisubstituted alkenes 1-14, generated from the condensation between benzaldehydes and 16, was determined based on the 3 J (between vicinal 1 H and 13 C) coupling constants of 1 H-coupled 13 C NMR spectra. Nair et al. previously reported that the double-bond geometries in trisubstituted enamides could be differentiated based on the observed vicinal 1 H and 13 C-coupling constant ( 3 J C,Hβ ) of the amide carbonyl carbon [44]. According to the report, vicinal coupling constants ranging from 3.6 to 7.0 Hz indicate that H β and the carbonyl carbon of the amide are on the same side. On the other hand, the vicinal coupling constants (≥7.5 Hz), which are usually greater than 10 Hz, indicate that H β and the carbonyl carbon are on the opposite sides. To confirm the double-bond geometry of the (Z)-BPT derivatives, the 13 C NMR of 2 was measured in the 1 H-coupled NMR mode. The amide carbonyl carbon had a vicinal 1 H and 13 C-coupling constant of 5.3 Hz (see S6. The proton-coupled 13 C NMR spectrum of compound 2 in the Supplementary Materials), suggesting that 2 has a (Z)-geometry.
x FOR PEER REVIEW 4 of 25 (see S6. The proton-coupled 13 C NMR spectrum of compound 2 in the Supplementary Materials), suggesting that 2 has a (Z)-geometry.

Tyrosinase Inhibition of (Z)-BPT Derivatives
Although melanogenesis is modulated by hormones, tyrosinase inhibition is the most effective method for controlling melanogenesis. Therefore, the inhibitory effect of 14 synthetic (Z)-BPT derivatives on mushroom tyrosinase was investigated according to a previously described procedure [39]. For evaluating tyrosinase inhibitors, L-tyrosine and kojic acid were utilized as a substrate and a positive control, respectively. All the 14 (Z)-BPT derivatives inhibited the mushroom tyrosinase activity in a concentration-dependent manner. Table 1 shows the IC50 values of (Z)-BPT derivatives and kojic acid.

Tyrosinase Inhibition of (Z)-BPT Derivatives
Although melanogenesis is modulated by hormones, tyrosinase inhibition is the most effective method for controlling melanogenesis. Therefore, the inhibitory effect of 14 synthetic (Z)-BPT derivatives on mushroom tyrosinase was investigated according to a previously described procedure [39]. For evaluating tyrosinase inhibitors, L-tyrosine and kojic acid were utilized as a substrate and a positive control, respectively. All the 14 (Z)- BPT derivatives inhibited the mushroom tyrosinase activity in a concentration-dependent manner. Table 1 shows the IC 50 values of (Z)-BPT derivatives and kojic acid. Tyrosinase inhibition experiments were performed independently in triplicates. Ltyrosine and kojic acid (KA) were used as the substrate and positive controls, respectively. The IC50 values are presented as the mean ± standard error of the mean (SEM). ChemDraw Ultra 12.0 was used to obtain the log P values.

Identifying the Inhibition Type of (Z)-BPT Derivatives 1-3
Since (Z)-BPT derivatives 1-3 exhibited more potent tyrosinase inhibitory activity than kojic acid, further studies were performed to determine their inhibition type. To investigate the type of inhibition of mushroom tyrosinase, kinetic studies of 1-3 were conducted using mushroom tyrosinase with various concentrations of L-DOPA as a substrate. As depicted in Figure 2, the mechanisms of tyrosinase inhibition were analyzed using Lineweaver-Burk plots. (Z)-BPT derivatives 1 and 2 showed similar Lineweaver-Burk plot patterns. Each of the four straight lines produced at each inhibitor concentration converged to a single point on the y-axis. Each Vmax value was constant, regardless of the inhibitor (1 and 2) concentration, and each KM value gradually increased as the concentration of each inhibitor increased. These results imply that the (Z)-BPT derivatives 1 and 2 bind to the same binding pocket as L-DOPA, a tyrosinase substrate, and inhibit mushroom tyrosinase activity in a competitive manner (Table 2). In contrast, in derivative 3, four straight lines generated at each inhibitor concentration intersected at a point in the second quadrant. According to the Lineweaver-Burk plot, as the inhibitor concentration increased, the Vmax value of 3 decreased, but the KM value increased. These results suggest that (Z)-BPT derivative 3 is a mixed-type inhibitor.
The (Z)-BPT derivatives 8-11, with no hydroxyl group, showed no potent tyrosinase inhibition (IC 50 > 300 µM). Although derivative 7 contained a 2-hydroxyl substituent on the β-phenyl ring of the PUSC template, it did not exhibit mushroom tyrosinase inhibitory activity. Derivative 1 contained a 4-hydroxyl substituent on the β-phenyl ring and exhibited a strong tyrosinase inhibitory activity (IC 50 = 6.4 ± 0.52 µM), and inhibited mushroom tyrosinase more strongly than kojic acid (IC 50 = 20.8 ± 1.34 µM). The introduction of a 3-bromo substituent on the β-phenyl ring of derivative 1 dramatically reduced the tyrosinase inhibitory activity (IC 50 value of derivative 13: 253.7 ± 1.87 µM). However, when an additional 5-bromo substituent was inserted into the β-phenyl ring of derivative 13, the tyrosinase inhibitory activity increased to some extent (IC 50 value of derivative 14: 87.2 ± 0.63 µM). The introduction of a 3-alkoxyl substituent on the β-phenyl ring of derivative 1 also decreased the tyrosinase inhibitory activity (IC 50 values of derivative 4 with a 3-methoxyl and derivative 5 with a 3-ethoxyl: 211.1 ± 1.84, and 73.8 ± 2.02 µM, respectively). The exchange of 3-methoxyl and 4-hydroxyl substituents in derivative 4 slightly increased the tyrosinase inhibitory activity (IC 50 value of derivative 6: 121.1 ± 2.72 µM). The introduction of an additional 3-hydroxyl substituent on the β-phenyl ring of derivative 1 did not influence the tyrosinase inhibitory activity, whereas the insertion of an additional 2-hydroxyl substituent greatly enhanced the inhibition of tyrosinase activity. Derivative 2 contained a 2,4-dihydroxyl substituent and had an IC 50 value of 0.1 ± 0.01 µM, indicating that 2 was a 190-fold stronger inhibitor than kojic acid. These findings imply that the 4-hydroxyl substituent on the β-phenyl ring of the PUSC template in (Z)-BPT derivatives plays a significant role in inhibiting tyrosinase activity, and that the type and position of additional substituents are closely correlated with the level of tyrosinase activity inhibition. As in previous studies [17,45], derivative 2 with 2,4-dihydroxyphenyl exhibited the strongest tyrosinase inhibition. Log P represents the logarithm of the partition coefficient (P). The greater the log P value of a compound, the better its absorption by the skin. The log P values of (Z)-BPT derivatives and kojic acid were obtained from ChemDraw Ultra 12.0. As shown in Table 1, kojic acid had a log P value of −2.45, whereas (Z)-BPT derivatives had log P values ranging from 3.16 to 4.38.
Tyrosinase inhibition experiments were performed independently in triplicates. Ltyrosine and kojic acid (KA) were used as the substrate and positive controls, respectively. The IC 50 values are presented as the mean ± standard error of the mean (SEM). ChemDraw Ultra 12.0 was used to obtain the log P values.

Identifying the Inhibition Type of (Z)-BPT Derivatives 1-3
Since (Z)-BPT derivatives 1-3 exhibited more potent tyrosinase inhibitory activity than kojic acid, further studies were performed to determine their inhibition type. To investigate the type of inhibition of mushroom tyrosinase, kinetic studies of 1-3 were conducted using mushroom tyrosinase with various concentrations of L-DOPA as a substrate. As depicted in Figure 2, the mechanisms of tyrosinase inhibition were analyzed using Lineweaver-Burk plots. (Z)-BPT derivatives 1 and 2 showed similar Lineweaver-Burk plot patterns. Each of the four straight lines produced at each inhibitor concentration converged to a single point on the y-axis. Each V max value was constant, regardless of the inhibitor (1 and 2) concentration, and each K M value gradually increased as the concentration of each inhibitor increased. These results imply that the (Z)-BPT derivatives 1 and 2 bind to the same binding pocket as L-DOPA, a tyrosinase substrate, and inhibit mushroom tyrosinase activity in a competitive manner (Table 2). In contrast, in derivative 3, four straight lines generated at each inhibitor concentration intersected at a point in the second quadrant. According to the Lineweaver-Burk plot, as the inhibitor concentration increased, the V max value of 3 decreased, but the K M value increased. These results suggest that (Z)-BPT derivative 3 is a mixed-type inhibitor.
Tyrosinase inhibition experiments were performed independently in triplicates. Ltyrosine and kojic acid (KA) were used as the substrate and positive controls, respectively. The IC50 values are presented as the mean ± standard error of the mean (SEM). ChemDraw Ultra 12.0 was used to obtain the log P values.

Identifying the Inhibition Type of (Z)-BPT Derivatives 1-3
Since (Z)-BPT derivatives 1-3 exhibited more potent tyrosinase inhibitory activity than kojic acid, further studies were performed to determine their inhibition type. To investigate the type of inhibition of mushroom tyrosinase, kinetic studies of 1-3 were conducted using mushroom tyrosinase with various concentrations of L-DOPA as a substrate. As depicted in Figure 2, the mechanisms of tyrosinase inhibition were analyzed using Lineweaver-Burk plots. (Z)-BPT derivatives 1 and 2 showed similar Lineweaver-Burk plot patterns. Each of the four straight lines produced at each inhibitor concentration converged to a single point on the y-axis. Each Vmax value was constant, regardless of the inhibitor (1 and 2) concentration, and each KM value gradually increased as the concentration of each inhibitor increased. These results imply that the (Z)-BPT derivatives 1 and 2 bind to the same binding pocket as L-DOPA, a tyrosinase substrate, and inhibit mushroom tyrosinase activity in a competitive manner ( Table 2). In contrast, in derivative 3, four straight lines generated at each inhibitor concentration intersected at a point in the second quadrant. According to the Lineweaver-Burk plot, as the inhibitor concentration increased, the Vmax value of 3 decreased, but the KM value increased. These results suggest that (Z)-BPT derivative 3 is a mixed-type inhibitor. Lineweaver-Burk plots for the inhibition of mushroom tyrosinase activity by (Z)-BPT derivatives 1-3. The inhibition modes of (Z)-BPT derivatives 1-3 were demonstrated by the analysis of Lineweaver-Burk plots. Three independent experiments with five different concentrations (0.5, 1, 2, 4, and 8 mM) of an L-DOPA substrate were used to measure the augmented dopachrome optical density per minute at 492 nm. The data are shown as mean values of 1/V (an inverse of that optical density per minute). The (Z)-BPT derivatives were used at concentrations of 0 (filled circles), 3 (unfilled circles), 6 (filled triangles), and 12 µM (unfilled triangles) for 1; 0 (filled circles), 0.05 (unfilled circles), 0.1 (filled triangles), and 0.2 µM (unfilled triangles) for 2; and 0 (filled circles), 2.5 (unfilled circles), 5 (filled triangles), and 10 µM (unfilled triangles) for 3. Each experiment was run independently in triplicates. Lineweaver-Burk plots for the inhibition of mushroom tyrosinase activity by (Z)-BPT derivatives 1-3. The inhibition modes of (Z)-BPT derivatives 1-3 were demonstrated by the analysis of Lineweaver-Burk plots. Three independent experiments with five different concentrations (0.5, 1, 2, 4, and 8 mM) of an L-DOPA substrate were used to measure the augmented dopachrome optical density per minute at 492 nm. The data are shown as mean values of 1/V (an inverse of that optical density per minute). The (Z)-BPT derivatives were used at concentrations of 0 (filled circles), 3 (unfilled circles), 6 (filled triangles), and 12 µM (unfilled triangles) for 1; 0 (filled circles), 0.05 (unfilled circles), 0.1 (filled triangles), and 0.2 µM (unfilled triangles) for 2; and 0 (filled circles), 2.5 (unfilled circles), 5 (filled triangles), and 10 µM (unfilled triangles) for 3. Each experiment was run independently in triplicates. Table 2. Inhibition type of (Z)-BPT derivatives 1-3 and kojic acid.

Compound
Inhibition Type a Compound Inhibition Type The inhibition type was analyzed using a Lineweaver-Burk plot, and L-DOPA was used as a substrate, b as previously reported in the presence of L-DOPA.
To determine the K i (inhibition constant) values of (Z)-BPT derivatives 1-3, Lineweaver-Burk plots for 1-3 were transformed into the corresponding Dixon plots (Figure 3). Each Dixon plot obtained from each Lineweaver-Burk plot produced four straight lines that intersected at one point in the second quadrant. The vertical lines from the merge point to the x-axis provide the K i value of each (Z)-BPT derivative. The K i values of (Z)-BPT derivatives 1-3 were 8.89, 0.10, and 4.80 µM, respectively, indicating that 2 generates the strongest enzyme-inhibitor complexes.  The inhibition type was analyzed using a Lineweaver-Burk plot, and L-DOPA was used as a substrate, b as previously reported in the presence of L-DOPA.
To determine the Ki (inhibition constant) values of (Z)-BPT derivatives 1-3, Lineweaver-Burk plots for 1-3 were transformed into the corresponding Dixon plots ( Figure  3). Each Dixon plot obtained from each Lineweaver-Burk plot produced four straight lines that intersected at one point in the second quadrant. The vertical lines from the merge point to the x-axis provide the Ki value of each (Z)-BPT derivative. The Ki values of (Z)-BPT derivatives 1-3 were 8.89, 0.10, and 4.80 µM, respectively, indicating that 2 generates the strongest enzyme-inhibitor complexes.

In Silico Docking of (Z)-BPT Derivatives 1-3 with Mushroom Tyrosinase
Using (Z)-BPT derivatives 1-3, which exhibit effective tyrosinase inhibition against mushroom tyrosinase, an in silico docking simulation was performed to examine their behavior in the mushroom tyrosinase active site. Schrodinger Suite (release 2021-2) and kojic acid were used as the docking software and a positive reference, respectively. The results of the in silico docking simulations are shown in Figure 4.
As in kojic acid, all the three (Z)-BPT derivatives 1-3 bound to the active site of tyrosinase and showed similar or greater binding affinities than kojic acid with a docking score of −4.5 kcal/mol (1: −4.3, 2: −6.2, and 3: −6.3 kcal/mol) ( Figure 4B). Interestingly, all compounds (kojic acid and 1-3) interacted with one or two copper ions in the active site. Kojic acid was found to interact with tyrosinase through three modes: metal coordination between Cu401 and a 2-hydroxymethyl, a hydrogen bond between a 5-hydroxyl and Met280, and π-π stacking between Ser282 and a pyranone ring. In (Z)-BPT derivative 1, a metal coordination between a phenolic 4-hydroxyl and Cu401 was formed, and three ππ stacking interactions (between a β-phenyl ring and His259 and His263, and a 2-phenyl ring and His244) were created. Unlike (Z)-BPT derivative 1, the phenolic 4-hydroxyl of The K i values of (Z)-BPT derivatives 1-3 were obtained from the absolute values of the x-coordinate where four straight lines converged at one point. Three independent experiments with four different concentrations of (Z)-BPT derivatives were conducted to measure the augmented dopachrome optical density per minute at 492 nm. The data are shown as the mean values of 1/V (an inverse of that optical density per minute). The (Z)-BPT derivatives were used at concentrations of 0, 3, 6, and 12 µM for 1; 0, 0.05, 0.1, and 0.2 µM for 2; and 0, 2.5, 5, and 10 µM for 3. The experiments were performed at four different concentrations (0.5, 1, 2, and 4 mM) of L-DOPA. Each experiment was run independently in triplicates.

In Silico Docking of (Z)-BPT Derivatives 1-3 with Mushroom Tyrosinase
Using (Z)-BPT derivatives 1-3, which exhibit effective tyrosinase inhibition against mushroom tyrosinase, an in silico docking simulation was performed to examine their behavior in the mushroom tyrosinase active site. Schrodinger Suite (release 2021-2) and kojic acid were used as the docking software and a positive reference, respectively. The results of the in silico docking simulations are shown in Figure 4.
As in kojic acid, all the three (Z)-BPT derivatives 1-3 bound to the active site of tyrosinase and showed similar or greater binding affinities than kojic acid with a docking score of −4.5 kcal/mol (1: −4.3, 2: −6.2, and 3: −6.3 kcal/mol) ( Figure 4B). Interestingly, all compounds (kojic acid and 1-3) interacted with one or two copper ions in the active site. Kojic acid was found to interact with tyrosinase through three modes: metal coordination between Cu401 and a 2-hydroxymethyl, a hydrogen bond between a 5-hydroxyl and Met280, and π-π stacking between Ser282 and a pyranone ring. In (Z)-BPT derivative 1, a metal coordination between a phenolic 4-hydroxyl and Cu401 was formed, and three π-π stacking interactions (between a βphenyl ring and His259 and His263, and a 2-phenyl ring and His244) were created. Unlike (Z)-BPT derivative 1, the phenolic 4-hydroxyl of (Z)-BPT derivative 2 interacted with two copper ions (Cu400 and Cu401), not a single copper ion, through a salt bridge rather than a metal coordination. In addition, as seen in 1, 2 also formed three π-π stacking interactions using the same interaction patterns as 1. Furthermore, the 2-hydroxyl group on the β-phenyl ring of 2 participated in hydrogen bonding interactions with Glu256. The phenolic 4-hydroxyl of derivative 3 generated two salt bridges with Cu400 and Cu401, as in 2, and the 2-phenyl ring formed one π-cation interaction with Arg268. These results suggest that (Z)-BPT derivatives 1-3 can bind to the active site of tyrosinase via various interactions, including interactions with copper ions, such as kojic acid. (Z)-BPT derivative 2 interacted with two copper ions (Cu400 and Cu401), not a single copper ion, through a salt bridge rather than a metal coordination. In addition, as seen in 1, 2 also formed three π-π stacking interactions using the same interaction patterns as 1. Furthermore, the 2-hydroxyl group on the β-phenyl ring of 2 participated in hydrogen bonding interactions with Glu256. The phenolic 4-hydroxyl of derivative 3 generated two salt bridges with Cu400 and Cu401, as in 2, and the 2-phenyl ring formed one π-cation interaction with Arg268. These results suggest that (Z)-BPT derivatives 1-3 can bind to the active site of tyrosinase via various interactions, including interactions with copper ions, such as kojic acid.

In Silico Docking of (Z)-BPT Derivatives 1-3 with the Human Tyrosinase Model
Since the X-ray crystal structure of human tyrosinase is not available, we previously designed a human tyrosinase (hTYR) homology model based on tyrosinase-related protein-1. To predict the binding ability to human tyrosinase, (Z)-BPT derivatives 1-3 were docked with a previously designed hTYR homology model [34], and kojic acid was used as a reference material. The docking simulation results are shown in Figure 5.

In Silico Docking of (Z)-BPT Derivatives 1-3 with the Human Tyrosinase Model
Since the X-ray crystal structure of human tyrosinase is not available, we previously designed a human tyrosinase (hTYR) homology model based on tyrosinase-related protein-1. To predict the binding ability to human tyrosinase, (Z)-BPT derivatives 1-3 were docked with a previously designed hTYR homology model [34], and kojic acid was used as a reference material. The docking simulation results are shown in Figure 5. a docking score of −7.2 kcal/mol. In contrast, although derivative 3 had two hydroxyl groups on the β-phenyl ring like 2, the hydroxyls of 3 did not interact with zinc ions (Figure 5E). Instead, the two hydroxyl groups formed two hydrogen bonds with Ser184 and Arg196. In addition, the 2-phenyl ring participated in π-π stacking with Phe347. These interactions provided a slightly lower docking score (−3.9 kcal/mol) than that of kojic acid. These results suggest that the (Z)-BPT derivative 2 has potent anti-hTYR activity.  Kojic acid interacted with Zn7, Met374, and Ser375 through metal coordination, π-π stacking, and hydrogen bonding, respectively ( Figure 5A). These interactions provided a docking score of −4.4 kcal/mole ( Figure 5B). Derivative 1 participated in two interactions ( Figure 5C). The 4-hydroxyl of the β-phenyl ring in 1 hydrogen-bonded with Arg196, and the 2-phenyl ring interacted with Phe347 via π-π stacking. The docking score was −4.1 kcal/mol. Derivative 2 showed similar interactions with mushroom tyrosinase ( Figure 5D). The 4-hydroxyl of the β-phenyl ring in 2 participated in the formation of salt bridges with two zinc ions (Zn6 and Zn7), the 2-hydroxyl of the β-phenyl ring generated a hydrogen bond with Met374, and the β-phenyl ring produced a π-π stacking interaction with His367. These interactions with the homology model afforded a high binding affinity with a docking score of −7.2 kcal/mol. In contrast, although derivative 3 had two hydroxyl groups on the β-phenyl ring like 2, the hydroxyls of 3 did not interact with zinc ions ( Figure 5E). Instead, the two hydroxyl groups formed two hydrogen bonds with Ser184 and Arg196. In addition, the 2-phenyl ring participated in π-π stacking with Phe347. These interactions provided a slightly lower docking score (−3.9 kcal/mol) than that of kojic acid. These results suggest that the (Z)-BPT derivative 2 has potent anti-hTYR activity.

Effects of (Z)-BPT Derivatives 1-3 on B16F10 Cell Viability
To examine whether (Z)-BPT derivatives 1-3 had the ability to exhibit inhibitory effects on melanogenesis and cellular tyrosinase activity, B16F10 murine melanoma cells were utilized for cell experiments. Before investigating their inhibitory effects on B16F10 cells, the viability of (Z)-BPT derivatives 1-3 was examined in B16F10 cells. Cell viability was monitored for 72 h using the EZ-Cytox assay, and the results are shown in Figure 6.
To examine whether (Z)-BPT derivatives 1-3 had the ability to exhibit inhibitory effects on melanogenesis and cellular tyrosinase activity, B16F10 murine melanoma cells were utilized for cell experiments. Before investigating their inhibitory effects on B16F10 cells, the viability of (Z)-BPT derivatives 1-3 was examined in B16F10 cells. Cell viability was monitored for 72 h using the EZ-Cytox assay, and the results are shown in Figure 6.
To determine the cell viability, six different concentrations (0, 1, 2, 5, 10, and 20 µM) of (Z)-BPT derivatives 1-3 were used, and the results were examined at two time points (48 and 72 h). Derivatives 1 and 2 did not show any perceptible cytotoxicity at concentrations ≤ 20 µM at either time point, but derivative 3 showed cytotoxic effects at concentrations ≥ 5 µM at 72 h and at concentrations ≥ 2 µM at 72 h. Based on the cell viability results, the concentrations for cell-based assays related to the melanin production level and cellular tyrosinase activity were determined.

Effects of (Z)-BPT Derivatives 1 and 2 on Melanin Production in B16F10 Cells
Derivative 3 showed an excellent inhibitory activity against mushroom tyrosinase but exhibited cytotoxicity in B16F10 cells. Thus, the experiments evaluating the effect of 3 on cellular tyrosinase activity and melanin production in B16F10 cells were excluded. To confirm whether (Z)-BPT derivatives 1 and 2, which show potent mushroom tyrosinase inhibitory activity, can exert anti-melanogenic effects in mammalian cells, the melanin levels were measured in B16F10 cells. Cells were exposed to four different concentrations (0, 5, 10, and 20 µM) of 1 and 2 prior to treatment with stimulators (3-isobutyl-1methylxanthine (IBMX, 200 µM), and α-melanocyte stimulating hormone (α-MSH, 1 µM)). After 72 h of incubation, the effects of 1 and 2 on melanin production were determined. Cellular melanin levels were assessed by measuring melanin levels in cell lysates. The results are shown in Figure 7. To determine the cell viability, six different concentrations (0, 1, 2, 5, 10, and 20 µM) of (Z)-BPT derivatives 1-3 were used, and the results were examined at two time points (48 and 72 h). Derivatives 1 and 2 did not show any perceptible cytotoxicity at concentrations ≤ 20 µM at either time point, but derivative 3 showed cytotoxic effects at concentrations ≥ 5 µM at 72 h and at concentrations ≥ 2 µM at 72 h. Based on the cell viability results, the concentrations for cell-based assays related to the melanin production level and cellular tyrosinase activity were determined.

Effects of (Z)-BPT Derivatives 1 and 2 on Melanin Production in B16F10 Cells
Derivative 3 showed an excellent inhibitory activity against mushroom tyrosinase but exhibited cytotoxicity in B16F10 cells. Thus, the experiments evaluating the effect of 3 on cellular tyrosinase activity and melanin production in B16F10 cells were excluded. To confirm whether (Z)-BPT derivatives 1 and 2, which show potent mushroom tyrosinase inhibitory activity, can exert anti-melanogenic effects in mammalian cells, the melanin levels were measured in B16F10 cells. Cells were exposed to four different concentrations (0, 5, 10, and 20 µM) of 1 and 2 prior to treatment with stimulators (3-isobutyl-1-methylxanthine (IBMX, 200 µM), and α-melanocyte stimulating hormone (α-MSH, 1 µM)). After 72 h of incubation, the effects of 1 and 2 on melanin production were determined. Cellular melanin levels were assessed by measuring melanin levels in cell lysates. The results are shown in Figure 7. content that was increased by the stimulators. Derivative 1 at 20 µM reduced the cellular melanin content more effectively than kojic acid at the same concentration. However, derivative 2 exhibited more remarkable results. The inhibitory effect of 2 at 5 µM on the cellular melanin content was similar to that of 20 µM kojic acid, and derivative 2 at 20 µM reduced the cellular melanin levels that were enhanced by stimulators to levels below those of the non-treated control.

Effects of (Z)-BPT Derivatives 1 and 2 on B16F10 Cellular Tyrosinase Activity
To examine the origin of the anti-melanogenic effects of (Z)-BPT derivatives 1 and 2, the anti-tyrosinase effects of these derivatives were investigated using B16F10 cells. Similar to the melanin production assay, four different (Z)-BPT derivative concentrations (0, 5, 10, and 20 µM of 1 and 2) were exposed to B16F10 cells prior to treatment with stimulators, IBMX (200 µM) and α-MSH (1 µM). After 72 h of incubation, the effects of 1 and 2 on cellular tyrosinase activity were analyzed. The results are depicted in Figure 8; kojic acid was used as the positive control.
Treatment of stimulators (IBMX and α-MSH) with B16F10 cells highly increased the cellular tyrosinase activity by 3.3-fold compared to that of the untreated control (100%). Kojic acid (20 µM) treatment diminished the enhanced cellular tyrosinase activity by 2.4fold. Treatment with derivatives 1 and 2 dose-dependently and significantly decreased the cellular tyrosinase activity, which was enhanced by treatment with stimulators. Derivative 1 exerted a stronger inhibitory activity than kojic acid at 20 µM. Derivative 2 had a more potent cellular anti-tyrosinase activity than derivative 1 and kojic acid. At a concentration of 5 µM, derivative 2 had cellular tyrosinase inhibitory activity comparable to that of 20 µM kojic acid, and at a concentration of 20 µM, it reduced the cellular tyrosinase activity to the levels observed in the untreated control. The results of derivatives 1 and 2 for the cellular tyrosinase activity were similar to those for cellular melanin content, implying that the anti-melanogenic capacity of 1 and 2 in B16F10 cells was mainly attributable to their ability to inhibit cellular tyrosinase. B16F10 cells treated with stimulators (IBMX and α-MSH) greatly increased the cellular melanin content by 3.3-fold compared to the untreated control (100%). However, exposure to 1 and 2 significantly and dose-dependently diminished the cellular melanin content that was increased by the stimulators. Derivative 1 at 20 µM reduced the cellular melanin content more effectively than kojic acid at the same concentration. However, derivative 2 exhibited more remarkable results. The inhibitory effect of 2 at 5 µM on the cellular melanin content was similar to that of 20 µM kojic acid, and derivative 2 at 20 µM reduced the cellular melanin levels that were enhanced by stimulators to levels below those of the non-treated control.

Effects of (Z)-BPT Derivatives 1 and 2 on B16F10 Cellular Tyrosinase Activity
To examine the origin of the anti-melanogenic effects of (Z)-BPT derivatives 1 and 2, the anti-tyrosinase effects of these derivatives were investigated using B16F10 cells. Similar to the melanin production assay, four different (Z)-BPT derivative concentrations (0, 5, 10, and 20 µM of 1 and 2) were exposed to B16F10 cells prior to treatment with stimulators, IBMX (200 µM) and α-MSH (1 µM). After 72 h of incubation, the effects of 1 and 2 on cellular tyrosinase activity were analyzed. The results are depicted in Figure 8; kojic acid was used as the positive control.
Treatment of stimulators (IBMX and α-MSH) with B16F10 cells highly increased the cellular tyrosinase activity by 3.3-fold compared to that of the untreated control (100%). Kojic acid (20 µM) treatment diminished the enhanced cellular tyrosinase activity by 2.4-fold. Treatment with derivatives 1 and 2 dose-dependently and significantly decreased the cellular tyrosinase activity, which was enhanced by treatment with stimulators. Derivative 1 exerted a stronger inhibitory activity than kojic acid at 20 µM. Derivative 2 had a more potent cellular anti-tyrosinase activity than derivative 1 and kojic acid. At a concentration of 5 µM, derivative 2 had cellular tyrosinase inhibitory activity comparable to that of 20 µM kojic acid, and at a concentration of 20 µM, it reduced the cellular tyrosinase activity to the levels observed in the untreated control. The results of derivatives 1 and 2 for the cellular tyrosinase activity were similar to those for cellular melanin content, implying that the anti-melanogenic capacity of 1 and 2 in B16F10 cells was mainly attributable to their ability to inhibit cellular tyrosinase. Data are presented as the mean ± standard error of the mean. ### p < 0.001; * p < 0.05, ** p < 0.01, and *** p < 0.001.

ABTS Cation Radical Scavenging Effects of (Z)-BPT Derivatives 1-14
It has been suggested that substances with antioxidant qualities may prevent the tyrosinase substrate (L-tyrosine and L-dopa) oxidation required for melanogenesis without interacting with the enzyme [46,47]. In addition, it has been reported that antioxidants that scavenge radicals may inhibit excess melanogenesis, because melanogenesis can be triggered by various radicals [48][49][50][51][52]. Therefore, the 2,2′-azino-bis(3-ethylbenzothiazoline-6-sulfonic acid) (ABTS) cation radical scavenging abilities of 14 (Z)-BPT derivatives were studied using Trolox as a positive control. The ABTS cation radical is formed by transferring one electron of ABTS to K2S2O8 (potassium persulfate). To obtain an ABTS cation radical solution, a 7.4 mM ABTS solution was mixed with 2.6 mM potassium persulfate and kept in the dark at room temperature for 24 h. After adjusting the absorbance of the ABTS cation radical solution to 0.700 ± 0.03 at 732 nm, the ABTS cation radical solution was mixed with test samples (1-14 and Trolox) to a final concentration of 100 µM and kept for 2 min in the dark. The radical scavenging effects were assessed based on the absorbance measured at 732 nm. Figure 9 shows the ABTS radical scavenging activity.

ABTS Cation Radical Scavenging Effects of (Z)-BPT Derivatives 1-14
It has been suggested that substances with antioxidant qualities may prevent the tyrosinase substrate (L-tyrosine and L-dopa) oxidation required for melanogenesis without interacting with the enzyme [46,47]. In addition, it has been reported that antioxidants that scavenge radicals may inhibit excess melanogenesis, because melanogenesis can be triggered by various radicals [48][49][50][51][52]. Therefore, the 2,2 -azino-bis(3-ethylbenzothiazoline-6-sulfonic acid) (ABTS) cation radical scavenging abilities of 14 (Z)-BPT derivatives were studied using Trolox as a positive control. The ABTS cation radical is formed by transferring one electron of ABTS to K 2 S 2 O 8 (potassium persulfate). To obtain an ABTS cation radical solution, a 7.4 mM ABTS solution was mixed with 2.6 mM potassium persulfate and kept in the dark at room temperature for 24 h. After adjusting the absorbance of the ABTS cation radical solution to 0.700 ± 0.03 at 732 nm, the ABTS cation radical solution was mixed with test samples (1-14 and Trolox) to a final concentration of 100 µM and kept for 2 min in the dark. The radical scavenging effects were assessed based on the absorbance measured at 732 nm. Figure 9 shows the ABTS radical scavenging activity.

DPPH Radical Scavenging Effects of (Z)-BPT Derivatives 1-14
The 2,2-diphenyl-1-picrylhydrazyl (DPPH) radical scavenging capacity of the (Z)-BPT derivatives (1 mM) was investigated using L-ascorbic acid (1 mM, vitamin C) as a positive control. Figure 10 shows the radical scavenging outcomes. A total of 7 of the 14 (Z)-BPT derivatives (1-5, 12, and 14) showed moderate to strong DPPH radical scavenging activity compared with L-ascorbic acid. While (Z)-BPT derivatives 8-11, with no hydroxyl substituents on the β-phenyl ring of the PUSC template, did not demonstrate any appreciable DPPH radical scavenging activity, those with hydroxyl substituents typically demonstrated moderate to strong radical scavenging activity. L-Ascorbic acid scavenged 97.1% of the DPPH radical, whereas (Z)-BPT derivative 1 with a 4-hydroxyl substituent on the β-phenyl ring scavenged 46.5% of the DPPH radical. Interestingly, when a 3bromo substituent was inserted into the β-phenyl ring of 1, the radical scavenging activity decreased to 14.8% (derivative 13), but when two bromo substituents were introduced at positions three and five of the β-phenyl ring of 1, the radical scavenging activity increased to 57.3% (derivative 14). The strongest radical scavenging activity was observed in derivative 3 with a β-3,4-dihydroxyphenyl moiety, which scavenged 96.2% of the DPPH radical, similar to L-ascorbic acid. Derivative 4, with a 4-hydroxy-3-methoxyl substituent, showed a strong DPPH radical scavenging activity of 81.0%, whereas derivative 6, in which the positions of the hydroxyl group and the methoxy groups on the β-phenyl ring of 4 were exchanged, showed a low radical scavenging activity of only 25.1%. Derivative 2, with a 2,4-dihydroxylphenyl ring, and derivative 5, with a 3-ethoxy-4-hydroxylphenyl ring, also showed good DPPH radical scavenging activities of 66.9 and 58.0%, respectively.

ROS Scavenging Effects of (Z)-BPT Derivatives 1-14
Owing to the association between antioxidant activity and anti-melanogenesis [46,47,53], the ROS scavenging effect of (Z)-BPT derivatives was evaluated. To assess the ROS scavenging activities of the (Z)-BPT derivatives, two assays were used to measure intracellular ROS and in vitro ROS generated by 3-morpholinosydnonimine (SIN-1).
The intracellular ROS scavenging activity was measured using the assay method reported by Bondy and Lebel [54] and Ali et al. [55]. We used 2′,7′-dichlorodihydrofluorescein diacetate (DCFH-DA) to measure intracellular ROS. DCFH-DA is converted to 2′,7′dichlorodihydrofluorescein (DCFH), a hydrolyzed form of DCFH-DA, which in turn, reacts with intracellular ROS to generate 2′,7′-dichlorofluorescein (DCF), a fluorescent material. B16F10 cells treated with 10 µM 3-morpholinosydnonimine (SIN-1), an ROS/RNS generator, greatly enhanced intracellular ROS levels ( Figure 11A). Of the fourteen (Z)-BPT derivatives, nine derivatives (1-6 and 12-14) at a concentration of 20 µM significantly decreased the intracellular ROS levels enhanced by SIN-1 treatment. In particular, five derivatives (3-5, 12 and 14) showed similar or more potent intracellular ROS scavenging activities than Trolox, a well-known antioxidant. These compounds have 3,4-dihydroxyl (catechol) (3), 3-alkoxy-4-hydroxyl (4 and 5), 3,5-dimethyl-4-hydroxyl (12), and 3,5-dibromo-4-hydroxyl (14) substituents on the β-phenyl ring, respectively. DCFH-DA, esterase, and SIN-1 were used to measure the scavenging activity of in vitro ROS generated by SIN-1. DCFH-DA reacts with esterase to generate DCFH, which in turn, reacts with the in vitro ROS generated by SIN-1 to produce DCF. To evaluate the in vitro ROS scavenging effects of the (Z)-BPT derivatives, a concentration of 40 µM was used. SIN-1 treatment greatly increased the ROS levels in vitro, and treatment with Trolox, a positive control, significantly decreased the enhanced ROS levels ( Figure 11B). As observed for Trolox, five (Z)-BPT derivatives (2)(3)(4)13, and 14) also significantly reduced the ROS levels that were increased by SIN-1 treatment to a similar or greater extent than did Trolox. Three derivatives, 2, 3, and 14, exerted strong scavenging abilities against both intracellular ROS and in vitro ROS. In particular, derivatives 2 and 3, which showed potent inhibitory activities against both tyrosinases in mushroom and B16F10 cells, had the most potent ROS scavenging activity.  DCFH-DA, esterase, and SIN-1 were used to measure the scavenging activity of in vitro ROS generated by SIN-1. DCFH-DA reacts with esterase to generate DCFH, which in turn, reacts with the in vitro ROS generated by SIN-1 to produce DCF. To evaluate the in vitro ROS scavenging effects of the (Z)-BPT derivatives, a concentration of 40 µM was used. SIN-1 treatment greatly increased the ROS levels in vitro, and treatment with Trolox, a positive control, significantly decreased the enhanced ROS levels ( Figure 11B). As observed for Trolox, five (Z)-BPT derivatives (2-4, 13, and 14) also significantly reduced the ROS levels that were increased by SIN-1 treatment to a similar or greater extent than did Trolox. Three derivatives, 2, 3, and 14, exerted strong scavenging abilities against both intracellular ROS and in vitro ROS. In particular, derivatives 2 and 3, which showed potent inhibitory activities against both tyrosinases in mushroom and B16F10 cells, had the most potent ROS scavenging activity.

Peroxynitrite (ONOO − ) Scavenging Effects of (Z)-BPT Derivatives 1-14
Reactive nitrogen species (RNS) have been reported to be involved in melanogenesis induction [56]. Thus, the peroxynitrite scavenging effects of the (Z)-BPT derivatives were examined using SIN-1, an RNS/ROS generator, and dihydrorhodamine (DHR123), an ROS indicator. SIN-1 generates peroxynitrite, which, in turn, oxidizes DHR123 to produce rhodamine 123, a fluorescent material. Therefore, the peroxynitrite scavenging activities of the (Z)-BPT derivatives 1-14 were determined at 50 µM by measuring the fluorescence of rhodamine 123. The peroxynitrite scavenging activity of the (Z)-BPT derivatives is shown in Figure 12.

Peroxynitrite (ONOO − ) Scavenging Effects of (Z)-BPT Derivatives 1-14
Reactive nitrogen species (RNS) have been reported to be involved in melanogenesis induction [56]. Thus, the peroxynitrite scavenging effects of the (Z)-BPT derivatives were examined using SIN-1, an RNS/ROS generator, and dihydrorhodamine (DHR123), an ROS indicator. SIN-1 generates peroxynitrite, which, in turn, oxidizes DHR123 to produce rhodamine 123, a fluorescent material. Therefore, the peroxynitrite scavenging activities of the (Z)-BPT derivatives 1-14 were determined at 50 µM by measuring the fluorescence of rhodamine 123. The peroxynitrite scavenging activity of the (Z)-BPT derivatives is shown in Figure 12.
Of the 14 derivatives, 13 significantly reduced peroxynitrite levels that were enhanced by treatment with SIN-1. Notably, eight derivatives (1-5 and 12-14), with a 4-hydroxyl substituent on the β-phenyl ring, showed potent peroxynitrite scavenging efficacy and were stronger peroxynitrite scavengers than penicillamine, which was used as the positive control. These results indicated that derivatives 2 and 3 are potent tyrosinase inhibitors with strong antioxidant capacities. Figure 12. Peroxynitrite scavenging effects of (Z)-BPT derivatives 1-14 using DHR123 (an ROS indicator) and SIN-1 (RNS/ROS generator). Penicillamine (PNC) was used as a positive reference material. PNC and test samples were used at a concentration of 50 µM and SIN-1 was used at a concentration of 10 µM. ### p < 0.001 for SIN-1-treated controls vs. non-treated controls; *** p < 0.001 and ** p < 0.01 for compound-treated groups vs. SIN-1-treated controls.

Effects of (Z)-BPT Derivatives 1-2 on Tyrosinase Expression
We investigated whether the derivatives 1 and 2 could influence the expression of tyrosinase protein using Western blotting. As shown in Figure 13, the tyrosinase expression levels in B16F10 cells were elevated via stimulator treatment (200 µM of IBMX plus 1 µM of α-MSH) by 6.6-fold for 1 and 5.7-fold for 2, compared to those in the untreated control. Treatment with derivatives 1 and 2 significantly reduced the tyrosinase expression levels in a dose-dependent manner. At 20 µM, derivatives 1 and 2 reduced the tyrosinase expression levels that increased during the stimulator treatment to 4.4-and 3.1- Figure 12. Peroxynitrite scavenging effects of (Z)-BPT derivatives 1-14 using DHR123 (an ROS indicator) and SIN-1 (RNS/ROS generator). Penicillamine (PNC) was used as a positive reference material. PNC and test samples were used at a concentration of 50 µM and SIN-1 was used at a concentration of 10 µM. ### p < 0.001 for SIN-1-treated controls vs. non-treated controls; *** p < 0.001 and ** p < 0.01 for compound-treated groups vs. SIN-1-treated controls.
Of the 14 derivatives, 13 significantly reduced peroxynitrite levels that were enhanced by treatment with SIN-1. Notably, eight derivatives (1-5 and 12-14), with a 4-hydroxyl substituent on the β-phenyl ring, showed potent peroxynitrite scavenging efficacy and were stronger peroxynitrite scavengers than penicillamine, which was used as the positive control. These results indicated that derivatives 2 and 3 are potent tyrosinase inhibitors with strong antioxidant capacities.

Effects of (Z)-BPT Derivatives 1-2 on Tyrosinase Expression
We investigated whether the derivatives 1 and 2 could influence the expression of tyrosinase protein using Western blotting. As shown in Figure 13, the tyrosinase expression levels in B16F10 cells were elevated via stimulator treatment (200 µM of IBMX plus 1 µM of α-MSH) by 6.6-fold for 1 and 5.7-fold for 2, compared to those in the untreated control. Treatment with derivatives 1 and 2 significantly reduced the tyrosinase expression levels in a dose-dependent manner. At 20 µM, derivatives 1 and 2 reduced the tyrosinase expression levels that increased during the stimulator treatment to 4.4-and 3.1-fold, respectively, implying that derivatives 1 and 2 may, in part, contribute to their anti-melanogenic actions by inhibiting tyrosinase expression. fold, respectively, implying that derivatives 1 and 2 may, in part, contribute to their antimelanogenic actions by inhibiting tyrosinase expression.

General Methods
Thermo Fisher Scientific (Carlsbad, CA, USA), Sigma-Aldrich (St. Louis, MO, USA), and SEJIN CI Co. (Seoul, Republic of Korea) provided the chemical reagents, which were all used without additional purification. All solvents requiring anhydrous conditions were distilled over CaH2 or Na/BP. All reactions were carried out under nitrogen, and the progress of the reactions were monitored via thin-layer chromatography (TLC) using Merck pre-coated 60F245 plates. Reaction mixture was purified via flash column chromatography using MP Silica 40-63 (60 Å). Low-resolution mass data were measured using the ESI positive or negative mode on an Expression CMS mass spectrometer (Advion, Ithaca, NY, USA). A Varian Unity AS500 unit (Agilent technologies, Santa Clara, CA, USA) was used to record the 1 H-NMR spectra at 500 MHz, and a Varian Unity INOVA 400 instrument (Agilent technologies, Santa Clara, CA, USA) was used to record the 1 H and 13 C-NMR spectra at 400 MHz and 100 MHz, respectively. Chloroform-d and dimethylsulfoxide-d6 were used as the NMR solvents. Parts per million (ppm) measurements of all chemical shifts were made in comparison to the corresponding residual solvent or deuterated peaks (δC 39.

General Methods
Thermo Fisher Scientific (Carlsbad, CA, USA), Sigma-Aldrich (St. Louis, MO, USA), and SEJIN CI Co. (Seoul, Republic of Korea) provided the chemical reagents, which were all used without additional purification. All solvents requiring anhydrous conditions were distilled over CaH 2 or Na/BP. All reactions were carried out under nitrogen, and the progress of the reactions were monitored via thin-layer chromatography (TLC) using Merck pre-coated 60F 245 plates. Reaction mixture was purified via flash column chromatography using MP Silica 40-63 (60 Å). Low-resolution mass data were measured using the ESI positive or negative mode on an Expression CMS mass spectrometer (Advion, Ithaca, NY, USA). A Varian Unity AS500 unit (Agilent technologies, Santa Clara, CA, USA) was used to record the 1 H-NMR spectra at 500 MHz, and a Varian Unity INOVA 400 instrument (Agilent technologies, Santa Clara, CA, USA) was used to record the 1 H and 13 C-NMR spectra at 400 MHz and 100 MHz, respectively. Chloroform-d and dimethylsulfoxide-d 6 were used as the NMR solvents. Parts per million (ppm) measurements of all chemical shifts were made in comparison to the corresponding residual solvent or deuterated peaks (δ C 39.8 and δ H 2.48 for DMSO-d 6 ; δ C 77.0 and δ H 7.27 for CDCl 3 ). The values for the coupling constant (J) are presented in hertz (Hz). For 1 H NMR, the following abbreviations were used: m, multiplet; brs, broad singlet; s, singlet; d, doublet; dd, doublet of doublets; t, triplet; and q, quartet.

General Procedure for the Synthesis of (Z)-BPT Derivatives 1-14
A solution of 16 (100 mg, 0.56 mmol) and an appropriate benzaldehyde (1.0 equiv.) in ethanol (3 mL) was refluxed in the presence of piperidine (0.02 mL, 0.17 mmol) as a base catalyst for 0.5-2.5 h. After cooling, water was added to the reaction mixture. The resultant solid was filtered and washed with water and ethanol for compounds 1, 2, 4-8, and 11-14,

Kinetic Studies on the Inhibition of Mushroom Tyrosinase by (Z)-BPT Derivatives 1-3
To determine the mode of action, Lineweaver-Burk plots of (Z)-BPT derivatives 1-3 were acquired. In brief, 10 µL of (Z)-BPT derivatives 1-3 (final concentrations:0, 3, 6, or 12 µM for 1; 0, 0.05, 0.1, or 0.2 µM for 2; 0, 2.5, 5, or 10 µM for 3) was added to a 96-well plate containing a substrate mixture (170 µL) consisting of an aqueous solution of L-dopa at final concentrations of 0.5, 1.0, 2.0, or 4.0 mM and 14.7 mM potassium phosphate buffer (pH 6.5), and a mushroom tyrosinase aqueous solution (20 µL, 20 units). Using a microplate reader (VersaMax TM , Molecular Devices, Sunnyvale, CA, USA) to measure the increase in absorbance at a wavelength of 492 nm (∆OD 492 /min), the initial rate of dopachrome synthesis in the reaction mixture was estimated. Using five different concentrations of L-dopa, Lineweaver-Burk plots (inverse of the response rate (1/V) vs. inverse of L-dopa concentration (1/(S)) were created. The convergence points of the four plot lines were used to define the mode of action of tyrosinase inhibition.

In Silico Study of Kojic Acid and (Z)-BPT Derivatives 1-3 on Mushroom Tyrosinase (mTYR)
The Schrödinger Suite (2021-1) was used for in silico analysis in accordance with the previously described procedures [35] with a few minor modifications. The Protein Preparation Wizard in Maestro12.4 was provided with the crystal structure of mTYR (PDB:2Y9X, Agaricus bisporus). The structure was then optimized by including hydrogen atoms, removing water molecules that were more than 3 Å away from the enzyme, and minimizing the structure. The binding site of tropolone was used to designate the glide grid of the enzyme. The structures of kojic acid and its derivatives 1-3 were imported into the CDXML format to the entry list of Maestro. The structures of kojic acid and its derivatives 1-3 were generated using LigPrep prior to ligand docking. Using glide from the task list [57], the ligand structures were docked to the glide grid of the enzyme. The glide extra precision (XP) approach was used to determine the binding affinity and interactions between ligands and proteins [58].

In Silico Study of Kojic Acid and (Z)-BPT Derivatives 1-3 on Human Tyrosinase Homology Model
To perform in silico docking simulations, the hTYR homology model previously created using the Swiss-Model online server and the Schrödinger Suite (2020-2) was used [34]. Kojic acid and (Z)-BPT derivatives 1-3 were docked using the hTYR homology model using the same procedures as those described above for mTYR docking.

B16F10 Murine Melanoma Cell Culture
B16F10 cells were purchased from the American Type Culture Collection (ATCC, Manassas, VA, USA). Penicillin, trypsin, streptomycin, phosphate-buffered solution (PBS), fetal bovine serum (FBS), and Dulbecco's modified Eagle's medium (DMEM) were purchased from Gibco/Thermo Fisher Scientific (Waltham, MA, USA). A DMEM solution containing penicillin-streptomycin solution (10,000 U/mL) and 10% heat-inactivated FBS was used to cultivate B16F10 cells at 37 • C in a humid environment with 5% CO 2 . These cells were cultivated in a 6-or 96-well culture plate, and cell viability, anti-tyrosinase activity, and anti-melanogenesis activity experiments were conducted.

Cytotoxicity Analysis of (Z)-BPT Derivatives 1-3 in B16F10 Melanoma Cells
Cytotoxicity studies of (Z)-BPT derivatives 1-3 were performed on B16F10 melanoma cells using the previously established EZ-Cytox (EZ-3000, DoGenBio, Seoul, Republic of Korea) test [33], and cytotoxicity studies of (Z)-BPT derivatives 1-3 were conducted on B16F10 melanoma cells. Briefly, B16F10 cells were seeded in a 96-well plate at a density of 1 × 10 4 cells/well and cultivated at 37 • C for 24 h in a humid environment with 5% CO 2 . Next day, the cells were treated with kojic acid or (Z)-BPT derivatives 1-3 at six different concentrations (0, 1, 2, 5, 10, and 20 µM) and incubated at 37 • C for 48 and 72 h, respectively, in a humid environment with 5% CO 2 . The cells were then incubated at 37 • C for 2 h after 10 µL of the EZ-Cytox solution had been added to each well. To determine the cell viability, the absorbance of each well was measured at 450 nm using a microplate reader (VersaMax TM , Molecular Devices, Sunnyvale, CA, USA). Each assay was performed in triplicates.
3.3.7. Anti-Melanogenesis Assay of Kojic Acid and (Z)-BPT Derivatives 1 and 2 in B16F10 Cells The melanin content assay standard procedure [37,59] was used to investigate the anti-melanogenic effect of (Z)-BPT derivatives 1 and 2 with a few minor modifications. Briefly, B16F10 melanoma cells were seeded in each well of a 6-well plate at a density of 1 × 10 5 cells/well and allowed to adhere to the well bottom under conditions that were identical to those used for cell culture. After 24 h, the cultured B16F10 cells were treated with kojic acid (20 µM) or four different concentrations (0, 5, 10, or 20 µM) of (Z)-BPT derivatives 1 and 2 for 1 h before stimulation with 200 µM IBMX and 1 µM α-MSH. For 72 h, B16F10 cells were incubated at 37 • C in a humid environment containing 5% CO 2 after being treated with IBMX and α-MSH. To determine the intracellular melanin content, B16F10 cells stimulated with IBMX and α-MSH for 72 h were rinsed twice with PBS and incubated in 200 µL 1N-NaOH solution including 10% dimethyl sulfoxide (DMSO) at 60 • C. The cell lysates were transferred to a 96-well plate after an hour of incubation, and the absorbance of the melanin content was calculated using a VersaMax TM microplate reader (Molecular Devices, Sunnyvale, CA, USA) to determine intracellular melanin levels. Each experiment was performed independently in triplicates.

Anti-Tyrosinase Activity Assay of Kojic Acid and (Z)-BPT Derivatives 1 and 2 in B16F10 Cells
The cellular tyrosinase activities of kojic acid and (Z)-BPT derivatives 1 and 2 were investigated using an anti-tyrosinase activity assay based on the oxidation rate of L-DOPA previously reported [60]. Briefly, B16F10 melanoma cells were seeded in a 6-well plate at a density of 1 × 10 5 cells/well and allowed to adhere to the bottom of each well under identical conditions to the cell culture. B16F10 cells were incubated for 24 h before treatment with kojic acid (20 µM) or various concentrations (0, 5, 10, or 20 µM) of (Z)-BPT derivatives 1 and 2. After incubation for 1 h, the tyrosinase activity of the cells was increased by treatment with 200 µM IBMX and 1 µM α-MSH. The cells were incubated for 72 h under identical conditions to cell culture before being rinsed twice with PBS, exposed to a 100 µL lysis buffer solution consisting of 90 µL phosphate buffer (50 mM) at pH 6.5, 5 µL PMSF (2 mM), and 5 µL Triton X-100 (20%), and then incubated for 30 min at −80 • C. The cell lysates were placed into microcentrifuge tubes after defrosting and they were centrifuged at 12,000 rpm at 4 • C for 30 min. The supernatants (80 µL) of the lysates were combined with 20 µL of 10 mM L-DOPA in a 96-well plate. Using a microplate reader used for the melanin content assay, the mixture's absorbances were measured every 10 min for 1 h at 37 • C and 492 nm. Each experiment was performed independently in triplicates.

ABTS Cation-Free Radical Scavenging Assay of (Z)-BPT Derivatives 1-14
The antioxidant properties of (Z)-BPT derivatives 1-14 were assessed using the ABTS cation-free radical scavenging assay, as previously described [61]. Potassium persulfate solution (10 mL, 2.45 mM) in distilled water was mixed with the ABTS solution (10 mL, 7 mM) in distilled water, and the mixture was incubated for 16 h at room temperature in the dark to produce ABTS cation-free radicals. Then, methanol was used to dilute the ABTS cation-free radical solution to an absorbance of 0.70 ± 0.02 at 734 nm. The test samples (10 µL; 90% EtOH + 10% DMSO) with a final concentration of 100 µM were added to the diluted ABTS cation-free radical solution (90 µL) and incubated in the dark at room temperature for 2 min. The absorbance of each well was measured at 734 nm using a VersaMax ™ microplate reader (Molecular Devices, Sunnyvale, CA, USA). The ABTS cation-free radical scavenging ability of the (Z)-BPT derivatives was determined using the following formula: ABTS cation radical scavenging activity (%) = (A C − A T ) × 100/A C (1) where A C is the absorbance of the untreated control, and A T is the absorbance of the tested compounds.
3.3.10. DPPH Radical Scavenging Assay of (Z)-BPT Derivatives 1-14 The antioxidant properties of (Z)-BPT derivatives 1-14 against DPPH (2,2-diphenyl-1-picrylhydrazyl) radicals were examined using a DPPH radical scavenging assay, as previously described with minor modifications [62]. DPPH solution (180 µL, 0.2 mM) in methanol was mixed with a (Z)-BPT derivative solution (1-14, 20 µL, 10 mM) in dimethyl sulfoxide in each well of a 96-well plate, and the mixture was kept for 30 min in the dark at room temperature. The absorbance at 517 nm was measured using a VersaMax ™ microplate reader (Molecular Devices, Sunnyvale, CA, USA). The capacity of the derivatives to scavenge DPPH radicals was compared with that of L-ascorbic acid, which was used as the positive control. Each experiment was run independently in triplicates.

Intracellular ROS Scavenging Activity Assay of (Z)-BPT Derivatives 1-14
The intracellular ROS scavenging activity was assessed as described by Lebel and Bondy [54] and Ali et al. [55]. The assay was evaluated using DCFH-DA as the ROSsensitive fluorescence probe. In brief, B16F10 melanoma cells at a density of 1 × 10 4 cells/well were incubated for 24 h in 96-well black plates and treated with 20 µM (Z)-BPT derivatives 1-14 for 2 h. Next, 10 µL of SIN-1 (20 µM in 50 mM sodium phosphate buffer, pH 7.4) was added for 1 h to stimulate B16F10 cells and induce ROS production. The cells were incubated with 20 µM DCFH-DA at 37 • C for 30 min. Fluorescence was measured every 5 min for 30 min using a microplate reader (Berthold Technologies GmbH & Co., Vienna, Austria) at 485 nm (excitation wavelength) and 535 nm (emission wavelength).
3.3.12. In Vitro ROS Scavenging Activity Assay of (Z)-BPT Derivatives 1-14 The ROS scavenging abilities of (Z)-BPT derivatives 1-14 and Trolox were assessed in accordance with methods previously reported by Lebel and Bondy [54]. Briefly, 50 mM phosphate buffer (pH 7.4) was dissolved with esterase (6 units/mL) and 1.25 mM DCFH-DA (2 ,7 -dichlorodihydrofluorescein diacetate) from Molecular Probes (Eugene, OR, USA), maintained at 37 • C for 30 min, and kept on ice in the dark until use. (Z)-BPT derivatives (final concentration: 40 µM; 10 µL) and SIN-1 (3-morpholinosydnonimine) (10 µM; 10 µL), dissolved in DMSO and phosphate buffer (180 µL), respectively, were added to each well of a black 96-well plate and kept in the dark for 30 min. The esterase-DCFH-DA mixture (50 µL) was then added to each well. Using a microplate reader (Berthold Advances GmbH & Co., Bad Wildbad, Germany), the fluorescence intensities of 2 ,7 -dichlorofluorescein (DCF or oxidized DCFH) were examined at excitation and emission wavelengths of 485 nm and 535 nm, respectively. Trolox was used as the reference control. The peroxynitrite scavenging activities were determined as described by Kooy et al. [63], with a few modifications. This technique measures the fluorescence of rhodamine 123, which is produced quickly when non-fluorescent DHR 123 (dihydrorhodamine 123) (Molecular Probes, Eugene, OR, USA) is used. Potassium chloride (5 mM), sodium chloride (90 mM), sodium phosphate (50 mM), and diethylenetriamine pentaacetic acid (DTPA, 100 µM) were used to prepare the pH 7.4 rhodamine buffer solution. DHR 123 at a final concentration of 5 µM was used. The buffer solution was prepared in advance for use in this test and was kept on ice. (Z)-BPT derivatives (10 µL), with a final concentration of 50 µM, were dissolved in DMSO and combined with SIN-1 (10 µL, 50 µM) and rhodamine buffer solution (180 µL) in each well of a black 96-well plate. The fluorescence intensities of the generated rhodamine 123 were measured using a fluorescent plate reader (Berthold Advances GmbH & Co., Bad Wildbad, Germany) at excitation and emission wavelengths of 485 and 535 nm, respectively. After deducting the background fluorescence from the final fluorescence intensity, the peroxynitrite scavenging activity was determined. L-penicillamine was used as a positive control, and the results are shown as the mean ± standard deviation.

Western Blot Assay of Tyrosinase Protein
After stimulating murine melanoma (B16F10) cells with α-MSH and IBMX for 24 h in the presence or absence of the tested derivatives at the indicated concentrations, the cells were lysed with lysis buffer containing protease and phosphatase inhibitors [38]. Protein concentrations were determined using a BCA protein assay kit, according to the manufacturer's instructions (Pierce, Rockford, IL, USA). Lysate proteins in each sample (10-20 µg) were subjected to 9% sodium dodecyl sulfate-polyacrylamide gel electrophoresis and electrophoretically transferred to a PVDF membrane (Millipore, Billerica, MA, USA) using a semi-dry system (Bio-Rad, Hercules, CA, USA). The membranes were incubated with 5% non-fat milk in TBST (0.1 M Tris-HCl, pH 7.5, 1.5 M NaCl, and 1% Tween 20) for 2 h. After blocking, the membranes were incubated with primary antibodies at 4 • C overnight. The following antibodies were used: tyrosinase and β-actin primary antibodies, as well as horseradish peroxidase-conjugated anti-mouse and anti-goat secondary antibodies (Santa Cruz Biotechnology, Dallas, TX, USA). The protein bands were visualized using the SuperSignal West Pico Chemiluminescence assay kit (Advansta, San Jose, CA, USA) and the Davinch-Chem program (Davinch-K, Seoul, Republic of Korea). The relative intensities were acquired using CS Analyzer 3.2 (Densitograph) image analysis software (http://www.attokorea.co.kr, accessed on 21 November 2022).

Statistical Analysis
The significance of the differences between treatment groups was determined using one-way analysis of variance (ANOVA) and the Bonferroni post hoc test. The statistical analysis was performed using GraphPad Prism 5 (La Jolla, CA, USA). All results are presented as the mean ± standard error of the mean (SEM). P-values with two sides less than 0.05 were regarded as statistically significant.

Conclusions
Novel anti-melanogenic compound (Z)-5-benzylidene-2-phenylthiazol-4(5H)-one ((Z)-BPT) derivatives were designed and synthesized from thiobenzamide and bromoacetic acid. Three (Z)-BPT derivatives (1-3) inhibited mushroom tyrosinase more effectively than kojic acid. Using kinetic studies, their modes of inhibition were demonstrated: 1 and 2 were competitive inhibitors, and 3 was a mixed-type inhibitor; the results were supported by docking simulation. Experiments using B16F10 cells revealed that derivatives 1 and 2 decreased the melanin production by inhibiting the tyrosinase activity and suppressing tyrosinase protein expression. Several (Z)-BPT derivatives exert strong antioxidant activities, scavenging DPPH and ABTS cation radicals, peroxynitrite, and intracellular and in vitro ROS.