Activity Prediction and Molecular Mechanism of Bovine Blood Derived Angiotensin I-Converting Enzyme Inhibitory Peptides

Abstract

Development of angiotensin I-converting enzyme (ACE, EC 3.4.15.1) inhibitory peptides from food protein is under extensive research as alternative for the prevention of hypertension. However, it is difficult to identify peptides released from food sources. To accelerate the progress of peptide identification, a three layer back propagation neural network model was established to predict the ACE-inhibitory activity of pentapeptides derived from bovine hemoglobin by simulated enzyme digestion. The pentapeptide WTQRF has the best predicted value with experimental IC₅₀ 23.93 μM. The potential molecular mechanism of the WTQRF / ACE interaction was investigated by flexible docking.

About the Authors:

Ting Zhang

Affiliation: Laboratory of Nutrition and Functional Food, Jilin University, Changchun, Jilin, China

Shaoping Nie

Affiliation: State Key Laboratory of Food Science and Technology, Nanchang University, Nanchang, Jiangxi, China

Boqun Liu

Affiliation: Laboratory of Nutrition and Functional Food, Jilin University, Changchun, Jilin, China

Yiding Yu

Affiliation: Laboratory of Nutrition and Functional Food, Jilin University, Changchun, Jilin, China

Yan Zhang

Affiliation: Laboratory of Nutrition and Functional Food, Jilin University, Changchun, Jilin, China

Jingbo Liu

* E-mail: [email protected]

Affiliation: Laboratory of Nutrition and Functional Food, Jilin University, Changchun, Jilin, China

Introduction

Hypertension is a risk factor for cardiovascular diseases including coronary heart disease, peripheral artery disease and stroke. [1, 2] Recently, several food-derived bioactive peptides have been found playing a significant role in decreasing blood pressure. Therefore, more and more attention has been paid to peptides from food sources with antihypertensive activity. [3]

Most of the antihypertension peptides regulate blood pressure by inhibiting the activity of angiotensin I-converting enzyme (ACE, EC.3.4.15.1). ACE is a zinc- and chloride- dependent metallopeptidase, which belongs to the M2 family of zinc metallopeptidases. [4, 5] It converts angiotensin I to angiotensin II (a potent vasoconstrictor) as well as inactivates the vasodilator bradykinin. [6] ACE plays a crucial role in the renin-angiotensin system (RAS), which is well known for its regulation of blood pressure and fluid homeostasis. [7, 8] Nowadays, inhibitors of ACE have been considered as first-line therapy for hypertension. [9, 10] It has been reported that a number of bioactive peptides, which derived from food sources, have ACE-inhibitory activity. Chibuike C. Udenigwe [11] summarized the major approaches in bioactive peptides research as the classical approach, the bioinformatics approach and the integrated approach. This classification is also suitable for ACE-inhibitory peptides.

The classic approach is the most widely used method for the discovery of ACE-inhibitory peptides from food proteins, involving peptides production (solvent extraction, enzyme hydrolysis, and microbial fermentation), purification (membrane-based separation and chromatography techniques) and identification (mass spectrometry methods). There are a number of ACE-inhibitory peptides derived from different food sources and obtained by the classic approach. For instance, ACE-inhibitory peptides derived from soy protein such as DLP, DG, IA, ILAGNQ, FFL, IYLL, VMDKPQG, IFL, WL, TPRVF, YVVFK, PNNKPFQ, EDENNPFYLR, NWGPLV, IPPGVPYWT, VLIVP, LAIPVNKP, LPHF, SPYP and WL, were found in published articles. [12–17] In wheat germ hydrolysates, 16 peptides [18] with the IC₅₀ value of less than 20 μM, composed of 2–7 amino acid residues were identified. And IAP [19] was identified in wheat gliadin hydrolysates. In pork meat hydrolysates digested by in vitro gastrointestinal digestion, 12 peptides were identified. [20] Also, in beef rump (biceps femoris) hydrolysates, Jang and Lee [21] identified VLAQYK.

In order to circumvent some challenges of the classical approach, the bioinformatics approach has been recently applied towards the discovery of ACE-inhibitory peptides encrypted in food proteins. This approach was recently used to study the distribution of ACE-inhibitory peptides within the primary structure of typical food proteins. [22]

Following the identification of bioactive peptides from protein sets by bioinformatics in databases populated following the classical approach, the remainder of the purportedly “inactive” peptides can be analyzed in silico to identify structural patterns that have previously been associated with known bioactivities. [11]

