Discovery of ACE Inhibitory Peptides Derived from Green Coffee Using In Silico and In Vitro Methods

Inhibition of angiotensin-I converting enzyme (ACE) is an important means of treating hypertension since it plays an important regulatory function in the renin-angiotensin system. The aim of this study was to investigate the ACE inhibitory effect of bioactive peptides from green coffee beans using in silico and in vitro methods. Alcalase and thermolysin were employed to hydrolyze protein extract from coffee beans. Bioactive peptides were identified by LC-MS/MS analysis coupled with database searching. The potential bioactivities of peptides were predicted by in silico screening, among which five novel peptides may have ACE inhibitory activity. In vitro assay was carried out to determine the ACE inhibitory degree. Two peptides (IIPNEVY, ITPPVMLPP) were obtained with IC50 values of 57.54 and 40.37 μM, respectively. Furthermore, it was found that two inhibitors bound to the receptor protein on similar sites near the S1 active pocket of ACE to form stable enzyme–peptide complexes through molecular docking, and the Lineweaver–Burk plot showed that IIPNEVY was a noncompetitive inhibitor, and ITPPVMLPP was suggested to be a mixed-type inhibitor. Our study demonstrated that two peptides isolated from coffee have potential applications as antihypertensive agents.


Introduction
Hypertension or high blood pressure is a major risk factor for several diseases, particularly stroke, heart disease, and chronic kidney disease [1].It is the leading cause of morality globally, affecting more than 1 billion people with the disease, and is currently experiencing a younger decline.The burden of high blood pressure is disproportionately high in low-and middle-income countries, where two-thirds of cases are present.
Angiotensin-I converting enzyme (ACE), a zinc metallopeptidase, plays a vital role in elevating blood pressure converting angiotensin-I into angiotensin-II, a powerful vasoconstrictor [2].Due to the evident side effects of synthesized ACE inhibitory medicine (captopril, enalapril), natural ACE inhibitors have become a hot topic of research interest.Therefore, bioactive peptides from natural sources such as animal products, plants, marine resources, microorganisms, and their fermentation products become safe and effective ACE inhibitors to treat hypertension.Animal products such as eggs and milk are rich in protein and are important sources of bioactive peptides.Some oligopeptides with ACE inhibitory activity were found in the egg whites of chicken [3] and ostrich eggs [4].Protein-rich plant-derived peptides such as beans, grains, and seeds have also been extensively studied in recent years [5][6][7].Moreover, a variety of bioactive peptides obtained from the protein of marine organisms including fish [8], mollusks [9,10], seaweed [11,12], and microbial fermentation products [13,14] have been confirmed to have an ACE inhibitory function.With the advance of peptidomics, it has been demonstrated that plant polypeptides exhibit different functions and may have significant bioactivities like being anti-microbial, anti-oxidative, anti-inflammatory, and antihypertensive [15].The inhibitory effect of peptides on ACE depends on the molecular weight of the active peptide and the hydrophobicity of the three amino acid residues at the C-terminal [16].Experimental results also suggested that the more hydrophobic amino acids the ACE inhibitory peptide contains, the easier it is to enter the hydrophobic region of the enzyme active center, thereby interacting with the amino acid residues to inhibit its activity [17].The main method for obtaining plant peptides is to hydrolyze proteins with different proteases [18].Alcalase and thermolysin are two commercialized microbial proteases, which are commonly used to hydrolyze proteins from food sources to prepare active oligopeptides.In the current study, these two enzymes were used due to their characteristics of enzymatic cleavage of hydrophobic amino acids and broad action sites.
Coffee is one of the most valuable traded commercial drinks in the world, especially Coffea arabica, which occupied over 55% of the world's coffee production in 2021 [19].Longterm research studies have shown that caffeine, chlorogenic acids (CGAs), and trigonelline are the main compounds with diverse bioactivities [20].In recent years, a variety of novel diterpenes have been demonstrated existing in green Arabica coffee beans in addition to cafestol and kahweol [21,22].Coffee beans, the "fruits and seeds of plants of the genus Coffea" ref. [23], containing about 10% protein [24] which plays an important role in the growing stage, can be considered as a natural source of bioactive peptides.Peptides from green coffee beans have been studied using simulated gastrointestinal digestion through in vitro hydrolysis and in silico prediction [25].Although both computer predictions and experimental results indicated that green coffee protein hydrolysates may have biological activities, there was limited information about the identification of peptides.Peptides isolated from spent coffee grounds (SCGs) have also been studied for their antioxidant and ACE inhibitory activity using enzymatic hydrolysis [26] and fermentation by bacteria [27].The Maillard reaction occurs during the roasting step, resulting in denaturation of protein and generation of melanoidins, which might cause interference during extraction and in vitro evaluation of activity [26].So far, there are few studies in the literature describing the oligopeptides with an ACE inhibitory effect obtained from the protein hydrolysates of green coffee.As a result, it is of great significance to discover bioactive peptides and investigate their bioactivity and inhibition mechanism.
This research was to explore active peptides with an ACE inhibitory effect in green coffee beans using the method of peptidomics.Identification, screening, and characterization of these active peptides were conducted, and molecular docking was also used to explore the kinetics and mechanism of inhibition.

