Next Article in Journal
Strategies for Improving Bioavailability, Bioactivity, and Physical-Chemical Behavior of Curcumin
Next Article in Special Issue
Application of Density Measurements for Discrimination and Evaluation of Chemical Composition of Different Types of Mechanically Separated Meat (MSM)
Previous Article in Journal
Maillard and Hydrolytic Reactions in Subcritical Water Extraction of Bioactive Compounds from Licorice
Previous Article in Special Issue
Identification and Quantification of Anthocyanin and Catechin Compounds in Purple Tea Leaves and Flakes
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

In Silico Approach of Glycinin and Conglycinin Chains of Soybean By-Product (Okara) Using Papain and Bromelain

1
Department of Food and Agricultural Product Technology, Faculty of Agricultural Technology, Universitas Gadjah Mada, Flora Street No. 1, Bulaksumur, Yogyakarta 55281, Indonesia
2
Research Centre for Appropriate Technology, National Research and Innovation Agency, KS. Tubun Street No.5, Subang 41213, Indonesia
3
Research Center for Food Technology and Processing, National Research and Innovation Agency, Jogja-Wonosari Street km 31,5 Playen, Gunungkidul, Yogyakarta 55861, Indonesia
4
Study Program of Chemistry, Department of Chemistry Education, Faculty of Mathematics and Science Education, Universitas Pendidikan Indonesia, Bandung 40154, Indonesia
*
Author to whom correspondence should be addressed.
Molecules 2022, 27(20), 6855; https://doi.org/10.3390/molecules27206855
Submission received: 29 September 2022 / Revised: 6 October 2022 / Accepted: 10 October 2022 / Published: 13 October 2022
(This article belongs to the Special Issue Food Chemistry: Food Quality and New Analytical Approaches)

Abstract

:
This study explores utilization of a sustainable soybean by-product (okara) based on in silico approach. In silico approaches, as well as the BIOPEP database, PeptideRanker database, Peptide Calculator database (Pepcalc), ToxinPred database, and AllerTop database, were employed to evaluate the potential of glycinin and conglycinin derived peptides as a potential source of bioactive peptides. These major protein precursors have been found as protein in okara as a soybean by-product. Furthermore, primary structure, biological potential, and physicochemical, sensory, and allergenic characteristics of the theoretically released antioxidant peptides were predicted in this research. Glycinin and α subunits of β-conglycinin were selected as potential precursors of bioactive peptides based on in silico analysis. The most notable among these are antioxidant peptides. First, the potential of protein precursors for releasing bioactive peptides was evaluated by determining the frequency of occurrence of fragments with a given activity. Through the BIOPEP database analysis, there are several antioxidant bioactive peptides in glycinin and β and α subunits of β-conglycinin sequences. Then, an in silico proteolysis using selected enzymes (papain, bromelain) to obtain antioxidant peptides was investigated and then analyzed using PeptideRanker and Pepcalc. Allergenic analysis using the AllerTop revealed that all in silico proteolysis-derived antioxidant peptides are probably nonallergenic peptides. We also performed molecular docking against MPO (myeloperoxidases) for this peptide. Overall, the present study highlights that glycinin and β and α subunits of β-conglycinin could be promising precursors of bioactive peptides that have an antioxidant peptide for developing several applications.

1. Introduction