Moreover, the strengths of each approach can be combined as deemed fit to enhance the discovery and use of ACE-inhibitory peptides. Bioinformatics software can be used to simulate proteolytic specificities of enzymes in order to establish the peptide database in silico. [11] Quantitative structure-activity relationship (QSAR) studies are widely undertaken for modeling the bioactivities such as the bioactivity of ACE-inhibitory peptides [22–26] and the sweetness of compounds [27]. Neural network, as a kind of artificial intelligence, has been applied to modeling non-linear systems, simulating the chaos bioprocess and predicting the results. It turns out to have higher modeling accuracy and generalization capacity [28] and becomes a potentially effective tool in modeling the QSAR.

Bovine blood, as a by-product generated in great volume in industrial abattoirs, gives rise to several possibilities for their recovery and use. The use of bovine blood as a food component has been widely reported due in part to their high nutritional value and there are many studies have recently demonstrated that bovine blood proteins can be used to obtain bioactive peptides. This is important because it gives an added value to bovine blood. A number of bioactive peptides released from bovine hemoglobin hydrolysates have been reported. [29–32] It is a highly desirable but difficult task to identify bovine blood derived peptides. In this aspect, prediction model would be a useful technique to highlight potential ACE-inhibitory peptides identification.

It was found that ACE has two homologous domains (the N-domain and the C-domain), each containing an active center. [33] The C-domain of ACE has been proved to be the dominant angiotensin-I converting site, which has a conserved HEXXH zinc-binding motif, for controlling blood pressure and serving cardiovascular functions. [34, 35] If inhibitory peptides occupied the active site of the C-domain of ACE and bound to specific amino acid residues, ACE will lose its activity. Therefore, it is possible to reveal the ACE inactivation mechanisms by analyzing the structural consequences of ACE-inhibitor interactions. [34]

In this study, we established a database of potential pentapeptides derived from bovine hemoglobin by simulated enzyme digestion, and a prediction model of ACE-inhibitory pentapeptides by back propagation neural network (BPNN). The peptide with best predictive value was synthesized and its IC₅₀ of ACE was measured. We also sought to elucidate the potential molecular mechanism of how the peptide with best predictive value exerts its ACE-inhibitory effects by automated molecular docking.

Materials and Methods

Materials and Chemicals

Angiotensin converting enzyme from rabbit lung, hippuryl-L-histidyl-L-leucine, hippuric acid and HPLC grade acetonitrile, trifluoroacetic acid were purchased from Sigma Chemical Co. (St. Louis, Mo., U.S.A.). All the other reagents were analytical grade.

Dataset

A total of 24 unique pentapeptides were collected from published articles with IC₅₀ of ACE from 0.00948μM to 848 μM in Table 1. [16, 36–48] Due to the wide range of IC₅₀ value, the negative logarithm of IC₅₀ (pIC₅₀) was used ranging from-2.928 to 2.023.

[Figure omitted. See PDF.]

Table 1. The IC₅₀ (μM) of ACE of 24 pentapeptides as taken from literature.

https://doi.org/10.1371/journal.pone.0119598.t001

Molecular Descriptors Calculation

Ten molecular descriptors were calculated for analogs of peptides including structural, spatial, thermodynamic and electronic. Structural descriptors included the number of rotatable bond (Rotbond), the number of hydrogen bond acceptor (H bond acceptor), the number of hydrogen bond donor (H bond donor), the number of ring (Ring) and the number of arcomatic ring (AR). [49–51] The molecular fractional polar surface area (PSA) [52] was used as spatial descriptor. The thermodynamic descriptors were taken to describe the hydrophobic character (ALog P: logarithm of the partition coefficient in octanol / water and Log D: logarithm of the partition coefficient in octanol / water in pH 7.4) and refractivity (MR: molar refractivity). [49] The atomic polarizabilities (Apol) were calculated as electronic descriptors. [50] Peptide structures were generated by Accelrys Discovery Studio 3.5 software (Accelrys Inc., San Diego, USA) and the energy was minimized with CHARMm program using steepest descent and conjugate gradient techniques. All descriptors were calculated by calculate molecular properties protocol.

Modeling by Back-Propagation Neural Network (BPNN)