Chemicals and Samples
The green coffee beans of coffea arabica were purchased from Baoshan, Yunnan Province of China.

Protein Extraction
The protein extraction procedure was similar to that described by Pérez-Míguez et al. [28] with some minor modifications; 25 g of grounded green coffee beans were defatted four times with 500 mL petroleum ether; 1 g of flour was extracted with 100 mL of 100 mM Tris-HCl buffer (pH 7.5) containing 1% SDS and 0.25% DTT under ultrasonic for 15 min at room temperature.The mixture was centrifuged at 4000 r/min at room temperature for 15 min and supernatants were collected.Protein in the supernatant was precipitated with 100 mL of cold acetone at 4 • C for 12 h and collected by centrifugation for 20 min at room temperature.The precipitate was dialyzed with distilled water overnight followed by lyophilization and stored at −20 • C.

Protein Digestions
Protein extracts were dissolved in the 5 mM phosphate buffer (pH 8.0) at the concentration of 2.5 mg/mL and submitted to enzymatic hydrolysis with alcalase (10,000 U/g substrate) and thermolysin (0.05 g enzyme/g substrate) separately at 50 • C for 2 h with slight stirring, according to the proposed procedure of Esteve et al. [29] The digestions were stopped by increasing the temperature to 90 • C for 20 min, and supernatants were collected by centrifugation at 10,000 g and 4 • C for 10 min.
Identification of the desalted sample was then conducted by LC-MS/MS using a Triple TOF 6600 + mass spectrometer in positive ion mode with a OptiFlow Turbo V ion source, coupled to a NanoLC 400 (Ekspert); 8 µL of hydrolysate was loaded on a trap Nano LC precolumn at 10 µL/min for 3 min.Thereafter, the sample was loaded at 5 µL/min onto a Micro LC Column (3 µm, ChromXP C18, 120 Å, 0.3 × 150 mm) at the column temperature of 40 • C for the separation of peptides.Mobile phase A was 2% ACN containing 0.1% formic acid (v/v), and mobile phase B was 98% ACN containing 0.1% formic acid (v/v).Elution gradient employed for separation was: 3-22% B in 40 min, 22-35% B in 10 min, 35-80% B in 1 min, 80% B for 3 min, 80-3% B in 1 min, and 3% B for 5 min.The following settings were applied: nebulizer (gas 1) 25 psi, heater (gas 2) 50 psi, curtain gas 35 psi, and ionization temperature 150 • C. TOF-MS scan range was set at 350 to 1500 m/z, followed by an MS/MS scan at 100 to 1500 m/z.

