Structural and molecular bases of angiotensin-converting enzyme inhibition by bovine casein-derived peptides: an in silico molecular dynamics approach

The angiotensin-converting enzyme (ACE) plays a key role in blood pressure regulation process, and its inhibition is one the main drug targets for the treatment of hypertension. Though various peptides from milk proteins are well-known for their ACE-inhibitory capacity, research devoted to understand the molecular bases of such property remain scarce, specifically for such peptides. Therefore, in this work, computational molecular docking and molecular dynamics calculations were performed to enlighten the intermolecular interactions involved in ACE inhibition by six different casein-derived peptides (FFVAPFPEVFGK, FALPQYLK, ALNEINQFYQK, YLGYLEQLLR, HQGLPQEVLNENLLR and NAVPITPTLNR). Two top ranked docking poses for each peptide (one with N- and the other C-terminal peptide extremity oriented towards the ACE active site) were selected for dynamic simulations (50 ns; GROMOS53A6 force field), and the results were correlated to in vitro ACE inhibition capacity. Two molecular features appeared to be essential for peptides to present high ACE inhibition capacity in vitro: i) to interact with the S1 active site residues (Ala354, Glu384, and Tyr523) by hydrogen bonds; ii) to interact with Zn2+ coordinated residues (His383, His387, and Glu411) by short-lenght hydrogen bonds, as observed in the cases of ALNEINQFYQK (IACE = 80.7%), NAVPITPTLNR (IACE = 80.7%), and FALPQYLK (IACE = 79.0%). Regardless of the temporal stability of these strong interactions, they promoted some disruption of Zn2+ tetrahedral coordination during the molecular dynamics trajectories, and were pointed as the main reason for the greatest ACE inhibition by these peptides. On the other hand, peptides with intermediate inhibition capacity (50% < I-ACE < 45%) interacted mainly by weaker interactions (e.g.: electrostatic and hydrophobic) with the Zn2+ coordinated residues, and were not able to change significantly its tetrahedral coordination structure. These findings may: i) assist the discrimination in silico of "good" and "bad" ACE-inhibitory peptides from other food sources, and/or ii) to aid in designing de novo new molecules with ACE-inhibitory capacity. Communicated by Ramaswamy Sarma