Next Article in Journal
Effect of Enyzmes on the Quality and Predicting Glycaemic Response of Chinese Steamed Bread
Next Article in Special Issue
A Prospective Review on the Research Progress of Citric Acid Modified Starch
Previous Article in Journal
Trade for Food Security: The Stability of Global Agricultural Trade Networks
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Communication

Virtual Screening of Soybean Protein Isolate-Binding Phytochemicals and Interaction Characterization

1
College of Food Science and Nutritional Engineering, China Agricultural University, Beijing 100083, China
2
China National Engineering Research Center for Fruit & Vegetable Processing, Beijing 100083, China
3
Beijing Key Laboratory for Food Non-Thermal Processing, Beijing 100083, China
*
Author to whom correspondence should be addressed.
These authors contributed equally to this work.
Foods 2023, 12(2), 272; https://doi.org/10.3390/foods12020272
Submission received: 22 November 2022 / Revised: 28 December 2022 / Accepted: 1 January 2023 / Published: 6 January 2023

Abstract

:
Soybean protein isolate (SPI) and small molecule interactions have drawn more and more attention regarding their benefits for both parts, while research on large-scale investigations and comparisons of different compounds is absent. In this study, a high throughput virtual screening was applied on a phytochemical database with 1130 compounds to pinpoint the potential SPI binder. Pentagalloylglucose, narcissoside, poliumoside, isoginkgetin, and avicurin were selected as the top-five ranking molecules for further validation. Fluorescence quenching assays illustrated that isoginkgetin has a significantly higher apparent binding constant (Ka) of (0.060 ± 0.020) × 106 L·mol−1, followed by avicularin ((0.058 ± 0.010) × 106 L·mol−1), pentagalloylglucose ((0.049 ± 0.010) × 106 L·mol−1), narcissoside ((0.0013 ± 0.0004) × 106 L·mol−1), and poliumoside ((0.0012 ± 0.0006) × 106 L·mol−1). Interface characterization by MD simulation showed that protein residues E172, H173, G202, and V204 are highly involved in hydrogen bonding with the two carbonyl oxygens of isoginketin, which could be the crucial events in SPI binding. Van der Waals force was identified as the major driven force for isoginketin binding. Our study explored SPI–phytochemical interaction through multiple strategies, revealing the molecular binding details of isoginkgetin as a novel SPI binder, which has important implications for the utilization of the SPI–phytochemical complex in food applications.

1. Introduction

Soybean protein isolate (SPI) is a widely used plant-based protein containing all the essential amino acids needed to fulfill adults’ nutritional requirements. In addition, SPI exhibits many biological activities including decreasing serum lipids, improving intestinal microbiota, modulating immunity, etc. [1,2,3]. Although widely used in meat products and others [4], SPI has some unfavorable characteristics affecting its application, such as low solubility, poor emulsifiability, etc. [5]. Physical, chemical, and enzymatic modification methods had been utilized to improve the functional properties of SPI. However, physical and chemical methods have low specificity and enzymatic methods tend to produce unfavorable bitter peptides [6]. Recently, molecular interaction between proteins and small molecules has been developed as a new strategy to improve the functional properties of proteins, especially for plant-based proteins [7]. Two types of interactions are demonstrated between proteins and small molecules, namely, covalent and non-covalent interactions. Non-covalent bonding can be achieved through hydrogen bonds, van der Waals forces, and hydrophobic interactions [8]. For example, (-)-epigallocatechin-3-gallate (EGCG) noncovalently binds SPI and increases the foaming, emulsifying, and antioxidant properties of SPI [6]. Rutin significantly increases the emulsifying properties of myofibrillar protein [9] and improves the foaming properties of soybean protein [10]. Alternatively, covalent bonding can be achieved by enzymatic, thermal, or alkaline methods. Covalent binding formed between catechins and whey protein changed the protein secondary and tertiary structure and the conjugate improved foaming and emulsifying properties [11]. However, due to undesirable side products formation such as quinone polymers and intense operating conditions such as high temperature and high pH in covalent interaction [12], non-covalent modification is more preferred and environmentally friendly for utilization in the food industry. For example, the composite of β-lactoglobulin and EGCG was used to manufacture stable Pickering emulsions carrying bioactive compounds [13]. EGCG was also used to alter the surface hydrophobicity of thermally induced SPI to improve the gel strength, water-holding capacity, and rheological properties of SPI-based gels [14]. Overall, small molecule non-covalent binding has become a promising strategy for protein modification.
To explore more compounds with potential SPI-binding capacity, SPI–small molecule interactions need to be investigated on a larger scale and with multiple strategies. Phytochemicals have drawn more and more attention in modern food science due to their unique bioactivity and structural diversity [15,16,17,18,19]. Many phytochemicals have been reported for their protein-binding capacities, in which polyphenolic compounds are the most studied, such as catechins, resveratrol, curcumin, et al. Until now, research on the characterization of SPI–phytochemical interactions is still limited [16]. In this study, a high throughput virtual screening was applied to a phytochemical database consisting of 1130 compounds to pinpoint the potential SPI-binding phytochemicals. Compounds with the lowest docking free energies were validated by in vitro fluorescence quenching assays. MD simulation was applied to the best binder to interpret the interaction details. Our work compared the binding preference of SPI to different phytochemicals, which has important implications for the utilization of natural small molecules in the fabrication of SPI conjugates utilized in the food industry.

