一种处理小分子和蛋白质相互作用的计算方法技术

本发明专利技术涉及一种处理小分子和蛋白质相互作用的计算方法，包括以下步骤：(1)读入蛋白质和小分子配体文件数据；(2)判断蛋白质和小分子配体中的原子类型；(3)随机初始化小分子配体的位置和角度；(4)运行SearchProtein算法，优化小分子配体的坐标和角度，Chemscore公式作为评价函数；(5)达到最大迭代次数后终止计算；(6)分析结果，RMSD、对接自由能ΔG和对接位点的坐标，输出结果。与现有技术相比，本发明专利技术具有做到对污染物与生物大分子之间时候有作用作出迅速判断等优点。

【技术实现步骤摘要】

本专利技术涉及一种计算毒理学领域，尤其是涉及。
技术介绍
随着人类社会的日益发展，人类所面对的环境问题愈发严重。各种各样的工、农业废水、废气、废渣排放入人类所居住的公共环境，这些“三废”中往往含有许多对人体有害的化合物。目前，世界上已经有许多国家和地区对“三废”的排放进行了限制和禁止，对环境中污染物的含量做出了严格的限定，制定出相关的法律规范以保护人体免受其伤害。污染物对人体造成损害的过程，是污染物进入体内，与人体中某些重要蛋白结合的过程。这些被结合的蛋白，往往是在人的生命过程中占有重要地位的蛋白。污染物分子与蛋白结合后，有两种行为，其一，使蛋白质的结构发生改变，结构的改变往往引起蛋白质功能改变，使其丧失或改变原有的功能。在通常情况下，污染物的存在会影响蛋白质的二级结构含量的变化；其二，蛋白质的活性位点被污染物占据，一些功能性的分子无法结合，从而导致一些对人体健康有关键作用的物质无法被合成或转化，人体最终将因缺乏这些应该被合成的重要物质而受到损害。污染物与蛋白质相互作用研究的重要性尤为突出。研究的重点在于首先从宏观的层面，污染物对蛋白质造成了何种影响，这种影响如何损害人体的健康。其次从微观层面，污染物与蛋白质造成影响的机理是如何的。微观层面与宏观层面相结合，将使我们获得污染物毒理的全面彻底的视角，使人们规避污染物，并增强人们的环保意识，促使社会各界重视环境污染的问题，加大环境整治和杜绝污染的力度，最终使地球环境更加健康。计算化学，是指通过计算方法解决化学领域中的问题。计算化学中的计算方法涵盖了各种方法。根据科学家研究的体系不同，可采用多种不同的计算方法。最常用的计算方法为分子力学和量子化学计算。计算化学的目的更明确的来说是为了获得研究对象的宏观物理或化学性质。这种宏观性质可分成两类其一，对象的静态平衡性质，例如两个物质作用的结合常数，一个体系的平均势能或是对象在液体中的径向分布函数；其二是动态、非平衡性质，例如对象的相变，大分子的构象变化，小分子在细胞膜或纳米管中的运动。事实上，一个大分子体系的各种性质，是组成该大分子的许许多多微粒表现出的量子态的平均，而获得微粒量子态的唯一途径便是解出薛定谔方程。这是因为薛定谔方程，暗含了微粒的全部性质。在一个有限时间、有限能力的情况下，用薛定谔方程对体系每一个微粒进行从头计算是不可能的，更是没有必要的，所以有必要提出一种更为简单有效地方法，获得体系较为可靠的性质。为了解决上述问题，提出了蒙特卡罗模拟法和分子动力学方法运用到计算化学的领域中。综上，计算方法能够很好地预测对象的静态平衡性质。进而获得可靠的污染物毒性数据
技术实现思路
本专利技术的目的就是为了克服上述现有技术存在的缺陷而提供。本专利技术的目的可以通过以下技术方案来实现，其特征在于，包括以下步骤(I)读入蛋白质和小分子配体文件数据；(2)判断蛋白质和小分子配体中的原子类型；(3)随机初始化小分子配体的位置和角度；(4)运行SearchProtein算法,优化小分子配体的坐标和角度，Chemscore公式作为评价函数；ChemScore评价函数是Eldridge等人在1997年提出的。在本文的程序中，就使用了 ChemScore函数作为能量评价的函数。ChemScore能量函数通过以下公式对配体与蛋白质的能量进行评分Δ Gbinding — Δ G0+ Δ GhbondShbond+ Δ GmetalSmetal+ Δ GlipoSlipo+ Δ GrotHrot(5)达到最大迭代次数后终止计算；(6)分析结果，RMSD (根均方偏移量)、对接自由能Λ G和对接位点的坐标，输出结果O所述的配体文件包括原子的类型、坐标以及键的类型信息。 所述的随机初始化小分子配体的位置和角度具体为小分子配体的位置和角度在活性中心周围ιοΑχιοΑχ IOA的范围内随机产生，并通过Presearch算法进行搜索和对接计算。所述的Presearch算法为将蛋白质的原配体与蛋白质对接,找到活性位点并得到原始对接自由能，读入蛋白质和原分子配体文件数据，判断蛋白质和原配体分子中的原子类型，随机初始化小分子配体的位置和角度；在蛋白质分子上加上极性氢，并计算库尔曼电荷，将其保存后缀名为pdbqs的文件,小分子经过加氢,计算电荷,并确定可旋转键后保存为Pdbq文件；利用mkdpf4命令生成对接所需要的格点参数文件，后缀名为.gpf，打开该文件，将网格中心修改为活性位点中心的坐标；用mkdpf4命令生成dpf文件。所述的SearchProtein算法为进行污染物小分子与蛋白质的分子对接,并得到对接自由能AG，读入蛋白质和原分子配体文件数据，判断蛋白质和原配体分子中的原子类型，随机初始化小分子配体的位置和角度；在蛋白质分子上加上极性氢，并计算库尔曼电荷，将其保存后缀名为pdbqs的文件,小分子经过加氢,计算电荷,确定可旋转键后保存为Pdbq文件；利用mkdpf4命令生成对接所需要的格点参数文件，后缀名为.gpf，打开该文件，将网格中心修改为活性位点中心的坐标；用mkdpf4命令生成dpf文件。与现有技术相比，本专利技术具有可以做到对污染物与生物大分子之间时候有作用做出迅速判断，有利于提升生产效率。该技术关键原理简单，操作建议，容易学习，可以在各个方面得到广泛应用。具体实施例方式下面结合具体实施例对本专利技术进行详细说明。实施例1本专利技术分为Presearch和SearchProtein两个功能模块。主要利用ANSI C++语言编写，可以运行在Windows和Linux操作系统下。程序包中有源程序和头文件18个,主要流程为受体和配体分子的读入，主要为mol2文件格式；受体和配体分子预处理，主要是原子类型，单轴可旋转键的指定等等。PreSearch程序在运行前需要准备presearch, txt文件；然后对小分子进行随机初始化，小分子初始化位置在活性中心周围IOAX IOAX IOA的范围内随机产生；然后利用Presearch进行搜索和对接计算。计算得到对接的RMSD(根均方偏移量)、对接自由能AG和对接位点的坐标。建立蛋白质数据库,利用SearchProtein模块对小分子和大分子进行逆向对接实验。当SearchProtein程序计算完毕后,会生成文件。文件包含了每次对接的能量以及对接位点的坐标。还另有文件分别给出AG与AGO以及原始配体对接评价值与小分子污染物的对接评价值。程序定义了输入文件，主要是定义受体分子文件名、配体分子文件名、初始化中心坐标、初始化参数、搜索算法种类、搜索算法参数等等。除了受体分子文件名、配体分子文件名必须指定外，其他参数可以指定，也可以使用默认值，程序的使用比较简便易用。程序的运行=ToxicitySearch程序里，主要优化的变量主要是小分子配体的中心的三维坐标x，y和z，配体分子的三个笛卡尔角度d，β和Y，其主要运算步骤为(1)读入蛋白质和小分子配体文件，主要是SYBYL mol2文件格式，文件里含有原子的类型、坐标以及键的类型；(2)判断蛋白质和小分子配体中的原子类型；(3)随机初始化小分子配体的位置和角度；(4)运行SearchProtein算法,优化小分子配体的坐标和角度,Chemscore公式作为评价函数；(5)达到本文档来自技高网...

