1、1化学信息学第十二讲第十二讲三维结构的实验数据资源三维结构的实验数据资源基于三维结构的计算方法基于三维结构的计算方法2剑桥晶体结构数据库（CSD）n有机物、有机金属和金属复合物n二维结构/三维结构n晶体学数据（空间群、晶胞参数）、X-ray/neutron 衍射数据等n参考文献3蛋白质数据库:PDBnHttp:/www.rcsb.org/pdbn大于10,000个大分子三维结构n包括蛋白质、碳水化合物、DNA、DNA/蛋白复合物nX-ray/Neutron衍射、NMR、结构计算4其它非实验数据计算得到的三维结构数据库w 商用化合物库nCAnBeilstennCMCnACD-3Dw 开放数据库和
2、专有数据库nNCI5其它非实验数据计算得到的三维结构数据库w 专业数据库n天然产物lDNPlNAPALERTn代谢库lMetabolite n毒性库lToxicity n药物分子库lMDDR-3DlWDIACD-3D数据库查询结果浏览界面6计算机辅助药物设计的一般过程计算机辅助药物设计的一般过程ACDCSDMDDRNIH中草药数据库(CHDD)构效关系模型计算组合化学优化 理想活性化合物a.原子生长法b.片段生长法c.模片裁剪法d.分子对接法数据库搜索a.构象搜索b.构象迭合c.药效团识别 药效团 QSAR先导化合物 合成 测活具有相同作用机制的系列化合物受体数据库(PDB等)受体模型.复合物
3、结构c.分子对接b.残基突变d.共振识别活性位点或结合位点计算机辅助药物设计的一般过程7基于三维结构的计算方法基于三维结构的计算方法8定量构效关系：QSARw 简介w 2D-QSARw 3D-QSARw 结构描述符分类w 建模方法9定量构效关系：QSARw 简介n1867年，Crum-Brow和Fraser提出构效关系的概念n1900年，Overton和Meyer等提出了麻醉作用的类脂学说，即化学结构各异的麻醉剂其活性随着脂水分配系数的增加而增加，这可能是最早提出的化合物生理活性和物理化学性质之间的定量关系模型。nHansch法和Free-Wilson法的提出，n多变量解析技术的引入使QSAR
4、方法逐渐发展和应用起来。nQSAR/QSPR10定量构效关系：QSARw 2D-QSARnHansch法:从盲目药物设计过渡到合理药物设计（Rational Drug Design）的重要标志 n定量Hansch法的原理l假定今采用合适的方法投药于生物体内，经过转运过程，药物（D）和靶点（T）之间产生相互作用。这个相互作用的具体过程我们可能并不清楚，但可假定形成的产物为（P），从而引起所期望的药物活性。PTDTDkK.11定量构效关系：QSARw 对同一药物的不同给药浓度nD0为给药部位浓度，logA表示与转运过程相关的量，G为药物靶点的结合自由能 w 对不同药物的相同给药浓度n相当于所测定的
5、效应值的对数 10434.0log1logCRTGAD2*434.0log)log(CRTGAP12定量构效关系：QSARw 抛物线模型n在多种情况下，药物在体内的扩散，主要是被动的转运过程。在此条件下，药物要交替穿过水相和类脂构成的体系，其移动难易程度和logP呈现函数关系，则经过一定时间后药物在最末一相中为浓度logC。CPbPaA)(log)(loglog213定量构效关系：QSARw 膜通透的双直线模型logPlogA抛物线模型双线性模型CPbPaA)1log(loglog14定量构效关系：QSARw Hansch-Fujita公式nEs是与立体效应有关的取代基参数。式中的a、b、和C
6、通过最小二乘法得到 CEsbaPD2*0log或1log15定量构效关系：QSARw Free-Wilson 法n假定分子的活性是由母体化合物和取代基的活性贡献之和n不论其它位置的取代基是否变化，每一取代基对生物活性的贡献是恒定的而且可加和的n如果对于同一母体构成的同系列化合物，当取代基位置为R个，引入各个位置上的取代基种类为Kj个时 n化合物的生物活性y是母体结构对活性的贡献加上取代基对活性的贡献之和 R1jK1kjkjkiay16定量构效关系：QSARw 3D-QSARn距离几何（distance geometry,DG）n分子形状分析（molecular shape analysis,M
7、SA）n比较分子场分析（comparative molecular field analysis,CoMFA）n虚拟受体（phesudo receptor）17定量构效关系：QSARw 3D-QSAR基本假设n生物活性是有待考察的化合物产生的而不是其代谢产物n使用的构象是具有生物活性的构想n所有参与计算的化合物的结合位点是相同的n生物活性主要是由焓变解释n熵效应对所有化合物来说都相似n只考虑平衡体系，不考虑反应的快慢n药理动力学中的溶剂效应、扩散、传输等过程都不考虑18定量构效关系：QSARw 分子场n一个分子的场效应可以使用格点的方法表示n每个单元代表目标分子与其相互作用的分子之间的吸引或排
8、斥作用n与目标分子相互作用的分子可以是水、醇及其它溶剂19定量构效关系：QSARw 常见分子场nMEP 静电场，（单位正电荷作探针）nMLP 疏水场，（不需要探针）nGRID 相互作用的总能量（包括空间效应(Lennard-Jones),氢键和静电势(可以使用任意探针)nCoMFA 标准的场包括:空间效应和静电势，此外，也可以考虑氢键等的作用20定量构效关系：QSARw COMFA基本原理：n能量场的计算是通过将探针原子放在每个单元进行计算产生的n分子场的标准类型包括l空间效应(Lennard-Jones)l静电势作用(库仑作用)l探针原子带有+1.0电荷的sp3杂化的碳原子21定量构效关系：
9、QSARwCOMFA使用基本步骤：1.确定药物分子的活性构象，再按照一定的规则（一般为骨架叠合或场叠合）进行药物分子的叠合。然后，在叠合好的分子周围定义分子场的空间范围，可以采用矩形或圆形的空间，把定义的空间按照一定的步长均匀划分产生格点 2.在每个格点上用一个探针离子来评价格点上的分子场特征（一般为静电场和立体场，有时也包括疏水场和氢键场），在CoMFA分子中可以采用不同的分子场势能函数 3.通过偏最小二乘（PLS）方法建立化合物活性和分子场特征之间的关系 22233D-QSAR/COMFA应用示例w 对15个CYP450复合物使用COMFA方法建模静电场：BLUEBLUE 带正电REDRE
10、D 带负电 立体场:GREENGREEN 空间填充匹配较好区域 YELLOWYELLOW 空间冲突区域24定量构效关系：QSARw 结构描述符的分类n疏水参数：n电子参数：n几何参数：n拓扑参数：n组成参数：n立体参数：n量化计算参数：25定量构效关系：QSARw 现代建模方法n输入:n 描述符P1,.Pn 和m各化合物的生物活性实验值(EC50,LD50,IC50等)BioP1P2.PnCpd10.73.7Cpd23.20.4Cpdm26定量构效关系：QSARn多元统计回归（MLR）n偏最小二乘（PLS）n神经网络n遗传算法27药效团模型w 简介n药效团的概念在1909年由P.Ehrilic
11、h提出，指载有活性必须的特征原子的分子框架。n现在泛指在生物活性分子（一般指配体小分子）中对活性起重要作用的“药效特征元素药效特征元素”的空间排列形式空间排列形式。这些“药效特征元素”可以认为是配体和受体发生相互作用时的活性部位。n这些“药效特征元素”可以是l某些具体的原子或原子团，比如氧原子、羟基、苯环等等l某些特定的化学功能结构，如疏水团、氢键给体、氢键受体等等28药效团模型w 简介n药效团的概念在1909年由P.Ehrilich提出，指载有活性必须的特征原子的分子框架。n现在泛指在生物活性分子（一般指配体小分子）中对活性起重要作用的“药效特征元素药效特征元素”的空间排列形式空间排列形式。
12、这些“药效特征元素”可以认为是配体和受体发生相互作用时的活性部位。n这些“药效特征元素”可以是l某些具体的原子或原子团，比如氧原子、羟基、苯环等等l某些特定的化学功能结构，如疏水团、氢键给体、氢键受体等等29药效团模型 OOHOHOHOCOHd1d2d3 NN+d1d2d3hydrophobic centerZn的配体氢键受体孤对电子方向阴离子阴离子中心原子d1=d2=d3=11.0 1.090 105.0 0.5d1=4.5 1.0d2=6.7 1.0d3=7.6 1.07.4 0.81.0 1.74.3 0.83.8 1.00.9 1.088-14830药效团模型w 与QSAR方法的比较n
13、相同：都是以配体效分子作为研究的起点n不同：lQSAR:一般基于母体相似的同系列化合物母体相似的同系列化合物，它得到的是化合物结构和活性之间的定量关系，这个定量关系一般只能要来指导先导化合物的改造指导先导化合物的改造；l药效团模型可以基于不同类的先导化合物基于不同类的先导化合物，它得到的是与生物活性有关的重要的药效团特征，只要分子拥有这些特征，就可能具有某种生物活性，而这些活性配体分子的结构未必需要相同，因此药效团模型方法可以用来寻找结构全新的先寻找结构全新的先导化合物导化合物 31药效团模型w基本步骤1.药效团模型的识别：药效团模型的识别仅仅是得到药效团模型，通过药效团模型，我们希望能够找到
14、新的先导化合物2.基于药效团模型的数据库搜索：通过数据库搜索，来寻找具有特定药效团特征的化合物，32药效团模型w 药效团模型的识别n在受体结构未知或者在作用机制不明确受体结构未知或者在作用机制不明确的情况下，对一系列化合物进行结构活性研究，并结合构象分析，分子叠合等手段，得到一些对化合物活性至关重要的原子、基团或化学功能结构，以及其空间关系。n在受体结构已知受体结构已知的情况下，分析受体的作用位点以及药物分子和受体之间的相互作用模式，根据预测的复合物结构或相互作用信息来推知可能的药效团结构。这里我们主要介绍第一种方法，即怎样仅仅根据配体结构得到可能的药效团模型。33药效团模型w药效团模型识别基
15、本步骤1.1.活性化合物的选择以及药效特征元素的定义活性化合物的选择以及药效特征元素的定义2.2.构象分析构象分析。对每个化合物都要进行构象分析，得到在某一能量范围内的构象。这个能量范围一般是根据经验来指定，3.3.分子叠合和药效团映射分子叠合和药效团映射（pharmacophore mapping）。得到药效特征元素以及构象以后，就需要对这些构象进行叠合以找到共同的药效团模型。4.4.药效团模型的修正药效团模型的修正。分子叠合得到的药效团模型往往不是最优的，可能还需要根据自己的经验和其它实验或计算的结果对药效团模型做进一步的修正，比如加入排除体积。34药效团模型w药效团模型数据库搜索的基本步
16、骤1.1.初筛初筛：把在药效团中出现频率较多的原子或基团作为约束条件，每个约束条件作为一个编码，如果在提问结构中出现的编码也出现在这些在这些被查找的分子中，则这个分子可以通过初筛，否则被淘汰。2.2.二维子结构匹配二维子结构匹配：这步的目的是确定测试分子中的药效模式间的连接方式是否和提问结构相同。3.3.三维结构搜索三维结构搜索：验证它在数据库中的构象和药效团中的药效特征元素是否满足空间的限制条件，如果满足，则该分子是一个命中结构。35药效团模型w 药效团模型搜索常用软件nDISCOnGASPnCATALYST36Muscarinic M3受体拮抗剂的设计wMuscarinic M3受体是一种
17、非常重要的药物设计靶点，它与过敏性肠道综合症（IBS），慢性阻塞性气管疾病（COAD）以及小便失禁有关w1999年，Marriott等采用基于药效团的数据库搜索的方法进行了Muscarinic M3受体拮抗剂的设计 37Muscarinic M3受体拮抗剂的设计 OHOONNN NOO OON 用来搭建药效团模型的Muscarinic M3受体拮抗剂的结构38Muscarinic M3受体拮抗剂的设计 DISCO产生的两个用于数据库搜索的药效团模型 N+HdonHaccHdon5.85.46.52.95.07.6N+HdonHaccHdon4.15.48.52.96.06.6N5.85.46.
18、52.95.07.6Hacceptor atomHdonor siteHdonor siteOONdonor sitedonor site39Muscarinic M3受体拮抗剂的设计 数据库搜索得到的三个活性化合物 NHNOO ONHNONH OON40Muscarinic M3受体拮抗剂的设计w 采用第一个药效团模型，数据库搜索给出了176个命中结构w 采用第二个药效团模型，数据库搜索给出了173个命中结构。w 在得到的这两组化合物中，有172个化合物是完全相同的。两次搜索结构包含177个完全不同的分子。w 在总共这177个命中结构中，对172个分子进行了生物活性测定，其中有三个分子显示了明显的活性，其中第三个化合物的活性最高，其pA2为6.670.07。
