Antitumor activity of a pexidartinib bioisostere inhibiting CSF1 production and CSF1R kinase activity in human hepatocellular carcinoma.

Macrophage colony-stimulating factor (M-CSF, also known as CSF1) in tumor tissues stimulates tumor growth and tumor-induced angiogenesis through an autocrine and paracrine action on CSF1 receptor (CSF1R). In the present study, novel bioisosteres of pexidartinib (1) were synthesized and evaluated their inhibitory activities against CSF1R kinase and tumor growth. Among newly synthesized bioisosteres, compound 3 showed the highest inhibition (95.1%) against CSF1R tyrosine kinase at a fixed concentration (1 µM). The half maximal inhibitory concentration (IC50) of pexidartinib (1) and compound 3 was 2.7 and 57.8 nM, respectively. Unlike pexidartinib (1), which cross-reacts to three targets with structural homology, such as CSF1R, c-KIT, and FLT3, compound 3 inhibited CSF1R, c-KIT, but not FLT3, indicating compound 3 may be a more selective CSF1R inhibitor than pexidartinib (1). The inhibitory effect of compound 3 on the proliferation of various cancer cell lines was the strongest in U937 cells followed by THP-1 cells. In the case of cancer cell lines derived from solid tumors, the anti-proliferative activity of compound 3 was weaker than pexidartinib (1), except for Hep3B. However, compound 3 was safer than pexidartinib (1) in terminally differentiated normal cells such as macrophages. Pexidartinib (1) and compound 3 suppressed the production of CSF1 in Hep3B liver cancer cells as well as in the co-culture of Hep3B cells and macrophages. Also, pexidartinib (1) and compound 3 decreased the population ratio of the M2/M1 phenotype and inhibited their migration. Importantly, compound 3 preferentially inhibited M2 phenotype over M1, and the effect was about 4 times greater than that of pexidartinib (1). In addition, compound 3 inhibited maintenance of cancer stem cell population. In a chick chorioallantoic membrane (CAM) tumor model implanted with Hep3B cells, tumor growth and tumor-induced angiogenesis were significantly blocked by compound 3 to a similar extent as pexidartinib (1). Overall, compound 3, a bioisostere of pexidartinib, is an effective dual inhibitor to block CSF1R kinase and CSF1 production, resulting in significant inhibition of tumor growth.

Hydroxamic Acid as a Potent Metal-Binding Group for Inhibiting Tyrosinase.

Tyrosinase, a metalloenzyme containing a dicopper cofactor, plays a central role in synthesizing melanin from tyrosine. Many studies have aimed to identify small-molecule inhibitors of tyrosinase for pharmaceutical, cosmetic, and agricultural purposes. In this study, we report that hydroxamic acid is a potent metal-binding group for interacting with dicopper atoms, thereby inhibiting tyrosinase. Hydroxamate-containing molecules, including anticancer drugs targeting histone deacetylase, vorinostat and panobinostat, significantly inhibited mushroom tyrosinase, with inhibitory constants in the submicromolar range. Of the tested molecules, benzohydroxamic acid was the most potent. Its inhibitory constant of 7 nM indicates that benzohydroxamic acid is one of the most potent tyrosinase inhibitors. Results from differential scanning fluorimetry revealed that direct binding mediates inhibition. The enzyme kinetics were studied to assess the inhibitory mechanism of the hydroxamate-containing molecules. Experiments with B16F10 cell lysates confirmed that the new inhibitors are inhibitory against mammalian tyrosinase. Docking simulation data revealed intermolecular contacts between hydroxamate-containing molecules and tyrosinase.

Discovery of Kinase and Carbonic Anhydrase Dual Inhibitors by Machine Learning Classification and Experiments.

