计算机辅助药物设计DPP-IV抑制肽的方法、DPP-IV抑制肽及其应用与流程

计算机辅助药物设计dpp
‑
iv抑制肽的方法、dpp
‑
iv抑制肽及其应用
技术领域
1.本发明属于生物活性肽领域，具体是指一种计算机辅助药物设计dpp
‑
iv抑制肽的方法、dpp
‑
iv抑制肽及其应用。


背景技术：

2.糖尿病(dm)是一种常见的以高血糖症为特征的慢性代谢紊乱。长期高血糖状态可能导致各组织器官，特别是眼、足、肝、肾、心脏等慢性损害与功能障碍，产生严重的并发症，危害生命健康。根据预测，到2035年ⅱ型糖尿病的患病人数将增至5.92亿。用于治疗ⅱ型糖尿病的药物主要有双胍类、磺脲类、糖苷酶抑制剂类、苯甲酸衍生物类和噻唑烷二酮类。但是这些被批准用于治疗t2dm的药物成本过高或对身体有潜在性危害。dpp
‑
iv酶是治疗糖尿病的靶点。在血糖调节的过程中，胰高血糖素样肽
‑
1(glucagon
‑
likepeptide1，glp
‑
1)能通过刺激胰岛素分泌、促进胰岛β细胞的增殖和分化与抑制胰高血糖素分泌等方式降低血糖。dpp
‑
iv酶可分解glp
‑
1。因此，抑制dpp
‑
iv活性从而减少glp
‑
1分解是糖尿病治疗的重要方向之一。
3.现已上市的dpp
‑
iv抑制药有西格列汀(sitagliptin)、维格列汀(vildagliptin)、沙格列汀(saxagliptin)等。这些都是化学合成药，长期使用要加大药量才能维持正常疗效，这有严重副作用或毒性。探索耐药性更低，安全性更高的抗糖尿病药物是药物研发趋势。dpp
‑
iv抑制肽是一种生物活性肽，能抑制dpp
‑
iv酶的活性，而且具有高效、无毒、价廉的优点。目前，对dpp
‑
iv抑制肽的许多研究还停留在表面阶段。获得dpp
‑
iv抑制肽主要通过酶解天然产物蛋白，但使用酶解法制备dpp
‑
iv抑制肽需要耗费大量的人力物力，且提取出的肽未必具有很好活性。
4.随着计算机技术的高速发展和设备快速更新，定量构效关系作为计算机辅助药物设计技术中的一种重要方法，逐渐成为科研人员的研究热点，三维定量构(3d
‑
qsar)效关系是目前定量构效关系中发展最为成熟的技术。使用3d
‑
qsar技术可以建立出化合物结构与生物活性之间的数学模型，为dpp
‑
iv抑制肽的设计、筛选和改造提供理论依据，且研发耗费的人力物力较小，效率更高，生产成本也更加低廉，对研发高活性和高安全性的新型降糖药物具有重要意义。


技术实现要素：