In Silico Prediction and Analysis of Bioactive Peptides
The bioactivities of all identified peptides were evaluated using in silico methods.The ACE inhibitory capacity of peptides was predicted using the AHTpin platform (https: //webs.iiitd.edu.in/raghava/ahtpin/(accessed on 1 March 2023)) in accordance with the support vector machine (SVM) classification model [33].Peptides with an SVM score greater than 0 were regarded as antihypertensive peptide inhibitors.In addition, the physiochemical properties (hydrophobicity and isoelectric point) of peptides were analyzed by the PepDraw server (http://pepdraw.com/(accessed on 28 February 2023)).

Peptide Synthesis and ACE Inhibitory Activity
The peptides predicted with potential antihypertensive capacity were synthesized by GenScript.The ACE inhibitory capacity of hydrolysates and selected peptides were measured as previously described with minor modifications [34]; 400 µL of samples (hydrolysates, synthesized peptides dissolved in the HEPES-NaCl buffer at a concentration of 1.25 mg/mL, 125 µM separately) were mixed with 100 µL of ACE (0.2 U/mL).After incubation of the mixture at 37 • C for 10 min, 500 µL of FAPGG was added and the absorbance value at 340 nm was detected instantly by a FlexStation 3 and recorded as A 0 min .After being kept at 37 • C for 30 min, the absorbance value at 340 nm was detected again and recorded: where A 0 min is the initial absorbance at 340 nm; A 30 min is the absorbance at 340 nm after keeping at 37 • C for 30 min, and ∆A is the reduction in the absorbance: A 0 min -A 30 min .
Referring to the methods in the literature, and combined with the above in vitro activity assay results, the design of the inhibition pattern test is as follows: various concentrations (0, 0.2, 0.4 mM) of ACE inhibitory peptides were added to each reaction system; meanwhile, the concentration of FAPGG were set as 0.25, 0.5, 1, 2 mM.According to the Michaelis-Menten kinetic equation, the Michaelis-Menten constant (K m ) and the maximum reaction rate (V max ) were calculated from the Lineweaver-Burk plot to determine the inhibition pattern of peptides.All evaluations were performed in triplicate.

Gastrointestinal Digestion Simulation Study
In vitro gastrointestinal digestion simulation was performed according to the method described by Fan et al. [3] with minor modifications.Two peptides (CP2, CP3) were dissolved in ultrapure water at 1.0 mg/mL, respectively.The solutions were adjusted to pH 2.0 with 0.1 HCl.Then, pepsin (1% E/S, w/w) was added to peptide solutions and the mixtures were incubated at 37 • C for 2 h.Half of the two digests were terminated in digestion by increasing the temperature to 90 • C.After the other half were adjusted to pH 7.5 with 0.1 NaOH, pancreatin (1% E/S, w/w) was added and the systems were digested at 37 • C for 2 h.The reactions were then stopped by boiling for 10 min.Subsequently, solvents were centrifuged at 10,000 rpm for 10 min.The supernatants were lyophilized to collect and analyzed using Agilent 6545 series Q-TOF mass spectrometer (Agilent, Waldbronn, Germany).The ACE inhibitory activity of the digested peptide solutions was determined under the same conditions.

Molecular Docking
Molecular docking was conducted using Autodock Vina software (version 1.5.7) to further study the mechanism of interaction between ACE (receptor) and peptides (ligands) according to methods previously reported with slight modification [35].The crystal structure of ACE was obtained from the RSCB protein data bank (https://www.rcsb.org/(accessed on 13 March 2023)) with the PDB code of 1O86.The structures of selected peptides were generated by ChemDraw.The number of points in the X-, Y-, and Z-dimensions were set as 75, 75, and 75, and spacing was set as 1Å.PyMol software (version 2.5.4) was used to visualize the pose of ligand from the docked complex.

Peptide Identification
The hydrolysates were identified with MaxQuant software (version 2.3.0.0) for each mass spectral data generated by LC-MS/MS.The results showed that 55 and 22 peptides were identified from alcalase and thermolysin hydrolysates, respectively.The identified peptides ranged in length from 7 to 15 amino acids and had molecular masses ranging from 600 to 1400 m/z.Tables S1 and S2 summarized the peptides after digestion, along with their experimental data.Figure 1 shows, as examples, the MS/MS spectra obtained for two of the identified peptides: GLPSGGAPSGY (Gly-Leu-Pro-Ser-Gly-Gly-Ala-Pro-Ser-Gly-Tyr) in the alcalase hydrolysate and ITPPVMLPP (Ile-Thr-Pro-Pro-Val-Met-Leu-Pro-Pro) in the thermolysin hydrolysate.It is worth noting that no common peptide fragments were found in the coffee protein hydrolysate after enzymatic hydrolysis by two proteases.mass spectral data generated by LC-MS/MS.The results showed that 55 and 22 peptides were identified from alcalase and thermolysin hydrolysates, respectively.The identified peptides ranged in length from 7 to 15 amino acids and had molecular masses ranging from 600 to 1400 m/z.Tables S1 and S2 summarized the peptides after digestion, along with their experimental data.Figure 1 shows, as examples, the MS/MS spectra obtained for two of the identified peptides: GLPSGGAPSGY (Gly-Leu-Pro-Ser-Gly-Gly-Ala-Pro-Ser-Gly-Tyr) in the alcalase hydrolysate and ITPPVMLPP (Ile-Thr-Pro-Pro-Val-Met-Leu-Pro-Pro) in the thermolysin hydrolysate.It is worth noting that no common peptide fragments were found in the coffee protein hydrolysate after enzymatic hydrolysis by two proteases.It has been reported that alcalase has a wide range of applications in food production and processing due to broad selectivity and specificity of the cleavage site of alcalase [36].The results have shown that among the products hydrolyzed by alcalase, the C-terminal amino acid residues of peptides were mostly Gln, Tyr, Met, Lys, and Leu.In addition, it was also found that in a certain proportion of peptide fragments, the hydrophobic amino acid was close to the N-terminal, which is very common in bioactive peptides [37].In contrast, thermolysin hydrolysate presented a smaller number of peptides, only 22 peptides were identified, and mostly produced large fragments under this hydrolysis condition.The results in Table S2 showed that a great number of peptides have Leu, Phe, Ile, and Val at the N-terminal position, which was consistent with literature [38], and the presence of large bulky amino acids in these peptide sequences suggested potential antihypertensive capacity.

In Silico Screening
Peptides with a score greater than 60 indicating the degree of matching between the secondary mass spectrum and their databases, were submitted to analyze hydrophobicity and isoelectric point and predict potential bioactivity using the in silico method.The physiochemical properties analyzed by PepDraw server and ACE inhibition prediction score of peptides hydrolyzed by alcalase are shown in Table 1, and those of peptides hydrolyzed by thermolysin are shown in Table 2. Light, dark shades correspond to peptides with SVM score higher than 0.0, 1.0, respectively.
The AHTpin platform was used to screen antihypertensive peptides in accordance with the SVM classification models.With the existence of the AHTpin platform, the computational prediction of potential antihypertensive peptides from the hundreds of thousands of peptide sequences can be performed in a short time.The 71 identified peptides were entered to predict probability in ACE inhibition; 16, 15 peptides were authenticated as antihypertensive peptides because their SVM scores were higher than 0 in alcalase, thermolysin hydrolysates, respectively [33].The higher the score, the higher the ACE inhibitory activity of peptides; when the SVM score is > 1, the peptide may have a relatively higher potential with ACE inhibitory activity [39].According to the result, five coffee bean peptide sequences with an SVM score of > 1.0 were chosen.The five selected sequences GLPSGGAPSGY, IIPNEVY, ITPPVMLPP, VLETPDGPL, and VKNPNPIPIP were named CP1, CP2, CP3, CP4, and CP5.

In Vitro Evaluation of ACE Inhibitory Capacity
The selected monomer peptides were synthesized by the solid-phase synthesis method with a purity of ≥95%.The ACE inhibitory capacity of five synthesized peptides predicted by AHTpin with potential antihypertensive activity were determined.The inhibition rate (%) of two enzyme hydrolysates and five peptides separately are shown in Figure 2.
Thermolysin was the enzyme yielding peptides with the higher ACE inhibitory capacity, and the hydrolysate of alcalase also had considerable ACE inhibitory activity, according to results.Generally, the antihypertensive peptides contained between 2 and 12 amino acids, and presented a certain amount of hydrophobic residues [40,41].Compared with captopril, the ACE inhibition rate of thermolysin hydrolysate reached 73.3%, indicating that there are very promising inhibitors among them, and these inhibitors are completely safe without any side effects.
The selected monomer peptides were synthesized by the solid-pha method with a purity of ≥95%.The ACE inhibitory capacity of five synthesiz predicted by AHTpin with potential antihypertensive activity were determi hibition rate (%) of two enzyme hydrolysates and five peptides separately a Figure 2. Thermolysin was the enzyme yielding peptides with the higher ACE in pacity, and the hydrolysate of alcalase also had considerable ACE inhibitory cording to results.Generally, the antihypertensive peptides contained betw amino acids, and presented a certain amount of hydrophobic residues [40,41 with captopril, the ACE inhibition rate of thermolysin hydrolysate reached cating that there are very promising inhibitors among them, and these inhibit pletely safe without any side effects.
The ACE inhibitory capacity of CP 1-5 was tested at a concentration of results showed that five peptides exhibited various degrees of inhibitory ca The ACE inhibitory capacity of CP 1-5 was tested at a concentration of 50 µM.The results showed that five peptides exhibited various degrees of inhibitory capacity from 6.0 to 53.3%.The inhibition rates of the two peptides of CP2 (IIPNEVY) and CP3 (ITPPVMLPP) were significantly higher than those of the other three peptides, which were 49.3% and 53.3%, respectively.Their IC 50 values were 57.54 ± 7.78 µM and 40.37 ± 5.56 µM, respectively.Among them, the inhibition rate of CP4 (VLETPDGPL) was the lowest, only 6.0%.Five peptides obtained by computational screening all showed certain activities, suggesting that the prediction results have reference value.

Inhibition Kinetics of Bioactive Peptides
The inhibition models of ACE inhibitory peptides mainly include competitive inhibitors, noncompetitive inhibitors, and mixed-type inhibitors and their inhibition mechanisms are different.To determine the inhibitory kinetic pattern of the ACE inhibitory peptides obtained from green coffee protein, Lineweaver-Burk plots were constructed, and the results are shown in Figure 3.As the concentration of IIPNEVY (CP2) increased, the Michaelis-Menten constant (K m ) remained unchanged, while the V max decreased, which is characteristic of the noncompetitive inhibition pattern indicating that this peptide bound to essential sites outside the active part of enzyme.The reported noncompetitive inhibitors include VFDGVLRPGQ from rice bran protein [34], QLDL from the mycelia of Ganoderma Lucidum [42], and VGLPPNSR and QAGLSPVR from tilapia skin gelatin [43], etc.As for the peptide ITPPVMLPP (CP3), Figure 3B showed that the fitted curve did not intersect the x-axis at the same point and the K m values of control and various inhibitory peptide concentration were not the same, suggesting that the inhibition of CP3 at different concentrations presented a mixed-type inhibition pattern.The V max values decreased as the peptide concentration increased and the K m increased with the addition of an inhibitory peptide.The result was similar to the ACEI peptides YGIKVGYAIP from palm kernel cake protein hydrolysates [44].
Ganoderma Lucidum [42], and VGLPPNSR and QAGLSPVR from tilapia skin gelatin [43], etc.As for the peptide ITPPVMLPP (CP3), Figure 3B showed that the fitted curve did not intersect the x-axis at the same point and the Km values of control and various inhibitory peptide concentration were not the same, suggesting that the inhibition of CP3 at different concentrations presented a mixed-type inhibition pattern.The Vmax values decreased as the peptide concentration increased and the Km increased with the addition of an inhibitory peptide.The result was similar to the ACEI peptides YGIKVGYAIP from palm kernel cake protein hydrolysates [44].

Molecular Docking Studies
To further explain the inhibitory mechanism between the identified peptides and ACE, the docking study between CP2 or CP3 and ACE was performed.Prior research have suggested that the binding pockets related to the activity of ACE are mainly distributed in S1(Ala354, Tyr523, and Glu384), S2(Gln281, Lys511, His353, Tyr520, and His513), and S1'(Glu162) [45].However, the docking results showed that CP2 or CP3 did not form any interactions with amino acid residues of ACE active pockets.As shown in Figure 4, two peptides bound to the receptor protein in the same cavity at the upper end of ACE

Molecular Docking Studies
To further explain the inhibitory mechanism between the identified peptides and ACE, the docking study between CP2 or CP3 and ACE was performed.Prior research have suggested that the binding pockets related to the activity of ACE are mainly distributed in S1(Ala354, Tyr523, and Glu384), S2(Gln281, Lys511, His353, Tyr520, and His513), and S1'(Glu162) [45].However, the docking results showed that CP2 or CP3 did not form any interactions with amino acid residues of ACE active pockets.As shown in Figure 4, two peptides bound to the receptor protein in the same cavity at the upper end of ACE close to active pocket S1 and tended to combine with the similar amino acid residues.Both peptides produced an H-bond very close to the critical residue Ala354 in an ACE active pocket.Moreover, CP3 bound with Arg522 in pocket S1 exhibiting a higher inhibition rate according to the in vitro assay.We hypothesized that the ability of CP2 was to generate interactions with nonactive sites, resulting in a distortion of the catalytic configuration, which is consistent with the inhibition kinetic pattern of bioactive peptides as ACE noncompetitive inhibitors.When CP3 bound to the receptor protein, it formed more complex interactions with the amino acid residues near the active sites (Ser355, Arg522), suggesting it may be a potent mixed-type inhibitor.

In Vitro Gastrointestinal Digestion Stability
As ACE inhibitory peptides exert physiological regulatory functions after oral administration, there is a risk of being degraded by proteases in the gastrointestinal tract.Therefore, the stability of the two inhibitory peptides against gastrointestinal digestion was analyzed using LC-MS characterization and variations in the ACE inhibition rate.The MS results of CP2 after pepsin hydrolysis were consistent with that of control (Figure 5A,B).After successive hydrolysis by pepsin and pancreatin, new peptide sequences were detected (Figure 5C).Compared with before digestion, the ACE inhibition rate of CP2 decreased after being digested by pepsin, and then after digestion by pancreatin, its inhibition rate decreased significantly (Table 3).The result was consistent with the LC-MS analysis, indicating that the CP2 is well resistant in gastric digestion, but may be affected by pancreatin.Among the products from in vitro digestion (pepsin, pepsin-pancreatin) of CP3, new peptides were determined.Only a small amount of CP3 was degraded in the digestion by pepsin (Figure 5E), and the result of the ACE inhibition assay also suggested that CP3 was partially degraded during pepsin hydrolysis as the inhibition rate did not change significantly.Great changes in the MS analysis and ACE inhibitory activity occurred after incubation with pepsin-pancreatin.Thus, pancreatin may be hydrolyzing the CP3 and its pepsin digestion product into smaller peptide fragments with a higher ACE inhibitory effect.The result of this simulation has suggested that CP3 may be a pro-drug peptide which could be converted to true inhibitors by gastrointestinal proteases, resulting in increased activity [46].
Foods 2023, 12, x FOR PEER REVIEW 10 of 14 close to active pocket S1 and tended to combine with the similar amino acid residues.Both peptides produced an H-bond very close to the critical residue Ala354 in an ACE active pocket.Moreover, CP3 bound with Arg522 in pocket S1 exhibiting a higher inhibition rate according to the in vitro assay.We hypothesized that the ability of CP2 was to generate interactions with nonactive sites, resulting in a distortion of the catalytic configuration, which is consistent with the inhibition kinetic pattern of bioactive peptides as ACE noncompetitive inhibitors.When CP3 bound to the receptor protein, it formed more complex interactions with the amino acid residues near the active sites (Ser355, Arg522), suggesting it may be a potent mixed-type inhibitor.

In Vitro Gastrointestinal Digestion Stability
As ACE inhibitory peptides exert physiological regulatory functions after oral administration, there is a risk of being degraded by proteases in the gastrointestinal tract.Therefore, the stability of the two inhibitory peptides against gastrointestinal digestion was analyzed using LC-MS characterization and variations in the ACE inhibition rate.The MS results of CP2 after pepsin hydrolysis were consistent with that of control (Figure 5A,B).After successive hydrolysis by pepsin and pancreatin, new peptide sequences were detected (Figure 5C).Compared with before digestion, the ACE inhibition rate of CP2 decreased after being digested by pepsin, and then after digestion by pancreatin, its inhibition rate decreased significantly (Table 3).The result was consistent with the LC-MS analysis, indicating that the CP2 is well resistant in gastric digestion, but may be affected by pancreatin.Among the products from in vitro digestion (pepsin, pepsin-pancreatin) of CP3, new peptides were determined.Only a small amount of CP3 was degraded in the digestion by pepsin (Figure 5E), and the result of the ACE inhibition assay also suggested that CP3 was partially degraded during pepsin hydrolysis as the inhibition rate did not  change significantly.Great changes in the MS analysis and ACE inhibitory activity occurred after incubation with pepsin-pancreatin.Thus, pancreatin may be hydrolyzing the CP3 and its pepsin digestion product into smaller peptide fragments with a higher ACE inhibitory effect.The result of this simulation has suggested that CP3 may be a pro-drug peptide which could be converted to true inhibitors by gastrointestinal proteases, resulting in increased activity [46].

Conclusions
In the current study, traditional methods coupled with emergent peptidomics methods were employed to discover potential antihypertensive peptides in green coffee.Five peptides derived from green coffee that may have ACE inhibitory activity were predicted through in silico screening.However, subsequent in vitro activity tests showed that the experimental results were slightly different from the prediction results.The hydrolysates obtained when using proteases both showed considerable ACE inhibitory capacity.Combined with the identification of peptides, it was considered that thermolysin and alcalase cleave peptide bonds near hydrophobic amino acid residues, resulting in peptides containing favorable amino acid residues for antihypertensive activity at the C-terminal.Based on the analysis of the above results, the peptide with the highest experimental inhibition rate was not the one with the highest score, indicating that the result predicted by the webserver AHTpin may be inaccurate.Additionally, some peptides with equivalent or stronger activities were missing, leading to limitations in finding active peptides effectively.In the screening results, two novel oligopeptides, IIPNEVY and ITPPVMLPP, exhibited excellent ACE inhibitory activity with IC 50 values of 57.54 and 40.37 µM.Lineweaver-Burk plots showed that the former peptide was a noncompetitive inhibitor and the latter might act in a mixed-mode type of inhibition.It was found that CP2 bound to the nonactive sites in ACE through molecular docking study, thereby destroying the natural conformation of ACE and affecting its activity.CP3 bound to the receptor protein in the same cavity to cause inhibition of ACE.However, future research should focus on active peptides that may be missed because of the high inhibition rate of coffee protein hydrolysate, and other potential peptides in coffee hydrolysate are worthy of further exploration.Moreover, in vitro gastrointestinal digestion simulation has been conducted to investigate biostability, and the CP3 can be considered to be a so-called pro-drug peptide, while the CP2 requires further experiments to explore its delivery method.In summary, peptides derived from the hydrolysates of Coffea arabica are a promising source for the treatment of high blood pressure.

Figure 2 .
Figure 2. ACE inhibitory capacities of the green coffee bean protein isolate digested thermolysin at the concentration of 500 µg/mL, five screened peptides identified fro drolysates at the concentration of 50 µM, and captopril at the concentration of 0.1 µM

Figure 2 .
Figure 2. ACE inhibitory capacities of the green coffee bean protein isolate digested with alcalase, thermolysin at the concentration of 500 µg/mL, five screened peptides identified from protein hydrolysates at the concentration of 50 µM, and captopril at the concentration of 0.1 µM.

Figure 4 .
Figure 4. Molecular docking of bioactive peptides binding on receptor protein: (A) 3D illustration of the interaction between CP2 and ACE; (B) 3D illustration of the interaction between CP3 and ACE.

Figure 4 .
Figure 4. Molecular docking of bioactive peptides binding on receptor protein: (A) 3D illustration of the interaction between CP2 and ACE; (B) 3D illustration of the interaction between CP3 and ACE.

Table 1 .
Physiochemical properties and SVM score of identified peptides from the alcalase hydrolysates.
Light, dark shades correspond to peptides with SVM score higher than 0.0, 1.0, respectively.

Table 2 .
Physiochemical properties and SVM score of identified peptides from the thermolysin hydrolysates.

Table 3 .
ACE inhibition rate of the peptides following digestion by gastrointestinal proteases.
a Peptides without any treatment.

Table 3 .
ACE inhibition rate of the peptides following digestion by gastrointestinal proteases.