A multi-target small molecule modulator is advantageous for treating complicated diseases such as cancers. However, the strategy and application for discovering a multi-target modulator have been less reported. This study presents the dual inhibitors for kinase and carbonic anhydrase (CA) predicted by machine learning (ML) classifiers, and validated by biochemical and biophysical experiments. ML trained by CA I and CA II inhibitor molecular fingerprints predicted candidates from the protein-specific bioactive molecules approved or under clinical trials. For experimental tests, three sulfonamide-containing kinase inhibitors, 5932, 5946, and 6046, were chosen. The enzyme assays with CA I, CA II, CA IX, and CA XII have allowed the quantitative comparison in the molecules' inhibitory activities. While 6046 inhibited weakly, 5932 and 5946 exhibited potent inhibitions with 100 nM to 1 µM inhibitory constants. The ML screening was extended for finding CAs inhibitors of all known kinase inhibitors. It found XMU-MP-1 as another potent CA inhibitor with an approximate 30 nM inhibitory constant for CA I, CA II, and CA IX. Differential scanning fluorimetry confirmed the direct interaction between CAs and small molecules. Cheminformatics studies, including docking simulation, suggest that each molecule possesses two separate functional moieties: one for interaction with kinases and the other with CAs.

Repositioning Food and Drug Administration-Approved Drugs for Inhibiting Biliverdin IXß Reductase B as a Novel Thrombocytopenia Therapeutic Target.

Biliverdin IXß reductase B (BLVRB) has recently been proposed as a novel therapeutic target for thrombocytopenia through its reactive oxygen species (ROS)-associated mechanism. Thus, we aim at repurposing drugs as new inhibitors of BLVRB. Based on IC50 (<5 µM), we have identified 20 compounds out of 1496 compounds from the Food and Drug Administration (FDA)-approved library and have clearly mapped their binding sites to the active site. Furthermore, we show the detailed BLVRB-binding modes and thermodynamic properties (ΔH, ΔS, and KD) with nuclear magnetic resonance (NMR) and isothermal titration calorimetry together with complex structures of eight water-soluble compounds. We anticipate that the results will serve as a novel platform for further in-depth studies on BLVRB effects for related functions such as ROS accumulation and megakaryocyte differentiation, and ultimately treatments of platelet disorders.

New Inhibitors of Laccase and Tyrosinase by Examination of Cross-Inhibition between Copper-Containing Enzymes.

Coppers play crucial roles in the maintenance homeostasis in living species. Approximately 20 enzyme families of eukaryotes and prokaryotes are known to utilize copper atoms for catalytic activities. However, small-molecule inhibitors directly targeting catalytic centers are rare, except for those that act against tyrosinase and dopamine-ß-hydroxylase (DBH). This study tested whether known tyrosinase inhibitors can inhibit the copper-containing enzymes, ceruloplasmin, DBH, and laccase. While most small molecules minimally reduced the activities of ceruloplasmin and DBH, aside from known inhibitors, 5 of 28 tested molecules significantly inhibited the function of laccase, with the Ki values in the range of 15 to 48 µM. Enzyme inhibitory kinetics classified the molecules as competitive inhibitors, whereas differential scanning fluorimetry and fluorescence quenching supported direct bindings. To the best of our knowledge, this is the first report on organic small-molecule inhibitors for laccase. Comparison of tyrosinase and DBH inhibitors using cheminformatics predicted that the presence of thione moiety would suffice to inhibit tyrosinase. Enzyme assays confirmed this prediction, leading to the discovery of two new dual tyrosinase and DBH inhibitors.

Novel Pyridine Bioisostere of Cabozantinib as a Potent c-Met Kinase Inhibitor: Synthesis and Anti-Tumor Activity against Hepatocellular Carcinoma.

