分子对接研究英文论文

Docking算法需要每个原子带有电荷并且需要标记原子的属性。这些信息通常未包含在PDB文件中。我们需要在对蛋白和小分子的PDB文件预处理，生成PDBQT文件同时包含以上信息和PDB文件中的原子坐标信息。进一步地对于“柔性配体docking”，我们还需要定义配体的柔性部分和刚性部分。所有这些都可以通过软件AutoDock Tools (adt)来完成。 =====准备受体蛋白=======    注：在windows下，我们可以手动选择，或者利用Excel的筛选功能。在Linux下，使用命令egrep "^(ATOM|TER)" >   2. 启动AutoDockTools   3. 依次点选加载蛋白分子。   按住左键拖动旋转分子结构；点击中键滚动缩放；按住右键移动晶体位置。  4. 更改展示方式：依次点选Color-By Atom Type-All Geometries-OK。  5. 加氢：晶体结构中通常缺少氢原子的坐标 (因为氢原子电子少，且质子核对电子吸引能力弱，因此很难定位，具体见)。但是在docking过程中，氢原子尤其是极性氢原子对计算静电作用是必须的。因此我们需要给蛋白加上氢原子，依次点选Edit-Hydrogen-Add-Polar only-OK。这时氢原子会以白色短线形式出现。  6. 存储对蛋白的每个原子所做的修改和原子类型判断：依次点选Grid-Macromolecule-Choose-1HSG_protein-Slect Molecule。ADT会弹出一个信息框包含程序所做的处理，比如合并非极性氢原子，计算原子局部电荷和判断原子类型，并提示保存。打开文件，查看最后两列，分别为每个原子的电量和类型 。如下图所示分别是原始的pdb文件和新生成的pdbqt文件。7. 在受体蛋白定义配体结合的3D搜索空间: 如果我们事先不知道结合位点，理论上可以定义一个长方体盒子包含整个蛋白或者随便一个特定区域 。 依次点选Grid-Grid box将会在蛋白上画出一个长方体，并且有一个弹出框。在弹出框中，拖拽刻度线查看长方体的变化，完成设置。在这个例子中，我们知道结合位点，就选取以其为中心的一个小空间。设置Spacing (angstrom)为1埃 (这实际是一个换算系数, 相当于步长; 默认为，是C-C单键长度的1/4，最大为1。spacing值与(各个维度上的点的数目+1)的乘机就是长方体Grid box的大小)。在我们调整的过程中，可以看到随着这个数值的变大，立方体也被放大了。另外我们设置x,y,z center为16,25,4,number of points in (x,y,z)-dimension为30,30,30(最大为126，必须为偶数，AutoDock会自动再每一维再加一个点)。记下我们设置的这些点，下面会用到。在刻度转盘处点击右键会弹出一个窗口，输入数字回车即可设置GRID的中心坐标和大小。较大的number of points in (xyz)-dimension和较小的Spacing会增加搜索的精度，同时需要花费更多的计算时间。   8. 设置受体的柔性残基：在ADT中依次点选Flexible Residues-Input-Choose Macromolecule-1hsg_prot; select-select from string-Residue: ARG8-Add-Dismiss, 8号ARG氨基酸残基就被选中了。  再依次点选Flexible Residues-Choose Torsions in Currently Selected Residues将选择的残基标记为柔性残基并设置可扭转的数量。在分子显示窗口中分别点击两个残基上CA和CB原子之间的建，使之变为非扭转的（紫色显示），这样两个残基中的32个键中有6个是可扭转的。这里设置配体的柔性残基或者使CA-CB的键为刚性都是可选操作。Flexible Residues-Output-Save Flexible PDBQT保存柔性残基文件。Flexible Residues-Output-Save Rigid PDBQT保存柔性残基文件。 9. 关掉grid和删除protein：Grid Options-File-Close w/out saving; Edit-Delete-Delete Molecule-1hsg_prot-Continue。 =====准备配体===== 蛋白结构类似，配体的结构也缺少氢原子，我们需要添加氢原子并且定义哪些键是可以旋转的以用于柔性docking。从PDB结构中提取配体的原子位置。indinavir的配体残基名字为MK1，以HETATM开头的行表示非核心多聚体的成分 (heteroatoms)。 Linux系统下，运行grep "^HETATM.*MK1" > Windows系统下，直接拷贝到文件   3. 将结构读入ADT；依次点选;Color-By Atom Type-All Geometreies-OK。  4. 晶体结构中通常缺少氢原子 (因为氢原子电子少，且质子核对电子吸引能力弱，因此很难定位)。但是在docking过程中，氢原子，尤其是极性氢原子对计算静电作用是必须的。因此我们需要给配体加上氢原子，Edit-Hydrogen-Add-Polar only-OK(之所以选择Polar only是因为vina的官方视频里面是这么选择的)。这时氢原子会以白色短线形式出现。 5. 在ADT中定义此化合物为配体，以便ADT为其计算局部电荷(partialcharges)和设置可旋转配体键。依次点选Ligand-input-Choose-indinavir-Select Molecule for AutoDock4。这时会有一个弹出框显示ADT所做的操作，包括合并非极性氢(只在添加了的情况下)、计算电荷电量和设置旋转键。然后点选Ligand-Output-Save as PDBQT存储结果。   6. 查看ADT检测出的旋转键，依次点选Ligand-TorsionTree-Choose Torsions，可以看到Number of rotatable bonds=14/32。====准备docking的配置文件==== docking配置文件包含了输入的受体(蛋白)、配体(化合物)和搜索参数的信息，为一个文本文件，名字任意，可以为，内容如下   receptor = ligand  = num_modes = 50 out = log = center_x = 16 center_y = 25 center_z = 4 size_x = 30 size_y = 30 size_z = 30 seed = 2009   receptor和ligand为输入文件的名字，与在同一目录下; out为输出文件的名字；log为输出日志文件的名字。centerhe和size定义搜索空间的位置和大小。num_modes设置最多显示的结合模型，鉴于只输出符合能量值要求的结果，最后输出的结合模型数量可能少于这一数值。seed设置随机数生成的种子，可以为任意整数。如果想重现之前的分析结果就需要使用相同的seed。 ====Docking 小分子化合物indinavir到HIV-1蛋白酶====== 使用AutoDockVina执行docking预测。  在windows命令行提示符或linux终端下运行命令  vina--config ，大概需要几分钟时间。  2. 输出结果包含两个文件，构象文件和日志文件。 : 包含所有docking的模式，通常第一个为结合最好的构象，但如果前几个能量值相差不大时也有例外。 : 日志文件，包含结合能量值(第一列，越低越稳定，默认由低到高排序，所以第一个为最好的构象)、每个构象与第一个构象的距离、每个构象与第一个构象的差别。 ====用PyMOL可视化docking结果====   打开PyMOL，依次点选File-Open文件类型选择All Files-选取结果文件、原始蛋白和配体的pdb文件、原教程的pdbqt文件。   这个是原始的蛋白受体和配体信息：受体蛋白信息：配体信息：预测的配体信息：

