Next Article in Journal
Copper-Catalyzed Ring-Opening Reactions of Alkyl Aziridines with B2pin2: Experimental and Computational Studies
Next Article in Special Issue
Naringenin Induces ROS-Mediated ER Stress, Autophagy, and Apoptosis in Human Osteosarcoma Cell Lines
Previous Article in Journal
Enhancement of Antibiotic Activity by 1,8-Naphthyridine Derivatives against Multi-Resistant Bacterial Strains
Previous Article in Special Issue
The Implication of Low Dose Dimethyl Sulfoxide on Mitochondrial Function and Oxidative Damage in Cultured Cardiac and Cancer Cells
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
Molecules
Volume 26
Issue 23
10.3390/molecules26237396
molecules-logo
Submit to this Journal Review for this Journal Propose a Special Issue
► Article Menu

Article Menu

share Share announcement Help format_quote Cite question_answer Discuss in SciProfiles thumb_up
...
Endorse textsms
...
Comment
first_page
settings
Order Article Reprints
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Computational Screening for the Anticancer Potential of Seed-Derived Antioxidant Peptides: A Cheminformatic Approach

by
Tsun-Thai Chai
1,2,*,
Jiun-An Koh
1,
Clara Chia-Ci Wong
1,
Mohamad Zulkeflee Sabri
3 and
Fai-Chu Wong
1,2
1
Department of Chemical Science, Faculty of Science, Universiti Tunku Abdul Rahman, Kampar 31900, Malaysia
2
Center for Agriculture and Food Research, Universiti Tunku Abdul Rahman, Kampar 31900, Malaysia
3
Green Chemistry and Sustainable Technology Cluster, Bioengineering Section, Universiti Kuala Lumpur Malaysian Institute of Chemical and Bioengineering Technology (UniKL MICET), Lot 1988, Bandar Vendor Taboh Naning, Alor Gajah 78000, Malaysia
*
Author to whom correspondence should be addressed.
Molecules 2021, 26(23), 7396; https://0-doi-org.brum.beds.ac.uk/10.3390/molecules26237396
Submission received: 18 November 2021 / Revised: 1 December 2021 / Accepted: 2 December 2021 / Published: 6 December 2021
(This article belongs to the Special Issue Redox Active Molecules in Cancer Treatments)
Browse Figures
Versions Notes

Abstract

:
Some seed-derived antioxidant peptides are known to regulate cellular modulators of ROS production, including those proposed to be promising targets of anticancer therapy. Nevertheless, research in this direction is relatively slow owing to the inevitable time-consuming nature of wet-lab experimentations. To help expedite such explorations, we performed structure-based virtual screening on seed-derived antioxidant peptides in the literature for anticancer potential. The ability of the peptides to interact with myeloperoxidase, xanthine oxidase, Keap1, and p47phox was examined. We generated a virtual library of 677 peptides based on a database and literature search. Screening for anticancer potential, non-toxicity, non-allergenicity, non-hemolyticity narrowed down the collection to five candidates. Molecular docking found LYSPH as the most promising in targeting myeloperoxidase, xanthine oxidase, and Keap1, whereas PSYLNTPLL was the best candidate to bind stably to key residues in p47phox. Stability of the four peptide-target complexes was supported by molecular dynamics simulation. LYSPH and PSYLNTPLL were predicted to have cell- and blood-brain barrier penetrating potential, although intolerant to gastrointestinal digestion. Computational alanine scanning found tyrosine residues in both peptides as crucial to stable binding to the targets. Overall, LYSPH and PSYLNTPLL are two potential anticancer peptides that deserve deeper exploration in future.

1. Introduction

The past decade has seen a surge in scientific interest towards the exploration of bioactive peptides for potential applications in health promotion and disease management. Bioactive peptides identified from plant food and other natural origins often range between 2 and 20 residues, although this is not a hard-and-fast definition as exceptions do exist [1,2,3]. Plant bioactive peptides, known to exhibit diverse bioactivities, such as antioxidant, antihypertensive, antimicrobial, and antitumor activities, are often purified and identified from enzymatic hydrolysates of edible plant sources and plant-based agricultural by-products. The bioactive potency of some such peptides have also been demonstrated in cellular and animal models [2,4,5]. To date, a growing body of research has shown that plant seeds are a good source of antioxidant peptides [2,5,6]. While such peptides could be developed into natural additive for food processing and nutraceuticals for health maintenance, they may also be therapeutically relevant as some could modulate cellular and/or in vivo antioxidant status [2]. Cellular redox homeostasis is connected to the initiation and/or progression of certain cancers [7,8]. Perturbation in reactive oxygen species (ROS) homeostasis resulting from unchecked ROS production is associated with carcinogenesis; scavenging of excessive ROS accumulation may prevent early neoplasia [9]. Significant reduction in the antioxidant activity of the blood serum of patients with malignant neoplasms has also been reported [10].
In the body, cellular redox status is regulated by oxidative and antioxidative enzymes, non-enzymatic antioxidants, and certain protein-protein interactions involved in regulating antioxidant gene expression. Myeloperoxidase (MPO), xanthine oxidase (XO), and nicotinamide adenine dinucleotide phosphate oxidase (NADPH oxidase) are three examples of such oxidative enzymes. MPO, an abundant heme-containing enzyme in the human neutrophils, catalyzes the reaction between hydrogen peroxide and chloride, generating hypochlorous acid, a potent oxidant. MPO-mediated oxidative burst has been linked to the initiation and progression of cancer, including tumor cell metastasis. Notably, downregulation of MPO gene expression is connected to reduction in the risk of lung, breast, and ovarian cancers [7,8]. XO is an enzyme that catalyzes the conversion of hypoxanthine to xanthine and ultimately to uric acid, producing ROS during the reaction. The importance of XO as an anticancer target is highlighted by the discovery that XO inhibitor febuxostat could repress breast cancer cell migration and the metastasis of breast cancer to the lung in animal models [11,12]. Six isoforms of NADPH oxidase are known to date. NADPH oxidase is a membrane-bound enzyme complex in phagocytes, whose primary function is the production of superoxide anion radicals. The assembly and activation of NADPH oxidase requires protein-protein interaction between the cytosolic factor p47phox and transmembrane component p22phox [13,14]. Due to the importance of p47phox-p22phox interaction in NADPH oxidase activation, the interaction can be targeted in structure-based virtual screening for NADPH oxidase inhibitors [15]. Notably, enhanced NADPH oxidase expression in multiple malignant diseases supports the recognition of the NADPH oxidase family as potential targets in cancer therapies [13,16]. The Kelch-like ECH-associated protein 1 (Keap1)-nuclear factor E2-related factor 2 (Nrf2) pathway is one of the major signaling cascades involved in protecting cells against oxidative stress. The Nrf2 transcription factor can activate the transcription of cytoprotective genes implicated in protection against cancer. However, Keap1-Nrf2 protein-protein interaction could trigger Nrf2 degradation mediated by the ubiquitin–proteasome pathway. Hence, there has been strong interest among researchers to discover inhibitors of Keap1-Nrf2 protein-protein interaction. Such inhibitors may preserve or enhance the transcription-activating role of Nrf2, counteracting ROS-mediated damage in cancers [17,18].
Although a growing number of seed-derived antioxidant peptides has been documented in the literature, knowledge of their ability to modulate cellular regulators of oxidative status (i.e., MPO, XO, NADPH oxidase, and Keap1-Nrf2), which are also promising targets of anticancer therapy, is still limited. A recent report of watermelon seed-derived antioxidant peptides targeting the Keap1-Nrf2 system [19] suggests that seed-derived antioxidant peptides should be explored more intensively as potential modulators of cellular regulators of ROS balance. Thus, this in silico study was undertaken to virtually screen the numerous seed-derived antioxidant peptides in the literature for their potential as anticancer peptides that can target two oxidative enzymes (MPO and XO) and two protein-protein interactions (Keap1-Nrf2 and p47phox-p22phox). In silico or virtual screening is a less costly and less time-consuming strategy to screen for desirable bioactive peptides and other compounds when compared with wet-lab screening [20]. In bioactive peptide screening, this approach can benefit from various freely available peptide databases (e.g., PlantPepDB [21], and other online tools, such as AntiCP 2.0 [22] and MLCPP [23], which are anticancer peptide and cell-penetrating potential prediction servers designed from machine learning models). Moreover, different molecular modelling and simulation methods [24,25,26] may also be used to clarify the mechanisms of action between the peptides and the protein targets of interest. Although virtual screening cannot replace wet-lab experimentation, the aforementioned benefits have driven increasing popularity of in silico research on bioactive peptides [20,27]. Notably, by narrowing down a large set of candidate peptides to a small number, in silico screening can facilitate a more focused research strategy in future wet-lab experimentation; this also allows more efficient use of limited research resources [20].
The goal of this in silico study was three-fold: (a) to compile a virtual library of seed antioxidant peptides from the literature, followed by screening for non-toxic, non-allergenic, and non-hemolytic anticancer peptides; (b) to perform structure-based screening of the predicted anticancer peptides for ability to target Keap1-Nrf2, MPO, XO, and p47phox-p22phox, followed by molecular dynamics validation of peptide-target interactions; and (c) to further characterize the predicted anticancer peptides based on computational alanine mutagenesis and prediction of cell- and blood-brain barrier penetrating potential, as well as plasma and gastrointestinal (GI) stability.

2. Results and Discussion

A virtual library consisting of 677 seed-derived antioxidant peptides was generated (Table S1), based on peptide sequences collected from Scopus and PlantPepDB databases, as outlined in Materials and Methods. The collection encompassed antioxidant peptides of 2–57 residues in length and 192–5338 Da in molecular mass. Seed sources in the virtual library included legumes, such as faba bean and soybean; cereals, such as wheat and rye; and seeds of plantation crop species, such as oil palm and coconut. The types of antioxidant activities reported for the seed-derived peptides included in vitro free radical scavenging activities, lipid peroxidation inhibitory activity, cellular antioxidant activity, and in vivo antioxidant activity (Table S1). Based on Figure 1a, 52% of the seed-derived antioxidant peptides contain five to ten residues. By contrast, seed-derived antioxidant peptides with more than 20 residues comprised only 0.15–0.59% of the virtual antioxidant peptide library. Among the 63 Scopus-indexed publications we examined for the preparation of the virtual library, 42 (67%) reported peptides of 5–10 residues. The prevalence of such peptide length could be accounted by many seed-derived antioxidant peptides being purified and identified from protein fractions of a relatively low molecular mass range, such as <3 kDa fractions [28,29,30].
Next, we proceeded to screening the virtual library for potential anticancer peptides. Only 592 peptides were screened since dipeptides, tripeptides, peptides with more than 50 residues, and peptides with unnatural or modified amino acid residues could not be analyzed by the AntiCP 2.0 tool. Among the 592 peptides, only 123 (21%) were predicted as anticancer (Figure 1b). We carefully examined the publications reporting the 123 peptides and found that none of the peptides had been tested for anticancer activity experimentally. The 123 predicted anticancer peptides averaged 6 residues in length and 750 Da in mass (data not shown). Safety is an important consideration in the design or discovery of anticancer peptides. A functional anticancer peptide should not exhibit toxicity, elicit immune response, and induce the lysis of erythrocytes [31,32,33]. Our screening found that at least 50% of the 592 seed-derived antioxidant peptides were predicted to be safe (i.e., non-toxic, non-allergenic, and non-hemolytic) (Figure 1b). Among the 592 peptides screened, only 0.7% (4 peptides) were predicted to be toxic (Figure 1b). An in silico study also predicted that all 253 antioxidant peptides liberated from the flaxseed proteome were non-toxic [34]. This agrees with our observation of high abundance (99%) of non-toxic peptides in our antioxidant peptide virtual library (Figure 1b). In comparison with toxicity prediction, 47–48% of our antioxidant peptide virtual library comprised allergenic and hemolytic peptides. In an in silico study of 26 antimicrobial peptides of rapeseed, 54% were predicted as non-allergenic and 46% allergenic [35]. This relative distribution of allergenicity and non-allergenicity resembles that observed in our virtual screening. Among the 308 non-hemolytic peptides (Figure 1b), the greatest proportion (41%) originated from legumes, which included soybean and chickpea (data not shown). Fourteen soybean-derived multifunctional cationic peptides were shown to have no hemolytic effect on sheep red blood cells [36]. Meanwhile, two chickpea-derived antioxidant peptides also did not cause any hemolysis in bovine red blood cells [37]. These findings support our observation of legumes being a potential source of non-hemolytic peptides.
Based on our in silico screening, five seed-derived antioxidant peptides were predicted to be anticancer, non-toxic, non-allergenic, and non-hemolytic. The two-dimensional (2D) structures and molecular weight of the five peptides are shown in Figure 2. The five peptides, identified from chickpeas, cherry seeds, and tomato seeds, are 5–9 residues in length and 615–1016 Da in mass. The five peptides each contain at least one imidazole functional group or one aromatic ring among their amino acid side chains. Notably, LPHFNS and LYSPH each contain both an imidazole functional group and an aromatic ring in their structures. This is characteristic of many food-derived antioxidant peptides; imidazole groups and aromatic rings are associated with the ability of the peptides to scavenge free radicals by electron transfer/proton donation [38]. On the other hand, among the five peptides (Figure 2), FGPEMEQ has Phe (F) at the N-terminus, whereas Leu (L), His (H), and Phe (F) are present in four, three, and two of the peptides, respectively. The N-terminal preference for Phe and the abundance of Leu, His, and Phe are both characteristics of experimentally validated anticancer peptides [22].
In this in silico study, to investigate whether the five predicted anticancer peptides (Figure 2) could modulate cellular targets of cancer treatments, we docked the five peptides on Keap1, MPO, XO, and p47phox. To the best of our knowledge, structure-based virtual screening of the five peptides on the four targets has not been reported. Molecular docking analysis found that LYSPH, a cherry seed peptide, had the strongest binding affinity to Keap1, whereas PSYLNTPLL, a tomato seed peptide, had the weakest (Table 1). LYSPH, LPHFNS, and AEHGSLH also had binding affinity values more negative than that of ETGE (−7.1 kcal/mol) (data not shown). ETGE is the key motif of the co-crystalized 16-mer Nrf2 peptide that is involved in Keap1-Nrf2 interaction [39]. Thus, LYSPH, LPHFNS, and AEHGSLH could form similarly stable or more stable binding to Keap1 when compared with Nrf2. Furthermore, all five peptides could bind to the key residues of Keap1 that are required for stable Keap1-Nrf2 complex formation, mostly accomplished via hydrogen bonds and hydrophobic interactions (Table 1). Two tripeptides (DKK and DDW) that could bind to the key residues of Keap1 have been shown experimentally to inhibit Keap1-Nrf2 interaction in vitro [40]. DKK, which possessed a stronger activity than DDW, was reported to bind to key residues Arg380 and Asn382 [40]. Similar to DKK, all five seed-derived peptides in Table 1 were predicted to bind to Arg380 and Asn382. Thus, our binding affinity and intermolecular interaction results suggest that LYSPH, LPHFNS, and AEHGSLH may serve as potential inhibitors of Keap1-Nrf2 interaction. Specifically, at the molecular level, LYSPH was predicted to bind with the same Keap1 residues as did ETGE, namely, Arg380, Arg415, Arg483, and Ser508 [39]. This observation, in addition to LYSPH having the most negative binding affinity to Keap1 among the five peptides, suggests that the peptide is the most promising for targeting Keap1-Nrf2 interaction. A graphical representation of a LYSPH-Keap1 docked model and the intermolecular interactions between LYSPH and Keap1 is shown in Figure 3.
LYSPH showed the strongest binding to MPO, whereas PSYLNTPLL the weakest (Table 2), similar to our observations when the five peptides were docked to Keap1 (Table 1). None of the peptides showed better binding affinity to MPO than did 7-benzyl-1H-[1–3]triazolo[4,5-b]pyridin-5-amine (7GD)(−7.1 kcal/mol) (data not shown), a co-crystalized inhibitor of MPO [41]. However, all five peptides could form hydrophobic interactions with one of the catalytic residues (Arg239) of MPO. Besides, all five peptides could interact with the heme group (Hec606) through hydrophobic interactions (Table 2); the heme group is a cofactor in the active site of MPO [42]. Based on the interactions with both catalytic residue Arg239 and the heme group of MPO, all the five peptides are potential MPO inhibitors. Supporting this possibility is the finding that two experimentally-validated anti-MPO peptides (TDY and FAPQY) could also bind to Arg239 and the heme group of MPO [43]. Analysis of intermolecular interactions revealed that LYSPH could form hydrophobic interactions with Phe99, Thr238, Arg239, Glu242, Phe366, Phe407, and Hec606 of MPO (Table 2), resembling to the binding pattern of 7GD [41]. Hence, LYSPH is the most promising MPO inhibitor among the five peptides as it showed the best binding affinity to MPO and could interact with MPO similarly as 7GD. A graphical representation of LYSPH-MPO docked model and the intermolecular interactions between LYSPH and MPO is depicted in Figure 4.
Comparison of binding affinities found LYSPH (−6.2 kcal/mol) to have the most stable binding to XO among the five peptides analyzed (Table 3). Nevertheless, all of the five peptides had less negative binding affinities to XO than quercetin (−8.2 kcal/mol) (data not shown), a co-crystalized inhibitor of XO [12]. This implies that none of the peptides could bind more stably to XO when compared with quercetin. On the other hand, analysis of intermolecular interactions showed that all five peptides could bind to at least nine of the XO residues known to bind to quercetin, mainly through hydrophobic interactions. Each of the peptides could also bind to at least one catalytic residue (Glu802 or Arg880) of XO [12] through hydrophobic interactions. FGPEMEQ and PSYLNTPLL could also hydrogen bond to Glu802 and Arg880, respectively. However, despite additional interactions with Glu802 and Arg880, FGPEMEQ-XO interaction was predicted to be slightly less favorable than LYSPH-XO interaction based on comparison of their binding affinity values. Meanwhile, PSYLNTPLL-XO interaction was likely non-favorable or non-spontaneous considering the positive value predicted for its binding affinity (Table 3). Previous studies found that experimentally-proven XO-inhibitory peptides, KGFP [45] and EEAK [46] could both bind to the catalytic residue Glu802. Thus, the aforementioned binding patterns of the five seed peptides to XO, particularly their binding to XO catalytic residues, suggest that the peptides are potential XO inhibitors. LYSPH could be the most promising XO inhibitor among the five peptides considering its strongest binding affinity and its binding to XO residues that known XO inhibitors bind to (Table 3). A graphical representation of a LYSPH-XO docked model and the intermolecular interactions between LYSPH and XO is shown in Figure 5.
In the molecular docking to p47phox, tomato seed-derived PSYLNTPLL had the best docking score, whereas the cherry seed-derived FGPEMEQ had the worst (Table 4). All peptides had docking score less negative than that of proline-rich peptide derived from p22phox (−309.862) (data not shown). Thus, none of the five peptides could bind more stably to p47phox than the p22phox-derived peptide. Nevertheless, the potential of the five peptides as inhibitors of p47phox-p22phox interaction could not be completely ruled out solely based on this. Supporting this proposition is the observation that four peptides (RRSSIRNAHSIHQRSRKRLS, ISNSESGPRGVHFIFNKENF, RSRKRLSQDAYRRNSVRFLQQR, and AGGPPGGPQVNPIPVTDEVV) that were experimentally demonstrated to inhibit p47phox-p22phox interaction [47,48] were also predicted to bind less strongly to p47phox than was p22phox (Table S4). In short, a peptide predicted to bind less strongly to p47phox than p22phox may still inhibit p47phox-p22phox interaction.
As shown in Table 4, each of the peptides could bind to at least six of the 17 p47phox residues known to bind to the p22phox-derived peptide. However, only PSYLNTPLL, LYSPH, and FGPEMEQ could interact with Phe209, a key residue of p47phox which accounts for high-affinity binding between p47phox and p22phox (Table 4). Furthermore, PSYLNTPLL could bind to p47phox in a similar manner as the co-crystalized p22phox-derived peptide, by binding to Trp193, Trp204, Pro206, Phe209, Tyr237, Trp263, Met278, and Tyr279. Hence, PSYLNTPLL is the most promising among the five peptides to target p47phox-p22phox interaction considering its docking score and pattern of binding to p47phox. A graphical representation of PSYLNTPLL-p47phox docked model and the intermolecular interactions between PSYLNTPLL and p47phox is shown in Figure 6.
Based on binding affinities and similarity of binding patterns to those of co-crystalized inhibitors/ligands and reported peptide-based inhibitors, our analyses found LYSPH and PSYLNTPLL to have the greatest potential as modulators of the four targets of cancer treatments that we investigated. Specifically, LYSPH may be a multi-target peptide which could bind to, thus inhibiting the activity of MPO and XO, as well as interrupting Keap1-Nrf2 complex formation. On the other hand, PSYLNTPLL is the most promising peptide that could bind to p47phox, thus precluding p47phox-p22phox interaction and the subsequent activation of NADPH oxidase. Inhibition of the four targets could potentially dampen ROS overproduction which is associated with the initiation and/or progression of certain cancers [11,49,50,51,52]. To our knowledge, the Keap1-, MPO-, XO-, and p47phox-binding activity of the two peptides have not been previously reported. Considering the in silico evidence presented here, future investigations of the effectiveness of LYSPH and PSYLNTPLL in modulating the four targets, thus repressing ROS production and even cancer initiation and/or progression are warranted.
The five anticancer peptides predicted from our virtual library were also screened for cell-penetrating potential, blood-brain barrier penetrating potential, plasma half-life, and tolerance to in silico GI digestion, which can shed light on their potential as anticancer agents. Among the five peptides, LYSPH and PSYLNTPLL had the top two best plasma half-life (Table 5). The two peptides predictably had cell- and blood-brain barrier penetrating potential, although both were susceptible to GI digestion. The predicted cell-penetrating potential of the two peptides supports their potential in entering body cells and binding to the four intracellular targets: MPO, XO, Keap1, and p47phox, modulating the functions of the four proteins. The predicted ability of LYSPH and PSYLNTPLL to cross the blood-brain barrier may also facilitate their development as brain-tumor targeting peptides. Plasma half-life and tolerance to in silico GI digestion are related to the bioavailability of a peptide. Our results suggests that the two peptides were similar in their level of susceptibility to plasma peptidases, thus not differing much in their stability during systemic circulation. When compared with other natural anticancer peptides, such as KENPVLSLVNGMF identified from the giant barrel sponge Xestospongia testudinaria (half-life of 3.2 h in human serum in vitro) [53], the half-life of LYSPH and PSYLNTPLL was relatively short (about 14 min). Meanwhile, LYSPH and PSYLNTPLL were similarly susceptible to degradation by GI proteases. So, poor stability in blood and susceptibility to GI digestion is a key potential weakness of the two peptides, despite their ability to target MPO, XO, Keap1, and p47phox, as well as predictably having cell- and blood-brain barrier penetrating potential. The stability issue may limit the potential effectiveness of LYSPH and PSYLNTPLL as anticancer agents in the body, whether introduced into the body through oral or non-oral routes. To enhance the in vivo bioavailability of LYSPH and PSYLNTPLL, structural modifications that could improve their resistance to plasma and GI peptidases, such as cyclization of peptides [54] could be considered in future research. Moreover, the application of innovative technology such as mucoadhesive nanoparticles [55] may also be explored for oral delivery of the peptides with reduced risk of GI degradation and enhanced bioavailability.
Based on computational alanine scanning, Tyr played the most significant role in the binding and stabilizing of peptide-protein complexes for Keap1, MPO, XO, and p47phox. This can be observed from the drastically elevated ∆∆G values after the substitution of Tyr to Ala in both LYSPH and PSYLNTPLL (Table 6). This suggests that the hydrophobic interactions between Tyr and the residues of Keap1 (Tyr572), of XO (Glu879, Thr1010, Phe1013), of MPO (Phe99, Glu102, Phe146, Leu415, Leu420) and of p47phox (Gly192, Asp261, Gly262, Met278) (Figure 3, Figure 4, Figure 5 and Figure 6) are critical to the formation of stable peptide-protein complexes. In line with our findings, Wu and co-workers [56] found that the only Tyr-containing peptide in their study had the highest XO inhibitory activity; Tyr in the peptide also interacted hydrophobically with Phe1013 of XO. On the other hand, Ala substitution of His in LYSPH also led to the second largest increase in ∆∆G by 14.2670 kJ/mol (Table 6) when the LYSPH-XO complex was analyzed. By contrast, Ala substitution of His in LYSPH led to only a minor increase in ∆∆G of the LYSPH-Keap1 and LYSPH-MPO complexes. A possible explanation is that the His residue of LYSPH could bind to more key residues in XO (Glu802, Phe914, Phe1009, and Leu1014) (Figure 5c). By contrast, the His residue of LYSPH interacted with only one key residue (Arg380) in Keap1 (Figure 3c) and with none in MPO (Figure 4c). Our analysis suggests that future research that considers re-designing LYSPH and PSYLNTPLL for enhanced interactions with Keap1, MPO, and p47phox should avoid replacing or removing the Tyr residue. For stable binding to XO, the Tyr and His residues of LYSPH both should not be replaced or removed.
Molecular dynamics (MD) is a simulation technique which applied to derive the statements about the structural, dynamical, and thermodynamic properties of a molecular system [57]. The approach is able to observe minor conformational changes corresponds to the residue side chains which affect the binding site of a protein and ligand complementarity [57]. In the current study MD was applied to observe the dynamic level stability of each peptide ligand against the targeted proteins, as the peptides can functions either as the receptor inhibitors [58,59] or as the mediator such as the peptide mediated interactions in cell signaling [60].
The MD simulations results in Figure 7a–d determines the protein-ligand complexes stability during the 50 ns duration. In Figure 7a, the all-atom averaged root mean square deviation (RMSD) value for protein target Keap1, XO and MPO (chain A and B) in the complex were shown to be low at 2.14 ± 0.11 Å, 3.15 ± 0.33 Å, 2.59 ± 0.18 Å and 3.17 ± 0.36Å, respectively. In comparison, Figure 7b shows the all-atom averaged RMSD value of ligand LYSPH docked on each Keap1, XO and MPO were 1.73 ± 0.26 Å, 3.07 ± 0.33 Å and 3.34 ± 0.47 Å, respectively. This shows that RMSD values of both the docked proteins and the ligands are below the allowed limit [61], confirming that the protein-ligand complexes are stable over time. The plotted RMSD graphs also shows that receptor p47phox took longer time to reach complex stability compared to the other docked proteins with the averaged RMSD value of 5.14 ± 0.13 Å, while all-atom averaged RMSD for its ligand, PSYLNTPLL was 4.36 ± 0.45 Å. The high p47phox RMSD value was contributed by the flexibility of the N- and also C-terminal residues of the protein which reached up to 6.00 Å due to the loop structure of both terminals, visible by the root mean square fluctuations RMSF plot (Figure S1). The ligand interacted residues, however, were not affected and gave relatively low fluctuations during the 50 ns duration. In addition, the RMSD of PSYLNTPLL was also similarly low with LYSPH docked on other protein target (Figure 7b) during the 50 ns duration.
The dynamic intermolecular hydrogen bonds formed between the docked peptide and receptor protein were summarized in Figure 7c. The figure shows that highest number of intermolecular hydrogen bonds formed was in between LYSPH-XO (ave: 7), followed by PSYLNTPLL-p47phox (ave: 5) and LYSPH-Keap1 (ave: 3). MPO protein was consists of chain A and chain B domain, where chain A formed only one intermolecular hydrogen bond with the ligand in average while chain B has the average of three intermolecular hydrogen bonds formed within the 50 ns duration. The polar group of the XO hot-spot region, and LYSPH peptide both contributed to the higher number of hydrogen bonds formed making the complex more stable [62]. Higher surface of interactions between PSYLNTPLL-p47phox due to the longer sequence of the peptide had stabilized its docking on the active site of p47phox [63]. The intermolecular hydrogen bonds formed in the complex had also correlated with the distance formed between each ligand and protein, as summarized in Figure 7d. The result shows that all complexes were tightly packed with the average protein-ligand distance of 1.43 Å–1.82 Å, except for LYSPH and the chain A of MPO which varied from 1.43 Å up to 5.31 Å. This was contributed by the binding site of the ligand which located closer to the chain B of MPO. Overall, the duration 50 ns were shown to be sufficient to evaluate the stability of protein-ligand complex formation where most of the ligand tends to reach conformational stability after 3 ns. The RMSF and radius of gyration (Rg) plots of each complex are available in Figures S1 and S2.

