Screening and Identification of Novel Soluble Epoxide Hydrolase Inhibitors from Corn Gluten Peptides

The objective of this study was to investigate the soluble epoxide hydrolase (sEH) inhibitory properties of corn gluten peptides. In total, 400 dipeptides and 8000 tripeptides were first virtually screened by molecular docking and 30 potential sEH inhibitory peptides were selected. Among them, WEY, WWY, WYW, YFW, and YFY showed the highest sEH inhibitory activities with IC50 values of 55.41 ± 1.55, 68.80 ± 7.72, 70.66 ± 9.90, 96.00 ± 7.5, and 94.06 ± 12.86 μM, respectively. These five peptides all behaved as mixed-type inhibitors and were predicted to form hydrogen bond interactions mainly with Asp333, a key residue located in the catalytic active site of sEH. Moreover, it was found that the corn gluten hydrolysates of Alcalase, Flavourzyme, pepsin and pancreatin all exhibited high sEH inhibitory activities, with IC50 values of 1.07 ± 0.08, 1.19 ± 0.24, and 1.46 ± 0.31 mg/mL, respectively. In addition, the sEH inhibitory peptides WYW, YFW, and YFY were successfully identified from the corn gluten hydrolysates by Alcalase using nano-LC-MS/MS. This study demonstrated the sEH inhibitory capacity of peptides for the first time and corn gluten might be a promising food protein source for discovering novel natural sEH inhibitory peptides.


Introduction
Soluble epoxide hydrolase (sEH) is a member of the α/β-hydrolase family and is widely expressed in tissues, such as the liver, kidney, heart, lung, brain, and intestine [1]. It is a homodimer and each monomer comprises a C-terminal hydrolase and an N-terminal phosphatase. The C-terminal hydrolase predominantly degrades epoxy-fatty acids, for example, epoxyeicosatrienoic acids (EETs), to their corresponding diols [2]. EETs are known to be generated from arachidonic acid via the cytochrome P450 epoxygenase pathway and have been shown to have anti-inflammatory, hypotensive, hypoglycemic, and neuroprotective properties. However, sEH can convert the active EETs to inactive dihydroxyeicosatrienoic acids (DHETs), causing them to lose their beneficial effects [3]. Therefore, sEH has become a new potential target for many metabolic, inflammatory, vascular, and neurodegenerative diseases [4]. Inhibiting the activity of sEH can effectively prevent the hydrolysis of EETs and increase the level of intracellular EETs so they can exert their physiological functions [5].
Discovering novel sEH inhibitors through chemical synthesis or from natural products has attracted a lot of attention. It has been found that most of the sEH inhibitors mainly include urea and amide groups [6]. Compared with synthesized drugs, natural inhibitors are considered safer and are thought to have fewer adverse effects. As a consequence, the interest in screening effective sEH inhibitors from natural sources has increased in the last decade [7]. Until now, many sEH inhibitors have been identified from traditional Chinese herbs, such as Alisma orientale [8], Scutellaria baicalensis [9], Leonurus japonicus [10], and Eucommia ulmoides [11]. However, these sEH inhibitors often have poor water solubility, poor intestinal absorption, and low bioavailability, resulting in undesirable pharmacokinetic properties and a decrease in potency after oral intake [12]. All this implies the importance of discovering possible sEH inhibitors from common nutritional components of foods.
Bioactive peptides prepared from food proteins by enzymatic hydrolysis or microbial fermentation have been considered promising nutraceutical treasures due to their excellent health-promoting benefits [13]. Corn gluten meal is a good food protein source for bioactive peptides and peptides generated from corn gluten have been demonstrated to have various biological activities, such as antioxidant, anti-inflammatory, antihypertensive, antihyperglycemic, hepatoprotective, alcohol metabolism-facilitating, anticancer, and antimicrobial activities [14]. In our previous study, it was found that peptides generated from corn gluten hydrolysis by the enzyme Alcalase had many free radicals and high intracellular reactive oxygen species scavenging capacity [15,16]. Moreover, corn gluten peptides could significantly suppress the expression of Epoxide Hydrolase 2 (EPHX2), the human gene encoding sEH, in hydrogen peroxide-induced HepG2 cells [17]. However, the effects of corn gluten peptides on sEH activity have never been reported.
Molecular docking has been widely used in the rapid screening of bioactive compounds by stimulating interactions between target proteins and candidate compounds. For example, Zhao et al. screened potential inhibitors against the main protease of SARS-CoV-2 from lactoferrin-derived peptides using molecular docking [18]. Yu et al. also successfully identified xanthine oxidase inhibitors from egg white-derived peptides using molecular docking [19]. Compared with traditional methods, the molecular docking-based in silico approach is highly efficient, time-saving, and low-cost. Therefore, in this study, the potential sEH inhibitory capacity of peptides was first explored through the computational docking of peptides into the active site of sEH. The peptides with a strong binding ability to sEH were selected and their sEH inhibitory activities were determined, as well as their inhibition modes and interactions with sEH. Furthermore, the inhibitory properties of corn gluten peptides from different proteases on sEH activity were studied and the corn gluten peptides were identified by nano-LC-MS/MS.