不可以。Jouralofmedicinechemistry药物化学杂志。分子对接的种类主要包括：1、刚体对接：指在对接过程中，研究体系（受体和配体）的构象不发生变化?。适合考察比较大的体系，如蛋白质和蛋白质间以及蛋白质和核酸之间的对接。2、半柔性对接：指在对接过程中，研究体系尤其是配体的构象允许在一定的范围内变化。适合处理大分子和小分子间的对接，对接过程中，小分子的构象一般是可以变化的，但大分子是刚性的。所以Jouralofmedicinechemistry的分子对接放在SI里不可以。3、柔性对接：指在对接过程中，研究体系的构象基本上可以自由变化的。一般用于精确考虑分子间的识别情况。由于计算过程中体系的构象可以变化，所以计算耗费最大。

2001年国内在药物设计和分子模拟领域取得的主要进展和成就。总的说来，三维定量构效关系方面的工作较多，并且与实验紧密结合；分子模拟领域的工作相对水平较高，既有方法方面的研究，也有在具体体系中的应用；而由于知识产权方面的原因，目前基于结构的药物设计与虚拟筛选等方面的工作发表较少。 关键词：三维定量构效关系，分子动力学，基于结构的药物设计，分子对接 计算机模拟已经成为继实验、理论之后的第三种研究手段，在科学的各个领域都得到了广泛的应用，并且发展迅速。药物设计与分子模拟也是如此，其总体趋势是计算精度越来越高，与实验结合日趋紧密。近年来，国内药物设计与分子模拟研究工作取得了长足的进展，其中，定量构效关系在药物研究中的应用广泛，已成为药物结构改造的一个有力工具；分子模拟研究发展较快，计算的方法多样化，模拟的时间、空间尺度也越来越接近真实体系，取得了一些水平较高的成果。 一、定量构效关系 定量构效关系是最早发展起来的一种分子设计方法，一旦确定已知的化合物的生物活性与其结构之间的定量关系，研究者就可以依据这一关系对化合物结构进行改造，以获得更高的活性。3D-QSAR（尤其是比较分子力场分析方法）由于其预测能力强，模型形象、直观，已成为最常用的药物设计方法之一。 噻唑烷二酮和芳酮酸类等PPAR激动剂，能够增强胰岛素与其受体结合的信号传导，促进脂肪代谢。利用比较分子力场分析方法对这两类化合物进行了三维定量构效关系研究，揭示了芳酮酸类化合物比噻唑烷二酮类化合物活性更高的原因，并推测PPAR与这一类化合物结合时，活性位点的形状可能和与罗咯列酮结合时不一样[1]。 双叔丁基苯基杂环化合物是COX-2的抑制剂，具有抗炎作用。对三类含双叔丁基苯基杂环化合物进行三维定量构效关系研究，在三种可能的叠合模型中，选择了交叉验证结果最好的一种。研究结果表明：静电作用对化合物的活性贡献稍大于范德华作用；增加杂环的正电荷，可以提高化合物的活性[2]。 采用比较分子力场分析方法对1－咪唑基苯并二氮杂卓（草头卓）类法尼基蛋白转移酶抑制剂进行研究，发现在苯并二氮杂卓（草头卓）中的苯环上增加取代基的体积会使化合物活性降低，在咪唑环上增大取代基的体积能够提高化合物的活性[3]。 吡里酮类化合物是一类很有药用前途的抗炎镇痛剂，利用计算机辅助药物设计专家系统（Apex-3D）研究了5,6－二芳基取代吡里酮类化合物的抗炎活性与其分子结构之间的定量关系。研究表明：化合物的活性与吡里酮环的1位取代基的电荷密度正相关，而与分子的疏水性负相关[4]。 采用比较分子力场分析方法分析了一组－羟基丁烯羟酸内酯类内皮激素拮抗剂的结构活性关系。在计算中，除考察静电场和立体场之外，还引入了氢键场。结果表明，引入氢键场可以明显地改善计算结果。另外，还分别考察了格点距离以及探针原子类型对计算结果的影响[5]。 采用一系列几何结构参数以及由ClogP软件计算得到的logP值作为结构参数，对21个钾离子通道开放剂进行了定量构效关系研究。研究结果表明，氢键受体原子与疏水部分中心距离HA-Hy，以及芳环中心垂线和芳环中心与疏水部分中心连线的夹角对动脉舒张活性影响较显著[6]。 利用CATALYST软件得到了酪氨酸激酶抑制剂的三维药效基团模型，所研究的两类化合物具有完全不同的骨架，但其药效基团相似，说明两者的结合模式相似。这两类化合物的药效基团包括一个氢键受体，一个氢键给体，一个脂肪类疏水团和一个芳香类疏水团。在此基础上，进行了三维构效关系研究[7]。 以黄酮为模板，对一系列黄酮类化合物进行了定量构效关系研究及从头算研究。根据CoMFA, CoMSIA和HQSAR分析以及静电势计算结果，提出了合理的药效基团模型，并提出了一种黄酮类化合物与GABAA受体结合的模型[8]。 二、基于结构的药物设计 基于结构的药物设计是寻找新结构的药物先导化合物的方法之一，国内运用基于结构的药物设计方法发现活性化合物的例子也越来越多，不过，因为知识产权等方面的原因，这方面的论文发表较少，因此本文仅能介绍方法的发展。 在对已知抑制剂与磷脂酶A2的复合物晶体结构的分析的基础上，设计了新的抑制剂，并利用SCORE程序对其生物活性进行了预测。随后进行了化学合成及药理测试，结果表明设计的两种化合物都有一定的活性，并且实验值与SCORE的预测值具有对应关系[9]。 利用同源模建方法、序列进化踪迹分析，研究了蛋白质前提加工酶furin/kexin与水蛭抑制剂（eglin C）突变体的相互作用模式，阐明了furin/kexin各个亚类的底物/抑制剂的共性和差异性。在此基础上，利用界面再设计策略进行了furin/kexin抑制剂的理性再设计[10]。 对DycoBlock(Dynamically Assembling Molecular Building Block)方法进行了改进，主要包括利用结合能和溶剂可及性表面积(SASA)对DycoBlock产生的分子进行筛选，并发展了一个基于分子相似性的分类方法。利用COX-2及其选择性抑制剂SC-558对其方法进行了测试，重现晶体结构的比率由以前的上升到[11]。在DycoBlock的基础上发展了F-DycoBlock，在进行多拷贝随机分子动力学的过程中，一方面考虑小分子配体柔性，同时考虑小分子配体的平均场对大分子的影响，因而还可以考虑蛋白质的柔性。利用COX-2以及HIV-1蛋白水解酶两个体系对其进行了测试[12]。 三、分子模拟 分子模拟通过对生物大分子的结构和动力学行为的模拟，为阐明其结构、生物功能提供了丰富的信息。一般来说，这些信息用实验方法获取是比较困难的，因此，分子模拟已成为研究生物大分子结构和功能的一种重要手段。 Threading方法对于根据同源性比较低的模板蛋白预测目标蛋白三维结构已经取得一定的成功，但仍难以处理那些模板蛋白以及目标蛋白三维结构相似，而二级结构在序列中的顺序不一致的情况。研究发现预先对目标序列中的交错的二级结构片段作适当的调整，可以显著地提高Threading方法的预测能力[13]。 利用蒙特卡罗方法对蛇毒神经毒素的八种不同Loop III序列在自由态以及与短链神经毒素埃拉布毒素b（Erabutoxin b）结合状态的构象性质进行研究，发现远程作用确实对蛋白质表面的Loop区的构象及稳定性有影响，而近程作用只是提供候选的构象，还须经过远程作用的过滤[14]。 BLIP是目前效果最好的TEM-1-内酰氨酶的蛋白类抑制剂，它的结构域1中的一个转角结构B1可以与TEM-1-内酰氨酶的活性部位结合，从而抑制其活性。利用圆二色谱、傅立叶红外光谱和核磁共振谱研究了B1在磷酸盐缓冲溶液中的构象，结果表明，B1在溶液中形成了转角结构 [15]。 采用分子动力学方法对猪胰岛素二聚体进行模拟，并利用连续介质模型有限差分方法计算二聚体的静电势，利用溶剂可及性表面分析二聚体的疏水性，还考察了不同的PH值对二聚体的静电作用和疏水作用的影响。分析表明，在弱酸弱碱条件下，二聚体的静电作用和疏水作用比较有利；在PH值为时，二聚体的疏水作用最佳[16]。 利用分子动力学方法对喜树碱进行构象搜寻，并利用AM1，PM3以及ab initio方法对与晶体结构比较接近的构象进行优化。ab initio方法优化得到的结构与晶体结构更接近，喜树碱E环内酯部分形成了分子内氢键，使其内酯键反应活性提高，使其易受水的影响开环而失去活性。研究还表明，喜树碱在与TopoI-DNA复合物作用时，要破坏分子内氢键，能量补偿较大[17, 18]。 采用快速退火进化算法对氨基酸与-环糊精的结合进行了模拟。对14种包合物的结合常数与作用能的线性回归分析表明，两者之间存在非常好的相关性。通过对作用能进行分析，发现范德华作用起主要作用。对不同手性的氨基酸与-环糊精的作用能进行了比较，解释了-环糊精的手性选择性[19]。 利用自由能微扰方法计算了两种hydroxamate类的抑制剂和MMP-3的相对结合自由能。自由能计算采用了慢增长和固定窗口增长两种方法，并且都采用了双向采样策略。计算得到的结合自由能与实验值吻合得很好，同时从动力学模拟的轨迹分析了抑制剂与MMP-3之间的重要相互作用[20]。 通过同源模建方法得到了枯草杆菌（Bacillus subtilis）色氨酸tRNA合成酶的三维结构，分析结果表明合成酶的二级结构含有16个螺旋和5个折叠，唯一的色氨酸Trp92位于两个亚基的界面上。另外，还对配基结合位点、活性位点以及与tRNATrp可能的结合方式进行了预测[21]。 蛋白质的折叠是一个自由能驱动的过程，天然的蛋白质构象处于一个热力学的低能量态。基于蛋白质结构域是折叠单位的设想，蛋白质结构域作为一个独立的折叠单元，应该不仅在结构上相对致密，而且在能量上也应该相对稳定。因此，用折叠自由能来划分结构域更为合理。利用去折叠自由能的方法对50种不同的双结构域蛋白质进行了结构域分析，大多数划分结果与文献一致；少数结果虽然与文献不一致，但去折叠自由能方法的划分结果可能更合理[22]。 采用布朗动力学对蝎毒Lq2与钾离子通道的结合进行模拟，蝎毒的三维结构取自NMR实验测定的22个结构，结果表明只有4个结构能够以高概率和较好的静电作用与钾离子通道结合。对上述4种结构进行三联体接触分析，构建了复合物的三维结构模型，并确定了Lq2与钾离子通道之间的主要相互作用[23]。 采用同源模建和分子动力学方法，模建了memapsin 2蛋白水解酶的三维结构。memapsin 2蛋白水解酶的活性部位为一个富含-折叠的结构域所支撑，活性部位通过一段螺旋与该结构域连接。-发夹结构中的Tyr132是天冬氨酸蛋白水解酶的共同特征。分析表明，虽然M2蛋白水解酶的结构比较独特，但仍属于天冬氨酸蛋白水解酶超家族[24]。 四、数据库与虚拟筛选 如果已经掌握了足够的有关靶标以及与靶标结合的化合物的信息，针对大的化合物数据库的虚拟筛选非常具有吸引力，并可作为一种与高通量筛选互为补充的寻找先导化合物的方法。最近，在超级计算机上对Kuntz发展的DOCK程序进行了并行化，大大加快了数据库筛选的速度。采用这种方法，从商用的化合物数据库以及中国天然产物数据库中找到了一系列对钾离子通道，-分泌酶以及PPAR等靶标具有生物活性的化合物。但由于要申请专利，结果还没有发表 一个包含范德华作用、金属配位作用、氢键作用、去溶剂化作用以及小分子构象变化等项的打分函数SCORE，被用于评价DOCK4找到的小分子与蛋白质的结合构象。用两者结合得到的程序SCOREDOCK对200个蛋白－配体复合物进行测试，发现SCORE能够明显地提高DOCK计算精度[25]。 通过Fortran PowerStation和Visual C++混合编程，将在DOS下运行的、用于研究分子结构和生物活性关系的计算机辅助筛选生物活性化合物软件（Computer-Aided Screening bioActive Compounds, CASAC）移植到Windows系统[26]。 基于XML（eXtensible Markup Language）,构建了海洋天然产物数据库，目前数据库中包含了大约11000种海洋天然化合物的分子结构及其生物来源、生物活性、药效、物理化学参数和参考文献等信息。与XML的集成有利于该数据库的管理、查询、Web发布、数据交换以及未来的扩展，目前已具备Web检索的功能[27, 28]。 提出了一种将遗传算法与随机搜寻相结合的分子对接算法，并采取分步考虑形状互补和能量最优的策略。对随机选取的8个复合物结构测试结果表明：计算所得构象及取向与晶体结构接近，大部分情况下，最小的均方根偏差小于[29]。 5,6-二芳基取代吡里酮类化合物是一类新型的抗炎镇痛剂，利用分子对接方法研究了这一类化合物的作用机理。从分子对接结果来看，ZZ-122与COX-2结合的比较好，静电作用和范德华作用都是负值，范德华作用是主要作用；而与COX-1结合构象严重扭曲，两者之间存在强烈的范德华斥力。这一结果表明，该类化合物可能是比较好的COX-2选择性抑制剂[30]。