Molecular Docking, Synthesis, and Tyrosinase Inhibition Activity of Acetophenone Amide: Potential Inhibitor of Melanogenesis

Department of Chemistry, Allama Iqbal Open University, Islamabad, Pakistan Faculty of Sciences, Department of Chemistry, University of Sialkot, Sialkot, Pakistan Department of Chemistry, University of Gujrat, Gujrat 50700, Pakistan Department of Zoology, University of Gujrat, Gujrat 50700, Pakistan Department of Physiology and Biochemistry, Faculty of Bio-Sciences, Cholistan University of Veterinary and Animal Sciences, Bahawalpur 63100, Pakistan Institute of Biochemistry, Biotechnology and Bioinformatics, Faculty of Science, The Islamia University of Bahawalpur, Bahawalpur 63100, Pakistan Department of Pharmacy, SBK Women University, Quetta, Pakistan Institute of Molecular Biology and Biotechnology (IMBB), University of Lahore, Punjab, Pakistan Department of Pharmacy, COMSATS University Islamabad, Abbottabad Campus, 22060, Pakistan


Introduction
A binuclear copper containing metalloenzyme tyrosinase has been detected in various species including fungi, bacteria, plants, and animals [1]. It is primarily involved in catalytic oxidation of its natural substrates L-tyrosine and L-DOPA (L-3,4-dihydroxyphenylalanine) to dopaquinone in a rate-limiting step during melanin biosynthesis which is associated with the reduction of oxygen (O 2 ) into water [2,3]. The dopaquinone synthesized in vivo in specialized organelles called melanosomes undergo further oxidation and several nonenzymatic polymeric reactions to form brown/black eumelanin and by reacting with cysteine and thiols to yellow/red pheomelanin [4][5][6]. Melanin is mainly responsible for skin, hair, and eye color and protects DNA degradation by ultraviolet (UV) radiation by acting as a photoprotectant in living organisms. The UV rays initiate melanogenesis for protecting skin cells from detrimental effects [7]. However, an abnormal production and accumulation of melanin in the dermal layer could cause several dermatological disorders including melasma, freckle, lentigo, and Riehl melanosis [8,9]. To date, numerous approaches including the tyrosinase inhibition, melanocyte to keratinocyte melanin transfer suppression, and cell signaling intervention for melanogenesis have been attempted to combat the increased amount of melanin [10][11][12][13][14]. Many natural and synthetic antimelanogenic substances with promising characteristics in both in vivo and in vitro have been reported [15,16], but the majority is associated with various side effects involving an in sufficient potency, permanent depigmentation and dermatitis [17,18], thyroid cancer [19], nephrotoxicity [20], genotoxicity [21], and toxicity to melanocytic cells [22]. Therefore, the search for novel inhibitors of melanogenesis is of valuable concern to both medicinal and cosmeceutical industries. The tyrosinase inhibitory activity of acetophenone moiety has only been rarely reported [23]. Nazir et al. and Rafiq et al. reported that hydroxyl-substituted benzoic/cinnamic acid analogues exhibited potent tyrosinase inhibition [24,25].
Over the past several years, our laboratory has reported various phenolic derivatives including tyramine, thymol, vanillin, carvacrol, and coumarin analogues which structurally mimic the natural substrates of the enzyme tyrosinase stronger than the L-tyrosine and L-DOPA [26][27][28][29][30]. We concluded from these investigations that substitution pattern of phenolic hydroxyls is directly related to degree of tyrosinase inhibition. Keeping in view these findings, the amide-based hydroxyl substituted benzoic/cinnamic acid analogues were synthesized to explore their tyrosinase inhibitory activity both in silico and in vitro so that more potent inhibitors are engineered. The synthesized benzoic 3(a-e) and cinnamic acids 5(a-e) derivatives passed the panassay interference screening except 3d [31]. The enzyme inhibitory kinetics and reversibility of the protein inhibitor complex were determined using Lineweaver-Burk and Dixon plots. The computational molecular docking was also performed against the crystal structure of target mushroom tyrosinase (PDB ID: 2Y9X) to predict the role of various amino acid residues in ligand-protein complex formation for further insights into biological properties.