5.为了寻找耐药性更低，安全性更高的dpp
‑
iv抑制肽，并减少传统dpp
‑
iv抑制药在研发过程中巨大的物资消耗，本发明利用计算机辅助药物设计的高效方法设计出6个新型dpp
‑
iv抑制三肽，氨基酸序列为lpi，lpt，lpe，mpl，cpw，wpv。经体外抑制dpp
‑
iv活性实验研究发现这6个肽具有很高的抑制活性，其活性能排在现已报导的dpp
‑
iv抑制肽中的前列，具有潜在的应用价值。
6.本发明的方法通过建立3d
‑
qsar模型来设计dpp
‑
iv抑制肽，主要的步骤包括：建立
数据库并划分成内部训练集和外部测试集；构建数据库中多肽的结构；多肽与dpp
‑
iv进行分子对接以选出评分最高的构象；以三肽ipi作为模板分子并以肽键为公共骨架进行叠合；手动调节各个肽的键角提高模型拟合度；训练集就是我们建立的目标模型，测试集用来测试模型的预测能力；用训练集算出测试集中的肽的pic
50
，以算出的肽的pic
50
与实验得出的pic
50
对比；然后导出训练集的等势图来指导设计dpp
‑
iv抑制肽；最后得到新型dpp
‑
iv抑制肽。
7.本发明通过下述技术方案实现：
8.计算机辅助药物设计dpp
‑
iv抑制肽的方法，包括如下步骤：
9.在已发表的文献中找出dpp
‑
iv抑制肽，筛选出用同一种实验方法得出的ic
50
数据的若干多肽建立数据库，并算出pic
50
值(pic
50
＝
‑
lg ic
50
)作为构建3d
‑
qsar模型的数据基础；
10.将筛选出来的多肽划分成内部训练集和外部测试集；
11.用sybyl2.1.1软件构建筛选出来的多肽的结构，并对其做进一步能量优化；
12.以dpp
‑
iv作为受体蛋白与构建的多肽进行分子对接，选出每个多肽评分最高的构象；
13.以活性最高的三肽ipi作为模板分子，并以多肽的共同结构
‑
肽键为基准进行叠合；采用偏最小二乘法建立3d
‑
qsar模型；采用比较分子场分析法和比较分子相似指数分析法研究结构与活性的关系；通过手动调整多肽的键角使交叉验证系数q2，非交叉验证系数r2，f检验值，标准估计误差都满足模型的标准；
14.训练集即目标模型，通过训练集设计和筛选活性肽，并通过测试集用来测试目标模型的预测能力；用训练集算出测试集中的肽的pic
50
，对算出的肽的pic
50
与实验得出的pic
50
算出差值，若差值的绝对值均小于0.6，说明这个训练集的预测能力好；
15.然后用sybyl
‑
x 2.1.1软件导出训练集的等势图来指导设计dpp
‑
iv抑制肽，最后得到新型dpp
‑
iv抑制肽，对新型dppd
‑
iv抑制肽进行活性预测，筛选出预测活性大的肽进行体外活性试验验证其真实活性。
16.本发明的原理：首先将在文献中找出已发表且是用相同实验方法得到了ic
50
数据的若干dpp
‑
iv抑制肽建立数据库；进行分子叠合，并建立相应的分子表单。采用comfa分析方法计算出立体场和静电场，采用pls将comfa方法中的两个力场以及相关的描述符与药物的生物活性值进行线性回归。然后采用内部样本集建立qsar模型，而利用外部测试集对模型做进一步验证分析。在进行pls建模研究时，先采用留一法进行交叉验证分析得到最佳的自变量数n、q
loo2
及立体场与静电场的贡献值比；然后将得到的最佳自变量数作为基础，通过非交叉验证法再进行分析，得到系列相关参数：s、r2、f，并建立相应的3d
‑
qsar模型。
17.随后利用stedev coeff方法，把立体场和静电场对dpp
‑
iv抑制肽生物活性值的贡献情况用等势图分别呈现出来，以此可以得到设计新的dpp
‑
iv抑制肽类似物的有利信息。采用comsia分析法计算出立体场、静电场以及疏水场，同comfa一样分析得到三维等势图，用于设计新的dpp
‑
iv抑制肽。
18.结合comfa和comsia分析，得出以下预测结果：dpp
‑
iv抑制三肽n端引入大基团，带负电荷的氨基酸有利于提高活性；第二位氨基酸处引入小基团，带负电荷和亲水的氨基酸有利于提高活性；c端为大基团有利于有利提高活性。根据该预测结果设计出了一系列dpp
‑
iv抑制肽，最终利用sybyl软件对新设计的肽进行活性预测，筛选出预测活性高的肽。
19.本发明经过体外实验验证，发现通过本发明的方法设计得到的新型肽具有很高的抑制dpp
‑
iv酶的活性。而且肽类本具有低耐药性，高安全性的优点。所以本发明的dpp
‑
iv抑制肽还可应用于制备治疗糖尿病的药物、或用于改善糖尿病病人症状的保健食品。
附图说明
20.图1是本发明comfa(se)模型的三维等势图；
21.图2是本发明comsia(she)模型的三维等势图；
22.图3是mpl，lpt，lpe的双倒数图；
23.图4是wpv，cpw，lpi的双倒数图。
具体实施方式
24.为使本发明要解决的技术问题、技术方案和优点更加清楚，下面将结合具体实施例进行详细描述。
25.计算机辅助药物设计dpp
‑
iv抑制肽的方法，包括如下步骤：
26.在已发表的文献中找出dpp
‑
iv抑制肽，筛选出用同一种实验方法得出的ic
50
数据的25个三肽和25个二肽，将pic
50
值(pic
50
＝
‑
lgic
50
)作为构建3d
‑
qsar模型的数据基础，按内部训练集内的多肽数是外部测试集内的多肽数的3倍左右这个要求，将50个多肽分子分为内部训练集(37个活性肽分子)和外部测试集(13个活性肽分子)：
27.训练集：ek、gl、sl、al、wrp、wri、wrq、wra、wrm、yp、ws、wry、wrr、wrt、wrw、wrs、wt、wrf、wrk、wrn、wrd、wre、wy、wm、lpl、va、wn、hl、wi、wa、wv、vpf、wrg、wl、wk、wr、ipi。
28.测试集：wrl、vr、lp、wrh、ww、wc、fl、ypy、ip、wq、wp、ypi、vpv。
29.用sybyl2.1.1软件的biopolymer模块构建50种多肽的结构，加入gasetiger
‑
hcukle电荷，采用powell共轭梯度法和tripos力场对分子的能量做进一步优化，其中将最大迭代系数设置为1000，能量收敛限定为0.005kj/mol，利用 1价碳原子作为探针计算出多肽药物的相关力场，其余的参数设置均采用默认值。
30.以dpp
‑
iv作为受体蛋白与构建的多肽序列进行分子对接，选出没种多肽评分最高的构象；
31.以活性最高的三肽ipi作为模板分子，并以上述50个肽的共同结构
‑
肽键为基准进行叠合。采用偏最小二乘法(pls)建立3d
‑
qsar模型；采用比较分子场分析法(comfa)和比较分子相似指数分析法(comsia)研究结构与活性的关系。通过手动调整多肽的键角使交叉验证系数q2，非交叉验证系数r2，f检验值，标准估计误差都满足模型的标准。模型的细节如表1所示。
32.表1comfa和comsia模式的细节
33.参数comfa(se field)comsia(she field)q20.8370.726r20.9910.990n66see0.0590.062
f值574.088518.800立体场贡献值0.6050.246静电场贡献值0.3950.343疏水场贡献值
‑
0.411
34.其中，n为最优主成分数，see为标准估计误差。
35.q2和r2由以下方程式计算：
[0036][0037]
其中y
obs
、y
cvpre
、y
pred
和y
mean
分别是数据库的观测、交叉验证、预测和平均活动值。
[0038]
为了获得最佳的comfa和comsia模型进行手动调整肽的键角：
[0039]
若算出q2小于0.6，则手动调整肽的键角，再算出q2和r2。然后重复叠合，调键角，算q2，r2这三个步骤，直到q2大于0.7，r2大于0.9为止。(模型的q2大于0.5，r2大于0.9即表示模型合格，为了提高模型的预测能力，本发明把comfa和comsia都调到了q2大于0.7，r2大于0.9)。
[0040]
训练集即目标模型，通过训练集设计和筛选活性肽，并通过测试集用来测试目标模型的预测能力。用训练集算出测试集中的肽的pic
50
，对算出的肽的pic
50
与实验得出的pic
50
算出差值，见表2，若差值均小于0.6，说明这个训练集的预测能力好。
[0041]
表2抑制肽分子comfa/comsia模型的pic
50
预测值与实验测量值之间的残差值表
[0042][0043][0044]
然后用sybyl
‑
x 2.1.1软件导出训练集的等势图来指导设计dpp
‑
iv抑制肽，将设
计的dpp
‑
iv抑制肽进行活性预测，筛选出预测活性最高的6个新型dpp
‑
iv抑制三肽，氨基酸序列分别为lpi(leu
‑
pro
‑
ile)，lpt(leu
‑
pro
‑
thr)，lpe(leu
‑
pro
‑
glu)，mpl(met
‑
pro
‑
leu)，cpw(cys
‑
pro
‑
trp)，wpv(trp
‑
pro
‑
val)。
[0045]
在comfa/comsia模型中，分子结构的相关区域和生物活性之间的关系可以通过三维等势线来可视化。本发明以模板分子ipi为参照分子，进行comfa/comsia等势图分析。图1中不同颜色的色块代表的不同的含义。绿色(g)表示引人大基团有利于提高活性，黄色(y)表示引人小基团有利于提高活性；红色(r)表示增加带负电基团有利提高活性，蓝色(b)表示带正电基团有利于提高活性。空间距离分子较远的色块不能提供有效的预测信息。
[0046]
图1中(a)图所示，根据模型得出以下预测结果。在第一位氨基酸残基(n端)处存在一个很大绿色色块(g)，这表明第一位残基处引入大基团有利提高活性。在第二位氨基酸残基处存在一个很大的黄色色块(y)，这表明第二位残基处引入小基团有利提高活性。如图1中b图所示，在第二个氨基酸残基处有一个中等大小的红色(r)色块，这表明第二位残基处引入带负电荷氨基酸也有利提高活性。在第三位氨基酸残基(c端)处存在一个很大绿色色块。这表明第三位残基处引入大基团有利提高活性。
[0047]
图1的(b)图中，显示在comfa的静电场中，蓝色(b)方块表示加入带正电荷的氨基酸有利于提高肽的活性，而红色(r)方块表示引入带负电荷的氨基酸有利于提高肽的活性。在第二个氨基酸残基处有一个中等大小的红色(r)色块，表明在第二个氨基酸残基处引入带负电荷的氨基酸有利于增强活性。
[0048]
在图2的三维等势图中，(a)图表示comsia的s场，得出的结果与comfa模型得出的结果相似。(b)图是comsia表示e场。第一位氨基酸残基(n端)处存在一个很小的红色(r)色块，这表明第一位残基处引入带负电荷氨基酸有利提高活性。在第二位氨基酸残基处存在一个中等大小红色(r)色块，这表明第二位残基处引入带负电荷氨基酸有利提高活性。(c)图是comsia的疏水场(h场)，黄色(y)表示增加疏水性基团有利于提高活性，白色(w)表示增加亲水性基团有利于提高活性。从图中可看出，在第二位氨基酸残基处存在一个很大的白色色块，这表明第二位残基处引入亲水氨基酸有利提高活性。
[0049]
表3为对照肽(ipi)和6新型肽在不同浓度下抑制ddp
‑
iv酶的百分比用probit分析计算方法算出这7个肽的ic
50
值。
[0050]
表3
[0051][0052]
从表中可见，对照肽(ipi)抑制dpp
‑
iv的ic
50
为3.66um，与许多文献报道的值一致。lpi、mpl和lpt抑制dpp
‑
iv的ic
50
分别为7.51、20.73和22.88um。这三种肽的活性在已报道的数百种dpp
‑
iv抑制肽中最高。lpi是三种多肽中活性最强的，其ic
50
水平达到1位。在已报道的数百种dpp
‑
iv抑制肽中，它排名第三。
[0053]
lpi的分子式为c
17
h
31
n3o4，化学结构简式为：
[0054][0055]
lpt的分子式为c
15
h
27
n3o5，化学结构简式为：
[0056][0057]
lpe的分子式为c
16
h
27
n3o6，化学结构简式为：
[0058][0059]
mpl的分子式为c
16
h
29
n3o4s，化学结构简式为：
[0060][0061]
cpw的分子式为c
21
h
27
n4o4，化学结构简式为：
[0062][0063]
wpv的分子式为c
19
h
23
n4o4s，化学结构简式为：
[0064][0065]
根据设计出的序列，由专业机构对该6个肽进行合成，纯度为99％，用于体外实验验证。
[0066]
实验材料
[0067]
ddp
‑
iv抑制肽lpi，lpt，lpe，mpl，cpw，wpv(纯度>99％)，由上海强耀生物科技有限公司合成；dpp
‑
iv/cd26,human recombinanty由biovision公司提供；甘氨酰
‑
脯氨酰
‑
对硝基苯胺(gly
‑
pro
‑
pna)和diprotin a(ipi)由cayman公司提供；盐酸三羟甲基氨基甲烷缓冲液(tris
‑
hcl buffer)由上海源叶生物科技有限公司提供。
[0068]
dpp
‑
iv抑制实验
[0069]
dpp
‑
iv抑制试验步骤如下：将测试样品肽以10至500um的浓度溶解在tris
‑
hcl缓冲液(100mm,ph 8.0)中。将测试样品(50μl)加入96孔微孔板上,再加入反应底物gly
‑
pro
‑
pna(终浓度0.2mm)。酶标仪的孵育器调至温度37℃，在酶标仪中震荡10分钟混匀。然后通过添加dpp
‑
iv(终浓度0.01u/ml)来引发反应。所有试剂和样品均在tris
‑
hcl缓冲液(100mm,ph 8.0)中稀释。diprotin a用作阳性对照。用tris
‑
hcl缓冲液(50μl)替代测试样品溶液作为阴性对照。每个样品进行三次重复实验。将微孔板在37℃下孵育15分钟，在405nm处测试释放的pna的吸光度。dpp
‑
iv的ic
50
值(使dpp
‑
iv活性下降一半时所需的肽的浓度)是通过将抑制百分比作为受试化合物浓度的函数作图来确定的。使用双倒数作图法(lineweaver
‑
burk作图法)研究了这六个肽的抑制模式。
[0070]
根据lineweaver
‑
burk作图法对mpl，lpt，lpe和wpv，cpw，lpi做双倒数图，得到图3和图4。原理是在一系列底物浓度(0.2至1mm)下，加入和不加入抑制剂对酶促反应的影响来判断竞争性抑制剂和非竞争性抑制剂。其中抑制剂的浓度是其ic
50
浓度。抑制肽mpl，lpt，lpe，cpw，lpi的直线与无抑制剂的直线相交于y轴。说明mpl，lpt，lpe，cpw，lpi是dpp
‑
iv的竞争性抑制剂。而抑制肽wpv的直线与无抑制剂的直线相交于x轴，说明wpv是dpp
‑
iv的非竞争性抑制剂。
[0071]
利用本发明的模型设计的多肽对dpp
‑
iv具有很高的抑制活性，如表4所示，预测的pic
50
值与实验结果基本一致。这表明我们的建模方法是非常正确的。而且所建立的模型预测能力强，可靠性高。
[0072]
表4本发明6个新型dpp
‑
iv肽的pic
50
[0073][0074]
以上所述是本发明的优选实施方式，应当指出，对于本技术领域的普通技术人员来说，在不脱离本发明所述原理的前提下，还可以做出若干改进和润饰，这些改进和润饰也应视为本发明的保护范围。

