Next Article in Journal
Bacterial Strains from Saline Environment Modulate the Expression of Saline Stress-Responsive Genes in Pepper (Capsicum annuum)
Previous Article in Journal
Influence of the Abiotic Elicitors Ag Salts of Aspartic Acid Derivatives, Self-Organized in Nanofibers with Monomeric and Dimeric Molecular Structures, on the Antioxidant Activity and Stevioside Content in Micropropagated Stevia rebaudiana Bert.
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Antioxidant Properties and Prediction of Bioactive Peptides Produced from Flixweed (sophia, Descurainis sophia L.) and Camelina (Camelina sativa (L.) Crantz) Seed Meal: Integrated In Vitro and In Silico Studies

by
Na Thi Ty Ngo
,
Tharindu R. L. Senadheera
and
Fereidoon Shahidi
*
Department of Biochemistry, Memorial University of Newfoundland, St. John’s, NL A1C 5S7, Canada
*
Author to whom correspondence should be addressed.
Plants 2023, 12(20), 3575; https://doi.org/10.3390/plants12203575
Submission received: 13 September 2023 / Revised: 5 October 2023 / Accepted: 9 October 2023 / Published: 14 October 2023

Abstract

:
Flixweed (sophia) seed meal and camelina, both by-products of oil processing, were employed to generate protein hydrolysates by applying Flavourzyme and Alcalase. This study aimed to integrate in vitro and in silico methods to analyze sophia and camelina protein hydrolysates for releasing potent antioxidative, dipeptidyl peptidase IV (DPP IV) inhibitors and angiotensin-converting enzyme (ACE) inhibitory peptides. In vitro methods were used to investigate the antioxidant potential of sophia/camelina protein hydrolysates. Bioinformatics techniques, including Peptideranker, BIOPEP, Toxinpred, AlgPred, and SwissADME, were employed to obtain the identification of bioactive peptides produced during the hydrolysis process. Protein hydrolysates produced from sophia and camelina seed meal exhibited higher ABTS and DPPH radical scavenging activities Ithan their protein isolates. Among the produced protein hydrolysates, Alcalase-treated samples showed the highest oxygen radical absorbance capacity and hydroxyl radical scavenging activity. In addition, sophia/camelina hydrolysates prevented hydroxyl and peroxyl radical-induced DNA scission and LDL cholesterol oxidation. In silico proteolysis was conducted on Alcalase-treated samples, and resultant peptides showed potential DPP IV and ACE-inhibitory activities. Identified peptides were further assessed for their toxicity and medicinal properties. Results indicate that all digestive-resistant peptides were non-toxic and had desirable drug-like properties. The findings of this study suggest that sophia/camelina protein hydrolysates are promising candidates for functional foods, nutraceuticals, and natural therapeutics.

1. Introduction

Natural plant-derived proteins, as opposed to those from animal sources, are currently attracting a lot of attention due to their sustainable availability and economic attraction [1]. Food-derived protein hydrolysates and their bioactive peptides are potentially active components suitable for formulating nutraceuticals and functional foods because they are safe and have a wide range of possible health benefits [2,3]. Thus, diversified sources of plant proteins increasingly address the global need for ingredients in the food industry. Furthermore, the identification of bioactive peptides encoded in oilseed processing waste has received much attention in recent years. It has also been highlighted as a viable biorefinery strategy for addressing the detrimental effects of unregulated waste disposal on the environment [3]. In this regard, camelina and sophia seed meals have emerged as potential functional food ingredients due to their nutritional value and beneficial health effects [4,5].
Protein hydrolysates are prepared by cleaving the peptide bonds, either enzymatically or chemically. The use of enzyme hydrolysis to manufacture protein hydrolysates has received much attention to create functional food ingredients and nutraceuticals for reducing disease risk and promoting health [6,7]. Many researchers have reported that peptides prepared by enzymatic hydrolysis from plant sources such as date seed, soybean, flaxseed, and quinoa exhibit a range of physiological and biological attributes, including antioxidative, dipeptidyl peptidase-IV (DPP IV) inhibition, and angiotensin-I-converting enzyme (ACE) inhibition [6,8,9,10]. According to the World Health Organization (WHO), the proportion of individuals affected by type 2 diabetes (T2D) and cardiovascular diseases (CVD) is growing tremendously. However, some synthetic therapeutic medications have been produced to treat these diseases, and many of them are considered harmful due to adverse effects such as inflammatory responses, taste disturbances, nausea, headaches, and allergic reactions [3,8]. Furthermore, oxidative damage is linked to several pathogenic problems, including cardiovascular disease, diabetes, cancer, and various chronic and inflammatory ailments [8]. For human disease prevention and management, natural antioxidants such as plant protein hydrolysates have arisen as an alternative to chemotherapy, as well as a strategy to maintain a healthy lifestyle.
Traditional methods have been employed for the discovery and creation of novel bioactive peptides, which include numerous processes for isolating, refining, and characterizing bioactive peptides. Nevertheless, this method is time-consuming and laborious. Currently, these problems have been solved by using bioinformatic methods that save time and work better than traditional methods, including the contributions of weakly polar interactions to structure and activity [11,12]. Generally, bioinformatics methods are designed to investigate bioactive peptides already present in protein sequences. This method may exclude potential precursor peptides as it only focuses on already studied protein sequences available in open-access databases. However, these virtual screening methods also detail the specificities and types of proteases used to generate bioactive peptide sequences [10]. In addition, the mass spectrometry (MS)-based proteomics method including sequence and conformation-specific structure determinations, a relatively recent technology, has been widely adopted and applied to the development of protein analysis to detect proteins and peptides from dietary sources [13,14]. After peptides are identified, the sequence can be further observed using bioinformatic tools, including UniProtKB, SwissProt, TrEMBL, BIOPEP, and PepBank, to assess the potential of proteins as sources for generating bioactive peptides [11,15]. A limited number of studies have reported the use of camelina and sophia meals as precursors for bioactive peptide. No study has so far been carried out to explore drug likeliness of camelina and sophia meal-derived bioactive peptides. This study marks the inaugural use of sophia and camelina seed meals for the creation of protein hydrolysates utilizing Flavourzyme and Alcalase, along with an examination of their antioxidant properties. This is also the first study to use in silico techniques to assess the antioxidative, ACE-inhibitory, and DPP IV-inhibitory activities of hydrolysates derived from camelina and sophia meals. The in silico techniques including PeptideRanker, BIOPEP, ToxinPred, and SwissADME, allergy FP tools have been used to screen these potential activities of camelina/sophia peptides to assess their toxicity and physicochemical characteristics, and drug likeliness were among variables tested. Therefore, the present study reports the importance of integrating in vitro and in silico techniques to identify the multifunctional peptides derived from underutilized camelina and sophia seed meals.

2. Results and Discussion

2.1. ABTS and DPPH Radical Scavenging Activities

The oxidative damage that occurs to the human body is caused by reactive oxygen species (ROS), namely superoxide anion (O2•−), peroxyl radicals (ROO), hydroxyl radicals (HO), hydrogen peroxide (H2O2), and peroxynitrite (ONOO). DPPH and ABTS are frequently employed for assessing the antioxidant potential of bioactive substances, including peptides and phenolic compounds [16,17,18,19,20].
Figure 1 and Figure 2 display the DPPH and ABTS radical scavenging activity of sophia/camelina hydrolysates and those of protein isolates; the hydrolysates exhibited significantly elevated values (p < 0.05) in comparison to their corresponding isolates. The significant difference observed between the hydrolyzed and non-hydrolyzed samples implies that sophia and camelina antioxidant peptides produced during enzyme digestion can donate electrons to decrease the content of ABTS•+. The ABTS radical scavenging activity of hydrolysates ranged between 26 and 46%, while DPPH radical scavenging activity varied from 17 to 36%. Among the tested samples, Alcalase-treated sophia protein isolate (SPI) exhibited the most notable (p < 0.05) capability in scavenging ABTS and DPPH radicals, whereas Flavourzyme hydrolyzed SPI showed the lowest radical scavenging activity. According to the previous study, the ability of food proteins and their hydrolysates to scavenge free radicals is influenced by various factors, including the amino acid composition, peptide sequences, and their size, as well as the protease’s specificity [21,22]. For example, aromatic amino acid residues (Tyr, His, Trp, and Phe) can donate hydrogen atoms to electron-deficient radicals through a procedure known as resonance stabilization. This enhances the radical-scavenging properties of the amino acid residues. Alkaline proteases, such as Alcalase and Flavourzyme, produce hydrolysates with more antioxidant activity compared to enzymes such as Neutrase, papain, bromelain, and pepsin [23]. This highlights that the observed trends in the radical scavenging abilities of sophia protein hydrolysates, particularly concerning DPPH and ABTS, are generated upon the treatment with Alcalase, which are comparable to those observed for camelina protein hydrolysates that were previously reported [21]. Therefore, our findings support the previously documented research into the radical scavenging potential of hydrolysates derived from diverse plant proteins [18,24].

2.2. Oxygen Radical Absorbance Capacity (ORAC) and Hydroxyl Radical Scavenging Activity

The most damaging forms of oxidative stress are caused by hydroxyl radicals, which interact with involving biomolecules such as amino acids, DNA, proteins, and membrane lipids to cause severe cell damage [21]. Moreover, oxygen radicals play a major role in oxidative stress, resulting from an imbalance between reactive oxygen species (ROS) and antioxidants, and leading to fatal diseases such as cancer, cardiovascular diseases, arthritis, diabetes, and aging [6]. In the present study, the hydroxyl radical scavenging and oxygen radical absorbance of camelina and sophia protein hydrolysates are shown in Table 1. The ability to scavenge hydroxyl radicals was highest in sophia and camelina protein hydrolysates prepared with Alcalase around 3.13 and 3.21 (µm of Histidine/mg protein), respectively, and lowest in those prepared with Flavourzyme at 1.91 and 2.5 (µm of Histidine/mg protein). There was a range of values for the oxygen radical absorption capacity of camelina/sophia protein hydrolysates, from 0.51 to 1.13 (µm of trolox per mg of protein). The results indicated that the ability to scavenge hydroxyl radicals and oxygen radical absorbance capacity of camelina hydrolysates are higher than that of sophia hydrolysates. According to the findings of the research, hydrolysates that were generated from Alcalase gave the greatest value of hydroxyl radical scavenging ability as well as oxygen radical absorbance capacity. A significant difference (p < 0.05) existed in oxygen radical absorbance capacity and hydroxyl radical scavenging ability of sophia protein hydrolysates, which decreased in the order of AL > AL + FL > FL. Similar results for Alcalase-assisted hydrolyzed products of date seed meals and camelina meals have been reported in previous studies [6,21].

2.3. Metal Chelation Activity

The Fenton reaction, which involves lipid oxidation, leads to the degradation of hydroperoxide into volatile compounds. The chelation of transition metal ions (Fe, Cu, and Co) by hydrolysates reduces the intensity of the pink Ferrozine-ferrous complex [19]. The chelation of metal ions by sophia protein hydrolysates revealed significantly better metal-chelating activity than their respective protein isolates, as shown in Figure 3, and the findings regarding camelina hydrolysates have been published [21]. Furthermore, the chelation of ferrous ions by sophia protein hydrolysates varied greatly in the order of AL> AL + FL > FL. The highest metal chelation activity of sophia hydrolysates was observed in hydrolysates prepared by Alcalase, while AL + FL design produced the highest activity in camelina hydrolysates. This capacity correlates with the existence of amino acid residues such as tryptophan, aspartate, cysteine, histidine, and glutamate on the surface of proteins and polypeptides, which have been found to bind divalent metal ions. Additionally, because of the imidazole group that it possesses, the N-terminal presence of histidine demonstrates significant potential for metal ion chelation [25]. The findings suggest that peptides can interact with transition metal ions, potentially slowing down the oxidation process.

