Next Article in Journal
Polyketide Derivatives from the Mangrove-Derived Fungus Penicillium sp. HDN15-312
Next Article in Special Issue
Preparation and Vasodilation Mechanism of Angiotensin-I-Converting Enzyme Inhibitory Peptide from Ulva prolifera Protein
Previous Article in Journal
Identification of Penexanthone A as a Novel Chemosensitizer to Induce Ferroptosis by Targeting Nrf2 in Human Colorectal Cancer Cells
Previous Article in Special Issue
Systematical Investigation on Anti-Fatigue Function and Underlying Mechanism of High Fischer Ratio Oligopeptides from Antarctic Krill on Exercise-Induced Fatigue in Mice
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

In Silico Identification and Molecular Mechanism of Novel Tyrosinase Inhibitory Peptides Derived from Nacre of Pinctada martensii

1
College of Food Science and Technology, National Research and Development Branch Center for Shellfish Processing (Zhanjiang), Guangdong Provincial Key Laboratory of Aquatic Products Processing and Safety, Guangdong Provincial Engineering Technology Research Center of Seafood, Guangdong Province Engineering Laboratory for Marine Biological Products, Guangdong Ocean University, Zhanjiang 524088, China
2
Shenzhen Institute of Guangdong Ocean University, Shenzhen 518108, China
3
Collaborative Innovation Center of Seafood Deep Processing, Dalian Polytechnic University, Dalian 116034, China
4
Guangdong Shaohe Pearl Co., Ltd., Shantou 515041, China
*
Author to whom correspondence should be addressed.
Mar. Drugs 2024, 22(8), 359; https://doi.org/10.3390/md22080359
Submission received: 15 July 2024 / Revised: 5 August 2024 / Accepted: 5 August 2024 / Published: 7 August 2024
(This article belongs to the Special Issue Marine Bioactive Peptides—Structure, Function, and Application 2.0)

Abstract

:
Pearl and nacre powders have been valuable traditional Chinese medicines with whitening properties for thousands of years. We utilized a high-temperature and high-pressure method along with compound enzyme digestion to prepare the enzymatic hydrolysates of nacre powder of Pinctada martensii (NP-PMH). The peptides were identified using LC–MS/MS and screened through molecular docking and molecular dynamics simulations. The interactions between peptides and tyrosinase were elucidated through enzyme kinetics, circular dichroism spectropolarimetry, and isothermal titration calorimetry. Additionally, their inhibitory effects on B16F10 cells were explored. The results showed that a tyrosinase-inhibitory peptide (Ala-His-Tyr-Tyr-Asp, AHYYD) was identified, which inhibited tyrosinase with an IC50 value of 2.012 ± 0.088 mM. The results of the in vitro interactions showed that AHYYD exhibited a mixed-type inhibition of tyrosinase and also led to a more compact enzyme structure. The binding reactions of AHYYD with tyrosinase were spontaneous, leading to the formation of a new set of binding sites on the tyrosinase. The B16F10 cell-whitening assay revealed that AHYYD could reduce the melanin content of the cells by directly inhibiting the activity of intracellular tyrosinase. Additionally, it indirectly affects melanin production by acting as an antioxidant. These results suggest that AHYYD could be widely used as a tyrosinase inhibitor in whitening foods and pharmaceuticals.

Graphical Abstract

1. Introduction

Pinctada martensii is found in mainland China and the South China Sea. Its primary use is in pearl production, making it a significant marine cultural shellfish resource in China [1]. Shells and pearls are similar in structure and material composition. Both are biominerals formed under the regulation of an organic matrix secreted by the outer coat membrane tissue [2]. The nacre is secreted by the outer mantle tissue and consists of a lustrous material that gradually envelops the foreign matter to form a pearl over time [3]. The Chinese Pharmacopoeia 2020 stipulates that the source of medicinal pearls is from the pearl shellfish family, which includes pearl shellfish and other bivalves stimulated by the formation of pearls, and it has been recorded to have the effect of calming the mind and nerves, detoxifying the muscles, moisturizing the skin, and removing blemishes [4]. Pearls and nacre are rich sources of calcium, as well as shell proteins, small amounts of trace metals, hydrophobic amino acids, and aromatic amino acids. Calcium constitutes more than 90% of the weight of the pearl and nacre [5]. Studies have indicated that peptides inhibiting tyrosinase typically consist of these amino acids [6]. Therefore, nacre peptides are considered to have the potential of being tyrosinase-inhibiting peptides. Studies have shown that the main components of nacre powder are extremely similar to those of pearl powder and exhibit a variety of biological activities. Sedative-hypnotic [5], anti-photoaging [6], antioxidant [7], trauma repair [8], are some of the functional activities of pearl powder as supported by modern medicine. However, the specific whitening ingredients and mechanisms remain unclear.
Excessive melanin deposition can lead to a range of skin problems. Melanin is a phenolic polymer widely distributed in plants and animals [9]. It is the main factor in determining the color of the skin and hair and can inhibit ultraviolet rays and free radicals that stimulate damage to the skin cells [10]. However, when melanin is over-synthesized, excessive deposition of melanin in the surface layer of the skin may result in a series of skin diseases such as dark spots, brown spots, and further promote skin aging; ultimately, this excessive melanin production could lead to melanoma and other skin cancers. Tyrosinase is a copper-containing metal oxidase and is the key enzyme that regulates melanin production. Tyrosinase oxidizes tyrosine to produce dopamine quinone, which reacts with amino acids or proteins, ultimately resulting in melanin production [11]. Therefore, tyrosinase inhibition and antioxidants are particularly crucial. There are two main ways to inhibit tyrosinase: one involves direct inhibition of tyrosinase, while the other involves indirect regulation of the antioxidant system to inhibit tyrosinase production [12,13]. Currently proven effective whitening ingredients, such as Kojic acid and its derivatives, Resveratrol and Arbutin, have been found to have a stronger inhibitory effect on tyrosinase [14,15,16]. However, they are also prone to side effects such as skin redness, itching, ulcers, and even toxicity to the body [17,18]. The current whitening products available on the market, while commonly used, still carry a higher risk of causing allergic reactions or other adverse effects on the user’s skin [19,20]. Therefore, safe and stable whitening preparations are receiving increasing attention. A large number of by-products, such as shells, are produced during the pearl harvest from the Pinctada martensii [21]. The main ways of utilization include decoration, building paint, feed, and other purposes. Shell meat, as well as other tissues, has been the main focus of existing studies, while the active peptides in the inner nacre of the shell have been less studied. Previous studies have shown that peptides extracted from the shell flesh tissues of pearl oysters and their related tissues have demonstrated a variety of biological activities, such as antioxidant [22], anti-photoaging [23], and blood pressure-lowering effects [24]. Zhou et al. simulated gastrointestinal hydrolysis to prepare peptides from the meat of Pinctada martensii and verified the antioxidant activity of the peptides using HepG2 cells [22]. A study successfully predicted potential tyrosinase-inhibitory peptides from abalone proteins using an anti-tyrosinase prediction tool [25]. However, there has been limited research on potential whitening peptides within the shell of pearl oysters.
In the present study, peptides derived from the nacre powder of Pinctada martensii after undergoing high temperature and pressure treatment. Neutral proteases have been used in the enzymatic digestion of shellfish. To break down the insoluble proteins in the pearl layer, we conducted a combination of neutral proteases and pineapple protease for enzymatic digestion, and then characterized them using LC–MS/MS. The potential tyrosinase-inhibitory peptides were then screened through molecular docking. Subsequently, they were synthesized in solid phase and examined for their tyrosinase-inhibitory mechanism and antioxidant capacity at both in vitro and cellular levels. These findings could contribute to the development of new cosmetic ingredients and natural food-based functional foods.

2. Results and Discussion

