一种壳聚糖水解酶可逆抑制保护剂的筛选方法与流程

本发明属于活性物质筛选技术领域，具体涉及一种壳聚糖水解酶可逆抑制保护剂的筛选方法。

背景技术：

酶的特点是用量少、催化效率高、专一性强。酶制剂产业在经历了半个多世纪的迅猛发展后，现已经形成一个富有活力的高技术产业，在食品、纺织、洗涤、造纸、饲料和生物能源等行业都有广泛应用。目前，液体酶制剂技术相比同类产品生产制造过程节约用水约30％～50％、节能50％以上,更符合低碳经济的发展需要；毫无疑问，高活性、高纯度、液体酶将成为今后酶制剂的发展方向。

酶的活性功能决定于其分子结构的完整和严格的构象，环境因素的变化极易导致酶的空间结构破坏，从而丧失其原有的生物活性。一些小分子的物质，可以与酶活中心的氨基酸残基相互作用并进行锁定，对酶的构象进行空间限位，使酶进入休眠的状态。之后在某种激活条件下，小分子物质与酶活性中心的作用力变弱，从而解除下来，活性中心的功能基团重新获得自由伸展，酶的活性得到恢复，这是酶可逆抑制剂作用的基本原理。20世纪60年代初，umezawa最早提出了酶抑制的概念，发展至今，主要用于药物筛选领域。

工业用酶可逆抑制剂没有现成的专一筛选模型和可参考的结构信息数据库，一直以来都是研究领域的难点问题，国内还未能有所突破。而在壳聚糖酶可逆抑制保护剂筛选技术领域还没有有效的方法。

技术实现要素：

本发明的目的是提供一种壳聚糖水解酶的可逆抑制保护剂的筛选方法，筛选的可逆抑制保护剂可提高壳聚糖水解酶的稳定性；从而弥补现有技术存在的不足。

本发明首先建立已知的可逆抑制保护剂训练集，根据已有的壳聚糖水解酶结构数据，进行分子对接并排序，并结合测试结果进行验证，确定方法的可靠性。后建立筛选规则和筛选模型，使用该规则和模型分别对小分子化合物数据库进行筛选、排序和打分，逐渐缩小范围，最终根据由上到下几次排序的结果，确定可以作为壳聚糖水解酶可逆抑制保护剂的化合物，并进行进一步的实验验证。

本发明所提供的壳聚糖水解酶的可逆抑制保护剂的筛选方法，包括以下步骤：

1)根据已知于壳聚糖水解酶的结构数据，确定可逆抑制保护剂的类别及对壳聚糖水解酶的抑制常数，确定训练集的化合物；

2)将步骤1)所述的训练集中的化合物与壳聚糖水解酶进行分子对接，并使用步骤1)中的可逆抑制保护剂对壳聚糖水解酶活性的抑制常数进行修正；

3)根据步骤2)获得的训练集中的化合物与壳聚糖水解酶的结合自由能和可逆抑制常数，结合步骤1)训练集中的化合物筛选结果，制定筛选规则，建立筛选模型；

4)根据步骤3)制定的筛选规则和建立的筛选模型进行筛选。

进一步，所述的步骤1)进入训练集的组成化合物有2,4-二硝基苯甲醚(2,4-dinitro-1-fluorbenzene)，2-羟-5硝基苄溴(2-hydroxy-5-nitrobenzylbromide)，乙酸钠(sodiumacetate)，苯扎氯铵(benzalkoniumchloride)，乙酰亚胺酸乙酯(ethylacetimidate)，乙二胺四乙酸二钠(edta)，葡糖胺(glucosamine)，盐酸胍(guanidiniumhydrochloride)，半胱氨酸(cysteine)，碘乙酸(monoiodoacetate)；这些化合物均是壳聚糖水解酶的抑制剂，包括可逆抑制剂和不可逆抑制剂，与蛋白的结合程度不同

进一步，所述的步骤2)中的分子对接方法为：利用reduce、autodocktools和autodoc4共同完成，首先为壳聚糖水解酶和抑制剂添加氢键和gaussian电荷，再采用分子和壳聚糖水解酶活性位点区域完全柔性的对接方法；

进一步，筛选的规则-6.000≤ei≤-5.000,50≤ki≤200，以此作为筛选规则。

进一步，所述的建立的筛选模型是一个有效的药效团模型，该模型具有2个疏水中心，以及3个氢键受体和5个给体；按照xscore打分，低于5.0分的被排除；

进一步，所述的筛选的方法为按照步骤3)制定的筛选规则从小分子数据库中筛选的化合物组成化合物池a，并进行排序打分，

将化合物池a中的化合物按照理论结合自由能和理论可逆抑制常数重新排序，取排在前2000位的分子使用x-score打分函数联合打分，计算与靶蛋白的亲和力，取打分值大于5.0的化合物再组成化合物池b；

获取化合物池b后，使用筛选模型进行拟合排序。

