Open Access Article
This Open Access Article is licensed under a
Creative Commons Attribution 3.0 Unported Licence

Understanding the inhibitory mechanism of tea polyphenols against tyrosinase using fluorescence spectroscopy, cyclic voltammetry, oximetry, and molecular simulations

Haifeng Tangace, Fengchao Cui*a, Haijuan Lib, Qingrong Huangd and Yunqi Li*a
aKey Laboratory of Synthetic Rubber, Changchun Institute of Applied Chemistry (CIAC), Chinese Academy of Sciences, Changchun, Jilin 130022, P. R. China. E-mail: fccui@ciac.ac.cn; yunqi@ciac.ac.cn
bState Key Laboratory of Electroanalytical Chemistry, Changchun Institute of Applied Chemistry (CIAC), Chinese Academy of Sciences, Changchun, Jilin 130022, P. R. China
cSchool of Life Science, Jilin University, Changchun, Jilin 130012, P. R. China
dDepartment of Food Science, Rutgers University, 65 Dudley Road, New Brunswick, NJ 08901, USA
eUniversity of Chinese Academy of Sciences, Beijing 100049, P. R. China

Received 24th November 2017 , Accepted 15th February 2018

First published on 22nd February 2018


Abstract

Inhibiting the activity of tyrosinase is a very effective and safe way to prevent enzymatic browning in food and to resist pests in agriculture. Tea polyphenols (TPs), regarded as safe and non-toxic food additives, have been reported due to their potential inhibitory capability against tyrosinase, but their ambiguous inhibitory mechanisms have severely limited their application. In the present work, fluorescence spectroscopy, cyclic voltammetry (CV), oximetry and molecular simulation approaches were employed to shed light on the underlying inhibitory mechanisms of TPs with different structures including (+)-catechin, (−)-epicatechin gallate (ECG) and (−)-epigallocatechin gallate (EGCG) against tyrosinase. Fluorescence spectra show that the three TPs are capable of binding tyrosinase with a molar proportion of 1[thin space (1/6-em)]:[thin space (1/6-em)]1. The analysis of CV curves and oxygen utilization suggests that these three TPs can be oxidized by tyrosinase, revealing that these three TPs are suicide inhibitors of tyrosinase. Furthermore, ECG and catechin make tyrosinase irreversibly inactivated due to their catechol group (ring B) being catalyzed by tyrosinase through a cresolase-like pathway, while EGCG inhibits the activity of tyrosinase by competing with or delaying the oxidation of substrate. Molecular simulations further confirm that ring B of ECG and catechin makes a significant contribution to tyrosinase inhibitory activities, and has a direct interaction with the coupled binuclear copper ions in the optimal orientation required by the cresolase-like pathway.


1. Introduction

Tyrosinase (EC 1.14.18.1) containing a catalytic center with bi-copper ions coordinated by six histidine residues, can catalyze the oxidation of both monophenol (cresolase or monophenolase activity) and diphenol (catecholase or diphenolase activity) to the corresponding ortho-quinones, followed by self-polymerizing to dark melanin or another pigment.1–4 In these oxidation reactions, the redox forms of the catalytic center include Cu2–O2 with both monophenolase and diphenolase activities, Cu2–O with diphenolase activity and Cu2 without any activity, according to the amount of coordinated oxygen.1 In the food industry, tyrosinase can expedite the enzymatic browning of post-harvest produce by catalyzing the oxidation of polyphenol compounds,4 leading to the loss of nutritional quality and the shortening of storage life.5 Moreover, the higher tyrosinase activity would adversely affect the pest control in agriculture due to promoting the physiological processes of larval maturity.6,7 Inhibiting tyrosinase bioactivity is therefore favorable to the preservation of agricultural product and pest prevention. Although a plenty of natural tyrosinase inhibitors have been discovered, the understanding of the inhibition mechanism of tyrosinase inhibitors is still deficient, which have severely limited their application in food industry and agriculture and the development of new inhibitors.8–11

Flavonoids, widely used as dietary supplements to functional foods,12 are the largest groups in natural tyrosinase inhibitors which have been discovered,13–17 up to now. Flavonoids were usually subdivided into six major groups including flavanols, flavones, flavonols, flavanones, isoflavones and anthocyanidins.18 The position and number of hydroxyl group on functional moiety of flavonoids have a significant effect on their inhibitory activity and associated mechanism. Some flavonols with a 3-hydroxyl-4-keto moiety, such as kaempferol, quercetin and morin, exhibit the inhibitory activity toward tyrosinase by chelating with the copper ions and occupying the catalytic center,15,19 but most of flavonols have weaker inhibitory activity than kojic acid.20 Two isoflavone metabolites (8-hydroxydaidzein and 8-hydroxygenistein), deemed as the suicide substrate of tyrosinase, can make tyrosinase irreversibly inactivated by substrate-like interaction.14 Tea polyphenols (TPs), such as (+)-catechin (catechin), (−)-epicatechin gallate (ECG) and (−)-epigallocatechin-3-O-gallate (EGCG), which belong to flavanols, frequently present in human diet due to their strong antioxidant activity. The report has shown that they can competitively inhibit the monophenolase activity of tyrosinase in vitro.21 Particularly, EGCG can penetrate the cell membrane and suppress the activity of tyrosinase in B16 murine melanoma cell.22 Although the inhibitory activity of these TPs against tyrosinase has been reported, their inhibitory mechanisms are still elusive.