2.1. Inhibition of Tyrosinase and Antioxidant Activity of the Enzymatic Hydrolysis Product

A preliminary investigation of tyrosinase activity and antioxidant activity was conducted on the enzymatic hydrolysis product obtained from nacre powder after high-temperature and high-pressure treatment and complex enzymatic digestion. The results, as presented in Table 1, indicate that the nacre peptide digest exhibits a certain inhibitory effect on tyrosinase and possesses antioxidant properties. These findings suggest that it holds promise for the development of tyrosinase-inhibiting peptides.

2.2. Identification and Molecular Docking

The mass spectrometry data were analyzed using the software PEAKS Studio 8.5. The identification results were compared using the UniProt database. For efficient screening, peptides with amino acid residues < 10 and score > 20 were selected to docking with tyrosinase (2Y9X) [26].
Molecular docking is a commonly used method to simulate the way molecules interact with each other and predict their binding modes and affinities using a computerized platform [27]. The method is widely used to probe the binding ability of small molecules interacting with large molecules [28]. A lower intermolecular binding energy indicates better peptide–protein affinity and theoretically more tyrosinase-inhibitory activity of the peptides.
As shown in Table 2, molecular docking of 32 peptides was performed using Autodock vina, and the receptor was tyrosinase [29]. In terms of amino acid composition, hydrophobic amino acids, aromatic ring-containing amino acids, and polar amino acids play a crucial role in protein–ligand interactions [30]. Additionally, the amino acid composition also influences the water solubility of peptides. Taking into account the peptide’s GRAVY score < −0.5, amino acid composition, and binding energy, three peptides—AHYYD, KPIWT, and TFSGNYP—were finally selected for the subsequent study. The predicted binding energies of AHYYD, KPIWT, and TFSGNYP were −8.0 kcal/mol, −7.4 kcal/mol, and −7.3 kcal/mol, respectively. Docking binding was slightly lower than the results for peptides extracted from Atrina pectinata Mantle, where a lower binding energy indicates a more favorable binding effect [31]. To visually demonstrate how the peptides bind to tyrosinase, a visual analysis was conducted using the Discovery Studio 2019 client tool. As shown in Table 3 and Figure 1, Docking results showed that the peptide formed interactions with residues such as HIS263, ASN260, and VAL283 near the active site pocket of tyrosinase [32]. AHYYD formed four hydrogen bonds with SER282, CYS83, GLU322, and ASN260 of tyrosinase. The bond lengths ranged from 2.10 Å to 3.18 Å. KPIWT formed two hydrogen bonds with MET280 and ASN81 of tyrosinase, with bond lengths of 2.62 Å and 3.36 Å. TFSGNYP formed six hydrogen bonds with ASN81, SER282, GLY281, VAL283, ALA246, and GLY245 of tyrosinase, and the bond lengths ranged from 2.13 to 5.26 Å. The shorter the hydrogen bond, the greater the bond energy, indicating that the ligand and receptor form a more stable bond [26]. In the hydrophobic interactions with tyrosinase, all three peptides form hydrophobic interactions with PRO284 and HIS263. KPIWT and AHYYD form hydrophobic interactions with VAL283. KPIWT and TFSGNYP form hydrophobic interactions with ALA286. On the basis of the amino acid residues where the three peptides form hydrophobic interactions, most of the amino acids they interact with are similar to each other. Due to the similarity of the hydrophobic interactions, we infer that the difference in binding energies is primarily due to variations in hydrogen bonding. Hydrogen bonding is identified as the main driving force.

2.3. Molecular Dynamics Simulation

The stability and binding affinity between tyrosinase and the three small molecules were further investigated using the GROMACS program, with the three complex docked systems as the initial structure. The 100 ns molecular dynamics (MD) simulation method was used to obtain stable complex systems for further comparison of the three small molecules at the molecular level. The trajectory data obtained from the 100 ns were analyzed to derive the kinetic properties of the three complexes. Root mean square deviation (RMSD), radius of gyration (Rg), solvent accessible surface area (SASA), root mean square fluctuation (RMSF), and Hbond number were utilized to assess each system in the MD study, and the results are presented. As shown in Figure 2, a graph of each parameter of the dynamics is plotted.
Root mean square deviation (RMSD) is a metric used to assess structural changes in proteins [33]. The results show that the RMSD between proteins and small molecules is large in the first 30 ns of the simulation. Tyrosinase with AHYYD stabilizes the fastest, the RMSD value remains within a small range, indicating that the binding of proteins to small molecules is relatively stable. Root mean square fluctuation (RMSF) is a metric used to assess the dynamics of proteins [34]. The results show that the residues in the key regions increase in flexibility, and the RMSF values are larger in the binding portion in the 300–400 region. Tyrosinase with AHYYD and tyrosinase with TFSGNYP were more variable. The RMSF of tyrosinase with AHYYD reaches a maximum value of 0.5 nm, while it is smaller in the unbound portion. The RMSF values are all below 1 Å. The minimal fluctuation indicates that these atoms create stable complexes with tyrosinase because of robust intermolecular interactions, restricting their movement in molecular dynamics simulations [35]. Gyrate is an indicator of the overall compactness of a protein and characterizes the distribution of system atoms along a specific axis [28]. The results indicate that after the binding of the protein with the three small molecules, the gyrate value of tyrosinase with AHYYD tends to stabilize after 20 ns, while the gyrate value of tyrosinase with TFSGNYP tends to decrease after 60 ns. This suggests that the binding of small molecules leads to a more compact protein structure. KPIWT, on the other hand, kept fluctuating within a small range repeatedly. SASA is a metric used to assess the surface area of proteins. The results indicate that the amplitude of the tyrosinase with AHYYD curve fluctuates slightly between 170 and 180 and then decreases slightly. This suggests that protein–protein interactions have little effect on the stability of protein molecules. The amplitudes of the tyrosinase with TFSGNYP and tyrosinase with KPIWT curves are similar, but they exhibit opposite trends. The SASA value of tyrosinase with TFSGNYP shows a decreasing trend after 60 ns, indicating that the binding of small molecules causes the protein molecule to become more compact, aligning with the results of gyrate analysis. A hydrogen bond is a metric used to assess hydrogen bonding between proteins and small molecules. The results show that many hydrogen bonds are formed between proteins and small molecules during the simulation. These bonds are dominated by interactions between key residues in tyrosinase and peptides in the molecular docking results.
Analyzing the RMSD, RMSF, Gyrate, SASA, and Hbond number results of the three peptide–enzyme complexes mentioned above, it can be concluded that AHYYD and TFSGNYP interact with tyrosinase more effectively than KPIWT, protein–small molecule interaction appears stable, and the binding of small molecules leads to a more compact protein structure.

2.4. Tyrosinase-Inhibitory Activity and Antioxidant Capacity of Synthetic Peptides