2.4. Inhibition against Copper-Induced LDL-Cholesterol Oxidation

LDL cholesterol plays a pivotal role in delivering triacylglycerols and cholesterol to body cells. The oxidation of plasma low-density lipoprotein (LDL), triggered by metal ions or reactive oxygen species (ROS), is a significant risk factor for atherosclerosis and coronary heart disease [26]. Cupric ion has been proven effective in initiating the oxidation of human LDL without EDTA, and the presence of cholesteryl linoleate and cholesteryl arachidonate in the LDL core supplies ample substrate for lipid peroxidation. According to the assay, Cu2+ is employed to induce in vitro human LDL-cholesterol oxidation, and the outcomes of conjugated dienes are measured at 234 nm. Carnosine is used as a positive control. Table 2 indicates that sophia, camelina hydrolysates, and carnosine displayed increased levels after 12 h of incubation, being 39.04−79.97 and 78.59%, respectively. Flavourzyme displayed the lowest amount of conjugated dienes production in comparison to all other treatments for protein hydrolysates. Camelina hydrolysates prepared by Alcalase have demonstrated their highest efficacy in inhibiting cupric ion-mediated LDL oxidation (79.97%). Camelina and sophia protein hydrolysates were found to have an inhibitory effect on human LDL cholesterol oxidation, which may be linked to their capacity to scavenge free radicals and chelate copper ions in LDL. In comparison to date seed protein, sophia/camelina protein hydrolysates were found to have a similar inhibitory impact on LDL oxidation [27].

2.5. Inhibition of Peroxyl and Hydroxyl Radical-Induced Supercoiled DNA Strand Scission

Free radicals such as hydroxyl and peroxyl radicals generated in living cells possess a more significant reduction potential and can engage with biomolecules, leading to DNA damage at both the phosphate backbone and nucleotide bases. This can lead to mutation, carcinogenesis, and other pathological processes. As a result, limiting DNA oxidation is critical for maintaining cell health. In this experiment, AAPH, a water-soluble azo compound with cationic properties, served as the source of alkoxyl and alkyl peroxyl radicals. It was employed to induce single-strand cleavage of supercoiled plasmid DNA and evaluate the protective effects of hydrolysates against DNA oxidation. This was carried out to find out if hydrolysates have a protective effect against DNA oxidation. Carnosine was chosen as a positive control due to its natural antioxidant properties. The results are shown in Table 3. Sophia and camelina protein hydrolysates were found to reduce hydroxyl radical-induced DNA scission by 19.85 to 53.01% and 28.25 to 59.52%, respectively. Additionally, these hydrolysates exhibited a decrease in peroxyl radical-induced DNA oxidation by 62.83 to 87.10% and 75.82 to 88.77%, respectively. In both hydroxyl radical- and peroxyl radical-induced oxidation, hydrolysates that were produced by Alcalase and the combination of Alcalase and Flavourzyme showed a more potent inhibitory effect than carnosine. Comparable results were observed in the case of date seed protein hydrolysates, with supercoiled DNA retention ranging from 13 to 33% for hydroxyl radicals and from 47 to 83% for peroxyl radicals [27]. As a result, sophia/camelina protein hydrolysates have the potential to act as an appropriate dietary source for the protection of DNA against oxidation.

2.6. Amino Acid Composition

The amino acid content of protein hydrolysates has a major impact on their antioxidant activity and plays an important part in the physiological benefits that these substances provide [24]. Table 4 summarizes the amino acid composition of sophia protein hydrolysates prepared by Alcalase, and these correspond with those of Sophia seeds in a previous study [28]. Furthermore, the amino acid profile of camelina and sophia protein hydrolysates almost meets the recommendations of the Agriculture Organization and the World Health Organization/Food (FAO/WHO) for the amounts of most essential amino acids in food products. However, the levels of methionine and histidine in sophia hydrolysates are quite low compared to those recommended by WHO/FAO. In addition, the amounts of most amino acids in camelina hydrolysates are significantly higher than that of sophia hydrolysates. The amino acid profile of camelina and sophia protein hydrolysates exhibited a high concentration of leucine, tryptophan, arginine, glutamic, and aspartic. The high amounts of leucine found in hydrolysates may promote scavenging of superoxide radicals. Furthermore, negatively charged amino acids including Asp and Glu have strong antioxidant properties because they are rich sources of electrons that can be donated to free radicals, while tyrosine residue can act as a robust hydrogen donor, and histidine can be acknowledged for its effective radical scavenging activity, attributed to the breakdown of its imidazole ring [21,25]. As a result, the amino acid profile has been proposed as a method for predicting the antioxidant activity of protein hydrolysates. According to the findings of this study, camelina/sophia protein hydrolysates have the potential to serve as a valuable contributor of protein to food.

2.7. Peptide Ranker Analysis

As shown in the research results decsribed above, the protein hydrolysate produced by Alcalase has high antioxidant activity. In this connection, peptideRanker (PepRank) is an in silico tool that has been used to predict the potential of bioactive peptides prepared using Alcalase. In general, any peptide with a threshold of more than 0.5 is considered bioactive. In addition, the false optimistic predictions reduce from 11 and 16% at a 0.5 threshold to 2 and 6% at a 0.8 threshold for long and short peptides, respectively. PepRank is more dependable that the peptide is bioactive when the expected probability is close to 1 [29]. Following the LC-MS/MS analysis, identified peptides from sophia/camelina protein hydrolysates were subjected to a peptide ranker tool to predict possible activity with a 0.8 threshold score. Table 5 shows 9 and 33 peptide sequences that were found to be potentially bioactive in sophia and camelina hydrolysates, respectively. These were with a threshold score of 0.8 to 0.96. FGFGPGL and GPPSGGGGGGGGGGGGGK peptides had the highest threshold score (0.96), while the lowest score was for SFPLPEL (0.8). Previous studies have shown that the biological activity of peptide sequence is determined by amino acid position and composition [30,31]. Therefore, the BIOPEP database was employed to investigate the antioxidative, ACE-inhibitory, and DPP IV-inhibitory effects of expected peptides [32].

2.8. In Silico Predictions of Potential Bioactive Peptides from Sophia/Camelina Protein Hydrolysate

Theprediction of the bioactive potential of selected peptides was demonstrated using the tool “Profiles of potential biological activity” from a BIOPEP-UWM database. This tool matches the predicted peptide sequences with antioxidative potential as well as ACE- and DPP IV-inhibitory activities to the available literature of diverse plant-, animal-, and other food-derived peptides [31,32]. Selected peptides were screened to investigate the type and locations of the bioactive fragments in peptide sequences.
According to BIOPEP-UWM database, Table 6 indicate profiles of potential biological activity of camelina/sophia protein hydrolysates, including antioxidative, ACE- and DPP IV-inhibitory activities. A total of 5 of the 9 peptide sequences and 12 of the 33 peptides sequences released from sophia and camelina hydrolysates, respectively, were shown to have the potential to be antioxidative, and all selected sophia/camelina-derived bioactive peptides were found to have the ability to inhibit ACE and DPP IV activities (Table 6). All active fragments are dipeptides or tripeptides except for GFGPGL, which is easily absorbed in the gastrointestinal tract. The total number of ACE-inhibitory and DPP IV-inhibitory dipeptides is significantly higher than the total of tripeptides.
Interestingly, tryptophan (W) was the most common amino acid residue present in antioxidative peptides derived from sophia hydrolysates, such as FVPVTGLWM, WYTICICIL, and RAPWLEPL. In addition, GPP was the most common bioactive sequence identified in peptides derived from camelina hydrolysates. A similar result was found in antioxidant peptides from protein hydrolysate of bluefin leatherjacket [33]. In addition, the antioxidant activities of IY and IR were also predicted for jack bean and canavalin protein [34]. Peptides with more hydrophobic and aromatic amino acids were shown to have high antioxidant activity, and their presence is considered a critical attribute to the antioxidative property of bioactive peptides. In addition, the presence of hydrophobic amino acids including alanine (A), valine (V), leucine (L), isoleucine (I), proline (P), phenylalanine (F), tyrosine (Y), methionine (M), histidine (H), and tryptophan (W) is responsible for the antioxidant activity [35]. For instance, tyrosine and phenylalanine residues can act as potent hydrogen donors, and histidine can be credited with strong radical scavenging activity due to its imidazole ring [3]. Furthermore, the active amino acid residues located at the C-terminal as tyrosine in peptide AWPDKNPFFPSDPY, methionine in peptide GGGGGGGGPPAMSM, and GLDPPDLPM contributed to their antioxidant activity. Several studies have indicated that peptides having a tyrosine residue located at the C-terminus strongly scavenged hydroxyl radicals and peroxyl radicals as well as hydrogen peroxide [36]. As a result, the type and position of amino acid residues of the peptide could be comprehensively applied to identify novel antioxidant sources and obtain desirable antioxidants. Previous research has shown that several peptides found in sophia/camelina protein hydrolysates are potent antioxidants, and these peptides come from a variety of animal and plant sources. According to the database, AY, PW, LW, WY, IY, and IR were identified from okara protein, buckwheat, marine bivalve and bovine beta-lactoglobulin, jack bean, and canavalin protein, respectively.
The highest number of potential ACE-inhibitory dipeptides and tripeptides (14) was found in the sequence SSTSGPAFNAGRSIWLPGWL. Numerous peptides found in sophia/camelina protein hydrolysates have previously been recognized as strong ACE inhibitors in diverse peptides from animal and plant origins. For example, GP, PL, GPL, and LPG were identified from Alaskan pollack skin, while CF, SF, GPP, and FVP were also previously confirmed from shark meat hydrolysate, garlic, wheat gliadin, and soya milk, respectively. Previous research has demonstrated that characteristics of peptides, such as chain length, composition, and sequence, can significantly affect their ACE-inhibitory activities. For example, the presence of hydrophobic amino acids such as proline (P), phenylalanine (F), isoleucine (I), valine (V), methionine (M), leucine (L), and alanine (A) close to the C-terminus of peptides promotes their binding to ACE and indeed inhibits the enzyme more effectively [6,37]. Additionally, the effect of proline on ACE-inhibitory activity is associated with its imidazole ring, which interacts strongly with the amino acid residues at the active centers of ACE. Proline was present in most ACE-inhibitory peptide sequences (Table 6). Thus, peptides derived from sophia/camelina meals are effective in inhibiting ACE. Alcalase was the most promising enzyme used for the isolation of both antioxidant and antihypertensive peptides, consistent with previous research on bioactive peptides derived from edible seeds [3].
According to Table 6, all the active segments are dipeptides or tripeptides, with dipeptides being the most abundant and effective DDP IV inhibitors. In addition, in active dipeptides and tripeptides, proline and alanine were the two most common amino acids such as GP, GA, PG, PA, AG, IP, LP, AP, PP, VP, APG, GPA, and PPL. This corresponds to some previous research results that peptides having proline or alanine residues in their sequences have potential DPP IV-inhibitory activity [38,39]. Furthermore, the potent DPP IV-inhibitory peptides have been found to possess a branched-chain amino acid (leucine, isoleucine, or valine) or an aromatic residue with a polar group in the sidechain (tryptophan) at their N-terminal position [40]. The interest in multifunctional peptides has been increasing as more of them have been discovered from food protein. According to the findings, it has been demonstrated that multifunctional peptides exhibit more than one major physiologically relevant bioactive characteristic. Because of this, there is more interest in multifunctional peptides now found in food proteins. Peptides, on the other hand, are susceptible to degradation during digestion in the gastrointestinal tract, and the biological activity of peptides can be either activated or inactivated depending on the structure and function of the peptides. Since in silico digestion is an appropriate method for predicting the release and resistance of bioactive peptides, the predicted peptides that were derived from the sophia/camelina protein hydrolysates were put through this process.