To construct a BPNN, the Levenberg-Marquardt algorithm [53] was used to train the network. With all data normalized, the descriptors and the pIC₅₀ of dataset were introduced as input and output values, respectively. In order to avoid overfitting, the dataset was randomly classified for training, validation, and test sets. To verify the suitable number of nodes in hidden layer and transfer function, the mean square error (MSE) and pearson correlation coefficient (R) of each BPNN model were calculated.

Simulated Enzyme Digestion and BPNN Prediction

Pepsin and trypsin degrade food proteins into peptides, which are critical in human digestion of protein. The online software Peptide Cutter (http://web.expasy.org/peptide_cutter/) can predict potential cleavage sites cleaved by proteases, including pepsin and trypsin, in a certain protein sequence. As hemoglobin is composed of four poly peptide subunits, two α and two β, we used Peptide Cutter to predict the potential pentapeptides of subunit α (accession number, P01966) and β (accession number, P02070), respectively. All pentapeptides were described by descriptors mentioned above. And then, put the normalized descriptors into the best BPNN model as input data to predict the ACE-inhibitory activity.

Peptide Synthesis

The peptide was synthesized (ChinaPeptides Co. Ltd., Shanghai, China) by the solid phase procedure peptide using FMOC protected amino acids synthesis. The purity (95%) and the molecular masses of the peptides were determined by HPLC and mass spectrometry.

Assay for ACE-Inhibitory Activity

The ACE-inhibitory activity of the peptide with best predict value was measured according to Liu, et al. [54]

Docking

The 3-D structure of human ACE was imported from Protein Data Bank (1O8A.pdb, a crystal structure of human ACE). The structure of peptide was generated by DS and the energy was minimized with CHARMm program. Automated molecular docking was performed using the flexible docking tool of DS software in the presence of cofactors (zinc and chloride ions). The binding site with a radius of 15 Å and a coordinates x: 45.0463, y: 38.6842, z: 45.8268. Evaluation of the molecular docking was performed according to the scores and total potential energy in order to obtain the best pose of peptide. The DS software was used to identify the hydrogen bonds as well as the hydrophobic, hydrophilic, electrostatic and coordination interactions between residues present within ACE active site.

Results and Discussion

Structural Features Analysis

Most of the ACE-inhibitory peptides derived from egg white hydrolysates in our lab previously [55–57] were pentapeptides. Therefore, we are interested in ACE-inhibitory pentapeptide and try to build a prediction model based on a pentapeptides database (Table 1). It shows that 75% of the pentapeptides taken from literature have 1 or 2 hydrophobic amino acid residues. A half of these pentapeptides contain aromatic or hydrophobic amino acid residues at the C-terminus, which are beneficial to ACE-inhibitory activity. [24] In addition, 66.67% of these pentapeptides contain hydrophilic amino acid residues in the second amino acid residues from the C-terminus. The aforementioned common structural features can be used to instruct ACE-inhibitory peptides’ design.

The molecular descriptor (Tables 2, 3) calculation is a logic and mathematical procedure which converts the chemical information of a molecule into some useful data. Most of the number ranges of descriptors of Pentapeptides derived from hemoglobin simulated enzyme digestion included in the scope of the reported pentapeptides, except the molar refractivity of GHGAK. The range of MR of reported pentapeptides is 111.291 to 186.874, while the GHGAK is 107.89. Nevertheless, the predicted IC₅₀ of GHGAK is the fourth lowest among the seven hemoglobin pentapeptides, because the BPNN is a nonlinearity model, leading to a nonlinear relationship between the input and the output.

[Figure omitted. See PDF.]

Table 2. Descriptors of 24 peptides.

https://doi.org/10.1371/journal.pone.0119598.t002

[Figure omitted. See PDF.]

Table 3. Pentapeptides derived from hemoglobin simulated enzyme digestion.

https://doi.org/10.1371/journal.pone.0119598.t003

BPNN model Building

There were three layers in the BPNN model (input layer, hidden layer and output layer), while the number of nodes in hidden layer was not certain. So that we trained BPNN models with different number of nodes in hidden layer and different transfer functions to obtain the suitable model structure (Fig. 1, Fig. 2 and S1 Table, S2 Table). For the log-sigmoid & purelin transfer function, the highest MSE appears when the number of hidden layers is twelve. However, in terms of the tan-sigmoid & tan-sigmoid transfer function, the MSE is relatively low when the number of hidden layers is between four and fifteen, and the lowest MSE (0.0587 ± 0.0351) is obtained when it reaches seven. Besides, seven hidden layers with the log-sigmoid & purelin transfer function bring the highest determination coefficient (square of R, 0.3819 ± 0.2781). Therefore, we decided to select 10–7–1 as the topological structure of BPNN and the tan-sigmoid transfer function as the transfer function between both input layer & hidden layer and hidden & output layer. To improve the accuracy of the final model, we set the aim as R > 0.9 and trained 10–7–1 BPNN with tan-sigmoid transfer function & tan-sigmoid transfer function several times. The MSE and correlation coefficients of this model were in acceptable ranges as shown in Table 4. These correlation coefficients indicate that there is a strong correlation between the predicted and experimental result and the MSE (0.162) is acceptable. The plots of experimental versus predicted values (Fig. 3) confirmed the discussed results. Hence, we chose this one as the final model. To our knowledge, it was the first time applying BPNN to predict the IC₅₀ of ACE of pentapeptide. In our previous work [58], a BPNN model was built to predict the IC₅₀ of ACE of tripeptide. The MSE of the tripeptide model (0.2148) is higher than this pentapeptide model (0.0162). Meanwhile, the R of the tripeptide model (0.854) was less than this new model (0.9176). These differences possibly resulted from the difference of descriptors.

[Figure omitted. See PDF.]

Fig 1. The MSE of different structure of BPNN models.

https://doi.org/10.1371/journal.pone.0119598.g001

[Figure omitted. See PDF.]

Fig 2. The coefficient of determination (square of pearson correlation coefficient) of different structure of BPNN models.

https://doi.org/10.1371/journal.pone.0119598.g002

[Figure omitted. See PDF.]

Fig 3. Predicted versus experimental plots for BPNN model.

https://doi.org/10.1371/journal.pone.0119598.g003

[Figure omitted. See PDF.]

Table 4. MSE and correlation coefficients of final BPNN model.

https://doi.org/10.1371/journal.pone.0119598.t004

Activity Prediction and Experimental Verification

After simulated enzyme digestion of hemoglobin, seven pentapeptides among all the cleavage fragments were obtained. The descriptor calculated value and the predicted IC₅₀ were shown in Table 3. The pentapeptide WTQRF and AAWGK contain two or three hydrophobic amino acid residues, which is in agreement with the common features of the 24 reported pentapeptides, and show a low predicted IC₅₀. The AHRYH, WTQRF and FTPVL contain hydrophpbic or aromatic amino acid residues in the C- terminus which contribute to ACE inhibition. For the second amino acid residues from C-terminus, the PTTKT, WTQRF and AAWGK contain hydrophilic amino acid residues, which also match the common features of the 24 reported pentapeptides. Obviously, WTQRF conforms to all the aforementioned common features. It include three hydrophilic amino acid residues (Thr, Gln, Arg) and two hydrophobic amino acid residues with aromatic (Trp) / heterocyclic (phe) ring in their side chain. It molecular weight is 736.83 g / mol. The hydrophobic amino acid residues in both of the N- and C- terminus as well as the hydrophilic amino acid residues in the second amino acid residues from C-terminus may make main contribution to ACE inhibition. The result of predicted IC₅₀ indicated that the peptide WTQRF has the best inhibitory activity of ACE. The experimental IC₅₀ of WTQRF is 23.93 μM, which is similar to a pentapeptide identified from egg white protein in our lab before (RVPSL, IC₅₀ 20 μM). [54] The difference between predicted and experimental IC50 may bring from data normalization.

Molecular Docking

Interactions between ACE and inhibitor are fundamental to ACE-inhibitory processes. In this work, docking is used to predict the preferred orientation of WTQRF to ACE when they bound to each other to form a stable complex. The free energy of binding represents the binding affinity between ACE and WTQRF. The molecular docking study of WTQRF presenting within the ACE catalytic site in the presence of the cofactor showed the best pose with potential binding energy:- 749.747 kj / mol. It is indicated that the ligand bind tightly to the receptor and the ACE-WTQRF complex is stable. As we know, ligand-protein affinities are influenced by non-covalent intermolecular interactions between the two molecules. The best pose of WTQRF was stabilized mainly by hydrophobic and hydrophilic interactions (Fig. 4), electrostatic interactions (Electrostatic energy:- 154.728 kj / mol; Van der Waals energy:- 898.743 kj / mol) and hydrogen bond (Table 5).

[Figure omitted. See PDF.]

Fig 4. WTQRF binding with the active site of ACE, the conformation extracted from docking result.

https://doi.org/10.1371/journal.pone.0119598.g004

[Figure omitted. See PDF.]

Table 5. Hydrogen bonds observed between ACE and the best docking pose of WTQRF.

https://doi.org/10.1371/journal.pone.0119598.t005

The hydroxyl group of Tyr 360, the guanidine and the ɛ-amino group of Arg 522 showed significant importance of the binding between WTQRF and ACE by hydrogen bond. The docking result suggested that WTQRF contact with residues Met 223, Ala 354, 356, Tyr 360, 394, Phe 391, Pro 407, 519, Val 518 by hydrophobic interactions and with residues His 353, 383, 387, 413, Ser 355, 517, Arg 402, 522, Glu 384, 403, 411, Asn 406, Val 518 by hydrophilic interactions. The active site of ACE was constituted of a zinc ion and a HEXXH…E motif, including His 383, Glu 384, His 387 on helix α 13 and Glu 411 on helix α 14. [34, 59] The WTQRF were positioned to interact with the HEXXH…E motif. Furthermore, the interactions between inhibitors and Zn²⁺ at the ACE active site also play a significant role in modulating catalysis. [60] It is believed that the shorter distance between the Zn²⁺ and the carbonyl oxygen of the peptide, the greater the degree of ACE inhibition. [34]

Conclusion

In this present work, a back propagation neural network model was built to predict the IC₅₀ of ACE of pentapeptides. The topological structure was 10–7–1 and the transfer function was tan-sigmoid transfer function & tan-sigmoid transfer function. We also built a database of potential pentapeptides derived from bovine hemoglobin by simulated enzyme digestion and found that WTQRF has the highest predictive value with experimental IC₅₀ 23.93 μM. The molecular docking result indicated that Tyr 360 and Arg 522 gave a significantly contribution to the stabilization between WTQRF and ACE. The result also demonstrated that the short distance between the Zn²⁺ and the carbonyl oxygen of the peptide are desirable to the ACE-inhibitory activity. The enzyme kinetics investigation as well as docking and molecular dynamics simulation analysis will be our further study.

Supporting Information

[Figure omitted. See PDF.]

S1 Table. The MSE of different structure of BPNN model.

https://doi.org/10.1371/journal.pone.0119598.s001

(DOC)

S2 Table. The coefficient of determination (square of pearson correlation coefficient) of different structure of BPNN models.

https://doi.org/10.1371/journal.pone.0119598.s002

(DOC)

Author Contributions

Conceived and designed the experiments: TZ BL YY YZ JL. Performed the experiments: TZ BL YY. Analyzed the data: TZ BL. Contributed reagents/materials/analysis tools: SN YZ JL. Wrote the paper: TZ.

Citation: Zhang T, Nie S, Liu B, Yu Y, Zhang Y, Liu J (2015) Activity Prediction and Molecular Mechanism of Bovine Blood Derived Angiotensin I-Converting Enzyme Inhibitory Peptides. PLoS ONE 10(3): e0119598. https://doi.org/10.1371/journal.pone.0119598

References

1. Hernandez-Ledesma B, del Mar Contreras M, Recio I. Antihypertensive peptides: production, bioavailability and incorporation into foods. Adv Colloid Interface Sci. 2011;165(1):23–35. pmid:21185549

2. Lu J, Sawano Y, Miyakawa T, Xue Y-L, Cai M-Y, Egashira Y, et al. One-Week Antihypertensive Effect of Ile-Gln-Pro in Spontaneously Hypertensive Rats. J Agric Food Chem. 2011;59(2):559–563. pmid:21182294

3. Vermeirssen V, Camp JV, Verstraete W. Bioavailability of angiotensin I converting enzyme inhibitory peptides. Br J Nutr. 2004;92:357–366. pmid:15469639

4. Sturrock ED, Natesh R, van Rooyen JM, Acharya KR. Structure of angiotensin I-converting enzyme. Cell Mol Life Sci. 2004;61(21):2677–2686. pmid:15549168

5. Turner AJ, Hooper NM. The angiotensin—converting enzyme gene family genomics and pharmacology. TRENDS PHARMACOL SCI. 2002;23(4):177–183. pmid:11931993

6. Fleming I. Signaling by the angiotensin-converting enzyme. Circ Res. 2006;98(7):887–896. pmid:16614314

7. Natesh R, Schwager SLU, Evans HR, Sturrock ED, Acharya KR. Structural details on the binding of antihypertensive drugs captopril and enalaprilat to human testicular angiotensin I-converting enzyme. Biochemistry. 2004;43(27):8718–8724. pmid:15236580

8. Lavoie JL, Sigmund CD. Minireview: Overview of the Renin-Angiotensin System—An Endocrine and Paracrine System. Endocrinology. 2003;144(6):2179–2183. pmid:12746271

9. Hanif K, Bid HK, Konwar R. Reinventing the ACE inhibitors: some old and new implications of ACE inhibition. Hypertens Res. 2010;33(1):11–21. pmid:19911001

10. Materson BJ, Preston RA. Angiotensin-converting enzyme inhibitors in hypertension. A dozen years of experience. Arch Intern Med. 1994;154(5):513–523. pmid:8122944

11. Udenigwe CC. Bioinformatics approaches, prospects and challenges of food bioactive peptide research. Trends Food Sci Technol. 2014;(36):137–143.

12. Wu J, Ding X. Characterization of inhibition and stability of soy-protein-derived angiotensin I-converting enzyme inhibitory peptides. FOOD RES INT. 2002;35(4):367–375.

13. Chen JR, Okada T, Muramoto K, Suetsuna K, Yang S-C. Identification of Angiotensin I-Converting enzyme inhibitory peptides derived from the peptic digest of soybean protein. J Food Biochem. 2002;26(6):543–554.

14. Kuba M, Tanaka K, Tawata S, Takeda Y, Yasuda M. Angiotensin I-converting enzyme inhibitory peptides isolated from tofuyo fermented soybean food. Biosci Biotech Bioch. 2003;67(6):1278–1283. pmid:12843654

15. Kodera T, Nio N. Identification of an Angiotensin I-converting Enzyme Inhibitory Peptides from Protein Hydrolysates by a Soybean Protease and the Antihypertensive Effects of Hydrolysates in 4 Spontaneously Hypertensive Model Rats. J Food Sci. 2006;71(3):C164–C173.

16. Mallikarjun Gouda K, Gowda LR, Rao AA, Prakash V. Angiotensin I-converting enzyme inhibitory peptide derived from glycinin, the 11S globulin of soybean (Glycine max). J Agric Food Chem. 2006;54(13):4568–4573. pmid:16786999

17. Kuba M, Tana C, Tawata S, Yasuda M. Production of angiotensin I-converting enzyme inhibitory peptides from soybean protein with Monascus purpureus acid proteinase. PROCESS BIOCHEM. 2005;40(6):2191–2196.

18. Matsui T, Li CH, Osajima Y. Preparation and characterization of novel bioactive peptides responsible for angiotensin I-converting enzyme inhibition from wheat germ. J Pept Sci. 1999;5(7):289–297. pmid:10442764

19. Motoi H, Kodama T. Isolation and characterization of angiotensin I-converting enzyme inhibitory peptides from wheat gliadin hydrolysate. Nahrung. 2003;47(5):354–358. pmid:14609094

20. Escudero E, Sentandreu MA, Arihara K, Toldra F. Angiotensin I-converting enzyme inhibitory peptides generated from in vitro gastrointestinal digestion of pork meat. J Agric Food Chem. 2010;58(5):2895–2901. pmid:20151679

21. Jang A, Lee M. Purification and identification of angiotensin converting enzyme inhibitory peptides from beef hydrolysates. Meat Sci. 2005;69(4):653–661. pmid:22063143

22. Gu Y, Majumder K, Wu J. QSAR-aided in silico approach in evaluation of food proteins as precursors of ACE inhibitory peptides. FOOD RES INT. 2011;44(8):2465–2474.

23. Wu J, Aluko RE. Quantitative structure-activity relationship study of bitter di- and tri-peptides including relationship with angiotensin I-converting enzyme inhibitory activity. J Pept Sci. 2007;13(1):63–69. pmid:17117396

24. Wu J, Aluko RE, Nakai S. Structural Requirements of Angiotensin I-Converting Enzyme Inhibitory Peptides: Quantitative Structure-Activity Relationship Modeling of Peptides Containing 4–10 Amino Acid Residues. QSAR Comb Sci. 2006;25(10):873–880.

25. Majumder K, Wu J. A new approach for identification of novel antihypertensive peptides from egg proteins by QSAR and bioinformatics. FOOD RES INT. 2010;43(5):1371–1378.

26. Sagardia I, Roa-Ureta H R, Bald C. A new QSAR model, for angiotensin I-converting enzyme inhibitory oligopeptides. Food Chem. 2013;130:1370–1376.

27. Zhong M, Chong Y, Nie X, Yan A, Yuan Q. Prediction of sweetness by multilinear regression analysis and support vector machine. J Food Sci. 2013;78(9):S1445–1450. pmid:23915005

28. Peng W, Zhong J, Yang J, Ren Y, Xu T, Xiao S, et al. The artificial neural network approach based on uniform design to optimize the fed-batch fermentation condition: application to the production of iturin A. Microb Cell Fact. 2014;13(1):54. pmid:24725635

29. Nedjar-Arroume N, Dubois-Delval V, Adje EY, Traisnel J, Krier F, Mary P, et al. Bovine hemoglobin: An attractive source of antibacterial peptides. Peptides. 2008;29:969–977. pmid:18342399

30. Nedjar-Arroume N, Dubois-Delval V, Miloudi K, Daoud R, Krier F, Kouach M, et al. Isolation and characterization of four antibacterial peptides from bovine hemoglobin. Peptides. 2006;27:2082–2089. pmid:16730859

31. Catiau L, Traisnel J, Delval-DuboiS V, Chihib N- E, Guillochon D, Nedjar-Arroume N. Minimal antimicrobial peptidic sequence from hemoglobin alpha-chain: KYR. Peptides. 2011;32:633–638. pmid:21262306

32. Catiau L, Traisnel J, Chihib N-E, Flem GL, Blanpain A, Melnyk O, et al. RYH A minimal peptidic sequence obtained from beta-chain hemoglobin exhibiting an antimicrobial activity. Peptides. 2011;32:1463–1468. pmid:21645568

33. Soubrier F, Alhenc-Gelas F, Hubert C, Allegrini J, John M, Tregear G, et al. Two putative active centers in human angiotensin I-converting enzyme revealed by molecular cloning. Proc Natl Acad Sci U S A. 1988;85(24):9386–9390. pmid:2849100

34. He R, Aluko RE, Ju X-R. Evaluating Molecular Mechanism of Hypotensive Peptides Interactions with Renin and Angiotensin Converting Enzyme. PloS one. 2014;9(3):e91051. pmid:24603692

35. Watermeyer JM. Structural determinants of the domain-selectivity of novel inhibitors of human testis angiotensin-converting enzyme: University of Cape Town; 2008.

36. Donkor ON, Henriksson A, Singh TK, Vasiljevic T, Shah NP. ACE-inhibitory activity of probiotic yoghurt. Int Dairy J. 2007;17(11):1321–1331.

37. Guang C, Phillips RD. Plant food-derived Angiotensin I converting enzyme inhibitory peptides. J Agric Food Chem. 2009;57(12):5113–5120. pmid:19449887

38. Gómez-Ruiz JÁ, Taborda G, Amigo L, Recio I, Ramos M. Identification of ACE-inhibitory peptides in different Spanish cheeses by tandem mass spectrometry. Eur Food Res Technol. 2006;223(5):595–601.

39. Majumder K, Wu J. Angiotensin I converting enzyme inhibitory peptides from simulated in vitro gastrointestinal digestion of cooked eggs. J Agric Food Chem. 2009;57(2):471–477. pmid:19154160

40. Hernández-Ledesma B, Quirós A, Amigo L, Recio I. Identification of bioactive peptides after digestion of human milk and infant formula with pepsin and pancreatin. Int Dairy J. 2007;17(1):42–49.

41. Hasan F, Kumada Y, Hashimoto N, Katsuda T, Terashima M, Katoh S. Fragmentation of angiotensin-I converting enzyme inhibitory peptides from bonito meat under intestinal digestion conditions and their characterization. Food Bioprod Process. 2006;84(2):135–138.

42. Gómez-Ruiz JÁ, Ramos M, Recio I. Identification of novel angiotensin-converting enzyme-inhibitory peptides from ovine milk proteins by CE-MS and chromatographic techniques. Electrophoresis. 2007;28(22):4202–4211. pmid:17948260

43. Hernandez-Ledesma B, Ramos M, Recio I, Amigo L. Effect of beta-lactoglobulin hydrolysis with thermolysin under denaturing temperatures on the release of bioactive peptides. J Chromatogr A. 2006;1116(1–2):31–37. pmid:16643928

44. Wu H, He H-L, Chen X-L, Sun C-Y, Zhang Y-Z, Zhou B-C. Purification and identification of novel angiotensin-I-converting enzyme inhibitory peptides from shark meat hydrolysate. PROCESS BIOCHEM. 2008;43(4):457–461.

45. Hyoung Lee D, Ho Kim J, Sik Park J, Jun Choi Y, Soo Lee J. Isolation and characterization of a novel angiotensin I-converting enzyme inhibitory peptide derived from the edible mushroom Tricholoma giganteum. Peptides. 2004;25(4):621–627. pmid:15165718

46. Contreras MdM, Carrón R, Montero MJ, Ramos M, Recio I. Novel casein-derived peptides with antihypertensive activity. Int Dairy J. 2009;19(10):566–573.

47. Muguruma M, Ahhmed AM, Katayama K, Kawahara S, Maruyama M, Nakamura T. Identification of pro-drug type ACE inhibitory peptide sourced from porcine myosin B: Evaluation of its antihypertensive effects in vivo. Food Chem. 2009;114(2):516–522.

48. Ren Y, Wan D-G, Lu X-m, Chen L, Zhang T-e, Guo J-L. Isolation and characterization of angiotensin I-converting enzyme inhibitor peptides derived from porcine hemoglobin. Sci Res Essays. 2011;6(30):6262–6269.

49. Lipinski CA, Lombardo F, Dominy BW, Feeney PJ. Experimental and computational approaches to estimate solubility and permeability in drug discovery and development settings. Adv Drug Deliv Rev. 2012;64:4–17.

50. Shen Q, Jiang J-H, Jiao C-X, Shen G-l, Yu R-Q. Modified particle swarm optimization algorithm for variable selection in MLR and PLS modeling: QSAR studies of antagonism of angiotensin II antagonists. EUR J PHARM SCI. 2004;22(2):145–152.

51. Ghose AK, Crippen GM. Atomic physicochemical parameters for three-dimensional-structure-directed quantitative structure-activity relationships. 2. Modeling dispersive and hydrophobic interactions. J Chem Inf Comput Sci. 1987;27(1):21–35. pmid:3558506

52. Ertl P, Rohde B, Selzer P. Fast calculation of molecular polar surface area as a sum of fragment-based contributions and its application to the prediction of drug transport properties. J Med Chem. 2000;43(20):3714–3717. pmid:11020286

53. Jalali-Heravi M, Asadollahi-Baboli M, Shahbazikhah P. QSAR study of heparanase inhibitors activity using artificial neural networks and Levenberg—Marquardt algorithm. Eur J Med Chem. 2008;43(3):548–556. pmid:17602800

54. Liu J, Yu Z, Zhao W, Lin S, Wang E, Zhang Y, et al. Isolation and identification of angiotensin-converting enzyme inhibitory peptides from egg white protein hydrolysates. Food Chem. 2010;122(4):1159–1163.

55. Yu Z, Liu B, Zhao W, Yin Y, Liu J, Chen F. Primary and secondary structure of novel ACE-inhibitory peptides from egg white protein. Food Chem. 2012;133(2):315–322.

56. Yu Z, Zhao W, Liu J, Lu J, Chen F. QIGLF, a novel angiotensin I-converting enzyme-inhibitory peptide from egg white protein. J Sci Food Agric. 2011;91(5):921–926. pmid:21384361

57. Yu Z, Yin Y, Zhao W, Yu Y, Liu B, Liu J, et al. Novel peptides derived from egg white protein inhibiting alpha-glucosidase. Food Chem. 2011;129(4):1376–1382. pmid:21214175

58. Liu B. The quantitative structure activity relationships of angiotensin converting enzyme inhibitor peptide from egg white. Changchun: Jilin University; 2013.

59. Ni H, Li L, Liu G, Hu S- Q. Inhibition mechanism and model of an angiotensin I-converting enzyme (ACE)-inhibitory hexapeptide from yeast (Saccharomyces cerevisiae). PloS one. 2012;7(5):e37077. pmid:22606330

60. Natesh R, Schwager SLU, Sturrock ED, Acharya KR. Crystal structure of the human angiotensin-converting enzyme-lisinopril complex. nature. 2003;421(30):551–554.