10分钟教你掌握分子对接模拟软件(医药向)

首先介绍一下自己吧本人毕业于南方某知名大学药学系目前于澳门科技大学攻读硕士研究生。从本科开始自己就在接触CADD（计算机辅助药物设计）方面的软件知识在此将分享一些自己的纯干货下面将以一个实例操作带大家迅速认识和掌握分子模拟对接望给各位从事医药行业和药物化学合成的同学带来帮助。话不多说，下面进入正题。首先我们搞清楚一个概念什么是分子模拟对接分子模拟对接简单来说就是利用电脑软件将受体蛋白与配体分子进行模拟对接，计算它们的结合能KJ/MOL）大小来判断结合是否紧密若结合效果比较理想那么该蛋白受体或配体则是我们理想的分子以进一步进行实验室操作免盲目实验带来的人力经济损失。接下来我将介绍一下本篇文章的主角是我们所要用到的软件Chemdraw、AutodockTools以及PyMol。为了便于理解，简要概括之Chemdraw为化合物分子绘图软件；PyRx为Vina算法搭载软件，能够调用其算法直接进行模拟对接；AutodockTools是为对接结果成像软件，可以进一步分析其结构。下面正式进入正题我将大致分为三个板块来进行推进受体配体的准备分子对接；结果分析。研究类型为：已知若干配体分子结构，通过受体蛋白测试配体分子活性。本次筛选意在以COMT酶为受体20种与常见氨基酸形成环二肽的目标化合物中筛选出与COMT酶受体结合最为紧密的一种环二肽结构，大大减少了随机筛选的盲目性利于进一步研究该类化合物分子的生物学活性与改造成抗帕金森疾病前药的可能。图1示了20种不同环二肽结构物质的统一结构，随着R基团的不同，所对应的氨基酸也不同。而表1则展示了种不同环二肽的分子式。图1Cycol[DOPA）-AA]

表1待选的20种配体分子配体名称（）-Ala]（）-Arg]（）-Asn]（）-Asp]（）-Cys]（）]（）]（）-Gly]（）-His]（）]（）-Leu（）]（）]（）-Phe]（）]（）-Ser]（）（）-Trp]（）-Tyr]（）]一受配的备

分子量307.079351.089350.104351.089339.145364.131365.116293.052373.141349.16349.16364.17365.199383.177333.117323.087337.105422.213399.176335.133首先谈谈受体分子的准备工作一般来讲受体分子作为目标靶点在我们人体内通常以大分子蛋白质的形式存在而PDB数据库则覆盖了世界上70%的人源鼠源等的蛋白数据我们打算在该网站上下载一个简单不含辅酶的蛋白我们选择了代码为4PYI的酶（基于配体分子靶标为酶从而选择之

图24PYI型COMT酶4PYI蛋白受体选择PDB下行另存为路径最好选择C:\Users\Administrator\Desktop以避免Windows操系统不兼容导致的程序运行崩溃下载文件导入AutodockTools进行受体前处理。1.除水：Autodocktools菜单栏Edit中点击除去蛋白分子中的水分子，避免其对分子对接的干扰。2.加氢AutodockTools菜单栏Edit中点击Hydrogens—Add参数设置PolarOnly；Nobongorder；Renumberatomstoincludenewhydrogens3.设置原子类型AutodockTools菜单栏中点击—AssignAD4type。4.保存：File—Save—WritePDBQT。保存为PDBQT格式便于分子对接。配体是分子对接的变量虚拟筛选的目标群体次筛选的20种（DOPA）需要手动绘制加以表现，本次用到的绘制工具为Chemdraw与Chem3D，绘制完毕后对配体进行后处理。1.绘制：利用20种Cyclo（DOPA（6-NO2）-AA）分别进行绘制，再将绘制结果直接导入Chem3D输出为mol2文件。2.能量最小化：分别将不同mol2文件导入PyRx软件，在右下角菜单栏选择MinimizationParameters”进行能量最小化处理，参数设置ForceField—uffOptimizationAlgorithm—ConjugateGradients；Totalnumberofsteps—200；Numberofstepsfor—1；Stopifenergyislessthan—0.001。最后将能量最小化的分子另存为原格式。

格式转换：将经过能量最小化后的mol2分子通过LigandInput导入AutodockTools，然后LigandOutput选项里选择SaveasPDBQT，最终将mol2文件格式转换为PDBQT格式（与受体格式相同二分对经过受体与配体分子的准备后可以得到两种分子的PDBQT格式的文件等待进一步分子对接，分子对接在软件中调用AutodockVina算法进行模拟对接。1.导入将20种配体分子导入PyRx软件的选项并全部选中另外将已进行处理的4PYI酶导入Macromolecules选项并选中。2.对接将20种配体与4PYI受体选中后选择计算网格大小参数设置CenterX：-20.757Y：-22.666Z:0.8213;Dimensions(Angstrom)X:55.0259Y:55.0259Z:55.0259；Exhaustiveness：30。最后RunVina等待计算结果。图3PyRx分对接操作界面

三虚筛的果析结表220种体分子结合能数据LigandBinding（kcal/mol）（）-Ala]（）-Arg]（）-Asn]（）-Asp]（）-Cys]（）]（）]（）-Gly]（）-His]（）]（）-Leu（）]（）]（）-Phe]（）]（）-Ser]（）（）-Trp]（）-Tyr]（）]

--8-8-9由表格数据可知Cycol[DOPA（6-NO2结合能为-9.2kcal/mol过PyMOL展现其氢键作用力良好，该分子共与五个氨基酸分子发生氢键相互作用，且在

surface模式下发现该分子与蛋白质受体凹槽部分结合紧密所以该分子属于最佳目标化合物。图4与COMT互作用图5的白结合位点本次虚拟筛选将二十种含常见氨基酸的配