Journal of Chemistry

Journal of Chemistry / 2013 / Article

Research Article | Open Access

Volume 2013 |Article ID 849782 | https://doi.org/10.1155/2013/849782

Ramesh L. Sawant, Prashant D. Lanke, Jyoti B. Wadekar, "Tyrosinase Inhibitory Activity, 3D QSAR, and Molecular Docking Study of 2,5-Disubstituted-1,3,4-Oxadiazoles", Journal of Chemistry, vol. 2013, Article ID 849782, 7 pages, 2013. https://doi.org/10.1155/2013/849782

Tyrosinase Inhibitory Activity, 3D QSAR, and Molecular Docking Study of 2,5-Disubstituted-1,3,4-Oxadiazoles

Academic Editor: Romdhane Karoui
Received19 Jan 2012
Revised15 May 2012
Accepted20 May 2012
Published01 Jul 2012

Abstract

In continuation with our research program, in search of potent enzyme tyrosinase inhibitor, a series of synthesized 2,5-disubstituted 1,3,4-oxadiazoles have been evaluated for enzyme tyrosinase inhibitory activity. Subsequently, 3D QSAR and docking studies were performed to find optimum structural requirements for potent enzyme tyrosinase inhibitor from this series. The synthesized 20 compounds of 2,5-disubstituted-1,3,4-oxadiazole series were screened for mushroom tyrosinase inhibitory activity at various concentrations by enzyme inhibition assay. The percentage enzyme inhibition was calculated by recording absorbance at 492 nm with microplate reader. 3D QSAR and docking studies were performed using VLife MDS 3.5 software. In the series 2,5-disubstituted-1,3,4-oxadiazoles enzyme tyrosinase inhibitory activity was found to be dose dependent with maximum activity for compounds 4c, 4h, 4m, and 4r. 3D QSAR and docking studies revealed that more electropositive and less bulky substituents if placed on 1,3,4-oxadiazole nucleus may result in better tyrosinase inhibitory activity in the series.

1. Introduction

Tyrosinase (E.C. 1.14.18.1), known as polyphenol oxidase (PPO), is a multifunctional, glycosylated, and copper-containing enzyme from the oxidase superfamily widely distributed in microorganisms, plants, and animals [1]. It is a nonessential amino acid made by the body that is a building block for several important neurotransmitters like epinephrine, norepinephrine, serotonin, and dopamine. It plays a key role in melanin biosynthesis catalyzing two divergent reactions: the hydroxylation of monophenols (cresolase or monophenolase activity) and the oxidation of o-diphenols (catecholase or diphenolase activity) into reactive o-quinones [24]. It is involved in human neuromelanin formation in the substantia nigra of the brain and dopamine neurotoxicity and contributes to Parkinson’s disease-related neurodegeneration. In insects, tyrosinase is uniquely associated with different physiologically important biochemical processes, including sclerotization of the insect cuticle, defensive encapsulation and melanization of foreign organisms, and wound healing [510].

Tyrosinase inhibitors have become increasingly important in the field of medicine, agriculture, and cosmetics [11]. Presently, tyrosinase inhibitors have been a great concern due to its role in both mammalian melanogenesis and fruit or fungi enzymatic browning. Although melanin has photoprotective function in human skin, the accumulation of an abnormal amount of melanin in different specific parts of the skin resulting in more pigmented patches (melasma, freckles, ephelides, senile lentigines, etc.) might become an esthetic problem [12, 13]. In addition, enzymatic browning in fruit and fungi is undesirable, for example, fresh fruits, beverages, vegetables, and mushrooms which decrease the commercial value of the products [14]. These phenomena prompted us to design novel tyrosinase inhibitors with higher bioactivities that could be useful in skin whitening and antibrowning of foods.

In continuation with our earlier work, here we report enzyme tyrosinase inhibitory activity, three-dimensional quantitative structure activity relationship (3D QSAR) analysis, and molecular docking study with enzyme tyrosinase of series of 2,5-disubstituted-1,3,4-oxadiazoles synthesized and reported by us [15].

2. Materials and Methods

2.1. Tyrosinase Inhibitory Activity