2. Materials and Methods

2.1. Materials

Soy protein isolate (90.0% protein) was purchased from Shandong Biological Products Co., Ltd. (Linyi, China). Pentagalloylglucose, Narcissoside, Poliumoside, Isoginkgetin, and Avicurin (purity > 90%) were purchased from Chengdu Push Bio-technology Co., Ltd. (Wuhou, China). All other chemicals and reagents were of analytical grade. Syringe filters were purchased from Sigma-Aldrich (St. Louis, MO, USA).

2.2. Virtual Screening of SPI-Binding Phytochemicals

The binding pocket on 11S glycinin (PDB ID: 1OD5) was determined in previous research [5]. Briefly, Schrödinger SiteMap (Schrödinger Release 2021-3, NewYork, NY, USA) was used to predict the binding pockets on soybean 11S glycinin. As a well-defined SPI binder, epigallocatechin gallate (EGCG) was used as a model compound and docked into the predicted pockets. The grid files were generated by the Receptor Grid Generation module and the docking was evaluated by the ligand docking function. The EGCG docking scores were used to select the optimal binding pocket. Then, a virtual screening was performed on a phytochemical database with 1130 compounds. We performed docking using Schrödinger Glide, where small molecules were placed separately into previously identified docking pockets [5], then calculated the docking score and root mean square deviation (RMSD). RMSD of each atom is used to resolve the molecular conformation change deviation. If RMSD values are less than 2, reproducing the protein–ligand interaction is considered to be successful.

2.3. Validation of SPI-Phytochemicals Binding by Fluorescence Quenching Assays

Binding of SPI and high-scored phytochemicals was characterized using a fluorescence quenching assay. Specifically, a series of concentration gradients of small molecules were mixed with protein solutions and added to the 96-well plate, with the final protein concentration maintained at 0.135 mg/mL. Excitation light at 280 nm was used and the emission spectra were scanned in the wavelength range of 310–450 nm.
The fluorescence in-filter effect is corrected using the following equation:
F c o r ( λ e x , λ e m ) = F o b s ( λ e x , λ e m )   C F p ( λ e x )   C F s ( λ e m ) F o b s ( λ e x , λ e m ) 10 ( A e m + A e x ) / 2
where CFP and CFS are the correction factors for excitation and emission light, respectively, Aex and Aem are the absorbances of emission and excitation light, Fobs is the observed fluorescence intensity, and Fcor is the corrected fluorescence intensity [20].
The Stern–Volmer equation was used to calculate the quenching constant:
F0/F = 1 + Ksv[Q].
F and F0 are fluorescence intensity with and without phytochemicals in the protein solution, respectively. Ksv is the quenching constant, and [Q] is the concentration of phytochemicals. Equation (2) was used to estimate the binding parameters:
lg[(F0− F)/F] = lgKa + nlg[Q].
Ka is the apparent binding constant, and n is the number of binding sites.

2.4. Molecular Dynamics Simulation

All MD simulations were performed using GROMACS 2021 [21] with Amberff99SBildn force field. The initial structures of 11S glycinin used in the MD simulation were derived from the PDB structure (ID:1OD5). The molecular structure parameters of isoginkgetin were obtained from ChemSpider and optimized by ORCA 5.0 [22] and Multiwfn [23], then Sobtop [24] was used to establish the GAFF force field. The complex was prepared by Pymol, using the conformation obtained by docking as the initial conformation [25]. First, SPI and isoginkgetin were placed in a dodecahedron box with a volume of 455.86 nm3 and the box was filled with TIP3P water molecules. To balance the charge, 13 Na+ were added to the system. Then, the system was energy optimized to convergence. Temperature was stabilized near 300 °C using NVT at 200 ps and pressure was stabilized near 1 bar using NPT at 1 ns. Finally, the limiting potential of the system was lifted and the simulation was performed for 50 ns.
Throughout simulations, NVT uses V-rescale coupling and NPT uses Parrinello–Rahman coupling with a time step of 2 fs. Van der Waals forces and short-range electrostatic interactions are truncated at a radius of 14 Å. The long-range electrostatic interactions are calculated using the particle mesh Ewald (PME) method. The root mean square deviation (RMSD) of the SPI skeleton, the radius of gyration (Rg) of the SPI, and the average distance between the SPI and ISO were calculated using the program that comes with gromacs 2021. Gmx_MMPBSA is a program based on Amber’s MMPBSA.py [26]. We can use it to calculate the complex’s binding free energy and the binding contribution of each amino acid residue using the molecular mechanics Poisson–Boltzmann surface area (MMPBSA) method.

2.5. Statistical Analysis

The statistical analysis software GraphPad (San Diego, CA, USA) was used to analyze the data. Comparisons of means were determined by Duncan’s test at the 5% significance level using a one-way analysis of variance (ANOVA).

3. Results and Discussion

3.1. Virtual Screening of Phytochemicals against Soybean Protein