2.9. In Silico Simulated Gastrointestinal (GI) Digestion of Sophia/Camelina Protein Hydrolysates

Peptides that make up the building blocks of the human body are created when protein products are digested. To determine whether peptides can provide a health benefit, they must be subjected to simulated digestive conditions [2,41]. In this perspective, simulating GI digestion in silico is a valuable tool for determining the bioactive potential of sophia/camelina peptide before assessing its bioavailability and bioaccessibility in vivo. The sequences of protein hydrolysates were subjected to in silico proteolysis using pepsin (EC 3.4.23.1), trypsin (EC 3.4.21.4), and chymotrypsin (EC 3.4.21.1) for the prediction of theoretically released peptides by the enzymatic action program available in BIOPEP tool. In the present study, in silico proteolysis demonstrates the ability of given enzymes to release the antioxidative and the ACE-inhibitory peptides from sophia/camelina protein hydrolysates.
After the simulated digestion process of sophia hydrolysates, none of the peptides displayed antioxidative potential, as predicted. This could be because there is insufficient information in the antioxidative peptide database, or it could be because gastrointestinal condition cleaves certain peptides [11]. The results showed that three out of the nine peptides identified from sophia hydrolysates were anticipated to release potent ACE-inhibitory and DPP IV-inhibitory fragments from their original sequences, respectively (Table 7). It is interesting to note that FGFGPGL possesses powerful ACE-inhibitory and DPP IV-inhibitory properties. GF is an active fragment sequence that can inhibit ACE as well as DPP IV. According to the findings given in Table 7, the bioactive potential upon GI digestion of peptides derived from camelina protein hydrolysate was reported. The results show that two peptides (AWPDKNPFFPSDPY, and SFPLPEL) have the potential antioxidant after the simulated digestion process. A total of 13 out of the 33 peptides identified from camelina hydrolysates were predicted to release potent ACE-inhibitory fragments from their original sequences. The ACE-inhibitory dipeptides are the most abundantly released from camelina protein hydrolysates, whereas few tri-peptides peptides sequences were found. Importantly, the frequency of release of potent ACE inhibitor peptides from each identified peptide is different because some peptides have more than one active fragment embedded in their sequence. For instance, SFPLPEL sequence was predicted to release three active ACE-inhibitory sequences, namely PL, SF, and PEL.
The “enzyme(s) action” option of the BIOPEP database is used to determine the degree of hydrolysis (DHt), the release frequency (AE), and the value of relative release frequency (W). AE and W are the two most important parameters to consider when estimating the potential release of bioactive peptides from protein sequences [31]. In silico digestion revealed that sophia and camelina hydrolysates have degrees of hydrolysis ranging from 12.50 (DIPPPRGPL) to 66.67% (FHWDLPQ), as indicated in Table 8. The highest AE value (0.29) was observed in peptide SFPLPEL, whereas the lowest (0.06) was FGGYAPGILSPSPAML. Regarding camelina hydrolysates, the highest relative frequency of the release of fragments with a given selected enzymes (W) was 1.0 in AWPDKNPFFPSDPY, whereas AAMGGFPGGGGGAHALG (0.06) is the lowest value. After in silico simulated digestion of camelina meals derived peptides, most of the bioactive fragments possess proline (P), glycine (G), leucine (L), and phenylalanine (F) in their dipeptide and tripeptide sequences. The bioactive motifs that result contain VF, GPL, AW, DF, GL, GR, PPL, PQ, GGY, SF, DM, IL, GW, PL, and PEL. Previous studies using in silico methods to identify biopeptides in proteins from a variety of plants and animals have shown similar results [31,42]. According to the findings of several studies, the presence of hydrophobic amino acids (such as leucine, serine, and methionine) or aromatic amino acids (such as phenylalanine) in this peptide can boost the inhibition of ACE and DPP IV [41,43]. In addition, peptides that contain proline in their sequences have the potential to be highly effective ACE and DPP IV inhibitors [44]. The multiple bioactivities displayed by these peptides can increase their impact on ameliorating more than one disease target or various symptoms since many human diseases are interrelated in etiology and progress. According to these findings, sophia/camelina protein hydrolysates may be employed as multifunctional components in nutraceutical or functional food products.

2.10. Toxicity and Allergenicity Prediction of Sophia-Derived Bioactive Peptide Fractions after In Silico Digestion

The toxicity of bioactive peptides is a major concern when developing peptide-based nutraceuticals or functional food ingredients. The ability to predict the toxicity of therapeutic peptides prior to their synthesis is critical for reducing the time and money spent developing peptide-based drugs [45]. In addition, the majority of allergens are proteins present in a variety of plant and animal sources, resulting in a global problem that is continuously deteriorating [31]. As a result, all peptides derived from sophia/camelina protein hydrolysates following in silico digestion must be assessed for potential toxicity and allergenicity using the ToxinPred and allergen FP v.1.0 tools, respectively [31,46,47]. Table 9 illustrates that the low molecular weight dipeptides (GF, SL, GM, VF, AL, VL, AW, PF, PY, DF, GL, PM, GR, IH, QF, QW, PQ, SF, DM, IL, GW, and PL) and tripeptides (PGL, GPL, PPL, GGY, and PEL) that are released from sophia and camelina protein hydrolysates in silico are non-toxic. Non-toxic peptides are primarily composed of amino acids such as Val, Thr, Arg, Gln, Met, Leu, Lys, Ile, Phe, and Ala. These findings are comparable to those obtained with peptides derived from other protein sources, such as tubers and quinoa, rice bran, and flaxseed protein, all of which were found to be non-toxic [31,45,48,49].
According to Table 9, 4 out of 5 sophia-derived bioactive peptide fractions were classified as probable non-allergens, while 16 of 23 camelina-derived bioactive peptide fractions were classified as probable allergens and 6 peptides were probable non-allergens. Moreover, only VF is considered as non-allergen. The isoelectric points (pI) of the predicted bioactive peptides were found to be in the pH range of 3.80–10.11. In addition, 18 peptides have an isoelectric point of 5.88 with 0 net charge. The dipeptides (DF and DM) and tripeptides (PEL) showed a negative charge (−1) with a pI of 3.8 and 4.0, respectively, whereas GR had a positive charge (+1) with a pI of 10.11. Thus, bioactive peptides with an alkali isoelectric point had strong water solubility, whereas peptides with an acidic isoelectric point had poor water solubility. Some authors have investigated the charge, peptide sequence, low molecular weight, hydrophobicity, and hydrophilicity as important structural properties for food-derived peptide bioavailability [50,51]. In the present study, the active fragment fraction had molecular weights ranging from 188.25 Da (GL) to 357.44 Da (PEL). The hydrophobicity of the bioactive peptides was shown to be in the range of −0.38 (PQ) to 0.57 (VF), whereas their hydrophilicity varied from −2 (VF) to 0.85 (DM). As a result, these peptides have the potential to be used as functional ingredients in the food industry. Allergen and toxic effects of bioactive peptides can be predicted quickly and appropriately using in silico methods. However, in vitro and in vivo investigations are required to verify their toxicity and allergenicity to confirm the safety of these peptides for human consumption.

2.11. In Silico Evaluation of Drug-Likeness Bioactive Peptides Derived from Protein Hydrolysates

In recent years, ADME computation, which stands for absorption, distribution, metabolism, and excretion, has been very helpful in the in silico prediction of the drug-likeness of bioactive peptides. SwissADME is a web-based tool for predicting and assessing small molecule pharmacokinetics, drug-likeness, and medicinal chemistry friendliness [42,52]. The assessment of drug-likeness, which determines the prospect of a component becoming an oral drug, is one of the SwissADME evaluation methods. Because of its flexibility, this method is considered a cost-effective alternative to time-consuming experimental methods [53,54]. SwissADME was used to study physicochemical features such as the number of rotatable bonds (ROTB), hydrogen bond acceptors (HBA), hydrogen bond donors (HBD), topological polar surface area (TPSA), water solubility, lipophilicity, drug-likeness (Lipinski filter and bioavailability score), and pharmacokinetics (GI absorption). The drug-like characteristics of five multifunctional dipeptides and tripeptides that were released from sophia hydrolysates were used to predict the ADME and pharmacokinetic aspects of these compounds. The in silico physicochemical characteristics, drug-likeness, and pharmacokinetics of sophia peptides are displayed in Table 10 and Figure 4, respectively. All these parameters were examined with reference to captopril, which serves as the standard medication. Similarities could be seen between the drug-likeness and pharmacokinetics of the peptides and those of inhibitor drugs such as captopril. According to the data presented in Table 10, all the predicted dipeptides and tripeptides derived from sophia protein hydrolysates qualified as described by Lipinski. The Lipnski rule of five is considered as the rule of thumb to determine the drug likeliness of a molecule [53]. For instance, the rotatable bonds (RORB) of these peptides that were observed had a range of 0 to 9, anforheir TPSA values were lower than 130 Å2, which led to high gastrointestinal absorption. In addition, Table 10 shows all peptides (GF, SL, PM, M, and PGL) as being highly soluble and having high values for GI absorption. In terms of camelina hydrolysates, Table 10 reveals all peptides without PEL, rotatable bonds observed ranged from 0 to 9. This is consistent with previous studies because compounds with more than 10 rotatable bonds show poor oral bioavailability [52]. Furthermore, TPSA values of two dipeptides (GR, DM, and QW) and two tripeptides (GGY and PEL) are higher than 130 Å2 leading to low gastrointestinal absorption. Table 10 shows that all peptides were highly soluble and 17 out of 23 bioactive di- and tripeptide released from camelina hydrolysates had high values for GI absorption, except for GR, QW, PQ, DM, GGY, and PEL. In addition, all peptides (GF, SL, PM, M, and PGL) released from sophia protein hydrolysates were highly soluble and had high levels of GI absorption. These results showed Lipinski filter and bioavailability score of all the identified peptides were similar to that of captopril. This means that most of the peptides derived from camelina/sophia protein hydrolysates have very good absorption.
In addition, bioavailability as reflected in lipophilicity (LIPO), molecular size (Size), polarity (Polar), insolubility (Insolu), flexibility (Flex), and instauration (Insatu) were observed when radar was displayed for a quick assessment of drug-likeness. Oral bioavailability is best achieved in the area colored in Figure 4. According to Figure 4, all selected di- and tripeptides derived from sophia protein hydrolysates have their optimal values in the highlighted zone. In addition, the optimal values of most di- and tripeptides derived from camelina protein hydrolysates, except GR, DM, QW, GGY, and PEL, are in the colored region. As a result, peptides obtained from sophia and camelina protein hydrolysates have the potential to exhibit drug-like properties which could be exploited in the pharmaceutical industry.

3. Materials and Methods

3.1. Materials

In this study, camelina seeds were sourced from Linnaeus Plant Sciences INC in Saskatoon, SK, Canada, while sophia seeds were acquired from Daghdagh Abad, near Hamedan city in Iran, and were purchased at the Tavazo store in Toronto, ON, Canada. Flavourzyme (1000 LAPU/g) and Alcalase (2.4 AU/g) were procured from Novozymes in Bagsvaerd, Denmark. All chemical reagents used in the study were obtained from Fisher Scientific Ltd. in Ottawa, ON, Canada, or Sigma-Aldrich Canada Ltd. in Oakville, ON, Canada.

3.2. Methods

3.2.1. Preparation of Sophia/Camelina Protein Isolatess

Sophia and camelina protein isolates were obtained from sophia and camelina seed meals using essentially the procedure outlined by Ngo and Shahidi [5] using defatted samples. The extracted proteins were freeze-dried and stored at −20 °C for further analyses.