In this work, we managed to reveal the inhibitory mechanisms of TPs (including catechin, ECG and EGCG as depicting in Fig. 1) using fluorescence spectroscopy, cyclic voltammetry (CV) and oximetry together with molecular simulations. Fluorescence spectra determined the nature of interactions between TPs and tyrosinase. CV identified the inhibitory mechanism of TPs with different functional moieties. Oxidation kinetics of tyrosinase discovered the inactivation of TPs against tyrosinase by Quintox mechanism. Molecular simulations provided a deep understanding of their inhibitory mechanism against tyrosinase at the molecular level.


image file: c7ra12749a-f1.tif
Fig. 1 Molecular structures of TPs.

2. Materials and methods

2.1 Chemical reagents and apparatus

Tyrosinase (TYR, EC1.14.18.1) from mushroom, L-3,4-dihydroxyphenylalanine (L-DOPA), kojic acid (KA), catechin, ECG, and EGCG were purchased from Sigma-Aldrich. Nitric acid and dimethyl formamide (DMF) were purchased from Beijing Chemical Works. Multiwalled carbon nanotubes (MWCNTs) were purchased from the Nanotech Port Co. Ltd (Shenzhen, China). Before use, the MWCNTs were acidified by nitric acid.23

Spectrophotometer reader (BioTek Instruments Inc., USA) and High Performance Liquid Chromatograph (HPLC: Shimadzu Instruments Inc., Japan) were used to measure the inhibitory activity of the tyrosinase inhibitors. CV experiments were performed in a standard three-electrode electrochemical cell with a CHI 660E electrochemical analyzer (CH Instruments, Chenhua Co., Shanghai, China). Oxygen consumptions of tyrosinase interacting with TPs were monitored with YSI 5000 dissolved oxygen meter (YSI Inc., American).

2.2 Inhibitory activity in vitro

The inhibitory activity of the TPs, i.e., catechin, ECG, and EGCG, and a positive control, i.e., KA, against diphenolase activity of tyrosinase was measured with L-DOPA as a substrate using the reported spectroscopy method.13,24 Assays were conducted in a 96-well micro-plate and a spectrophotometer reader was used to determine the absorbance.25 Briefly, 6 μL of 1000 unit per mL tyrosinase solution was combined with potassium phosphate buffer (pH 6.5) and 1 mM KA or TPs (both dissolved in potassium phosphate buffer, pH 6.5). After pre-incubating at 25 °C for 5 min, 150 μL of 1 mM L-DOPA aqueous solution was added, followed by incubation at 25 °C for 10 min, again, before the concentration determination of DOPA quinone at 475 nm absorbance. The total volume of each reaction system was 300 μL. As an exception, the tyrosinase inhibitory activity of catechin was measured with HPLC rather than spectroscopy, as it possessed an overlap adsorption band with DOPA quinone which will interfere the analysis of inhibitory activity.26 The inhibitory activity was expressed by IC50 value, which is equal to the concentration of inhibitor at 50% inhibitory rate according to the following equation:25
 
image file: c7ra12749a-t1.tif(1)
where, A is the optical density (OD475) of potassium phosphate buffer; B, C and D are the OD475 of potassium phosphate buffer with tyrosinase, with both TPs and tyrosinase, and with TPs, respectively.

2.3 Fluorescence spectroscopy

The fluorescence spectrum of the TP-TYR complex was measured according to our previous work.27 Briefly, 10 mM TP was titrated into 3 mL potassium phosphate buffer (pH 6.5) containing 10 units tyrosinase under continuous stirring and nitrogen exists. After 5 minutes, tyrosinase was excited at 274 nm and the emission spectrum over the range of 280 nm to 400 nm was detected through a 3 nm slit. The emission spectrum of the tyrosinase solution was also directly measured as a background.

2.4 Cyclic voltammetry

The acidified MWCNT was suspended in DMF at a concentration of 200 mg mL−1, and then 5 μL of this dispersion was dropped onto the polished surface of the GC electrode and allowed to dry naturally for 12 h. The dried electrode was submerged in 3 mg mL−1 tyrosinase dissolved in 10 mM potassium phosphate buffer (pH 4.5) for 24 h.

CVs were performed in an electrolyte of 10 mM TP and 10 mM potassium phosphate buffer (pH 6.5) in the range from −0.2 V to 0.4 V at a scan rate of 0.05 V s−1.

2.5 Oximetry

The oxygen consumption in the process of tyrosinase interacting with TPs was monitored with the reported method.28 Briefly, 0.4 mL of 8.2 mM TP was added to 19.2 mL phosphate buffer (pH 6.5). After equilibrating at 30 °C for 5 min, 0.4 mL of 500 unit per ml tyrosinase was injected to TP solution under continuous stirring. The decrease of oxygen was monitored with the dissolved oxygen meter. In order to evaluate the remaining activity of tyrosinase, L-DOPA was added to the above mixture of tyrosinase and TP after 5 minutes if no oxygen decreasing was observed, followed by monitoring the secondary oxygen utilization in the oxidation reaction of L-DOPA catalyzed by the remaining tyrosinase. The concentration of L-DOPA was 194 μM.

2.6 Molecular simulations

