Unmasking the efficacy of Skipjack Tuna (Katsuwonus Pelamis) dark muscle hydrolyzate in lipid regulation: a promising component for functional food

ABSTRACT To effectively utilize the by-products of skipjack tuna (Katsuwonus Pelamis) for the preparation of peptides with lipid-lowering activity, alkaline proteinase was selected from six commercial proteases to hydrolyze dark muscle. An L16 (54) orthogonal design was employed to optimize the hydrolysis conditions, resulting in the production of a high-peptide-content hydrolyzate named KPHs-AL. This hydrolyzate exhibited the ability to significantly reduce TG content. KPHs-AL also showed higher persuasive pancreatic lipase (PL) and cholesterol esterase (CEase) inhibitory potential than various molecular weight (Mw) fractions and other enzymatic hydrolyzates, with IC50 values of 5.14 and 8.51 mg/mL, respectively. Furthermore, KPHs-AL exerted a combined effect of regulating lipid accumulation in FAA-induced steatosis HepG2 cells, leading to decreased total cholesterol (TC) and total triglyceride (TG). Based on the high binding likelihood and interaction with specific binding sites predicted by Pepsite2, 38 peptide sequences were identified as potential PL and CEase inhibitors. Hence, KPHs-AL holds promise as a beneficial component of functional food and exhibits the potential for anti-obesity effects through the alleviation of lipid accumulation


Introduction
Obesity is defined as a chronic, complex disease characterized by abnormal lipid metabolism and excessive fat accumulation, which may impair health.It affects various system in the body, including cardiovascular, aerodigestive, circulatory and genitourinary systems, and can lead to a range of non-communicable diseases.Recent statistics from the WHO Global Health Observatory show that globally, more than 1 billion children, adolescents and adults are obese.The prevalence rate has nearly tripled since 1975, resulting in a significant rise in healthcare costs associated with related disease.Consequently, obesity remains an ongoing global public health problem. [1,2]C-1: References should be number format?Pancreatic lipase (PL) catalyzes the lipid hydrolysis of triacylglycerols into glycerol monoesters and fatty acids, playing an essential role in the digestion and absorption of dietary fats in the digestive tract.Similarly, cholesterol esterase (CEase) primarily hydrolyzes dietary cholesterol esters, liberating cholesterol and free fatty acids tofacilitate the absorption process of the body. [3]Inhibiting the activities of both enzymes can effectively control the hydrolysis of dietary fats and cholesterol esters, subsequently regulating the plasma levels of triglycerides and cholesterol.This approach provides an effective strategy for treating metabolic disorders and managing obesity-associated complications. [4,5]In this regard, various synthetic drugs based on the principle of PL and CEase inhibition have been developed and proven useful in managing obesity-associated complications. [6]For instance, Orlistat has been clinically proven to be effective in significantly and sustainably treating obesity by inhibiting PL's enzymatic action.However, these drugs come with severe gastrointestinal adverse effects, reduced levels of fat-soluble vitamins, and high long-term costs, imposing certain limitations. [7]Given these concerns with synthetic drugs, there is a growing need to discover and understand safer, more economical natural alternatives.
In recent years, the exploration of bioactive peptides has opened up exciting new avenues for addressing lipid disorders. [8]Peptides, short chains of amino acids, have demonstrated a wide range of biological activities, including the ability to modulate lipid metabolism.These bioactive peptides have garnered significant attention due to their potential in improving lipid profiles and, consequently, their implications for health and medicine.Sajid Maqsood's group systematically investigated the potential anti-obesity activities of camel milk-derived bioactive peptides, including in vitro inhibition of PL [9] and CEase [10] of camel milk protein hydrolyzate, the antiobesity potential and interactions mechanism [11] of camel whey protein hydrolyzate, [12] and the inhibitory potential of camel casein hydrolyzate against key lipid-digesting enzymes. [13]Several proteins or by-products from the food industry are considered high-quality sources of peptides or hydrolyzates that interfere with fat digestion and absorption, and alleviate lipid accumulation, such as Cuminum cyminum seeds, [14] olive seeds, [15] Phoenix dactylifera L. seeds, [16] amaranth seeds [17] and brewer's spent grain. [18]kipjack Tuna (Katsuwonus Pelamis) is a species of significant commercial importance, primarily sought after for its use in canned food products [.19]As consumer demand for this species steadily increases, a substantial portion of the fish, including not just the desirable white meat but also the head, tail, viscera, bones, and dark muscle, inevitably becomes available as byproducts. [20]To date, a wide range of bioactive peptides has been identified from those processing by-products of tun, specifically, antioxidant and angiotensin-I-converting enzyme inhibitory peptides from dark muscle [21,22] or blood [23] or roe [24] protein hydrolyzates, antioxidant peptides from the in vitro gastrointestinal digestion products of the head [25] or milt protein hydrolyzate [26] or bone collagen, [27,28] high fischer ratio oligo-peptides hydrolyzate from minced meat. [28]eptides with lipid-lowering properties can be incorporated into functional foods and nutraceuticals, providing individuals with dietary options to naturally manage their lipid levels.These products offer a convenient and appealing way to promote heart health through everyday consumption.Nonetheless, there are still gaps in the research on the potential of skipjack tuna dark muscle protein hydrolyzate to alleviate lipid accumulation, especially the PL and CEase inhibitory activities.
The objective of this work was to explore the potential of skipjack tuna dark muscle protein hydrolyzates, generated using commercial enzymes, to inhibit two key metabolic enzymes and reduce lipid accumulation in HepG 2 cells.L16 (5 4 ) orthogonal design was used to optimize the hydrolysis process for the release of potential inhibitory peptides derived from dark muscle.The obtained hydrolyzate underwent characterization, including analysis of its molecular weight (Mw) distribution, amino acid (AA) composition, and identification of peptide sequences.In addition, the effects of these potent hydrolyzates were assessed through in vitro assays and in-silico methods to gain insights into the possible interactions of inhibitory peptides with pancreatic lipase (PL) and cholesterol esterase (CEase).