技术特征：
1.计算机辅助药物设计dpp
‑
iv抑制肽的方法，其特征在于，包括如下步骤：在已发表的文献中找出dpp
‑
iv抑制肽，筛选出用同一种实验方法得出的ic
50
数据的若干多肽建立数据库，并算出pic
50
值作为构建3d
‑
qsar模型的数据基础；将筛选出来的多肽划分成内部训练集和外部测试集；用sybyl2.1.1软件构建筛选出来的多肽的结构，并对其做进一步能量优化；以dpp
‑
iv作为受体蛋白分别与构建的多肽进行分子对接，并选出每个多肽评分最高的构象；以活性最高的三肽ipi作为模板分子，并以多肽的共同结构
‑
肽键为基准进行叠合；采用偏最小二乘法建立3d
‑
qsar模型；采用比较分子场分析法和比较分子相似指数分析法研究结构与活性的关系；通过手动调整多肽的键角使交叉验证系数q2，非交叉验证系数r2，f检验值，标准估计误差都满足模型的标准；训练集即目标模型，通过训练集设计和筛选活性肽，并通过测试集用来测试目标模型的预测能力；用训练集算出测试集中多肽的pic
50
，并用此pic
50
与实验得出的pic
50
算出差值；若差值的绝对值均小于0.6，说明这个训练集的预测能力好；然后用sybyl
‑
x 2.1.1软件导出训练集的等势图来指导设计dpp
‑
iv抑制肽，最后得到新型dpp
‑
iv抑制肽。2.如权利要求1所述的计算机辅助药物设计dpp
‑
iv抑制肽的方法，其特征在于，筛选出来的用于建立模型的多肽数据库包括25个三肽和25个二肽：ek、gl、sl、al、wrp、wri、wrq、wra、wrm、yp、ws、wry、wrr、wrt、wrw、wrs、wt、wrf、wrk、wrn、wrd、wre、wy、wm、lpl、va、wn、hl、wi、wa、wv、vpf、wrg、wl、wk、wr、ipi、wrl、vr、lp、wrh、ww、wc、fl、ypy、ip、wq、wp、ypi、vpv。3.如权利要求2所述的计算机辅助药物设计dpp
‑
iv抑制肽的方法，其特征在于，训练集有37个多肽：ek、gl、sl、al、wrp、wri、wrq、wra、wrm、yp、ws、wry、wrr、wrt、wrw、wrs、wt、wrf、wrk、wrn、wrd、wre、wy、wm、lpl、va、wn、hl、wi、wa、wv、vpf、wrg、wl、wk、wr、ipi；测试集有13个多肽：wrl、vr、lp、wrh、ww、wc、fl、ypy、ip、wq、wp、ypi、vpv。4.如权利要求1所述的计算机辅助药物设计dpp
‑
iv抑制肽的方法，其特征在于，能量优化中，将最大迭代系数设置为1000，能量收敛限定为0.005kj/mol，利用 1价碳原子作为探针计算出多肽药物的相关力场，其余的参数设置均采用默认值。5.如权利要求1所述的计算机辅助药物设计dpp
‑
iv抑制肽的方法，其特征在于，在叠合后，为了获得最佳的comfa和comsia模型，还包括手动调整肽键角的操作：若算出q2小于0.6，则手动调整肽的键角，再算出q2和r2；然后重复叠合，调键角，算q2，r2这三个步骤，直到q2大于0.7，r2大于0.9为止。6.如权利要求1～5任一项所述的计算机辅助药物设计dpp
‑
iv抑制肽的方法，其特征在于，还包括对得到的新型dpp
‑
iv抑制肽进行活性预测，筛选出预测活性高的六个三肽，其氨基酸序列分别为lpi，lpt，lpe，mpl，cpw，wpv。7.一种如权利要求1～5任一项所述的计算机辅助药物设计dpp
‑
iv抑制肽的方法得到的dpp
‑
iv抑制肽，包括lpi，lpt，lpe，mpl，cpw，wpv。8.一种如权利要求1～5任一项所述的计算机辅助药物设计dpp
‑
iv抑制肽的方法得到
的dpp
‑
iv抑制肽在制备治疗糖尿病的药物上的应用。9.一种如权利要求1～5任一项所述的计算机辅助药物设计dpp
‑
iv抑制肽的方法得到的dpp
‑
iv抑制肽在制备改善糖尿病病人症状的食品上的应用。
技术总结
本发明提供一种计算机辅助药物设计DPP


技术研发人员：张荣 梁振 严文利 李红梅 李颖
受保护的技术使用者：广东药科大学
技术研发日：2021.04.08
技术公布日：2021/6/25