Two novel bioisosteres of cabozantinib, 3 and 4, were designed and synthesized. The benzene ring in the center of the cabozantinib structure was replaced by trimethylpyridine (3) and pyridine (4), respectively. Surprisingly, the two compounds showed extremely contrasting mesenchymal-epithelial transition factor (c-Met) inhibitory activities at 1 µM concentration (4% inhibition of 3 vs. 94% inhibition of 4). The IC50 value of compound 4 was 4.9 nM, similar to that of cabozantinib (5.4 nM). A ligand-based docking study suggested that 4 includes the preferred conformation for the binding to c-Met in the conformational ensemble, but 3 does not. The anti-proliferative activity of compound 4 against hepatocellular carcinoma (Hep3B and Huh7) and non-small-cell lung cancer (A549 and H1299) cell lines was better than that of cabozantinib, whereas 3 did not show a significant anti-proliferative activity. Moreover, the tumor selectivity of compound 4 toward hepatocellular carcinoma cell lines was higher than that of cabozantinib. In the xenograft chick tumor model, compound 4 inhibited Hep3B tumor growth to a much greater extent than cabozantinib. The present study suggests that compound 4 may be a good therapeutic candidate against hepatocellular carcinoma.

Synthesis and anti-hepatocellular carcinoma activity of aminopyridinol-sorafenib hybrids.

Sorafenib is recommended as the primary therapeutic drug for patients with hepatocellular carcinoma. To discover a new compound that avoids low response rates and toxic side effects that occur in sorafenib therapy, we designed and synthesized new hybrid compounds of sorafenib and 2,4,5-trimethylpyridin-3-ols. Compound 6 was selected as the best of 24 hybrids that inhibit each of the four Raf kinases. The anti-proliferative activity of 6 in HepG2, Hep3B, and Huh7 cell lines was slightly lower than that of sorafenib. However, in H6c7 and CCD841 normal epithelial cell lines, the cytotoxicity of 6 was much lower than that of sorafenib. In addition, similar to sorafenib, compound 6 inhibited spheroid forming ability of Hep3B cells in vitro and tumour growth in a xenograft tumour model of the chick chorioallantoic membrane implanted with Huh7 cells. Compound 6 may be a promising candidate targeting hepatocellular carcinoma with low toxic side effects on normal cells.

Structural Deformation of MTX Induced by Nanodrug Conjugation Dictate Intracellular Drug Transport and Drug Efficacy.

BACKGROUND: Understanding structural interactions between the active drug and conjugated nanoparticles is critical for optimizing intracellular drug transport and for increasing nano drug efficacy. In this regard, analyzing the conformational deformation of conjugated drugs surrounding nanoparticles is essential to understand the corresponding nanodrug efficacy. PURPOSE: The objective of this study is to present an optimal synthesis method for efficient drug delivery through a clear structural analysis of nanodrugs according to the type of conjugation. METHODS AND RESULTS: In this study, the structural variation of methotrexate (MTX) surrounding carbon nanotubes, depending on the type of conjugation style, such as covalent and non-covalent (PEGylation) bonds, was investigated. Specifically, covalent bonds of MTX surrounding CNTs induced greater structural deformation compared to non-covalent bonds (ie, PEGylated CNT). CONCLUSION: Greater changes in the structural variations of MTX analyzed by nuclear magnetic resonance (NMR) significantly improved the anti-inflammatory drug efficacy of human fibroblast-like synovial cells (FLS) via stable drug release in the extracellular environment and burst drug release under intracellular conditions.

X-ray Crystal Structure-Guided Design and Optimization of 7H-Pyrrolo[2,3-d]pyrimidine-5-carbonitrile Scaffold as a Potent and Orally Active Monopolar Spindle 1 Inhibitor.

Triple-negative breast cancer (TNBC) is an aggressive breast-cancer subtype associated with poor prognosis and high relapse rates. Monopolar spindle 1 kinase (MPS1) is an apical dual-specificity protein kinase that is over-expressed in TNBC. We herein report a highly selective MPS1 inhibitor based on a 7H-pyrrolo[2,3-d]pyrimidine-5-carbonitrile scaffold. Our lead optimization was guided by key X-ray crystal structure analysis. In vivo evaluation of candidate (9) is shown to effectively mitigate human TNBC cell proliferation.