Preparation of Corn Gluten Peptides
Corn gluten peptides were prepared through enzymatic hydrolysis according to the method by You et al. with minor modifications [20]. The corn gluten meal was dissolved in deionized water to obtain the aqueous solution of corn gluten (5%) and heated at 90 • C for 10 min to denature the protein. After cooling down (Alcalase: 50 • C; Flavourzyme: 50 • C; papain: 55 • C; Neutrase proteases: 45 • C), the pH was adjusted with 1M NaOH or HCl (Alcalase: 8; Flavourzyme: 7; Papain: 7; Neutrase proteases: 7) followed by addition of different enzymes at an enzyme: substrate of 4%. The corn gluten was hydrolyzed for 180 min in a thermally controlled incubator under constant stirring. In addition, an in vitro simulated gastrointestinal digestion experiment of corn gluten was carried out. The corn gluten (2%) was first hydrolyzed by pepsin at pH 2 and 37 • C for 90 min and then by pancreatin at pH 7 and 37 • C for 240 min. Both pepsin and pancreatin were added at an enzyme: substrate of 2%. At the end of hydrolysis, all the resulting mixtures were heated at 90 • C for 10 min to inactivate the enzymes.
Then, all the corn gluten hydrolysates were centrifuged at 4000 r/min for 10 min. The supernatant was fractionated using ultrafiltration with a molecular weight cut-off membrane of 1 kDa (Shanghai Laungy Membrane Separation Equipment Engineering Co., Ltd., Shanghai, China). The corn gluten hydrolysate fractions with a molecular weight less than 1 kDa were selected, lyophilized, and stored at −20 • C.

