北京大学化学信息学course-13三维结构的实验数据资源基于三维结构的计算方法课件

化学信息学第十二讲 三维结构的实验数据资源 基于三维结构的计算方法1剑桥晶体结构数据库（CSD）有机物、有机金属和金属复合物二维结构/三维结构晶体学数据（空间群、晶胞参数）、X-ray/neutron衍射数据等参考文献2蛋白质数据库:PDB大于10,000个大分子三维结构包括蛋白质、碳水化合物、DNA、DNA/蛋白复合物X-ray/Neutron衍射、NMR、结构计算3其它非实验数据计算得到的三维结构数据库专业数据库天然产物DNPNAPALERT代谢库Metabolite毒性库Toxicity药物分子库MDDR-3DWDIACD-3D数据库查询结果浏览界面5计算机辅助药物设计的一般过程ACDCSDMDDRNIH中草药数据库(CHDD)构效关系模型计算组合化学优化理想活性化合物a.原子生长法b.片段生长法c.模片裁剪法d.分子对接法数据库搜索a.构象搜索b.构象迭合c.药效团识别

药效团QSAR先导化合物合成测活具有相同作用机制的系列化合物受体数据库(PDB等)受体模型.复合物结构c.分子对接b.残基突变d.共振识别活性位点或结合位点计算机辅助药物设计的一般过程6基于三维结构的计算方法7定量构效关系：QSAR简介1867年，Crum-Brow和Fraser提出构效关系的概念1900年，Overton和Meyer等提出了麻醉作用的类脂学说，即化学结构各异的麻醉剂其活性随着脂－水分配系数的增加而增加，这可能是最早提出的化合物生理活性和物理化学性质之间的定量关系模型。Hansch法和Free-Wilson法的提出，多变量解析技术的引入使QSAR方法逐渐发展和应用起来。QSAR/QSPR9定量构效关系：QSAR2D-QSARHansch法:从盲目药物设计过渡到合理药物设计（RationalDrugDesign）的重要标志定量Hansch法的原理假定今采用合适的方法投药于生物体内，经过转运过程，药物（D）和靶点（T）之间产生相互作用。这个相互作用的具体过程我们可能并不清楚，但可假定形成的产物为（P），从而引起所期望的药物活性。10定量构效关系：QSAR对同一药物的不同给药浓度D0为给药部位浓度，logA表示与转运过程相关的量，G为药物–靶点的结合自由能对不同药物的相同给药浓度相当于所测定的效应值的对数11定量构效关系：QSAR膜通透的双直线模型logPlogA抛物线模型双线性模型13定量构效关系：QSARHansch-Fujita公式Es是与立体效应有关的取代基参数。式中的a、b、、和C通过最小二乘法得到14定量构效关系：QSARFree-Wilson法假定分子的活性是由母体化合物和取代基的活性贡献之和不论其它位置的取代基是否变化，每一取代基对生物活性的贡献是恒定的而且可加和的如果对于同一母体构成的同系列化合物，当取代基位置为R个，引入各个位置上的取代基种类为Kj个时化合物的生物活性y是母体结构对活性的贡献加上取代基对活性的贡献之和

15定量构效关系：QSAR3D-QSAR基本假设生物活性是有待考察的化合物产生的而不是其代谢产物使用的构象是具有生物活性的构想所有参与计算的化合物的结合位点是相同的生物活性主要是由焓变解释熵效应对所有化合物来说都相似只考虑平衡体系，不考虑反应的快慢药理动力学中的溶剂效应、扩散、传输等过程都不考虑17定量构效关系：QSAR分子场一个分子的场效应可以使用格点的方法表示每个单元代表目标分子与其相互作用的分子之间的吸引或排斥作用与目标分子相互作用的分子可以是水、醇及其它溶剂18定量构效关系：QSAR常见分子场MEP–静电场，（单位正电荷作探针）MLP–疏水场，（不需要探针）GRID–相互作用的总能量（包括空间效应(Lennard-Jones),氢键和静电势(可以使用任意探针)CoMFA–标准的场包括:空间效应和静电势，此外，也可以考虑氢键等的作用19定量构效关系：QSARCOMFA使用基本步骤：确定药物分子的活性构象，再按照一定的规则（一般为骨架叠合或场叠合）进行药物分子的叠合。然后，在叠合好的分子周围定义分子场的空间范围，可以采用矩形或圆形的空间，把定义的空间按照一定的步长均匀划分产生格点