Almost all the world’s soybean production is used directly for human food [1]. Soybean has been commercialized in nature and soy-derived products include textured soy, soymilk, tofu, and fermented products (miso, shoyo, and tempeh). Soybean is a very important commodity in several areas as a protein source [1,2]. Okara is first developed by the production of soymilk. Soymilk can be made from whole SB or fat soy flour. Its production usually consists of five main stages: (1) SB washing to remove impurities; (2) SB soaking/hydration for 12 h at 25 °C, then draining and rinsing with water; (3) cooking at 98 °C for 5 min, with the aim of both sterilizing and improving aroma and nutritional value via the inactivation of trypsin inhibitors; (4) grinding in a blender with distilled water (1:10 ratio w/v SBs/water) for the preparation of a slurry; (5) separation of the slurry into soymilk (water-soluble SB extract) and okara by mechanical means (usually filtration) [3]. The production of soymilk made from soybean results in an insoluble by-product called okara with limited market value which is usually used as animal feed [4]. Okara is the waste produced during the production of soy milk, tofu, and soybean protein isolate in the food industry [5]. Okara has a high nutritive value due to its high-quality protein, fat, carbohydrates, and fiber that can still be utilized further. Concerning the proteins initially present in soybeans, approximately 23% are retained in the okara [3]. These proteins have a high nutritive value and superior protein efficiency ratio, indicating that okara is a potential source of low-cost vegetable protein for human consumption [6]. Besides protein, soybean and okara also contain bioactive compounds such as isoflavones [7].
In silico analysis is a common approach of bioinformatics that is widely used in analyzing the amino acid sequence, protein domain, and protein structure, through a computational approach [6,8,9]. In silico methods of peptide analyses could include different approaches such as homology modeling, molecular dynamics, protein docking, and PPI targeting. Structural characterization of the peptides could be carried out by x-ray crystallography, NMR spectroscopy, and cryo-electron microscopy instruments. The obtained structural data are stored and available in structural deposition databases like the Protein Data Bank (PDB). The advantages of computational in silico methods over empirical methods are their low cost, faster procedure speed, simple process, and reliability in target PPIs using peptides [10]. It is an advancement for the simulation of biological experiments in the early stage of research, especially for food proteins, to predict the proteins and peptide sequences, as well as their bioactivities [11,12]. In addition, potentially harmful substances that may be present within the proteins or peptides, such as toxins or allergens, can be discovered by using an in silico approach. Simulation of theoretical enzymatic hydrolysis, prediction of peptide sequences, and their potential biological activities can be conducted using the BIOPEP-UWM database, a freely accessible bioinformatics database that is widely used in research on biologically active peptides [13]. We already successfully applied this in silico approach analysis to identify bioactive peptides from several food sources such as yellowfin, tilapia, and jack bean [14,15].
Nowadays, the exploratory study of bioactive peptides from food proteins is gaining momentum, especially in finding functional properties that are beneficial for human health [16,17]. The bioactive peptides are commonly generated when their parent proteins are hydrolyzed by endogenous or exogenous enzymes into smaller fragments with their functional properties. Bioactive peptides derived from food proteins, containing 2–20 amino acid residues, have been reported to possess immunomodulatory, anti-cancer, anti-hypertensive, antioxidant, anti-proliferative, hypocholesterolemic, metal chelating, and anti-inflammatory activities [16,17]. For example, some health benefits of soybean bioactive peptides have been reported in various studies, including antioxidant activities [18]. They are potential food functional ingredients and can be used for different purposes, such as aroma enhancers, food additives, nutritional supplements, or therapeutic drugs. Much of the interest in plant-derived bioactive peptides present in foodstuffs lies in their potential pharmaceutical and/or nutraceutical benefits [19]. Well-known examples are soymetide (from soybean) and oryzatensin (from rice). In both cases, the bioactive peptide is a degradation product of an abundant protein present in the consumed part of the plant. The present bioinformatics or in silico approach has identified large numbers of proteins in food-producing plants from which bioactive peptides could be generated [19].
The antioxidant properties of soy protein have been ascribed to the cooperative effect of several properties including their ability to scavenge free radicals, act as a metal-ion chelator, oxygen quencher, or hydrogen donor, and the possibility of preventing the penetration of lipid oxidation initiators by forming a membrane around oil droplets [16]. The operational conditions employed in the processing of protein isolates, the type of protease, and the degree of hydrolysis affect the antioxidant and other biological activities [20]. The alternative substitution of synthetic antioxidants with natural ones is gaining interest due to consumer preferences and health concerns associated with the use of synthetic food additives [21]. Additionally, the importance of several of these natural antioxidants in the prevention or control of certain chronic diseases has been increasing. The use of proteins or their hydrolysates for food and/or cosmetic applications presents additional advantages over other antioxidants since they also confer nutritional and functional properties [22].
Glycinin (32 and 20 kDa respectively) and the β and α subunits of β-conglycinin (52 and 68 kDa) are major proteins in okara [23]. In this study, in silico analysis was performed on Glycinin and the β and α subunits of β-conglycinin to determine their effectiveness in releasing potential bioactive peptides with antioxidant activities using two different proteases, papain and bromelain, in BIOPEP-UWM database. We used two different proteases, papain and bromelain, in silico since these two proteases are mainly used and can be released into the bioactive peptides with antioxidant activities [24,25]. Previous research has already measured the physical treatment such as HPH and Alcalase to obtain the antioxidant peptide in Okara [5]. In our study, the antioxidant peptides will be analyzed to elucidate their physicochemical properties, gastrointestinal stability, toxicity, and allergenicity. In this study, we hypothesized that soy glycinin and the β and α subunits of β-conglycinin are proteins from which gastrointestinal-stable, non-toxic, and non-allergenic bioactive peptides can be liberated. The outcome of this study, although preliminary, could add to a better understanding of the potential of protein-rich food wastes, such as the soybean protein by-product okara, as a source of health-promoting peptides.

2. Results and Discussion

2.1. Profile of Potential Biological Activity

