0% READ
FULL TEXT
INTRODUCTION
Tyrosinase (EC1.14.18.1) is a key enzyme in the process of melanogenesis. Tyrosinase catalyses the rate-limiting steps of hydroxylation of L-tyrosine to 3,4-dihydroxyphenylalanine (L-DOPA) and further oxidization of L-DOPA to produce melanin.1 Melanin is responsible for skin, hair and eye colour and plays a critical role in protecting the skin against ultraviolet radiation. Over accumulation of melanin is associated with disorders, such as freckles, melasma, solar lentigo, post-inflammatory hyperpigmentation and neurodegeneration like Parkinson’s disease.2 Tyrosinase inhibitors have found application in cosmetics and medicine. However, many known tyrosinase inhibitors are found to be cytotoxic, unstable and carcinogenic.3 Hence, there is a need for novel non-toxic, effective and stable tyrosinase inhibitory agents.
Silk cocoon shell is composed of two main proteins, sericin (20-30% w/w) and fibroin (70-80% w/w) and some trace amounts of non-protein components (1-5% w/w). The non-protein components include pigments, waxes, sugars, metabolites and certain inorganic compounds. The sericin-containing component has reported to have several valuable properties like antioxidant, antitumor and antimicrobial activities.4 However, the composition, preparation and biological activities of the non-protein component of the sericin solution haven’t been looked into yet.5 In the present investigation, the non-protein component of the sericin extract from the cocoon shells of Antheraea proylei J was determined using GC-MS chromatography. In silico and in vitro tyrosinase inhibition was performed to check the possibility of furthering the non-sericin component as a potential tyrosinase inhibitor.
MATERIALS AND METHODS
Materials
Cocoon shells of Antheraea proylei J. were acquired from tasar farm, Mortbung, Manipur, India. Ethanol from Jebsen and Jessen GmbH and Co. (Germany). Mushroom tyrosinase, kojic acid were from Sigma Chemical Co. (USA).3,4-dihydroxy-L-phenylalanine (L-DOPA) was from HiMedia (India).
Methods
Preparation of the non-protein component of sericin extract
The sericin extract from the cocoon shells of oak tasar silkworm was prepared using hot water extraction method.6 The cocoon shells were autoclaved and the aqueous solutions obtained were centrifuged at 8000 g for 15 min to remove any suspended insoluble materials. The resulting final aqueous solution was the sericin extract. The non-protein component in the sericin extract was fractionated using ethanol precipitation. Sericin extract and ethanol was mixed in the ratio of 1:8. The extract was fractionated using ethanol precipitation. The mixture was kept overnight at -20°C, followed by centrifugation at 15000 g for 30 min at 40°C. The resulting supernatant contained the non-protein component which was concentrated by rotary evaporation.
GC-MS analysis
Gas Chromatography-Mass Spectrometry (GC-MS) analysis of the ethanolic extract was carried out using the GC-MS (Shimadzu QP-2010 Plus with Thermal Desorption System TD 20). RTx-5MS column (30 m in length, 0.25 mm internal diameter and 0.25 μm film thickness) was used. The analysis was done by injecting 1 μL of the sample with a split ratio of 10.0. The column oven temperature was set at 80°C for 0.5 min then increased from 80°C to 300°C at a rate of 10°C/min. Helium gas was used as the carrier gas. The column flow rate was maintained at 1.2 mL/min and the injection temperature was 260°C. Compounds were identified by comparing the spectral data of peaks with the corresponding standard mass spectra from the library database.7
In silico analysis for tyrosinase inhibition
Preparation of ligands and protein for docking
The structures of the compounds identified from GS-MS were obtained from PubChem database. The Open Babel and Autodock formats were used to obtain the structures in.pdbqt format.8,9 The tyrosinase protein (PDB ID- 4OUA) was converted to its.pdbqt format.
Screening of compounds for binding with tyrosinase
The Autodock Vina software was utilized to screen the compounds identified by GC-MS analysis.10 The compound showing the best binding affinity was further subjected to Autodock 4.0 analysis. The resultant protein-ligand interactions were analysed using the Discovery Studio software.9
Molecular dynamics simulation studies
The 1,3-Bis(cinnamoyloxymethyl)adamantine-tyrosinase complex was subjected to molecular dynamics simulation to check for stability in docking. Molecular dynamics simulation was performed by GROMACS (Groningen MAchine for Chemical Simulations) version 2022.3. The initial structure is prepared using pdb2gmxby adding hydrogens. CHARMM36 force field is used to describe the interactions and ligand topology.11 Trajectory analysis is performed using the inbuilt GROMACS programs. All visual analysis is performed using Visual Molecular Dynamics (VMD) software.12
Binding free energy analysis by MM/PBSA
g_mmpbsa program is used to evaluate the binding energy profile of protein and ligand.13 Residue-wise free energy decomposition is derived using the g_mmpbsa associated python program, MmPbSaDecomp.py script.
In vitro antityrosinase assay
The antityrosinase activity of the ethanolic extract of the non-protein component was determined according to Kim et al., 2005.14 The reaction mixture was constituted by mixing 2 U of mushroom tyrosinase, 8 mM DOPA and the sample in 100 mM phosphate buffer, pH 6.8, followed by incubation for 5 min at room temperature. The amount of dopachrome produced was measured at 480 nm. Kojic acid was taken as a positive control. The enzyme activity was expressed as per cent enzyme activity.
IC50determination
The concentration of the sample required for 50% inhibition of the enzyme (IC50) was determined using OriginPro 8. The results were expressed as Mean±Standard Deviation (SD) of three independent experiments.
RESULTS
The composition of the non-protein component of the sericin extract from the cocoon of Antheraea proylei J. was identified by GC-MS analysis as described in GC-MS analysis. The GC-MS spectra are shown in Figure 1. The compounds identified by the GC-MS analysis are mentioned in Table 1. Thirty-four metabolites were identified and it includes polyamines, phenolics, cholesterol derivatives, pyrazole, pyrazole derivatives, fatty acids and fatty acid derivatives.
| Sl. No. | Compound name | Binding energy (kcal/mol) |
|---|---|---|
| 1 | 3,3’-Iminobispropylamine | -4.4 |
| 2 | Phenol, 2,4-Bis(1,1-Dimethylethyl)- | -6.7 |
| 3 | 3-Methyl-1,4-diazabicydo[4.3.0]nonan-2,5-dione, N-acetyl | -6.7 |
| 4 | 4-(1,1-Dioxidotetrahydro-3-thienyl)morpho | -5.4 |
| 5 | Pyrrolo[1,2-a]pyrazine-1,4-dione, hexahydro- | -7.1 |
| 6 | 8-Octadecanone | -5.4 |
| 7 | Pyrrolo[1,2-a]pyrazine-1,4-dione, hexahydro-3-(2-methylpr | -6.4 |
| 8 | n-Hexadecanoic acid | -5.1 |
| 9 | N,N’-Dibutylidene-hydrazine | -4.2 |
| 10 | 9-Octadecenoic acid (Z)- | -5.6 |
| 11 | Hexadecanamide | -5.6 |
| 12 | 9,12-Octadecadienoic Acid (Z,Z)- | -5.6 |
| 13 | 9-Octadecenamide, (Z)- | -5.4 |
| 14 | Octadecanamide | -5.1 |
| 15 | octadecanoic acid, 3-oxo-, ethyl ester | -5.8 |
| 16 | 2-Ethylbutyric acid, tridecyl ester | -5.2 |
| 17 | 9-t-Butyltricyclo[4.2.1.1(2,5)]decane-9,10-diol | -6.8 |
| 18 | 1-Cyclohexyldimethylsilyloxybutane | -4.6 |
| 19 | 4-(2-Tert-Butyl-5-oxo-1,3-Dioxolan-4-yl)butylformamide | -6.1 |
| 20 | Octadecanoic acid, 2,3-dihydroxypropyl ester | -5.2 |
| 21 | 4,7,10,13,16,19-Docosahexaenoic acid, methyl ester, (all-Z | -6.8 |
| 22 | 9-Octadecenamide | -5.4 |
| 23 | 1,3-Bis(cinnamoyloxymethyl)adamantine | -8.4 |
| 24 | 5,8,11,14-Eicosatetraenoic acid, ethyl ester | -6.7 |
| 25 | Cholesta-3,5-diene | -7.9 |
| 26 | 9-Octadecenamide | -4.8 |
| 27 | 1,37-Octatriacontadien | -5.1 |
| 28 | Cholesterol | -8 |
| 29 | 8-Nonenoic acid, 9-(1,3-nonadienyloxy)-, methyl ester | -5.6 |
| 30 | Z-11-Tetradecen-1-ol trifluoroacetate | -7.2 |
| 31 | Betulinaldehyde | -8 |
| 32 | Palmitic acid vinyl ester | -5 |
| 33 | 1-Pentatriacontanol | -5.2 |
| 34 | Tetrahydrosarsasapogenin | -8.1 |
The Autodock Vina-based screening results for the compounds against the tyrosinase enzyme showed that the compound 1,3-Bis(cinnamoyloxymethyl)adamantine showed the best binding affinity for the active site of tyrosinase (Table 1). This compound was subjected to Autodock 4.0 analysis. 1,3-Bis(cinnamoyloxymethyl)adamantine showed a binding energy of -7.18 kcal/mol (Table 2) and hydrogen bonding with THR268 amino acid residue (Figure 2a).
| Compound | Binding energy (kcal/mol) | Inhibition constant (μM) | Ligand efficiency | Number of Hydrogen bonds |
|---|---|---|---|---|
| 1,3-Bis(cinnamoyloxymethyl)adamantine | -7.18 | 5.42 | -0.27 | 1 |
Molecular dynamics simulation of the 1,3-Bis(cinnamoyloxymethyl)adamantine and tyrosinase complex showed an increase in the Root Mean Square Deviation (RMSD) value compared to the protein structure alone (Figure 2b). The Radius of gyration (Rg) stayed steady for tyrosinase both in the presence and absence of the ligand (Figure 2c). The Root Mean Square Fluctuation (RMSF) study of the amino acids in the tyrosinase enzyme shows a restriction in the motion of the amino acids between positions 75 to 85 and 250 to 270. The RMSF values (Figure 2d) correspond with the 2D image (Figure 2a) of the interaction between the ligand and HIS82 and the amino acids between positions 250 to 273 of the protein. Solvent Accessible Surface Area (SASA) of the protein alone and in the presence of the ligand didn’t show any difference (Figure 2e). The number of hydrogen bonds between the protein and ligand stayed at a maximum of 2 (Figure 3a). The total number of pairs of hydrogen bonds in the entire structure has significantly reduced in the protein in the presence of ligand (Figure 3b).
The Molecular Mechanics/Poisson-Boltzmann Surface Area (MM/PBSA) analysis reveals that van der Waals non-bonded interactions have contributed most to the ligand-protein interactions (Table 3). From Figure 3c we observed that the most dominant contributions are residues GLU309, GLU248, GLU307 and the residues in the positions of 75 to 85 and 250 to 270.
| Energy components | Energy (kJ/mol) |
|---|---|
| van der Waal energy (ΔEvdw) | -78.370±64.666 kJ/mol |
| Electrostatic energy (ΔEvdw) | -20.103±27.626 |
| Polar solvation energy (ΔGsol) | 49.750± 67.177 |
| SASA energy (ΔGSASA) | -9.691±7.989 |
| Binding energy (ΔGbinding) | -58.415±43.331 |
The ethanolic extract also caused a gradual reduction in the tyrosinase activity with increasing doses of the extract (Figure 3d). The extract showedin vitro antityrosinase activity with an IC50 value of 3.808±0.920 μg/mL (Table 4).
| Sample | IC50 (μg/mL) |
|---|---|
| Kojic acid | 0.006±0.002 |
| Ethanolic extract | 3.808±0.920 |
DISCUSSION
The identity and composition of these secondary metabolites found associated with the cocoon shell depend on the nature of the leaf of the food plant that the silkworm ingests.15 Hence different types of silk cocoon shells exhibited different metabolite composition. These metabolites are derived in the body of the silkworm. Hence the cocoon shell contained metabolites not present in the food plants which could be a source of finding novel bioactive agents.16
The non-protein fraction from the cocoon shell of Antheraea proylei J. was found to exhibit in vitro antityrosinase activity (IC50=3.808±0.920 μg/mL) (Table 4). An attempt to identify metabolites contributing to in vitro antityrosinase activity was made through Autodock Vina screening. Out of all the metabolites, the compound 1,3-Bis(cinnamoyloxymethyl) adamantine was found to exhibit the best binding affinity for the active site of the enzyme tyrosinase (Table 1). 1,3-Bis(cinnamoyloxymethyl)adamantine formed bonds with conserved residues like HIS57 and HIS255.17 Molecular dynamics simulation of the 1,3-Bis(cinnamoyloxymethyl)adamantine and tyrosinase complex showed an increase in the Root Mean Square Deviation (RMSD) value compared to the protein structure alone (Figure 2b). The Radius of gyration (Rg) stayed steady for tyrosinase both in the presence and absence of the ligand (Figure 2c). The increase in the RMSD value indicates the increased interaction between the protein and the ligand. Rg values showed the overall compactness of the structure.18 The lowering of RMSF values of the active site amino acids (Figure 2a and 2d) in proteins corresponded with the decreased protein activity.19 The total number of pairs of hydrogen bonds in the entire structure has significantly reduced in the protein in the presence of ligand (Figure 3a and 3b) which shows reduced structural stability of the protein induced by the presence of ligand.20 The in vitro results lend further credence to our in silico analysis of the potential of the non-protein component of silkworm extract to act as a tyrosinase inhibitor.
CONCLUSION
Our results demonstrated that the ethanolic extract of the non-protein component of the sericin extract was composed of several metabolites. Of those components, 1,3-Bis(cinnamoyloxymethyl)adamantine showed the best binding with the conserved amino acids of the tyrosinase active site. The stability of the binding was further demonstrated by the molecular dynamics simulation studies. In vitro inhibition of the activities of the tyrosinase enzyme by the ethanolic extracts of the non-sericin component showed further validation of the in silico studies.