As an extension of our work, synthesized 20 analogues of 2,5-disubstituted-1,3,4-oxadiazole series (Table 1) were used to study their enzyme tyrosinase inhibitory activity. Tyrosinase inhibition activity was determined by the method described earlier [16]. The title compounds were dissolved in DMSO. The 96 well plate was prepared by applying 140 μL of phosphate buffer (pH 6.8), 20 μL of mushroom tyrosinase (48 units/mL), 20 μL of sample, and 20 μL of DL-DOPA (0.85 mM) to make the concentrations of 10, 20, 40, 80, and 100 μg/mL. After incubation for 10 minutes, the enzyme activity was determined by measuring the absorbance at 492 nm using the microplate reader (Biotek Company, US). Kojic acid (1 mg/mL) was used as positive control. The percentage of tyrosinase inhibition was calculated as follows: where, and are absorbance at 492 nm with and without test sample, respectively. The percentage tyrosinase inhibition at various concentrations of title compounds was calculated (Table 2), and logarithmic value of percentage inhibition at 100 μg (Log % ) was used as a dependant variable for the development of valid 3D-QSAR models.


Serial numberCompoundStructure

1.4a849782.tab.001a
2.4b849782.tab.001c
3.4c849782.tab.001e
4.4d849782.tab.001g
5.4e849782.tab.001i
6.4f849782.tab.001k
7.4g849782.tab.001m
8.4h849782.tab.001o
9.4i849782.tab.001q
10.4j849782.tab.001s
11.4k849782.tab.001b
12.4l849782.tab.001d
13.4m849782.tab.001f
14.4n849782.tab.001h
15.4o849782.tab.001j
16.4p849782.tab.001l
17.4q849782.tab.001n
18.4r849782.tab.001p
19.4s849782.tab.001r
20.4t849782.tab.001t


CompoundPercentage inhibition of enzyme tyrosinase
10 μg20 μg40 μg80 μg100 μg

Kojic acid7.8615.2228.9542.5254.01
4a2.714.025.387.2115.21
4b2.613.627.5110.128.91
4c4.3817.8925.3935.4140.51
4d4.9115.6723.7732.1836.13
4e5.0312.3820.6128.4138.11
4f2.715.547.8110.0212.02
4g2.576.918.2712.1417.61
4h6.5718.3728.1634.5939.71
4i4.7110.6218.7128.7235.78
4j3.6815.1221.9127.8333.81
4k2.276.9210.6112.1818.62
4l3.1112.6220.6825.5132.56
4m3.9815.6828.1937.7641.14
4n2.5110.1612.7216.2222.16
4o2.978.9115.7422.0228.34
4p4.1815.7117.8122.6131.57
4q3.388.6220.6326.5134.62
4r4.6115.7628.6336.4447.92
4s4.0111.8820.8126.6131.54
4t4.6210.8115.4220.9833.02

2.2. Computational Details
2.2.1. Geometry Optimization

The 3D QSAR studies of 2,5-disubstituted-1,3,4-oxadiazole derivatives were carried out using VLife molecular design suite software version 3.5 (MDS 3.5) running on Pentium IV processor. Three-dimensional structures of title compounds were constructed and optimized their geometries to make the conformations with least potential energy using merck molecular force field (MMFF) and MMFF charge for the atom [17] followed by considering distance-dependent dielectric constant of 1.0 and convergence criteria (rms gradient) of 0.01 kcal/mol.

2.2.2. Alignment of Molecules

All molecules in the data set are aligned by template-based method where a template is built by considering common substructures in the series (Figure 1). Highly bioactive energetically stable conformation in the series is chosen as a reference molecule (Figure 2) on which other molecules are aligned (Figure 3).

2.2.3. Activity Prediction

The predictability of the QSAR model would be good if the values of biological activity predicted by the QSAR model do not appreciably differ from the observed results of biological activity for the given data set. Quality of selected models was further ascertained by , , and test. The models were cross-validated by “leave one out” scheme and cross-validation corelation coefficient was calculated. The model with high value is said to have high predictability.

2.2.4. 3D QSAR