In terms of amino acid composition, AHYYD contains two repetitive benzene-ring-containing amino acid residues, Tyr, and a hydrophobic amino acid, Ala. It has been shown that the repetitive amino acid sequence enhances tyrosinase inhibition [36]. TFSGNYP contains the polar amino acid Ser, a benzene-ring-containing amino acid residue, Phe, and a hydrophobic amino acid, Pro. KPIWT has the benzene-ring-containing amino acid residue Trp and the hydrophobic amino acid Pro. It has been reported that hydrophobic amino acids, benzene-ring-containing amino acid residues, and polar amino acids contribute to tyrosinase inhibition. As shown in Table 4, the synthetic peptides AHYYD, KPIWT, and TFSGNYP derived from Pinctada martensii exhibited varying degrees of inhibitory activity against tyrosinase, with Kojic acid serving as a positive control. From the molecular docking visualization results, it can be seen that although the hydrophobic interactions between the three peptides and the enzyme are relatively similar, their distinct hydrogen-bonding interactions may affect the inhibition of tyrosinase activity in vitro. Among these, AHYYD demonstrated a significant inhibitory effect, its IC50 value is 2.012 ± 0.088 mM, consistent with the findings of molecular docking screening; this result is similar to the peptides activity identified in the mantle of Atrina pectinata [31].
As depicted in Figure 3, there was no significant difference in DPPH clearance between AHYYD and TFSGNYP at any of the five concentrations (p > 0.05). However, KPIWT showed a significant difference from both AHYYD and TFSGNYP at concentrations of 0.1–2 mg/mL; AHYYD and TFSGNYP exhibited a stronger DPPH clearance effect than KPIWT (p < 0.05) at the corresponding concentrations. There was no significant difference in the clearance of DPPH by the three peptides at a concentration of 4 mg/mL. The IC50 values for ABTS clearance by AHYYD, KPIWT, and TFSGNYP were 1.416 mg/mL, 0.7955 mg/mL, and 0.2785 mg/mL, respectively, with a significant difference between the three peptides (p < 0.05). A smaller IC50 value indicates a better scavenging effect. It can be observed that TFSGNYP exhibits the strongest scavenging ability for ABTS radicals, followed by KPIWT and AHYYD. Vitamin C (Vc) was utilized as a positive control.

2.5. Inhibition Dynamics Analysis

Based on the observation that the inhibition of tyrosinase by the synthetic peptide AHYYD was reversible, the type of inhibition was further elucidated by constructing a Lineweaver–Burk double inverse plot of the enzymatic reaction rate against the concentration of the substrate (L-tyrosine). In this plot, the intercept with the Y-axis represents 1/Vmax, while the intercept with the X-axis represents −1/Km. As depicted in Figure 4, when the concentration of tyrosinase, the concentration of synthetic peptide (0, 0.3, 0.6, 1.2 mM, respectively), and the concentration of substrate L-tyrosine remained constant, the reciprocal of the enzymatic reaction rate exhibited a linear correlation with the reciprocal of the substrate concentration. The slope of the line increased with the rise in synthetic peptide concentration, indicating a mixed type of inhibition [37]. By using the quadratic plotting method, the slopes of the four straight lines in the Lineweaver–Burk double inverse plot were plotted against the synthetic peptide concentration [38]. The inhibition constants of AHYYD on the free enzyme were calculated as Ki = 0.601 mmol/L, and as Kis = 38.375 mmol/L on the enzyme–substrate complex. The smaller the inhibition constants, the higher the binding affinity of the inhibitor to the enzyme, Ki << Kis indicates that AHYYD clearly favors the binding site with the free enzyme.

2.6. CD Spectra Analysis

CD spectroscopy is a precise method used to study the secondary structure of biomolecules, including proteins. It allows for the rapid analysis of protein conformational changes caused by the addition of ligands [39]. The addition of AHYYD caused changes in the secondary structure of tyrosinase (Figure 5). In the spectrum of free tyrosinase, there are two negative bands near 208 nm and 222 nm, which are characteristic peaks of the α-helical structure in the secondary structure [40]. The content of α-helix, β-folding, and β-turning increased from 32.53%, 16.27%, and 18.20% to 36.8%, 23.77%, and 21.93% and the random coil content decreased from 30.83% to 18.33%, respectively. The rise in the α-helix content suggests that the interaction between the peptide and tyrosinase caused the tyrosinase structure to become more compact. This finding contrasts with the interaction between tretinoin and tyrosinase, suggesting that it leads to a change in the tyrosinase conformation to some extent [37].

2.7. Isothermal Titration Calorimetry Analysis

ITC can measure the binding strength of protein–protein interactions by quantifying the thermodynamic changes in complex reactions [41]. The titration of tyrosinase by AHYYD exhibited significant heat changes, indicating binding (Figure 6). The data were analyzed using the independent model, and the thermodynamic binding parameters were calculated. The results showed that the enthalpy (∆H) was −7.282 kJ/mol, the binding constant (Ka) was 2.495 × 104 M−1, the entropy change (∆S) was 60.7 J/mol·K, and the standard free energy change (∆G) was −1.883 × 104 kJ/mol. Negative values of ∆G indicate that the bimolecular reaction is spontaneous. The change in enthalpy indicates that the process is essentially exothermic, which is consistent with the findings of Yu et al [34]. The titration curves are more consistent with the “one set of binding sites” model [42]. The mechanism of inhibition categorizes AHYYD as competitive inhibition, indicating that AHYYD binds to the active site of tyrosinase. This finding further confirms the proposed mechanism of inhibition based on inhibition kinetics.

2.8. Study on the Whitening Effect of Peptides Based on B16F10 Cells

2.8.1. Cell Viability

As depicted in Figure 7A, there was no significant difference in the viability of B16F10 cells within the concentration range of 15–1200 μM concentration of the synthetic peptide AHYYD compared to the control group (p > 0.05). This suggests that the synthetic peptide AHYYD did not exhibit any cytotoxic effects. Consequently, the concentrations of 15, 37.5, and 75 μM were chosen for the subsequent experiments.

2.8.2. Effect of AHYYD on the Tyrosinase Activity and Melanin Content of B16F10 Cells

As shown in Figure 7B,C, compared with the control group, AHYYD significantly inhibited the tyrosinase activity of B16F10 cells in the concentration range of 15–75 μM (p < 0.05). Additionally, the intracellular melanin content was significantly reduced by AHYYD in the concentration range of 15–75 μM (p < 0.05). The trend in tyrosinase activity and melanin content mirrored that of the peptides extracted from pearl shell meat [26]. At the maximum treatment concentration, the intracellular tyrosinase activity and melanin content were reduced by 9.08% and 13.48%, respectively. Kojic acid caused a reduction of 19.51% in tyrosinase activity and 18.54% in melanin content.

2.8.3. Effect of AHYYD on the Antioxidant Enzyme Activity and ROS Content of B16F10 Cells

SOD, CAT, and GSH-Px are antioxidant enzymes that collaborate to regulate antioxidant processes and safeguard cells from oxidative stress damage [43]. As can be seen from Figure 8, the activities of GSH-Px and SOD enzymes exhibited a concentration-dependent relationship with increasing concentrations of AHYYD. The activities of CAT, GSH-Px, and SOD increased by 73.05%, 139.08%, and 41.25%, respectively, under the maximum treatment concentration, and by 24.55%, 52.85%, and 68.92%, respectively. At the maximum treatment concentration, the CAT and GSH-Px activities were superior (p < 0.05) compared to Kojic acid. However, both of them were still far away from Kojic acid in terms of SOD activity (p < 0.05). Therefore, AHYYD exhibited a protective effect on the antioxidant system by enhancing the activity of antioxidant enzymes in B16F10 cells, thereby preventing melanin production. Reactive oxygen species (ROS) serve as a crucial indicator of antioxidant connectivity in the organism. As shown in Figure 8G, at a concentration of 75 μM, the intracellular ROS levels were reduced to 52.61% and 65.44% of the untreated group after AHYYD and tretinoin treatment, respectively.
GSH, GSSG, and MDA are crucial metabolic regulators in the cell [44]. GSH helps maintain normal immune system function and has antioxidant effects. It exists in two forms: reduced glutathione (GSH) and oxidized glutathione (GSSG). MDA content is a significant parameter that reflects the body’s antioxidant potential, indicating the rate and intensity of lipid peroxidation and indirectly reflecting the extent of tissue peroxidation damage. Under the maximum treatment concentration, the content of GSH increased by 44.07%, while the content of GSSG and MDA decreased by 52.86% and 56.56%, respectively. These changes were not significantly different from those observed in the positive control group treated with Kojic acid (p > 0.05). This suggests that AHYYD demonstrates a superior ability in regulating GSH, GSSG, and MDA levels. Along with the antioxidant enzymes mentioned above, it effectively protects B16F10 cells from oxidative damage.