3.2.2. Preparation of Sophia/Camelina Protein Hydrolysates

The hydrolysis of camelina and sophia protein isolates was performed using Alcalase (0.3 AU/g) and Flavourzyme (50 LAPU/g) under various conditions in a batch-wise manner [6,21]. In the case of enzyme combination treatments, the hydrolysis process began with Alcalase for 2 h (pH 8, 50 °C), followed by the addition of Flavourzyme for an additional 2 h (pH 7, 50 °C). The enzymatic reactions were terminated by heat inactivation at 90 °C for 10 min. The resulting hydrolysates were collected by centrifugation.

3.2.3. DPPH Radical Scavenging Assay

The 2,2-diphenyl-1-picrylhydrazyl (DPPH) radical scavenging activity of both protein hydrolysates and their peptide fractions was assessed using a previously established method [16,55], with slight modifications as reported elsewhere using trolox as a reference [21]. The DPPH radical scavenging activity (%) was calculated in terms of micromoles (μM) of trolox equivalents (TE) per milligram of protein hydrolysate, using the following equation:
D P P H   r a d i c a l   s c a v e n g i n g   a c t i v i t y   ( % )   =   ( A b s o r b a n c e   o f   t h e   c o n t r o l A b s o r b a n c e   o f   t h e   s a m p l e ) A b s o r b a n c e   o f   t h e   c o n t r o l × 100

3.2.4. ABTS Radical Cation Scavenging Assay

The radical scavenging activity of protein hydrolysates was assessed using the 2,2′-azinobis (3-ethylbenzothiazoline-6-sulfonate) radical cation (ABTS•+) method with slight modifications as reported previously [21]. The ABTS radical scavenging activity was expressed in µmol of Trolox equivalents (TE) per milligram (mg) of protein hydrolysates, and calculated using the following equation:
A B T S   r a d i c a l   s c a v e n g i n g   a b i l i t y   ( % )   =   ( A b s o r b a n c e   o f   t h e   b l a n k A b s o r b a n c e   o f   t h e   s a m p l e   a f t e r   6   m i n ) A b s o r b a n c e   o f   t h e   b l a n k × 100

3.2.5. Hydroxyl Radical Scavenging Activity

The hydroxyl radical scavenging capacity was determined using EPR spectroscopy (Bruker E-scan, Bruker Biospin Co., Billericia, MA, USA) as reported by Ngo and Shahidi [21]. The hydroxyl radical scavenging capacity was quantified in terms of micromoles (μM) of histidine equivalents per milligram of protein hydrolysate.

3.2.6. Oxygen Radical Absorbance Capacity (ORAC)

The ORAC assay was conducted utilizing a Fluostar Optima plate reader (BMG Labtech, Durham, NC, USA), following a previously described method with slight adjustments [6]. All chemicals and materials were diluted in a 75 mM/L phosphate buffer at pH 7.0 prior to the experiment. In Costar® 3695 flat-bottomed 96-well black microplates (Corning Incorporated, Corning, NY, USA), camelina and sophia protein hydrolysates (1 mg/mL, 20 μL) were combined with 200 μL of fluorescein (0.11 μM in PBS). The mixture was pre-incubated for 15 min at 37 °C in the built-in incubator. The equipment was programmed to introduce 75 μL of AAPH (17.2 mg/mL in PBS) at the end of the incubation period during the first cycle. Changes in fluorescence were recorded every minute for 25 cycles, with each cycle lasting 210 s. Excitation and emission wavelengths were set at 485 nm and 520 nm, respectively. The ORAC values of protein hydrolysates were quantified as micromoles of trolox equivalents per milligram of protein using the trolox standard curve (ranging from 6.25 to 50 μM).

3.2.7. Metal Chelation Activity

Ferrous ion chelating activity was assessed following a previously established method [6,21] using the trisodium salt of ethylenediaminetetraacetic acid (Na3EDTA) as a reference. The metal ion chelating ability (%) was calculated using the following formula:
M e t a l   c h e l a t i n g   a c t i v i t y   ( % )   =   ( 1 A b s o r b a n c e   o f   t h e   s a m p l e A b s o r b a n c e   o f   t h e   c o n t r o l ) × 100

3.2.8. Cupric Ion-Induced Human Low-Density Lipoprotein (LDL) Peroxidation

The inhibitory activity of camelina/sophia protein hydrolysates against copper-induced LDL cholesterol oxidation was measured according to the reported method [27]. LDL (5 mg/mL) was dialyzed in 10 mM phosphate buffer (pH 7.4, 0.15 M NaCl) using a dialysis tube with a molecular weight cut-off of 12−14 kDa (Fisher Scientific, Nepean, ON, Canada) at 4 °C under a nitrogen blanket in the dark for 12 h. Diluted LDL cholesterol (0.04 mg LDL/mL) was mixed with protein hydrolysate solutions (0.1 mg/mL). The positive control used was carnosine. The reaction was initiated by adding 0.1 mL of 100 μM solution of CuSO4. A blank containing only a sample without LDL or CuSO4 was prepared for each test compound. After incubation of the reaction mixture at 37 ℃ for 12 h, the conjugated dienes formed were recorded at 234 nm using a diode array spectrophotometer (Agilent, Palo Alto, CA, USA). The inhibitory effect of tested samples was expressed as percentage inhibition of conjugated diene formation according to the following equation:
Inhibition (%) = [(Acontrol − Asample)/(Acontrol − Anative)] × 100
Acontrol = (A°0 − A°t)
Asample = (A0 − At)
Acontrol: the absorbance of LDL + CuSO4 + PBS
A sample: the absorbance of LDL + CuSO4 + sample/standard
Anative: the absorbance of LDL + PBS.
A0, At, and A°0, A°t are absorbance values for test samples and control, respectively, measured at zero time and at time t after incubation.

3.2.9. Supercoiled DNA Strand Scission Inhibition

The inhibitory activity, induced by peroxyl and hydroxyl radical, of protein hydrolysates from sophia/camelina isolates against DNA strand scission was determined according to previous studies, with some modifications [27]. The existing bands were visualised using the Alpha-Imager gel documentation system (Cell Biosciences, Santa Clara, CA, USA) under trans-illumination of UV light, and data processing was achieved using the Gel Analyzer 19.1 software. The percentage oxidation inhibition of the supercoiled DNA strand was calculated according to Equation (7).
Inhibition (%) = (DNA samples/DNAblank) × 100

3.2.10. Amino Acid Composition

The amino acid composition analysis of each protein hydrolysate was conducted at the Analytics, Robotics, and Chemical Biology Centre (SPARC BioCentre), located at the Hospital for Sick Children in Toronto, ON, Canada, according to Mohan and Udenigwe [56]. All amino acids with the exception of tryptophan, cysteine, and methionine, were hydrolysed in 6 M HCl. Tryptophan, cysteine, and methionine analyses were carried out separately [21].

3.2.11. Simulated Gastrointestinal Digestion

LC-MS/MS Analysis

The LC-MS/MS analysis was performed at the Analytics, Robotics, and Chemical Biology Centre (SPARC BioCentre) at The Hospital for Sick Children in Toronto, ON, Canada. This analysis utilized a Q_Exactive Orbitrap analyzer equipped with a nanospray source and an EASY-nLC nano-LC system from Thermo Fisher (San Jose, CA, USA), as previously described elsewhere [25]. Data analysis was performed using PEAKS+ studio 11 software (Bioinformatic Solutions, Waterloo, ON, Canada).

Peptide Ranking

