Virtual Screening and Mechanism of Action of DPP-IV Inhibitory Peptides from Pinctada fucata (P. fucata)
-
摘要: 为了快速发现马氏珍珠贝来源二肽基肽酶IV（dipeptidyl peptidase IV，DPP-IV）抑制活性肽，本研究利用数据库中20个已报道的DPP-IV抑制肽组成训练集，构建了药效团模型并通过测试集分子和Fisher随机验证法对模型进行了评估。利用在线网站PeptideCutter，以珍珠贝肉蛋白为原料，进行虚拟酶解。使用最优药效团模型（Hypo1）对虚拟酶解获得的192个低分子肽（氨基酸数小于或等于5）进行初步筛选，接着以分子对接的方法进一步筛选，并且对潜在的DPP-IV抑制活性肽进行固相合成，体外验证其活性并分析其作用机制。结果表明，药效团结合分子对接技术筛选了4个理论上可能具有高活性的DPP-IV抑制肽，即LPIY、VQDR、PIY和APSL。其中，VQDR、PIY没有表现出抑制活性，而LPIY和APSL具有较高的DPP-IV抑制活性，其IC50值分别为521.19 和258.67 µmol/L。与DPP-IV相互作用的机理表明，肽LPIY和APSL与DPP-IV活性口袋中的氨基酸残基形成了多个氢键，且N-端第二个位置的脯氨酸（P）也均与活性口袋中的残基形成了Pi-Alkyl作用，因此促进了LPIY和APSL的DPP-IV抑制作用。本研究为高效筛选珍珠贝来源的DPP-IV抑制肽提供了理论方法。Abstract: In order to discover the dipeptidyl peptidase IV (DPP-IV) peptides quickly from Pinctada fucata (P. fucata), in this research, 20 reported DPP-IV inhibitory peptides in the database were used to form a training set and construct pharmacophore models, which were then valuated by test set molecules and fisher random validation. The online website PeptideCutter was employed to perform the virtual hydrolysis with P. fucata meat protein as raw materials. The optimal pharmacophore model (Hypo1) was used to initially screen the 192 low-molecular peptides (the number of amino acids was no more than 5), which were obtained by the virtual hydrolysis, and then the molecular docking was performed for further screening. Furthermore, the potential DPP-IV inhibitory peptides were synthesized using a solid-phase method, and their activities were verified in vitro and the interaction mechanism was analyzed. The results showed that 4 potential DPP-IV inhibitory peptides, namely LPIY, VQDR, PIY, and APSL were screened through pharmacophore combined with molecule docking. Among them, VQDR and PIY did not show inhibitory activity, while LPIY and APSL had strong DPP-IV inhibitory activity, with IC50 values of 521.19 and 258.67 µmol/L respectively. The interaction mechanism indicated both peptides formed multiple hydrogen bonds with the amino acid residues within the active pocket of DPP-IV, and the proline (P) at the second position of the N-terminal constructed Pi-Alkyl interactions with the residues of the active pocket, which promoted the DPP-IV inhibitory activity for LPIY and APSL. The current work would provide a theoretical method for screening DPP-IV inhibitory peptides derived from P. fucata efficiently.
-
近年来，活性肽作为一种天然安全的食品补充剂得到了迅速发展，其重要的生理调节功能致使其在人类健康保健预防上占有一席之地[1]。目前的研究已证实牛奶，大豆等食品蛋白的水解物具有一系列生理活性，如抗氧化、抗炎、降血压、降血糖血脂等[2-6]。糖尿病是一种高发高危性的慢性疾病，据报道，到2030年，全球糖尿病患病人数将突破4亿，其中2型糖尿病占绝大多数，约占比90%[7]。在糖尿病的发生发展中，胰岛β细胞分泌的胰岛素多肽（GIP）和胰高血糖素样肽1（GLP-1）起着关键的调控作用，而二肽基肽酶IV（DPP-IV）则可以降解这些激素，致使其在体内的含量降低，从而导致血糖的升高[8-10]。因此，抑制DPP-IV的活性，成为控制血糖的关键步骤。已知上市的具有DPP-IV 抑制活性的药物有西格列汀、沙格列汀、阿格列汀等[11-12]。然而，其副作用使这些化学药物的使用受到限制。因此，寻找安全有效的DPP-IV活性抑制替代物成为近年来的研究热点。经研究发现，某些具有2~8个氨基酸的肽可抑制DPP-IV，其活性受氨基酸组成，序列长短，疏水性等各种因素影响[13]。
在活性肽的制备过程中，传统的制备方法是以活性为导向的一系列分离纯化，然后鉴定出活性肽单体，其过程繁琐，需要耗费大量的人力物资，且收率低[14]，并且因其分段的活性测定方法导致错过某些活性极强的单肽。因此需要寻找一种高效寻找活性肽的方法，以实现大量制备。
随着计算机技术的发展，其在活性药物筛选，结构优化设计上近年来有了明显的突破，使用计算机对小分子化合物进行虚拟筛选，结构设计在医药领域已得到广泛的应用。然而，使用这一技术对活性肽进行全面系统的虚拟筛选并不多见。常见的虚拟筛选方法包括基于受体蛋白的分子对接以及基于配体小分子的药效团模型。Nongonierma等[15]运用分子对接技术预测了二肽抑制DPP-IV的潜力，然而，利用药效团模型对DPP-IV抑制活性多肽进行筛选并不多见。Wang等[16]采用已知的ACE抑制肽构建了药效团模型，并进行了初步验证，然而并未利用构建的药效团进行肽分子的筛选。
在前期的研究中，利用Biopep数据库采用在线虚拟酶解方法，发现马氏珍珠贝蛋白来源的活性肽具有较高的DPP-IV抑制潜力。因此本研究继续以马氏珍珠贝蛋白为研究对象，利用PeptideCutter平台在线酶解，产生小分子肽库，通过建立DPP-IV抑制肽的药效团模型，采用最优药效团模型及分子对接方法对肽库进行虚拟筛选，以鉴定出具有高活性的DPP-IV抑制肽。
1. 材料与方法
1.1 材料与仪器
二肽基肽酶IV（DPP-IV）、抑二肽素A（IPI）、甘氨酰-脯氨酰-对硝基苯胺对甲苯磺酸（Gly-Pro-pNA） 美国Sigma公司；合成肽LPIY、VQDR、PIY和APSL（纯度>98%） 吉尔生化（上海）有限公司。
Discovery Studio（DS）软件 美国Accelrys公司；K3 Plus酶标仪 上海宝予德科学仪器有限公司。
1.2 实验方法
1.2.1 药效团模型的产生
如图1所示，从Biopep数据库中，收集了40个已知的具有DPP-IV抑制活性小肽作为模型建立和验证的数据集[17]，其中20个作为训练集用于产生药效团模型，剩余的20个组成测试集用作药效团模型的验证。采用ChemOffice2018构建活性肽结构，并赋与CHARMM力场进行能量优化。采用Discovery Studio软件中“3D QSAR Pharmacophore Generation”生成药效团模型。在构象产生时，设置“Best”构象生成方法，并限定产生的最大构象数为250，新产生的构象与最低能量构象之间的能量差不大于20 kcal·mol−1，用于产生药效团模型的化学特征有：氢键受体（A）、氢键给体（D）、疏水中心（H）和阴离子中心（N）。将不确定值（Uncertainty）值设为2，其余参数为软件默认值。
![]()
图 1 用于构建药效团模型的训练集的肽结构Figure 1. Structure of peptides in the training set used to construct the pharmacophore model1.2.2 药效团模型的验证
为检验药效团是否具有预测能力，采用测试集分子，Fisher随机验证对药效团模型进行评价。
1.2.2.1 测试集分子验证
为检测该药效团模型的预测能力，本研究选取除训练集以外已报道的20个不同结构的DPP-IV抑制活性肽分子作为测试集，采用DS软件中的“Ligand Pharmacophore Mapping”模块对模型进行验证，测试集肽分子的构象产生方式与训练集相同。
1.2.2.2 Fisher 随机验证
Fisher随机验证是通过对训练集中的分子重新分配活性值，从而产生药效团模型，如果随机产生的药效团模型优于之前产生的模型，则证明先前产生的模型是随机产生的，不具有统计学意义，反之则具有统计学意义，通过设置置信度为95%，对模型进行Fisher随机验证。
1.2.3 肽库的构建
从UniProt数据库（http://www.uniprot.org/）中选取了十一条来源于马氏珍珠贝肉的蛋白，其编码为M5AKJ3，M5AJM7，M5AKK3，M5AJP7，B3GRV7，Q2AAD2，Q967M3，Q0ZVW7，A0A510EA34，A0A0N7J6H6，A0A0P0LDM5。参考Yang等[18]的做法，采用PeptideCutter在线酶切所选的蛋白序列，所用到的蛋白酶有胰蛋白酶、胃蛋白酶、糜蛋白酶以及脯氨酸内切蛋白酶，在线虚拟酶解完后，选取氨基酸个数小于或等于5的小分子多肽组成多肽库。
1.2.4 药效团模型对肽库的筛选
使用构建好的药效团模型对建立的小分子多肽库进行搜索，整个搜索过程在“Search 3D Database”中完成，设定搜索方法为‘BEST’，其余参数保持不变。
1.2.5 分子对接
参考Nongonierma等[19]的方法进行分子对接，DPP-IV 3D结构从PDB生物大分子数据库（PDB code: 1WCY）检索获得，使用AutoDock Tools （version 1.5.6）对其去水，加氢并设置搜索空间。格点坐标及盒子体积设置如下：size_x=18，size_y=16，size_z=18，center_x=56.133，center_y=63.602，center_z=39.577。对接前将蛋白晶体结构中原始配体Diprotin A提取出来重新对接，计算RMSD值，如果RMSD<2 Å，证明对接程序可信。将经药效团模型筛选后理论活性值较高的前23个多肽（匹配值FitValue>6.0）与DPP-IV进行Autodock Vina分子对接，进一步筛选得到理论上活性高的多肽。肽结构均由ChemOffice 2018软件画的，能量最小化后与蛋白对接，产生的结合能用于表示小分子与蛋白结合的稳定性，结合能越低，表明小分子与蛋白结合得越稳定，其理论上的活性值越大。
1.2.6 肽的固相合成
肽由上海吉尔生化有限公司合成的，纯度大于98%，其分子量采用LC-MS确定，纯度采用HPLC测定[20]。
1.2.7 体外DPP-IV抑制活性测定
DPP-IV抑制实验参考Ibrahim等[21]的方法稍作修改。将25 µL的DPP-IV（总浓度为0.0035 U/mL）与100 µL的样品溶液混合并保温10 min，控制组以Tris-HCl缓冲液代替样品加入。随后加50 µL浓度为1 mM底物（Gly-Pro-pNA溶液）分别到样品组和控制组，继续保温反应1 h，使用酶标仪测量405 nm的吸光度。每个样品进行三次重复实验，DPP-IV酶活性抑制率根据以下公式计算：
DPP-IV抑制率（%）=（1−AS/AC）×100
式中：As组为多肽样品存在时的吸光度；Ac为控制组的吸光度。
1.3 数据处理
采用GraphPad Prism 8.0.1的单因素方差分析进行统计分析。所有实验均重复三次，数据表示为平均值±标准偏差，以P<0.05时表示差异显著。
2. 结果与分析
2.1 药效团模型的产生
由20个多肽组成的训练集产生了10个药效团模型，其相应的统计评分如表1所示，一般来说，评差（Δcost），config值，均差（RMS），以及相关度可用于综合评价药效团是否具有统计学意义。在理想的药效团模型中，总评值（total cost）应位于固定值（fixed cost）和空值（null cost）之间，接近前者，并远离后者，并且，null cost和total cost差值（Δcost）应大于60，才表明该模型统计学意义大于90%，即说明模型具有良好的预测能力[22-24]。另外，模型的config值应小于17，RMS值应小于1.5，相关度应大于0.9。从下表1可以看出，产生的十个药效团的Δcost，RMS，相关度都在合理的范围内，均具有统计学意义。同时，Hypo1因具有最大的Δcost和相关度，最小的RMS，被认为是理论上最优药效团模型。Hypo1对训练集的预测结果如表2所示，从表2可以看出，其对训练集中各分子的预测IC50值与真实IC50值较为接近，误差值（Error）保持在±4以内，预测结果良好。
表 1 建立的药效团的各项参数Table 1. Parameters of the established pharmacophore药效团（Hypo） 总评 评差 均差 相关度 特征 1 86.6076 67.036 1.029 0.943 HBA，HBA，HBD，HY 2 87.6027 66.041 1.061 0.939 HBD，HBD，HY，HY，Ne 3 88.5037 65.140 1.116 0.933 HBA，HBD，HY，Ne 4 88.8108 64.833 1.124 0.932 HBD，HY，HY，Ne 5 89.1022 64.542 1.114 0.933 HBA，HBA，HBD，HY 6 89.1124 64.532 1.124 0.932 HBD，HBD，HY，HY 7 89.4906 64.153 1.130 0.931 HBA，HBD，HY，HY 8 90.3185 63.326 1.184 0.924 HBA，HY，HY，Ne 9 90.6879 62.956 1.201 0.922 HBA，HY， HY，Ne 10 91.4894 62.155 1.236 0.917 HBA，HBD，HBD，HY 注：Null cost=153.644; Fixed cost=76.0183; Config=16.8324; △cost=null cost−total cost。 表 2 Hypo1对训练集分子的预测Table 2. Prediction of Hypo1 for training set molecules多肽序列 真实IC50 预测IC50 误差 不确定值 Ile-Pro-Ile 3.5 5.7 1.6 2 Val-Pro-Leu 16 19 1.2 2 Ile-Pro-Ala 49 71 1.5 2 Leu-Leu-Ala-Pro 54 94 1.7 2 Ile-Pro-Met 69 21 −3.3 2 Thr-Trp 84 130 1.5 2 Ile-Pro 150 500 3.3 2 Val-Ala 170 440 2.6 2 Pro-Pro-Leu 390 640 1.6 2 Gly-Pro-Met 420 250 −1.7 2 Tyr-Pro 660 630 −1.1 2 Ala-Leu 880 1800 2 2 Leu-Trp 990 420 −2.4 2 Leu-Gln-Pro 1200 420 −2.8 2 Pro-Pro-Gly 2300 950 −2.4 2 Gly-Leu 2600 1700 −1.5 2 Ser-Pro 6000 1700 −3.4 2 Lys-Ala 6300 6500 1 2 Ala-Pro 8000 10000 1.3 2 Gly-Pro 9700 16000 1.6 2 Hypo1含有的药效团特征有氢键受体，氢键给体，疏水特征，各药效团之间的空间距离如下图2A所示。Hypo1对训练集中的活性最强，即真实IC50值最小的肽分子（Ile-Pro-Ile）以及活性最弱的肽分子（Gly-Pro）预测结果分别如图2B和2C所示，活性最强的分子同药效团完全匹配上，而活性最弱的分子则有一个疏水特征没有匹配上。
![]()
图 2 最优药效团模型及其与训练集分子的匹配注：A为DPP-IV抑制肽产生的最优药效团模型Hypo1及其相互间距离；B为Hypo1对训练集中活性最高肽分子的匹配；C为Hypo1对训练集中活性最低肽分子的匹配；浅蓝色代表疏水中心，紫色代表氢键给体，绿色代表氢键受体。Figure 2. Optimal pharmacophore model and its aligning with training set molecules2.2 药效团模型的验证
为了验证Hypo1对除训练集以外的其他肽分子也有预测作用，以文献中20个具有DPP-IV抑制作用的其他分子组成测试集，对Hypo1进行验证，结果如图3A所示，得到其真实活性与预测活性值之间的决定系数R2为0.81，相关性良好。这表明Hypo1不仅能预测训练集中的肽，对除训练集以外的多肽也具有一定活性预测能力。采用Fisher随机验证以判断药效团的产生是否具有偶然性，结果如图4B，可以看出，随机产生的19个药效团模型相比Hypo1（original），其总评值（cost值）都要更大，这表明本研究中产生的Hypo1不是偶然产生的，具有一定的统计学意义。
![]()
图 3 药效团模型的验证注：A为测试集对最优药效团Hypo1的预测；B为置信度为95%的Fisher随机验证。Figure 3. Verification of the pharmacophore model2.3 肽库的筛选及活性验证
将珍珠贝肉虚拟酶解产生的192个小分子多肽组成多肽库，能量最小化后，使用Hypo1对多肽库进行筛选，选择Fit value值大于6.0的肽进行分子对接，以进一步筛选高活性DPP-IV抑制肽，其结果如下表3所示。从表3可以看出，多肽LPIY，VQDR，PIY及APSL与DPP-IV的结合能最低，分别为−8.0，−7.8，−7.6，−7.6 kcal/mol。因此选择这四个肽分别固相合成，测定2.5 mg/mL时它们的DPP-IV抑制率，结果如图4所示，VQDR和PIY并没有展示活性，而LPIY及APSL却表现了极好的DPP-IV抑制率，分别为79.43%±0.61%和80.3%±0.24%（P>0.05），而阳性对照IPI在浓度为1.5 mg/mL时，DPP-IV抑制率为80.06%±0.44%。进一步测得LPIY和APSL的IC50值，分别为521.19和258.67 µmol/L。通过观察这两个肽分子的结构特点，可以发现它们的N-端第二位均为脯氨酸（Pro），且两端均有疏水性氨基酸。Nongonierma等[25]研究报道，肽氨基酸序列中含有脯氨酸可以极大地促进DPP-IV抑制活性，Nong等[26]也指出，N-端第二个位置为Pro或Ala是DPP-IV抑制肽的明显特征。并且，两端的疏水性氨基酸也可以促进酶抑制活性的提高，本研究与以往的报道相符。
表 3 23个预测的DPP-IV抑制肽及分子对接结合能Table 3. 23 predicted DPP-IV inhibitory peptides and molecular docking binding energy编号 肽序列 匹配值 结合能（kcal/mol） 1 TVR 6.81 −7.4 2 VQDR 6.73 −7.8 3 TVPIY 6.66 −6.5 4 NAPAM 6.60 −6.9 5 IDL 6.54 −7.0 6 IE 6.49 −6.6 7 IQL 6.48 −7.2 8 PIY 6.45 −7.6 9 NSIM 6.41 −7.0 10 LPIY 6.34 −8.0 11 ETT 6.34 −7.0 12 ISGAE 6.30 −7.1 13 AAS 6.28 −6.3 14 AVF 6.21 −7.3 15 IGGSI 6.21 −7.4 16 APSL 6.20 −7.6 17 IEDM 6.19 −7.3 18 VAVD 6.17 −7.0 19 VTY 6.12 −7.2 20 ISK 6.08 −6.1 21 NPPK 6.08 −7.3 22 SEAE 6.03 −7.0 23 GSL 6.01 −6.2 2.4 分子作用机理
在多肽分子对接前，首先将原配体重新对接后到晶体结构中，其RMSD值为1.2Å，表明对接算法程序适用于本次对接。图5展示了两个高活性的四肽LPIY及APSL与DPP-IV之间的相互作用图，从图5A可以看出，两个肽均很好地嵌合在DPP-IV的活性口袋中。DPP-IV一共包含有三个活性口袋，即S1、S2及S3，S1口袋包括氨基酸残基Ser630、Asn710和His740，S2口袋由Arg125、Glu205和Glu206组成，S3口袋由Tyr547、Arg358和Phe357组成[27-28]。从图5B可以看出，肽LPIY与DPP-IV中的Asn710形成了两个氢键，与Glu206，Arg669各形成了一个氢键，其中Asn710为S1中的残基，Glu206为S2中的氨基酸残基。其中的脯氨酸（Pro）与Tyr547形成了Pi-Alkyl作用，Tyr547位于S3活性口袋中。图5C表示四肽APSL与氨基酸残基Glu205形成了两个氢键，与Glu206、His740、Asn710、Arg125和Tyr631各形成了一个氢键，其中His740和Asn710属于S1口袋，Glu205、Glu206和Arg125属于S2口袋，中间的脯氨酸与Phe357形成了Pi-Alkyl作用，Phe357也是S3活性口袋中的氨基酸残基。可见，两个肽均与DPP-IV活性口袋中的氨基酸残基形成了相互作用，并且，它们N-端第二位的脯氨酸均与S3口袋中的氨基酸残基形成了Pi-Alkyl作用，这些作用对于肽-酶结合的稳定性起到了极为关键的作用。另外，肽与DPP-IV之间还有范德华作用，Pi-Sigma作用等。Xu等[29]发现他们鉴定的所有的肽均与S1，S2或S3口袋残基有着相互作用，如PAGPF（Asn740和Ser630）、KTMPGP（Tyr547和Arg358）、IPQVS（Asn740，Arg125和Glu206）、ELHQEEPL（Arg125，Glu206和Arg358），且这些肽均含有一个或多个脯氨酸，这说明脯氨酸及其与活性口袋残基之间的相互作用可能是肽具有DPP-IV抑制作用的重要原因之一。
![]()
图 5 活性肽与DPP-IV的对接注：A为LPIY和APSL在DPP-IV活性口袋中的3D图；B为LPIY与DPP-IV对接时的2D图；C为APSL与DPP-IV对接时的2D图。Figure 5. Molecular docking between the peptides and DPP-IV3. 结论
本研究构建了用于DPP-IV抑制肽筛选的3D-QSAR药效团模型，其中Hypo1经评估后为最优药效团模型，其Δcost为67.036，相关度为0.943，RMS为1.029，含有的药效团特征为两个氢键受体，一个氢键给体，一个疏水基团。采用在线网站PeptideCutter对珍珠贝肉蛋白进行虚拟酶解，产生的192个氨基酸数小于或等于5的小分子肽组成肽库，采用Hypo1以及分子对接对肽库进行筛选，得到四个理论上具有高潜力的DPP-IV抑制活性肽（LPIY，VQDR，PIY和APSL）。活性验证后发现VQDR和PIY没有活性，而LPIY和APSL具有较高的DPP-IV抑制活性，其IC50值分别为521.19 和258.67 µmol/L。经文献调研后发现肽序列中含有脯氨酸并且脯氨酸在N-端第二位时其DPP-IV抑制活性明显，本研究筛选的两个高活性肽符合这一特征。这一结果也表明，以后的研究工作中，虚拟筛选时若结合构效关系研究将进一步提高筛选的成功率。两个肽与DPP-IV的结合模式表明脯氨酸与活性口袋中氨基酸残基形成的相互作用对于维持肽-酶稳定性起到了关键作用，下一步的研究将集中在体内活性抑制机理中。
-
![]()
图 1 用于构建药效团模型的训练集的肽结构
Figure 1. Structure of peptides in the training set used to construct the pharmacophore model
![]()
图 2 最优药效团模型及其与训练集分子的匹配
注：A为DPP-IV抑制肽产生的最优药效团模型Hypo1及其相互间距离；B为Hypo1对训练集中活性最高肽分子的匹配；C为Hypo1对训练集中活性最低肽分子的匹配；浅蓝色代表疏水中心，紫色代表氢键给体，绿色代表氢键受体。
Figure 2. Optimal pharmacophore model and its aligning with training set molecules
![]()
图 3 药效团模型的验证
注：A为测试集对最优药效团Hypo1的预测；B为置信度为95%的Fisher随机验证。
Figure 3. Verification of the pharmacophore model
![]()
图 5 活性肽与DPP-IV的对接
注：A为LPIY和APSL在DPP-IV活性口袋中的3D图；B为LPIY与DPP-IV对接时的2D图；C为APSL与DPP-IV对接时的2D图。
Figure 5. Molecular docking between the peptides and DPP-IV
表 1 建立的药效团的各项参数
Table 1 Parameters of the established pharmacophore
药效团（Hypo） 总评 评差 均差 相关度 特征 1 86.6076 67.036 1.029 0.943 HBA，HBA，HBD，HY 2 87.6027 66.041 1.061 0.939 HBD，HBD，HY，HY，Ne 3 88.5037 65.140 1.116 0.933 HBA，HBD，HY，Ne 4 88.8108 64.833 1.124 0.932 HBD，HY，HY，Ne 5 89.1022 64.542 1.114 0.933 HBA，HBA，HBD，HY 6 89.1124 64.532 1.124 0.932 HBD，HBD，HY，HY 7 89.4906 64.153 1.130 0.931 HBA，HBD，HY，HY 8 90.3185 63.326 1.184 0.924 HBA，HY，HY，Ne 9 90.6879 62.956 1.201 0.922 HBA，HY， HY，Ne 10 91.4894 62.155 1.236 0.917 HBA，HBD，HBD，HY 注：Null cost=153.644; Fixed cost=76.0183; Config=16.8324; △cost=null cost−total cost。 
下载: 导出CSV
表 2 Hypo1对训练集分子的预测
Table 2 Prediction of Hypo1 for training set molecules
多肽序列 真实IC50 预测IC50 误差 不确定值 Ile-Pro-Ile 3.5 5.7 1.6 2 Val-Pro-Leu 16 19 1.2 2 Ile-Pro-Ala 49 71 1.5 2 Leu-Leu-Ala-Pro 54 94 1.7 2 Ile-Pro-Met 69 21 −3.3 2 Thr-Trp 84 130 1.5 2 Ile-Pro 150 500 3.3 2 Val-Ala 170 440 2.6 2 Pro-Pro-Leu 390 640 1.6 2 Gly-Pro-Met 420 250 −1.7 2 Tyr-Pro 660 630 −1.1 2 Ala-Leu 880 1800 2 2 Leu-Trp 990 420 −2.4 2 Leu-Gln-Pro 1200 420 −2.8 2 Pro-Pro-Gly 2300 950 −2.4 2 Gly-Leu 2600 1700 −1.5 2 Ser-Pro 6000 1700 −3.4 2 Lys-Ala 6300 6500 1 2 Ala-Pro 8000 10000 1.3 2 Gly-Pro 9700 16000 1.6 2
下载: 导出CSV
表 3 23个预测的DPP-IV抑制肽及分子对接结合能
Table 3 23 predicted DPP-IV inhibitory peptides and molecular docking binding energy
编号 肽序列 匹配值 结合能（kcal/mol） 1 TVR 6.81 −7.4 2 VQDR 6.73 −7.8 3 TVPIY 6.66 −6.5 4 NAPAM 6.60 −6.9 5 IDL 6.54 −7.0 6 IE 6.49 −6.6 7 IQL 6.48 −7.2 8 PIY 6.45 −7.6 9 NSIM 6.41 −7.0 10 LPIY 6.34 −8.0 11 ETT 6.34 −7.0 12 ISGAE 6.30 −7.1 13 AAS 6.28 −6.3 14 AVF 6.21 −7.3 15 IGGSI 6.21 −7.4 16 APSL 6.20 −7.6 17 IEDM 6.19 −7.3 18 VAVD 6.17 −7.0 19 VTY 6.12 −7.2 20 ISK 6.08 −6.1 21 NPPK 6.08 −7.3 22 SEAE 6.03 −7.0 23 GSL 6.01 −6.2
下载: 导出CSV
-
[1] Priya S. Therapeutic perspectives of food bioactive peptides: A mini review[J]. Protein and Peptide Letters,2019,26(9):664−675. doi: 10.2174/0929866526666190617092140
[2] Singh B P, Vij S. In vitro stability of bioactive peptides derived from fermented soy milk against heat treatment, pH and gastrointestinal enzymes[J]. LWT-Food Science and Technology,2018,91:303−307. doi: 10.1016/j.lwt.2018.01.066
[3] Parmar H, Hati S, Sakure A. In vitro and in silico analysis of novel ace-inhibitory bioactive peptides derived from fermented goat milk[J]. International Journal of Peptide Research and Therapeutics,2018,24:441−453. doi: 10.1007/s10989-017-9630-4
[4] Akan E. An evaluation of the in vitro antioxidant and antidiabetic potentials of camel and donkey milk peptides released from casein and whey proteins[J]. Journal of Food Science and Technology,2020:1−9.
[5] Chatterjee C, Gleddie S, Xiao C W, et al. Soybean bioactive peptides and their functional properties[J]. Nutrients,2018,10(9):16.
[6] Lammi C, Arnoldi A, Aiello G. Soybean peptides exert multifunctional bioactivity modulating 3-hydroxy-3-methylglutaryl-CoA reductase and dipeptidyl peptidase-IV targets in vitro[J]. Journal of Agricultural and Food Chemistry,2019,67:4824−4830. doi: 10.1021/acs.jafc.9b01199
[7] Power O, Nongonierma A B, Jakeman P, et al. Food protein hydrolysates as a source of dipeptidyl peptidase IV inhibitory peptides for the management of type 2 diabetes[J]. Proceedings of the Nutrition Society,2014,73(1):34−46. doi: 10.1017/S0029665113003601
[8] Jia C L, Hussain N, Obaroakpo J U, et al. Generation and characterization of dipeptidyl peptidase-IV inhibitory peptides from trypsin-hydrolyzed α-lactalbumin-rich whey proteins[J]. Food Chemistry,2020,318:126333. doi: 10.1016/j.foodchem.2020.126333
[9] Jin R T, Teng X Y, Shang J Q, et al. Identification of novel DPP-IV inhibitory peptides from Atlantic salmon (Salmo salar) skin[J]. Food Research International,2020,133:10.
[10] Nongonierma A B, Fitzgerald R J. Susceptibility of milk protein-derived peptides to dipeptidyl peptidase IV (DPP-IV) hydrolysis[J]. Food Chemistry,2014,145:845−852. doi: 10.1016/j.foodchem.2013.08.097
[11] Rivero-Pino F, Espejo-Carpio F J, Guadix E M. Production and identification of dipeptidyl peptidase IV (DPP-IV) inhibitory peptides from discarded Sardine pilchardus protein[J]. Food Chemistry,2020,328:7.
[12] Liu R, Zhou L, Zhang Y, et al. Rapid identification of dipeptidyl peptidase-IV (DPP-IV) inhibitory peptides from ruditapes philippinarum hydrolysate[J]. Molecules,2017,22(10):12.
[13] Kska P, Stadnik J. Structure-activity relationships study on biological activity of peptides as dipeptidyl peptidase IV inhibitors by chemometric modeling[J]. Chemical Biology & Drug Design,2020,95(2):291−301.
[14] Zhao W Z, Zhang D, Yu Z P, et al. Novel membrane peptidase inhibitory peptides with activity against angiotensin converting enzyme and dipeptidyl peptidaseIV identified from hen eggs[J]. Journal of Functional Foods,2020,64:11.
[15] Nongonierma A B, Mooney C, Shields D C, et al. In silico approaches to predict the potential of milk protein-derived peptides as dipeptidyl peptidase IV (DPP-IV) inhibitors[J]. Peptides,2014,7:43−51.
[16] Wang W, Shen S R, Feng F Q, et al. Pharmacophore-based structure optimization of angiotensin converting enzyme inhibitory peptide[J]. Science in China (Series B: Chemistry),2008,51(8):786−793. doi: 10.1007/s11426-008-0085-5
[17] Iwaniak A, Minkiewicz P, Darewicz M, et al. BIOPEP database of sensory peptides and amino acids[J]. Food Research International,2016,85:155−161. doi: 10.1016/j.foodres.2016.04.031
[18] Yang X X, Chen K N, Liu H G, et al. Purification and identification of peptides with high angiotensin-I converting enzyme(ACE) inhibitory activity from honeybee pupae(Apis mellifera) hydrolysates with in silico gastrointestinal digestion[J]. European Food Research and Technology,2019,245(3):535−544. doi: 10.1007/s00217-018-03223-7
[19] Nongonierma A B, Mooney C, Shields D C, et al. Inhibition of dipeptidyl peptidase IV and xanthine oxidase by amino acids and dipeptides[J]. Food Chemistry,2013,141(1):644−653. doi: 10.1016/j.foodchem.2013.02.115
[20] Li J, Su J L, Chen M, et al. Two novel potent ACEI peptides isolated from Pinctada fucata meat hydrolysates usingin silico analysis: Identification, screening and inhibitory mechanisms[J]. RSC Advances,2021,11:12172−12182. doi: 10.1039/D0RA10476K
[21] Ibrahim M A, Serem J C, Bester M J, et al. The dipeptidyl peptidase IV inhibitory activity and multifunctional antidiabetic properties of SQSPA: Structure-activity analysis evaluated with alanine scanning[J]. International Journal of Biological Macromolecules,2020,160:1220−1229. doi: 10.1016/j.ijbiomac.2020.05.250
[22] Ghayas S, Masood M A, Parveen R, et al. 3D QSAR pharmacophore-based virtual screening for the identification of potential inhibitors of tyrosinase[J]. Journal of Biomolecular Structure & Dynamics,2020,38:2916−2927.
[23] Kandakatla N, Ramakrishnan G. Ligand based pharmacophore modeling and virtual screening studies to design novel HDAC2 inhibitors[J]. Advances in Bioinformatics,2014:e1−11.
[24] 邵鑫, 衣倩颖, 杨春蕾. 基于药效团模型及虚拟筛选方法发现EphB4全新抑制剂[J]. 中国药科大学学报,2016,47(1):38−47. [25] Nongonierma A B, Fitzgerald R J. Inhibition of dipeptidyl peptidase IV (DPP-IV) by tryptophan containing dipeptides[J]. Food & Function,2013,4(12):1843−1849.
[26] Nong N T P, Chen Y K, Shih W L, et al. Characterization of novel dipeptidyl peptidase-IV inhibitory peptides from soft-shelled turtle yolk hydrolysate using orthogonal bioassay-guided fractionations coupled with in vitro and in silico study[J]. Pharmaceuticals,2020,13(10):21.
[27] Nongonierma A B, Fitzgerald R J. Features of dipeptidyl peptidase IV (DPP-IV) inhibitory peptides from dietary proteins[J]. Journal of Food Biochemistry,2019,43(1):11.
[28] Zhu Q S, Chen X J, Wu J J, et al. Dipeptidyl peptidase IV inhibitory peptides from Chlorella vulgaris: in silico gastrointestinal hydrolysis and molecular mechanism[J]. European Food Research & Technology,2017,243:1739−1748.
[29] Xu F R, Yao Y J, Xu X Y, et al. Identification and quantification of DPP-IV-inhibitory peptides from hydrolyzed-rapeseed-protein-derived napin with analysis of the interactions between key residues and protein domains[J]. Journal of Agricultural and Food Chemistry,2019,67:3679−3690. doi: 10.1021/acs.jafc.9b01069
下载:
计量
- 文章访问数:
- HTML全文浏览量:
- PDF下载量:
