2. Purification of 1-h Trypsin Hydrolysate by Ultrafiltration
1-h Trypsin hydrolysate (T1H) was separated by ultrafiltration (UF), yielding two UF fractions. The >3 kDa fraction exhibited stronger antioxidant activity than the <3 kDa fraction based on their EC
50 for 2,2′-azino-bis(3-ethylbenzothiazoline-6-sulfonic acid) diammonium salt radical cation (ABTS
•+) (28.4 and 45.6 µg dry mass (DM)/mL, respectively) and H
2O
2 (174.5 and 461.7 µg DM/mL, respectively) scavenging activities (
Figure 1C,D). Thiobarbituric acid reactive species (TBARS) value of the negative control increased from 0.7 to 1.6 µM malondialdehyde equivalents after incubation of the lecithin liposomes from 24 to 48 h, indicating an increase in lipid peroxidation. As revealed by the TBARS values, treatment with 0.5 mg DM/mL of >3 kDa and <3 kDa fractions inhibited lipid peroxidation by 22% and 16% after 24 h, and by 64% and 50% after 48 h, respectively (
Figure 1E). Our results are concordant with previous observations on the antioxidant activities of the UF fractions of fennel seed hydrolysate
[10]. In theirs and this study, the superiority of the >3 kDa fraction in scavenging radicals relative to the <3 kDa fraction was observed. This may be attributed, in part, to the presence of large compounds, such as long peptides, partially degraded proteins, or other components with antioxidant properties in the >3 kDa fraction. In fact, the protein content of the >3 kDa fraction was 8-fold higher than that of the <3 kDa fraction (
Figure 1A). By contrast, the <3 kDa fraction was peptide-rich, with peptide content 5-fold greater than the >3 kDa fraction (
Figure 1B).
Figure 1. (A) Protein content, (B) peptide content, (C) 2,2′-azino-bis(3-ethylbenzothiazoline-6-sulfonic acid) diammonium salt radical cation (ABTS•+) scavenging activity and (D) hydrogen peroxide (H2O2) scavenging activity of ultrafiltration (UF) fractions. (E) Thiobarbituric acid reactive species (TBARS) values of negative control and UF fractions tested at 0.5 mg dry mass (DM)/mL. For each bar chart, data are mean ± standard errors (n = 3). Mean values denoted by different superscript letters are significantly different (p < 0.05) according to Tukey’s test.
Notably, the <3 kDa fraction was comparable or superior to T1H in scavenging ABTS
•+ and repressing lipid peroxidation, based on activity data reported in our previous study
[4]. This implies that at least part of the antioxidant activities exhibited by T1H could be attributed to the presence of antioxidant peptides in the hydrolysate. The <3 kDa fraction was weaker than antioxidant tripeptide GSH (EC
50 6.7 µg DM/mL) as an ABTS
•+ scavenger but 3-fold stronger than GSH (EC
50 1378.7 µg DM/mL) as H
2O
2 scavenger. H
2O
2 is an ROS molecule that can diffuse through biological membranes and convert to highly reactive hydroxyl radicals in the body cells
[5]. Hence, the <3 kDa UF fraction likely contained antioxidant peptides with the ability to scavenge biologically relevant ROS. In this study, a time-dependent increase in the lipid peroxidation inhibitory activity of the UF fractions was observed (
Figure 1E). This observation agrees with a previous study that compared the ability of corn gluten meal hydrolysates to inhibit lipid peroxidation in a ground pork system following 8 h and 16 h of treatment
[11]. Our result showed that the lipid peroxidation inhibitory effects of the CS peptides could persist up to 48 h. Low molecular weight (MW) peptides are desirable because they could be more easily absorbed by the body compared to high MW peptides
[5]. Thus, the peptide-enriched <3 kDa fraction was selected for further purification.
3. Purification of <3 kDa Fraction by Gel Filtration Chromatography
Purification by gel filtration chromatography (GFC) resulted in three pooled fractions: GF-I, GF-II, and GF-III (
Figure 2A). Among the three, GF-I exhibited the strongest effects in ABTS
•+ scavenging activity (
Figure 2C). GF-I possibly comprised more non-aromatic peptide residues with radical scavenging activity, such as Leu and Pro
[12] than the other two fractions. Meanwhile, GF-III showed the highest H
2O
2 scavenging activity among the three pooled fractions (
Figure 2D). All three pooled GFC fractions could dampen the time-dependent increase in lipid peroxidation in the liposome model at 0.1 mg peptide/mL, with 26–35% inhibition of TBARS formation after 48 h (
Figure 2E), although their activities were significantly lower than that of GSH (57% inhibition after 48 h). The peptide content of GF-III (0.68 mg peptide/mg DM) was 3.5-fold greater than those of GF-I and GF-II (
Figure 2B). When compared to GF-III, 3.5-fold higher DM of GF-I and GF-II was required to achieve the standardized peptide concentration used for evaluating the lipid peroxidation inhibitory activity depicted in
Figure 2E. It can be anticipated that when expressed in terms of DM, the lipid peroxidation inhibitory activity of GF-III may exceed that of GF-I. Thus, we also analyzed the lipid peroxidation inhibitory activity of the three pooled fractions at 0.5 mg DM/mL. As expected, among the three fractions, GF-III showed the strongest inhibition of TBARS formation, with 26% and 51% inhibition after 24 and 48 h, respectively (
Figure 2F). The antioxidant activity of GF-III could be owing to its relatively high absorbance at 280 nm (
Figure 2A), which suggests an abundance of aromatic amino acid residues (e.g., Phe, Tyr, and Trp) in the pooled fraction. In keeping with this study, a previous study on Chinese chestnut also found that the GFC fraction with the most prominent absorbance at 280 nm had the highest antioxidant activity among all GFC fractions
[13].
Figure 2. Purification of the <3 kDa fraction by gel filtration chromatography (GFC). (A) Elution profile, (B) peptide content, (C) ABTS•+ scavenging activity (at 30 µg peptide/mL), and (D) H2O2 scavenging activity (at 150 µg peptide/mL) of pooled GFC fractions. TBARS values of negative control and GFC fractions (E) at 0.1 mg peptide/mL and (F) at 0.5 mg DM/mL. Data are mean ± standard errors (n = 3). For each bar chart, mean values denoted by different superscript letters are significantly different (p < 0.05) according to Tukey’s test.
In our GFC experiment, GF-II (intermediate molecular size) had the lowest antioxidant activity, whereas GF-I (greatest molecular size) and GF-III (smallest molecular size) had relatively higher antioxidant activities. Hence, our results suggest that the antioxidant potential of the GFC fractions is not directly related to their molecular size. This is in agreement with the lack of explicit relationship between antioxidant activity and molecular size among 81 corn gluten meal peptide fractions collected in a GFC experiment
[14]. On the other hand, peptide content data and the different trends in the lipid peroxidation inhibitory activities of the pooled fractions when tested based on DM and peptide mass pointed to the presence of non-peptide constituents in the GFC fractions. Thus, further purification was desirable. Considering that GF-III was the richest in peptide content and to discover peptides containing aromatic amino acid residues from T1H, we proceeded to perform purification on GF-III.
4. Purification of GF-III by Strong-Cation-Exchange Solid-Phase Extraction
GF-III was further purified by strong-cation-exchange solid-phase extraction (SCX-SPE), producing six SPE fractions. As shown in Figure 3A, most of the peptide constituents of GF-III were found in 50 mM KCl fraction (0.9 mg peptide/mL), which was 22–118 times greater than the other five SPE fractions. The 50 mM KCl fraction showed relatively low or no ABTS•+ and H2O2 scavenging activities at the peptide concentrations tested (Figure 3B,C). Our result suggests that SCX-SPE has partitioned most of the non-antioxidant peptides and/or peptides with weak antioxidant activity into the 50 mM KCl fraction.
Figure 3. Purification of the GF-III fraction by strong-cation-exchange solid-phase extraction (SPE). (A) Peptide concentration, (B) ABTS•+ scavenging activity (at 0.7 µg peptide/mL), and (C) H2O2 scavenging activity (at 10 µg peptide/mL) of SPE fractions. Data are mean ± standard errors (n = 3). For each bar chart, mean values denoted by different superscript letters are significantly different (p < 0.05) according to Tukey’s test.
The effectiveness of SCX-SPE in concentrating the antioxidant peptides into single fractions were also evident, as affirmed by the enhancement in antioxidant activity after the SPE step. Briefly, 30 µg peptide/mL of GF-III scavenged 57% of ABTS
•+ (
Figure 2C). After the purification of GF-III by SCX-SPE, the resultant 0 and 20 mM KCl fractions scavenged 35% and 22% ABTS
•+ at a 43-fold lower concentration (0.7 µg peptide/mL), respectively (
Figure 3B). By estimation, the 0 and 20 mM KCl fractions may be 26-fold and 17-fold stronger than GF-III as ABTS
•+ scavengers, respectively. Other studies also showed 33%
[15] and 44%
[10] improvement in the ABTS
•+ scavenging activity of a peptide fraction purified by the SCX-SPE. The potency of the 0 and 20 mM KCl fractions over the other SPE fractions as ABTS
•+ scavenger may be attributed to the presence of negatively charged amino acids (e.g., Glu)
[16] or proton-donating amino acids (e.g., Trp and Gln)
[14]. The presence of such residues may impart antioxidant activity to peptides by transferring electrons or protons to free radicals
[12]. For instance, rapeseed peptides predominantly made up of Glu (19.5%) were reported to have potent radical scavenging activity
[16]. Corn gluten meal-derived peptides made up of 67% of Trp and Gln displayed high ABTS
•+ scavenging activity
[14].
A similar improvement in the antioxidant activity of SPE fractions following SCX-SPE was revealed by the H2O2 scavenging assay. For example, 150 µg peptide/mL of GF-III scavenged 65% H2O2 (Figure 2D). In contrast, at a 15-fold lower concentration (10 µg peptide/mL), the 0 and 200 mM KCl fractions scavenged more than 80% of H2O2 (Figure 3C). The 20 mM KCl fraction also scavenged 56% H2O2 when tested at 10 µg peptide/mL. Thus, based on theoretical calculations, our results imply a 13–21-fold improvement in the H2O2 scavenging activity of the 0, 20, and 200 mM KCl fractions resulting from the purification of GF-III by using SCX-SPE. Altogether, the 0, 20, and 200 mM KCl fractions potentially contained potent antioxidant peptides; thus, they were taken to peptide sequencing.
5. Identification and Characterization of Antioxidant Peptides
Liquid chromatography-tandem mass spectrometry analysis identified 29 peptide sequences comprising 6–14 residues (633.33 to 1517.81 Da) from the 0, 20, and 200 mM KCl fractions (
Table 1). This range of peptide masses agrees with the observation that the molecular masses of food-derived antioxidant peptides commonly range between 500–1800 Da
[12]. Twenty-three of the 29 peptides contain 11–56% aliphatic amino acid residues (
Table 1). Such residues are responsible for the thermal stability of proteins
[17]. Two thermal-stable antioxidant peptides WAFAPA and MYPGLA that were identified from the blue-spotted stingray, for instance, are composed of 50% and 33% of aliphatic residues, respectively
[18]. Based on the comparison of the aliphatic index, 10 of the 29 peptides were likely superior to both WAFAPA and MYPGLA in terms of thermal stability (
Table 1). The discovery of such peptides also supports our previous observation of the thermal stability of T1H, the protein hydrolysate from which the 29 peptides were purified. T1H retained its radical scavenging and ferric reducing activity at temperatures up to 100 °C
[4]. These CS peptides that are likely to be thermal-stable can thus be utilized as alternatives for food additives to address the concerns regarding the food processing heat treatment.
Table 1. Physicochemical properties of the peptides identified from the 0, 20, and 200 mM KCl fractions.
SPE Fractions |
Peptides |
Measured m/z [M + 2H]2 |
Molecular Mass (Da) a |
Aromatic Residues (%) b |
Basic Residues (%) b |
Hydrophobic Residues (%) b |
Aliphatic Residues (%) b |
Aliphatic Index b |
0 mM KCl |
KRYFKR |
449.28 |
896.57 |
33 |
67 |
33 |
0 |
0 |
PRVRVAGR |
455.79 |
909.58 |
0 |
38 |
63 |
38 |
85 |
PVVWAAKR |
463.79 |
925.57 |
13 |
25 |
75 |
50 |
98 |
QVASGPLQR |
478.28 |
954.55 |
0 |
11 |
56 |
33 |
87 |
MAPRTPRK |
478.78 |
955.57 |
0 |
38 |
50 |
13 |
13 |
NKVVKLMR |
494.31 |
986.62 |
0 |
38 |
50 |
38 |
121 |
KVPLAVFSR |
508.82 |
1015.64 |
11 |
22 |
67 |
44 |
119 |
LKKGSPLKR |
513.84 |
1025.69 |
0 |
44 |
44 |
22 |
87 |
FQLKPVFR |
517.82 |
1033.63 |
25 |
25 |
63 |
25 |
85 |
THAVKGVVHK |
538.34 |
1074.67 |
20 |
40 |
50 |
40 |
97 |
YTWKFKGR |
543.31 |
1084.61 |
38 |
38 |
50 |
0 |
0 |
ARVPQQSYR |
552.80 |
1103.61 |
11 |
22 |
44 |
22 |
43 |
VHFNKGKKR |
557.34 |
1112.69 |
22 |
56 |
33 |
11 |
32 |
TAPLSSKALKR |
586.37 |
1170.73 |
0 |
27 |
45 |
36 |
89 |
FSCPLVMKGPNGLR |
759.91 |
1517.81 |
7 |
14 |
71 |
21 |
76 |
20 mM KCl |
RHGSGR |
335.18 |
668.37 |
17 |
50 |
33 |
0 |
0 |
NMVPGR |
337.17 |
672.34 |
0 |
17 |
67 |
17 |
48 |
FMFFVYK |
491.25 |
980.50 |
57 |
14 |
86 |
14 |
41 |
MCFHHHFHK |
612.27 |
1222.53 |
67 |
56 |
44 |
0 |
0 |
200 mM KCl |
DFPGAK |
317.66 |
633.33 |
17 |
17 |
67 |
17 |
17 |
NDGPSR |
323.15 |
644.29 |
0 |
17 |
33 |
0 |
0 |
AGFPLGK |
345.20 |
688.41 |
14 |
14 |
86 |
29 |
70 |
AMQQDK |
360.66 |
719.32 |
0 |
17 |
33 |
17 |
17 |
NLEGYR |
376.19 |
750.38 |
17 |
17 |
50 |
17 |
65 |
YETLNR |
398.20 |
794.41 |
17 |
17 |
33 |
17 |
65 |
MPPKSTR |
408.72 |
815.43 |
0 |
29 |
43 |
0 |
0 |
TAGASLVAR |
423.25 |
844.49 |
0 |
11 |
67 |
56 |
109 |
SSPATGGSLR |
466.74 |
931.49 |
0 |
10 |
50 |
20 |
49 |
NANSLAGPQR |
514.27 |
1026.55 |
0 |
10 |
50 |
30 |
59 |
Amino acid composition is a key factor influencing the antioxidant activity of peptides
[5]. Hydrophobic residues made up 33–86% of the compositions of the 29 CS peptides (
Table 1). The presence of hydrophobic residue-containing peptides may account for the lipid peroxidation inhibitory effect of T1H
[4], as well as that of <3 kDa UF (
Figure 1E) and GF-III (
Figure 2E,F), the peptide fractions which the 29 peptides were purified from. Hydrophobic residues may enhance the interaction of antioxidant peptides with lipid-soluble free radicals, thus attenuating the progression of lipid peroxidation
[5]. Lipid oxidation is one of the major contributors to the deterioration of food quality during food processing and storage
[5]. Thus, these CS peptides, as well the CS hydrolysate and partially purified fractions containing them, may be useful in the preservation of lipid-rich foods. On the other hand, it has been reported that peptides with C-terminal Lys could act as potent H
2O
2 scavengers
[19]. Thus, MAPRTPRK and THAVKGVVHK in the 0 mM KCl fraction, FMFFVYK and MCFHHHFHK in the 20 mM KCl fraction, as well as DFPGAK, AGFPLGK, and AMQQDK in the 200 mM KCl fraction may have contributed to the H
2O
2 scavenging activity of the three SPE fractions.
The 0 mM KCl fraction had the highest proportion of peptides containing aromatic residues (9 peptides), followed by the 200 mM KCl (4 peptides) and the 20 mM KCl fractions (3 peptides). Aromatic residues may promote the antioxidant activity of peptides by donating protons to the electron-deficient radicals
[12]. Notably, the 0 mM KCl fraction had the highest number of basic residue-containing peptides and the highest percentage of basic residues in peptides (
Table 1). Thus, the strongest ABTS
•+ scavenging activity of the 0 mM KCl fraction may be attributed to its richness in peptides comprising aromatic and basic amino acids. Our result agrees with the finding that abundance in basic amino acids may account for the strong ABTS
•+ scavenging activity of a Chinese chestnut peptide fraction
[13].
As revealed by the AnOxPePred analysis, 10 CS peptides (MCFHHHFHK, NLEGYR, AGFPLGK, FMFFVYK, NMVPGR, PVVWAAKR, DFPGAK, FSCPLVMKGPNGLR, RHGSGR, and VHFNKGKKR) had comparable or higher free radical scavenger (FRS) scores relative to the four reference peptides VGPWQK, MYPGLA, FPLPSF and WAFAPA (
Table 2). The four reference peptides were empirically proven as ABTS
•+ scavengers
[18][15][20]. Our results are in accordance with the finding that His, Trp, Tyr, and Pro are common in free-radical-scavenging peptides
[21]. In this study, the four amino acids account for 11–44% of the residues making up the 10 CS peptides and the four reference peptides (
Table 2). Notably, His-containing peptides were only found in the 0 mM KCl (THAVKGVVHK and VHFNKGKKR) and 20 mM KCl (RHGSGR and MCFHHHFHK) fractions, but none in the 200 mM KCl fraction (
Table 1). This is in keeping with our observation that the 0 and 20 mM KCl fractions had at least 4-fold greater ABTS
•+ scavenging activity than the 200 mM KCl (
Figure 3B). Our results, therefore, support the role of His residues in imparting radical scavenging activity to peptides
[12].
Table 2. Free radical scavenger (FRS) scores of corn silk peptides were identified from the three SPE fractions, in comparison with reference peptides.
Peptides |
SPE Fractions |
FRS Scores |
MCFHHHFHK |
20 mM KCl |
0.68068 |
VGPWQK * |
- |
0.52254 |
MYPGLA * |
- |
0.49386 |
NLEGYR |
200 mM KCl |
0.48158 |
AGFPLGK |
200 mM KCl |
0.44866 |
FMFFVYK |
20 mM KCl |
0.44397 |
NMVPGR |
20 mM KCl |
0.44319 |
PVVWAAKR |
0 mM KCl |
0.43744 |
DFPGAK |
200 mM KCl |
0.43574 |
FPLPSF * |
- |
0.43352 |
FSCPLVMKGPNGLR |
0 mM KCl |
0.41864 |
WAFAPA * |
- |
0.41519 |
RHGSGR |
20 mM KCl |
0.41088 |
VHFNKGKKR |
0 mM KCl |
0.41055 |
NANSLAGPQR |
200 mM KCl |
0.40415 |
QVASGPLQR |
0 mM KCl |
0.40213 |
MAPRTPRK |
0 mM KCl |
0.39973 |
NDGPSR |
200 mM KCl |
0.38760 |
KRYFKR |
0 mM KCl |
0.38352 |
YETLNR |
200 mM KCl |
0.37938 |
FQLKPVFR |
0 mM KCl |
0.37599 |
ARVPQQSYR |
0 mM KCl |
0.37580 |
YTWKFKGR |
0 mM KCl |
0.36769 |
AMQQDK |
200 mM KCl |
0.36324 |
SSPATGGSLR |
200 mM KCl |
0.35382 |
THAVKGVVHK |
0 mM KCl |
0.35200 |
MPPKSTR |
200 mM KCl |
0.33529 |
LKKGSPLKR |
0 mM KCl |
0.32957 |
PRVRVAGR |
0 mM KCl |
0.32698 |
KVPLAVFSR |
0 mM KCl |
0.32525 |
TAGASLVAR |
200 mM KCl |
0.32285 |
TAPLSSKALKR |
0 mM KCl |
0.29320 |
NKVVKLMR |
0 mM KCl |
0.27437 |
6. Molecular Docking between CS Peptides and ABTS•+
Docking simulation was performed to clarify the interactions between ABTS
•+ and the 10 CS peptides with the best FRS scores. All seven peptides originating from the 0 and 20 mM KCl fractions had higher binding affinities towards ABTS
•+ than the two peptides (AGFPLGK and DFPGAK) from the 200 mM KCl fraction (
Table 3). Five of the seven peptides (from 0 and 20 mM KCl fractions) were also stronger than NLEGYR (from 200 mM KCl fraction) in binding to ABTS
•+. The overall trend is in accordance with the relative levels of in vitro ABTS
•+ scavenging activity of the three SPE fractions (
Figure 3B). Notably, the binding affinities of all seven peptides from the 0 and 20 mM KCl fractions were more negative than that of reference peptides MYPGLA. The binding energy of three peptides MCFHHHFHK, VHFNKGKKR, and PVVWAAKR was up to 21% more negative than all four reference peptides (
Table 3). Taken together, the seven peptides originating from the 0 and 20 mM KCl fractions could bind to ABTS
•+ similarly or more stably than could the four reference peptides. Peptides that bind stably to free radicals can neutralize them. For instance, FPLPSF that was predicted to bind to ABTS
•+ has been experimentally demonstrated to quench ABTS
•+ in vitro
[20]. Furthermore, our prediction of WAFAPA binding to ABTS
•+ more stably than could MYPGLA (
Table 3) is also consistent with their relative in vitro antioxidant activity
[18]. Altogether, our results suggest that the stronger ABTS
•+ scavenging activities of the 0 and 20 mM KCl fractions, relative to that of the 200 mM KCl fraction, could be accounted for, at least in part, by the affinity of their seven peptides to ABTS
•+.
Table 3. Binding affinities and types of interactions between 10 corn silk peptides and ABTS•+, in comparison with four reference peptides.
Peptides |
SPE Fractions |
Binding Affinity (kcal/mol) |
Peptide Residues Interacting with ABTS•+ a |
Hydrogen Bond |
Hydrophobic Interaction |
MCFHHHFHK |
20 mM KCl |
−4.8 |
- |
Phe3, His6, Phe7 |
VHFNKGKKR |
0 mM KCl |
−4.7 |
Lys7, Arg9 |
Val1, His2, Gly6, Lys7, Arg9 |
PVVWAAKR |
0 mM KCl |
−4.7 |
Arg8 (2) |
Val2, Trp4, Ala5, Ala6, Arg8 |
FMFFVYK |
20 mM KCl |
−4.4 |
Lys7 |
Phe1, Phe3, Phe4, Lys7 |
FSCPLVMKGPNGLR |
0 mM KCl |
−4.2 |
Arg14 (2) |
Leu5, Lys8, Gly9, Pro10, Gly12, Arg14 |
NMVPGR |
20 mM KCl |
−4.1 |
Asn1, Arg6 (2) |
Asn1, Pro4, Gly5, Arg6 |
NLEGYR |
200 mM KCl |
−4.1 |
- |
Tyr5, Arg6 |
RHGSGR |
20 mM KCl |
−3.9 |
Arg1, Arg6 |
Arg1, Gly5, Arg6 |
AGFPLGK |
200 mM KCl |
−3.7 |
- |
Phe3, Pro4, Leu5 |
DFPGAK |
200 mM KCl |
−3.6 |
- |
Pro3, Gly4, Lys6 |
FPLPSF * |
- |
−4.6 |
Phe1, Ser5 |
Phe1, Pro2, Leu3, Pro4, Ser5 |
WAFAPA * |
- |
−4.3 |
- |
Trp1, Ala4, Pro5 |
VGPWQK * |
- |
−3.9 |
- |
Pro3, Trp4, Lys6 |
MYPGLA * |
- |
−3.8 |
Pro3 |
Pro3, Leu5, Ala6 |
Our LigPlot+ analysis indicates significant participation of aromatic residues in CS peptide-ABTS
•+ interactions (
Table 3). For instance, aromatic residues in 67% of the aromatic residue-containing peptides could bind to ABTS
•+ through hydrophobic interactions. Aromatic residue-ABTS
•+ interactions made up 14–100% of the total number of interactions between individual aromatic residue-containing peptides and ABTS
•+. Remarkably, all interactions formed between the best-binding-affinity MCFHHHFHK and ABTS
•+ were contributed by the aromatic residues Phe and His (
Table 3). Four CS peptides, namely VHFNKGKKR, PVVWAAKR, FSCPLVMKGPNGLR, and NMVPGR, were predicted to have the highest number of interactions with ABTS
•+. Notwithstanding, the lack of participation of aromatic residues in the interactions between ABTS
•+ and the two peptides FSCPLVMKGPNGLR and NMVPGR apparently made their binding to ABTS
•+ 12–15% less stable relative to VHFNKGKKR and PVVWAAKR. Similarly, both WAFAPA and MYPGLA were predicted to form comparable numbers of interactions with ABTS
•+. However, the lack of participation of aromatic residues in MYPGLA-ABTS
•+ interaction may explain the reported weaker ABTS
•+ scavenging activity of MYPGLA when compared with WAFAPA
[18]. Moreover, we also observed the significant participation of basic residues in peptide-ABTS
•+ interactions. Briefly, basic residues in 90% of basic residue-containing CS peptides were involved in the interactions with ABTS
•+. Such interactions account for 33–80% of the total number of interactions between individual basic residue-containing peptides and ABTS
•+. To further verify the role of the basic residues of peptides in binding to ABTS
•+, in silico alanine substitution was performed on those that were involved in ABTS
•+-peptide interactions, followed by docking of the alanine-substituted peptides to ABTS
•+. The binding affinities of all CS peptides, except FSCPLVMKGPNGLR, were diminished upon alanine mutagenesis of selected basic residues in the peptides (
Table 3 and
Table 4). For instance, alanine substitution of His6 in MCFHHHFHK has resulted in a 15% reduction in the binding affinity of peptides towards ABTS
•+. Besides, a decline of 24% in binding affinity of NMVPGR towards ABTS
•+ was observed upon alanine substitution of Arg6. Our results suggest that the basic residues are likely to be critical in binding and stabilizing ABTS
•+. This observation further reinforces our finding that the 0 mM KCl fraction with the highest number of basic residue-containing peptides displayed the strongest ABTS
•+ scavenging activity (
Table 1 and
Figure 3B). Besides, Leu-ABTS
•+ interactions were observed in 67% of Leu-containing CS peptides (
Table 3). Our observation agrees with a previous report of the participation of Leu in the binding between antioxidant peptides and ABTS
•+ [20].
Table 4. Binding affinities of corn silk peptides toward ABTS•+ upon alanine substitution of the basic residues that were involved in ABTS•+-peptide interactions.
Peptides a |
Basic Residues |
Mutant Peptides |
Binding Affinity (kcal/mol) |
MCFHHHFHK |
His6 |
MCFHHAFHK |
−4.1 |
VHFNKGKKR |
His2 |
VAFNKGKKR |
−4.8 |
Lys7 |
VHFNKGAKR |
−5.0 |
Arg9 |
VHFNKGKKA |
−4.4 |
PVVWAAKR |
Arg8 |
PVVWAAKA |
−4.3 |
FMFFVYK |
Lys7 |
FMFFVYA |
−4.3 |
FSCPLVMKGPNGLR |
Lys8 |
FSCPLVMAGPNGLR |
−4.2 |
Arg14 |
FSCPLVMKGPNGLA |
−4.7 |
NMVPGR |
Arg6 |
NMVPGA |
−3.1 |
NLEGYR |
Arg6 |
NLEGYA |
−3.7 |
RHGSGR |
Arg1 |
AHGSGR |
−3.8 |
Arg6 |
RHGSGA |
−4.2 |
DFPGAK |
Lys6 |
DFPGAA |
−3.4 |
7. Molecular Docking of Peptides on Keap1
Food-derived bioactive peptides, in addition to scavenging free radicals, can confer cellular protection by modulating the gene expression and activities of antioxidant and oxidant enzymes
[5]. Given this, we conducted a docking-based screening experiment to unravel the potential of the 29 CS peptides identified in this study in interacting with cellular protein targets that can regulate the endogenous oxidant status: Keap1, MPO, and XO. Soy-derived DEQIPSHPPR was predicted in molecular docking study to interact stably with Keap1, in keeping with its demonstrated ability to disrupt Keap1-Nrf2 binding and increase Nrf2 levels in the nucleus
[6]. Hence, DEQIPSHPPR was used as a reference peptide for comparison with CS peptides. Our docking results show that 13 of the 29 CS peptides could bind to Keap1 similarly or more stably than DEQIPSHPPR (
Table S2 (could be found in
https://www.mdpi.com/2076-3921/10/11/1822#supplementary)). Further in silico screening for low toxicity and allergenicity as well as high cell-penetrating potential narrowed down the 13 CS peptides to five, namely NDGPSR, NLEGYR, NMVPGR, SSPATGGSLR, and NANSLAGPQR (
Table 5). Screening based on these parameters allows the search for CS peptides that might be able to cross the cell membrane barrier and block the Keap1-Nrf2 interaction in cells with minimal or no harmful effects. Unlike the five CS peptides, the reference peptide DEQIPSHPPR may elicit allergy (
Table 5), hence it is less desirable for the application of functional food ingredients.
Table 5. Toxicity, allergenicity, and cell-penetrating ability predicted for selected corn silk peptides that have the same or higher affinity to Kelch-like ECH-associated protein 1, myeloperoxidase, and xanthine oxidase, in comparison with reference peptides.
Peptides |
Toxicity |
Allergenicity |
CPP Prediction |
NDGPSR |
Non-toxin |
Probable non-allergen |
CPP |
NLEGYR |
Non-toxin |
Probable non-allergen |
CPP |
NMVPGR |
Non-toxin |
Probable non-allergen |
CPP |
SSPATGGSLR |
Non-toxin |
Probable non-allergen |
CPP |
NANSLAGPQR |
Non-toxin |
Probable non-allergen |
CPP |
KRYFKR |
Non-toxin |
Probable non-allergen |
CPP |
RHGSGR |
Non-toxin |
Probable non-allergen |
CPP |
YETLNR |
Non-toxin |
Probable non-allergen |
Non-CPP |
AGFPLGK |
Non-toxin |
Probable non-allergen |
Non-CPP |
KVPLAVFSR |
Non-toxin |
Probable non-allergen |
Non-CPP |
TAGASLVAR |
Non-toxin |
Probable allergen |
Non-CPP |
YTWKFKGR |
Non-toxin |
Probable allergen |
CPP |
AMQQDK |
Non-toxin |
Probable allergen |
CPP |
MPPKSTR |
Non-toxin |
Probable allergen |
CPP |
PVVWAAKR |
Non-toxin |
Probable allergen |
CPP |
DFPGAK |
Non-toxin |
Probable allergen |
Non-CPP |
FMFFVYK |
Non-toxin |
Probable allergen |
Non-CPP |
QVASGPLQR |
Non-toxin |
Probable allergen |
Non-CPP |
DEQIPSHPPR * |
Non-toxin |
Probable allergen |
Non-CPP |
DTETGVPT * |
Non-toxin |
Probable non-allergen |
Non-CPP |
VPY * |
Non-toxin |
Probable allergen |
CPP |
ACECD * |
Non-toxin |
Probable allergen |
CPP |
The five aforementioned CS peptides successfully docked into Keap1 and interacted with 3–6 of the seven key residues known to be involved in Keap1-Nrf2 interactions (
Table 6). Our analysis on the reference peptide DEQIPSHPPR agrees with a previous report of its interaction with the key residues Arg380, Asn382, and Arg415 of Keap1
[6]. The participation of the three residues was also observed in the interactions between Keap1 and CS peptides (NANSLAGPQR and SSPATGGSLR) (
Table 6). Besides, the binding of NLEGYR to Keap1 via hydrogen bond, hydrophobic interaction, and the salt bridge was also found in the DEQIPSHPPR-Keap1 interaction (
Table 6). Thus, the five CS peptides apparently have comparable Keap1-binding properties as DEQIPSHPPR and could potentially activate the Keap1/Nrf2 pathway, triggering cellular antioxidant defense.
Table 6. Binding affinities and types of interactions between Kelch-like ECH-associated protein 1 (Keap1) and five corn silk peptides predicted as non-toxic, non-allergenic and cell-penetrating peptides, in comparison with a reference peptide.
Peptides |
Binding Affinity (kcal/mol) |
Interaction with Keap1 a |
|
Hydrogen Bond |
Hydrophobic Interaction |
Salt Bridge |
NLEGYR |
−8.7 |
Arg415, Arg483, Ser508, Gln530, Ser555 |
Tyr334, Ser363, Gly364, Leu365, Ala366, Arg415, Ile416, Gly417, Gly462, Phe478, Arg483, Ser508, Gly509, Ala510, Tyr525, Gln530, Ser555, Ala556, Leu557, Tyr572, Phe577, Ser602, Gly603, Val604 |
Arg415 |
NANSLAGPQR |
−8.2 |
Arg415 (3), Val418, Val465, Arg483 |
Ser363, Gly364, Leu365, Arg380, Asn382, Asn414, Arg415, Ile416, Gly417, Ile461, Gly462, Val463, Val465, Phe478, Arg483, Ser508, Gly509, Tyr525, Gln530, Ser555, Ala556, Ile559, Phe577, Gly603 |
- |
NMVPGR |
−8.1 |
Ser363, Leu365, Asn382, Ser602 |
Tyr334, Ser363, Gly364, Leu365, Ala366, Asn382, Arg415, Ile416, Ile461, Gly462, Ser508, Gly509, Ala510, Tyr525, Gln530, Ser555, Ala556, Ser602 |
- |
SSPATGGSLR |
−8.1 |
Ser363, Arg380, Asn414, Arg415, Ser431, Ser602 |
Tyr334, Gly364, Leu365, Arg380, Asn382, Asn414, Arg415, Ile416, Ser431, Gly433, His436, Gly462, Phe478, Arg483, Ser508, Gly509, Ala556, Ser602, Gly603 |
- |
NDGPSR |
−8.0 |
Arg415 (2), Ala510 |
Tyr334, Gly364, Leu365, Arg415, Ile461, Gly462, Phe478, Ser508, Gly509, Tyr525, Ala556, Ser602, Gly603, Val604 |
- |
DEQIPSHPPR * |
−8.0 |
Tyr334, Asn414, Arg415 (4), Ser431, Arg483 (3), Ser555 |
Tyr334, Ser363, Arg380, Asn382, Asn414, Arg415, Ser431, Gly433, His436, Gly462, Phe478, Arg483, Ser508, Gly509, Tyr525, Ser555, Ala556, Tyr572, Phe577, Ser602 |
Arg483 (2) |