Several 3D QSAR techniques such as comparative molecular field analysis (COMFA), comparative molecular similarity analysis (COMSIA), and k-nearest neighbor (kNN) [18] are being used in modern QSAR research. In the present study, molecular field analysis coupled with partial least squares (PLS) was applied to obtain a 3D QSAR model [19]. The calculated steric and electrostatic field descriptors were used as independent variables and log % values were used as dependent variables to derive the 3D QSAR models.

2.2.5. Docking

The crystal structure of mushroom tyrosinase (1wx2) was obtained from protein data bank and water molecules in the crystal structure were deleted. The optimized receptor was then saved as mol file and used for docking simulation.

2.2.6. Ligand Preparation

The 2D structures of the compounds were built and then converted into the 3D. The 3D structures were then energetically minimized up to the rms gradient of 0.01 using MMFF.

2.2.7. Identification of Cavities

By using cavity determination option of software, cavities of enzyme were determined. The cavities in the receptor were mapped to assign an appropriate active site. The basic feature used to map the cavities were the surface mapping of the receptor and identifying the geometric voids as well as scaling the void for its hydrophobic characteristics. Hence, all the cavities that are present in receptor are identified and ranked based on their size and hydrophobic surface area. Considering the dimensions and hydrophobic surface area, cavity-1 with volume 219.531250 ų and surface area 340.284790 Ų was found to be the best void as an active site.

2.2.8. Scoring Function

Distinction of good or bad docked conformation is based on scoring or fitness function. MDS uses fitness functions on only electrostatic and both steric and electrostatic interactions between receptor ligand as well as dock score scoring function. The dock score compute binding affinity of a given protein-ligand complex with known 3D structure.

3. Results and Discussion

3.1. Tyrosinase Inhibition Activity

The percentage inhibition of enzyme tyrosinase by compounds 4a4t is estimated at concentrations of 10, 20, 40, 80, and 100 μg and observed to be dose dependent. The compounds 4c, 4h, 4m, and 4r shows maximum where compounds 4b, 4f, 4g, and 4k shows minimum inhibition of enzyme tyrosinase from the series.

3.2. 3D QSAR

The model 1 describes the optimum structural feature for the tyrosinase inhibition activity as shown in Table 3. The training set of 16 molecules and test set of 4 molecules were used as described earlier. S_793, S_757, E_550 are the steric and electrostatic field energy of interactions between probe (CH3) and compounds at their corresponding spatial grid points of 793, 757, and 550 as shown in Figure 4. Low residual value observed in case of training as well as test set (Table 4 and Figure 5) confirms reliability of the QSAR model 1. The selected descriptors have good correlation with biological activity and a low intercorrelation among them depicted in Table 5. As per model 1, steric descriptor S_793 with negative coefficients represent less bulky substituent is favourable; S_757 with positive coefficients represent more bulky substituent is favourable at 2, 5-disubstituted oxadiazole nucleus while electrostatic field descriptor E_550 with positive coefficient indicate electropositive groups is favourable for activity. Compound 4r in the dataset bearing less bulky and electropositive chloro(hyroxy)phenyl rings on oxadiazole nucleus is most active, whereas less electropositive substituted oxadiazole nucleus in compound 4b is least active which indicates that developed 3D QSAR model satisfies model requirements.


Model numberEquation 𝑛 𝑟 2 𝑞 2 𝐹 test

1.Log % 𝐼 1 0 0 𝜇 g = −0.0005 − 147.0670 (±10.4994)
S_793 + 0.0183 (±0.0028) S_757 + 0.0339 (±0.0074) E_550
160.89200.818733.0329


CompoundDescriptorsActivity
S_793S_757E_550ActualPredictedResidue

** 4a−0.0063117.920552.0790881.1821291.325233−0.1431
4b−0.0055710.99283−0.79360.9498780.992882−0.043
4c−0.0074320.208072.1512231.6075621.5335830.073979
4d−0.0074320.208072.1512231.5578681.5335830.024285
4e−0.00651301.9273361.5810391.5698360.011203
4f−0.0057419.12609−1.805361.0799041.131727−0.05182
4g−0.0062918.094683.1323911.2457591.361143−0.11538
** 4h−0.0057410.46793−0.807661.59891.0075270.591373
** 4i−0.0074419.835724.2671131.553641.600056−0.04642
4j−0.0074419.835724.2671131.5290451.600056−0.07101
4k−0.0065430−5.242571.269981.331729−0.06175
4l−0.007630−5.366611.5126841.4838610.028823
4m−0.00607302.5703451.6142641.5273440.08692
4n−0.00696.3836753.1568991.345571.237790.10778
** 4o−0.0077217.788551.6607191.45241.516313−0.06391
4p−0.0077217.788551.6607191.4992751.516313−0.01704
4q−0.00673301.1708051.5393271.577306−0.03798
4r−0.00787301.0068011.6805171.73882−0.0583
4s−0.0063730−2.434341.8525541.8059040.097768
4t−0.0069430−2.224931.7826291.6943160.025586

**Test set compounds.

log % 𝐼 1 0 0 𝜇 g E_550S_757S_793

log % 𝐼 1 0 0 𝜇 g 1.000000
E_5500.2563601.000000
S_7570.403622−0.3200181.000000
S_793−0.46893−0.372145−0.0664911.000000

3.3. Docking

Docking studies of the title compounds with enzyme tyrosinase yielded docking score ranging from −5.8169 to −4.7232 (Table 6) indicating stable enzyme-substrate interactions. The compound 4r with best docking score (−5.8169) has shown highest inhibition of enzyme tyrosinase in the dataset. The compounds interact with enzyme tyrosinase by binding with SER142, SER146, GLY145, and GLY198 amino acid residues as shown in Figure 6. The tyrosinase inhibitory assay and molecular docking study reveals that 2,5-disubstituted-1,3,4-oxadiazoles inhibits enzyme tyrosinase may be because of its interaction with SER142, SER146, GLY145, and GLY198 amino acid residues.


Sr. no.CompoundDocking score
(Kcal/mol)

14a−4.8032
24b−4.9936
34c−5.2345
44d−5.3241
54e−5.1800
64f−4.7232
74g−4.8539
84h−5.2711
94i−5.3660
104j−5.2453
114k−5.3106
124l−5.4441
134m−5.0124
144n−4.9321
154o−5.4969
164p−5.4867
174q−5.3207
184r−5.8169
194s−5.1214
204t−5.1318

4. Conclusion

It can be concluded from the whole study that electropositivity is crucial for the inhibition of enzyme tyrosinase by 2,5-disubstituted-1,3,4-oxadiazoles. It would be worthwhile to synthesize a novel 1,3,4-oxadiazole analogues with less bulky and more electropositivity group as a better enzyme tyrosinase inhibitors.