本发明所建立的壳聚糖水解酶的可逆抑制保护剂的筛选方法，具有以下优点：(1)筛选方法比较快速，可以迅速减少候选分子的数量，减少时间和试验材料的浪费；(2)方法可靠，因为筛选规则和模型的建立是在对已有抑制剂训练集的基础上进行的，最终的筛选结果也使用实验数据验证，理论与实验结果想互补充，可靠性高；(3)针对性好，该方法中使用的分子对接，是建立在壳聚糖水解酶已有结构数据模建基础上(junsheng,xiaofengji,yuanzheng,等.improvementinthethermostabilityofchitosanasefrombacillusehimensisbyintroducingartificialdisulfidebonds[j].biotechnologyletters,2016.)的，针对性强。这种筛选方法，适用于壳聚糖水解酶在工业液体酶制剂上的应用。

附图说明

图1：壳聚糖水解酶结构模型图；

图2：壳聚糖水解酶与抑制保护剂结合后模建结构模型图。

具体实施方式

本发明用下列实施例来进一步说明本发明，但本发明保护范围并非限于下列实施例。

壳聚糖水解酶抑制常数的计算方法：

壳聚糖水解酶催化的米氏常数(km)使用lineweaver–burk双倒数作图法完成，抑制常数ki值通过dixon作图法得到。lineweaver–burk方程可以写作：

二次求导后得到：

从中可以求得抑制常数ki

实施例1壳聚糖水解酶的可逆抑制保护剂的筛选方法

壳聚糖水解酶的可逆抑制保护剂的筛选方法，包括以下步骤：

(1)根据已知壳聚糖水解酶的序列结构、抑制保护剂的结构数据，以及酶与小分子化合物反应的试验数据(表1)，抑制剂包含有乙酸钠(sodiumacetate)，苯扎氯铵(benzalkoniumchloride)，乙酰亚胺酸乙酯(ethylacetimidate)，乙二胺四乙酸二钠(edta)，葡糖胺(glucosamine)，盐酸胍(guanidiniumhydrochloride)，半胱氨酸(cysteine)，碘乙酸(monoiodoacetate)；

表1：壳聚糖水解酶抑制剂反应的抑制常数表

抑制常数(ki)可以反映出化合物与酶结合的能力大小，在这18种化合物中，与壳聚糖水解酶的结合能力太强(ki<0.3)，容易导致酶分子完全变性，被排除掉，其余进入可逆抑制保护剂训练集。

(2)将训练集中的化合物与壳聚糖水解酶进行分子对接

使用的分子对接方法为：利用reduce、autodocktools和autodoc4来完成。首先通过reduce程序添加氢键，然后用autodocktools脚本添加gaussian电荷。选取以坐标点(12.403，-20.749，24.022)为中心，的范围内的空间区域为对接的格点区，利用autogrid来划定格点。采用抑制剂和酶活性位点区域完全柔性的对接方法，利用autodock4来进行对接。打分采用autodock4自带打分和xscore打分相结合的方法。

对接后的结果如表2所示。

表2：训练集分子对接后的理论结合自由能(e)和理论可逆抑制常数(ki)

(3)筛选规则和筛选模型的建立及完善；

获取训练集中的化合物与壳聚糖水解酶的理论结合自由能(ei，kcal/mol)和可逆抑制常数(ki，μm)后，参考表1的实验结果数据，这些化合物对接能量值和理论可逆抑制常数(ki)需要满足选择条件(-6.000≤ei≤-5.000,50≤ki≤200)，因此以此作为筛选规则，在这个范围内的，可逆抑制结果为阳性。

据此建立的筛选模型是一个有效的药效团模型，该模型具有该模型具有2个疏水中心，以及3个氢键受体和5个给体；按照xscore打分，低于5.0分的被排除；

(4)筛选,

首先使用获得的筛选规则(-6.000≤ei≤-5.000,50≤ki≤200)对zinc数据库中大约1900万个小分子化合物进行初步筛选，获得了13000种化合物，组成化合物池a。

将化合物池a中的化合物按照理论结合自由能(ei，)和理论可逆抑制常数(ki)重新排序，取排在前2000位的分子使用x-score打分函数联合打分，计算与靶蛋白的亲和力，取打分值大于5.0的，共计283种化合物再组成化合物池b。

获取化合物池b后，进行筛选模型的拟合。使用以上获得的筛选模型对283种化合物进行匹配，其中与模型中的集团不匹配的或者比模型分子集团的数量多了2倍以上的被去除，最终获得了排名在前3位的结构最为匹配的化合物(表3)：

表3：十种与筛选模型匹配的化合物

这些化合物可以作为壳聚糖水解酶的可逆抑制保护剂的化合物分子。

实施例2化合物zinc01006243与壳聚糖水解酶的对接模拟

选取打分值第一的化合物5,7-diphenyl-1,7-dihydro-pyrrolo[2,3-d]pyrimidine-2,4-dithione进行研究，将其溶解在乙醇中配置成1mm溶液，按照25μm的量加入到浓缩壳聚糖水解酶溶液(8000u/ml)中，酶活收到抑制，保留活性18.76％。40℃保存30d后，按照1：10倍稀释后检测，结果活性为652u/ml，活性得到恢复，相对酶活超过80％。

上述结果表明本发明所建立的筛选方法的理论与实验结果想互补充，可靠性高。