The antioxidant biological activities peptides released by glycinin and the β and α subunits of β-conglycinin are reported in Table 1. Hence, the time-saving and more economical computer-simulated method or in silico methods can be used to predict the generation of the bioactive peptides from glycinin and the β and α subunits of β-conglycinin. The biological activities of the glycinin and the β and α subunits of β-conglycinin were evaluated using the BIOPEP analysis tool. It was predicted by the tool that glycinin and the β and α subunits of β-conglycinin have the potential to release peptides with several antioxidant activities (Table 1). The antioxidant activities of bioactive peptides are mainly due to the presence of some aromatic amino acids, and histidine [11].

2.2. In Silico Proteolysis

The potential glycinin and the β and α subunits of β-conglycinin derived peptide sequences displaying the antioxidant peptide were predicted using the BIOPEP analysis tool. The lists of the theoretically released antioxidant peptide from glycinin and the β and α subunits of β-conglycinin are summarized in Table 2. The glycinin and the β and α subunits of β-conglycinin and protein sequences were also subjected to in silico proteolysis using the BIOPEP tool. The supposed antioxidant peptide derived from the hydrolysis of glycinin and the β and α subunits of β-conglycinin by papain and bromelain proteases are summarized and presented in (Table 2).
The BIOPEP analysis method was used to predict possible glycinin and the β and α subunits of β-conglycinin-derived peptide sequences expressing the antioxidant activity. Table 2 summarizes the lists of antioxidant peptides released from glycinin as well as the β and α subunits of β-conglycinin using papain and bromelain enzymes.
According to Table 2, the majority of peptides produced by hydrolysis with papain and bromelain are dipeptides. This finding is consistent with prior studies employing these two distinct enzymes but with a different raw material [25,26,27]. The antioxidant peptides produced by the enzyme combination were primarily dipeptides and tripeptides with antioxidant activity, as predicted by the BIOPEP tool. The combination of papain and bromelain can generate several antioxidant dipeptides and tripeptides.

2.3. Physiochemical Properties of Peptides

The antioxidant activity of peptides is affected by their structural characteristics, especially their molecular masses, amino acid composition, and peptide sequence [28]. Peptides composed of two to three amino acids have previously been shown to have a high antioxidant potential [29,30]. The calculated molecular weight of peptides derived from glycinin and conglycinin can be seen in Table 2.
Furthermore, several properties of peptides are mentioned in Table 3. Peptides with low molecular weight are generally non-toxic and tend to be non-allergenic. The solubility, net charge, and iso-electric point are important physical parameters that are supposed to be considered in the design of novel drugs since these play a role in the distribution in the human body and in targeting specific bacteria, cells, proteins, or viruses. Consequently, further physical chemistry parameters of the peptides proposed were determined by the web server INNOVAGEN’s peptide calculator (PEPCALC) and the result is shown in Table 3. Peptides IR and EL were soluble in aqueous media. Moreover, identified peptides’ charges are varied (Table 3) at pH 7, which is a desirable property for MPO inhibitors (|charge| < −11), which is the net charge of MPO to avoid a possible cytotoxic effect.

2.4. Peptide Ranking of Antioxidant Peptides

The identified antioxidant peptide profile released from glycinin and conglycinin was subjected to activity prediction by in silico methods. PeptideRanker has been used to determine the antioxidant peptides’ potential. The PeptideRanker server can rank the peptide sets according to structure–function patterns [29]. The peptide generated from the papain and bromelain hydrolysis with the maximum peptide ranking scores can be seen in Table 4. The antioxidant dipeptides PHF and YYL showed the maximum peptide rank scores of 0.93 and 0.59, respectively (Table 4). Papain and bromelain mostly generate several active antioxidant peptides [2,14]. The toxicity of peptides derived from glycinin and conglycinin is one of the factors that can inhibit peptides’ development as functional food ingredients [31,32]. Based on Table 4, the sensory evaluation of the peptides that can be released from glycinin and conglycinin by papain and bromelain enzymes are mostly bitter and umami peptides. On the other hand, allergenic effect characteristics showed that most peptides can also cause allergenic effects. In silico toxicity and physicochemical predictions not only save time and money, but also facilitate the design of better therapeutic peptides, with low toxicity without compromising functionality [33,34]. Based on the result shown in Table 4, all the peptides derived from glycinin and conglycinin are not toxic.

2.5. Molecular Docking of the Peptide with the Antioxidant Enzyme In Silico (Binding Affinity and Interaction)

Molecular docking is an in silico method of mimicking the interaction model between receptors and small molecules [35]. Multiple molecular docking techniques were employed to choose MPO inhibitors with antioxidant activity, demonstrating the efficacy of this approach. Consequently, molecular docking is a helpful method to identify antioxidants and reveal underlying molecular mechanisms. Based on the energy binding affinity obtained through Pyrix, all the peptides derived from glycinin and conglycinin successfully docked to MPO spontaneously. Peptides derived from glycinin and conglycinin exhibit high affinity against MPO, and molecular interaction between peptides and MPO involved hydrogen bonds, van der walls bond, and several protein-ligand interactions (Figure 1). The presence of hydrogen bond interactions shows the stability of the complexes formed between peptide and protein. One of the most significant interactions between molecules is the hydrogen bond, which stabilizes the three-dimensional structure of proteins and nucleic acids. Moreover, all the peptides derived from glycinin and conglycinin show interaction with the important residue of MPO. The present finding confirmed the possibilities for peptides, especially dipeptides and tripeptides derived from glycinin and conglycinin, to act as antioxidants, hindering the substrate (such as H2O2) from accessing the heme-containing core of MPO. Compared to small compounds and proteins, peptides are much more flexible and do not have stable conformation before binding to a receptor [36]. To efficiently consider peptide flexibility in this study, we employed energy affinity binding dock the peptides derived from glycinin and conglycinin to MPO. The molecular docking result was expressed by (-) energy binding activity (Figure 2). The value of (-) energy indicates the bonding affinity between MPO and peptides, and a higher value means a more favorable combination [37].
The primary enzyme responsible for the production of ROS in vivo is MPO because heme peroxidase is known to occur in neutrophils and monocytes [34] MPO can produce highly active ROS (like HOCl) under peroxide H2O2 and halide conditions, leading to cell and tissue damage [34]. Consequently, molecular docking was performed to assess whether the various peptides were effective. According to Figure 2, the energy binding activity of peptides suggested that the peptide might spontaneously bind to the MPO.
Inhibitory peptides have the potential to inactivate an enzyme by occupying its active sites or obstructing access to the active site cavity of the enzyme [34]. As can be seen in Figure 3, the various peptides’ MPO docking postures were comparable to the earlier findings of antioxidant peptides from hemp seed, and both prevented entrance to the MPO active cavity, demonstrating their efficacy as antioxidants [33]. Hydrogen bonds are one sort of interaction between MPO and peptides. By preventing substances (such H2O2) from accessing the heme-containing core of MPO, peptides that attach to these locations can act as antioxidants in vivo. Therefore, like an earlier finding, their inhibitory effects on MPO may contribute to the antioxidant action of hemp seed protein-derived peptides in vivo [34].

3. Material and Methods

3.1. Profiles of the Potential Biological Activity