References

  1. Y. J. Kima and H. Uyamab, “Tyrosinase inhibitors from natural and synthetic sources: structure, inhibition mechanism and perspective for the future,” Cellular and Molecular Life Sciences, vol. 62, no. 15, pp. 1707–1723, 2005. View at: Publisher Site | Google Scholar
  2. F. G. Cánovas, F. García-Carmona, J. V. Sánchez, J. L. Iborra-Pastor , and J. A. Lozano-Teruel, “The role of pH in the melanin biosynthesis pathway,” The Journal of Biological Chemistry, vol. 257, no. 15, pp. 8738–8744, 1982. View at: Google Scholar
  3. J. N. Rodríguez-López, J. Tudela, R. Varón, and F. García-Cánovas, “Kinetic study on the effect of pH on the melanin biosynthesis pathway,” Biochimica et Biophysica Acta, vol. 1076, no. 3, pp. 379–386, 1991. View at: Google Scholar
  4. C. J. Cooksey, P. J. Garratt, E. J. Land et al., “Evidence of the indirect formation of the catecholic intermediate substrate responsible for the autoactivation kinetics of tyrosinase,” The Journal of Biological Chemistry, vol. 272, no. 42, pp. 26226–26235, 1997. View at: Publisher Site | Google Scholar
  5. J. M. Pawelek and A. M. Korner, “The biosynthesis of mammalian melanin,” American Scientist, vol. 70, no. 2, pp. 136–145, 1982. View at: Google Scholar
  6. A. M. Mayer, “Polyphenol oxidases in plants-recent progress,” Phytochemistry, vol. 26, no. 1, pp. 11–20, 1986. View at: Google Scholar
  7. J. R. Whitaker, “Polyphenol oxidase,” in Food Enzymes: Structure and Mechanism, D. W. S. Wong, Ed., vol. 271, Chapman and Hall, New York, NY, USA, 1995. View at: Google Scholar
  8. M. Friedman, “Food browning and its prevention: an overview,” Food Chemistry, vol. 44, no. 3, pp. 631–653, 1996. View at: Publisher Site | Google Scholar
  9. F. M. Barrett, “Wound-healing phenoloxidase in larval cuticle of Calpodes ethlius (Lepidoptera: Hesperiidae),” Canadian Journal of Zoology, vol. 62, no. 5, pp. 834–838, 1984. View at: Google Scholar
  10. M. Sugumaran, “Molecular mechanisms for cuticular sclerotization,” Advances in Insect Physiology, vol. 21, pp. 179–231, 1988. View at: Publisher Site | Google Scholar
  11. K. Maeda and M. Fukuda, “In vitro effectiveness of several whitening cosmetic components in human melanocytes,” Journal of the Society of Cosmetic Chemists, vol. 42, no. 6, pp. 361–368, 1991. View at: Google Scholar
  12. T. S. Chang, “An updated review of tyrosinase inhibitors,” International Journal of Molecular Sciences, vol. 10, no. 6, pp. 2440–2475, 2009. View at: Publisher Site | Google Scholar
  13. M. T. H. Khan, “Molecular design of tyrosinase inhibitors: a critical review of promising novel inhibitors from synthetic origins,” Pure and Applied Chemistry, vol. 79, no. 12, pp. 2277–2295, 2007. View at: Publisher Site | Google Scholar
  14. F. Artes, M. Castaner, and M. I. Gil, “Review: enzymatic browning in minimally processed fruit and vegetables,” Food Science and Technology International, vol. 4, no. 6, pp. 377–389, 1998. View at: Google Scholar
  15. R. L. Sawant and P. D. Lanke, “Microwave assisted synthesis and 3D QSAR analysis of analgesic oxadiazoles,” International Journal of Drug Design and Discovery, vol. 1, no. 4, pp. 336–344, 2010. View at: Google Scholar
  16. S. Tadtong, A. Viriyaroj, S. Vorarat, S. Nimkulrat, and S. Suksamrarn, “Antityrosinase and antibacterial activities of mangosteen pericarp extract,” Journal of Health Research, vol. 23, no. 2, pp. 99–102, 2009. View at: Google Scholar
  17. T. A. Halgren, “Merck molecular force field. III. Molecular geometries and vibrational frequencies for MMFF94,” Journal of Computational Chemistry, vol. 17, no. 5-6, pp. 553–586, 1996. View at: Google Scholar
  18. S. Ajmani, K. Jadhav, and S. A. Kulkarni, “Three-dimensional QSAR using the k-nearest neighbor method and its interpretation,” Journal of Chemical Information and Modeling, vol. 46, no. 1, pp. 24–31, 2006. View at: Publisher Site | Google Scholar
  19. M. S. Bhatia, K. B. Ingale, P. B. Choudhari, R. L. Sawant, and C. R. Patil, “Two- and three-dimensional quantitative structure-activity relationships studies on a series of diuretics,” Latin American Journal of Pharmacy, vol. 28, no. 6, pp. 927–931, 2009. View at: Google Scholar

Copyright © 2013 Ramesh L. Sawant et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.


More related articles

 PDF Download Citation Citation
 Download other formatsMore
 Order printed copiesOrder
Views2662
Downloads1580
Citations

Related articles

We are committed to sharing findings related to COVID-19 as quickly as possible. We will be providing unlimited waivers of publication charges for accepted research articles as well as case reports and case series related to COVID-19. Review articles are excluded from this waiver policy. Sign up here as a reviewer to help fast-track new submissions.