A phytochemical library containing 1130 compounds was used for virtual screening against the binding pocket on 11S glycinin (PDB ID: 1OD5) determined in previous research [5]. The top 20 compounds are listed in Table 1. The results show that pentagalloylglucose exhibits the lowest binding energy (−9.91 kcal/mol), followed by narcissoside, poliumoside, isoginkgetin, avicularin, which are −8.222, −8.183, −7.984, −7.984 kcal/mol, respectively (Figure 1). It is worth noticing that the top 5 compounds are all glycoside-substituted polyphenols except for isoginkgetin. The stereo views of the selected 11S-phytochemical complex are shown in Figure 1. Detailed interactions are shown in Figure 2. The major amino acids involved in pentagalloylglucose binding include Arg161, His173, Met177, Glu200, and Gly202 (Figure 2A). For narcissoside, amino acids including Glu172, Thr176, Glu200, and Gly202 dominate the binding (Figure 2B). For poliumoside (Figure 2C), Glu172, Thr176, Glu200, and Gly202 are highly involved. For isoginkgetic, Glu172, His173, and Ser206 are involved in the interaction (Figure 2D). As for avicularin, it mainly interacts with Arg161, Glu172, His173, and Gly202 (Figure 2E). The results indicate that electrically charged amino acids (His and Arg with positive charge and Glu with negative charge) and polar amino acids (Thr and Ser) are highly involved in the phytochemical binding. It indicates that electrostatic interaction and hydrogen bonding could be the driving forces for SPI–phytochemical binding.

3.2. Characterization of SPI-Phytochemical Interaction by Fluorescence Quenching Assays

Fluorescence quenching assay is widely used in characterizing protein–ligand interactions. Fluorescence quenching is a process that lowers the fluorescence intensity of protein because of interaction with quencher molecules. A lot of small molecules have been proven to be fluorescence quenchers. In this study, the binding of top phytochemicals including pentagalloylglucose, narcissoside, poliumoside, isoginkgetin, and avicularin to SPI was investigated based on a fluorescence quenching assay. Fluorescence quenching assay was also carried out on SPI and n-butanol interaction, as a negative control (Figure S1).
1,2,3,4,6-O-pentagalloylglucose (PGG) exhibited the lowest binding free energy in virtual screening (−9.91 kcal/mol, see Table 1). It is found in a wide variety of herbals, including Punica granatum, Mangifera indica, and Elaeocarpus sylvestris. It has been used for cancer prevention because of its antioxidant and antiinflammation properties [27]. Figure 3B shows the fluorescence spectra of SPI at different concentrations of PGG. Fluorescence intensity significantly decreases with increasing PGG concentration (from 10 to 60 μM), indicating the interaction between PGG and SPI. It is worth noting that a red shift (from 337 nm to 340 nm) was observed with the addition of PGG (Figure 3B), implying that the fluorescent group of SPI was in a more hydrophobic environment in the presence of PGG. This could be due to a slight conformational change of SPI induced by PGG [5].
As shown in Figure 3C, F0/F and PGG concentration were fitted well by the Stern-Volmer equation (Equation (1)). The linear Stern–Volmer diagram illustrates that only one type of quenching mechanism occurs (dynamic or static). The quenching constant Ksv was calculated to be 0.047 ± 0.001 × 106 L·mol−1 (Table 2). The apparent binding constant Ka, and the average number of binding sites n were obtained by fitting experimental data to Equation (2). PGG–SPI interaction has a Ka of (0.049 ± 0.010) ×106 L·mol−1 and an average number of binding sites (n) as 1.00 ± 0.04.
Narcissoside is the phytochemical with the second lowest binding energy in virtual screening (−8.222 kcal/mol, see Table 1). It is a monomethoxy flavone derivative (Figure 4A), which has been reported to alleviate mitochondrial oxidative stress and have anti-acute myeloid leukemia effects [28]. Like PGG, adding narcissoside also leads to a decrease in the fluorescence intensity of SPI. With increasing concentration, red shift of fluorescence peak happens (from 340 to 345 nm), which indicates that the bonding model of narcissoside is similar to PGG. As shown in Figure 4 and Table 2, the quenching constant Ksv for narcissoside was fitted to be (0.020 ± 0.001) ×106 L·mol−1. The apparent binding constant for narcissoside Ka was obtained to be (0.0013 ± 0.0004) × 106 L·mol−1, which is 38 times lower than PGG, indicating a significantly lower binding capacity. The average number of binding sites (n) was fitted to be 1.71 ± 0.16.
Figure 5 shows the SPI fluorescence quenching by poliumoside, which has a docking energy of −8.183 kcal/mol in the virtual screening. Poliumoside is a phenylethanoid glycosides isolated from Brandisia hancei which is an advanced glycation product formation inhibitor and has anti-inflammatory and antioxidant activities [29,30,31,32]. Fluorescence intensity of SPI decreased with the increase in poliumoside concentration (from 10 to 60 µM; see Figure 5Β). Like PGG, A slight red shift (from 340 to 345 nm) was observed with increasing poliumoside concentration, indicating a conformational change. The binding parameters are shown in Table 2. Poliumoside showed a Ksv of (0.0094 ± 0.0010) × 106 L·mol−1. The apparent binding constant Ka was calculated to be (0.0012 ± 0.0006) × 106 L·mol−1, which is close to poliumoside and much lower than 1,2,3,4,6-O-pentagalloylglucose. The average number of binding sites n is 1.52 ± 0.18.
The highest apparent binding constant was observed for isoginkgetin-SPI interaction (Figure 6). Isoginkgetin is a bis-flavonoid compound isolated from Ginkgo biloba, attenuates lipopolysaccharide induced monoamine neurotransmitter deficiency and depression-like behaviors by downregulating the p38/NF-κB signaling pathway [33]. It shows a binding energy of −7.984 kcal/mol in molecular docking. A quenching constant Ksv of (0.042 ± 0.002) × 106 L·mol−1 was obtained and the apparent binding constant Ka was calculated to be (0.060 ± 0.020) ×106 L·mol−1, much higher than other compounds (Table 2). Besides, it did not perturb the red or blue shift of SPI fluorescence, representing a different binding mode. The average number of binding sites n is 0.98 ± 0.12.
Avicularin, a glycoside of quercetin from the leaves of guava, has been reported to possess a variety of biological properties such as anti-inflammatory, anti-allergic, antioxidant, hepatoprotective, anti-tumor activities, and anti-obesity [34,35]. Figure 7B shows the fluorescence emission spectra of SPI with the addition of avicularin. Unlike other compounds, a blue shift was observed, suggesting that the structure of SPI is more compacted in the presence of avicularin, indicating a different mode of conformational change caused by avicularin. As shown in Figure 7C and Table 2, the quenching constant Ksv for avicularin was fitted to be (0.040 ± 0.001) × 106 L·mol−1. The apparent binding constant Ka and the average number of binding sites n were calculated to be (0.058 ± 0.010) × 106 L·mol−1 and 0.90 ± 0.05, respectively.
Overall, fluorescence quenching assays have indicated distinct SPI-binding parameters and protein structural change among the top 5 phytochemicals. Isoginkgetin exhibited the biggest binding constant ((0.060 ± 0.020) × 106 L·mol−1), indicating that isoginkgetin has the highest binding affinity to SPI, which is followed by avicularin, pentagalloylglucose, narcissoside, and poliumoside. Thus, to investigate why SPI has such a preference for isoginkgetin, we further carried out an MD simulation to characterize the molecular details in SPI–isoginkgetin interaction.

3.3. Molecular Dynamics Simulation on SPI-isoginkgetin interaction

A 50 ns MD simulation was performed on isoginkgetin and soybean 11S glycinin. As shown in Figure 8A, the RMSD reflecting the magnitude of the protein motion equilibrated at 20 ns, indicating a more stable state after binding. The decrease in the radius of gyration (Rg) implies a more compacted complex structure. Overall, the simulation reaches equilibrium after 20 ns. Free energies of the amino acid residues with a distance less than 6 Å from the protein were plotted in Figure 8A, among which residues I171, M177, and V204 exhibited the lowest binding energy, indicating that they are highly involved in the binding. Besides, hydrogen bonding occupancy (Figure 8A) showed that residues E172, H173, G202, and V204 are highly involved in hydrogen bonding to isoginkgetin [36]. In the representative stereo structure of the 11S–isoginkgetin complex shown in Figure 8B, the backbone amide nitrogen of E172 formed a hydrogen bond with the carbonyl oxygen 1 in the C ring of isoginkgetin. The backbone amide nitrogen of G202 formed a hydrogen bond with the carbonyl oxygen 2 in C’ ring. This result indicated the crucial role of carbonyl groups in 11S binding. Among the top five compounds, isoginkgetin is the only one with two carbonyl groups except for pentagalloylglucose. The low binding of pentagalloylglucose to SPI could be due to the steric hindrance caused by the high molecular weight. To further investigate the major driven forces for isoginkgetin–11S interaction, we calculated 200 frames at 30–50 ns using gmx_MMPBSA. Binding free energies for the complexes are shown in Table 3. The total binding energy is −34.64 kcal/mol. Van der Waals interaction is the main driven force of isoginkgetin–11S binding, which is up to −64.15 kcal/mol. Electrostatic interactions and nonpolar solvation are also involved in the binding.

4. Conclusions

SPI interaction is highly dependent on the chemical structure, with a high tendency for polyphenol binding. We found that isoginkgetin has the highest apparent binding constant (Ka) of (0.060 ± 0.020) × 106 L·mol−1, followed by avicularin ((0.058 ± 0.010) × 106 L·mol−1), pentagalloylglucose ((0.049 ± 0.010) × 106 L·mol−1), narcissoside ((0.0013 ± 0.0004) × 106 L·mol−1) and poliumoside ((0.0012 ± 0.0006) × 106 L·mol−1). The interaction between soybean 11S glycinin and isoginkgetin is mostly driven by van der Waals force. Protein residues E172, H173, G202, and V204 are highly involved in hydrogen bonding with carbonyl oxygens of isoginketin, which could be the crucial events in SPI binding. Our study located isoginkgetin as a novel SPI binder, which could be a potential modifier used in SPI-based emulsions or gels, likely in the manufacture of plant-based drinks and meat analogs. Our study also provides a research strategy combining virtual screening, fluorescence quenching, and MD simulation, to characterize plant protein and phytochemicals interaction, which has important implications for the utilization of SPI–phytochemical complex in the food matrix.

Supplementary Materials

The following supporting information can be downloaded at: https://www.mdpi.com/article/10.3390/foods12020272/s1, Figure S1: Characterization of SPI-n-butanol interaction by fluorescence quenching assay, as a negative control. (A) Molecular chemical structure of n-butanol. (B) The fluorescence emission spectra of SPI titrated with different concentrations of n-butanol. (C) The Stern-Volmer plots for binding parameter calculation. The obtained F0/F and lg[(F0 − F)/F] were not fitted well with Equation (2) (Figure S1C, R2 = −0.9734) and Equation (3) (Figure S1C, R2 = 0.7234) due to the poor concentration dependence.

Author Contributions

Conceptualization, J.Z. and Y.S.; methodology, J.Z., P.L. and A.W.; software, P.L. and A.W.; formal analysis, P.L. and A.W.; investigation, P.L. and A.W.; writing—original draft preparation, P.L. and A.W.; writing—review and editing, P.L., A.W., Y.S. and J.Z.; funding acquisition, J.Z. and Y.S. All authors have read and agreed to the published version of the manuscript.

Funding

This research was supported by the Young Elite Scientists Sponsorship Program by China Association for Science and Technology (CAST) to J.Z.

Data Availability Statement

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

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Sirtori, C.R.; Eberini, I.; Arnoldi, A. Hypocholesterolaemic effects of soya proteins: Results of recent studies are predictable from the Anderson meta-analysis data. Br. J. Nutr. 2007, 97, 816–822. [Google Scholar] [CrossRef] [Green Version]
  2. An, C.; Kuda, T.; Yazaki, T.; Takahashi, H.; Kimura, B. Caecal fermentation, putrefaction and microbiotas in rats fed milk casein, soy protein or fish meal. Appl. Microbiol. Biotechnol. 2013, 98, 2779–2787. [Google Scholar] [CrossRef] [PubMed]
  3. Zhang, J.; Li, H.; Xiao, L. Immune Modulation of Soybean Protein and Soybean Oligopeptides on Epidermal Trauma Infection in Aged Mice with Negative Nitrogen Balance. Food Sci. 2018, 39, 145–151. [Google Scholar] [CrossRef]
  4. Saz, J.M.; Marina, M.L. High performance liquid chromatography and capillary electrophoresis in the analysis of soybean proteins and peptides in foodstuffs. J. Sep. Sci. 2007, 30, 431–451. [Google Scholar] [CrossRef] [PubMed]
  5. Ao, L.; Liu, P.; Wu, A.; Zhao, J.; Hu, X. Characterization of Soybean Protein Isolate-Food Polyphenol Interaction via Virtual Screening and Experimental Studies. Foods 2021, 10, 2813. [Google Scholar] [CrossRef]
  6. You, Y.; Yang, L.; Chen, H.; Xiong, L.; Yang, F. Effects of (−)-Epigallocatechin-3-gallate on the Functional and Structural Properties of Soybean Protein Isolate. J. Agric. Food Chem. 2021, 69, 2306–2315. [Google Scholar] [CrossRef]
  7. Vesic, J.; Stambolic, I.; Apostolovic, D.; Milcic, M.; Stanic-Vucinic, D.; Velickovic, T.C. Complexes of green tea polyphenol, epigalocatechin-3-gallate, and 2S albumins of peanut. Food Chem. 2015, 185, 309–317. [Google Scholar] [CrossRef]
  8. Meng, Y.; Li, C. Conformational changes and functional properties of whey protein isolate-polyphenol complexes formed by non-covalent interaction. Food Chem. 2021, 364, 129622. [Google Scholar] [CrossRef]
  9. Cheng, J.; Zhu, M.; Liu, X. Insight into the conformational and functional properties of myofibrillar protein modified by mulberry polyphenols. Food Chem. 2019, 308, 125592. [Google Scholar] [CrossRef]
  10. Ye, J.; Deng, L.; Wang, Y.; McClements, D.J.; Luo, S.; Liu, C. Impact of rutin on the foaming properties of soybean protein: Formation and characterization of flavonoid-protein complexes. Food Chem. 2021, 362, 130238. [Google Scholar] [CrossRef]
  11. Jia, Z.; Zheng, M.; Tao, F.; Chen, W.; Huang, G.; Jiang, J. Effect of covalent modification by (−)-epigallocatechin-3-gallate on physicochemical and functional properties of whey protein isolate. LWT 2016, 66, 305–310. [Google Scholar] [CrossRef]
  12. Yildirim-Elikoglu, S.; Erdem, Y.K. Interactions between milk proteins and polyphenols: Binding mechanisms, related changes, and the future trends in the dairy industry. Food Rev. Int. 2017, 34, 665–697. [Google Scholar] [CrossRef]
  13. Su, J.; Guo, Q.; Chen, Y.; Mao, L.; Gao, Y.; Yuan, F. Utilization of β-lactoglobulin- (−)-Epigallocatechin- 3-gallate(EGCG) composite colloidal nanoparticles as stabilizers for lutein pickering emulsion. Food Hydrocoll. 2019, 98, 105293. [Google Scholar] [CrossRef]
  14. Xu, J.; Guo, S.; Li, X.; Jiang, S.; Zhong, X.; Zheng, Z. Gel properties of transglutaminase-induced soy protein isolate–polyphenol complex: Influence of epigallocatechin-3-gallate. J. Sci. Food Agric. 2021, 101, 3870–3879. [Google Scholar] [CrossRef] [PubMed]
  15. Hassan, S.; Azim, S. A Therapeutic Connection between Dietary Phytochemicals and ATP Synthase. Curr. Med. Chem. 2017, 24, 3894–3906. [Google Scholar] [CrossRef]
  16. Shahidi, F.; Pan, Y. Influence of food matrix and food processing on the chemical interaction and bioaccessibility of dietary phytochemicals: A review. Crit. Rev. Food Sci. Nutr. 2021, 62, 6421–6445. [Google Scholar] [CrossRef] [PubMed]
  17. Zembyla, M.; Murray, B.S.; Radford, S.J.; Sarkar, A. Water-in-oil Pickering emulsions stabilized by an interfacial complex of water-insoluble polyphenol crystals and protein. J. Colloid Interface Sci. 2019, 548, 88–99. [Google Scholar] [CrossRef]
  18. Le Bourvellec, C.; Renard, C. Interactions between Polyphenols and Macromolecules: Quantification Methods and Mechanisms. Crit. Rev. Food Sci. Nutr. 2012, 52, 213–248. [Google Scholar] [CrossRef]
  19. Müller, L.; Caris-Veyrat, C.; Lowe, G.; Böhm, V. Lycopene and Its Antioxidant Role in the Prevention of Cardiovascular Diseases—A Critical Review. Crit. Rev. Food Sci. Nutr. 2014, 56, 1868–1879. [Google Scholar] [CrossRef]
  20. Panigrahi, S.K.; Mishra, A.K. Inner filter effect in fluorescence spectroscopy: As a problem and as a solution. J. Photochem. Photobiol. C Photochem. Rev. 2019, 41, 100318. [Google Scholar] [CrossRef]
  21. Van Der Spoel, D.; Lindahl, E.; Hess, B.; Groenhof, G.; Mark, A.E.; Berendsen, H.J.C. GROMACS: Fast, flexible, and free. J. Comput. Chem. 2005, 26, 1701–1718. [Google Scholar] [CrossRef] [PubMed]
  22. Neese, F. Software update: The ORCA program system—Version 5.0. WIREs Comput. Mol. Sci. 2022, 12, e1606. [Google Scholar] [CrossRef]
  23. Lu, T.; Chen, F. Multiwfn: A multifunctional wavefunction analyzer. J. Comput. Chem. 2012, 33, 580–592. [Google Scholar] [CrossRef] [PubMed]
  24. Lu, T. Sobtop, Version [1.0(dev3)]. Available online: https://sobereva.com/soft/Sobtop (accessed on 26 March 2022).
  25. Liang, L.; Duan, W.; Zhang, J.; Huang, Y.; Zhang, Y.; Sun, B. Characterization and molecular docking study of taste peptides from chicken soup by sensory analysis combined with nano-LC-Q-TOF-MS/MS. Food Chem. 2022, 383, 132455. [Google Scholar] [CrossRef] [PubMed]
  26. Valdés-Tresanco, M.S.; Valiente, P.A.; Moreno, E. gmx_MMPBSA: A New Tool to Perform End-State Free Energy Calculations with GROMACS. J. Chem. Theory Comput. 2021, 17, 6281–6291. [Google Scholar] [CrossRef] [PubMed]
  27. Yoon, S.; Kim, D.; Roh, K.M.; Ahn, D.; Kang, H.J.; Chung, S.J. Identification of Vaccinia -H1 Related Phosphatase as an Anticancer Target for 1,2,3,4,6- O -Pentagalloylglucose. Chem. Biodivers. 2019, 17, e1900414. [Google Scholar] [CrossRef]
  28. Liu, T.; Cao, L.; Zhang, T.; Fu, H. Molecular docking studies, anti-Alzheimer’s disease, antidiabetic, and anti-acute myeloid leukemia potentials of narcissoside. Arch. Physiol. Biochem. 2020, 19, 1–11. [Google Scholar] [CrossRef]
  29. Deng, R.; Xu, Y.; Feng, F.; Liu, W. Identification of poliumoside metabolites in rat feces by high performance liquid chromatography coupled with quadrupole time-of-flight tandem mass spectrometry. J. Chromatogr. B 2014, 969, 285–296. [Google Scholar] [CrossRef]
  30. Yu, S.Y.; Lee, I.-S.; Jung, S.-H.; Lee, Y.M.; Lee, Y.-R.; Kim, J.-H.; Sun, H.; Kim, J.S. Caffeoylated Phenylpropanoid Glycosides from Brandisia hancei Inhibit Advanced Glycation End Product Formation and Aldose Reductase in Vitro and Vessel Dilation in Larval Zebrafish in Vivo. Planta Medica 2013, 79, 1705–1709. [Google Scholar] [CrossRef] [Green Version]
  31. He, Z.-D.; Lau, K.-M.; Xu, H.-X.; Li, P.-C.; But, P.P.-H. Antioxidant activity of phenylethanoid glycosides from Brandisia hancei. J. Ethnopharmacol. 2000, 71, 483–486. [Google Scholar] [CrossRef]
  32. Wu, A.; Yang, Z.; Huang, Y.; Yuan, H.; Lin, C.; Wang, T.; Zhao, Z.; Zhou, Y.; Zhu, C. Natural phenylethanoid glycosides isolated from Callicarpa kwangtungensis suppressed lipopolysaccharide-mediated inflammatory response via activating Keap1/Nrf2/HO-1 pathway in RAW 264.7 macrophages cell. J. Ethnopharmacol. 2020, 258, 112857. [Google Scholar] [CrossRef] [PubMed]
  33. Li, P.; Zhang, F.; Li, Y.; Zhang, C.; Yang, Z.; Zhang, Y.; Song, C. Isoginkgetin treatment attenuated lipopolysaccharide-induced monoamine neurotransmitter deficiency and depression-like behaviors through downregulating p38/NF-κB signaling pathway and suppressing microglia-induced apoptosis. J. Psychopharmacol. 2021, 35, 1285–1299. [Google Scholar] [CrossRef] [PubMed]
  34. Vo, V.A.; Lee, J.-W.; Chang, J.-E.; Kim, J.-Y.; Kim, N.-H.; Lee, H.J.; Kim, S.-S.; Chun, W.; Kwon, Y.-S. Avicularin Inhibits Lipopolysaccharide-Induced Inflammatory Response by Suppressing ERK Phosphorylation in RAW 264.7 Macrophages. Biomol. Ther. 2012, 20, 532–537. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  35. Amadi, J.A.; Amadi, P.; Njoku, U.C.; Onitchi, C.L. Lettuce-avicularin treatment reverses insulin resistance through stimulation of glycolytic kinases and insulin signaling molecules. Iran. J. Basic Med. Sci. 2021, 24, 232–239. [Google Scholar] [CrossRef] [PubMed]
  36. Liang, L.; Zhou, C.; Zhang, J.; Huang, Y.; Zhao, J.; Sun, B.; Zhang, Y. Characteristics of umami peptides identified from por-cine bone soup and molecular docking to the taste receptor T1R1/T1R3. Food. Chem. 2022, 387, 132870. [Google Scholar] [CrossRef]
