计算机辅助药物分子设计

计算机辅助药物分子设计第1页/共58页药物：是人类用来预防、治疗、诊断疾病，或为了调节人体功能、提高生活质量、保持身体健康的化学品。新药需求：对生活水平要求提高；对新疾病（HIV、SARS、H1N1）的诊断与治

疗；耐药性。药物的设计与研发：成为一门新的学科。第2页/共58页我国制药业目前状况：以仿制

为主，自己创制的新药仅占2％-3％。药物研发的时间和费用：一种新药，花费10~12年，耗资10亿美元。计算机辅助药物设计：发展新的理论计算方法，从传统偶然发现到计算机辅助设计。第3页/共58页第4页/共58页新药研究的三个重要阶段第5页/共58页先导(导向)化合物(Lead

compound)的发现;先导化合物的优化研究;临床与开发研究。第6页/共58页第7页/共58页新药物研究的主要途径第8页/共58页天然产物的获得(先导药物主要来源);现有药物的开发(包括部分结构改造);生理学机理研究(包括代谢研究);90年代诞生了组合化学(CombinatorialChemistry)；偶然发现与随机筛选药物合理设计与高通量筛选(MassScreening

or

High

Throughput

Screening)第9页/共58页青霉素的发现第10页/共58页1928年的一天，弗莱明和往常一样，一到实验室，便观察培养皿里的葡萄球菌的生长情况。突然，他发现一个培养皿里生长了一团青绿色的霉。他注意到，青霉周围呈现出一片清澈。弗莱明立刻意识到，这是葡萄球菌被杀死的迹象。为了证实自己的判断，弗莱明用吸管从培养皿中吸取了一滴溶液，涂在干净的玻璃上，然后放在高倍显微镜下观察。结果，溶液里一个葡萄球菌也没有。发现，青霉所分泌的物质，对白喉菌、炭疽菌等，都有强效的杀菌效果。第11页/共58页计算机辅助药物设计第12页/共58页传统药物设计具有很大的盲目性，一般平均要筛选10000种化合物以上才能得到一种新药，因此开发效率很低。随着计算机技术及计算化学、分子生物学和药物化学的发展，药物设计进入了理性阶段，它是依据生物化学、酶学、分子生物学以及遗传学等生命科学的研究成果，针对潜在的药物设计靶点，并参考其它类源性配体或天然产物的化学结构特征，利用计算化学和分子模拟方法设计出合理的药物分子。第13页/共58页第14页/共58页第15页/共58页第16页/共58页第17页/共58页白血病青光眼第18页/共58页第19页/共58页药物设计目标：设计一个与受体结合能力强的化合物第20页/共58页依据酶催化反应物的结构设计第21页共58页计算受体与配体相互作用的方法第22页/共58页量子化学法QM分子力学MM量子化学与分子力学结合方法QM/MM分子力学方法第23页/共58页分子力学优缺点第24页/共58页优点：计算速度快概念清楚，便于理解及应用缺点：不能处理电荷的极化现象不能很好处理金属原子不能处理化学键的断裂第25页/共58页第26页/共58页0.5~1.0第27页/共58页第28页/共58页第29页/共58页第30页/共58页第31页/共58页第32页/共58页第33页/共58页第34页/共58页第35页/共58页第36页/共58页第37页/共58页类药5规则第38页/共58页类药五原则，是辉瑞公司资深药物化学家Christopher

A.Lipinski在1997年提出的筛选类药分子的五条基本法则，符合Lipinski规则的化合物会有更好的药代动力学性质，在生物体内代谢过程中会有更高的生物利用度，因而也更有可能成为口服药物。类药五原则（rule

of

five）也称为Lipinski规则，其内容如下：一个小分子药物中要具备以下性质分子量小于500；氢键给体数目小于5；氢键受体数目小于10；脂水分配系数小于5；可旋转键的数量不超过10个。类药5规则第39页/共58页在长期的实践过程中，药物化学家们对

Lipinski规则作出简化，形成“四规则”和

“三规则”，但是四规则和三规则有时仍

然被称作“五规则”，这里的五指的是各

条规则的判别值均为5或500。简化后的四

规则去掉了关于可旋转键的数量限制；三规则进一步去掉了对氢键受体数量的限制。这是根据口服药物总结出的经验型规律，是小分子药物设计的有效指导原则。脂水分配系数脂水分配系数（lipo-hydro

partition

coefficient）为化合物在脂相和水相间达到平衡时的浓度比值，通常是以化合物在有机相中的浓度为分子，在水相中的浓度为分母。脂水分配系数越大，越易溶于脂，反之则越易溶于水。易溶于脂的物质在机体内呈现亲脂性或疏水性，而易溶于水的现象称为亲水性脂水分配系数以P表示，化合物在有机相中浓度Co，在水相的浓度Cw第40页/共58页脂水分配系数(cont.)Since

the

differences

areusually

on

a

very

large

scale,Log10(P)

is

used.OH1-octanolOH第41页/共58页Hwater第42页/共58页基于密度泛函理论与蛋白模型的流感病毒抑制剂的设计第43页/共58页Neu

minidases第44页/共58页cleave

host

receptorshelp

release

of

new

virions流感病毒神经氨酸酶与抑制剂的复合物第45页/共58页第46页/共58页基于密度泛函与蛋白模型的方法第47页/共58页传统力场计算方法：缺点：精度低，不具有普遍性；优点：速度快。传统量子力学方法：缺点：速度慢，体系小，没有被广泛应用；优点：精

度高；基于蛋白模型与密度泛函的方法：简化计算体系，提高计算精度；蛋白模型方法的应用－流感病毒抑制剂的设计第48页/共58页

流感病毒的危害性：传染性强；威胁生命安全，例如H5N1、H1N1。

神经氨酸酶流感药物：扎那米韦

(Zanamivir)、奥塞米韦(Oseltamivir)。

新型流感病毒抑制剂的必要性：抗药性的新变异病毒的出现；提高药效。蛋白质三维晶体结构（PDPBDBc:o1dFeK:G1）FKG）受体蛋白质模型（~400个原子）DFT方法优化各项构型（蛋白质模型，配体，蛋白质-配体复合体）计算各项的能量受体蛋白与配体的相互作用能判断配体与受体的结合能力强弱配体设计研究方法第49页/共58页FKBP12与抑制剂的复合物the

binding

pocket

of

FKBP12第50页/共58页FKBP12与抑制剂的复合物

DFT结果第51页/共58页FKBP12与抑制剂的复合物

Autodock结果第52页/共58页第53页/共58页基于奥塞米韦(Oseltamivir)的抑制剂设计第54页/共58页已设计了5个化合物，其活性约为奥塞米韦的100倍。磺酰脲除草剂氯嘧磺隆（CE)与拟南芥乙酰乳酸合成酶复合物的晶体结构第55页/共58页配体分子CE结论：量子力学与模型蛋白相结合的方法可以作为除草剂设计的一种可靠方