PepRank (http://bioware.ucd.ie/~compass/biowareweb/ (accessed on 19 March 2022)) was used to screen peptides for their potential bioactivity which predicts the probability (between 0 and 1) of the peptide being bioactive. To limit the number of false positives, a threshold of 0.8 was selected [29]. The most potent camelina/ sophia-derived bioactive peptides were selected for further analysis.

Prediction of Potential Biological Activity Profile

The bioactive properties of selected peptides, including their antioxidant and antihypertensive activities, as well as their dipeptidyl peptidase IV inhibitory effects, were explored using the “profiles of potential biological activity” feature within the BIOPEP-UWM™ database, accessible at http://www.uwm.edu.pl/biochemia/index.php/en/biopep (accessed on 9 April 2022).

Prediction of the Bioactivities of the Released Peptides after Simulated Digestion

The identified sequences of bioactive peptides underwent in silico hydrolysis using the enzyme action functionality of BIOPEP-UWM™. This simulation employed pepsin (EC 3.4.23.1), trypsin (EC 3.4.21.4), and chymotrypsin (EC 3.4.21.1) as representative digestive enzymes to replicate the process of in vivo digestion. After that, the option to “search for active fragments” was chosen to get a list of peptides with potential activities such as antioxidant and antihypertensive activities and dipeptidyl peptidase IV inhibitory effects [14,52].

Toxicity and Allergenicity Prediction of Camelina/Sophia Bioactive Peptides Released after Simulated Digestion

One of the significant barriers to the sustainable utilization of camelina/ sophia proteins for the development of nutraceuticals and functional food ingredients is their toxicity and allergenicity. Therefore, the potential toxicity and allergenicity of the peptides released after simulated digestive proteolysis were investigated using ToxinPred online tool (https://webs.iiitd.edu.in/raghava/toxinpred/index.html (accessed on 25 April 2022)) [46] and allergen FP v.1.0 tool (http://www.ddg-pharmfac.net/AllergenFP (accessed on 25 April 2022)), respectively [57].

In Silico Physicochemical Properties and Drug-Likeness of Released Peptides after Simulated Digestion

The in silico assessment of the drug-like properties of the identified peptides derived from camelina and sophia was conducted, considering parameters related to absorption, distribution, metabolism, and excretion. This evaluation was carried out using the SwissADME tool, as previously outlined [53], and can be accessed at http://www.swissadme.ch/index.php# (accessed on 8 May 2022).

3.3. Statistical Analysis

All experiments were replicated three times, and the results were presented as the mean  ±  standard deviation. Statistical analysis was carried out using one-way ANOVA, and the means were compared using Tukey’s HSD test (p  <  0.05) with SPSS 16.0 for Windows (SPSS Inc., Chicago, IL, USA)

4. Conclusions

The current study explored the antioxidant potential of protein hydrolysates sourced from sophia and camelina using in vitro analysis. Based on in vitro results, Alcalase assisted protein hydrolysates were selected for further evaluation of their bioactive potential using peptidomics and bioinformatics tools. Findings indicate that protein hydrolysates produced from camelina and sophia seed meals could be utilized as precursors for releasing peptides with dual functions in ACE and DPP IV inhibition. Furthermore, the in silico predictions support the fact that these peptides were non-toxic and have drug-like properties. These findings revealed that camelina and sophia-derived peptides have the potential to serve as natural therapeutic ingredients as well as formulating functional food and nutraceuticals. However, future work should be focused on validating the optimal functionality of these bioactive peptides in vivo.

Author Contributions

F.S. conceptualized and oversaw the project, with N.T.T.N. conducting the experiments and analyzing the data. T.R.L.S. provided support and assistance in the analysis of in-silico data. The manuscript was initially drafted by N.T.T.N., and it underwent review and editing by F.S. The final document was reviewed and approved by all the authors. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Data Availability Statement

Please contact Dr. Fereidoon Shahidi for data requests.

Acknowledgments

We express our gratitude to the Natural Science and Engineering Research Council (NSERC) of Canada for their generous financial support. Additionally, we would like to extend our thanks to the Vietnam International Education Development for providing a scholarship to Ngo.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Deng, Y.; Huang, L.; Zhang, C.; Xie, P.; Cheng, J.; Wang, X.; Li, S. Physicochemical and functional properties of Chinese quince seed protein isolate. Food Chem. 2019, 283, 539–548. [Google Scholar] [CrossRef]
  2. Udenigwe, C.C.; Aluko, R.E. Food protein-derived bioactive peptides: Production, processing, and potential health benefits. J. Food Sci. 2012, 77, R11–R24. [Google Scholar] [CrossRef]
  3. Samtiya, M.; Acharya, S.; Pandey, K.K.; Aluko, R.E.; Udenigwe, C.C.; Dhewa, T. Production, Purification, and Potential Health Applications of Edible Seeds’ Bioactive Peptides: A Concise Review. Foods 2021, 10, 2696. [Google Scholar] [CrossRef] [PubMed]
  4. Russo, R.; Reggiani, R. Seed Protein in Camelina sativa (L.) Crantz var. Calena. Int. J. Plant Soil. Sci. 2015, 8, 1–6. [Google Scholar] [CrossRef]
  5. Ngo, N.T.T.; Shahidi, F. Functional properties of protein isolates from camelina (Camelina sativa (L.) Crantz) and flixweed (sophia, Descurainis sophia L.) seed meals. Food Prod. Process. Nutr. 2021, 3, 31. [Google Scholar] [CrossRef]
  6. Ambigaipalan, P.; Al-Khalifa, A.S.; Shahidi, F. Antioxidant and angiotensin I converting enzyme (ACE) inhibitory activities of date seed protein hydrolysates prepared using Alcalase, Flavourzyme and Thermolysin. J. Funct. Foods 2015, 18, 1125–1137. [Google Scholar] [CrossRef]
  7. Tavano, O.L.; Berenguer-Murcia, A.; Secundo, F.; Fernandez-Lafuente, R. Biotechnological Applications of Proteases in Food Technology. Compr. Rev. Food Sci. Food Saf. 2018, 17, 412–436. [Google Scholar] [CrossRef] [PubMed]
  8. Mune Mune, M.A.; Minka, S.R.; Henle, T. Investigation on antioxidant, angiotensin converting enzyme and dipeptidyl peptidase IV inhibitory activity of Bambara bean protein hydrolysates. Food Chem. 2018, 250, 162–169. [Google Scholar] [CrossRef] [PubMed]
  9. Langyan, S.; Khan, F.N.; Yadava, P.; Alhazmi, A.; Mahmoud, S.F.; Saleh, D.I.; Zuan, A.T.K.; Kumar, A. In silico proteolysis and analysis of bioactive peptides from sequences of fatty acid desaturase 3 (FAD3) of flaxseed protein. Saudi J. Biol. Sci. 2021, 28, 5480–5489. [Google Scholar] [CrossRef]
  10. Mudgil, P.; Kilari, B.P.; Kamal, H.; Olalere, O.A.; FitzGerald, R.J.; Gan, C.-Y.; Maqsood, S. Multifunctional bioactive peptides derived from quinoa protein hydrolysates: Inhibition of α-glucosidase, dipeptidyl peptidase-IV and angiotensin I converting enzymes. J. Cereal Sci. 2020, 96, 103130. [Google Scholar] [CrossRef]
  11. Udenigwe, C.C. Bioinformatics approaches, prospects and challenges of food bioactive peptide research. Trends Food Sci. Technol. 2014, 36, 137–143. [Google Scholar] [CrossRef]
  12. Chass, G.A.; Lovas, S.; Murphy, R.F.; Csizmadia, I.G. The role of enhanced aromatic -electron donating aptitude of the tyrosyl sidechain with respect to that of phenylalanyl in intramolecular interactions. Eur. Phys. J. D 2002, 20, 481–497. [Google Scholar] [CrossRef]
  13. Chass, G.A.; Marai, C.N.J.; Setiadi, D.H.; Csizmadia, I.G.; Harrison, A.G. A Hartree–Fock, MP2 and DFT computational study of the structures and energies of ″b2 ions derived from deprotonated peptides. A comparison of method and basis set used on relative product stabilities. J. Mol. Struct. THEOCHEM 2004, 675, 149–162. [Google Scholar] [CrossRef]
  14. Panjaitan, F.C.A.; Gomez, H.L.R.; Chang, Y.W. In silico Analysis of Bioactive Peptides Released from Giant Grouper (Epinephelus lanceolatus) Roe Proteins Identified by Proteomics Approach. Molecules 2018, 23, 2910. [Google Scholar] [CrossRef] [PubMed]
  15. Agyei, D.; Tsopmo, A.; Udenigwe, C.C. Bioinformatics and peptidomics approaches to the discovery and analysis of food-derived bioactive peptides. Anal. Bioanal. Chem. 2018, 410, 3463–3472. [Google Scholar] [CrossRef]
  16. Chandrasekara, A.; Shahidi, F. Inhibitory activities of soluble and bound millet seed phenolics on free radicals and reactive oxygen species. J. Agric. Food Chem. 2011, 59, 428–436. [Google Scholar] [CrossRef] [PubMed]
  17. Udenigwe, C.C.; Aluko, R.E. Chemometric analysis of the amino acid requirements of antioxidant food protein hydrolysates. Int. J. Mol. Sci. 2011, 12, 3148–3161. [Google Scholar] [CrossRef] [PubMed]
  18. He, R.; Girgih, A.T.; Malomo, S.A.; Ju, X.; Aluko, R.E. Antioxidant activities of enzymatic rapeseed protein hydrolysates and the membrane ultrafiltration fractions. J. Funct. Foods 2013, 5, 219–227. [Google Scholar] [CrossRef]
  19. Shahidi, F.; Zhong, Y. Measurement of antioxidant activity. J. Funct. Foods 2015, 18, 757–781. [Google Scholar] [CrossRef]
  20. Senadheera, T.R.L.; Dave, D.; Shahidi, F. Antioxidant potential and physicochemical properties of protein hydrolysates from body parts of North Atlantic sea cucumber (Cucumaria frondosa). Food Prod. Process. Nutr. 2021, 3, 3. [Google Scholar] [CrossRef]
  21. Ngo, N.T.T.; Shahidi, F. Antioxidant properties of camelina (Camelina sativa (L.) Crantz) protein hydrolysates. J. Food Bioact. 2021, 16, 75–82. [Google Scholar] [CrossRef]
  22. Girgih, A.T.; Udenigwe, C.C.; Aluko, R.E. In vitro Antioxidant Properties of Hemp Seed (Cannabis sativa L.) Protein Hydrolysate Fractions. J. Am. Oil Chem. Soc. 2011, 88, 381–389. [Google Scholar] [CrossRef]
  23. Hu, R.; Chen, G.; Li, Y. Production and Characterization of Antioxidative Hydrolysates and Peptides from Corn Gluten Meal Using Papain, Ficin, and Bromelain. Molecules 2020, 25, 4091. [Google Scholar] [CrossRef] [PubMed]
  24. Alashi, A.M.; Blanchard, C.L.; Mailer, R.J.; Agboola, S.O.; Mawson, A.J.; He, R.; Girgih, A.; Aluko, R.E. Antioxidant properties of Australian canola meal protein hydrolysates. Food Chem. 2014, 146, 500–506. [Google Scholar] [CrossRef]
  25. Udenigwe, C.C.; Udechukwu, M.C.; Yiridoe, C.; Gibson, A.; Gong, M. Antioxidant mechanism of potato protein hydrolysates against in vitro oxidation of reduced glutathione. J. Funct. Foods 2016, 20, 195–203. [Google Scholar] [CrossRef]
  26. Pan, M.; Huo, Y.; Wang, C.; Zhang, Y.; Dai, Z.; Li, B. Positively charged peptides from casein hydrolysate show strong inhibitory effects on LDL oxidation and cellular lipid accumulation in Raw264.7 cells. Int. Dairy J. 2019, 91, 119–128. [Google Scholar] [CrossRef]
  27. Ambigaipalan, P.; Shahidi, F. Antioxidant potential of date (Phoenix dactylifera L.) seed protein hydrolysates and carnosine in food and biological systems. J. Agric. Food Chem. 2015, 63, 864–871. [Google Scholar] [CrossRef]
  28. Mohamed, N.H.; Mahrous, A.E. Chemical Constituents of Descurainia sophia L. and its Biological Activity. Acad. Chem. Globe Publ. 2009, 3, 58–67. [Google Scholar]
  29. Mooney, C.; Haslam, N.J.; Pollastri, G.; Shields, D.C. Towards the improved discovery and design of functional peptides: Common features of diverse classes permit generalized prediction of bioactivity. PLoS ONE 2012, 7, e45012. [Google Scholar] [CrossRef] [PubMed]
  30. Garg, S.; Apostolopoulos, V.; Nurgali, K.; Mishra, V.K. Evaluation of in silico approach for prediction of presence of opioid peptides in wheat. J. Funct. Foods 2018, 41, 34–40. [Google Scholar] [CrossRef]
  31. Pooja, K.; Rani, S.; Prakash, B. In silico approaches towards the exploration of rice bran proteins-derived angiotensin-I-converting enzyme inhibitory peptides. Int. J. Food Prop. 2017, 20, 2178–2191. [Google Scholar] [CrossRef]
  32. Minkiewicz, P.; Iwaniak, A.; Darewicz, M. BIOPEP-UWM Database of Bioactive Peptides: Current Opportunities. Int. J. Mol. Sci. 2019, 20, 5978. [Google Scholar] [CrossRef]
  33. Chi, C.F.; Wang, B.; Wang, Y.M.; Zhang, B.; Deng, S.G. Isolation and characterization of three antioxidant peptides from protein hydrolysate of bluefin leatherjacket (Navodon septentrionalis) heads. J. Funct. Foods 2015, 12, 1–10. [Google Scholar] [CrossRef]
  34. Harvian, Z.A.; Ningrum, A.; Anggrahini, S.; Setyaningsih, W. In silico Approach in Evaluation of Jack Bean (Canavalia ensiformis) Canavalin Protein as Precursors of Bioactive Peptides with Dual Antioxidant and Angiotensin I-Converting Enzyme Inhibitor. Mater. Sci. Forum 2019, 948, 85–94. [Google Scholar] [CrossRef]
  35. Acquah, C.; Stefano, E.D.; Udenigwe, C.C. Role of hydrophobicity in food peptide functionality and bioactivity. J. Food Bioact. 2018, 4, 88–98. [Google Scholar] [CrossRef]
  36. Zou, T.B.; He, T.P.; Li, H.B.; Tang, H.W.; Xia, E.Q. The Structure-Activity Relationship of the Antioxidant Peptides from Natural Proteins. Molecules 2016, 21, 72. [Google Scholar] [CrossRef]
  37. Daskaya-Dikmen, C.; Yucetepe, A.; Karbancioglu-Guler, F.; Daskaya, H.; Ozcelik, B. Angiotensin-I-Converting Enzyme (ACE)-Inhibitory Peptides from Plants. Nutrients 2017, 9, 316. [Google Scholar] [CrossRef]
  38. Kartal, C.; Kaplan Türköz, B.; Otles, S. Prediction, identification and evaluation of bioactive peptides from tomato seed proteins using in silico approach. J. Food Meas. Charact. 2020, 14, 1865–1883. [Google Scholar] [CrossRef]
  39. Ding, L.; Ma, R.; You, H.; Li, J.; Ge, Q.; Yu, Z.; Wang, L. Identification and characterization of dipeptidyl peptidase IV inhibitory peptides from wheat gluten proteins. J. Cereal Sci. 2022, 103, 103396. [Google Scholar] [CrossRef]
  40. Dugardin, C.; Cudennec, B.; Tourret, M.; Caron, J.; Guerin-Deremaux, L.; Behra-Miellet, J.; Lefranc-Millot, C.; Ravallec, R. Explorative Screening of Bioactivities Generated by Plant-Based Proteins after In vitro Static Gastrointestinal Digestion. Nutrients 2020, 12, 3746. [Google Scholar] [CrossRef]
  41. Wu, J.; Aluko, R.E.; Nakai, S. Structural Requirements of Angiotensin I-Converting Enzyme Inhibitory Peptides: Quantitative Structure−Activity Relationship Study of Di- and Tripeptides. J. Agric. Food Chem. 2006, 54, 732–738. [Google Scholar] [CrossRef]
  42. Iwaniak, A.; Minkiewicz, P.; Pliszka, M.; Mogut, D.; Darewicz, M. Characteristics of Biopeptides Released In silico from Collagens Using Quantitative Parameters. Foods 2020, 9, 965. [Google Scholar] [CrossRef] [PubMed]
  43. Zhu, Q.; Chen, X.; Wu, J.; Zhou, Y.; Qian, Y.; Fang, M.; Xie, J.; Wei, D. Dipeptidyl peptidase IV inhibitory peptides from Chlorella vulgaris: In silico gastrointestinal hydrolysis and molecular mechanism. Eur. Food Res. Technol. 2017, 243, 1739–1748. [Google Scholar] [CrossRef]
  44. Cermeno, M.; Stack, J.; Tobin, P.R.; O’Keeffe, M.B.; Harnedy, P.A.; Stengel, D.B.; FitzGerald, R.J. Peptide identification from a Porphyra dioica protein hydrolysate with antioxidant, angiotensin converting enzyme and dipeptidyl peptidase IV inhibitory activities. Food Funct. 2019, 10, 3421–3429. [Google Scholar] [CrossRef]
  45. 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]
  46. Gupta, S.; Kapoor, P.; Chaudhary, K.; Gautam, A.; Kumar, R.; Open Source Drug Discovery, C.; Raghava, G.P. In silico approach for predicting toxicity of peptides and proteins. PLoS ONE 2013, 8, e73957. [Google Scholar] [CrossRef]
  47. Gupta, S.; Kapoor, P.; Chaudhary, K.; Gautam, A.; Kumar, R.; Raghava, G.P. Peptide toxicity prediction. Methods Mol. Biol. 2015, 1268, 143–157. [Google Scholar] [CrossRef]
  48. Gao, J.; Li, T.; Chen, D.; Gu, H.; Mao, X. Identification and molecular docking of antioxidant peptides from hemp seed protein hydrolysates. LWT 2021, 147, 111453. [Google Scholar] [CrossRef]
  49. Ibrahim, M.A.; Bester, M.J.; Neitz, A.W.; Gaspar, A.R.M. Tuber Storage Proteins as Potential Precursors of Bioactive Peptides: An In silico Analysis. Int. J. Pept. Res. Ther. 2018, 25, 437–446. [Google Scholar] [CrossRef]
  50. Peredo-Lovillo, A.; Hernández-Mendoza, A.; Vallejo-Cordoba, B.; Romero-Luna, H.E. Conventional and in silico approaches to select promising food-derived bioactive peptides: A review. Food Chem. X 2022, 13, 100183. [Google Scholar] [CrossRef] [PubMed]
  51. Sun, X.; Acquah, C.; Aluko, R.E.; Udenigwe, C.C. Considering food matrix and gastrointestinal effects in enhancing bioactive peptide absorption and bioavailability. J. Funct. Foods 2020, 64, 103680. [Google Scholar] [CrossRef]
  52. Ji, D.; Xu, M.; Udenigwe, C.C.; Agyei, D. Physicochemical characterisation, molecular docking, and drug-likeness evaluation of hypotensive peptides encrypted in flaxseed proteome. Curr. Res. Food Sci. 2020, 3, 41–50. [Google Scholar] [CrossRef]
  53. Daina, A.; Michielin, O.; Zoete, V. SwissADME: A free web tool to evaluate pharmacokinetics, drug-likeness and medicinal chemistry friendliness of small molecules. Sci. Rep. 2017, 7, 42717. [Google Scholar] [CrossRef]
  54. Daina, A.; Zoete, V. A BOILED-Egg To Predict Gastrointestinal Absorption and Brain Penetration of Small Molecules. ChemMedChem 2016, 11, 1117–1121. [Google Scholar] [CrossRef] [PubMed]
  55. Aluko, R.E.; Monu, E. Functional and Bioactive Properties of Quinoa Seed Protein Hydrolysates. J. Food Sci. 2003, 68, 1254–1258. [Google Scholar] [CrossRef]
  56. Mohan, A.; Udenigwe, C.C. Towards the design of hypolipidaemic peptides: Deoxycholate binding affinity of hydrophobic peptide aggregates of casein plastein. J. Funct. Foods 2015, 18, 129–136. [Google Scholar] [CrossRef]
  57. 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]