Materials and Methods
2.1. General. Dry dichloromethane (DCM) was dried following the standard method. The IR-spectra (cm -1 ) were taken with a FT-IR spectrophotometer (Perkin Elmer, USA) using attenuated total reflection sampling technique (ATR). NMR spectra were taken using a Bruker spectrometer ( 1 H, 400 MHz, 13 C, 100 MHz, DMSO-d 6 ), and chemical shifts (δ) are reported in parts per million (ppm). The Flash elemental analyzer connected with thermal conductivity detector (Flash-TCD 2000 Series, USA) was used to perform elemental analysis (C, H) and reported uncorrected with ±0.3%. The melting points determined by Digimelt MPA 160 (USA) are reported uncorrected. Silica gel thin layer column chromatography-TLC was used to check the purity of title derivatives (petroleum spirit : ethyl acetate, 2 : 1). All other reagents such as the enzyme mushroom tyrosinase (Mtyr, EC 1.14.18.1), L-DOPA, and 1-(3-acetylphenyl)ethan-1-one were analytical grade (Sigma Aldrich, USA) and are used without further purification.  [32], was condensed with various hydroxyl-substituted benzoic acids 2(a-e) and cinnamic acids 4(a-e) (0.01 mol) with equimolar triethylamine and potassium iodide (KI) in dimethylformamide (DMF) (25 mL) under nitrogen at 25°C and stirred overnight. The reaction mixture was then concentrated using rotary evaporator and extracted using ethyl acetate (3 × 25 mL). Finally, the combined ethyl acetate layer was treated with HCl (5%) and brine followed by drying over magnesium sulphate MgSO 4 , filtered and concentrated by rotary evaporator to get the crude title products 3(a-e) and 5(a-e) (Schemes 1 and 2) and further purified by normal phase column chromatography (petroleum spirit : ethyl acetate, 3 : 1).

The Spectroscopic Characterization of Synthesized
Compounds 3(a-e) and 5(a-e). The 1 H-NMR of final compounds 3(a-e) and 5(a-e) showed an overlap of peaks in the alkyl region for homologues of this series. The singlet of terminal ketonic methyl proton of peripheral phenyl ring occupies the most upfield region of the spectrum encompassing the range of 2.54-2.56 ppm. In 13 C NMR, the carbon of this terminal ketonic methyl group resonates between chemical shifts 26.73 and 26.80 ppm. The most shielded protons in the peripheral chain were the methylene (-CH 2 ) proton which are directly linked to ester groups. This methylene group has a chemical shift value of 4.67-4.94 ppm while the carbon shift of this group in 13 Scheme 1: Synthetic pathway for title derivatives 3(a-e); reaction conditions and reagents; (i) dichloromethane (CH 2 Cl 2 )/triethylamine (C 2 H 5 ) 3 N, 0-5°C, 5 hours stirring (ii) dimethylformamide (DMF)/(C 2 H 5 ) 3 N/KI, stirring for 24 hours at room temperature.    2.3. Mushroom Tyrosinase Inhibitory Activity. In vitro antityrosinase potential was determined following the method previously described with little modifications [33,34]. Briefly, 140 μL of sodium phosphate buffer, 20 μL (30 U/mL) of the enzyme mushroom tyrosinase, and 20 μL of the tested compound were added in a well of a 96-well microplate. After a preincubation of 10 minutes at room temperature of 20 μL (0.85 mM), the substrate L-DOPA was added followed by an incubation of 20 minutes at the same temperature. Subsequently, the absorbance of the intermediate dopachrome was recorded by using an Opti Max tunable microplate reader (Sunnyvale, USA) at 472 nm. The standard inhibitors were kojic acid and arbutin. The negative inhibitor was a phosphate buffer without an inhibitor. All the compounds were dissolved in DMSO, and blank tests were performed with DMSO. IC 50 values were calculated by Origin (8.6, 64-bit).

Mushroom Tyrosinase Inhibitory
Kinetics. Based upon the antityrosinase assay, the most potent derivative (5c) was chosen for kinetic inhibition studies using enzyme and L-DOPA as a substrate [35][36][37]. Various concentrations of the inhibitor 5c were among 0.00 to 0.013 μM and that of L-DOPA was 0.00 to 0.0064 μM during all kinetic assays. For preincubation and calculations, the same protocol was followed as described in in vitro tyrosinase inhibition assay and continuous monitoring at 30 s interval for 5 minutes after adding the enzyme tyrosinase, and the formation of dopachrome was measured at 475 nm. The enzyme

Retrieval of Mushroom Tyrosinase in Maestro.
For this purpose, the crystal structure of target protein (PDB ID 2Y9X) was retrieved from RCSB and prepared in "Protein Preparation Wizard" workflow in Maestro Schrödinger for adding hydrogens, adjusting protonation and bond orders appropriate for pH 7. The water molecules beyond 5 Å from het groups were removed and minimized the prepared protein structure for further grid generation and docking analysis.

Receptor Grid Generation and Docking.
For grid generation, the catalytic pocket is selected from its cocrystallized ligand and literature [38]. The grid was generated by specifying the cocrystallized ligand tropolone of the active site of the target protein. The receptor grid box was defined as 20 Å box. After grid preparation, Glide_dock_XP precision docking experiment was performed with default docking setup parameters reporting the 15 top-ranked poses per ligand [39]. The predicted binding scores (binding energies) and proper orientation of ligands within catalytic region of tyrosinase was also performed. Finally, the most favorable binding mode of active compounds within the binding pocket was investigated in terms of docking score and 3D graphical images of the binding pose of the best docked score was also done with Maestro (Schrödinger).

Results and Discussions
3.1. In Vitro Mushroom Tyrosinase Inhibition Assay. Benzoic and cinnamic acid analogues 3(a-e) and 5(a-e) were designed and synthesized to evaluate in vitro tyrosinase inhibition activity. Table 1 demonstrates that title amide derivatives bearing a cinnamic acid moiety showed better tyrosinase inhibition than reference drugs kojic acid and arbutin and benzoic acid 3(a-e) analogues. The cinnamic acid analogues 5a, 5b, 5c, and 5d have different substitutions on the phenyl ring. The analogue 5a is an unsubstituted cinnamic acid analogue, and 5b and 5c contain 4-hydroxyl and 2,4-dihydroxyl substituted cinnamic acid moieties,  5 BioMed Research International respectively, whereas the inhibitor 5d bears a chloro substituent at the para position. The hydroxyl substitution at the second and fourth positions of the cinnamic acid phenyl ring in 5c resulted in a considerable rise in the enzyme inhibitory activity ( Table 1). The phenolic hydroxyls play a crucial role in tyrosinase inhibition as the natural substrates of tyrosinase, L-tyrosine, and L-DOPA also bear the phenolic hydroxyls [40]. A naturally occurring cinnamic acid has been reported with extensive physiological activities including tyrosinase inhibitory action [41]. The compounds 5a (IC 50 46.43 μM), 5b (IC 50 309.20 μM), and 5d (IC 50 61.62 μM) displayed lower inhibitory activity than 5c (IC 50 0.0020 μM) compared to reference inhibitor kojic acid (IC 50 16.69 μM) and arbutin (IC 50 191.17 μM). It is reported here that the ortho-and parahydroxyls of the phenyl ring in 5c appear to show an increased tyrosinase inhibitory activity ( Figure 1).   BioMed Research International ing site of mushroom tyrosinase was determined using Michaelis-Menten kinetic studies. The inhibitors exhibit a dose-dependent inhibition of the enzyme tyrosinase. A striking reduction in the reaction rate in the presence of inhibitors refers to a decrease in the final absorbance in comparison to control without inhibitors. The inhibitory index of all these tested compounds varied depending upon the position of various substituents and classes of compounds [27]. Inhibition kinetics was analyzed by the Lineweaver-Burk and Dixon plots. The kinetics analyses of 1/V max vs. 1/[S] at different doses of 5c reveal that the Michaelis-Menten constant (K m ) changes while that of 1/V max remained the same representing the competitive nature of the most potent inhibitor (5c) (Figure 2(a)). The dissociation constant K i for 5c was 0.0072 μM calculated by Dixon plots as shown in Figure 2(b) [27,42]. showed a group of straight lines with varying slopes and intersecting on the same point on the y-axis preferring the reversible effect of 5c on tyrosinase ( Figure 3). The results suggested that 5c effectively inhibits the enzyme tyrosinase by binding to its active site reversibly [35].

In Silico Studies
3.3.1. Binding Pocket Analysis of Analogues 3(a-e) and 5(a-e) against Mushroom Tyrosinase. A computational molecular docking study was also conducted to examine the binding conformations of all the synthesized compounds within the catalytic pocket of enzyme tyrosinase. The docked ligandprotein complexes were investigated based on docking score (kcal/mol) and the hydrophobic/philic bonding interaction pattern. The docking scores had little fluctuations, and the comparison depicted no significant energy difference among all docked molecules due to similar basic skeleton of the ligands. Therefore, majority of the ligands showed efficient docking energy values. From docking results and in vitro enzyme inhibitory assay, the most active compound 5c (-6.568 kcal/mol) was visualized to determine its interactions in the catalytic site of the protein tyrosinase ( Figure 4). Two

Chemoinformatics and Lipinski's Rule of Five (RO5).
The Molsoft tool used to predict pharmacokinetics properties such as the number of hydrogen bond acceptor/donor, solubility, polar surface area, and molar volume. The standard range for molar molecular weight (160-500), polarizability ≤ 5, and solubility ≤ 4 has been reported in literature [44]. The computational results predicted that the molecular weight of 5c is 355.11 g/mol, and it possesses 6 HBA (≤10), 3 HBD (≤5), 2.72 LogP (≤5), -3.11 mol/L LogS, and 91.64 A 2 PSA, and the molar volume of 355.34 A 3 proved to have a drug-like behavior within the standard range (Table 2). Furthermore, Lipinski's rule of five (RO5) states nothing about structural features found in drugs and nondrugs for their medicinal potential.

Data Availability
Data could be provided upon request from the shared corresponding author Dr. Zaman Ashraf, mzchem@yahoo.com

Conflicts of Interest
The authors declare that they have no known conflict of interest to influence the work reported in this research.