Quantitative Structure–Activity Relationship Modeling Coupled with Molecular Docking Analysis in Screening of Angiotensin I-Converting Enzyme Inhibitory Peptides from Qula Casein Hydrolysates Obtained by Two-Enzyme Combination Hydrolysis

In this study, Qula casein derived from yak milk casein was hydrolyzed using a two-enzyme combination approach, and high angiotensin I-converting enzyme (ACE) inhibitory activity peptides were screened by quantitative structure-activity relationship (QSAR) modeling integrated with molecular docking analysis. Hydrolysates (<3 kDa) derived from combinations of thermolysin + alcalase and thermolysin + proteinase K demonstrated high ACE inhibitory activities. Peptide sequences in hydrolysates derived from these two combinations were identified by liquid chromatography-tandem mass spectrometry (LC-MS/MS). On the basis of the QSAR modeling prediction, a total of 16 peptides were selected for molecular docking analysis. The docking study revealed that four of the peptides (KFPQY, MPFPKYP, MFPPQ, and QWQVL) bound the active site of ACE. These four novel peptides were chemically synthesized, and their IC50 was determined. Among these peptides, KFPQY showed the highest ACE inhibitory activity (IC50 = 12.37 ± 0.43 μM). Our study indicated that Qula casein presents an excellent source to produce ACE inhibitory peptides.