在每个格点上用一个探针离子来评价格点上的分子场特征（一般为静电场和立体场，有时也包括疏水场和氢键场），在CoMFA分子中可以采用不同的分子场势能函数

通过偏最小二乘（PLS）方法建立化合物活性和分子场特征之间的关系21223D-QSAR/COMFA应用示例对15个CYP450复合物使用COMFA方法建模静电场：

BLUE–带正电

RED–带负电立体场:

GREEN–空间填充匹配较好区域

YELLOW–空间冲突区域23定量构效关系：QSAR现代建模方法输入:n描述符P1,..Pn和m各化合物的生物活性实验值(EC50,LD50,IC50等)BioP1P2...PnCpd10.73.7Cpd23.20.4…Cpdm25定量构效关系：QSAR多元统计回归（MLR）偏最小二乘（PLS）神经网络遗传算法26药效团模型

d1=d2=d3=11.01.090105.00.5d1=4.51.0d2=6.71.0d3=7.61.07.40.81.0–1.74.30.83.81.00.9–1.088-14829药效团模型与QSAR方法的比较相同：都是以配体效分子作为研究的起点不同：QSAR:一般基于母体相似的同系列化合物，它得到的是化合物结构和活性之间的定量关系，这个定量关系一般只能要来指导先导化合物的改造；药效团模型可以基于不同类的先导化合物，它得到的是与生物活性有关的重要的药效团特征，只要分子拥有这些特征，就可能具有某种生物活性，而这些活性配体分子的结构未必需要相同，因此药效团模型方法可以用来寻找结构全新的先导化合物

30药效团模型基本步骤药效团模型的识别：药效团模型的识别仅仅是得到药效团模型，通过药效团模型，我们希望能够找到新的先导化合物基于药效团模型的数据库搜索：通过数据库搜索，来寻找具有特定药效团特征的化合物，

31药效团模型药效团模型的识别在受体结构未知或者在作用机制不明确的情况下，对一系列化合物进行结构－活性研究，并结合构象分析，分子叠合等手段，得到一些对化合物活性至关重要的原子、基团或化学功能结构，以及其空间关系。在受体结构已知的情况下，分析受体的作用位点以及药物分子和受体之间的相互作用模式，根据预测的复合物结构或相互作用信息来推知可能的药效团结构。这里我们主要介绍第一种方法，即怎样仅仅根据配体结构得到可能的药效团模型。32药效团模型药效团模型识别基本步骤活性化合物的选择以及药效特征元素的定义构象分析。对每个化合物都要进行构象分析，得到在某一能量范围内的构象。这个能量范围一般是根据经验来指定，分子叠合和药效团映射（pharmacophoremapping）。得到药效特征元素以及构象以后，就需要对这些构象进行叠合以找到共同的药效团模型。药效团模型的修正。分子叠合得到的药效团模型往往不是最优的，可能还需要根据自己的经验和其它实验或计算的结果对药效团模型做进一步的修正，比如加入排除体积。33药效团模型药效团模型数据库搜索的基本步骤初筛：把在药效团中出现频率较多的原子或基团作为约束条件，每个约束条件作为一个编码，如果在提问结构中出现的编码也出现在这些在这些被查找的分子中，则这个分子可以通过初筛，否则被淘汰。二维子结构匹配：这步的目的是确定测试分子中的药效模式间的连接方式是否和提问结构相同。三维结构搜索：验证它在数据库中的构象和药效团中的药效特征元素是否满足空间的限制条件，如果满足，则该分子是一个命中结构。34药效团模型药效团模型搜索常用软件DISCOGASPCATALYST35MuscarinicM3受体拮抗剂的设计MuscarinicM3受体是一种非常重要的药物设计靶点，它与过敏性肠道综合症（IBS），慢性阻塞性气管疾病（COAD）以及小便失禁有关1999年，Marriott等采用基于药效团的数据库搜索的方法进行了MuscarinicM3受体拮抗剂的设计

36MuscarinicM3受体拮抗剂的设计

用来搭建药效团模型的MuscarinicM3受体拮抗剂的结构37MuscarinicM3受体拮抗剂的设计

DISCO产生的两个用于数据库搜索的药效团模型

38MuscarinicM3受体拮抗剂的设计

数据库搜索得到的三个活性化合物

39MuscarinicM3受体拮抗剂的设计