The TP structures were generated using Molinspiration Galaxy 3D Structure Generator v. 2013.02 beta (https://www.molinspiration.com) followed by optimization using B3LYP29–31 with a 6-31G (d, p) basis set implanted in Gaussian 09. The crystal structure of tyrosinase with tropolone (PDB ID: 2Y9X)32 was taken from the Protein Data Bank (PDB).33 Ligands and crystal water molecules were removed and missing hydrogen atoms were added using Schrödinger software.34 The protonation and disassociation of chargeable residues (His, Asp, Glu, Arg, and Lys) and N- and C-terminus residues were set at pH 7.0.35,36 Partial charges in the six histidine residues (60, 84, 93, 258, 262, and 295) and two copper ions in the catalytic center were set according to the report of Choi et al.37 All residues of the protein were parameterized with the AMBER-03 force field.38 The cluster of bi-copper ions coordinated by six histidine residues were firstly optimized using B3LYP with a 6-31G (d, p) basis set, and then the parameters of bonds, angles and dihedral formed between the histidine residues and copper ions were evaluated using the Visual Force Field Derivation Toolkit (VFDT).39 A semi-flexible docking was performed by AutoDock Vina40 to identify binding poses of TPs with tyrosinase, where tyrosinase was treated as a rigid body and all rotatable bonds in the ligand were sampled. Optimal binding sites were searched in a grid box of 30 × 30 × 30 Å3 centered on the bi-copper center of tyrosinase.

The optimal binding modes were further refined with fully flexible atomic molecular dynamics (MD) simulation using NAMD (v. 2.9).41 The binding free energy of all TPs with tyrosinase was evaluated using the Molecular Mechanics/Poisson–Boltzmann Surface Area (MM/PBSA) method42 and was further decomposed into the van der Waals (ΔEvdW), electrostatic (ΔEele), nonpolar solvation (ΔGnonpolar) and polar solvation (ΔGpolar). The conformational entropy (−TΔS) was calculated by a normal-mode analysis.43 More detailed simulation information can be found in our previous papers.44,45

3. Results and discussion

3.1 Tyrosinase inhibitory activities of TPs

The IC50 value of the TPs and KA (acted as positive control) was measured with L-DOPA as a substrate and listed in Table 1 and the profile was depicted by ESI Fig. 1. These natural compounds can inhibit diphenolase activity of tyrosinase with dose-dependent manner and their inhibitory activities against tyrosinase ranked by IC50 values were: ECG > KA > catechin > EGCG. The structural difference between ECG and EGCG is less of hydroxyl group on ring B of ECG, the tyrosinase inhibitory activity of ECG was 6-fold stronger than that of EGCG. This suggested that the catechol group (with two hydroxyl group) on the TPs was crucial in enhancing its tyrosinase inhibitory activity. We speculated that the catechol group (ring B) should have a direct interaction with a coupled binuclear copper site in the catalytic center. Moreover, we found that the pyrogallol group (ring G) would be favorable to the tyrosinase inhibitory activity of TPs by comparing ECG and catechin. The important role of the pyrogallol group in tyrosinase inhibitory activity of compounds has also been reported by Kim and Uyama.18
Table 1 The measured IC50 value and Stern–Volmer quenching constant (KSV), bimolecular quenching constant (κq), binding constant (KA) and binding number (n) obtained from fluorescence spectra
Inhibitors IC50 (μM) KSV (×104 M−1) κq (×1012 M−1 s−1) KA (×105) n
Catechin 57.12 0.48 0.48 0.04 0.99
ECG 22.63 4.10 4.10 0.53 1.03
EGCG 142.40 3.13 3.13 3.07 1.23
KA 35.70 NA NA NA NA


3.2 Fluorescence quenching

Fluorescence spectroscopy was used to further detect the interaction between tyrosinase and TPs by monitoring the quenching of the intrinsic residues in tyrosinase. The fluorescence emission spectra of tyrosinase excited at 274 nm wavelength were collected with different concentration gradients of TPs and depicted in Fig. 2. It can be observed that the fluorescence intensity of tyrosinase decreased gradually with the increase of the concentration of TPs, indicating that TPs could interact with tyrosinase and quench the fluorescence of residues. Fig. 2 shows a strong 305 nm fluorescence emission peak of tyrosinase in each TP aqueous solution, and there is no significant shift in the fluorescence emission peak of tyrosinase with either the increase of TPs concentration or the addition of other TP molecules.
image file: c7ra12749a-f2.tif
Fig. 2 Emission spectra of tyrosinase excited at 274 nm wavelength in the presence of different concentrations of catechin (A), ECG (B) and EGCG (C). Inset: Stern–Volmer plots describe the tyrosinase quenching caused by binding with TPs.

Fluorescence quenching mechanism is usually divided into static or dynamic manner.46 To clarify the quenching mechanism of the interaction between tyrosinase and TPs, we analyzed the fluorescence quenching data using the Stern–Volmer equation.47

 
image file: c7ra12749a-t2.tif(2)
where F0 and F are the relative fluorescence intensities before and after the addition of inhibitor, respectively; κq is the bimolecular quenching rate constant; τ0 is the fluorescence lifetime at 10−8 s for most biological macromolecules,48 [Q] is the concentration of inhibitor and KSV is the Stern–Volmer quenching constant. The Stern–Volmer quenching constant (KSV) determined by linear regression of a curve of F0/F versus [Q] was plotted in Fig. 2 (inset). As listed in Table 1, the KSV value of ECG with less one hydroxyl group on ring B is much higher than that of EGCG. Good linear relationships in three Stern–Volmer plots of TPs are generally the indication of a single type of quenching manner. The calculated κq values by KSV/τ0 are 0.48 × 1012, 4.10 × 1012 and 3.13 × 1012 M−1 for catechin, ECG and EGCG, respectively. The κq values of tyrosinase quenching procedure caused by three TPs were higher than the limiting diffusion collision quenching constant of the biomolecules (2.0 × 1010 M−1 s−1) by two order of magnitude, suggesting that the quenching process of tyrosinase did not follow dynamic manner,49 but was a static one.48,50 For static quenching mechanism, the binding site number (n) and association constant (KA) can be obtained by the following equation:49
 
image file: c7ra12749a-t3.tif(3)

The binding sites number of each TP with tyrosinase approximately equals to 1 (ESI Fig. 2), which indicates that one tyrosinase could only be bound by one TP molecule. However, the calculated association constants (KA) are in the order of EGCG (3.07 × 105) > ECG (0.53 × 105) > catechin (0.04 × 105), being significant inverse relation with their inhibitory activities. This may suggest that the inhibitory mechanism of TPs on tyrosinase is not apparently competitive with its substrate binding into the catalytic site of tyrosinase.

3.3 Inhibitory mechanism of TPs

A few of investigations have reported51–55 that flavonoids inhibited tyrosinase activity either by chelating with the coupled bi-copper ions or by acting as substrate of tyrosinase. To deeply understand tyrosinase inhibitory mechanism of TPs, belonging to flavonoids, we traced the fluctuation of oxidized current accompanied by the addition of TPs into tyrosinase solution using the cyclic voltammetry (CV) approach, which has been extensively applied to determine the antioxidant potential of flavonoids.56–59 Fig. 3 depicts the CV curves of catechin, ECG and EGCG. The presence of significant current peaks in the CV curves indicates that three TPs should be oxidized by tyrosinase as its substrate. The catechin has always been used as a substrate of tyrosinase.51,60–62 Thus, we speculated that three TPs could interact with tyrosinase with substrate-like manner and inhibit tyrosinase activity by reducing tyrosinase. As shown in Fig. 3, two peaks were detected for ECG and EGCG, while for catechin only one peak was found. The peak at 0.26 V in the CV curve of catechin should arise from the oxidization of the catechol group (ring B with ortho-diphenols) catalyzed by tyrosinase, because meta-diphenols, i.e. ring A, have been found to have higher oxidation potentials than ortho-diphenols.56 Kilmartin et al.57 reported that the oxidization peak of meta-diphenols on ring A of flavonoids (such as catechin) was higher than 0.7 V. By comparing the CV curves of ECG with the one of catechin, it can be concluded that the first peak of ECG (0.27 V) corresponds to the oxidation of the catechol group (ring B), contained by both catechin and ECG. The difference between the CV curves of ECG and EGCG can identify that the second peak originates from the oxidation of the pyrogallol group (ring G) of both ECG and EGCG, while the first peak of EGCG (0.21 V) is caused by the oxidation of its pyrogallol group (ring B with three hydroxyl groups). The lower oxidization peak potential of the pyrogallol group compared with the catechol group has also been found by Furuno.63 It should be noted that the first peak of ring B of EGCG has a much lower oxidation potential than that of ECG, indicating it should possess much stronger reductive activity. Thus, the order of ability in reducing tyrosinase should be EGCG > catechin ≈ ECG. However, this order is inversely correlated with their order of tyrosinase inhibitory activity. One possible explanation is that the oxidation of catechin, ECG and EGCG catalyzed by tyrosinase should undergo different reaction pathways, leading to the different inhibitory mechanisms.
image file: c7ra12749a-f3.tif
Fig. 3 Cyclic voltammetry (CV) curves of 20 μM catechin (A), ECG (B) and EGCG (C) measured on glassy carbon electrode modified by MWCNT and tyrosinase at 50 mV s−1 in pH 6.5 phosphate buffer. The oxidation peaks are marked.

To further reveal tyrosinase inhibitory mechanisms of TPs, we monitored oxygen utilization against time starting from the injecting of TPs into tyrosinase solution using oximetry. The corresponding time-dependent curves of oxygen utilization are depicted in Fig. 4A. In the initial stage, the increase of oxygen utilization for TPs further corroborates that TPs could be oxidized by tyrosinase like-as substrate manner. The time-dependent curves of oxygen utilization can be better fitted with the following equation:28

 
image file: c7ra12749a-t4.tif(4)
where, Ut is the oxygen consumption at time t; E0 is the initial amount of tyrosinase; k1 and k2 are the oxidation rate and the inactivation rate of tyrosinase, respectively. As listed in Table 2, ECG and catechin can inactivate tyrosinase with inactivation rate (k2) of 0.018 and 0.016, respectively, while the inactivation rate of EGCG is ∼4 times lower than catechin. The order of tyrosinase inactivation rate by TPs is in according with their inhibitory activity, indicating that whether inactivating tyrosinase more or not is the vital factor for the inhibitory activity of TPs. Furthermore, the catecholase activities of residual tyrosinase were evaluated by the oxygen utilization of L-DOPA. Fig. 4B displays the curves of oxygen utilization against time in the oxidation reaction of L-DOPA catalyzed by the residual tyrosinase in each TP solution. The slope (K) of oxygen consumption curve is positive proportional to the amount of residual tyrosinase.64 Obviously, the remaining quantity of tyrosinase for EGCG is largest, while the residual amount is least for ECG. This indicates that the catechol group (ring B with ortho-dihydroxy groups) on ECG or catechin plays crucial roles in inhibiting tyrosinase activity. Therefore, we proposed that the catechol group (ring B) of ECG and catechin resulted in the inactivation of tyrosinase by reducing tyrosinase with Quintox mechanism through cresolase-like catalytic pathway, as reported by Land et al.28 In this inactivation pathway, the catechol is converted to an intermediate by deprotonation and reductive elimination along with the leaving of forming Cu(0) from catalytic center, leading to the irreversible inactivation of tyrosinase.28 However, the oxidation of ring B on EGCG or ring G on both ECG and EGCG (the pyrogallol group) catalyzed by tyrosinase should be competitive with the natural substrate (such as L-DOPA) through the diphenolase-like catalytic pathway (as shown in Fig. 5), because the three adjacent hydroxyl groups on the pyrogallol group having larger steric hindrance and worse orientation were unfavorable to the insert of additional oxygen on the aromatic ring.28 Tyrosinase at last can be reactivated after catalyzing the oxidation of two molecules of EGCG, accounting for the lower inhibitory activity of EGCG than ECG and catechin. The partition ratio (r), which is the molar proportion of inhibitors and inactive enzyme, is usually used to describe the inactivation capability of substrate with inactivating ability.14,65 As shown in ESI Fig. 3 (the intercept on abscissa), inactivating one molecule of tyrosinase requires 14.7 molecules of ECG or 38.7 molecules of catechin, being better consistent with their inactivation rate. The lower inactivation efficiency against tyrosinase explained that catechin has weaker inhibitory activity than ECG. Overall, three TPs inhibit the activity of tyrosinase with different inhibitory mechanisms: the catechin inactivates tyrosinase irreversibly with lower efficiency; EGCG inhibits tyrosinase activity by competing with or delaying the oxidation of substrate; ECG possesses both two inhibitory mechanisms and higher inactivation efficiency against tyrosinase.


image file: c7ra12749a-f4.tif
Fig. 4 The time-dependent curves of oxygen utilization in the oxidation of (A) TPs catalyzed by tyrosinase and (B) L-DOPA catalyzed by the remaining tyrosinase. The initial reaction velocities of tyrosinase in oxidizing L-DOPA are illustrated by the slope of oxygen utilization curves at the starting reaction time.
Table 2 The oxidation rate (k1), inactivation rate (k2) of TPs against tyrosinase and the slope (K) of the residual tyrosinase oxidizing L-DOPA substrate
Inhibitors k1 k2 K
Catechin 0.063 0.016 0.009
ECG 0.068 0.018 0.007
EGCG 0.007 0.004 0.390



image file: c7ra12749a-f5.tif
Fig. 5 Cresolase-like and diphenolase-like catalytic pathways of tyrosinase, which catalyze the monophenols (black) and diphenols (blue), respectively. Both catalytic mechanisms were reported by Solomon et al. and Land et al. and re-depicted here.

3.4 Binding behaviors of TPs with tyrosinase

To further confirm the different inhibitory mechanisms of TPs at the molecular level, we investigated the binding behaviors of TPs with tyrosinase using molecular docking, molecular dynamics simulations and MM/PBSA calculations. The binding mode of TPs was determined by molecular docking. As depicted in ESI Fig. 4, all functional groups (ring A, B and G) of three TPs can approach to the catalytic center of tyrosinase containing bi-copper ions. Based on the optimal binding conformation for each functional group, we performed 30 ns MD simulations to refine the structure of TP–TYR complex. All TPs survive the binding pose during the whole simulation time. Then, the binding free energies of each functional group with tyrosinase were evaluated using MM/PBSA method, followed by decomposing into four energy components, i.e., ΔEvdW, ΔEele, ΔGnonpolar and ΔGpolar. The entropic contributions (−TΔS) were calculated using normal mode analysis43 based on five conformations. Table 3 lists the results of free energy calculations. It can be observed that the ring B for three TPs with the lower binding free energies has higher proneness to closely contact with the catalytic center of tyrosinase than other functional groups. The order of binding free energies of catechin, ECG and EGCG is highly consistent with their associated constants (KA) determined by fluorescence quenching. Although EGCG possesses the highest binding ability with tyrosinase, the weakest ability of inactivating tyrosinase than both catechin and ECG results in its lowest tyrosinase inhibitory activity. It was obvious that ECG with the pyrogallol group (ring G) exhibited much stronger binding with tyrosinase than catechin due to the lower binding free energies. This is consistent with the measured inhibitory activity. The larger van der Waals contribution for binding free energy of ECG compared with catechin also explained why the pyrogallol group could enhance the tyrosinase inhibitory activity of TPs very well. The weaker binding ability of catechin with tyrosinase than ECG may be the main reason of its lower inactivation efficiency against tyrosinase. In addition, it can be noted that the pyrogallol group (ring G) also enhances the binding free energy or associated constants (KA) of EGCG with tyrosinase, but EGCG has weaker tyrosinase inhibitory activity than catechin due to the lack of the catechol group on EGCG, which gives rise to the inactivation of tyrosinase. Moreover, we also found that ring B of ECG and EGCG had higher possibility to closely contact with the catalytic center of tyrosinase than ring G, indicating that their tyrosinase inhibitory activities mainly originated from the contributions of ring B. Further, geometrically, the orientation of ring B of ECG and catechin was approximately orthogonal to the line connecting the coupling copper ions with ∼90° included angles (as shown in Fig. 6), indicating both of them could present in the optimal orientation required by Quintox mechanism,66 which results in the inactivation of tyrosinase, while the ring B of EGCG with lowest inhibitory activity possesses unfavorable geometric orientation for the cresolase-like pathway. This further confirms that EGCG inhibits tyrosinase activity by mainly competing with or delaying the oxidation of substrate.
Table 3 Binding free energy (kcal mol−1) and the energy components for catechin, ECG and EGCG with tyrosinase and their orientationsa
Complexes ΔEvdW ΔEele ΔGpolar ΔGnonpolar TΔS ΔGbind
TPs Binding group
a The van der Waals (ΔEvdW), electrostatic interaction (ΔEele), polar solvation free energy (ΔGpolar) and non-polar solvation energy (ΔGnonpolar) are calculated and decomposed using MM/PBSA method and the change of conformational entropy (−TΔS) is evaluated by normal mode analysis. The binding free energy (ΔGbind) is the total of five energy items.
Catechin Ring A −26.63 11.16 −15.43 −3.66 22.21 −12.35
Ring B −22.97 6.99 −14.01 −3.66 21.26 −12.39
ECG Ring A −40.63 4.76 −9.48 −4.91 26.19 −24.07
Ring B −45.58 8.56 −12.05 −5.03 25.25 −28.85
Ring G −35.53 5.63 −11.65 −5.25 29.83 −16.97
EGCG Ring A −34.07 3.99 −4.48 −4.95 25.15 −14.36
Ring B −41.52 3.06 −9.49 −5.02 23.43 −29.54
Ring G −36.11 3.60 −8.56 −5.28 34.03 −12.32



image file: c7ra12749a-f6.tif
Fig. 6 The optimal binding structures of catechin (A), ECG (B) and EGCG (C) with tyrosinase clustered out from MD simulation trajectories. The dicopper ions in the catalytic center are represented as orange spheres, tyrosinase and the skeleton of TPs are shown using cartoon and ball-and-stick model, respectively. The included angles between the ring B of TPs and the connected line of two copper ions (illustrated in D) and the distances between the oxygen atoms of ring B and copper ions were labeled.

4. Conclusion

In this work, fluorescence spectroscopy, cyclic voltammetry, oximetry and molecular simulations were employed to shed light on the tyrosinase inhibitory mechanisms of three TPs, including catechin, ECG, and EGCG. Their inhibitory activities are in the order of ECG > catechin > EGCG. The catechol group (ring B) is necessary to inhibit tyrosinase activity with high efficient. TPs have potent quenching effect on the intrinsic fluorescence of tyrosinase with static quenching manner and interact with tyrosinase with the molar proportion of 1[thin space (1/6-em)]:[thin space (1/6-em)]1. Cyclic voltammetry profiles and oxidization kinetics of tyrosinase suggest that catechin, ECG and EGCG inhibit the activity of tyrosinase, but with different inhibitory mechanisms. The high inhibitory activities of ECG and catechin were accomplished mainly by making tyrosinase inactivate through cresolase-like catalytic pathway, while EGCG is by competing with or delaying the oxidation of substrate mainly through diphenolase-like catalytic pathway. Molecular simulations further verified the speculated two types of inhibitory mechanisms.

Conflicts of interest

There are no conflicts to declare.

Acknowledgements

This work was supported by National Natural Science Foundation of China (21504092 and 21374117), and One Hundred Person Project of the Chinese Academy of Sciences. We are also grateful to Computing Center of Jilin Province for essential support.

References

  1. E. I. Solomon, U. M. Sundaram and T. E. Machonkin, Chem. Rev., 1996, 96, 2563–2605 CrossRef CAS PubMed.
  2. I. Kubo, Q. Chen and K. Nihei, Food Chem., 2003, 81, 241–247 CrossRef CAS.
  3. W. Yi, R. Cao, W. Peng, H. Wen, Q. Yan, B. Zhou, L. Ma and H. Song, Eur. J. Med. Chem., 2010, 45, 639–646 CrossRef CAS PubMed.
  4. M. Sikora and M. Swieca, Food Chem., 2018, 239, 1160–1166 CrossRef CAS PubMed.
  5. M. Friedman, J. Agric. Food Chem., 1996, 44, 631–653 CrossRef CAS.
  6. Z. Pan, H. Li, X. Yu, Q. Zuo, G. Zheng, Y. Shi, X. Liu, Y. Lin, G. Liang, Q. Wang and Q. Chen, J. Agric. Food Chem., 2011, 59, 6645–6649 CrossRef CAS PubMed.
  7. W. Chai, Y. Shi, H. Feng, L. Xu, Z. Xiang, Y. Gao and Q. Chen, J. Agric. Food Chem., 2014, 62, 6382–6389 CrossRef CAS PubMed.
  8. L. Zhou, Z. Xiong, W. Liu and L. Zou, Food Chem., 2017, 234, 445–454 CrossRef CAS PubMed.
  9. H. Satooka, P. Cerda, H. Kim, W. F. Wood and I. Kubo, Biochem. Biophys. Res. Commun., 2017, 487, 840–846 CrossRef CAS PubMed.
  10. G. P. Li, D. Zhou, L. N. Kan, Y. W. Wu, J. F. Fan and J. Ouyang, Acta Aliment., 2017, 46, 100–108 CrossRef.
  11. X. Chen, L. Ren, M. Li, J. Qian, J. Fan and B. Du, Food Chem., 2017, 214, 432–439 CrossRef CAS PubMed.
  12. B. Yang, J. Yang, Y. Zhao, H. Liu and Y. Jiang, Curr. Med. Chem., 2016, 23, 3078–3115 CrossRef CAS PubMed.
  13. S. H. Jeong, Y. B. Ryu, M. J. Curtis Long, H. W. Ryu, Y. S. Baek, J. E. Kang, W. S. Lee and K. H. Park, J. Agric. Food Chem., 2009, 57, 1195–1203 CrossRef CAS PubMed.
  14. T. Chang, J. Agric. Food Chem., 2007, 55, 2010–2015 CrossRef CAS PubMed.
  15. L. Xie, Q. Chen, H. Huang, H. Wang and R. Zhang, Biochemistry-Moscow, 2003, 68, 487–491 CrossRef CAS.
  16. Z. Zheng, K. Cheng, J. T. To, H. Li and M. Wang, Mol. Nutr. Food Res., 2008, 52, 1530–1538 CAS.
  17. Y. J. Kim and T. Yokozawa, Biol. Pharm. Bull., 2009, 32, 1155–1159 CAS.
  18. Y. J. Kim and H. Uyama, Cell. Mol. Life Sci., 2005, 62, 1707–1723 CrossRef CAS PubMed.
  19. I. Kubo, I. KinstHori, S. K. Chaudhuri, Y. Kubo, Y. Sanchez and T. Ogura, Bioorg. Med. Chem., 2000, 8, 1749–1755 CrossRef CAS PubMed.
  20. T. Chang, Int. J. Mol. Sci., 2009, 10, 2440–2475 CrossRef CAS PubMed.
  21. J. K. No, D. Y. Soung, Y. J. Kim, K. H. Shim, Y. S. Jun, S. H. Rhee, T. Yokozawa and H. Y. Chung, Life Sci., 1999, 65, PL241–PL246 CrossRef CAS PubMed.
  22. Y. Liang, S. Kang, L. Deng, L. Xiang and X. Zheng, Trop. J. Pharm. Res., 2014, 13, 1825–1831 CrossRef CAS.
  23. J. Liu, A. G. Rinzler, H. J. Dai, J. H. Hafner, R. K. Bradley, P. J. Boul, A. Lu, T. Iverson, K. Shelimov, C. B. Huffman, F. Rodriguez-Macias, Y. S. Shon, T. R. Lee, D. T. Colbert and R. E. Smalley, Science, 1998, 280, 1253–1256 CrossRef CAS PubMed.
  24. M. J. Matos, C. Varela, S. Vilar, G. Hripcsak, F. Borges, L. Santana, E. Uriarte, A. Fais, A. Di Petrillo, F. Pintus and B. Era, RSC Adv., 2015, 5, 94227–94235 RSC.
  25. C. Chen, P. Kuo, Y. Chen, J. Huang, M. Ho, R. Lin, J. Chang and H. Wang, J. Taiwan Inst. Chem. Eng., 2010, 41, 129–135 CrossRef CAS.
  26. B. Gasowska-Bajger and H. Wojtasek, J. Agric. Food Chem., 2016, 64, 5417–5427 CrossRef CAS PubMed.
  27. C. Shi, H. Tang, J. Xiao, F. Cui, K. Yang, J. Li, Q. Zhao, Q. Huang and Y. Li, J. Agric. Food Chem., 2017, 65, 656–665 CrossRef CAS PubMed.
  28. E. J. Land, C. A. Ramsden and P. A. Riley, Tohoku J. Exp. Med., 2007, 212, 341–348 CrossRef CAS PubMed.
  29. A. D. Becke, J. Chem. Phys., 1992, 96, 2155–2160 CrossRef CAS.
  30. A. D. Becke, J. Chem. Phys., 1992, 97, 9173–9177 CrossRef CAS.
  31. A. D. Becke, J. Chem. Phys., 1993, 98, 5648–5652 CrossRef CAS.
  32. W. T. Ismaya, H. J. Rozeboom, A. Weijn, J. J. Mes, F. Fusetti, H. J. Wichers and B. W. Dijkstra, Biochemistry, 2011, 50, 5477–5486 CrossRef CAS PubMed.
  33. H. M. Berman, J. Westbrook, Z. Feng, G. Gilliland, T. N. Bhat, H. Weissig, I. N. Shindyalov and P. E. Bourne, Nucleic Acids Res., 2000, 28, 235–242 CrossRef CAS PubMed.
  34. G. M. Sastry, M. Adzhigirey, T. Day, R. Annabhimoju and W. Sherman, J. Comput.-Aided Mol. Des., 2013, 27, 221–234 CrossRef PubMed.
  35. D. Liebschner, M. Dauter, A. Brzuszkiewicz and Z. Dauter, Acta Crystallogr., Sect. D: Biol. Crystallogr., 2013, 69, 1447–1462 CAS.
  36. S. J. Fisher, M. P. Blakeley, M. Cianci, S. McSweeney and J. R. Helliwell, Acta Crystallogr., Sect. D: Biol. Crystallogr., 2012, 68, 800–809 CAS.
  37. J. Choi, K. Choi, S. J. Park, S. Y. Kim and J. Jee, J. Chem. Inf. Model., 2016, 56, 354–367 CrossRef CAS PubMed.
  38. Y. Duan, C. Wu, S. Chowdhury, M. C. Lee, G. M. Xiong, W. Zhang, R. Yang, P. Cieplak, R. Luo, T. Lee, J. Caldwell, J. M. Wang and P. Kollman, J. Comput. Chem., 2003, 24, 1999–2012 CrossRef CAS PubMed.
  39. S. Zheng, Q. Tang, J. He, S. Du, S. Xu, C. Wang, Y. Xu and F. Lin, J. Chem. Inf. Model., 2016, 56, 811–818 CrossRef CAS PubMed.
  40. O. Trott and A. J. Olson, J. Comput. Chem., 2010, 31, 455–461 CAS.
  41. J. C. Phillips, R. Braun, W. Wang, J. Gumbart, E. Tajkhorshid, E. Villa, C. Chipot, R. D. Skeel, L. Kalé and K. Schulten, J. Comput. Chem., 2005, 26, 1781–1802 CrossRef CAS PubMed.
  42. J. Srinivasan, T. E. Cheatham, P. Cieplak, P. A. Kollman and D. A. Case, J. Am. Chem. Soc., 1998, 120, 9401–9409 CrossRef CAS.
  43. B. R. Miller III, T. D. McGee Jr, J. M. Swails, N. Homeyer, H. Gohlke and A. E. Roitberg, J. Chem. Theory Comput., 2012, 8, 3314–3321 CrossRef PubMed.
  44. F. Cui, K. Yang and Y. Li, PLoS One, 2015, 10, e0125848 Search PubMed.
  45. F. Cui, L. Liu, H. Tang, K. Yang and Y. Li, Chem. Biol. Drug Des., 2017, 89, 482–491 CAS.
  46. J. Xiao, X. Chen, X. Jiang, M. Hilczer and M. Tachiya, J. Fluoresc., 2008, 18, 671–678 CrossRef CAS PubMed.
  47. D. Kim, J. Park, J. Kim, C. Han, J. Yoon, N. Kim, J. Seo and C. Lee, J. Agric. Food Chem., 2006, 54, 935–941 CrossRef CAS PubMed.
  48. M. R. Eftink, Fluorescence Quenching Reactions, in Biophysical and Biochemical Aspects of Fluorescence Spectroscopy, ed. T. G. Dewey, Springer, Boston, MA, 1991, pp. 1–41 Search PubMed.
  49. M. Jiang, M. X. Xie, D. Zheng, Y. Liu, X. Y. Li and X. Chen, J. Mol. Struct., 2004, 692, 71–80 CrossRef.
  50. W. R. Ware, J. Phys. Chem., 1962, 66, 455–458 CrossRef CAS.
  51. Y. J. Kim, J. E. Chung, M. Kurisawa, H. Uyama and S. Kobayashi, Biomacromolecules, 2004, 5, 474–479 CrossRef CAS PubMed.
  52. P. Maisuthisakul and M. H. Gordon, Food Chem., 2009, 117, 332–341 CrossRef CAS.
  53. Y. Xue, T. Miyakawa, Y. Hayashi, K. Okamoto, F. Hu, N. Mitani, K. Furihata, Y. Sawano and M. Tanokura, J. Agric. Food Chem., 2011, 59, 6011–6017 CrossRef CAS PubMed.
  54. M. R. Loizzo, R. Tundis and F. Menichini, Compr. Rev. Food Sci. Food Saf., 2012, 11, 378–398 CrossRef CAS.
  55. Y. Si, Z. Wang, D. Park, H. Chung, S. Wang, L. Yan, J. Yang, G. Qian, S. Yin and Y. Park, Int. J. Biol. Macromol., 2012, 50, 257–262 CrossRef CAS PubMed.
  56. P. A. Kilmartin, H. L. Zou and A. L. Waterhouse, J. Agric. Food Chem., 2001, 49, 1957–1965 CrossRef CAS PubMed.
  57. P. A. Kilmartin, H. L. Zou and A. L. Waterhouse, Am. J. Enol. Vitic., 2002, 53, 294–302 CAS.
  58. S. Gunckel, P. Santander, G. Cordano, J. Ferreira, S. Munoz, L. J. Nunez-Vergara and J. A. Squella, Chem.-Biol. Interact., 1998, 114, 45–59 CrossRef CAS PubMed.
  59. E. Diego, L. Agui, A. Gonzalez-Cortes, P. Yanez-Sedeno, J. M. Pingarron and J. M. Kauffmann, Electroanalysis, 1998, 10, 33–38 CrossRef CAS.
  60. M. Tse, S. Kermasha and A. Ismail, J. Mol. Catal. B: Enzym., 1997, 2, 199–213 CrossRef CAS.
  61. S. Kermasha, H. H. Bao and B. Bisakowski, J. Mol. Catal. B: Enzym., 2001, 11, 929–938 CrossRef CAS.
  62. M. Y. Moridani, H. Scobie, P. Salehi and P. J. O'Brien, Chem. Res. Toxicol., 2001, 14, 841–848 CrossRef CAS PubMed.
  63. K. Furuno, T. Akasako and N. Sugihara, Biol. Pharm. Bull., 2002, 25, 19–23 CAS.
  64. M. R. L. Stratford, C. A. Ramsden and P. A. Riley, Bioorg. Med. Chem., 2013, 21, 1166–1173 CrossRef CAS PubMed.
  65. Z. Wang, J. Theor. Biol., 1990, 147, 497–508 CrossRef CAS PubMed.
  66. M. Goldfeder, M. Kanteev, S. Isaschar-Ovdat, N. Adir and A. Fishman, Nat. Commun., 2014, 5, 4505 CAS.

Footnote

Electronic supplementary information (ESI) available. See DOI: 10.1039/c7ra12749a

This journal is © The Royal Society of Chemistry 2018