Preparation & optimization of katsuwonus pelamis dark muscle hydrolysates (KPHs)
The pre-homogenized Katsuwonus pelamis dark muscle was hydrolyzed as the working parameters listed in Table 1.A constant-flow pump system was equipped with 0.5 M sodium hydroxide and 0.5 M hydrochloric acid was equilibrated to the desired pH (Lohand, DG150, Hangzhou, China).After that, the reaction system was heated to 100°C to inactivate these enzymes.Dark muscle from Katsuwonus pelamis was enzymatically hydrolyzed using different proteases to obtain distinct hydrolyzates.Specifically, hydrolysis with neutral protease resulted in KPHs-NP, alkaline protease yielded KPHs-AL, papain produced KPHs-PP, acid protease generated KPHs-ACP, bromelain led to KPHs-BM, and pancreatin produced KPHs-PC.These hydrolyzates were designated based on the respective proteases used for the hydrolysis process.The molecular weight fractions of KPHs-AL were categorized as follows: KPHs-AL-I for those with a molecular weight greater than 10 kDa, KPHs-AL-II for those ranging between 5 kDa and 10 kDa, KPHs-AL-III for those within the 5 kDa to 3 kDa range, and KPHs-AL-Ⅳ for those spanning between 3 kDa and 0.15 kDa.The supernatant was collected by centrifugation (Himac, CF16RN, Tokyo, Japan) and lyophilized for further use (Marin Christ, α1-4 LD plus, Germany).A five-factor-four level orthogonal design [28] was applied to optimize the hydrolysis parameters to discriminate the lipid accumulation inhibitory effect and peptide content of different hydrolyzates.The peptides fractions were separated using a Cogent Scale TFF ultrafiltration system (Merck Millipore, Darmstadt, Germany) equipped with Mw cutoff membranes of 10, 5, 3 and 0.15 kDa, which are used for further bioactivity assay after dialyzation and lyophilization.

Peptide content
The peptide content of the hydrolyzates was determined by the trichloroacetic acid precipitationbiuret colorimetric method. [29]Specifically, the hydrolyzates were pre-precipitated with 10% trichloroacetic acid, the centrifuged supernatant was mixed with biuret reagent, and the absorbance was measured at 540 nm to calculate the peptide concentrations (mg/mL) with the standard curve of Gly-Gly-Tyr-Arg tetrapeptide.In addition, the peptide content (mg/mg) was normalized according to the mass of the lyophilized powder.