3. Materials and Methods

3.1. Compilation of a Virtual Library of Seed-Derived Antioxidant Peptides

Seed-derived antioxidant peptides were compiled from the publications in the Scopus database by using the search words listed in Table S2 (Accessed: 5–7 October 2021). A total of 63 publications were carefully examined to find antioxidant peptides identified from different seed sources. In addition, peptides were compiled from the PlantPepDB database (http://14.139.61.8/PlantPepDB/index.php) [21] by using “Simple Search”, searching “Antioxidant” and selecting “Peptide Activity” as search field (Accessed: 5–6 October 2021). Following the exclusion of redundant sequences, the resulting collection of seed-derived antioxidant peptide sequences was used in subsequent screening and molecular modelling analyses, as depicted in Figure 8.

3.2. Virtual Screening for Anticancer Potential, Toxicity, Allergenicity and Hemolyticity

Anticancer potential was predicted by using AntiCP 2.0 (https://webs.iiitd.edu.in/raghava/anticp2/index.html) [22] with the default SVM threshold of 0.45. Seed-derived antioxidant peptide sequences that were predicted as anticancer peptides by both Model 1 and Model 2 in the AntiCP 2.0 tool were noted. Toxicity was predicted by using ToxinPred (https://webs.iiitd.edu.in/raghava/toxinpred/index.html) [33] with the SVM threshold of 0.0, by using two methods: (i) SVM (Swiss-Prot) + Motif, and (ii) SVM (TrEMBL) + Motif. Only peptide sequences that were predicted as non-toxic by both of the aforementioned methods are regarded as non-toxic. Allergenicity was predicted by using AllerTOP v. 2.0 (https://www.ddg-pharmfac.net/AllerTOP/index.html) [32]. Hemolyticity was predicted by using HemoPI (https://webs.iiitd.edu.in/raghava/hemopi/index.php) [31] with the SVM + Motif (HemoPI-1) method. The aforementioned tools were accessed on 8–9 October 2021. Only seed-derived antioxidant peptides of 4–50 residues were screened as the peptides outside this range cannot be analyzed by AntiCP 2.0. The 2D structures of selected anticancer peptides were drawn by using the ACD/ChemSketch freeware (ACD/ChemSketch, version 2019.2.1, Advanced Chemistry Development, Inc., Toronto, ON, Canada, www.acdlabs.com, 2019). The molecular masses of peptides were calculated by using PepDraw (https://pepdraw.com/) (Accessed: 9 October 2021).

3.3. Molecular Docking Analysis

The 3D structures of peptides predicted to be anticancer, non-toxic, non-allergenic, and non-hemolytic were constructed by using PEP-FOLD 3 (https://bioserv.rpbs.univ-paris-diderot.fr/services/PEP-FOLD3/) [64,65,66] (Accessed: 9 October 2021). Two hundred simulations were run and the resulting models were sorted by the sOPEP method. The best output model of each peptide was downloaded and used in molecular docking.
The crystal structures of human Keap1 complexed with 16-mer peptide of Nrf2 (PDB ID: 2FLU) [39], human MPO complexed with 7GD (PDB ID: 6WYD) [41], bovine XO complexed with quercetin (PDB ID: 3NVY) [12], and p47phox complexed with p22phox-derived proline-rich peptide (PDB ID: 1WLP) [14] were downloaded from the RCSB Protein Data Bank (https://www.rcsb.org/) [67] (Accessed: 9 October 2021). The separation of proteins and ligands from the crystals were performed by using BIOVIA Discovery Studio Visualizer (BIOVIA, Dassault Systèmes, BIOVIA Discovery Studio Visualizer, Version 20.1.0.192, San Diego: Dassault Systèmes, CA, USA, 2020).
For Keap1, MPO, and XO, the proteins were prepared as receptors in the PDBQT format after deleting water, adding polar hydrogen, and adding Kollman charges by using the AutoDock Tools 1.5.6 [25]. The prepared receptors were used for molecular docking on Webina 1.0.3 (https://durrantlab.pitt.edu/webina/) [24]. The co-crystalized ligands and the selected seed-derived peptides to be docked to the proteins were prepared as ligands and saved in the PDBQT format by using AutoDock Tools 1.5.6. Redocking of the co-crystalized ligands to Keap1, MPO and XO was performed using Webina 1.0.3 and RMSD was predicted by using LigRMSD v1.0 (https://ligrmsd.appsbio.utalca.cl/) [68]. For Keap1, redocking was performed by using tetrapeptide ETGE, the key motif of the 16-mer peptide of Nrf2, as recommended previously [40]. The coordinates of box center and box size used in molecular docking on Webina 1.0.3 were tabulated in Table S3. Molecular docking between p47phox and the peptide ligands was performed by using the HPEPDOCK Server (http://huanglab.phys.hust.edu.cn/hpepdock/) [69,70,71,72,73]. Redocking of the co-crystalized p22phox-derived proline-rich peptide (GPLGSKQPPSNPPPRPPAEARKKPS) to p47phox was also performed on HPEPDOCK and RMSD was predicted by using LigRMSD v1.0. Webina 1.0.3 and HPEPDOCK were accessed between 9 and 12 October 2021. Intermolecular interactions between proteins and peptides in selected docked models were analyzed and 2D interaction diagrams were generated by using LigPlot+ v.2.2.4 [74,75].

3.4. Prediction of Cell-Penetrating Potential, Blood-Brain Barrier Penetrating Potential, Plasma Half-Life, and Tolerance to In Silico GI Digestion

Cell-penetrating potential was predicted by using MLCPP (http://www.thegleelab.org/MLCPP/MLCPP.html) [23]. Blood-brain barrier penetrating potential was predicted by using B3Pred (https://webs.iiitd.edu.in/raghava/b3pred/index.html) [76] with Random- Forest (RF)-based prediction model and RF probability threshold of 0.1. Plasma half-life was predicted by using PlifePred (https://webs.iiitd.edu.in/raghava/plifepred/index.php) [77]. Tolerance to in silico GI digestion was predicted using the “enzyme(s) action” tool on BIOPEP-UWM (http://www.uwm.edu.pl/biochemia/index.php/en/biopep) [78] as previously reported [29]. The aforementioned tools were accessed on 18 October 2021.

3.5. Computational Alanine Scanning Mutagenesis

To assess the energetic contribution of individual residues in the selected seed-derived antioxidant peptides in conferring stability of binding to Keap1, MPO, XO, and p47phox, computational alanine scanning was performed by using BUDE Alanine Scan (https://pragmaticproteindesign.bio.ed.ac.uk/balas/) [79,80] (Accessed: 18–19 October 2021) as previously reported [27].

3.6. Molecular Dynamics Simulation

For a comprehensive analysis of the biomolecular dynamics, molecular dynamics (MD) simulation has evolved as the most powerful technique [26]. The detailed MD simulations of the complexes were conducted in GROMACS 2020 using the GROMOS96 54a7 force field [81]. The 54a7 force field was shown to improve the stability of α-helical structures in proteins and widely used in peptide simulations [82]. Molecular dynamics simulation was performed on each peptide ligand and protein complex of LYSPH-Keap1, LYSPH-MPO, LYSPH-XO, and PSYLNTPLL-p47phox for 50 ns duration. In the MD, each complex was solvated in a cubic box with the distance of 1.2 nm between the complex and each side of the solvated box [83]. Sodium and chloride ions were added to neutralize the total charge of the system. The complex then was energy-minimized using the steepest descent algorithm [84]. The simulation condition was set at the room temperature (300 K) and the atmospheric pressure (1 bar) to closely mimic the general experiment conditions. The NVT thermal equilibration was carried out with a constrained structure and a velocity rescale thermostat specific to GROMACS, followed by NPT pressure equilibration was applied with the same velocity-rescale temperature coupling in addition to the Parrinello−Rahman pressure coupling [85]. The fully temperature and pressure equilibrated system was then used as the initial configuration for the MD production dynamic analysis. All simulations were conducted using a 2 fs time step [86]. The results were then analyzed using GROMACS functions such as RMSD and RMSF, while the formation of hydrogen bonds between each peptide and target proteins were analyzed using GROMACS “gmx_hbond” functions. Additionally, the distance between each protein and its ligand peptide was measured using the “gmx_pairdist” function.

4. Conclusions

Our computational study narrowed down the 677 peptides in the virtual library to five candidates predicted to have anticancer potential, in addition to non-toxicity, non-allergenicity and non-hemolyticity. Structure-based virtual screening found that LYSPH was the most promising peptide in targeting MPO, XO, and Keap1. On the other hand, PSYLNTPLL was the candidate that interacted most stably with p47phox. LYSPH and PSYLNTPLL were predicted to have cell- and blood-brain barrier penetrating potential. Taken together, LYSPH and PSYLNTPLL are two potential candidates of anticancer peptides that deserve more in-depth explorations, particularly wet-lab experimental validations, in future.

Supplementary Materials

Table S1: Seed-derived antioxidant peptides compiled from Scopus and PlantPepDB databases; Table S2: Search words used in Scopus to compile seed-derived antioxidant peptides; Table S3: Coordinates of box center and box size for different targets in molecular docking, and RMSD values; Table S4: Docking scores for peptides that were experimentally demonstrated to inhibit p47phox-p22phox interaction and NADPH oxidase, in comparison with p22phox; Figure S1: RMSF plots for (a) PSYLNTPLL-p47phox, (b) LYSPH-Keap1, (c) LYSPH-XO, and (d) LYSPH-MPO; Figure S2: Gyration (Rg) plots of each complex.

Author Contributions

Conceptualization, T.-T.C. and F.-C.W.; methodology, T.-T.C. and F.-C.W.; software, T.-T.C. and M.Z.S.; validation, T.-T.C., F.-C.W., and M.Z.S.; formal analysis, C.C.-C.W., M.Z.S., and J.-A.K.; investigation, T.-T.C., M.Z.S., C.C.-C.W., and J.-A.K.; resources, T.-T.C., F.-C.W.; data curation, C.C.-C.W., J.-A.K., T.-T.C., and F.-C.W.; writing—original draft preparation, T.-T.C., J.-A.K., and M.Z.S.; writing—review and editing, T.-T.C., J.-A.K., F.-C.W., C.C.-C.W., and M.Z.S.; visualization, C.C.-C.W., J.-A.K., F.-C.W., and M.Z.S.; supervision, T.-T.C. and F.-C.W.; project administration, T.-T.C. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

The data presented in this study are available on request from the corresponding author.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Jakubczyk, A.; Karaś, M.; Rybczyńska-Tkaczyk, K.; Zielińska, E.; Zieliński, D. Current trends of bioactive peptides—New sources and therapeutic effect. Foods 2020, 9, 846. [Google Scholar] [CrossRef]
  2. Wong, F.-C.; Xiao, J.; Wang, S.; Ee, K.-Y.; Chai, T.-T. Advances on the antioxidant peptides from edible plant sources. Trends Food Sci. Technol. 2020, 99, 44–57. [Google Scholar] [CrossRef]
  3. Apostolopoulos, V.; Bojarska, J.; Chai, T.-T.; Elnagdy, S.; Kaczmarek, K.; Matsoukas, J.; New, R.; Parang, K.; Lopez, O.P.; Parhiz, H.; et al. A global review on short peptides: Frontiers and perspectives. Molecules 2021, 26, 430. [Google Scholar] [CrossRef]
  4. Chai, T.-T.; Tan, Y.-N.; Ee, K.-Y.; Xiao, J.; Wong, F.-C. Seeds, fermented foods, and agricultural by-products as sources of plant-derived antibacterial peptides. Crit. Rev. Food Sci. Nutr. 2019, 59, S162–S177. [Google Scholar] [CrossRef]
  5. Díaz-Gómez, J.L.; Castorena-Torres, F.; Preciado-Ortiz, R.E.; García-Lara, S. Anti-cancer activity of maize bioactive peptides. Front. Chem. 2017, 5, 44. [Google Scholar] [CrossRef]
  6. Ramkisson, S.; Dwarka, D.; Venter, S.; Mellem, J.J. In vitro anticancer and antioxidant potential of Amaranthus cruentus protein and its hydrolysates. Food Sci. Technol. 2020, 40, 634–639. [Google Scholar] [CrossRef]
  7. Mika, D.; Guruvayoorappan, C. Myeloperoxidase: The yin and yang in tumour progression. J. Exp. Ther. Oncol. 2011, 9, 93–100. [Google Scholar]
  8. Kargapolova, Y.; Geißen, S.; Zheng, R.; Baldus, S.; Winkels, H.; Adam, M. The enzymatic and non-enzymatic function of myeloperoxidase (MPO) in inflammatory communication. Antioxidants 2021, 10, 562. [Google Scholar] [CrossRef]
  9. Wang, Y.; Qi, H.; Liu, Y.; Duan, C.; Liu, X.; Xia, T.; Chen, D.; Piao, H.-L.; Liu, H.-X. The double-edged roles of ROS in cancer prevention and therapy. Theranostics 2021, 11, 4839–4857. [Google Scholar] [CrossRef]
  10. Kazakov, Y.; Tarasov, A.; Alyoshina, L.; Brainina, K. Interplay between antioxidant activity, health and disease. Biointerface Res. Appl. Chem. 2020, 10, 4893–4901. [Google Scholar] [CrossRef]
  11. Oh, S.-H.; Choi, S.-Y.; Choi, H.-J.; Ryu, H.-M.; Kim, Y.-J.; Jung, H.-Y.; Cho, J.-H.; Kim, C.-D.; Park, S.-H.; Kwon, T.-H.; et al. The emerging role of xanthine oxidase inhibition for suppression of breast cancer cell migration and metastasis associated with hypercholesterolemia. FASEB J. 2019, 33, 7301–7314. [Google Scholar] [CrossRef]
  12. Cao, H.; Pauff, J.M.; Hille, R. X-ray crystal structure of a xanthine oxidase complex with the flavonoid inhibitor quercetin. J. Nat. Prod. 2014, 77, 1693–1699. [Google Scholar] [CrossRef]
  13. Vermot, A.; Petit-Härtlein, I.; Smith, S.M.E.; Fieschi, F. NADPH oxidases (NOX): An overview from discovery, molecular mechanisms to physiology and pathology. Antioxidants 2021, 10, 890. [Google Scholar] [CrossRef] [PubMed]
  14. Ogura, K.; Nobuhisa, I.; Yuzawa, S.; Takeya, R.; Torikai, S.; Saikawa, K.; Sumimoto, H.; Inagaki, F. NMR solution structure of the tandem Src homology 3 domains of p47phox complexed with a p22phox-derived proline-rich peptide. J. Biol. Chem. 2006, 281, 3660–3668. [Google Scholar] [CrossRef] [Green Version]
  15. Macías Pérez, M.E.; Hernández Rodríguez, M.; Cabrera Pérez, L.C.; Fragoso-Vázquez, M.J.; Correa-Basurto, J.; Padilla-Martínez, I.I.; Méndez Luna, D.; Mera Jiménez, E.; Flores Sandoval, C.; Tamay Cach, F.; et al. Aromatic regions govern the recognition of NADPH oxidase inhibitors as diapocynin and its analogues. Arch. Der Pharm. 2017, 350, 1700041. [Google Scholar] [CrossRef]
  16. Konaté, M.M.; Antony, S.; Doroshow, J.H. Inhibiting the activity of NADPH oxidase in cancer. Antioxid Redox Signal. 2020, 33, 435–454. [Google Scholar] [CrossRef] [PubMed]
  17. Deshmukh, P.; Unni, S.; Krishnappa, G.; Padmanabhan, B. The Keap1–Nrf2 pathway: Promising therapeutic target to counteract ROS-mediated damage in cancers and neurodegenerative diseases. Biophys. Rev. 2017, 9, 41–56. [Google Scholar] [CrossRef]
  18. Tascioglu Aliyev, A.; Panieri, E.; Stepanić, V.; Gurer-Orhan, H.; Saso, L. Involvement of NRF2 in breast cancer and possible therapeutical role of polyphenols and melatonin. Molecules 2021, 26, 1853. [Google Scholar] [CrossRef] [PubMed]
  19. Wen, C.; Zhang, J.; Zhang, H.; Duan, Y.; Ma, H. Study on the structure–activity relationship of watermelon seed antioxidant peptides by using molecular simulations. Food Chem. 2021, 364, 130432. [Google Scholar] [CrossRef]
  20. Koh, J.-A.; Ong, J.-H.; Abd Manan, F.; Ee, K.-Y.; Wong, F.-C.; Chai, T.-T. Discovery of bifunctional anti-DPP-IV and anti-ACE peptides from housefly larval proteins after in silico gastrointestinal digestion. Biointerface Res. Appl. Chem. 2022, 12, 4929–4944. [Google Scholar] [CrossRef]
  21. Das, D.; Jaiswal, M.; Khan, F.N.; Ahamad, S.; Kumar, S. PlantPepDB: A manually curated plant peptide database. Sci. Rep. 2020, 10, 2194. [Google Scholar] [CrossRef] [Green Version]
  22. Agrawal, P.; Bhagat, D.; Mahalwal, M.; Sharma, N.; Raghava, G.P.S. AntiCP 2.0: An updated model for predicting anticancer peptides. Brief. Bioinform. 2020, 22, bbaa153. [Google Scholar] [CrossRef] [PubMed]
  23. Manavalan, B.; Subramaniyam, S.; Shin, T.H.; Kim, M.O.; Lee, G. Machine-learning-based prediction of cell-penetrating peptides and their uptake efficiency with improved accuracy. J. Proteome Res. 2018, 17, 2715–2726. [Google Scholar] [CrossRef]
  24. Kochnev, Y.; Hellemann, E.; Cassidy, K.C.; Durrant, J.D. Webina: An open-source library and web app that runs AutoDock Vina entirely in the web browser. Bioinformatics 2020, 36, 4513–4515. [Google Scholar] [CrossRef] [PubMed]
  25. Morris, G.M.; Ruth, H.; Lindstrom, W.; Sanner, M.F.; Belew, R.K.; Goodsell, D.S.; Olson, A.J. Software news and updates AutoDock4 and AutoDockTools4: Automated docking with selective receptor flexibility. J. Comput. Chem. 2009, 30, 2785–2791. [Google Scholar] [CrossRef] [Green Version]
  26. Shukla, R.; Tripathi, T. Molecular dynamics simulation of protein and protein-ligand complexes. In Computer-Aided Drug Design; Springer: Singapore, 2020; Chapter 7; pp. 133–161. [Google Scholar]
  27. Ong, J.-H.; Liang, C.-E.; Wong, W.-L.; Wong, F.-C.; Chai, T.-T. Multi-target anti-SARS-CoV-2 peptides from mealworm proteins: An in silico study. Malays. J. Biochem. Mol. Biol. 2021, 24, 83–91. [Google Scholar]
  28. Mohana, D.S.; Chai, T.-T.; Wong, F.-C. Antioxidant and protein protection potentials of fennel seed-derived protein hydrolysates and peptides. Mod. Food Sci. Technol. 2019, 35, 22–29. [Google Scholar] [CrossRef]
  29. Chai, T.-T.; Soo, Z.-Y.; Hsu, K.-C.; Li, J.-C.; Abd Manan, F.; Wong, F.-C. Antioxidant activity of semen cassiae protein hydrolysate: Thermal and gastrointestinal stability, peptide identification, and in silico analysis. Mod. Food Sci. Technol. 2019, 35, 38–48. [Google Scholar]
  30. Chai, T.-T.; Xiao, J.; Mohana Dass, S.; Teoh, J.-Y.; Ee, K.-Y.; Ng, W.-J.; Wong, F.-C. Identification of antioxidant peptides derived from tropical jackfruit seed and investigation of the stability profiles. Food Chem. 2021, 340, 127876. [Google Scholar] [CrossRef] [PubMed]
  31. Chaudhary, K.; Kumar, R.; Singh, S.; Tuknait, A.; Gautam, A.; Mathur, D.; Anand, P.; Varshney, G.C.; Raghava, G.P.S. A web server and mobile app for computing hemolytic potency of peptides. Sci. Rep. 2016, 6, 22843. [Google Scholar] [CrossRef]
  32. Dimitrov, I.; Bangov, I.; Flower, D.R.; Doytchinova, I. AllerTOP v.2—A server for in silico prediction of allergens. J. Mol. Model. 2014, 20, 2278. [Google Scholar] [CrossRef]
  33. Gupta, S.; Kapoor, P.; Chaudhary, K.; Gautam, A.; Kumar, R.; Open Source Drug Discovery, C.; Raghava, G.P.S. In silico approach for predicting toxicity of peptides and proteins. PLoS ONE 2013, 8, e73957. [Google Scholar] [CrossRef] [Green Version]
  34. Ji, D.; Udenigwe, C.C.; Agyei, D. Antioxidant peptides encrypted in flaxseed proteome: An in silico assessment. Food Sci. Hum. Wellness 2019, 8, 306–314. [Google Scholar] [CrossRef]
  35. Duan, X.; Zhang, M.; Chen, F. Prediction and analysis of antimicrobial peptides from rapeseed protein using in silico approach. J. Food Biochem. 2021, 45, e13598. [Google Scholar] [CrossRef]
  36. Taniguchi, M.; Aida, R.; Saito, K.; Kikura, T.; Ochiai, A.; Saitoh, E.; Tanaka, T. Identification and characterization of multifunctional cationic peptides from enzymatic hydrolysates of soybean proteins. J. Biosci. Bioeng. 2020, 129, 59–66. [Google Scholar] [CrossRef]
  37. Ghribi, A.M.; Sila, A.; Przybylski, R.; Nedjar-Arroume, N.; Makhlouf, I.; Blecker, C.; Attia, H.; Dhulster, P.; Bougatef, A.; Besbes, S. Purification and identification of novel antioxidant peptides from enzymatic hydrolysate of chickpea (Cicer arietinum L.) protein concentrate. J. Funct. Foods 2015, 12, 516–525. [Google Scholar] [CrossRef]
  38. Chai, T.-T.; Law, Y.-C.; Wong, F.-C.; Kim, S.-K. Enzyme-assisted discovery of antioxidant peptides from edible marine invertebrates: A review. Mar. Drugs 2017, 15, 42. [Google Scholar] [CrossRef] [PubMed]
  39. Lo, S.-C.; Li, X.; Henzl, M.T.; Beamer, L.J.; Hannink, M. Structure of the Keap1:Nrf2 interface provides mechanistic insight into Nrf2 signaling. EMBO J. 2006, 25, 3605–3617. [Google Scholar] [CrossRef] [Green Version]
  40. Li, L.; Liu, J.; Nie, S.; Ding, L.; Wang, L.; Liu, J.; Liu, W.; Zhang, T. Direct inhibition of Keap1–Nrf2 interaction by egg-derived peptides DKK and DDW revealed by molecular docking and fluorescence polarization. RSC Adv. 2017, 7, 34963–34971. [Google Scholar] [CrossRef] [Green Version]
  41. Shaw, S.A.; Vokits, B.P.; Dilger, A.K.; Viet, A.; Clark, C.G.; Abell, L.M.; Locke, G.A.; Duke, G.; Kopcho, L.M.; Dongre, A.; et al. Discovery and structure activity relationships of 7-benzyl triazolopyridines as stable, selective, and reversible inhibitors of myeloperoxidase. Bioorganic Med. Chem. 2020, 28, 115723. [Google Scholar] [CrossRef]
  42. Maiocchi, S.L.; Ku, J.; Thai, T.; Chan, E.; Rees, M.D.; Thomas, S.R. Myeloperoxidase: A versatile mediator of endothelial dysfunction and therapeutic target during cardiovascular disease. Pharmacol. Ther. 2021, 221, 107711. [Google Scholar] [CrossRef]
  43. Kang, N.; Kim, E.-A.; Kim, J.; Lee, S.-H.; Heo, S.-J. Identifying potential antioxidant properties from the viscera of sea snails (Turbo cornutus). Mar. Drugs 2021, 19, 567. [Google Scholar] [CrossRef] [PubMed]
  44. Davey, C.A.; Fenna, R.E. 2.3 Å resolution x-ray crystal structure of the bisubstrate analogue inhibitor salicylhydroxamic acid bound to human myeloperoxidase:  A model for a prereaction complex with hydrogen peroxide. Biochemistry 1996, 35, 10967–10973. [Google Scholar] [CrossRef] [PubMed]
  45. Hu, X.; Zhou, Y.; Zhou, S.; Chen, S.; Wu, Y.; Li, L.; Yang, X. Purification and identification of novel xanthine oxidase inhibitory peptides derived from round scad (Decapterus maruadsi) protein hydrolysates. Mar. Drugs 2021, 19, 538. [Google Scholar] [CrossRef]
  46. Yu, Z.; Kan, R.; Wu, S.; Guo, H.; Zhao, W.; Ding, L.; Zheng, F.; Liu, J. Xanthine oxidase inhibitory peptides derived from tuna protein: Virtual screening, inhibitory activity, and molecular mechanisms. J. Sci. Food Agric. 2021, 101, 1349–1354. [Google Scholar] [CrossRef]
  47. Nakanishi, A.; Imajoh-Ohmi, S.; Fujinawa, T.; Kikuchi, H.; Kanegasaki, S. Direct evidence for interaction between COOH-terminal regions of cytochrome b558 subunits and cytosolic 47-kDa protein during activation of an O(2−)-generating system in neutrophils. J. Biol. Chem. 1992, 267, 19072–19074. [Google Scholar] [CrossRef]
  48. Huang, J.; Kleinberg, M.E. Activation of the phagocyte NADPH oxidase protein p47phox: Phosphorylation controls SH3 domain-dependent binding to p22phox *. J. Biol. Chem. 1999, 274, 19731–19737. [Google Scholar] [CrossRef] [Green Version]
  49. Arora, R.; Sawney, S.; Saini, V.; Steffi, C.; Tiwari, M.; Saluja, D. Esculetin induces antiproliferative and apoptotic response in pancreatic cancer cells by directly binding to KEAP1. Mol. Cancer 2016, 15, 64. [Google Scholar] [CrossRef] [Green Version]
  50. Zhang, H.; Jing, X.; Shi, Y.; Xu, H.; Du, J.; Guan, T.; Weihrauch, D.; Jones, D.W.; Wang, W.; Gourlay, D.; et al. N-acetyl lysyltyrosylcysteine amide inhibits myeloperoxidase, a novel tripeptide inhibitor1[S]. J. Lipid Res. 2013, 54, 3016–3029. [Google Scholar] [CrossRef] [Green Version]
  51. Wang, Z.; Tang, T.; Wang, S.; Cai, T.; Tao, H.; Zhang, Q.; Qi, S.; Qi, Z. Aloin inhibits the proliferation and migration of gastric cancer cells by regulating NOX2-ROS-mediated pro-survival signal pathways. Drug Des. Devel. 2020, 14, 145–155. [Google Scholar] [CrossRef] [Green Version]
  52. Rymaszewski, A.L.; Tate, E.; Yimbesalu, J.P.; Gelman, A.E.; Jarzembowski, J.A.; Zhang, H.; Pritchard, K.A., Jr.; Vikis, H.G. The role of neutrophil myeloperoxidase in models of lung tumor development. Cancers 2014, 6, 1111–1127. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  53. Quah, Y.; Mohd Ismail, N.I.; Ooi, J.L.S.; Affendi, Y.A.; Abd Manan, F.; Wong, F.-C.; Chai, T.-T. Identification of novel cytotoxic peptide KENPVLSLVNGMF from marine sponge Xestospongia testudinaria, with characterization of stability in human serum. Int. J. Pept. Res. Ther. 2018, 24, 189–199. [Google Scholar] [CrossRef]
  54. Amigo, L.; Hernández-Ledesma, B. Current evidence on the bioavailability of food bioactive peptides. Molecules 2020, 25, 4479. [Google Scholar] [CrossRef]
  55. Sharifi, S.; Samani, A.A.; Ahmadian, E.; Eftekhari, A.; Derakhshankhah, H.; Jafari, S.; Mokhtarpour, M.; Vahed, S.Z.; Salatin, S.; Dizaj, S.M. Oral delivery of proteins and peptides by mucoadhesive nanoparticles. Biointerface Res. Appl. Chem. 2019, 9, 3849–3852. [Google Scholar] [CrossRef]
  56. Wu, Y.; He, H.; Hou, T. Purification, identification, and computational analysis of xanthine oxidase inhibitory peptides from kidney bean. J. Food Sci. 2021, 86, 1081–1088. [Google Scholar] [CrossRef]
  57. Salo-Ahen, O.M.H.; Alanko, I.; Bhadane, R.; Bonvin, A.M.J.J.; Honorato, R.V.; Hossain, S.; Juffer, A.H.; Kabedev, A.; Lahtela-Kakkonen, M.; Larsen, A.S.; et al. Molecular dynamics simulations in drug discovery and pharmaceutical development. Processes 2021, 9, 71. [Google Scholar] [CrossRef]
  58. Mangold, M.; Gütschow, M.; Stirnberg, M. A short peptide inhibitor as an activity-based probe for matriptase-2. Pharmaceuticals 2018, 11, 49. [Google Scholar] [CrossRef] [Green Version]
  59. Baig, M.S.; Alagumuthu, M.; Rajpoot, S.; Saqib, U. Identification of a potential peptide inhibitor of SARS-CoV-2 targeting its entry into the host cells. Drugs R D 2020, 20, 161–169. [Google Scholar] [CrossRef]
  60. Zhou, P.; Miao, Q.; Yan, F.; Li, Z.; Jiang, Q.; Wen, L.; Meng, Y. Is protein context responsible for peptide-mediated interactions? Mol. Omics 2019, 15, 280–295. [Google Scholar] [CrossRef] [PubMed]
  61. Azizian, H.; Nabati, F.; Sharifi, A.; Siavoshi, F.; Mahdavi, M.; Amanlou, M. Large-scale virtual screening for the identification of new Helicobacter pylori urease inhibitor scaffolds. J. Mol. Model. 2012, 18, 2917–2927. [Google Scholar] [CrossRef] [PubMed]
  62. Keskin, O.; Ma, B.; Nussinov, R. Hot regions in protein–protein interactions: The organization and contribution of structurally conserved hot spot residues. J. Mol. Biol. 2005, 345, 1281–1294. [Google Scholar] [CrossRef]
  63. London, N.; Movshovitz-Attias, D.; Schueler-Furman, O. The structural basis of peptide-protein binding strategies. Structure 2010, 18, 188–199. [Google Scholar] [CrossRef] [Green Version]
  64. Thévenet, P.; Shen, Y.; Maupetit, J.; Guyon, F.; Derreumaux, P.; Tufféry, P. PEP-FOLD: An updated de novo structure prediction server for both linear and disulfide bonded cyclic peptides. Nucleic Acids Res. 2012, 40, W288–W293. [Google Scholar] [CrossRef] [Green Version]
  65. Lamiable, A.; Thévenet, P.; Rey, J.; Vavrusa, M.; Derreumaux, P.; Tufféry, P. PEP-FOLD3: Faster de novo structure prediction for linear peptides in solution and in complex. Nucleic Acids Res. 2016, 44, W449–W454. [Google Scholar] [CrossRef] [Green Version]
  66. Shen, Y.; Maupetit, J.; Derreumaux, P.; Tufféry, P. Improved PEP-FOLD approach for peptide and miniprotein structure prediction. J. Chem. Theory Comput. 2014, 10, 4745–4758. [Google Scholar] [CrossRef]
  67. Burley, S.K.; Berman, H.M.; Bhikadiya, C.; Bi, C.; Chen, L.; Di Costanzo, L.; Christie, C.; Dalenberg, K.; Duarte, J.M.; Dutta, S.; et al. RCSB Protein Data Bank: Biological macromolecular structures enabling research and education in fundamental biology, biomedicine, biotechnology and energy. Nucleic Acids Res. 2018, 47, D464–D474. [Google Scholar] [CrossRef] [Green Version]
  68. Velázquez-Libera, J.L.; Durán-Verdugo, F.; Valdés-Jiménez, A.; Valdés-Jiménez, A.; Núñez-Vivanco, G.; Caballero, J. LigRMSD: A web server for automatic structure matching and RMSD calculations among identical and similar compounds in protein-ligand docking. Bioinformatics 2020, 36, 2912–2914. [Google Scholar] [CrossRef]
  69. Zhou, P.; Jin, B.; Li, H.; Huang, S.Y. HPEPDOCK: A web server for blind peptide-protein docking based on a hierarchical algorithm. Nucleic Acids Res. 2018, 46, W443–W450. [Google Scholar] [CrossRef]
  70. Huang, S.-Y.; Zou, X. Ensemble docking of multiple protein structures: Considering protein structural variations in molecular docking. Proteins Struct. Funct. Bioinform. 2007, 66, 399–421. [Google Scholar] [CrossRef]
  71. Huang, S.-Y.; Zou, X. An iterative knowledge-based scoring function for protein–protein recognition. Proteins: Struct. Funct. Bioinform. 2008, 72, 557–579. [Google Scholar] [CrossRef]
  72. Yan, Y.; Zhang, D.; Huang, S.-Y. Efficient conformational ensemble generation of protein-bound peptides. J. Cheminform. 2017, 9, 59. [Google Scholar] [CrossRef]
  73. Zhou, P.; Li, B.; Yan, Y.; Jin, B.; Wang, L.; Huang, S.Y. Hierarchical flexible peptide docking by conformer generation and ensemble docking of peptides. J. Chem. Inf. Model. 2018, 58, 1292–1302. [Google Scholar] [CrossRef]
  74. Laskowski, R.A.; Swindells, M.B. LigPlot+: Multiple ligand-protein interaction diagrams for drug discovery. J. Chem. Inf. Model. 2011, 51, 2778–2786. [Google Scholar] [CrossRef]
  75. Wallace, A.C.; Laskowski, R.A.; Thornton, J.M. LIGPLOT: A program to generate schematic diagrams of protein-ligand interactions. Protein Eng. Des. Sel. 1995, 8, 127–134. [Google Scholar] [CrossRef] [PubMed]
  76. Kumar, V.; Patiyal, S.; Dhall, A.; Sharma, N.; Raghava, G.P.S. B3pred: A random-forest-based method for predicting and designing blood–brain barrier penetrating peptides. Pharmaceutics 2021, 13, 1237. [Google Scholar] [CrossRef]
  77. Mathur, D.; Singh, S.; Mehta, A.; Agrawal, P.; Raghava, G.P.S. In silico approaches for predicting the half-life of natural and modified peptides in blood. PLoS ONE 2018, 13, e0196829. [Google Scholar] [CrossRef] [Green Version]
  78. Minkiewicz, P.; Iwaniak, A.; Darewicz, M. BIOPEP-UWM database of bioactive peptides: Current opportunities. Int. J. Mol. Sci. 2019, 20, 5978. [Google Scholar] [CrossRef] [Green Version]
  79. Ibarra, A.A.; Bartlett, G.J.; Hegedüs, Z.; Dutt, S.; Hobor, F.; Horner, K.A.; Hetherington, K.; Spence, K.; Nelson, A.; Edwards, T.A.; et al. Predicting and experimentally validating hot-spot residues at protein–protein interfaces. ACS Chem. Biol. 2019, 14, 2252–2263. [Google Scholar] [CrossRef] [PubMed]
  80. Wood, C.W.; Ibarra, A.A.; Bartlett, G.J.; Wilson, A.J.; Woolfson, D.N.; Sessions, R.B. BAlaS: Fast, interactive and accessible computational alanine-scanning using BudeAlaScan. Bioinformatics 2020, 36, 2917–2919. [Google Scholar] [CrossRef]
  81. Schmid, N.; Eichenberger, A.P.; Choutko, A.; Riniker, S.; Winger, M.; Mark, A.E.; Van Gunsteren, W.F. Definition and testing of the GROMOS force-field versions 54A7 and 54B7. Eur. Biophys. J. 2011, 40, 843–856. [Google Scholar] [CrossRef]
  82. Huang, W.; Lin, Z.; van Gunsteren, W.F. Validation of the GROMOS 54A7 force field with respect to β-peptide folding. J. Chem. Theory Comput. 2011, 7, 1237–1243. [Google Scholar] [CrossRef]
  83. Leherte, L.; Vercauteren, D.P. Reduced point charge models of proteins: Effect of protein–water interactions in molecular dynamics simulations of ubiquitin systems. J. Phys. Chem. B 2017, 121, 9771–9784. [Google Scholar] [CrossRef] [PubMed]
  84. Nunes-Alves, A.; Ormersbach, F.; Wade, R.C. Prediction of the drug–target binding kinetics for flexible proteins by comparative binding energy analysis. J. Chem. Inf. Model. 2021, 61, 3708–3721. [Google Scholar] [CrossRef]
  85. Sabri, M.Z.; Hamid, A.A.A.; Hitam, S.M.S.; Rahim, M.Z.A. The assessment of three dimensional modelling design for single strand DNA aptamers for computational chemistry application. Biophys. Chem. 2020, 267, 106492. [Google Scholar] [CrossRef] [PubMed]
  86. Sabri, M.Z.; Abdul Hamid, A.A.; Sayed Hitam, S.M.; Abdul Rahim, M.Z. In silico screening of aptamers configuration against hepatitis b surface antigen. Adv. Bioinform. 2019, 2019, 6912914. [Google Scholar] [CrossRef] [Green Version]
Molecules 26 07396 g001 550
Figure 1. (a) The distribution of peptide length in the virtual library of seed-derived antioxidant peptides. (b) Outcome of the screening of 592 seed-derived antioxidant peptides of 4–50 residues for anticancer potential, toxicity, allergenicity, and hemolyticity. Numbers in the pie charts represent numbers of peptides with the predicted properties.
Figure 1. (a) The distribution of peptide length in the virtual library of seed-derived antioxidant peptides. (b) Outcome of the screening of 592 seed-derived antioxidant peptides of 4–50 residues for anticancer potential, toxicity, allergenicity, and hemolyticity. Numbers in the pie charts represent numbers of peptides with the predicted properties.
Molecules 26 07396 g001
Molecules 26 07396 g002 550
Figure 2. 2D structures of five selected seed-derived antioxidant peptides predicted to be anticancer, non-toxic, non-allergenic, and non-hemolytic.
Figure 2. 2D structures of five selected seed-derived antioxidant peptides predicted to be anticancer, non-toxic, non-allergenic, and non-hemolytic.
Molecules 26 07396 g002
Molecules 26 07396 g003 550
Figure 3. (a) Three-dimensional (3D) diagram of LYSPH; (b) 3D diagram of LYSPH-Keap1 docked model; (c) 2D LYSPH-Keap1 interaction diagram. In (a,b), LYSPH is displayed in a blue-stick style. In (b), Keap1 is displayed as red ribbon. In (c), green dashed lines and red spoked arcs represent hydrogen bonds and hydrophobic interactions, respectively. Residues of LYSPH are shown in purple bonds, whereas residues of Keap1 are shown in brown bonds and also represented by the red spoked arcs.
Figure 3. (a) Three-dimensional (3D) diagram of LYSPH; (b) 3D diagram of LYSPH-Keap1 docked model; (c) 2D LYSPH-Keap1 interaction diagram. In (a,b), LYSPH is displayed in a blue-stick style. In (b), Keap1 is displayed as red ribbon. In (c), green dashed lines and red spoked arcs represent hydrogen bonds and hydrophobic interactions, respectively. Residues of LYSPH are shown in purple bonds, whereas residues of Keap1 are shown in brown bonds and also represented by the red spoked arcs.
Molecules 26 07396 g003
Molecules 26 07396 g004 550
Figure 4. (a) 3D diagram of LYSPH; (b) 3D diagram of LYSPH-MPO docked model; (c) 2D LYSPH-MPO interaction diagram. In (a,b), LYSPH is displayed in a blue-stick style. In (b), MPO is displayed as red ribbon. In (c), red spoked arcs represent hydrophobic interactions. Residues of LYSPH are shown in purple bonds. Residues of MPO are represented by red spoked arcs.
Figure 4. (a) 3D diagram of LYSPH; (b) 3D diagram of LYSPH-MPO docked model; (c) 2D LYSPH-MPO interaction diagram. In (a,b), LYSPH is displayed in a blue-stick style. In (b), MPO is displayed as red ribbon. In (c), red spoked arcs represent hydrophobic interactions. Residues of LYSPH are shown in purple bonds. Residues of MPO are represented by red spoked arcs.
Molecules 26 07396 g004
Molecules 26 07396 g005 550
Figure 5. (a) 3D diagram of LYSPH; (b) 3D diagram of LYSPH-XO docked model; (c) 2D LYSPH-XO interaction diagram. In (a,b), LYSPH is displayed in a blue-stick style. In (b), XO is displayed as red ribbon. In (c), green dashed line and red spoked arcs represent hydrogen bond and hydrophobic interactions, respectively. Residues of LYSPH are shown in purple bonds, whereas residues of XO are shown in brown bonds and also represented by the red spoked arcs.
Figure 5. (a) 3D diagram of LYSPH; (b) 3D diagram of LYSPH-XO docked model; (c) 2D LYSPH-XO interaction diagram. In (a,b), LYSPH is displayed in a blue-stick style. In (b), XO is displayed as red ribbon. In (c), green dashed line and red spoked arcs represent hydrogen bond and hydrophobic interactions, respectively. Residues of LYSPH are shown in purple bonds, whereas residues of XO are shown in brown bonds and also represented by the red spoked arcs.
Molecules 26 07396 g005
Molecules 26 07396 g006 550
Figure 6. (a) 3D diagram of PSYLNTPLL; (b) 3D diagram of PSYLNTPLL-p47phox docked model; (c) 2D PSYLNTPLL-p47phox interaction diagram. In (a,b), PSYLNTPLL is displayed in a blue-stick style. In (b), p47phox is displayed as red ribbon. In (c), green dashed lines and red spoked arcs represent hydrogen bonds and hydrophobic interactions, respectively. Residues of PSYLNTPLL are shown in purple bonds, whereas residues of p47phox are shown in brown bonds and also represented by the red spoked arcs.
Figure 6. (a) 3D diagram of PSYLNTPLL; (b) 3D diagram of PSYLNTPLL-p47phox docked model; (c) 2D PSYLNTPLL-p47phox interaction diagram. In (a,b), PSYLNTPLL is displayed in a blue-stick style. In (b), p47phox is displayed as red ribbon. In (c), green dashed lines and red spoked arcs represent hydrogen bonds and hydrophobic interactions, respectively. Residues of PSYLNTPLL are shown in purple bonds, whereas residues of p47phox are shown in brown bonds and also represented by the red spoked arcs.
Molecules 26 07396 g006
Molecules 26 07396 g007 550
Figure 7. (a) all-atom RMSD value for proteins; (b) all-atom RMSD value for ligands; (c) number of intermolecular hydrogen bonds; and (d) minimum distance between the proteins and each peptide ligand for the 50 ns duration.
Figure 7. (a) all-atom RMSD value for proteins; (b) all-atom RMSD value for ligands; (c) number of intermolecular hydrogen bonds; and (d) minimum distance between the proteins and each peptide ligand for the 50 ns duration.
Molecules 26 07396 g007
Molecules 26 07396 g008 550
Figure 8. An overview of the computational approach used in this study.
Figure 8. An overview of the computational approach used in this study.
Molecules 26 07396 g008
Table
Table 1. Intermolecular interactions between seed-derived antioxidant peptides and Keap1.
Table 1. Intermolecular interactions between seed-derived antioxidant peptides and Keap1.
PeptideBinding Affinity (kcal/mol)Interaction with Keap1 a
Hydrogen BondHydrophobic InteractionSalt Bridge
LYSPH−7.6Arg380, Arg415(2), Arg483, Ser508Tyr334, Ser363, Arg380, Asn382, Arg415, Ser508, Tyr525, Gln530, Ala556, Tyr572, Phe577, Ser602-
LPHFNS−7.4Ser363, Arg380, Asn414, Arg415(2), Gln530Tyr334, Ser363, Arg380, Asn382, Asn414, Arg415, Arg483, Tyr525, Gln530, Ser555, Ala556, Tyr572, Phe577, Ser602-
AEHGSLH−7.2Tyr334(2), Arg380, Asn382(2), Arg415Tyr334, Arg336, Ser363, Arg380, Asn382, Ser383, Pro384, Arg415, Ile461, Arg483, Ser508, Tyr525, Gln530, Ala556, Tyr572, Ser602-
FGPEMEQ−7.0Ser363, Arg380(2), Asn382(2), Asn387, Asn414Tyr334, Ser363, Arg380, Asn382, Asn387, Asp389, Arg415, Gly433, Ile461, Ser555, Ala556, Tyr572, Phe577Arg380(2)
PSYLNTPLL−6.4Arg380(2), Asn382, Arg415, Arg483, Ser555, Tyr572Tyr334, Ser363, Gly364, Arg380, Asn382, Arg415, Arg483, Tyr525, Gln530, Ser555, Ala556, Tyr572, Phe577-
a Number in brackets indicates the number of hydrogen bonds or salt bridges formed with the same residue of Keap1. Keap1 residues that were reported to bind to ETGE (the key motif of Nrf2 peptide) [39] are marked in boldface type. Residues in the Keap1 binding pocket that were reported to contribute to stability of the Keap1:Nrf2 complex as evidenced by mutagenesis studies [39] are underlined.
Table
Table 2. Intermolecular interactions between seed-derived antioxidant peptides and MPO.
Table 2. Intermolecular interactions between seed-derived antioxidant peptides and MPO.
PeptideBinding Affinity (kcal/mol)Interaction with MPO
Hydrogen BondHydrophobic Interaction
LYSPH−7.1-Phe99, Glu102, Phe146, Pro220, Thr238, Arg239, Glu242, Phe366, Phe407, Leu415, Leu420, Hec606
LPHFNS−6.9-Phe99, Glu102, Glu116, Phe147, Pro220, Thr238, Arg239, Phe366, Phe407, Met411, Hec606
FGPEMEQ−6.8Glu102Phe99, Thr100, Glu102, Glu116, Pro145, Phe147, Leu216, Pro220, Thr238, Arg239, Glu242, Phe366, Phe407, Met411, Leu415, Arg424, Hec606
AEHGSLH−6.3-Phe99, Thr100, Glu102, Glu116, Pro145, Phe146, Phe147, Thr238, Arg239, Glu242, Phe366, Phe407, Val410, Met411, Leu415, Arg424, Hec606
PSYLNTPLL−3.0Thr100Phe99, Glu102, Glu116, Pro145, Phe147, Leu216, Pro220, Thr238, Arg239, Phe366, Phe407, Val410, Met411, Arg412, Leu415, Hec606
MPO residues that were observed to interact with 7GD (co-crystalized inhibitor) based on LigPlot+ analysis of the crystal (PDB ID: 6WYD) are marked in boldface type. MPO residues that were reported to be involved in catalysis [44] are underlined.
Table
Table 3. Intermolecular interactions between seed-derived antioxidant peptides and XO.
Table 3. Intermolecular interactions between seed-derived antioxidant peptides and XO.
PeptideBinding Affinity (kcal/mol)Interaction with XO
Hydrogen Bond aHydrophobic InteractionSalt BridgeExternal Bond
LYSPH−6.2Ser876Leu648, Phe649, Glu802, Leu873, His875, Ser876, Glu879, Phe914, Phe1009, Thr1010, Val1011, Pro1012, Phe1013, Leu1014, Ala1078, Ala1079--
FGPEMEQ−5.9Glu802, Ser876(2), Ala1079Leu648, Phe649, Gln767, Phe798, Gly799, Glu802, Thr803, Leu873, His875, Ser876, Glu879, Arg880, Ala910, Phe911, Arg912, Phe914, Phe1009, Thr1010, Val1011, Pro1012, Leu1014, Pro1076, Ala1078, Ala1079, Ser1080, Glu1261His875-
LPHFNS−4.7His875, Ser876Leu648, Phe649, Glu802, Leu873, His875, Ser876, Glu879, Arg880, Phe914, Phe1009, Thr1010, Val1011, Pro1012, Phe1013, Leu1014, Ala1078, Ala1079, Glu1261--
AEHGSLH−3.4Glu879Leu648, Phe649, Leu712, Glu802, Leu873, His875, Ser876, Glu879, Phe914, Phe1009, Thr1010, Val1011, Pro1012, Phe1013, Leu1014, Pro1076, Tyr1140, Phe1142His875-
PSYLNTPLL3.0Asn768, Asp872, Ser876(2), Arg880, Thr1010(2)Leu648, Phe649, Leu712, Asn768, Glu802, Thr803, Arg871, Asp872, Leu873, Ser874, His875, Ser876, Glu879, Arg880, Phe914, Ser1008, Phe1009, Thr1010, Val1011, Pro1012, Phe1013, Leu1014, Pro1076, Ala1079, Tyr1140, Phe1142, Glu1261-Ala1079
a Number in brackets indicates the number of hydrogen bonds formed with the same residue of XO. XO residues that were reported to bind to quercetin (co-crystalized inhibitor in the crystal PDB ID 3NVY) [12] are marked in boldface type. XO residues that were reported to be involved in catalysis [12] are underlined.
Table
Table 4. Intermolecular interactions between seed-derived antioxidant peptides and p47phox.
Table 4. Intermolecular interactions between seed-derived antioxidant peptides and p47phox.
PeptideDocking ScoreInteraction with p47phox
Hydrogen BondHydrophobic Interaction
PSYLNTPLL−216.493Trp193, Trp263Thr170, Ser171, Ser173, Glu174, Gly192, Trp193, Trp204, Pro206, Ser208, Phe209, Lys235, Tyr237, Asp261, Gly262, Trp263, Met278, Tyr279
LPHFNS−195.377Tyr279Tyr167, Thr170, Ser191, Trp193, Pro206, Tyr237, Asp243, Glu244, Asp261, Gly262, Trp263, Tyr274, Pro276, Tyr279
LYSPH−185.715Thr170, Ser208Tyr167, Thr170, Ser171, Ser173, Glu174, Trp193, Trp204, Pro206, Ser208, Phe209, Trp263, Met278
AEHGSLH−181.729Trp263Tyr167, Thr170, Ser171, Glu174, Trp193, Glu241, Asp243, Trp263, Tyr274, Pro276, Met278
FGPEMEQ−175.680Thr170, Trp193Thr170, Ser173, Glu174, Gly192, Trp193, Trp204, Pro206, Phe209, Asp261, Met278, Tyr279
p47phox residues that were reported to bind to the ligand p22phox-derived proline-rich peptide in the crystal (PDB ID: 1WLP) [14] are marked in boldface type. Key residues that were reported for high-affinity binding between p47phox and p22phox as evidenced by mutagenesis studies [14] are underlined.
Table
Table 5. The predictions of the cell-penetrating potential, blood-brain barrier penetrating potential, plasma half-life, and tolerance to in silico GI digestion of the five selected seed-derived antioxidant peptides.
Table 5. The predictions of the cell-penetrating potential, blood-brain barrier penetrating potential, plasma half-life, and tolerance to in silico GI digestion of the five selected seed-derived antioxidant peptides.
PeptideCell-Penetrating PotentialBlood-Brain Barrier Penetrating PotentialPlasma Half-Life (Seconds)Tolerance to In Silico GI Digestion
AEHGSLHNoNo828.91No
FGPEMEQNoYes796.21No
LPHFNSYesYes823.51No
LYSPHYesYes832.41No
PSYLNTPLLYesYes833.41No
Table
Table 6. Changes in the binding free energies (∆∆G) of the LYSPH- and PSYLNTPLL-protein complexes as revealed by the computational alanine scanning of the peptide residues.
Table 6. Changes in the binding free energies (∆∆G) of the LYSPH- and PSYLNTPLL-protein complexes as revealed by the computational alanine scanning of the peptide residues.
PeptideResidue∆∆G (kJ/mol)
Keap1MPOXOp47phox
LYSPHLeu8.45967.77646.3704-
Tyr14.323223.045816.8951-
Ser2.5540−0.1909−0.1396-
Pro5.03705.27646.3965-
His0.63993.101414.2670-
PSYLNTPLLPro---0.0856
Ser---0.5437
Tyr---27.1615
Leu---−0.4367
Asn---0.7800
Thr---1.1903
Pro---2.9660
Leu---2.7811
Leu---0.4836
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Chai, T.-T.; Koh, J.-A.; Wong, C.C.-C.; Sabri, M.Z.; Wong, F.-C. Computational Screening for the Anticancer Potential of Seed-Derived Antioxidant Peptides: A Cheminformatic Approach. Molecules 2021, 26, 7396. https://0-doi-org.brum.beds.ac.uk/10.3390/molecules26237396

AMA Style

Chai T-T, Koh J-A, Wong CC-C, Sabri MZ, Wong F-C. Computational Screening for the Anticancer Potential of Seed-Derived Antioxidant Peptides: A Cheminformatic Approach. Molecules. 2021; 26(23):7396. https://0-doi-org.brum.beds.ac.uk/10.3390/molecules26237396

Chicago/Turabian Style

Chai, Tsun-Thai, Jiun-An Koh, Clara Chia-Ci Wong, Mohamad Zulkeflee Sabri, and Fai-Chu Wong. 2021. "Computational Screening for the Anticancer Potential of Seed-Derived Antioxidant Peptides: A Cheminformatic Approach" Molecules 26, no. 23: 7396. https://0-doi-org.brum.beds.ac.uk/10.3390/molecules26237396

Article Metrics

Back to TopTop