Tubulosine selectively inhibits JAK3 signalling by binding to the ATP-binding site of the kinase of JAK3.

Gain- or loss-of-function mutations in Janus kinase 3 (JAK3) contribute to the pathogenesis of various haematopoietic malignancies and immune disorders, suggesting that aberrant JAK3 signalling is an attractive therapeutic target to treat these disorders. In this study, we performed structure-based computational database screening using the 3D structure of the JAK3 kinase domain and the National Cancer Institute diversity set and identified tubulosine as a novel JAK3 inhibitor. Tubulosine directly blocked the catalytic activity of JAK3 by selective interacting with the JAK3 kinase domain. Consistently, tubulosine potently inhibited persistently activated and interleukin-2-dependent JAK3, and JAK3-mediated downstream targets. Importantly, it did not affect the activity of other JAK family members, particularly prolactin-induced JAK2/signal transducer and activator of transcription 5 and interferon alpha-induced JAK1-TYK2/STAT1. Tubulosine specifically decreased survival and proliferation of cancer cells, in which persistently active JAK3 is expressed, by inducing apoptotic and necrotic/autophagic cell death without affecting other oncogenic signalling. Collectively, tubulosine is a potential small-molecule compound that selectively inhibits JAK3 activity, suggesting that it may serve as a promising therapeutic candidate for treating disorders caused by aberrant activation of JAK3 signalling.

Sulfisoxazole inhibits the secretion of small extracellular vesicles by targeting the endothelin receptor A.

Inhibitors of the secretion of cancer exosomes, which promote cancer progression and metastasis, may not only accelerate exosome biology research but also offer therapeutic benefits for cancer patients. Here we identify sulfisoxazole (SFX) as an inhibitor of small extracellular vesicles (sEV) secretion from breast cancer cells through interference with endothelin receptor A (ETA). SFX, an FDA-approved oral antibiotic, showed significant anti-tumor and anti-metastatic effects in mouse models of breast cancer xenografts, the reduced expression of proteins involved in biogenesis and secretion of sEV, and triggered co-localization of multivesicular endosomes with lysosomes for degradation. We demonstrate the important role of ETA, as target of SFX, by gain- and loss-of-function studies of the ETA protein, through a direct binding assay, and pharmacological and genetic approaches. These findings may provide a foundation for sEV-targeted cancer therapies and the mechanistic studies on sEV biology.

The inhibitory potential of Broussochalcone A for the human cytochrome P450 2J2 isoform and its anti-cancer effects via FOXO3 activation.

BACKGROUND: Broussonetia papyrifera (L.) Ventenat, a traditional medicinal herb, has been applied as a folk medicine to treat various diseases. Broussochalcone A (BCA), a chalcone compound isolated from the cortex of Broussonetia papyrifera (L.) Ventenat, exhibits several biological activities including potent anti-oxidant, antiplatelet, and cytotoxic effects. PURPOSE: The purpose of this study is to elucidate the inhibitory effect of BCA against CYP2J2 enzyme which is predominantly expressed in human tumor tissues and carcinoma cell lines. STUDY DESIGN: The inhibitory effect of BCA on the activities of CYP2J2-mediated metabolism were investigated using human liver microsomes (HLMs), and its anti-cancer effect against human hepatoma HepG2 cells was also evaluated. METHODS: Two representative CYP2J2-specific probe substrates, astemizole and ebastine, were incubated in HLMs with BCA. After incubation, the samples were analyzed using liquid chromatography-tandem mass spectrometry. To investigate the binding model between BCA and CYP2J2, we carried out structure-based docking simulations by using software and scripts written in-house. RESULTS: BCA inhibited CYP2J2-mediated astemizole O-demethylation and ebastine hydroxylase activities in a concentration dependent manner with Ki values of 2.3 and 3.7 µM, respectively. It also showed cytotoxic effects against human hepatoma HepG2 cells in a dose-dependent manner with activation of apoptosis related proteins. CONCLUSION: Overall, this was the first report of the inhibitory effects of BCA on CYP2J2 in HLMs. The present data suggest that BCA is a potential candidate for further evaluation for its CYP2J2 targeting anti-cancer activities.

4-Hydroxynonenal-induced GPR109A (HCA2 receptor) activation elicits bipolar responses, Gαi -mediated anti-inflammatory effects and Gßγ -mediated cell death.

BACKGROUND AND PURPOSE: In this study, we examined the possibility that 4-hydroxynonenal (4-HNE) acting as a ligand for the HCA2 receptor (GPR109A) elicits both anti-inflammatory and cell death responses. EXPERIMENTAL APPROACH: Agonistic activity of 4-HNE was determined by observing the inhibition of cAMP generation in CHO-K1-GPR109A-Gi cell line, using surface plasmon resonance (SPR) binding and competition binding assays with [3 H]-niacin. 4-HNE-mediated signalling pathways and cellular responses were investigated in cells expressing GPR109A and those not expressing these receptors. KEY RESULTS: Agonistic activity of 4-HNE was stronger than that of niacin or 3-OHBA at inhibiting forskolin-induced cAMP production and SPR binding affinity. In ARPE-19 and CCD-841 cells, activation of GPR109A by high concentrations of the agonists 4-HNE (≥10 µM), niacin (≥1000 µM) and 3-OHBA (≥1000 µM) induced apoptosis accompanied by elevated Ca2+ and superoxide levels. This 4-HNE-induced cell death was blocked by knockdown of GPR109A or NOX4 genes, or treatment with chemical inhibitors of Gßγ (gallein), intracellular Ca2+ (BAPTA-AM), NOX4 (VAS2870) and JNK (SP600125), but not by the cAMP analogue 8-CPT-cAMP. By contrast, low concentrations of 4-HNE, niacin and 3-OHBA down-regulated the expression of pro-inflammatory cytokines IL-6 and IL-8. These 4-HNE-induced inhibitory effects were blocked by a cAMP analogue but not by inhibitors of Gßγ -downstream signalling molecules. CONCLUSIONS AND IMPLICATIONS: These results revealed that 4-HNE is a strong agonist for GPR109A that induces Gαi -dependent anti-inflammatory and Gßγ -dependent cell death responses. Moreover, the findings indicate that specific intracellular signalling molecules, but not GPR109A, can serve as therapeutic targets to block 4-HNE-induced cell death.

A targeted ferritin-microplasmin based thrombolytic nanocage selectively dissolves blood clots.

The use of thrombolytic therapies is limited by an increased risk of systemic hemorrhage due to lysis of hemostatic clots. We sought to develop a plasmin-based thrombolytic nanocage that efficiently dissolves the clot without causing systemic fibrinolysis or disrupting hemostatic clots. Here, we generated a double chambered short-length ferritin (sFt) construct that has an N-terminal region fused to multivalent clot targeting peptides (CLT: CNAGESSKNC) and a C-terminal end fused to a microplasmin (µPn); CLT recognizes fibrin-fibronectin complexes in clots, µPn efficiently dissolves clots, and the assembly of double chambered sFt (CLT-sFt-µPn) into nanocage structure protects the activated-µPn from its circulating inhibitors. Importantly, activated CLT-sFt-µPn thrombolytic nanocage showed a prolonged circulatory life over activated-µPn and efficiently lysed the preexisting clots in both arterial and venous thromboses models. Thus, CLT-sFt-µPn thrombolytic nanocage platform represents the prototype of a targeted clot-busting agent with high efficacy and safety over existing thrombolytic therapies.

Structural ensemble-based docking simulation and biophysical studies discovered new inhibitors of Hsp90 N-terminal domain.

Heat shock protein 90 (Hsp90) is one of the most abundant cellular proteins and plays a substantial role in the folding of client proteins. The inhibition of Hsp90 has been regarded as an attractive therapeutic strategy for treating cancer because many oncogenic kinases are Hsp90 client proteins. In this study, we report new inhibitors that directly bind to N-terminal ATP-binding pocket of Hsp90. Optimized structure-based virtual screening predicted candidate molecules, which was followed by confirmation using biophysical and cell-based assays. Among the reported crystal structures, we chose the two structures that show the most favourable early enrichments of true-positives in the receiver operating characteristic curve. Four molecules showed significant changes in the signals of 2D [1H, 15N] correlation NMR spectroscopy. Differential scanning calorimetry analysis supported the results indicating direct binding. Quantified dissociation constant values of the molecules, determined by a series of 2D NMR experiments, lie in the range of 0.1-33 µM. Growth inhibition assay with breast and lung cancer cells confirmed the cellular activities of the molecules. Cheminformatics revealed that the molecules share limited chemical similarities with known inhibitors. Molecular dynamics simulations detailed the putative binding modes of the inhibitors.

Thiopurine Drugs Repositioned as Tyrosinase Inhibitors.

Drug repositioning is the application of the existing drugs to new uses and has the potential to reduce the time and cost required for the typical drug discovery process. In this study, we repositioned thiopurine drugs used for the treatment of acute leukaemia as new tyrosinase inhibitors. Tyrosinase catalyses two successive oxidations in melanin biosynthesis: the conversions of tyrosine to dihydroxyphenylalanine (DOPA) and DOPA to dopaquinone. Continuous efforts are underway to discover small molecule inhibitors of tyrosinase for therapeutic and cosmetic purposes. Structure-based virtual screening predicted inhibitor candidates from the US Food and Drug Administration (FDA)-approved drugs. Enzyme assays confirmed the thiopurine leukaemia drug, thioguanine, as a tyrosinase inhibitor with the inhibitory constant of 52 µM. Two other thiopurine drugs, mercaptopurine and azathioprine, were also evaluated for their tyrosinase inhibition; mercaptopurine caused stronger inhibition than thioguanine did, whereas azathioprine was a poor inhibitor. The inhibitory constant of mercaptopurine (16 µM) was comparable to that of the well-known inhibitor kojic acid (13 µM). The cell-based assay using B16F10 melanoma cells confirmed that the compounds inhibit mammalian tyrosinase. Particularly, 50 µM thioguanine reduced the melanin content by 57%, without apparent cytotoxicity. Cheminformatics showed that the thiopurine drugs shared little chemical similarity with the known tyrosinase inhibitors.

Interaction between human angiogenin and the p53 TAD2 domain and its implication for inhibitor discovery.

Interaction between angiogenin and the p53 TAD2 domain in cancer cells can inhibit the function of the p53 tumor suppressor and promote cell survival. Based on a model structure using NMR and mutational analysis, positively charged 31 RRR33 and 50 KRSIK54 motifs of human angiogenin were identified as p53-binding sites that could interact with negatively charged D48/E51 and E56 residues of the p53 TAD2 domain, respectively. These results suggest that 31 RRR33 and 50 KRSIK54 motifs of human angiogenin might play a critical role in the regulation of p53-mediated apoptosis and angiogenesis in cancer cells. This study identifies potential target sites for screening angiogenin-specific inhibitors that could not only inhibit p53 binding but could also simultaneously inhibit cell binding, internalization, DNA binding, and nuclear translocation of human angiogenin.

Antithrombotic properties of JJ1, a potent and novel thrombin inhibitor.

The development of new anticoagulants is an important goal for the improvement of thrombosis treatment. Recent studies have suggested the importance of thrombin inhibitors in the modulation of thromboembolic disorders. The aim of this study was to discover a new small-molecule thrombin inhibitor. In this study, the compound JJ1, which has a novel scaffold, was selected by structure-based docking simulation to determine its potential inhibitory activity against thrombin. JJ1 was shown to inhibit the catalytic activity of human α-thrombin with a K i of 0.019 µM by direct binding to the active site and with at least 10,000-fold selectivity relative to that reported for the inhibition of other biologically important serine proteases. JJ1 prolonged clotting times (activated partial thromboplastin time and prothrombin time) and inhibited the activity and production of thrombin. Furthermore, it inhibited thrombin-catalyzed fibrin polymerization and platelet aggregation. Similar to its in vitro antithrombotic activities, JJ1 showed enhanced antithrombotic effects in an in vivo pulmonary embolism and arterial thrombosis model. It also exhibited anticoagulant effects in mice. Collectively, these results demonstrated that JJ1 was a potent, direct, and selective thrombin inhibitor that may be useful in the management of various thrombotic disorders.

Ligand-Mediated Folding of the OmpA Periplasmic Domain from Acinetobacter baumannii.

The periplasmic domain of OmpA from Acinetobacter baumannii (AbOmpA-PD) binds to diaminopimelate and anchors the outer membrane to the peptidoglycan layer in the cell wall. Although the crystal structure of AbOmpA-PD with its ligands has been reported, the mechanism of ligand-mediated folding of AbOmpA remains elusive. Here, we report that in vitro refolded apo-AbOmpA-PD in the absence of ligand exists as a mixture of two partially folded forms in solution: mostly unfolded (apo-state I) and hololike (apo-state II) states. Binding of the diaminopimelate or glycine ligand induced complete folding of AbOmpA-PD. The apo-state I was highly flexible and contained some secondary structural elements, whereas the apo-state II closely resembled the holo-state in terms of both structure and backbone dynamics, except for the ligand-binding region. 15N-relaxation-dispersion analyses for apo-state II revealed substantial motion on a millisecond timescale of residues in the H3 helix near the ligand-binding site, with this motion disappearing upon ligand binding. These results provide an insight into the ligand-mediated folding mechanism of AbOmpA-PD in solution.

Ensemble-Based Virtual Screening Led to the Discovery of New Classes of Potent Tyrosinase Inhibitors.

In this study, we report new classes of potent tyrosinase inhibitors identified by enhanced structure-based virtual screening prediction; the enzyme and melanin content assays were also confirmed. Tyrosinase, a type-3 copper protein, participates in two distinct reactions, hydroxylation of tyrosine to DOPA and conversion of DOPA to dopaquinone, in melanin biosynthesis. Although numerous inhibitors of this reaction have been reported, there is a lag in the discovery of the new functional moieties. In order to improve the performance of virtual screening, we first produced an ensemble of 10,000 structures using molecular dynamics simulation. Quantum mechanical calculation was used to determine the partial charges of catalytic copper ions based on the met and deoxy states. Second, we selected a structure showing an optimal receiver operating characteristic (ROC) curve with known direct binders and their physicochemically matched decoys. The structure revealed more than 10-fold higher enrichment at 1% of the ROC curve than those observed in X-ray structures. Third, high-throughput virtual screening with DOCK 3.6 was performed using a library consisting of approximately 400,000 small molecules derived from the ZINC database. Fourth, we obtained the top 60 molecules and tested their inhibition of mushroom tyrosinase. The extended assays included 21 analogs of the 21 initial hits to test their inhibition properties. Here, the moieties of tetrazole and triazole were identified as new binding cores interacting with the dicopper catalytic center. All 42 inhibitors showed inhibitory constant, Ki, values ranging from 11.1 nM and 33.4 µM, with a tetrazole compound exhibiting the strongest activity. Among the 42 molecules, five displayed more than 30% reduction in melanin production when treated in B16F10 melanoma cells; cell viability was >90% at 20 µM. Particularly, a thiosemicarbazone-containing compound reduced melanin content by 55%.