Molecular weight distribution
The molecular weight (Mw) distribution of the hydrolyzates was determined by size exclusion chromatography as the same process described in our previous methods. [30]Specifically, samples and the standard mix were eluted with acetonitrile/H 2 O/TFA (45:55:0.1)at a flow rate of 1.0 mL/min by using Agilent 1260 Infinity II HPLC System (Agilent Technologies, Inc., CA, USA) equipped with a TSK gel G2000 SWXL analytical column (7.8 × 300 mm, 5 µm, TOSOH, Tokyo, Japan).The function of standard Mw and retention time was used to calculate the Mw distribution of KPHs.

Amino acid composition
The amino acid (AA) composition of the hydrolyzates was determined by three processes [31] : ion exchange chromatography separation, ninhydrin derivatization and absorbance monitoring at 570 nm/440 nm (HITACHI L-8900, Tokyo, Japan).Prior to the 60 min standard protein hydrolyzate analysis, the samples and high-purity hydrochloric acid (containing 2% phenol, m/v) were maintained at 110°C, 100% nitrogen atmosphere for 22 h to fully hydrolyze.

Identification of peptide sequences
Peptide sequences identification for samples was conducted using a Q Exactive HF LC/MS system (Thermo Scientific, CA, USA).The sample was injected into an 80-min gradient system equipped with a Nano Chrom XP C18 column (250 mm × 7.5 μm, 3 μm, Eksigent, CA, USA).The binary mobile phase consisting of 0.1% formic acid in acetonitrile (solvent A) and water (solvent B) was used at a constant flow rate of 600 nL/min.The elution procedure was performed as follows: 0 min, 94% solvent A, 5 min, 88% solvent A, 57 min, 71% solvent A, 71 min, 60% solvent A, 72 min, 5% solvent A, 80 min, 5% solvent A. Data analysis and library search for the MS/MS information collected in positive ion mode were performed with PEAKS studio version 6.0.Only peptides with a Denovo resolution confidence >80% or contained in the protein sequences from the species were identified in the following step. [8]

Pancreatic lipase inhibition assay
Pancreatic lipase (PL) inhibitory assay was carried out according to the method previously described by with minor modifications. [32]The preincubation system (37°C, 10 min) consisted of 50 μL of KPHs or fractions, 50 μL of PL (1 mg/mL) in 0.1 M Tris-HCl, pH 7.4.Then, 50 μL p-nitrophenyl butyrate substrate (5 mM) was added to initiate the incubation reaction for 20 min, and the Abs 405 was measured to calculate the inhibition rate.
where Aa is the absorbance obtained without inhibitor (KPHs) and with PL, Ab is the absorbance obtained without (KPHs) and PL, Ac is the absorbance obtained with inhibitor and PL, and Ad is the absorbance obtained with inhibitor and without PL.

Cholesterol esterase inhibition assay
Cholesterol esterase (CEase) inhibition assay was performed following the previous method [33] with slight modification.The KPHs or fractions (50 μL) at different concentrations were mixed with 100 μL sodium phosphate buffer (100 mM, pH 7.0) containing taurocholic acid (5.16 mM) and p-nitrophenyl butyrate substrate (5 mM).Then, 50 μL CEase (5 μg/mL) initiated the reaction system and the Abs 405 after 30 min was measured to calculate the inhibition rate.
Where Aa is the absorbance obtained without inhibitor (KPHs) and with CEase, Ab is the absorbance obtained without (KPHs) and CEase, Ac is the absorbance obtained with inhibitor and CEase, and Ad is the absorbance obtained with inhibitor and without CEase.

Cell viability assay
HepG 2 cells were seeded into 96-well plates (1 × 10 4 /well) overnight and treated with FFAs and/or KPHs for 24 h.After treatment, an additional 10 μL CCK-8 reagent was added to each well and incubated at 37°C for 1 h.The cell viability was measured by Cell Counting Kit-8 (CCK-8) (Beyotime Biotechnology, Shanghai, China).The absorbance at 450 nm was measured to calculate cell viability using the following equation [34] :