The BIOPEP analysis was used to evaluate the biological potential of glycinin and the β and α subunits of β-conglycinin (http://www.uwm.edu.pl/biochemia/index.php/en/biopep (accessed on 1 August 2022)). The sequences of glycinin and the β and α subunits of β-conglycinin obtained from the NCBI database were evaluated for potential biological activity profiles. The potential antioxidant peptide sequences of glycinin and the β and α subunits of β-conglycinin were also investigated. The frequency of occurrence of bioactive fragments in glycinin, as well as the β and α subunits of β-conglycinin, was also estimated. The sequences of glycinin and the β and α subunits of β-conglycinin were also subjected to in silico proteolysis to predict the theoretical peptide sequences cleaved by the enzymes papain and bromelain. Finally, a list of potential antioxidant peptides was compiled for further analysis.

3.2. Peptide Ranking

PeptideRanker (http://bioware.ucd.ie/compass/biowareweb/ (accessed on 1 August 2022)) was used to predict the ability of glycinin and the β and α subunits of β-conglycinin to generate antioxidant peptides. The peptide score of the chosen glycinin and conglycinin-derived antioxidant peptides was calculated using the PeptideRanker program. PeptideRanker assigned each peptide a score between 0 and 1. The highest score (1) indicated the most active peptides, while the lowest score (0) indicated the least active peptides.

3.3. Sensory Characteristics Prediction

Peptides and amino acids are examples of molecules that can modify the flavor of food. The BIOPEP analysis was used to predict the occurrence frequency of sensory features in hydrolyzed glycinin and the β and α subunits of β-conglycinin. The occurrence frequencies of the various sensory characteristics were also predicted for the generated peptides from glycinin and the β and α subunits of β-conglycinin using papain and bromelain. Furthermore, the sensory characteristics of the selected antioxidant were also predicted.

3.4. The Physicochemical Characteristics of the Antioxidant Peptides

Online peptide calculators were used to analyzing the potential antioxidant peptides. The online Pepcalc program (http://pepcalc.com/ (accessed on 1 August 2022)) was used to calculate the theoretical molecular weight, isoelectric point, peptide charge at pH 7, estimated solubility, and extinction coefficient of the screened antioxidative derived peptides from glycinin and β and α subunits of β-conglycinin peptides.

3.5. Allergenicity Prediction

One of the key problems in the development of soybean protein-derived functional food components is the allergenicity of the peptides. Therefore, AllerTop (https://ddg-pharmfac.net/AllerTOP/index.html (accessed on 1 August 2022)) was used to predict the allergenicity of the discovered glycinin and conglycinin-derived antioxidant peptides.

3.6. Molecular Docking against Peptide

AutoDock Vina- 4.2.6 was used to conducting a molecular docking investigation. The Crystal Structure of Human MPO (PDB ID: 3F9P) was retrieved from the RCSB Protein Data Bank (https://www.rcsb.org/ (accessed on 1 August 2022)), and the peptides derived were retrieved from PubChem. Both protein and peptides were prepared using PyRix. The activity was generated through the global docking method.

3.7. Visualization

The visualization of docking results was carried out through PyMol (3D), and Biovia discovery (3D).

4. Conclusions

We can conclude that glycinin and conglycinin as major protein fractions in okara can be a precursor for several bioactive peptides with antioxidant activity using an in silico approach. The current study showed that using a papain and bromelain-based in silico method with glycinin and conglycinin chains as protein precursors successfully generated antioxidant peptides, the majority in the form of dipeptides and tripeptides. These results can be applied further for the use of a soybean by-product called okara as a precursor to antioxidant peptides in the development of functional food and nutraceuticals. Therefore, the use of okara, a by-product of soybean production, for producing bioactive peptides is in line with the application idea of valorization and can be viewed as a waste management strategy.

Author Contributions

Conceptualization, A.N., A.S., R.K., R.E., D.K. and W.S.; methodology, A.N.; software, A.N.; validation, A.N.; formal analysis, A.N.; investigation, A.N., N.V. and D.W.W.; resources, A.N. and A.S.; data curation, A.N., N.V. and D.W.W.; writing—original draft preparation, A.N., W.S. and H.S.H.M.; writing—review and editing, A.N., W.S. and H.S.H.M.; visualization, A.N.; supervision, A.N. and A.S.; project administration, A.N. and W.S.; funding acquisition, A.N. All authors have read and agreed to the published version of the manuscript.

Funding

The authors wish to thank Kementerian Pendidikan, Kebudayaan, Riset, dan Teknologi, Republik Indonesia under the Project of Penelitian Penugasan (PDKN) Number 3799/UN1/DITLIT/Dit-Lit/PT.01.03/2022.

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.

Sample Availability

Samples of the compounds are not available from the authors.

References

  1. Arunrat, N.; Sereenonchai, S.; Chaowiwat, W.; Wang, C. Climate change impact on major crop yield and water footprint under CMIP6 climate projections in repeated drought and flood areas in Thailand. Sci. Total Environ. 2022, 807, 150741. [Google Scholar] [CrossRef]
  2. Šarčević, H.; Iljkić, D.; Andrijanić, Z.; Sudarić, A.; Gunjača, J.; Varnica, I.; Rastija, M.; Godena, S.; Maričević, M.; Stepinac, D.; et al. Stability of Protein and Oil Content in Soybean across Dry and. Agronomy 2022, 12, 915. [Google Scholar] [CrossRef]
  3. Colletti, A.; Attrovio, A.; Boffa, L.; Mantegna, S.; Cravotto, G. Valorisation of by-products from soybean (Glycine max (L.) Merr.) processing. Molecules 2020, 25, 2129. [Google Scholar] [CrossRef]
  4. Yang, L.C.; Fu, T.J.; Yang, F.C. Biovalorization of soybean residue (okara) via fermentation with Ganoderma lucidum and Lentinus edodes to attain products with high anti-osteoporotic effects. J. Biosci. Bioeng. 2020, 129, 514–518. [Google Scholar] [CrossRef] [PubMed]
  5. Fang, J.; Lu, J.; Zhang, Y.; Wang, J.; Wang, S.; Fan, H.; Zhang, J.; Dai, W.; Gao, J.; Yu, H. Structural properties, antioxidant and immune activities of low molecular weight peptides from soybean dregs (Okara). Food Chem. X 2021, 12, 100175. [Google Scholar] [CrossRef] [PubMed]
  6. Sbroggio, M.F.; Montilha, M.S.; de Figueiredo, V.R.G.; Georgetti, S.R.; Kurozawa, L.E. Influence of the degree of hydrolysis and type of enzyme on antioxidant activity of okara protein hydrolysates. Food Sci. Technol. 2016, 36, 375–381. [Google Scholar] [CrossRef] [Green Version]
  7. Feng, J.Y.; Thakur, K.; Ni, Z.J.; Zhu, Y.Y.; Hu, F.; Zhang, J.G.; Wei, Z.J. Effects of okara and vitamin B2 bioenrichment on the functional properties and in vitro digestion of fermented soy milk. Food Res. Int. 2021, 145, 110419. [Google Scholar] [CrossRef]
  8. Suleman, M.; Riaz, S. 3D in silico study of magnetic fluid hyperthermia of breast tumor using Fe3O4 magnetic nanoparticles. J. Therm. Biol. 2020, 91, 102635. [Google Scholar] [CrossRef] [PubMed]
  9. Liu, D.; Chen, X.; Huang, J.; Huang, M.; Zhou, G. Generation of bioactive peptides from duck meat during post-mortem aging. Food Chem. 2017, 237, 408–415. [Google Scholar] [CrossRef] [PubMed]
  10. Hashemi, Z.S.; Zarei, M.; Fath, M.K.; Ganji, M.; Farahani, M.S.; Afsharnouri, F.; Pourzardosht, N.; Khalesi, B.; Jahangiri, A.; Rahbar, M.R.; et al. In silico Approaches for the Design and Optimization of Interfering Peptides Against Protein–Protein Interactions. Front. Mol. Biosci. 2021, 8, 282. [Google Scholar] [CrossRef]
  11. Gu, Y.; Majumder, K.; Wu, J. QSAR-aided in silico approach in evaluation of food proteins as precursors of ACE inhibitory peptides. Food Res. Int. 2011, 44, 2465–2474. [Google Scholar] [CrossRef]
  12. Iwaniak, A.; Dziuba, J. Animal and Plant Proteins as Precursors of Peptides with ACE Inhibitory Activity—An in silico Strategy of Protein Evaluation. Food Technol. Biotechnol. 2009, 47, 441–449. [Google Scholar]
  13. Minkiewicz, P.; Dziuba, J.; Iwaniak, A.; Dziuba, M.; Darewicz, M. BIOPEP database and other programs for processing bioactive peptide sequences. J. AOAC Int. 2008, 91, 965–980. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  14. Siti, H.; Munawaroh, H.; Gun, G.; Dea, J.; Aisyah, S.; Afina, V.; Ningrum, A.; Susanto, E.; Martha, L.; Kurniawan, I.; et al. In silico proteolysis and molecular interaction of tilapia (Oreochromis niloticus) skin collagen-derived peptides for environmental remediation. Environ. Res. 2022, 212, 113002. [Google Scholar]
  15. Ningrum, A.; Anggrahini, S.; Anggraeni, Z. In silico approach for identification of bioactive peptides in jack bean. Qual. Assur. Saf. Crops Foods 2018, 10, 20. [Google Scholar]
  16. López-Barrios, L.; Gutiérrez-Uribe, J.A.; Serna-Saldívar, S.O. Bioactive Peptides and Hydrolysates from Pulses and Their Potential Use as Functional Ingredients. J. Food Sci. 2014, 79, R273–R283. [Google Scholar] [CrossRef] [PubMed]
  17. Tonolo, F.; Moretto, L.; Grinzato, A.; Fiorese, F.; Folda, A.; Scalcon, V.; Ferro, S.; Arrigoni, G.; Bellamio, M.; Feller, E.; et al. Fermented soy-derived bioactive peptides selected by a molecular docking approach show antioxidant properties involving the keap1/nrf2 pathway. Antioxidants 2020, 9, 1306. [Google Scholar] [CrossRef]
  18. Yokomizo, A.; Takenaka, Y.; Takenaka, T. Antioxidative Activity of Peptides Prepared from Okara Protein. Food Sci. Technol. Res. 2002, 8, 357–359. [Google Scholar] [CrossRef] [Green Version]
  19. Maestri, E.; Marmiroli, M.; Marmiroli, N. Bioactive peptides in plant-derived foodstuffs. J. Proteom. 2016, 147, 140–155. [Google Scholar] [CrossRef]
  20. Wang, R.; Thakur, K.; Feng, J.Y.; Zhu, Y.Y.; Zhang, F.; Russo, P.; Spano, G.; Zhang, J.G.; Wei, Z.J. Functionalization of soy residue (okara) by enzymatic hydrolysis and LAB fermentation for B2 bio-enrichment and improved in vitro digestion. Food Chem. 2022, 387, 132947. [Google Scholar] [CrossRef]
  21. Montserrat-de la Paz, S.; Villanueva, A.; Pedroche, J.; Millan, F.; Martin, M.E.; Millan-Linares, M.C. Antioxidant and anti-inflammatory properties of bioavailable protein hydrolysates from lupin-derived agri-waste. Biomolecules 2021, 11, 1458. [Google Scholar] [CrossRef]
  22. Han, R.; Álvarez, A.J.H.; Maycock, J.; Murray, B.S.; Boesch, C. Comparison of alcalase- and pepsin-treated oilseed protein hydrolysates—Experimental validation of predicted antioxidant, antihypertensive and antidiabetic properties. Curr. Res. Food Sci. 2021, 4, 141–149. [Google Scholar] [CrossRef]
  23. Stanojevic, S.P.; Barac, M.B.; Pesic, M.B.; Jankovic, V.S.; Vucelic-Radovic, B.V. Bioactive proteins and energy value of okara as a byproduct in hydrothermal processing of soy milk. J. Agric. Food Chem. 2013, 61, 9210–9219. [Google Scholar] [CrossRef]
  24. Ahmad Nadzri, F.N.; Tawalbeh, D.; Sarbon, N.M. Physicochemical properties and antioxidant activity of enzymatic hydrolysed chickpea (Cicer arietinum L.) protein as influence by alcalase and papain enzyme. Biocatal. Agric. Biotechnol. 2021, 36, 102131. [Google Scholar] [CrossRef]
  25. Firmansyah, M.; Abduh, M.Y. Production of protein hydrolysate containing antioxidant activity from Hermetia illucens. Heliyon 2019, 5, e02005. [Google Scholar] [CrossRef] [Green Version]
  26. Atef, M.; Ojagh, S.M. Health benefits and food applications of bioactive compounds from fish byproducts: A review. J. Funct. Foods 2017, 35, 673–681. [Google Scholar] [CrossRef]
  27. Baehaki, A.; Nopianti, R.; Anggraeni, S. Antioxidant activity of skin and bone collagen hydrolyzed from striped catfish (Pangasius pangasius) with papain enzyme. J. Chem. Pharm. Res. 2015, 7, 131–135. [Google Scholar]
  28. Udenigwe, C.C. Bioinformatics approaches, prospects and challenges of food bioactive peptide research. Trends Food Sci. Technol. 2014, 36, 137–143. [Google Scholar] [CrossRef]
  29. 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]
  30. Fu, Y.; Young, J.F.; Løkke, M.M.; Lametsch, R.; Aluko, R.E.; Therkildsen, M. Revalorisation of bovine collagen as a potential precursor of angiotensin I-converting enzyme (ACE) inhibitory peptides based on in silico and in vitro protein digestions. J. Funct. Foods 2016, 24, 196–206. [Google Scholar] [CrossRef]
  31. Lo, W.M.Y.; Li-Chan, E.C.Y. Angiotensin I Converting Enzyme Inhibitory Peptides from In Vitro Pepsin−Pancreatin Digestion of Soy Protein. J. Agric. Food Chem. 2005, 53, 3369–3376. [Google Scholar] [CrossRef] [PubMed]
  32. Roy, A.; Nair, S.; Sen, N.; Soni, N.; Madhusudhan, M.S. In silico methods for design of biological therapeutics. Methods 2017, 131, 33–65. [Google Scholar] [CrossRef] [PubMed]
  33. Barati, M.; Javanmardi, F.; Jabbari, M.; Mokari-Yamchi, A.; Farahmand, F.; Eş, I.; Farhadnejad, H.; Davoodi, S.H.; Mousavi Khaneghah, A. An in silico model to predict and estimate digestion-resistant and bioactive peptide content of dairy products: A primarily study of a time-saving and affordable method for practical research purposes. LWT 2020, 130, 109616. [Google Scholar] [CrossRef]
  34. 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]
  35. Herowati, R.; Widodo, G.P. Molecular Docking Studies of Chemical Constituents of Tinospora cordifolia on Glycogen Phosphorylase. Procedia Chem. 2014, 13, 63–68. [Google Scholar] [CrossRef] [Green Version]
  36. Nongonierma, A.B.; FitzGerald, R.J. Enhancing bioactive peptide release and identification using targeted enzymatic hydrolysis of milk proteins. Anal. Bioanal. Chem. 2018, 410, 3407–3423. [Google Scholar] [CrossRef]
  37. Kallenberg, C.G.M. Antineutrophil Cytoplasmic Autoantibodies with Specificity for Myeloperoxidase, 3rd ed.; Elsevier: Amsterdam, The Netherlands, 2013; ISBN 9780444563781. [Google Scholar]
Figure 1. 2D Ligand Protein interaction. (a) PHF; (b)YYL; (c) YNL; (d) HL; (e) IR; (f) YYV; (g) EL.
Figure 1. 2D Ligand Protein interaction. (a) PHF; (b)YYL; (c) YNL; (d) HL; (e) IR; (f) YYV; (g) EL.
Molecules 27 06855 g001aMolecules 27 06855 g001b
Figure 2. Binding free energy (binding activity) of ligand (peptides).
Figure 2. Binding free energy (binding activity) of ligand (peptides).
Molecules 27 06855 g002
Figure 3. Example of 3D interaction of ligand-protein -YYV.
Figure 3. Example of 3D interaction of ligand-protein -YYV.
Molecules 27 06855 g003
Table 1. Profile of biological activity as an antioxidative peptide of glycinin and to the β and α subunits of β-conglycinin.
Table 1. Profile of biological activity as an antioxidative peptide of glycinin and to the β and α subunits of β-conglycinin.
NoPrecursorActivityAB
1Glycininantioxidative0.0543
2Conglycinin alpha subunitantioxidative0.05293.01898 × 10−6
3Conglycinin beta subunitantioxidative0.09344.16056 × 10−6
Table 2. In silico assay using papain and bromelain release peptides with antioxidant activity.
Table 2. In silico assay using papain and bromelain release peptides with antioxidant activity.
NoProtein PrecursorPeptide IDSequenceLocationMonoisotopic MassChemical Mass
1Glycinin3317HL(84–850268.142268.302
2 10019YNL(508–510)408.1899408.434
beta conglycinin alpha subunit
1 7888EL(190–191)260.126260.276
2 7888EL(513–514)260.126260.276
3 7942YYV(301–303)443.194443.473
4 8025PHF(465–467)399.178399.436
5 8215IR(220–221)287.185287.348
beta conglycinin beta subunit
1 7888EL(218–219)260.126260.276
2 7888EL(341–342)260.126260.276
3 7941YYL(134–136)457.21457.5
4 8025PHF(295–297)399.178399.436
Table 3. Physicochemical characterization of peptides derives from glycinin and conglycinin with antioxidant activities.
Table 3. Physicochemical characterization of peptides derives from glycinin and conglycinin with antioxidant activities.
NoSequenceSolubilityNet-ChargeIso-Electric Point (pH)Extinction Coefficient
1PHFPoor0.18.260
2YYLPoor03.322560 M−1cm−1
3YNLPoor03.341280 M−1cm−1
4HLPoor0.17.560
5IRGood110.850
6YYVPoor03.352560 M−1cm−1
7ELGood−10.920
Table 4. PeptideRanker and sensory evaluation and allergenic prediction of peptides from glycinin and conglycinin as antioxidant peptides.
Table 4. PeptideRanker and sensory evaluation and allergenic prediction of peptides from glycinin and conglycinin as antioxidant peptides.
NumberPeptidePeptide RankerSensory EvaluationAllergenicityToxicity
1PHF0.938016bitterProbable Non AllergenNon Toxic
2YYL0.598987bitterProbable Non AllergenNon Toxic
3YNL0.375314bitterProbable Non AllergenNon Toxic
4HL0.374865bitterProbable Non AllergenNon Toxic
5IR0.332363bitterProbable Non AllergenNon Toxic
6YYV0.188713bitterNon AllergenicNon Toxic
7EL0.0728272umamiProbable Non AllergenNon Toxic
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Ningrum, A.; Wardani, D.W.; Vanidia, N.; Sarifudin, A.; Kumalasari, R.; Ekafitri, R.; Kristanti, D.; Setiaboma, W.; Munawaroh, H.S.H. In Silico Approach of Glycinin and Conglycinin Chains of Soybean By-Product (Okara) Using Papain and Bromelain. Molecules 2022, 27, 6855. https://doi.org/10.3390/molecules27206855

AMA Style

Ningrum A, Wardani DW, Vanidia N, Sarifudin A, Kumalasari R, Ekafitri R, Kristanti D, Setiaboma W, Munawaroh HSH. In Silico Approach of Glycinin and Conglycinin Chains of Soybean By-Product (Okara) Using Papain and Bromelain. Molecules. 2022; 27(20):6855. https://doi.org/10.3390/molecules27206855

Chicago/Turabian Style

Ningrum, Andriati, Dian Wahyu Wardani, Nurul Vanidia, Achmat Sarifudin, Rima Kumalasari, Riyanti Ekafitri, Dita Kristanti, Woro Setiaboma, and Heli Siti Halimatul Munawaroh. 2022. "In Silico Approach of Glycinin and Conglycinin Chains of Soybean By-Product (Okara) Using Papain and Bromelain" Molecules 27, no. 20: 6855. https://doi.org/10.3390/molecules27206855

APA Style

Ningrum, A., Wardani, D. W., Vanidia, N., Sarifudin, A., Kumalasari, R., Ekafitri, R., Kristanti, D., Setiaboma, W., & Munawaroh, H. S. H. (2022). In Silico Approach of Glycinin and Conglycinin Chains of Soybean By-Product (Okara) Using Papain and Bromelain. Molecules, 27(20), 6855. https://doi.org/10.3390/molecules27206855

Article Metrics

Back to TopTop