Figure 1. DPPH radical scavenging activity of sophia and camelina protein hydrolysates and their protein isolates (PI). FL, Flavourzyme; AL, Alcalase; and AL + FL, Combination of Alcalase and Flavourzyme.
Figure 1. DPPH radical scavenging activity of sophia and camelina protein hydrolysates and their protein isolates (PI). FL, Flavourzyme; AL, Alcalase; and AL + FL, Combination of Alcalase and Flavourzyme.
Plants 12 03575 g001
Figure 2. ABTS radical scavenging of sophia and camelina protein hydrolysates and their protein isolates (PI). FL, Flavourzyme; AL, Alcalase; and AL + FL, Combination of Alcalase and Flavourzyme.
Figure 2. ABTS radical scavenging of sophia and camelina protein hydrolysates and their protein isolates (PI). FL, Flavourzyme; AL, Alcalase; and AL + FL, Combination of Alcalase and Flavourzyme.
Plants 12 03575 g002
Figure 3. Metal chelation activity of sophia and camelina protein hydrolysates and their protein isolates (PI). FL, Flavourzyme; AL, Alcalase; and AL + FL, Combination of Alcalase and Flavourzyme.
Figure 3. Metal chelation activity of sophia and camelina protein hydrolysates and their protein isolates (PI). FL, Flavourzyme; AL, Alcalase; and AL + FL, Combination of Alcalase and Flavourzyme.
Plants 12 03575 g003
Figure 4. Bioavailability radar of sophia-derived bioactive peptides and camelina-derived bioactive peptides and inhibitor drug (captopril) based on physicochemical indices ideal for oral bioavailability. INSOLU, Insolubility: 0 < Log S (ESOL) < 6; INSATU, Instauration: 0.25 < Fraction Csp3 < 1; FLEX, Flexibility: 0 < Number of rotatable bonds < 9; LIPO, Lipophilicity: −0.7 < XLOGP3 < +5.0; SIZE, Molecular size: 150 g/mol < mol. wt. < 500 g/mol; POLAR, Polarity: 20 Å2 < TPSA <130 Å2. The colored region represents the appropriate physicochemical parameters for oral bioavailability.
Figure 4. Bioavailability radar of sophia-derived bioactive peptides and camelina-derived bioactive peptides and inhibitor drug (captopril) based on physicochemical indices ideal for oral bioavailability. INSOLU, Insolubility: 0 < Log S (ESOL) < 6; INSATU, Instauration: 0.25 < Fraction Csp3 < 1; FLEX, Flexibility: 0 < Number of rotatable bonds < 9; LIPO, Lipophilicity: −0.7 < XLOGP3 < +5.0; SIZE, Molecular size: 150 g/mol < mol. wt. < 500 g/mol; POLAR, Polarity: 20 Å2 < TPSA <130 Å2. The colored region represents the appropriate physicochemical parameters for oral bioavailability.
Plants 12 03575 g004aPlants 12 03575 g004bPlants 12 03575 g004cPlants 12 03575 g004d
Table 1. Oxygen radical absorbance capacity (ORAC) and hydroxyl radical scavenging activity of camelina/sophia protein hydrolysates.
Table 1. Oxygen radical absorbance capacity (ORAC) and hydroxyl radical scavenging activity of camelina/sophia protein hydrolysates.
SamplesProtein HydrolysatesHydroxyl Radical Scavenging ActivityOxygen Radical Absorbance Capacity
(µm of Histidine/mg Protein)(µm of Trolox/mg Protein)
CamelinaAlcalase3.21 ± 0.03 a1.13 ± 0.05 a
Flavourzyme2.50 ± 0.05 b0.62 ± 0.04 b
Alcalase + Flavourzyme2.71 ± 0.16 c0.93 ± 0.03 c
SophiaAlcalase3.13 ± 0.08 a1.03 ± 0.05 a
Flavourzyme1.91 ± 0.05 b0.51 ± 0.01 b
Alcalase + Flavourzyme2.63 ± 0.16 c0.93 ± 0.04 c
Data are presented as mean ± SD (n = 3). Results with identical lowercase letters in the column do not exhibit significant differences at p < 0.05.
Table 2. The inhibition against human LDL cholesterol oxidation by carnosine and sophia/camelina protein hydrolysates.
Table 2. The inhibition against human LDL cholesterol oxidation by carnosine and sophia/camelina protein hydrolysates.
Protein HydrolysatesLDL Oxidation Inhibition (%)
SophiaCamelina
Carnosine78.59 ± 0.11 a78.59 ± 0.11 a
Alcalase73.66 ± 1.19 b79.97 ± 0.29 b
Flavourzyme39.04 ± 0.28 c45.94 ± 0.44 c
Alcalase + Flavourzyme72.39 ± 0.59 b62.68 ± 0.87 d
Data are presented as mean ± SD (n = 3). Results with identical lowercase letters in the column do not exhibit significant differences at p < 0.05.
Table 3. Inhibition of peroxyl and hydroxyl radical-induced supercoiled DNA strand scission by sophia/camelina protein hydrolysates.
Table 3. Inhibition of peroxyl and hydroxyl radical-induced supercoiled DNA strand scission by sophia/camelina protein hydrolysates.
SamplesProtein HydrolysatesDNA Scission Inhibition (%)
Hydroxyl RadicalPeroxyl Radical
SophiaCarnosine30.44 ± 1.36 d22.63 ± 1.15 d
Alcalase53.01 ± 0.50 a87.10 ± 2.35 a
Flavourzyme19.85 ± 0.87 b62.83 ± 1.37 b
Alcalase + Flavourzyme46.35 ± 0.39 c86.52 ± 0.32 a
CamelinaAlcalase59.62 ± 2.21 a88.77 ± 0.75 a
Flavourzyme28.25 ± 0.73 b75.82 ± 1.55 b
Alcalase + Flavourzyme47.02 ± 0.35 c80.07 ± 1.36 c
Data are presented as mean ± SD (n = 3). Results with identical lowercase letters in the column do not exhibit significant differences at p < 0.05.
Table 4. Amino acid composition of sophia/camelina protein hydrolysates (Alcalase) (g/100 g).
Table 4. Amino acid composition of sophia/camelina protein hydrolysates (Alcalase) (g/100 g).
Amino AcidSophia HydrolysatesCamelina Hydrolyses [9]WHO/FAO
Essential amino acid (EAA)
Histidine (His)0.92.11.6
Isoleucine (Ile)1.42.91.3
Leucine (Leu)2.24.91.9
Lysine (Lys)1.42.31.6
Methionine (Met)0.61.51.7
Phenylalanine (Phe)1.53.2
Threonine (Thr)1.32.80.9
Tryptophan (Trp)2.11.30.5
Valine (Val)1.63.81.3
Non-EAA (NEEA)
Alanine (Ala)1.12.8
Arginine (Arg)2.87.0
Aspartic acid +Asparagine (Asp + Asn)3.06.4
Cystine (Cys)0.81.4
Glutamic acid + Glutamine (Glu + Gln)5.812.7
Glycine (Gly)1.83.2
Proline (Pro)1.53.4
Serine (Ser)1.33.0
Tyrosine (Tyr)1.22.4
Phenylalanine + Tyrosine (Phe + Tyr)2.73.31.9
Table 5. Ranking of the predicted bioactive peptide sequences prepared sophia and camelina protein obtained by PepRank.
Table 5. Ranking of the predicted bioactive peptide sequences prepared sophia and camelina protein obtained by PepRank.
SamplesPepRankPeptide Sequence
Sophia hydrolysates0.96FGFGPGL
0.89SSTSGPAFNAGRSIWLPGWL
0.85CAYDVAPGGLL
0.85SLCGIPPL
0.83PMITGFM
0.82FVPVTGLWM
0.82WYTICICIL
0.82RAPWLEPL
0.81LGMLPGL
Camelina hydrolysates0.96GPPSGGGGGGGGGGGGGK
0.94IDLFFVFL
0.94AAMGGFPGGGGGAHALGVL
0.92LNPCFTGGPLM
0.92NGGGGGGGGGGGPPKMVL
0.91PPPPGAL
0.9GGGGGGGFGGGAGGGLGGGGGL
0.9GGSPGIGGGL
0.89DIPPPRGPL
0.89GGGIGGGGGGGGGFGGGSGSGGGAGGGAGGGL
0.89LLGNGIGSGGGHGGKGGRVCY
0.88FGGGNLPAFVL
0.88AWPDKNPFFPSDPY
0.88FVPPFNPY
0.86DFSIFSPL
0.85FGGGNIPAFVL
0.84GGGGGGGGPPAMSM
0.84GLDPPDLPM
0.83VIGPGLGRDPFLL
0.83GGGAGGGGGL
0.83VFGSGLLGAFL
0.83FNAPIYL
0.83GPFGVIRPPL
0.82IHPIPPL
0.82AGAATGGFL
0.82NSGGGGGGWKGGGGQGGGWKGGGGQ
0.81QFQWIEFK
0.81FHWDLPQ
0.81FGGYAPGILSPSPAML
0.81FSFSPTVFDMILK
0.81FGWDKDL
0.8SFPLPEL
Table 6. In silico predictions of potential bioactive peptides from sophia and camelina protein hydrolysate.
Table 6. In silico predictions of potential bioactive peptides from sophia and camelina protein hydrolysate.
SamplesPeptide SequenceBioactive Segment Sequence
AntioxidativeACE-Inhibitory PeptidesDPP IV-Inhibitory Peptides
Peptides
Sophia hydrolysatesFGFGPGLGFGPGLPGL, GP, GF, GL, FG, PGGP, GL, GF, PG
SSTSGPAFNAGRSIWLPGWL-GPA, LPG, P, IW, GW, AF, AG, GR, SG, PG, WL, SGP, ST, LPGP, PA, LP, GPA, WL, AF, AG, FN, GW, IW, NA, PG, SI, TS
CAYDVAPGGLLAYVAP, AY, AP, GL, GG, PGVA, AP, LL, APG, GL, AY, GG, PG, YD
SLCGIPPL-IPP, PL, IP, GI, PP, PPLPP, IP, SL, PL, PPL, GI
PMITGFM-GF, TGGF, MI, PM, TG
FVPVTGLWMLWM, LWLW, VP, GL, TGVP, GL, WM, LW, PV, TG, VT
WYTICICILWY, WYTFVP, WMWY, YT, IL, TI
RAPWLEPLPWL, PWPL, AP, RA, WLAP, RA, WP, PL, WL, PW
LGMLPGL-PGL, LPG, GM, GL, LG, PG, LPLP, GL, ML, PG
Camelina hydrolysatesGPPSGGGGGGGGGGGGGKGPPGP, GK, GG, SG, GPP, PPGP, PP, GG, PS
IDLFFVFL-VF, LF, FFFL, VF, FF
AAMGGFPGGGGGAHALGVLAH, GAHFP, AA, GF, GA, GV, MG, GG, LG,HA, FP, GA, AL, AA, AH, GF, GG,
PG, AH, LGVGV, MG, PG
LNPCFTGGPLM-LNP, GPL, GP, PL, GG, TG, CF, LNGP, NP, PL, GG, LM, LN, TG
NGGGGGGGGGGGPPKMVLGPPGP, GG, NG, PPK, GPP, PPGP, PP, GG, MV, NG, PK, VL
PPPPGAL-GA, PG, PP, PPPPPPP, PP, GA, AL, PPG, PG
GGGGGGGFGGGAGGGLGGGGGL-GF, GA, GL, AG, FG, GG, LG, FGGGA, GL, AG, GF, GG
GGSPGIGGGL-IG, GI, GL, GS, GG, PGSP, GL, GG, GI, PG
DIPPPRGPL-IPP, PR, GPL, GP, PL, IP, PP, PPP, RGGP, PP, IP, PL, RG
GGGIGGGGGGGGGFGGGSGS-GF, IG, GA, GL, AG, FG, GS,GA, GL, AG, GF, GG, GI
-GGGAGGGAGGGLGG, SG, FGG
LLGNGIGSGGGHGGKGGRVCY-IG, GI, GH, GR, KG, GS, GK, HG, GG,LL, GG, GH, GI, KG, NG
SG, LG, NG, GHG
FGGGNLPAFVL-AF, FG, GG, FGG, LPPA, LP, AF, GG, NL, VL
AWPDKNPFFPSDPYAWFP, AW, FFWP, FP, NP, AW, DP, PF, PS, PY, FF
FVPPFNPY-VPP, VP, PP, FVPPP, VP, NP, FN, PF, PY
DFSIFSPL-PL, IF, DFSP, PL, SI
FGGGNIPAFVL-IPA, IP, AF, FG, GG, FGGPA, IPA, IP, AF, GG, VL
GGGGGGGGPPAMSMGPPGP, GG, GPP, PPGP, PP, PA, GG
GLDPPDLPMLPMDLP, GL, PP, LPPP, LP, GL, DP, PM
VIGPGLGRDPFLL-PGL, GP, IG, GL, GR, LG, PGGP, LL, FL, GL, DP, PF, PG, VI
VFGSGLLGAFL-VF, AF, GA, GL, FG, GS, SG, LG, AFLGA, GL, AG, GG
FNAPIYLIYIY, YL, APAP, FN, NA, PI, YL
GPFGVIRPPLIRGP, PL, IRP, RP, FG, GV, PP, RPPGP, PP, RP, PL, PPL, GV, IR, PF, VI
IHPIPPL-IPP, PL, IP, PP, HP, PPLPP, IP, HP, PL, PPL, IH, PI
AGAATGGFLGAAAA, GF, GA, AG, GG, TGGA, FL, AA, AG, AT, GF, GG, TG
NSGGGGGGWKGGGGQGG-GW, KG, GQ, GG, QG, SGWK, GG, GW, KG, QG
-GWKGGGGQ
QFQWIEFK-IE, FQ, EFWI, FQ, QF, QW
FHWDLPQ-DLP, PQ, LPLP, HW, PQ, WD, LPQ
FGGYAPGILSPSPAML-GGY, GY, LSP, AP, YA, GI, FG, GG,AP, PA, APG, SP, GG, GI, GY, IL, ML,
PG, IL, FGGPG, PS, YA
FSFSPTVFDMILKLKVF, SF, PT, DM, IPSP, IL, MI, PT, SF, TV, VF
FGWDKDLKDGW, FGGW, WD
SFPLPELEL, PELPLP, FP, PL, SF, LPLP, FP, PL, SF
SGLGGGQGIGGGSGTGM-IG, GI, GM, GL, GS, GQ, GT, GG, QG,GL, GG, GI, QG, TG
SG, LG, TG, GTG
Table 7. Remaining bioactive properties after in silico simulated gastrointestinal digestion of peptides prepared sophia and camelina protein hydrolysate.
Table 7. Remaining bioactive properties after in silico simulated gastrointestinal digestion of peptides prepared sophia and camelina protein hydrolysate.
SamplePeptideResults of Enzyme ActionLocation of Released PeptidesActive Fragment SequenceLocationBioactivity of Identified
Peptide
Sophia hydrolysatesFGFGPGLF-GF-GPGL[1-1], [2-3], [4-7]GF[2-3]ACE inhibitor, DPP IV inhibitor
SSTSGPAFNAGRSIWLPGWLSSTSGPAF-N-AGR-
SIW-L-PGW-L
[1-8], [9-9], [10-12], [13-15], [16-16], [17-19], [20-20]---
CAYDVAPGGLLCAY-DVAPGGL-L[1-3], [4-10], [11-11]---
SLCGIPPLSL-CGIPPL[1-2], [3-8]SL[1-2]DPP IV inhibitor
PMITGFMPM-ITGF-M[1-2], [3-6], [7-7]PM[1-2]DPP IV inhibitor
FVPVTGLWMF-VPVTGL-W-M[1-1], [2-7], [8-8], [9-9]---
WYTICICILW-Y-TICICIL[1-1], [2-2], [3-9]---
RAPWLEPLR-APW-L-EPL[1-1], [2-4], [5-5], [6-8]---
LGMLPGLL-GM-L-PGL[1-1], [2-3], [4-4], [5-7]PGL[5-7]ACE inhibitor
GM[2-3]ACE inhibitor
Camelina hydrolysatesGPPSGGGGGGGGGGGGGK-----
IDLFFVFLIDL-F-F-VF-L[1-3], [4-4], [5-5],
[6-7], [8-8]
VF[6-7]ACE inhibitor; DPP IV inhibitor
AAMGGFPGGGGGAHALGVLAAM-GGF-PGGGGGAH-AL-GVL[1-3], [4-6], [7-14],
[15-16], [17-19]
AL[15-16]DPP IV inhibitor
LNPCFTGGPLML-N-PCF-TGGPL-M[1-1], [2-2], [3-5],
[6-10], [11-11]
---
NGGGGGGGGGGGPPKMVLN-GGGGGGGGGGGPPK -M-VL[1-1], [2-15],
[16-16], [17-18]
VL[17-18]DPP IV inhibitor
PPPPGAL-----
GGGGGGGFGGGAGG-GLGGGGGLGGGGGGGF-GGGAGGGL -GGGGGL[1-8], [9-16], [17-22]---
GGSPGIGGGL-----
DIPPPRGPLDIPPPR-GPL[1-6], [7-9]GPL[7-9]ACE inhibitor
GGGIGGGGGGGGGFGGGSG-SGGGAGGGAGGGLGGGIGGGGGGGGGF -GGGSGSGGGAGGGAGGGL[1-14], [15-32]---
LLGNGIGSGGGHGGKGGRVCYL-L-GN-GIGSGGGH
-GGK-GGR-VCY
[1-1], [2-2], [3-4], [5-12],---
FGGGNLPAFVLF-GGGN-L-PAF-VL[1-1], [2-5], [6-6],
[7-9], [10-11]
VL[10-11]DPP IV inhibitor
AWPDKNPFFPSDPYAW-PDK-N-PF
-F-PSDPY
[1-2], [3-5], [6-6],
[7-8], [9-9], [10-14]
AW[1-2]ACE inhibitor;
DPP IV inhibitor; antioxidative
PF[7-8]DPP IV inhibitor
FVPPFNPYF-VPPF-N-PY[1-1], [2-5], [6-6], [7-8]PY[7-8]DPP IV inhibitor
DFSIFSPLDF-SIF-SPL[1-2], [3-5], [6-8]DF[1-2]ACE inhibitor
FGGGNIPAFVLF-GGGN-IPAF-VL[1-1], [2-5], [6-9], [10-11]VL[10-11]DPP IV inhibitor
GGGGGGGGPPAMSMGGGGGGGGPPAM-SM[1-12], [13-14]---
GLDPPDLPMGL-DPPDL-PM[1-2], [3-7], [8-9]GL[1-2]ACE inhibitor;
DPP IV inhibitor
PM[8-9]DPP IV inhibitor
VIGPGLGRDPFLLVIGPGL-GR-DPF-L-L[1-6], [7-8], [9-11],
[12-12], [13-13]
GR[7-8]ACE inhibitor
GGGAGGGGGL-----
VFGSGLLGAFLVF-GSGL-L-GAF-L[1-2], [3-6], [7-7],
[8-10], [11-11]
VF[1-2]ACE inhibitor;
DPP IV inhibitor
FNAPIYLF-N-APIY-L1-1], [2-2], [3-6], [7-7]---
GPFGVIRPPLGPF-GVIR-PPL[1-3], [4-7], [8-10]PPL[8-10]ACE inhibitor;
DPP IV inhibitor
IHPIPPLIH-PIPPL[1-2], [3-7]IH[1-2]DPP IV inhibitor
AGAATGGFLAGAATGGF-L[1-8], [9-9]---
NSGGGGGGWKGGGGQGG-GWKGGGGQN- SGGGGGGW-K -GGGGQGGGW-K-GGGGQ[1-1], [2-9], [10-10],
[11-19], [20-20], [21-25]
---
QFQWIEFKQF-QW-IEF-K[1-2], [3-4], [5-7], [8-8]QF[1-2]DPP IV inhibitor
QW[3-4]DPP IV inhibitor
FHWDLPQF-H-W-DL-PQ[1-1], [2-2], [3-3],
[4-5], [6-7]
PQ[6-7]ACE inhibitor;
DPP IV inhibitor
FGGYAPGILSPSPAMLF-GGY-APGIL-SPSPAM-L[1-1], [2-4], [5-9],
[10-15], [16-16]
GGY[2-4]ACE inhibitor
FSFSPTVFDMILKF-SF-SPTVF-DM-IL-K1-1], [2-3], [4-8],
[9-10], [11-12], [13-13]
SF[2-3]ACE inhibitor;
DPP IV inhibitor
DM[9-10]ACE inhibitor
IL[11-12]ACE inhibitor;
DPP IV inhibitor
FGWDKDLF-GW-DK-DL[1-1], [2-3], [4-5], [6-7]GW[2-3]ACE inhibitor;
DPP IV inhibitor
SFPLPELSF-PL-PEL[1-2], [3-4], [5-7]SF[1-2]ACE inhibitor;
DPP IV inhibitor
PL[3-4]ACE inhibitor;
DPP IV inhibitor
PEL[5-7]antioxidative
Table 8. In silico proteolysis of sophia and camelina protein hydrolysate for release bioactive peptides sequences.
Table 8. In silico proteolysis of sophia and camelina protein hydrolysate for release bioactive peptides sequences.
SamplePeptideActive Fragment
Sequence
DHt (%)AEWActivity
Sophia hydrolysatesFGFGPGLGF33.330.140.14ACE inhibitor
0.140.25DPP IV inhibitor
SLCGIPPLSL14.290.130.17DPP IV inhibitor
PMITGFMPM33.330.140.25DPP IV inhibitor
LGMLPGLPGL, GM500.290.29ACE inhibitor
Camelina hydrolysatesIDLFFVFLVF57.140.130.33ACE inhibitor
0.130.33DPP IV inhibitor
AAMGGFPGGGGGAHALGVLAL22.220.050.06DPP IV inhibitor
NGGGGGGGGGGGPPKMVLVL17.650.060.06DPP IV inhibitor
DIPPPRGPLGPL12.50.110.1ACE inhibitor
FGGGNLPAFVLVL400.090.14DPP IV inhibitor
AWPDKNPFFPSDPYAW, PF38.460.070.33ACE inhibitor
0.140.22DPP IV inhibitor
0.071Antioxidative
FVPPFNPYPY42.850.130.17DPP IV inhibitor
DFSIFSPLDF28.570.130.33ACE inhibitor
FGGGNIPAFVLVL300.090.14DPP IV inhibitor
GLDPPDLPMGL, PM250.110.25ACE inhibitor
0.220.4DPP IV inhibitor
VIGPGLGRDPFLLGR33.330.080.14ACE inhibitor
VFGSGLLGAFLVF400.090.11ACE inhibitor
0.090.17DPP IV inhibitor
GPFGVIRPPLPPL22.220.10.1ACE inhibitor
0.10.11DPP IV inhibitor
IHPIPPLIH16.670.140.14DPP IV inhibitor
QFQWIEFKQF, QW42.860.250.5DPP IV inhibitor
FHWDLPQPQ66.670.140.33ACE inhibitor
0.140.2DPP IV inhibitor
FGGYAPGILSPSPAMLGGY26.670.060.09ACE inhibitor
FSFSPTVFDMILKSF, DM, IL41.670.230.6ACE inhibitor
0.150.29DPP IV inhibitor
FGWDKDLGW500.140.5ACE inhibitor
0.140.5DPP IV inhibitor
SFPLPELPL, SF, PEL33.330.290.4ACE inhibitor
0.290.5DPP IV inhibitor
0.140.5Antioxidative
Table 9. Prediction of toxicity and allergenicity of the potential bioactive fragments of sophia and camelina protein hydrolysate obtained after the simulation in silico digestion.
Table 9. Prediction of toxicity and allergenicity of the potential bioactive fragments of sophia and camelina protein hydrolysate obtained after the simulation in silico digestion.
SamplesActive FragmentHydrophobicityHydrophilicityChargepIMolecularToxinAllergenicity
SequenceWeight (Da)PredictionPrediction
Sophia hydrolysatesGF0.39−1.2505.88222.26Non-ToxinProbable allergen
SL0.14−0.7505.88218.27Non-ToxinProbable non-allergen
PM0.1−0.6505.88246.34Non-ToxinProbable non-allergen
GM0.21−0.6505.88206.28Non-ToxinProbable non-allergen
PGL0.21−0.605.88285.38Non-ToxinProbable non-allergen
Camelina hydrolysatesVF0.57−205.88264.34Non-ToxinNon-allergen
AL0.39−1.1505.88202.27Non-ToxinProbable allergen
VL0.54−1.6505.88230.33Non-ToxinProbable non-allergen
GPL0.21−0.605.88285.38Non-ToxinProbable allergen
AW0.31−1.9505.88275.32Non-ToxinProbable allergen
PF0.27−1.2505.88262.32Non-ToxinProbable allergen
PY−0.03−1.1505.88278.32Non-ToxinProbable allergen
DF−0.050.25−13.8280.29Non-ToxinProbable allergen
GL0.35−0.905.88188.25Non-ToxinProbable allergen
PM0.1−0.6505.88246.34Non-ToxinProbable non-allergen
GR−0.81.5110.11231.27Non-ToxinProbable allergen
PPL0.13−0.605.88325.44Non-ToxinProbable allergen
IH0.16−1.150.57.1268.34Non-ToxinProbable allergen
QF−0.04−1.1505.88293.34Non-ToxinProbable allergen
QW−0.16−1.605.88332.38Non-ToxinProbable allergen
PQ−0.380.105.88243.28Non-ToxinProbable allergen
GGY0.11−0.7705.88295.33Non-ToxinProbable non-allergen
SF0.17−1.105.88252.28Non-ToxinProbable non-allergen
DM−0.230.85−13.8264.31Non-ToxinProbable allergen
IL0.63−1.805.88244.36Non-ToxinProbable non-allergen
GW0.27−1.705.88261.3Non-ToxinProbable non-allergen
PL0.23−0.905.88228.31Non-ToxinProbable allergen
PEL−0.050.4−14357.44Non-ToxinProbable allergen
Table 10. In silico physicochemical properties and absorption, distribution, metabolism, excretion (ADME) profile of sophia and camelina bioactive peptides.
Table 10. In silico physicochemical properties and absorption, distribution, metabolism, excretion (ADME) profile of sophia and camelina bioactive peptides.
SamplesActive Fragment
Sequence
Physicochemical PropertiesLipophilicityDrug LikelinessPharmacokinetics
ROTBHBAHBDTPSA (Ų)ESOLC LogPLipinski
Filter
Bioavailability
Score
GI Absorption
Sophia hydrolysatesCaptopril43196.41−1.140.62Yes (0)0.56High
GF64392.420.32−0.24Yes (0)0.55High
SL754112.650.84−0.86Yes (0)0.55High
PM743103.731.04−0.32Yes (0)0.55High
GM743117.721.32−0.82Yes (0)0.55High
PGL954107.530.91−0.51Yes (0)0.55High
Camelina hydrolysatesVF74392.420.110.6Yes (0)0.55High
AL64392.421.05−0.41Yes (0)0.55High
VL74392.420.50.27Yes (0)0.55High
GPL853112.730.99−0.6Yes (0)0.55High
AW644108.21−0.090.11Yes (0)0.55High
PF64378.4300.33Yes (0)0.55High
PY65498.660.140Yes (0)0.55High
DF864129.720.83−0.59Yes (0)0.56High
GL64392.420.86−0.49Yes (0)0.55High
PM743103.731.04−0.32Yes (0)0.55High
GR855156.822.31−2.16Yes (0)0.55Low
PPL85398.740.090.06Yes (0)0.55High
IH854121.10.2−0.33Yes (0)0.55High
QF954134.511.01−0.68Yes (0)0.55High
QW953151.30.49−0.48Yes (0)0.55Low
PQ754121.522.19−1.6Yes (0)0.55Low
GGY965141.751.14−1.22Yes (0)0.55Low
SF754112.651.03−0.77Yes (0)0.55High
DM964155.021.67−1.2Yes (0)0.11Low
IL84392.420.250.63Yes (0)0.55High
GW644108.210.23−0.21Yes (0)0.55High
PL64378.430.660.04Yes (0)0.55High
PEL1275144.830.92−0.52Yes (0)0.56Low
Abbreviations: ROTB, number of rotatable bonds; HBA, number of hydrogen bond acceptors; HBD, number of hydrogen bond donors; TPSA, topological polar surface area; ESOL, estimate solubility; C LogP, logarithm of compound partition coefficient between n-octanol and water; GI absorption, gastrointestinal absorption.
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Ngo, N.T.T.; Senadheera, T.R.L.; Shahidi, F. Antioxidant Properties and Prediction of Bioactive Peptides Produced from Flixweed (sophia, Descurainis sophia L.) and Camelina (Camelina sativa (L.) Crantz) Seed Meal: Integrated In Vitro and In Silico Studies. Plants 2023, 12, 3575. https://doi.org/10.3390/plants12203575