Oil red O staining
HepG 2 cells were seeded into 6-well plates (2 × 10 5 /well) overnight and treated with FFAs and/or KPHs for 24 h.Cells were fixed with 4% paraformaldehyde for 10 min at room temperature.Following fixation, the slides underwent a series of carefully executed steps.They were rinsed with distilled water twice and subsequently immersed in 60% isopropanol for a duration of 5 minutes.Subsequently, freshly prepared Oil Red O reagent was applied to the slides and allowed to incubate for precisely 15 minutes.After the incubation period, the slides underwent thorough washing with distilled water until complete removal of Oil Red O dye was achieved.For further staining, hematoxylin stain was applied to the slides for exactly 1 minute, followed by a subsequent rinse with distilled water, adhering to the manufacturer's guidelines (Solarbio, Beijing, China).The slices were acquired images with a Biomicroscope C×31inverted microscope (OLYMPUS, Tokyo, Japan).Orlistat was used as the positive control.

Biochemical parameters assay
HepG 2 cells were seeded into 6-well plates (2 × 10 5 /well) overnight and treated with FFAs and/or KPHs (5 mg/mL and 10 mg/mL) for 24 h.Orlistat (10 µg/mL) was used as the positive control.Intracellular levels of TG and TC were measured by commercial kits according to the manufacturer's instructions, and normalized by quantification of total protein concentration.

