《比较分子力场分析法研究苯基烷胺类致幻剂的定量构效关系》由会员分享，可在线阅读，更多相关《比较分子力场分析法研究苯基烷胺类致幻剂的定量构效关系（9页珍藏版）》请在金锄头文库上搜索。
1、1比较分子力场分析法研究苯基烷胺类致幻剂的定量构效关系【摘要】 本文运用比较分子力场分析（CoMFA）方法对 89 个苯基脘胺类致幻剂进行了三维定量构效关系研究。考察了不同的步长、探针原子以及电荷的计算方法等对构建模型的影响。对于研究苯烷基胺类致幻剂的作用机理以及指导发现新的致幻剂均具有重要意义。【关键词】 致幻剂； 苯基烷胺； 三维定量构效关系； 比较分子力场分析法CoMFA applied to QSAR study of phenylalkylamine hallucinogens【Abstract】 The three dimensional structure-activity re
2、lationship has been studied on a series of phenylalkylamine hallucinogens by comparative molecular field analysis (CoMFA). The contributions from steric and electrostatic fields were 0.536 and 0.464 respectively.It is very helpful to study the mechanism and find new compounds of such kind 2of halluc
3、inogens. 【Key words】 hallucinogens； phenylalkylamine； 3D-QSAR； CoMFA致幻剂（hallucinogens）又名拟精神病类药物，这类药物能严重改变知觉、情绪及主体认知过程,引起幻觉或精神病样症状。这类化合物在某些国家和地区当作毒品的成分使用。特别是现在许多青少年为寻找刺激，大量滥用这类药物，严重阻碍了他们的健康成长。致幻剂很有可能成为恐怖分子利用的武器，因此，为了维护国际社会的安定，保障人的健康发展，加强对这类化合物的严格管理，快速注册具有依赖性及潜在致幻性的药物，建立这类物质的定量构效关系模型及设计新的化合物都是很必要的。按照化
4、学结构，致幻剂主要分为两大类，一类是苯基烷胺类（包括苯乙胺及苯异丙胺衍生物等） ，另一类是吲哚基烷胺类（包括色胺及麦角灵 LSD 及其衍生物等） 。为了在分子水平上理解致幻剂性物质的活性，人们曾建立了一些 QSAR 模型14 ，这些模型多是基于量子化学参数建立的，而量化计算一般速度慢、耗时长。3本文利用三维定量构效关系研究中常用的比较分子力场分析法（CoMFA）研究了 89 个苯基烷胺类化合物的结构活性关系，建立了三维定量构效关系模型。该法较量化计算速度快、简便、直观，而且它充分考虑分子的三维结构和药物与受体的结合模式，建立的模型有很好的预测能力。1 研究方法1.1 CoMFA 的原理 比较分
5、子力场分析（comparative molecular field analysis, CoMFA）是目前应用最广泛的一种 3D-QSAR 研究方法。其基本假设是：所研究的化合物作用于同一受体作用位点，在分子水平上，影响生物活性的相互作用通常是非共价性的，分子力场侧重处理立体和静电等非共价性相互作用，它通过计算分子周围一定范围内三维空间中假想网格上探针与分子之间的立体场能和静电场能，运用偏最小二乘法（PLS）寻找最佳主成分数建立预测模型。并将结果以直观的三维图形方式表示出来。由于CoMFA 充分考虑了分子的三维结构和药物与受体的结合模式，在合理药物设计工作方面取得了较大的成功58 。41.2
6、化合物活性数据 本文研究的 89 个苯基烷胺类致幻剂活性数据与结构都来自文献3 ，活性数据是由 A.T.Shugin 小组通过人体试验获得的。数据以阈剂量与最大剂量的平均值表示，并以Mescaline 做参考，得到的是各个化合物的相对效价，在建立模型时采用其对数形式（在表中以 logA 表示） 。随机选取其中的 15 个化合物作为预测集，其他 74 个化合物为训练集。1.3 CoMFA 分析 采用 sybyl7.1 中的 QSAR 模块中的CoMFA 方法，将所有化合物放置于一个三维网格中，计算过程中网格区域应包含所有分子,并且在 X、Y、Z 轴方向至少延伸0.4nm。选择合适的步长及探针原子
7、计算立体场与静电场能，立体场一般采用 Lennard Jones 势的计算公式 ,静电场采用 Coulomb势，能量截断值(cutoff) 设为 126kJ/mol。利用偏最小二乘法 (PLS)建立模型，首先做 LOO(Leave-one-out)交叉验证，得到最佳主成分数。然后在最佳主成分数下做非交叉验证。为检验模型的预测能力，对预测组化合物做外部检验，由下面的公式得到 R2pre，最后将结果以直观的三维等高图显示。R2pre=SD-PRESS SD 其中 PRESS 是预测组化合物活性预测值与实验值残差平方和，SD 是预测组化合物活预测值与训练组化合物活性平均值的残差平方和。52 结果与讨
8、论2.1 空间网格及步长的选择 一般不同的网格设置对结果会产生一定的影响，采用系统自动生成的网格效果较差（q20.489，R20.851） ，自动生成的网格范围为 X -1.140.92，Y -0.901.20，Z -0.850.81（单位:nm）在此基础上，以 0.4nm 为单位调整网格大小并建模，最终确定网格为 X -1.000.80，Y -0.900.90，Z -0.800.80（单位:nm） ，并在此网格范围内研究了步长取值对结果的影响。结果见表 1。表 1 不同的步长值对结果的影响略2.2 探针的选择 探针用来模拟化合物分子与受体之间的非共价性相互作用，一般常用的是 sp3 杂化的
9、C,电荷为+1。实验表明不同的探针原子会对 CoMFA 模型产生一定的影响。本文尝试用不同的探针建立模型，结果表明，以 sp3 杂化的 O 作为探针得到的模型 q2 最大，而且非交叉验证时误差最小，因此我们最终选用了 sp3杂化的 O 电荷为+1 作为探针建立了模型。62.3 不同的电荷计算方法 采用不同的电荷计算方法在一定程度上也会对结果产生影响，本文比较了MMFF94、Gasteiger、Gasteiger-Hckle 、Hckle 等四种方法对计算结果的影响。从结果可以看出采用 MMFF94 和 Gasteiger-Hckle 电荷比用 Gasteiger 或 Hckle 所得模型效果要
10、好，这可能是因为后两种方法只考虑 或 电子，而这些化合物与受体发生作用时两种电子都会产生重要作用。前两种方法相比，MMFF94 电荷计算得到的模型标准偏差较大。因此，我们选用了 Gasteiger-Hckle 作为最佳的电荷计算方法。2.4 CoMFA 结果 图 1 是训练集和预测集中苯烷基胺类致幻剂的实验值与测定值的相关关系图。从图 1 中可以看出，实验值与预测值均匀分布在直线 Y=X 两侧，这表明实验值与预测值比较接近。与参考文献 5 的结果相比，本文所建立的模型得到的最大残差为0.422,标准偏差为 0.105，而文献 5 中分别为 0.73 和 0.192。相比之下，本文利用 CoMF
11、A 方法所建立的模型预测能力更好，结果更可靠。7图 1 致幻剂的活性实验值与预测值相关图略3 结论本文采用 CoMFA 方法对 89 个苯烷基胺类致幻剂进行了构效关系研究。考察了不同步长、探针及电荷计算方法对构建模型的影响, 最后得到了具有较强预测能力的模型,并对随机选取的预测组化合物活性进行了预测，得到了满意的结果。【参考文献】1 Gupta S P, Singh P, Bindal M C.Chem Rev，1983, 83(6):633 649.2 Clare B W.J Med Chem，1990,33 (2): 687702.83 Clare B W.Computer-Aided M
12、ol Des， 2002,16(8/9): 611 633.4 HU Wen-Xiang(胡文祥),WANG Jian-Ying(王建营). Synergetic Combinatorial Chemistry（协同组合化学） .Beijing: Science Press, 2003， 54.5 Kulkarni S S, Grundt P. J Med Chem， 2004, 47(13): 3388-3398.6 GU Yan(谷妍),DONG Xi-Cheng(董喜城),CHEN Hai-Feng(陈海峰)，et al. Acta Chimica Sinica(化学学报),2000,58(12):15401545.7 MA Yi (马翼),JIANG Lin (姜林), LI Zheng-M ing (李正名)，et al. Chem J Chinese Universities (高等学校化学学报),2004,25(11):20312033.98 Cruz V L, Ramos J. Organometallics, 2005, 24(21), 5095-5102.