3. Materials and Methods

3.1. Materials

The nacre powder of the Pinctada martensii was purchased from Guangdong Shaohe Pearl Co., Ltd. (Shantou, China). The tyrosinase (from Mushroom, 500 U/mg) and PBS buffer were purchased from Shanghai Yuanye Biotechnology Co., Ltd. (Shanghai, China). The L-tyrosine was obtained from Shanghai Macklin Biochemical Technology Co., Ltd. (Shanghai, China). The fetal bovine serum (FBS) was purchased from Serana Europe GmbH (Dorfstrasse 17A, Pessin, Brandenburg, Germany), and the penicillin–streptomycin mixture was purchased from Beijing Soleberg Biotechnology Co., Ltd. (Beijing, China). The DMEM/F−12 medium and trypsin-EDTA digest were purchased from Thermo Fisher Scientific (Shanghai, China).

3.2. Preparation of Enzymatic Digests of Nacre Peptides

A total of 50 g of nacre powder from Pinctada martensii was dissolved in 250 mL (v:w = 5:1) of water and heated at high temperature (121 °C) and pressure for 20 min. Then, 0.1250 g of pineapple protease (0.25% by mass of nacre powder) and 0.2500 g of neutral protease (0.5% by mass of nacre powder) were added to the pre-treated nacre powder solution, and the enzyme digestion was conducted at a temperature of 50 °C for 4 h. Finally, the mixture was heated at 100 °C for 20 min. The supernatant was centrifuged at 8000 rpm for 20 min at 4 °C to obtain the pearl peptide solution after vacuum concentration.

3.3. Tyrosinase-Inhibitory Activity and Antioxidant Activity Assay of Enzyme Digests

We followed the method of Yu [37], with a slight modification. L-tyrosine solution (0.5 mg/mL), sample solution, and PBS buffer were sequentially added into the 96-well plate, thoroughly mixed, and then incubated at 37 °C for 10 min, then 20 μL of tyrosinase solution (500 U/mL) was added into each well, followed by immediate placement into the enzyme marker after the reaction at 37 °C for (10 min ± 5 s). The results were analyzed at 475 nm. DPPH and ABTS free-radical-scavenging capacity assays were conducted using kits from Grace Biotechnology Co., Ltd. (Jiangsu, China).
T y r o s i n a s e   i n h i b i t i o n   r a t e % = 1 A 4 A 3 A 2 A 1 × 100 %
where A4 is the absorption value of the reaction well, A3 is the absorption value of the Sample Background well, A2 is the absorption value of the Solvent reaction well, and A1 is the absorption value of the Solvent base well.

3.4. Protein Sequence Identification

The amino acid sequence was identified by using the Bio-Tech Pack Technology Co. (Beijing, China). The nacre peptides were reduced and alkylated separately as samples, an Easy-nLC 1200 system coupled with a Q Exactive™ Hybrid Quadrupole-Orbitrap™ Mass Spectrometer (Thermo Fisher Scientific, Waltham, MA, USA) with an ESI nanospray source. Then, they were analyzed by liquid chromatography–mass spectrometry (LC–MS/MS) to generate a raw file of the mass spectrometry results. The total ion chromatograms were obtained using Xcalibur 4.0 software. The mass spectrometry data were analyzed using the software PEAKS Studio 8.5, and the identification results were compared with the UniProt database.

3.5. Molecular Docking Studies