【技术保护点】
一种处理小分子和蛋白质相互作用的计算方法，其特征在于，包括以下步骤：(1)读入蛋白质和小分子配体文件数据；(2)判断蛋白质和小分子配体中的原子类型；(3)随机初始化小分子配体的位置和角度；(4)运行SearchProtein算法，优化小分子配体的坐标和角度，Chemscore公式作为评价函数；(5)达到最大迭代次数后终止计算；(6)分析结果，RMSD、对接自由能ΔG和对接位点的坐标，输出结果。

【技术特征摘要】
1.一种处理小分子和蛋白质相互作用的计算方法，其特征在于，包括以下步骤(1)读入蛋白质和小分子配体文件数据；(2)判断蛋白质和小分子配体中的原子类型；(3)随机初始化小分子配体的位置和角度；(4)运行SearchProtein算法,优化小分子配体的坐标和角度，Chemscore公式作为评价函数；(5)达到最大迭代次数后终止计算；(6)分析结果，RMSD、对接自由能AG和对接位点的坐标，输出结果。2.根据权利要求1所述的一种处理小分子和蛋白质相互作用的计算方法，其特征在于，所述的配体文件包括原子的类型、坐标以及键的类型信息。3.根据权利要求1所述的一种处理小分子和蛋白质相互作用的计算方法，其特征在于，所述的随机初始化小分子配体的位置和角度具体为小分子配体的位置和角度在活性中心周围ιοΑχ ιοΑχ ιοΑ的范围内随机产生，并通过 Presearch算法进行搜索和对接计算。4.根据权利要求3所述的一种处理小分子和蛋白质相互作用的计算方法，其特征在于，所述的Presearch算法为将蛋白质的原配体与蛋白质对接，找到活性位点并得到原始对接自由能，读...

【专利技术属性】
技术研发人员：管翀，赵国华，曹同成，黄晓峰，张莹莹，
申请(专利权)人：同济大学，
类型：发明
国别省市：