Virtual Screening of Potential sEH Inhibitory Peptides
Molecular docking was performed to screen potential sEH inhibitory peptides. The X-ray structure of human sEH (PDB code: 1ZD3) was used as the receptor [21]. The structures of the native ligand (4-{[(cyclohexylamino)carbonyl]amino}butanoic acid, NC4) in the X-ray structure of human sEH and 400 dipeptides and 8000 tripeptides were defined as ligands in Molecular Operating Environment (MOE) (Chemical Computing Group Inc., Montreal, QC, Canada) after energy minimization. The binding site of sEH was set with the coordinates: x = −17.400, y = −7.960, and z = 66.430. The 2D structures of peptides were generated by the Protein Builder module in MOE. Three docking tools, MOE, AutoDock Vina, and CDOCKER were used to verify the better suitable docking method of redocking the native ligand NC4 to the binding site of sEH. The position of the native ligand NC4 in the X-ray structure of sEH was defined as the binding site.
For MOE, before docking, the force field of AMBER10: EHT and the implicit solvation model of the Reaction Field (R-field) were selected. The 2D structures of peptides were converted to 3D structures through an energy minimization module under the force field of AMBER10: EHT in MOE. The docking workflow followed the "induced fit" protocol, in which the side chains of the receptor pocket were allowed to move according to ligand conformations, with a constraint on their positions. The weight used for tethering side chain atoms to their original positions was 10. For each ligand, all docked poses of which were ranked by London dG scoring first, and then a force field refinement was carried out on the top 30 poses followed by a rescoring of GBVI/WSA dG. The conformation with the lowest binding free energy was finally identified as the best probable binding mode.
For AutoDock Vina, the X-ray structure of sEH was converted to a PDBQT file that contains a protein structure with hydrogens in all polar residues in AutoDock4. Additionally, the ligands were also converted to a PDBQT file and all the bonds were set to be rotatable in AutoDock4. All bonds of ligands were set to be rotatable. Lamarckian Genetic Algorithm (LGA) is used to calculate the flexible docking of protein-fixed ligands. AutoDock Vina Extended automatically calculates the grid maps and clusters the results in a way transparent to the user. The best conformation was chosen with the lowest binding energy after the docking search was completed.
For CDOCKER, the water molecules and the native ligands were removed from the complex crystal structure and then prepared by Prepare Protein module in Discovery Studio (DS). The 2D structures of the ligands were prepared (including adding hydrogen atoms, ionizing at the pH range from 6.5 to 8.5, generating tautomers and isomers, fixing bad valencies, and generating 3D coordinates) by Prepare Ligands module in DS. The peptides were prepared based on their sequence, and then modification and energy minimization were performed by the Energy Minimize procedure. Finally, the Prepare Ligands module was carried out to make the peptides ready before docking. The CDOCKER protocol was conducted between the ligands and the receptor. The conformation with the lowest CDOCKER Energy was finally recognized as the most stable conformation.
After the molecular docking, the interactions between the peptides and sEH including the hydrogen bonds and the bond lengths were analyzed using PyMol [22].

Validation of the Reaction Conditions between the sEH and Substrate
The effects of the concentrations of the enzyme, substrate, and the hydrolysis time of sEH on PHOME were determined according to the methods reported previously [11,23]. Briefly, 130 µL of 12.5 or 62.5 ng/mL human sEH solution (final concentration, dissolved in 25 mM bis-Tris-HCl containing 0.1% BSA, pH 7.0) and 20 µL of 10% methanol were added to a 96-well plate and incubated at room temperature for 15 min. Then, 50 µL of PHOME solutions at final concentrations of 0, 2.5, 5, 7.5, 10, 12.5, 15, 17.5, 20, and 22.5 µM (dissolved in 25 mM bis-Tris-HCl containing 0.1% BSA, pH 7.0) was added to initiate the reaction. After incubating at 37 • C for 5 to 50 min, the fluorescence intensity was measured at an emission wavelength of 465 nm and an excitation wavelength of 330 nm. Bis-Tris-HCl instead of PHOME was used as the blank control. The reaction rate of the system was calculated as follows: where S t (RFU) and C t (RFU) represent the fluorescence intensity of the sample and blank group at t min, respectively. The t (min) represents the reaction time.

Determination of the sEH Inhibitory Activities of Peptides
The sEH inhibitory activities of peptides were determined as above described. Briefly, 130 µL of 62.5 ng/mL human sEH solution (final concentrations, dissolved in 25 mM bis-Tris-HCl containing 0.1% BSA, pH 7.0) and 20 µL of peptides (final concentrations: synthesized peptides, 20, 40, 60, 80, and 100 µM; corn gluten peptides, 1, 2, 3, 4, and 5 mg/mL, respectively) were added to 96-well plate followed by adding 50 µL of 10 µM PHOME solution (final concentration, dissolved in 25 mM bis-Tris-HCl containing 0.1% BSA, pH 7.0). After incubating at 37 • C for 5 to 40 min, the fluorescence intensity of the mixture at an emission wavelength of 465 nm and an excitation wavelength of 330 nm was determined. AUDA and bis-Tris-HCl instead of peptide samples were used as the positive and blank control, respectively. The sEH inhibitory activity was calculated as follows: where C 40 (RFU) and C 0 (RFU) represented the fluorescence intensity of the blank control group at 40 and 0 min, respectively. S 40 (RFU) represented the fluorescence intensity of the sample group at 40 min. The median inhibitory concentration (IC 50 )-the concentration inhibiting 50% of sEH enzyme activity was calculated by through linear regression according to the percent inhibitory activity versus the sample concentration curve.
The sEH inhibition mode of peptides was determined using the Lineweaver−Burk double-reciprocal plot with PHOME and peptide concentrations ranging from 1 to 30 µM and 20 to 200 µM (final concentration), respectively.

Identification of Corn Gluten Peptides Using Nano LC-MS/MS
The peptide sequences of corn gluten hydrolysate were identified using nano-LC-MS/MS. Briefly, 5 µL of the sample was first loaded onto a nanocolumn (Acclaim PepMap RPLC C18, 150 µm × 150 mm, 3 µm, 100 Å; Dr. Maisch GmbH, Ammerbuch, Germany) in an Easy-nLC 1200 system (Thermo Fisher Scientific, Waltham, MA, USA). The sample was eluted by a linear gradient with mobile phase B (20% 0.1% formic acid in water-80% acetonitrile) from 4% to 8% for 2 min, from 8% to 28% for 43 min, from 28% to 40% for 10 min, from 40% to 95% for 1 min, and from 95% to 95% for 10 min. The mobile phase A was 0.1% formic acid in water. The flow rate was set at 600 nL/min. Then, the sample was analyzed using Q Exactive™ Hybrid Quadrupole-Orbitrap™ Mass Spectrometer (Thermo Fisher Scientific). The electrospray voltage was set at 2.2 kV and the capillary temperature was set at 270 • C. The m/z for the full scan was 300 to 1000. The activation type was selected as HCD and the normalized collision energy was set at 28. Up to the top 20, the most intense peptide ions from the preview scan in the Orbitrap were used for the data-dependent acquisition. The raw MS file was analyzed and searched against the UniProt Zea mays L. (corn) database using Byonic.

Physicochemical Property and Toxicity Analysis
In addition, the sEH inhibitory peptides were analyzed for physicochemical properties and toxicity. The physicochemical properties of the identified peptides, including isoelectric point (pI), net charge at neutral pH, and water solubility, were analyzed by the peptide property calculator, which is available at http://www.innovagen.com/proteomics (accessed on 10 October 2022). The prediction of peptide toxicity using ToxinPred is available at https://webs.iiitd.edu.in/raghava/toxinpred/multisubmit.php (accessed on 10 October 2022).

Statistical Analysis
All data were presented as mean ± SD (n = 3). The difference in significance between the two groups was analyzed by Student's t-test, and that between three and more groups was analyzed by one-way analysis of variance (ANOVA) followed by the Duncan multiple comparisons. Statistical analysis was performed by SPSS 19.0 with the significance level at p < 0.05.

The Optimal Reaction Conditions between the sEH and Substrate
PHOME is widely used in the high-throughput screening of sEH inhibitors. As a substrate, PHOME is highly sensitive, hydrolytically stable, and shows a large change in the fluorescence in the hydrolysis reaction [24]. When hydrolyzed by sEH, the nonfluorescent PHOME was converted to fluorescent 6-methoxy-2-naphthaldehyde. However, the inhibitor can competitively bind to the active site of sEH and lead to a decrease in the formation of the fluorescent products [25].
In the present study, the reaction conditions, including sEH and PHOME concentrations and reaction time, were optimized. A high reaction rate was observed when the sEH enzyme concentration of 62.5 ng/mL was used. The reaction rate increased significantly with the substrate concentration from 0 to 10 µM and decreased from 12.5 to 22.5 µM ( Figure 1A). In addition, the reaction rate increased from 5 to 40 min and was maintained from 40 to 50 min with substrate concentrations ranging from 7.5 to 17.5 µM ( Figure 1B). Therefore, sEH and PHOME concentrations at 62.5 ng/mL and 10 µM, respectively, and the reaction time of 40 min were selected for further determination of the sEH inhibitory activity.

Virtual Screening of sEH Inhibitory Peptides by Molecular Docking
In the past few years, the speed and accuracy of structural predictions and docking models have been significantly improved. Therefore, in this study, we carried out molecular docking of the sEH enzyme to screen potential sEH inhibitory peptides. It is known

Virtual Screening of sEH Inhibitory Peptides by Molecular Docking
In the past few years, the speed and accuracy of structural predictions and docking models have been significantly improved. Therefore, in this study, we carried out molecular docking of the sEH enzyme to screen potential sEH inhibitory peptides. It is known that the active side of sEH is located in a cavity with two narrow openings [6]. Therefore, in the present study, 400 dipeptides and 8000 tripeptides were selected and docked into the active site of sEH. Before the screening, three different docking tools-MOE, AutoDock Vina, and CDOCKER were validated by redocking the native ligand NC4 into the binding site of sEH. The root means square deviation (RMSD) between the redocked pose and the native pose in MOE, AutoDock Vina, and CDOCKER were calculated as 1.062, 1.939, and 2.445 Å, respectively. The lower RMSD means better redock performance. It was found that MOE showed the lowest RMSD, indicating that it had a good ability to reproduce the native conformation. The alignment between the docked pose and the native pose was as shown in Figure S1.
The 400 dipeptides and 8000 tripeptides were docked to sEH by MOE. After being docked, all peptides were ranked according to the score, namely, the binding free energy. The peptide-sEH complex with the lowest binding free energy was recognized as the most stable conformation and this peptide might have a strong binding ability to the active site of sEH. Finally, the top 20 tripeptides and the top 10 dipeptides were selected. It was noted that the tripeptides generally showed a lower binding free energy than dipeptides (Table 1). These 30 peptides were synthesized and their sEH inhibitory activities were determined. It was interesting to find that, at the final concentration of 100 µM, all the peptides showed sEH inhibitory activities ranging from 19.63% to 70.01% (Table 1). Among these peptides, WEY had the highest sEH inhibitory activity of 70.01% ± 6.25% followed by WYW, YFW, WWY, and YFY with 60.30% ± 5.60%, 57.38% ± 7.62%, 55.15%± 2.38%, and 52.96% ± 3.09%, respectively (p < 0.05). In addition, YMW, YRW, RV, FRY, WF, WLR, and YHY also exhibited sEH inhibitory activities of more than 40% at 100 µM. All this suggests that the small peptides might have potential as sEH inhibitors. Recently, a number of natural phytochemicals have been demonstrated to display sEH inhibitory activities and their IC 50 values mainly range from 1 to 100 µM. For instance, quercetin and kaempferol showed sEH inhibitory activities with IC 50 values of 22.5 ± 0.9 and 31.3 ± 2.6 µM, respectively [11]. In addition, quercetin-3-O-arabinoside, ursolic acid, corosolic acid, and 2-oxopomolic acid were found to exhibit sEH inhibitory activities with IC 50 values of 39.3 ± 3.4, 84.5 ± 9.5, 51.3 ± 4.9, and 11.4 ± 2.7 µM, respectively [26].
In our present study, it was found that the peptides WEY, WYW, YFW, WWY, and YFY all exhibited sEH inhibitory activities with a significant dose-effect from 20 to 100 µM (final concentration) (p < 0.05), as shown in Figure 2. Peptide WEY showed the lowest IC 50 value of 55.41 ± 1.55 µM, followed by WWY and WYW of 68.80 ± 7.72 and 70.66 ± 9.90 µM, and then YFW and YFY of 96.00 ± 7.5 and 94.06 ± 12.86 µM, respectively (p < 0.05) ( Table 1). While the IC 50 value of the positive control of AUDA was determined to be 63.53 ± 9.03 nM in this study. It is known that the peptide bond in peptides is one of the amide groups. In addition, the tripeptides and dipeptides have small linear structures which may be appropriate for the narrow cavity of the active site of sEH. All this may contribute to the sEH inhibitory activity of peptides. However, what should be noted is that the sEH inhibitory activity of peptides was determined for the first time in our present study. This implies the potential of peptides as natural sEH inhibitors.

Inhibition Mode of sEH Inhibitory Peptides
The inhibition kinetics of peptides WYW, WWY, YFW, YFY, and WEY on sEH determined. The reaction rates of the five peptides all increased and then remained high level with substrate PHOME ranging from 1 to 30 μM ( Figure S2). To determin inhibition mode of peptides, the Lineweaver−Burk double-reciprocal plots were ducted and the Michaelis constant (Km) and the maximum rate (Vmax) were calculate was found that, as shown in Figure 3, with the increase in peptide concentrations, the increased while the Vmax decreased, indicating that the five peptides all acted as com itive and non-competitive mixed inhibitors. What should be noted is that the adequa the Lineweaver−Burk double-reciprocal plots in this study may be due to the low bin affinity of peptides to sEH.

Inhibition Mode of sEH Inhibitory Peptides
The inhibition kinetics of peptides WYW, WWY, YFW, YFY, and WEY on sEH were determined. The reaction rates of the five peptides all increased and then remained at a high level with substrate PHOME ranging from 1 to 30 µM ( Figure S2). To determine the inhibition mode of peptides, the Lineweaver−Burk double-reciprocal plots were conducted and the Michaelis constant (K m ) and the maximum rate (V max ) were calculated. It was found that, as shown in Figure 3, with the increase in peptide concentrations, the Km increased while the Vmax decreased, indicating that the five peptides all acted as competitive and non-competitive mixed inhibitors. What should be noted is that the adequacy of the Lineweaver−Burk double-reciprocal plots in this study may be due to the low binding affinity of peptides to sEH.

Interactions between Peptides and sEH Enzyme
To further understand the inhibition mechanisms of peptides on sE lecular interactions between peptides and the sEH enzyme were anal known that the hydrolase catalytic pocket of sEH is primarily composed o Asp333-Asp495-His523 and two tyrosine residues, Tyr381 and Tyr465. W enters the active site, the oxygen atom of the epoxide group will be recog and Tyr465, which in turn activate the epoxide ring-opening by Asp333 [ key residue responsible for the hydrolase activity. The NH of the urea or many sEH inhibitors can act as a hydrogen bond donor and form a hydr the oxygen atom of Asp333. Therefore, various urea and amide-contain have been shown to be as competitive inhibitors for sEH [6].
In the present study, it was found that the peptides WEY, WYW, Y YFY were predicted to form multiple interactions mainly including hydr π stackings with the residues of the active site of sEH, as shown in Fig  mation on the interactions between peptides and the sEH enzyme is sum 2. The hydrogen atom of NH3 + of the N-terminal amino acid residue of W YFW were all predicted to behave as hydrogen bond donors to form hydr

Interactions between Peptides and sEH Enzyme
To further understand the inhibition mechanisms of peptides on sEH, the intermolecular interactions between peptides and the sEH enzyme were analyzed. It is well-known that the hydrolase catalytic pocket of sEH is primarily composed of a catalytic triad Asp333-Asp495-His523 and two tyrosine residues, Tyr381 and Tyr465. When an epoxide enters the active site, the oxygen atom of the epoxide group will be recognized by Tyr381 and Tyr465, which in turn activate the epoxide ring-opening by Asp333 [7]. Asp333 is the key residue responsible for the hydrolase activity. The NH of the urea or amide group of many sEH inhibitors can act as a hydrogen bond donor and form a hydrogen bond with the oxygen atom of Asp333. Therefore, various urea and amide-containing compounds have been shown to be as competitive inhibitors for sEH [6].
In the present study, it was found that the peptides WEY, WYW, YFW, WWY, and YFY were predicted to form multiple interactions mainly including hydrogen bonds and π stackings with the residues of the active site of sEH, as shown in Figure 4. The information on the interactions between peptides and the sEH enzyme is summarized in Table 2. The hydrogen atom of NH 3 + of the N-terminal amino acid residue of WEY, WYW, and YFW were all predicted to behave as hydrogen bond donors to form hydrogen bonds with the oxygen atom of Asp333. Two hydrogen bonds between the oxygen atom of Asp333 and the hydrogen atoms of the carbon backbone of the N-terminal amino acid residue and the NH of the middle amino acid residue of YFY were also observed. In addition, the peptides were predicted to form hydrogen bonds with other residues that belonged to the active site of sEH. It was noted that hydrogen bonds between WWY and the residues of Ser413 and Met418 not Asp333 of sEH were predicted. The hydrogen bonds between peptides and sEH all had satisfactory bond distances [30]. Moreover, various π stacking interactions, such as arene-cation stacking, arene-arene stacking, and arene-H stacking between peptides and the residues of sEH were formed. All of this might contribute to the binding and inhibition of peptides on sEH.   [21]. Molecular graphics were generated by PyMOL (http://www.pymol.org, accessed on 10 October 2022) [22].  [21]. Molecular graphics were generated by PyMOL (http://www.pymol.org, accessed on 10 October 2022) [22].
The amide group of chemically synthesized sEH inhibitors can form more hydrogen bond interactions with the key residues of the active site of sEH. For example, it was found that not only the NH of the amide group formed a hydrogen bond with the carboxyl of Asp333 (Asp335) but also the CO of the amide group formed hydrogen bonds with the hydroxy of Tyr381 (Tyr383) and Tyr465 (Tyr 466) of sEH [31,32]. In the present study, only a hydrogen bond between the NH of the amide group of YFY and the carboxyl of Asp333 of sEH was observed. For WEY, WYW, and YFW, it was the free NH 3 + not the NH of the amide group that formed hydrogen bond interactions with Asp333. In addition, different from the amide-based synthesized compounds, no hydrogen bond interactions between peptides Tyr381 and Tyr465 were observed. This might be associated with the relatively low sEH inhibitory activity of peptides than the amide-based synthetic inhibitors. Generally, the structure-activity relationships among these peptides can be used to further refine structural models of the sEH improving the power of future computational approaches.

Physicochemical Property and Toxicity Analysis of sEH Inhibitory Peptides
The physicochemical properties of the identified peptides were analyzed by the peptide property calculator (Table S1). The peptide WEY had a negative charge at pH 7 and good solubility. The other peptides showed poor water solubility and electroneutrality. All peptides were predicted without toxicity. These results provided the basis for further analysis of the pharmacokinetic properties including the absorption in the intestine, bioavailability, and half-time in animals and humans in the future.

The sEH Inhibitory Activities of Corn Gluten Peptides
Different enzymes were used to determine the sEH inhibitory activities of corn gluten hydrolysates. The results showed that corn gluten hydrolysates by Alcalase, papain, Flavourzyme, and pepsin and pancreatin all displayed good sEH inhibitory activities with a significant dose effect (p < 0.05) ( Figure 5). The corn gluten hydrolysates by Alcalase had the lowest IC 50 value of 1.07 ± 0.08 mg/mL. This is slightly lower than Flavourzyme, pepsin and pancreatin levels of 1.19 ± 0.24 and 1.46 ± 0.31 mg/mL, respectively, but no significance was observed (p > 0.05). However, the IC 50 values of the corn gluten hydrolysates produced by Alcalase, Flavourzyme, and pepsin and pancreatin were all significantly lower than that of papain of 2.94 ± 0.26 mg/mL (p < 0.05). Different from Alcalase, papain, Flavourzyme, and pepsin and pancreatin, the corn gluten hydrolysates produced by Neutrase proteases had a low sEH inhibitory activity of less than 30% even at the concentration of 5 mg/mL (data not shown). Corn gluten hydrolysates by Alcalase, Flavourzyme, and pepsin and pancreatin were thought to be promising sources of natural sEH inhibitors. What should be noted is that although the sEH inhibitory activity of corn glute drolysate was not as high as most of the traditional urea and amide-containing syn chemicals and some natural compounds [7], the corn gluten hydrolysates by Alc Flavourzyme, and pepsin and pancreatin could inhibit more than 90% of the sEH en activity at the concentration of 5 mg/mL. Moreover, food protein hydrolysates gen have good water solubility and few adverse effects. They can be given safely at far h concentrations compared with synthetic sEH inhibitors. This shows that food prote rived sEH inhibitory peptides have potential as biological ingredients in functional or can be considered as a food with additional benefits.
The peptide sequences from the corn gluten hydrolysates by Alcalase were iden by nano-LC-MS/MS to further understand the sEH inhibitory activity of corn glute drolysates. It was interesting to find that three of the above five peptides, WYW, and YFY, were successfully identified in the corn gluten hydrolysates by Alcalase (F 6). This may contribute to the superior sEH inhibitory activity of the corn gluten hy ysates by Alcalase. Nevertheless, more sEH inhibitory peptides should be identifie verified in ongoing studies. What should be noted is that although the sEH inhibitory activity of corn gluten hydrolysate was not as high as most of the traditional urea and amide-containing synthetic chemicals and some natural compounds [7], the corn gluten hydrolysates by Alcalase, Flavourzyme, and pepsin and pancreatin could inhibit more than 90% of the sEH enzyme activity at the concentration of 5 mg/mL. Moreover, food protein hydrolysates generally have good water solubility and few adverse effects. They can be given safely at far higher concentrations compared with synthetic sEH inhibitors. This shows that food proteinderived sEH inhibitory peptides have potential as biological ingredients in functional foods or can be considered as a food with additional benefits.
The peptide sequences from the corn gluten hydrolysates by Alcalase were identified by nano-LC-MS/MS to further understand the sEH inhibitory activity of corn gluten hydrolysates. It was interesting to find that three of the above five peptides, WYW, YFW, and YFY, were successfully identified in the corn gluten hydrolysates by Alcalase ( Figure 6). This may contribute to the superior sEH inhibitory activity of the corn gluten hydrolysates by Alcalase. Nevertheless, more sEH inhibitory peptides should be identified and verified in ongoing studies.

Conclusions
In this study, 8400 dipeptides and tripeptides were screened for potential sEH inhib itors using molecular docking. Finally, 30 peptides were selected and tested for sEH in hibitory activities. Peptides WEY, WWY, WYW, YFW, and YFY showed the highest activ ities, and all five synthetic peptides were competitive and non-competitive mixed inhibi tors. Furthermore, sEH inhibitory activities were found in corn gluten hydrolysates pro duced by Alcalase, Papain, Flavourzyme, pepsin and pancreatin. In addition, three sEH inhibitory peptides WYW, YFW, and YFY were successfully identified in the corn gluten hydrolysates by Alcalase. Our results suggest that corn gluten peptides have potential a natural sEH inhibitors and corn gluten may be a prominent protein source for sEH inhib itory peptides.
Supplementary Materials: The following supporting information can be downloaded at www.mdpi.com/xxx/s1, Figure S1: The alignment between MOE docked pose (cyan) and native pose (green); Figure

Conclusions
In this study, 8400 dipeptides and tripeptides were screened for potential sEH inhibitors using molecular docking. Finally, 30 peptides were selected and tested for sEH inhibitory activities. Peptides WEY, WWY, WYW, YFW, and YFY showed the highest activities, and all five synthetic peptides were competitive and non-competitive mixed inhibitors. Furthermore, sEH inhibitory activities were found in corn gluten hydrolysates produced by Alcalase, Papain, Flavourzyme, pepsin and pancreatin. In addition, three sEH inhibitory peptides WYW, YFW, and YFY were successfully identified in the corn gluten hydrolysates by Alcalase. Our results suggest that corn gluten peptides have potential as natural sEH inhibitors and corn gluten may be a prominent protein source for sEH inhibitory peptides.