Studying the peptides and tyrosinase using molecular docking methods, we referred to the method of Wang [31]. The ligands and proteins required for molecular docking were prepared using AutoDock Vina 1.1.2 software. Tyrosinase (PDB: 2Y9X) was utilized as the receptor, and a peptide was employed as the ligand. The peptides were designed using ChemDraw and energy-minimized using Chem3D. The crystal structure of the target protein was obtained from the PDB database (https://www.rcsb.org/structure/2Y9X, accessed on 25 November 2023). Autodock Vina 1.1.2 software was used to perform molecular docking with default parameters to calculate the binding energy of the peptide sequence to the enzyme. The binding-free energy (kcal/mol) of the target structure represents the binding ability of the two molecules, and the lower the binding-free energy, the more stable the binding between the ligand and the receptor. Pymol (https://pymol.org/2/, accessed on 3 December 2023) was used for visual analysis, while Discovery Studio 2019 was utilized for further exploration. Discovery Studio 2019 (https://www.3ds.com/products/biovia/discovery-studio/visualization, accessed on 3 December 2023) was used for visual analysis of the 2D images. The tyrosinase docking centers were X = −10.09, Y = −28.03, and Z = −43.14.

3.6. Molecular Dynamics Simulations

GROMACS 2022.3 software was used for the molecular dynamics simulation [45,46]. For small molecule preprocessing, AmberTools22 was employed to apply the GAFF force field to small molecules, while Gaussian 16W was used for hydrogenating small molecules and calculating the RESP potential. Potential data will be added to the topology file of molecular dynamics system. Following the completion of the simulation, the software’s built-in tool was utilized to analyze the trajectory. The root mean square variance (RMSD), root mean square fluctuation (RMSF), solvent accessible surface area (SASA), radius of gyration (Rg), and hydrogen bond were calculated for each amino acid trajectory.

3.7. Peptides Synthesis

The peptide was obtained through solid-phase synthesis and synthesized by Aminolink Biotechnology Co., Ltd. (Shanghai, China). Three peptides with purity above 98% (w/w) were detected by HPLC: Ala-His-Tyr-Tyr-Asp (AHYYD), Thr-Phe-Ser-Gly-Asn-Tyr-Pro (TFSGNYP), and Lys-Pro-Ile-Trp-Thr (KPIWT). Mobile phase A, 0.1% trifluoroacetic in 100% water; mobile phase B, 0.1% trifluoroacetic in 100% acetonitrile; flow rate, 1 mL/min; column, Inertsil ODS-SP (4.6 × 250 mm × 5 µm), with a detection wavelength of 220 nm. Eventually, the purified peptides were identified by ESI–MS spectroscopy.

3.8. Determination of Inhibition Type

Keeping the L-tyrosine concentration constant at 0.5 mg/mL, the concentrations of the synthetic peptide and tyrosinase were varied. The change in the OD475 absorbance value was continuously monitored at 37 °C after mixing. The enzyme concentration is plotted on the horizontal axis, while the rate of absorbance change is plotted on the vertical axis. Under the conditions of varying the peptide concentration and substrate concentration while keeping the tyrosinase solution constant, the reciprocal of the substrate concentration was plotted on the horizontal axis, and the reciprocal of the rate of absorbance change was plotted on the vertical ax.

3.9. CD Measurements

We used circular dichroism to study peptide–enzyme interactions, referring to the method of Yu [47]. The concentration of AHYYD was 0.6 mM, and the concentration of tyrosinase was 0.4 mg/mL. The reaction system comprised 100 μL of tyrosinase solution mixed with 20 μL of AHYYD solution, which reacted for over 30 min at 37 °C. All data were collected three times. The changes in the secondary structure were analyzed using CDNN 2.1 software.

3.10. Isothermal Titration Calorimetry Analysis

Peptide (1 mM) was used to titrate tyrosinase (0.001 mM) to study thermodynamic changes. The initial volume of the peptide solution for titration was 50 μL (inhaled into the syringe); the initial volume of the tyrosinase solution for titration was 300 μL (added to the cuvette). The reaction conditions included adding 25 drops of 2 μL each with a 120 s interval between drops at 37 °C and a stirring speed of 350 r/min; the experimental group involved titrating the peptide with tyrosinase, while the control group involved titrating the peptide with a solvent buffer. We calculated the binding constant (Ka), enthalpy change (∆H), and entropy change (∆S) during the Gibbs free energy (∆G) changes [39].
G = H T · S

3.11. Cytotoxicity of AHYYD to B16F10 Cells

The B16F10 cells were obtained from FengHui Biotechnology Co., Ltd. (Hunan, China). B16F10 cell cultures from generations 10–25 were grown in medium (DMEM supplemented with 10% fetal bovine serum, 1% antibiotic–antifungal solution) in a 5% CO2 incubator at 37 °C.
Cell viability was assessed using the Cell Counting Kit-8 (CCK-8) assay. During the logarithmic growth phase, 100 µL of cell suspension was dispensed into 96-well plates, with 3000 cells per well, and the supernatant was removed after incubation at 37 °C for 24 h. Wells containing 100 µL of medium served as blank controls. Determination of viability after 48 h of sample treatment.

3.12. Effect of AHYYD on Melanin and Tyrosinase Synthesis in B16F10 Melanoma Cells

The cells in the logarithmic growth phase were seeded at a density of 2 × 105 cells per well in 6-well plates. After 24 h of culture, the supernatant was removed. The control group, positive control group for Kojic acid (350 μM) [48], and peptide-treated groups (15, 37.5, 75 μM) were established. After 48 h of incubation, the cells were lysed, and the cell lysates were centrifuged at 12,000 rpm for 10 min. The reaction mixture consisted of 50 µL of cell lysate supernatant, 50 µL of 1 mg/mL L-DOPA, and 50 µL of PBS buffer. The mixture was incubated at 37 °C for 1 h. The dopamine formation level was then measured spectrophotometrically at 475 nm. Tyrosinase activity was expressed as a percentage of the blank control [49].
The cell culture stage for melanin content determination was as follows: the cells were washed twice with PBS buffer after 48 h of peptide treatment. Trypsin digestion was performed, and the cells were collected. The lower precipitate was removed by centrifugation for 5 min and the precipitate was lysed by dissolving it in 1 mol/L NaOH containing 10% DMSO for 1 h at 80 °C. Subsequently, the cell lysate was centrifuged at 12,000 rpm for 10 min. A total of 200 µL of the supernatant was transferred to a 96-well plate, and the absorbance was measured at 405 nm. The melanin content was then expressed as a percentage of the control group [49].

3.13. Effect of AHYYD on Antioxidant Enzyme Activity and ROS in B16F10 Cells

The cell culture and treatment methods were the same as those described in Section 3.12. The activities of catalase (CAT), glutathione peroxidase (GSH-Px), and superoxide dismutase (SOD) were determined using the assay kit from Nanjing Jianjian Bioengineering Research Institute (Nanjing, China). The ROS content was determined using the method provided by Biyuntian Bio-Technology Company Limited (Shanghai, China).

3.14. Statistical Analysis

All data are expressed as the mean ± standard deviation (SD) of at least three different experiments. Multiple group comparisons were conducted using one-way analysis of variance (ANOVA) and Duncan’s multiple-range test in IBM SPSS Statistics 27. A p-value less than 0.05 was considered statistically significant. The molecular docking was performed using AutoDock Vina, and the molecular dynamics simulations were performed by GROMACS 2022.3. The graphical abstracts were created using Figdraw 2.0.

4. Conclusions

In the present study, we found that peptides derived from the nacre powder of the Pinctada martensii exhibit significant potential for inhibiting tyrosinase. The inhibitory peptide AHYYD was identified from nacre powder peptide by LC–MS/MS and molecular docking screening. It exhibited tyrosinase inhibition with an IC50 value of 2.012 ± 0.088 mM. Inhibition kinetics and isothermal titration results showed that AHYYD was a reversible competitive inhibitor. Molecular docking and molecular dynamics simulations indicated that AHYYD had a strong binding affinity, with a binding energy of −8.0 kcal/mol. The effects of AHYYD on tyrosinase activity, melanin production, and antioxidant enzyme activities were investigated in mouse B16F10 melanoma cells. The results showed that AHYYD significantly inhibited intracellular tyrosinase activity and melanin content, and had a positive effect on the intracellular antioxidant enzyme system. It can be concluded that the peptide AHYYD extracted from Pinctada martensii has good therapeutic efficacy. This can be further confirmed through animal experiments and utilized for the development of food or cosmetic products.

Author Contributions

Conceptualization, F.L.; methodology, F.L. and H.L.; software, validation, and formal analysis, F.L. and H.L.; investigation, F.L. and H.L.; resources, H.L. and S.X.; data curation and writing—original draft preparation, F.L. and H.L.; writing—review and editing, X.Q., J.G., Z.C., W.C. and H.Z.; visualization, F.L.; supervision, project administration, and funding acquisition, H.L., X.Q., H.Z. and S.X. All authors have read and agreed to the published version of the manuscript.

Funding

This work was supported by the earmarked fund for China Agriculture Research System (CARS-49), Doctoral Startup Project of Guangdong Ocean University (R17082), the Guangdong Province Modern Agricultural Industry Technology System Innovation Team Construction Project (Grant No. 2023KJ146), the Innovative Team Program of High Education of Guangdong Province (2021KCXTD021).

Data Availability Statement

The data shown in this study are contained within the article.

Conflicts of Interest

Shaohe Xie is employed by Guangdong Shaohe Pearl Co., Ltd., and the other authors declare that there are no potential conflicts of interest. Guangdong Shaohe Pearl Co., Ltd. has no role in the study design, collection, analysis, interpretation of data, the writing of this article or the decision to submit it for publication.

References

  1. Deng, Z.; Sun, J.; Wei, H.; Zhao, W.; Chen, M.; Li, Y.; Yu, G.; Wang, Y. Shell colors and microstructures of four pearl oyster species in the South China Sea. Aquac. Rep. 2022, 25, 101214. [Google Scholar] [CrossRef]
  2. Marie, B.; Joubert, C.; Tayale, A.; Zanella-Cleon, I.; Belliard, C.; Piquemal, D.; Cochennec-Laureau, N.; Marin, F.; Gueguen, Y.; Montagnani, C. Different secretory repertoires control the biomineralization processes of prism and nacre deposition of the pearl oyster shell. Proc. Natl. Acad. Sci. USA 2012, 109, 20986–20991. [Google Scholar] [CrossRef] [PubMed]
  3. Mariom; Take, S.; Igarashi, Y.; Yoshitake, K.; Asakawa, S.; Maeyama, K.; Nagai, K.; Watabe, S.; Kinoshita, S. Gene expression profiles at different stages for formation of pearl sac and pearl in the pearl oyster Pinctada fucata. BMC Genom. 2019, 20, 240. [Google Scholar] [CrossRef] [PubMed]
  4. Committee, C.P. Pharmacopoeia of the People’s Republic of China; China Medical Science and Technology Press: Beijing, China, 2020. [Google Scholar]
  5. Zhang, J.-X.; Li, S.-R.; Yao, S.; Bi, Q.-R.; Hou, J.-J.; Cai, L.-Y.; Han, S.-M.; Wu, W.-Y.; Guo, D.-A. Anticonvulsant and sedative-hypnotic activity screening of pearl and nacre (mother of pearl). J. Ethnopharmacol. 2016, 181, 229–235. [Google Scholar] [CrossRef] [PubMed]
  6. Song, Y.; Chen, S.; Li, L.; Zeng, Y.; Hu, X. The Hypopigmentation Mechanism of Tyrosinase Inhibitory Peptides Derived from Food Proteins: An Overview. Molecules 2022, 27, 2710. [Google Scholar] [CrossRef] [PubMed]
  7. Zhang, Z.; Xu, Y.; Lai, R.; Deng, H.; Zhou, F.; Wang, P.; Pang, X.; Huang, G.; Chen, X.; Lin, H.; et al. Protective Effect of the Pearl Extract from Pinctada fucata martensii Dunker on UV-Induced Photoaging in Mice. Chem. Biodivers. 2022, 19, e202100876. [Google Scholar] [CrossRef] [PubMed]
  8. Chiu, H.-F.; Hsiao, S.-C.; Lu, Y.-Y.; Han, Y.-C.; Shen, Y.-C.; Venkatakrishnan, K.; Wang, C.-K. Efficacy of protein rich pearl powder on antioxidant status in a randomized placebo-controlled trial. J. Food Drug Anal. 2018, 26, 309–317. [Google Scholar] [CrossRef] [PubMed]
  9. Loh, X.J.; Young, D.J.; Guo, H.; Tang, L.; Wu, Y.; Zhang, G.; Tang, C.; Ruan, H. Pearl Powder-An Emerging Material for Biomedical Applications: A Review. Materials 2021, 14, 2797. [Google Scholar] [CrossRef] [PubMed]
  10. Cordero, R.J.B.; Casadevall, A. Melanin. Curr. Biol. 2020, 30, R142–R143. [Google Scholar] [CrossRef] [PubMed]
  11. Slominski, R.M.; Sarna, T.; Plonka, P.M.; Raman, C.; Brozyna, A.A.; Slominski, A.T. Melanoma, Melanin, and Melanogenesis: The Yin and Yang Relationship. Front. Oncol. 2022, 12, 842496. [Google Scholar] [CrossRef] [PubMed]
  12. Wang, Y.; Xiong, B.; Xing, S.; Chen, Y.; Liao, Q.; Mo, J.; Chen, Y.; Li, Q.; Sun, H. Medicinal Prospects of Targeting Tyrosinase: A Feature Review. Curr. Med. Chem. 2023, 30, 2638–2671. [Google Scholar] [CrossRef] [PubMed]
  13. Fu, W.; Wu, Z.; Zheng, R.; Yin, N.; Han, F.; Zhao, Z.; Dai, M.; Han, D.; Wang, W.; Niu, L. Inhibition mechanism of melanin formation based on antioxidant scavenging of reactive oxygen species. Analyst 2022, 147, 2703–2711. [Google Scholar] [CrossRef] [PubMed]
  14. Liu, F.; Qu, L.; Li, H.; He, J.; Wang, L.; Fang, Y.; Yan, X.; Yang, Q.; Peng, B.; Wu, W.; et al. Advances in Biomedical Functions of Natural Whitening Substances in the Treatment of Skin Pigmentation Diseases. Pharmaceutics 2022, 14, 3154. [Google Scholar] [CrossRef] [PubMed]
  15. Saeedi, M.; Eslamifar, M.; Khezri, K. Kojic acid applications in cosmetic and pharmaceutical preparations. Biomed. Pharmacother. 2019, 110, 582–593. [Google Scholar] [CrossRef] [PubMed]
  16. He, M.; Fan, M.; Yang, W.; Peng, Z.; Wang, G. Novel kojic acid-1,2,4-triazine hybrids as anti-tyrosinase agents: Synthesis, biological evaluation, mode of action, and anti-browning studies. Food Chem. 2023, 419, 136047. [Google Scholar] [CrossRef] [PubMed]
  17. Shaito, A.; Posadino, A.M.; Younes, N.; Hasan, H.; Halabi, S.; Alhababi, D.; Al-Mohannadi, A.; Abdel-Rahman, W.M.; Eid, A.H.; Nasrallah, G.K.; et al. Potential Adverse Effects of Resveratrol: A Literature Review. Int. J. Mol. Sci. 2020, 21, 2084. [Google Scholar] [CrossRef] [PubMed]
  18. Zilles, J.C.; dos Santos, F.L.; Kulkamp-Guerreiro, I.C.; Contri, R.V. Biological activities and safety data of kojic acid and its derivatives: A review. Exp. Dermatol. 2022, 31, 1500–1521. [Google Scholar] [CrossRef] [PubMed]
  19. Masub, N.; Khachemoune, A. Cosmetic skin lightening use and side effects. J. Dermatol. Treat. 2022, 33, 1287–1292. [Google Scholar] [CrossRef]
  20. Hu, S.; Laughter, M.R.; Anderson, J.B.; Sadeghpour, M. Emerging topical therapies to treat pigmentary disorders: An evidence-based approach. J. Dermatol. Treat. 2022, 33, 1931–1937. [Google Scholar] [CrossRef] [PubMed]
  21. Reboucas, J.S.A.; Oliveira, F.P.S.; Araujo, A.C.D.S.; Gouveia, H.L.; Latorres, J.M.; Martins, V.G.; Prentice Hernandez, C.; Tesser, M.B. Shellfish industrial waste reuse. Crit. Rev. Biotechnol. 2023, 43, 50–66. [Google Scholar] [CrossRef] [PubMed]
  22. Zhou, J.; Wei, M.; You, L. Protective Effect of Peptides from Pinctada Martensii Meat on the H2O2-Induced Oxidative Injured HepG2 Cells. Antioxidants 2023, 12, 535. [Google Scholar] [CrossRef] [PubMed]
  23. Wei, M.; Qiu, H.; Zhou, J.; Yang, C.; Chen, Y.; You, L. The Anti-Photoaging Activity of Peptides from Pinctada martensii Meat. Mar. Drugs 2022, 20, 770. [Google Scholar] [CrossRef] [PubMed]
  24. Sasaki, C.; Tamura, S.; Tohse, R.; Fujita, S.; Kikuchi, M.; Asada, C.; Nakamura, Y. Isolation and identification of an angiotensin I-converting enzyme inhibitory peptide from pearl oyster (Pinctada fucata) shell protein hydrolysate. Process Biochem. 2019, 77, 137–142. [Google Scholar] [CrossRef]
  25. Kongsompong, S.; E-Kobon, T.; Taengphan, W.; Sangkhawasi, M.; Khongkow, M.; Chumnanpuen, P. Computer-Aided Virtual Screening and In Vitro Validation of Biomimetic Tyrosinase Inhibitory Peptides from Abalone Peptidome. Int. J. Mol. Sci. 2023, 24, 3154. [Google Scholar] [CrossRef] [PubMed]
  26. Huang, P.; Miao, J.; Liao, W.; Huang, C.; Chen, B.; Li, Y.; Wang, X.; Yu, Y.; Liang, X.; Zhao, H.; et al. Rapid screening of novel tyrosinase inhibitory peptides from a pearl shell meat hydrolysate by molecular docking and the anti-melanin mechanism. Food Funct. 2023, 14, 1446–1458. [Google Scholar] [CrossRef] [PubMed]
  27. Jin, Z.; Wei, Z. Molecular simulation for food protein-ligand interactions: A comprehensive review on principles, current applications, and emerging trends. Compr. Rev. Food Sci. Food Saf. 2024, 23, e13280. [Google Scholar] [CrossRef] [PubMed]
  28. Li, X.; Guo, J.; Lian, J.; Gao, F.; Khan, A.J.; Wang, T.; Zhang, F. Molecular Simulation Study on the Interaction between Tyrosinase and Flavonoids from Sea Buckthorn. Acs Omega 2021, 6, 21579–21585. [Google Scholar] [CrossRef] [PubMed]
  29. Trott, O.; Olson, A.J. Software News and Update AutoDock Vina: Improving the Speed and Accuracy of Docking with a New Scoring Function, Efficient Optimization, and Multithreading. J. Comput. Chem. 2010, 31, 455–461. [Google Scholar] [CrossRef] [PubMed]
  30. Skolnick, J.; Zhou, H. Implications of the Essential Role of Small Molecule Ligand Binding Pockets in Protein-Protein Interactions. J. Phys. Chem. B 2022, 126, 6853–6867. [Google Scholar] [CrossRef] [PubMed]
  31. Wang, W.; Lin, H.; Shen, W.; Qin, X.; Gao, J.; Cao, W.; Zheng, H.; Chen, Z.; Zhang, Z. Optimization of a Novel Tyrosinase Inhibitory Peptide from Atrina pectinata Mantle and Its Molecular Inhibitory Mechanism. Foods 2023, 12, 3884. [Google Scholar] [CrossRef] [PubMed]
  32. Xiong, S.-L.; Lim, G.T.; Yin, S.-J.; Lee, J.; Si, Y.-X.; Yang, J.-M.; Park, Y.-D.; Qian, G.-Y. The inhibitory effect of pyrogallol on tyrosinase activity and structure: Integration study of inhibition kinetics with molecular dynamics simulation. Int. J. Biol. Macromol. 2019, 121, 463–471. [Google Scholar] [CrossRef] [PubMed]
  33. Zhao, W.; Tan, L.; Zhang, Q.; Chen, F.; Yu, Z. In silico identification and mechanistic evaluation of novel tyrosinase inhibitory peptides derived from coconut proteins. Food Biosci. 2024, 61, 104595. [Google Scholar] [CrossRef]
  34. Yu, Z.; Fu, L.; Zhang, Q.; Zhao, W. In silico identification and molecular mechanism of novel egg white-derived tyrosinase inhibitory peptides. Food Biosci. 2024, 57, 103567. [Google Scholar] [CrossRef]
  35. Najafi, Z.; Haramabadi, M.Z.; Chehardoli, G.; Ebadi, A.; Iraji, A. Design, synthesis, and molecular dynamics simulation studies of some novel kojic acid fused 2-amino-3-cyano-4H-pyran derivatives as tyrosinase inhibitors. BMC Chem. 2024, 18, 41. [Google Scholar] [CrossRef] [PubMed]
  36. Lee, Y.-C.; Hsiao, N.-W.; Tseng, T.-S.; Chen, W.-C.; Lin, H.-H.; Leu, S.-J.; Yang, E.-W.; Tsai, K.-C. Phage Display-Mediated Discovery of Novel Tyrosinase-Targeting Tetrapeptide Inhibitors Reveals the Significance of N-Terminal Preference of Cysteine Residues and Their Functional Sulfur Atom. Mol. Pharmacol. 2015, 87, 218–230. [Google Scholar] [CrossRef] [PubMed]
  37. Yu, S.; He, M.; Zhai, Y.; Xie, Z.; Xu, S.; Yu, S.; Xiao, H.; Song, Y. Inhibitory activity and mechanism of trilobatin on tyrosinase: Kinetics, interaction mechanism and molecular docking. Food Funct. 2021, 12, 2569–2579. [Google Scholar] [CrossRef] [PubMed]
  38. Yoshino, M.; Murakami, K. A graphical method for determining inhibition constants. J. Enzym. Inhib. Med. Chem. 2009, 24, 1288–1290. [Google Scholar] [CrossRef] [PubMed]
  39. Whitmore, L.; Mavridis, L.; Wallace, B.A.; Janes, R.W. DichroMatch at the protein circular dichroism data bank (DM@PCDDB): A web-based tool for identifying protein nearest neighbors using circular dichroism spectroscopy. Protein Sci. 2018, 27, 10–13. [Google Scholar] [CrossRef] [PubMed]
  40. Paudyal, S.; Sharma, S.K.; da Silva, R.L.C.G.; Mintz, K.J.; Liyanage, P.Y.; Al-Youbi, A.O.; Bashammakh, A.S.; El-Shahawi, M.S.; Leblanc, R.M. Tyrosinase enzyme Langmuir monolayer: Surface chemistry and spectroscopic study. J. Colloid Interface Sci. 2020, 564, 254–263. [Google Scholar] [CrossRef] [PubMed]
  41. Ott, F.; Rabe, K.S.; Niemeyer, C.M.; Gygli, G. Toward Reproducible Enzyme Modeling with Isothermal Titration Calorimetry. Acs Catal. 2021, 11, 10695–10704. [Google Scholar] [CrossRef]
  42. Liu, J.; Liu, Y.; He, X.; Teng, B.; McRae, J.M. Valonea Tannin: Tyrosinase Inhibition Activity, Structural Elucidation and Insights into the Inhibition Mechanism. Molecules 2021, 26, 2747. [Google Scholar] [CrossRef] [PubMed]
  43. Hu, Z.; Sha, X.; Zhang, L.; Huang, S.; Tu, Z. Effect of Grass Carp Scale Collagen Peptide FTGML on cAMP-PI3K/Akt and MAPK Signaling Pathways in B16F10 Melanoma Cells and Correlation between Anti-Melanin and Antioxidant Properties. Foods 2022, 11, 391. [Google Scholar] [CrossRef] [PubMed]
  44. Soares de Castro, R.J.; Sato, H.H. Biologically active peptides: Processes for their generation, purification and identification and applications as natural additives in the food and pharmaceutical industries. Food Res. Int. 2015, 74, 185–198. [Google Scholar] [CrossRef] [PubMed]
  45. Van der Spoel, D.; Lindahl, E.; Hess, B.; Groenhof, G.; Mark, A.E.; Berendsen, H.J.C. GROMACS: Fast, flexible, and free. J. Comput. Chem. 2005, 26, 1701–1718. [Google Scholar] [CrossRef] [PubMed]
  46. Abraham, M.J.; Murtola, T.; Schulz, R.; Páll, S.; Smith, J.C.; Hess, B.; Lindahl, E. GROMACS: High performance molecular simulations through multi-level parallelism from laptops to supercomputers. SoftwareX 2015, 1, 19–25. [Google Scholar] [CrossRef]
  47. Yu, Q.; Fan, L.; Ding, Z. The inhibition mechanisms between asparagus polyphenols after hydrothermal treatment and tyrosinase: A circular dichroism spectrum, fluorescence, and molecular docking study. Food Biosci. 2022, 48, 101790. [Google Scholar] [CrossRef]
  48. Wang, W.; Gao, Y.; Wang, W.; Zhang, J.; Yin, J.; Le, T.; Xue, J.; Engelhardt, U.H.; Jiang, H. Kojic Acid Showed Consistent Inhibitory Activity on Tyrosinase from Mushroom and in Cultured B16F10 Cells Compared with Arbutins. Antioxidants 2022, 11, 502. [Google Scholar] [CrossRef]
  49. Zhang, X.; Li, J.; Li, Y.; Liu, Z.; Lin, Y.; Huang, J.-A. Anti-melanogenic effects of epigallocatechin-3-gallate (EGCG), epicatechin-3-gallate (ECG) and gallocatechin-3-gallate (GCG) via down-regulation of cAMP/CREB/MITF signaling pathway in B16F10 melanoma cells. Fitoterapia 2020, 145, 104634. [Google Scholar] [CrossRef] [PubMed]
Figure 1. The 3D and 2D visualizations of molecular docking of AHYYD, KPIWT, and TFSGNYP with tyrosinase (2Y9X).
Figure 1. The 3D and 2D visualizations of molecular docking of AHYYD, KPIWT, and TFSGNYP with tyrosinase (2Y9X).
Marinedrugs 22 00359 g001
Figure 2. Molecular dynamics results of tyrosinase with AHYYD, KPIWT, and TFSGNYP. (A) RMSD; (B) RMSF; (C) Rg; (D) SASA; (E) Hbond number.
Figure 2. Molecular dynamics results of tyrosinase with AHYYD, KPIWT, and TFSGNYP. (A) RMSD; (B) RMSF; (C) Rg; (D) SASA; (E) Hbond number.
Marinedrugs 22 00359 g002
Figure 3. Scavenging capacity of peptides for DPPH free radicals (A); Scavenging capacity of peptides for ABTS free radicals (B).
Figure 3. Scavenging capacity of peptides for DPPH free radicals (A); Scavenging capacity of peptides for ABTS free radicals (B).
Marinedrugs 22 00359 g003
Figure 4. (A) Plots of enzymatic reaction rate versus tyrosinase concentration. (B) Plots of Lineweaver–Burk; the secondary plots of slope and Y-intercept versus concentration of AHYYD are shown in the inset.
Figure 4. (A) Plots of enzymatic reaction rate versus tyrosinase concentration. (B) Plots of Lineweaver–Burk; the secondary plots of slope and Y-intercept versus concentration of AHYYD are shown in the inset.
Marinedrugs 22 00359 g004
Figure 5. The CD spectra of tyrosinase in the absence and presence of AHYYD. c (tyrosinase) = 0.4 mg/mL, and c (AHYYD) = 0.6 mM.
Figure 5. The CD spectra of tyrosinase in the absence and presence of AHYYD. c (tyrosinase) = 0.4 mg/mL, and c (AHYYD) = 0.6 mM.
Marinedrugs 22 00359 g005
Figure 6. Interaction of AHYYD with tyrosinase studied by ITC at 37 °C showed the thermogram and binding isotherm (tyrosinase) = 0.001 mM, and c (AHYYD) = 1 mM.
Figure 6. Interaction of AHYYD with tyrosinase studied by ITC at 37 °C showed the thermogram and binding isotherm (tyrosinase) = 0.001 mM, and c (AHYYD) = 1 mM.
Marinedrugs 22 00359 g006
Figure 7. (A) Effects of AHYYD on cell viability. (B) Effects of AHYYD and Kojic acid on tyrosinase activity in B16F10 cells. (C) Effects of AHYYD and Kojic acid on melanin production in B16F10 cells. c (Kojic acid) = 350 μM. Different letters indicate that there are significant differences in data (p < 0.05).
Figure 7. (A) Effects of AHYYD on cell viability. (B) Effects of AHYYD and Kojic acid on tyrosinase activity in B16F10 cells. (C) Effects of AHYYD and Kojic acid on melanin production in B16F10 cells. c (Kojic acid) = 350 μM. Different letters indicate that there are significant differences in data (p < 0.05).
Marinedrugs 22 00359 g007
Figure 8. Effect of AHYYD on the intracellular antioxidant capacity in B16F10 cells; (A) CAT; (B) GSH-Px; (C) SOD; (D) GSH; (E) GSSG; (F) MDA; (G) ROS. c (Kojic acid) = 350 μM. Different letters indicate that there are significant differences in data (p < 0.05).
Figure 8. Effect of AHYYD on the intracellular antioxidant capacity in B16F10 cells; (A) CAT; (B) GSH-Px; (C) SOD; (D) GSH; (E) GSSG; (F) MDA; (G) ROS. c (Kojic acid) = 350 μM. Different letters indicate that there are significant differences in data (p < 0.05).
Marinedrugs 22 00359 g008
Table 1. Tyrosinase-inhibitory activity and antioxidant activity of Enzymatic hydrolysis product NP-PMH.
Table 1. Tyrosinase-inhibitory activity and antioxidant activity of Enzymatic hydrolysis product NP-PMH.
Enzymatic Hydrolysis ProductTyrosinase Inhibitory Activity
(IC50, mg/mL)
DPPH Free Radical Scavenging
(1 mg/mL, %)
ABTS Free Radical Scavenging
(1 mg/mL, %)
NP-PMH6.743 ± 0.0671.139 ± 0.23285.050 ± 2.770
Table 2. Potential tyrosinase-inhibitory peptide sequences: screening by mass spectrometry and molecular docking.
Table 2. Potential tyrosinase-inhibitory peptide sequences: screening by mass spectrometry and molecular docking.
NO.SequencesLengthScoreToxicityAffinityGravy
1AHYYD527.12Non−8.0−1.5
2GGFGNW630.94Non−7.8−0.47
3KPIWT520.49Non−7.4−0.52
4TFSGNYP721.74Non−7.3−0.79
5ATFDAI627.86Non−7.21.12
6NRIPN525.39Non−7.2−1.72
7KRSLE521.85Non−7.0−1.78
8HKDGY526.84Non−6.9−2.46
9ERHLGY624.12Non−6.9−1.52
10SIIDEVVA823.32Non−6.91.43
11KDLFF520.31Non−6.90.4
12GHSLTQF735.06Non−6.8−0.29
13FGSLSF624.8Non−6.81.23
14GGSFSVR726Non−6.80.01
15TNNFT520.4Non−6.7−1.12
16LPEEV523Non−6.7−0.12
17SASTTLEE822.84Non−6.6−0.55
18VTANPANT821Non−6.6−0.28
19HSSAHS620.36Non−6.5−1.17
20TNTSNP623.65Non−6.5−1.8
21SDLGGI621.62Non−6.50.53
22MVSLEG620.69Non−6.50.87
23LKGHEDL720.01Non−6.4−0.99
24SIDLYK620.7Non−6.4−0.2
25MDLSHA621.73Non−6.30
26DYQLP521.75Non−6.2−1.22
27KEMQGG620.62Non−6.2−1.63
28HTLESKPNPD1022.99Non−5.9−1.85
29KEPNK520.77Non−5.9−3.28
30TDIIDG624.25Non−5.80.15
31KKQLM520.39Non−5.8−1.12
32MQVTPASA822.86Non−5.80.39
Table 3. Molecular docking sites of AHYYD, KPIWT, and TFSGNYP with tyrosinase (2Y9X).
Table 3. Molecular docking sites of AHYYD, KPIWT, and TFSGNYP with tyrosinase (2Y9X).
PeptidesHydrogen BondsHydrophobic InteractionElectrostatic Interaction
AHYYDSer282,Cys83 Glu322,Asn260Val283,Pro284,His263
KPIWTMet280,Asn81Ala286,His263Val283,Pro284
TFSGNYPAsn81,Ser282,Gly281,Val283,Ala246,Gly245His263,Ala286,Pro284
Table 4. Effects of AHYYD, TFSGNYP, KPIWT, and positive control Kojic acid on tyrosinase activity.
Table 4. Effects of AHYYD, TFSGNYP, KPIWT, and positive control Kojic acid on tyrosinase activity.
SamplesIC50 (mM)
AHYYD2.012 ± 0.088
TFSGNYP>5
KPIWT>10
Kojic acid0.01 ± 0.003
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Li, F.; Lin, H.; Qin, X.; Gao, J.; Chen, Z.; Cao, W.; Zheng, H.; Xie, S. In Silico Identification and Molecular Mechanism of Novel Tyrosinase Inhibitory Peptides Derived from Nacre of Pinctada martensii. Mar. Drugs 2024, 22, 359. https://doi.org/10.3390/md22080359

AMA Style

Li F, Lin H, Qin X, Gao J, Chen Z, Cao W, Zheng H, Xie S. In Silico Identification and Molecular Mechanism of Novel Tyrosinase Inhibitory Peptides Derived from Nacre of Pinctada martensii. Marine Drugs. 2024; 22(8):359. https://doi.org/10.3390/md22080359

Chicago/Turabian Style

Li, Fei, Haisheng Lin, Xiaoming Qin, Jialong Gao, Zhongqin Chen, Wenhong Cao, Huina Zheng, and Shaohe Xie. 2024. "In Silico Identification and Molecular Mechanism of Novel Tyrosinase Inhibitory Peptides Derived from Nacre of Pinctada martensii" Marine Drugs 22, no. 8: 359. https://doi.org/10.3390/md22080359

APA Style

Li, F., Lin, H., Qin, X., Gao, J., Chen, Z., Cao, W., Zheng, H., & Xie, S. (2024). In Silico Identification and Molecular Mechanism of Novel Tyrosinase Inhibitory Peptides Derived from Nacre of Pinctada martensii. Marine Drugs, 22(8), 359. https://doi.org/10.3390/md22080359

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop