Tyrosinase is a tetrameric enzyme that plays an important role in pigment production. Overproduction of melanin, which may lead to several skin disorders, is a result of tyrosinase activity. Hence, tyrosinase inhibitors are of key importance in the treatment of these disorders. In the present study, four flavonoid inhibitors, namely chrysin, naringin, quercetin, and kaempferol, were evaluated physiochemically, and the inhibitory effects of these compounds on tyrosinase activity were evaluated using the molecular dynamics (MD) simulation method. To create the best conformation of the enzyme-substrate/inhibitor, the docking process for enzyme-substrate, i.e., enzymechrysin, enzyme-quercetin, enzyme-naringin, and enzyme-kaempferol, was performed. The complexes with the best binding energies were selected as the models for the MD simulation process. Furthermore, the structural (RMSD, Rg, RMSF, and Distance) and the thermodynamics properties of the complexes were evaluated. Additionally, the PMF was conducted to calculate the binding free energies. The results showed that chrysin, quercetin and the substrate were at similar distances to the amino acids of the active site, but naringin and kaempferol were closer to the active site of the enzyme than the substrate. Moreover, the analysis of the binding energy revealed that the substrates, chrysin, kaempferol, quercetin, and naringin bound to the enzyme with binding energies of !7.8, !3.1, !7.1, !3.9, and !8.4 kcal/mol, respectively, which confirms that naringin has the highest inhibitory effect on tyrosinase among other inhibitors, which makes it an appropriate candidate as a whitening agent in skin disorders.