In silico analysis of KPHs
The PeptideRanker Server (http://distilldeep.ucd.ie/PeptideRanker/,accessed February 22, 2022) based on the N-to-1 neural network returns the predicted bioactivity probability ranking report of KPHs, and the peptides predicted over the threshold of 0.8 were labeled.Subsequently, the ToxinPred server (https://webs.iiitd.edu.in/raghava/toxinpred/design.php,accessed February 23, 2022) was used to predict the toxicity of given peptides to help identify the particular peptide without development restrictions.Pepsite2(http:// pepsite2.russelllab.org/help,accessed February 25, 2022) was used to validate the interactions between PL (PDB code: 1ETH) or CEase (PDB code: 1AQL) and the above output peptide sequences to identify the promising peptide fragments.Predicted matches with statistical significance (P < .05)were highlighted as the potential binding sites and the candidate peptides that interacted with amino acids at the catalytic and substrate binding sites.Within the limits of the Pepsite2 server for input peptide sequences, the sequences containing more than ten amino acid residues were bisected to obtain p-values.

Statistical analysis
Results were presented as mean ± standard deviation.Data were analyzed by one-way variance method (ANOVA), followed by multiple tests using the Prism 9.0 software (GraphPad, La Jolla, CA) to assure the significant difference at P < .05.

Preparation of lipid-lowering hydrolysates
Katsuwonus pelamis dark muscle hydrolyzates (KPHs) were prepared using six commercially available proteases.The enzyme inhibition assay and cellular TG accumulation assay were performed to screen the potent lipid-lowering peptides fraction (Table 2 & Figure 1).The IC 50 values of the protein hydrolyzates produced by alkaline proteinase (KPHs-AL), neutral protease (KPHs-NP) and papain (KPHs-PP) on PL were 5.14, 7.55 and 14.56 mg/mL, respectively, which were significantly lower than other hydrolyzates.In addition, KPHs-AL showed the lowest IC 50 value (8.51 mg/mL) for inhibiting CEase activity.KPHs-AL exhibited the most competitive TG content reduction performance in the hepatic steatosis cell model (The results of FFA-induced steatosis in HepG 2 cells were shown in Figure S1), which was 38.13 ± 0.03%, suggesting that Alkaline protease may be a high-quality manufacturer of lipid-lowering peptides from Katsuwonus pelamis dark muscle.Alkaline proteinase has been reported to have differential effects on hydrolyzate activity due to its different cleavage patterns. [33]Several reports demonstrated that alkaline proteinase hydrolyzates obtained from porcine liver or scallops had better functional activity than other commercial enzymatic hydrolyzates [35] and. [36]ith regards to this, the hydrolysis parameters of KPHs-AL were comprehensively designed and analyzed based on the five-factor-four level orthogonal experiment to obtain a more reliable combination (Table 3).The results of range (R) analysis revealed that the levels change of temperature (A) and hydrolysis time (D) had a great influence on the lipid-lowering performance, which were the optimal factors, but the effects of liquid to material (C), pH (B) and enzyme concentration (E) were all insignificant.In terms of peptides content, the influential order of the four factors within the extreme value R was A > D > C > B > E, which also indicated the significant effect of factors A and D. In combination with the result of decrease rate of TG content and peptides content, the superior parameters combination was determined to A 1 B 1 C 4 D 2 E 3 , which were 45°C, pH 8, liquid to material 5:1 mL/g, hydrolysis time 4 h and enzyme concentration 5000 U/g.

Separation and characterization of lipid-lowering hydrolysates
The in vitro inhibitory effects of KPHs-AL fractions on PL or CEase were elucidated with the concentration-response curve Figure 2(a-b).The different Mw fractions of KPHs-AL presented a concentration-response percentage inhibition activity, and KPHs-AL-Ⅳ (3 kDa-0.15kDa) exhibited a general inhibition percentage, which was 88.07% to PL and 78.33% to CEase at the concentration of 25 mg/mL, which was comparable to KPHs-AL-III (5 kDa-3 kDa).The percentage inhibition activity of KPHs-AL-I (>10 kDa) and KPHs-AL-II (10 kDa-5 kDa) is lower than the KPHs-AL and other fractions.However, the KPHs-AL fraction showed the best half inhibitory concentration value of 5.14 mg/mL and 8.51 mg/mL compared with different Mw fractions.In addition, a significant decrease (P < .05vs Model) in TG content was also observed in FFA-induced HepG 2 cells with KPHs-AL intervention, which was only about 45.06% of the lipid accumulation in the FFA model group Figure 2(c).Heretofore, KPHs-AL fractions failed to present a significant cytotoxic effect below the test range of 5-30 mg/mL Figure 2(d) The Mw distribution of KPHs-AL was determined by the peak area normalization method using the calibration equation fitted to the standard Mw mix (Figure 3).For KPHs-AL, the percentage of Mw in the range of 3 kDa-0.15kDa was the highest (91.32%) among other Mw fractions, which were 1.6% (>10 kDa).Previous research indicated that low Mw peptides generated by enzymatic hydrolysis had shown great potential for bioactivity. [37]In general, the absorption of high molecular weight (HMW, >1,000) peptides within the human body is relatively limited.The human body possesses the capacity to absorb peptides of various molecular sizes.However, smaller peptides are generally more readily absorbed.Nevertheless, it is crucial to note that the absorption characteristics of peptides are influenced by multiple factors beyond just their molecular weight.Factors such as peptide structure, gastrointestinal conditions, and the presence of synergistic compounds can also impact their absorption within the body.Therefore, further research is warranted to gain a deeper understanding of the absorption and metabolism of KPHs-AL within the human body.These studies will contribute to a comprehensive assessment of the feasibility of KPHs-AL as a potential therapeutic approach or dietary supplement The amino acid composition of KPHs-AL is summarized in Table 4, and 16 amino acids including glutamic acid (11.26 ± 0.05 g/100 g pro.), lysine (5.59 ± 0.01 g/100 g pro.) and aspartic acid (5.24 ± 0.04 g/100 g pro.) were detected.Furthermore, the branched-chain amino acids (BCAAs, Ile, Leu and Val) accounted for 17.24% of total amino acids, and glutamic and aspartic acids accounted for 18.39% and 8.49% in KPHs-AL, respectively.Previous reports indicated that the BCAAs, [34] Asp and Glu [15] in the peptide sequences were considered to have regulatory functions in lipid metabolism.The study conducted by D. Ruan, et al showed that Dietary Lys supplementation influenced the activities of coenzymes in connection with lipid metabolism and promoted fatty acid β-oxidation. [34]It had been reported that leucine (Leu) decreased fat deposition .

Inhibitory effects of KPHs-AL on lipid accumulation
Combined with the Oil Red O staining results, the co-treatment of KPHs-AL significantly decreased the lipid droplet deposition in FFA-induced HepG 2 cells Figure 4(a).Simultaneously, the effects of KPHs-AL on TC and TG in steatosis HepG 2 cells were determined, as shown in Figure 4. Orlistat co-treated (positive control) could extremely decreased  intracellular TC (P < .01)and TG (P < .01)levels.Compared with the model group, similar effects to positive control were observed in KPHs-AL (5 mg/mL and 10 mg/mL), which significantly attenuates FFA-induced lipid accumulation in HepG 2 cells.

Identification and molecular interaction of KPHs-AL
Forasmuch as the potent inhibitory activities in vitro and cellular lipid-lowering effect, the fraction was used for database search and de novo sequencing of peptide sequences.In the sequencing phase, the peptide sequences with an average local confidence (ALC) value of ≥ 80% were screened for subsequent procedures (data not shown).Overall, a total of 1938 nontoxic and potentially bioactive  peptides were identified from the KPHs-AL fraction (Table S1), as PeptideRanker predicted probability >0.8 and ToxinPred SVM score <0.Moreover, the protein-peptide interaction in-silico prediction of potential bioactive peptides was performed by the Pepsite 2 server.Table S2-5 shows the P-values and interacting acceptor residues of the best match between the identified peptides in KPHs-AL fraction and PL or CEase.The peptide fragments with foul residue numbers were equally divided and submitted to the server, with statistical significance in both portions retained.Among the identified peptides, 1845 and 662 peptide fragments were confirmed to bind with PL or CEase (P < .05).In particular, both portions of 58 and 11 peptide fragments (up to 16 residues) tend to interact with enzymes.Apart from the high binding likelihood, interactions with specific binding sites of PL (PDB code: 1ETH) were taken into consideration, as these activity residues were regarded as potential indicators of inhibitory properties.The interaction of the active site maintained almost all peptide fragments that significantly binding to PL, and 1766 peptide fragments could bind to 3 catalytic sites (SER153、PHE216 and HIS264) simultaneously.Hence, these peptides were proposed to be potential PL inhibitors.Previously study showed that many bioactive peptides involved in the binding of CEase active sites, including catalytic triads (SER194, ASP320, andHIS435) and oxygen anion holes (GLY107, ALA108, and ALA195), may inhibit the process of cholesterol ester catalytic hydrolysis to release free cholesterol and reduce dietary cholesterol digestion/absorption. Next, the peptide fragments with high interaction likelihood (P < .05)were evaluated with putative CEase binding profiles.However, only 38 peptide sequences were found to be characterized by interaction residues favorable for the development of reductase inhibitors and showed reasonable interactions with the above active sites.
Consequently, the peptide fragments with favorable binding profiles with PL and CEase were cross-checked, emphasized, and shown in Table 5 & 6.Selected peptide sequences were confirmed to have a high statistical significance of binding (P < .05),further supported by the abundant interactions with specific binding sites.In terms of interaction with PL, the association of the catalytic triad with peptide residues predominated, while other peptide sequences, including AFF, CTF, FAL, FFF, FTF, LAF and VFF, also cooperated with SER153 and PHE216.As shown in Table 6, it was found that peptide residues frequently bound on CEase include catalytic residues SER194 and HIS435, and oxygen anion hole residues GLY107 and ALA108, but no interaction of peptide sequences with ASP320 was noted.Selected peptide sequences included a number of Phe (F), Leu (L) and Pro (P), as seen in Table 5 and 6.Hydrophobic residues were considered involved in restricting the functionality of cholesterol esterase and lipase.The presence of hydrophobic residues at the N-terminus may have followed this assumption.

Conclusion
In the present study, the qlkaline proteinase protein hydrolyzate from skipjack tuna dark muscle may affect intracellular lipid accumulation by interacting with PL and CEase.Moreover, 38 peptide fragments were identified and predicted as potential inhibitors of PL and CEase that can bind to multiple active sites on key metabolic enzymes.These potential peptides could be candidates for further activity validation and mechanistic studies during more challenging in vitro and animal experiments.This research not only enhances our understanding of the economic potential of utilizing dark muscle from fishery processing by-products but also contributes valuable insights into novel peptides with the capacity to modulate lipid metabolism-related enzymes.These peptide sequences, with their inhibitory potential, will undergo comprehensive investigation in subsequent studies.The subsequent research stemming from this study is expected to build upon these significant findings and explore various aspects that can enhance our understanding and potentially unlock practical applications of the alkaline proteinase protein hydrolyzate from skipjack tuna dark muscle.

Figure 2 .
Figure 2. Effects of KPHs-AL fractions on PL (A) or CE (B) inhibition, cellular TG contents (C) and cell viability (D) in HepG2 cells.# P < .05 vs.Control group.

Table 2 .
Inhibitory effects of Katsuwonus pelamis dark muscle hydrolyzates (KPHs) prepared using different commercial proteases on PL and CEase.

Table 5 .
Potential binding sites between multienzyme inhibitory peptides and PL using PepSite 2 server and their P-values.

Table 6 .
Potential binding sites between multienzyme inhibitory peptides and CEase using PepSite 2 server and their P-values.