AMA Style

Ngo NTT, Senadheera TRL, Shahidi F. Antioxidant Properties and Prediction of Bioactive Peptides Produced from Flixweed (sophia, Descurainis sophia L.) and Camelina (Camelina sativa (L.) Crantz) Seed Meal: Integrated In Vitro and In Silico Studies. Plants. 2023; 12(20):3575. https://doi.org/10.3390/plants12203575

Chicago/Turabian Style

Ngo, Na Thi Ty, Tharindu R. L. Senadheera, and Fereidoon Shahidi. 2023. "Antioxidant Properties and Prediction of Bioactive Peptides Produced from Flixweed (sophia, Descurainis sophia L.) and Camelina (Camelina sativa (L.) Crantz) Seed Meal: Integrated In Vitro and In Silico Studies" Plants 12, no. 20: 3575. https://doi.org/10.3390/plants12203575

APA Style

Ngo, N. T. T., Senadheera, T. R. L., & Shahidi, F. (2023). Antioxidant Properties and Prediction of Bioactive Peptides Produced from Flixweed (sophia, Descurainis sophia L.) and Camelina (Camelina sativa (L.) Crantz) Seed Meal: Integrated In Vitro and In Silico Studies. Plants, 12(20), 3575. https://doi.org/10.3390/plants12203575

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

Article Metrics

Back to TopTop