Figure 1. The stereo views of leading phytochemical–11S complexes from virtual screening. (A) Pentagalloylglucose, (B) narcissoside, (C) poliumoside, (D) isoginkgetin, (E) avicularin.
Figure 1. The stereo views of leading phytochemical–11S complexes from virtual screening. (A) Pentagalloylglucose, (B) narcissoside, (C) poliumoside, (D) isoginkgetin, (E) avicularin.
Foods 12 00272 g001
Figure 2. Molecular interaction diagram of top-ranking phytochemical binding with soybean 11S globulin. (A) Pentagalloylglucose, (B) narcissoside, (C) poliumoside, (D) isoginkgetin, (E) avicularin.
Figure 2. Molecular interaction diagram of top-ranking phytochemical binding with soybean 11S globulin. (A) Pentagalloylglucose, (B) narcissoside, (C) poliumoside, (D) isoginkgetin, (E) avicularin.
Foods 12 00272 g002
Figure 3. Characterization of SPI–pentagalloylglucose interaction by fluorescence quenching assay. (A) Molecular chemical structure of pentagalloylglucose. (B) The fluorescence emission spectra of SPI titrated with different concentrations of pentagalloylglucose (10, 20, 40, 50, 60 μΜ). (C) The Stern–Volmer plots for binding parameter calculation.
Figure 3. Characterization of SPI–pentagalloylglucose interaction by fluorescence quenching assay. (A) Molecular chemical structure of pentagalloylglucose. (B) The fluorescence emission spectra of SPI titrated with different concentrations of pentagalloylglucose (10, 20, 40, 50, 60 μΜ). (C) The Stern–Volmer plots for binding parameter calculation.
Foods 12 00272 g003
Figure 4. Characterization of SPI–narcissoside interaction by fluorescence quenching assay. (A) Molecular chemical structure of narcissoside. (B) The fluorescence emission spectra of SPI titrated with different concentrations of narcissoside (10, 20, 40, 50, 60 μΜ). (C) The Stern–Volmer plots for binding parameter calculation.
Figure 4. Characterization of SPI–narcissoside interaction by fluorescence quenching assay. (A) Molecular chemical structure of narcissoside. (B) The fluorescence emission spectra of SPI titrated with different concentrations of narcissoside (10, 20, 40, 50, 60 μΜ). (C) The Stern–Volmer plots for binding parameter calculation.
Foods 12 00272 g004
Figure 5. Characterization of SPI–poliumoside interaction by fluorescence quenching assay. (A) Molecular chemical structure of poliumoside. (B) The fluorescence emission spectra of SPI titrated with different concentrations of poliumoside (10, 20, 40, 50, 60, μΜ). (C) The Stern–Volmer plots for binding parameter calculation.
Figure 5. Characterization of SPI–poliumoside interaction by fluorescence quenching assay. (A) Molecular chemical structure of poliumoside. (B) The fluorescence emission spectra of SPI titrated with different concentrations of poliumoside (10, 20, 40, 50, 60, μΜ). (C) The Stern–Volmer plots for binding parameter calculation.
Foods 12 00272 g005
Figure 6. Characterization of SPI–isoginkgetin interaction by fluorescence quenching assay. (A) Molecular chemical structure of isoginkgetin. (B) The fluorescence emission spectra of SPI titrated with different concentrations of isoginkgetin (10, 20, 25, 30, 35 μΜ). (C) The Stern–Volmer plots for binding parameter calculation.
Figure 6. Characterization of SPI–isoginkgetin interaction by fluorescence quenching assay. (A) Molecular chemical structure of isoginkgetin. (B) The fluorescence emission spectra of SPI titrated with different concentrations of isoginkgetin (10, 20, 25, 30, 35 μΜ). (C) The Stern–Volmer plots for binding parameter calculation.
Foods 12 00272 g006
Figure 7. Fluorescence spectroscopy of SPI–avicularin interaction. (A) Molecular structure of avicularin and fluorescence spectra of SPI at different concentrations of avicularin (10, 20, 40, 50, 60, 80 μΜ). (B) The Stern–Volmer curve chart. (C) The Stern-Volmer plots for binding parameter calculation.
Figure 7. Fluorescence spectroscopy of SPI–avicularin interaction. (A) Molecular structure of avicularin and fluorescence spectra of SPI at different concentrations of avicularin (10, 20, 40, 50, 60, 80 μΜ). (B) The Stern–Volmer curve chart. (C) The Stern-Volmer plots for binding parameter calculation.
Foods 12 00272 g007
Figure 8. MD simulation of SPI and isoginkgetin interaction. (A) RMSD, radius of gyration, hydrogen bonding occupancy and amino acid energy contribution during simulation. (B) The representative stereo structure of 11S–isoginkgetin complex.
Figure 8. MD simulation of SPI and isoginkgetin interaction. (A) RMSD, radius of gyration, hydrogen bonding occupancy and amino acid energy contribution during simulation. (B) The representative stereo structure of 11S–isoginkgetin complex.
Foods 12 00272 g008
Table 1. Top 20 phytochemical compounds from virtual screening.
Table 1. Top 20 phytochemical compounds from virtual screening.
CompoundsCASMolecular Weight (Da)Docking Energy (kcal/mol)Ranking
Pentagalloylglucose14937-32-7940.7−9.911
Narcissoside604-80-8624.5−8.2222
Poliumoside94079-81-9770.7−8.1833
Isoginkgetin548-19-6566.5−7.9844
Avicularin572-30-5434.3−7.9845
Isoliensinine6817-41-0610.7−7.9386
Isoquercitrin482-35-9464.4−7.9137
Thonningianin A271579-11-4874.7−7.8238
2′-O-galloylhyperin53209-27-1616.5−7.7499
Forsythoside B81525-13-5756.7−7.74310
Crotonoside1818-71-9283.24−7.72711
Gnetol86361-55-9244.24−7.712
Oxyresveratrol29700-22-9244.24−7.68213
(-)-swainsonine136997-64-3712.7−7.63914
Amentoflavone1617-53-4538.5−7.58415
Brazilin474-07-7286.28−7.53816
Liquiritin Apioside74639-14-8550.5−7.49117
Rhodiosin86831-54-1610.5−7.48918
Chelidonine476-32-4353.4−7.45219
Rhapontigenin500-65-2258.269−7.44220
Table 2. Binding parameters of selected SPI–phytochemical compound interaction.
Table 2. Binding parameters of selected SPI–phytochemical compound interaction.
CompoundsKsv (×106 L·mol−1)Ka (×106 L·mol−1)n
Pentagalloylglucose0.047 ± 0.001 a0.049 ± 0.010 a1.00 ± 0.04 b
Narcissoside0.020 ± 0.001 c0.0013 ± 0.0004 b1.71 ± 0.16 a
Poliumoside0.0094 ± 0.0010 d0.0012 ± 0.0006 b1.52 ± 0.18 a
Isoginkgetin0.042 ± 0.002 b0.060 ± 0.020 a0.98 ± 0.12 b
Avicularin0.040 ± 0.001 b0.058 ± 0.010 a0.90 ± 0.05 b
a, b, c, d Significant difference in values with different letters in the same column.
Table 3. Binding energy contribution of isoginkgetin–11S complex.
Table 3. Binding energy contribution of isoginkgetin–11S complex.
Title 1Title 2
Van der Waals−64.15
Electrostatic−8.13
Polar solvation42.46
Nonpolar solvation−4.82
Ggas−72.28
Gsolv37.64
Total−34.64
Notes: Ggas = van der Waals + electrostatic, Gsolv = polar solvation + nonpolar solvation, Total = Ggas + Gsolv.
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

Liu, P.; Wu, A.; Song, Y.; Zhao, J. Virtual Screening of Soybean Protein Isolate-Binding Phytochemicals and Interaction Characterization. Foods 2023, 12, 272. https://doi.org/10.3390/foods12020272

AMA Style

Liu P, Wu A, Song Y, Zhao J. Virtual Screening of Soybean Protein Isolate-Binding Phytochemicals and Interaction Characterization. Foods. 2023; 12(2):272. https://doi.org/10.3390/foods12020272

Chicago/Turabian Style

Liu, Panhang, Annan Wu, Yi Song, and Jing Zhao. 2023. "Virtual Screening of Soybean Protein Isolate-Binding Phytochemicals and Interaction Characterization" Foods 12, no. 2: 272. https://doi.org/10.3390/foods12020272

APA Style

Liu, P., Wu, A., Song, Y., & Zhao, J. (2023). Virtual Screening of Soybean Protein Isolate-Binding Phytochemicals and Interaction Characterization. Foods, 12(2), 272. https://doi.org/10.3390/foods12020272

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