中国药科大学cadd重点概括.docx
《中国药科大学cadd重点概括.docx》由会员分享,可在线阅读,更多相关《中国药科大学cadd重点概括.docx(13页珍藏版)》请在冰豆网上搜索。
![中国药科大学cadd重点概括.docx](https://file1.bdocx.com/fileroot1/2022-10/11/f374b373-4eee-43e2-8ba5-f484f95f2601/f374b373-4eee-43e2-8ba5-f484f95f26011.gif)
中国药科大学cadd重点概括
第1章概论
一、药物发现一般过程
新药的研究有三个决定阶段:
先导化合物的发现,新药物的优化研究,临床与开发研究。
计算机辅助药物设计的主要任务就是先导化合物的发现与优化。
二、合理药物设计
1、合理药物设计(rationaldrugdesign)是依据与药物作用的靶点,即广义上的受体,如酶、受体、离子通道、病毒、核酸、多糖等,寻找和设计合理的药物分子。
通过对药物和受体的结构在分子水平甚至电子水平的全面准确了解进行基于结构的药物设计和通过对靶点的结构、功能、与药物作用方式及产生生理活性的机理的认识基于机理的药物设计。
CADD通过内源性物质或外源性小分子作为效应子作用于机体的靶点,考察其形状互补,性质互补(包括氢键、疏水性、静电等),溶剂效应及运动协调性等进行分子设计。
2、方法分类
(1)合理药物设计有基于靶点结构的三维结构搜索和全新药物设计等方法。
后者分为模板定位法、原子生长法、分子碎片法(碎片连接法和碎片生长法)。
(2)根据受体是否已知分为直接药物设计和间接药物设计。
前者即通过结构测定已知受体或受体-配体复合物的三维结构,根据受体的三维结构要求设计新药的结构。
受体结构测定方法:
同源模建(知道氨基酸序列不知道空间结构时),X射线衍射(可结晶并得到晶体时),多维核磁共振技术(溶液状态)。
后者通过一些配体的结构知识(SAR,计算机图形显示等)推测受体的图像,提出家乡受体,采用建立Pharmacophore模型或3D-QSAR和基于药效团模型的三维结构搜索等方法,间接进行药物设计。
三、计算化学
计算化学包括分子模型、计算方法、计算机辅助分子设计(CAMD)、化学数据库及有机合成设计。
计算方法包括很多种,但基本上可以分为两大类:
分子力学和量子力学(分为从头计算方法和半经验方法)。
常用的计算应用有:
(1)单点能计算:
根据模型中原子的空间位置给出相应原子坐标的势能;
(2)几何优化:
系统的修改原子坐标使原子的三维构象能量最小化;(3)性质计算:
预测某些物理化学性质,如电荷、偶极矩、生成热等;(4)构象搜索:
寻找能量最低的构象;(5)分子动力学模拟:
模拟分子的构象变化。
四、计算化学中的基本概念
1、坐标系统笛卡尔坐标和内坐标
2、原子类型3、势能面4、面积
VanderWaals面积:
空间填充模型或CPK模型。
分子面积:
试探分子(常为半径Å的水分子)在VanderWaals面积上滚动的面积(接触面积+悬空面积)。
可接近面积:
试探分子原点产生的面积。
5、单位:
键长多用Å(埃,angstroms),键能多用kcal/mol表示。
五、计算机辅助药物设计软件及限制
目前CADD存在的问题:
蛋白质结构三维结构的真实性和可用性问题(细胞膜上的受体或跨膜蛋白离开原先环境,空间排列会发生很大变化,难以得到真实的三维空间结构;大量受体结构未知;很多受体只有一级结构,获得的三维结构有限);受体-配体相互作用的方式问题;设计的分子能否进行化学合成;药物体内转运、代谢和体内毒副作用问题等。
第2章分子力学
分子力学(molecularmechanics)是基于原子间存在化学键、非键原子之间的范德华及静电相互作用这一经典理论,通过分子集合、能量、振动光谱及其他物理性质的计算寻求分子的平衡构型及能量,确定有机分子的结构、构象、能量及动力学模型。
其计算忽略了电子的贡献,只考虑核。
一般分子动力学软件提供三种位能面采样算法:
单点(singlepoint):
位能面上某一点计算,给出系统能量和梯度
几何优化(geometryoptimization):
单点采样,梯度为零,局部最小
分子动力学(moleculardynamics):
增加动能,按Newton定律运动
主要用于能量最小化和构象搜索。
受体结构已知,计算R-L结合能;受体未知通过已知配体导出药效团模型(pharmacophoremapping)。
一、理论简介
分子力学基本思想是通过选择一套势函数和从实验中得到的一套力常数,从给定的分子体系原子的空间坐标的初值,用分子力场描述的体系总能量对于原子坐标的梯度,通过多次迭代的数值算法来得到合理的分子体系的结构。
分子力学优化只能是局部优化,所以选择初始构象非常关键。
若为了找到全局能量最低构象,须将所有可能的初始构象分别进行优化,最后进行比较确定分子体系的最优构象。
二、分子力场
如果一个解析表达式能拟合位能面,则此解析表达式就成为分子力场(forcefield),亦即一个力场的确定就是选择解析函数形式和确定参数。
分子力学用几个典型结构参数和作用力来描述结构的变化,由分子内相互作用和分子间相互作用(静电、H-bond、vdW)构成。
常用的分子力场:
其中,适用于大分子的力场:
AMBER,OPLS/AMBER(Schrödinger),CHARMM,CHARMm(Accelrys),GROMOS,ECEPP(freeenergyforcefield)。
适用于小分子的力场:
MM3,CFF95(Biosym/Accelrys),MMFF94,UFF
MM3:
MM力场对静电的相互作用采用键的偶极方法,对于极性或电荷系统不能充分模拟,适用于非极性小分子的结构和热动力学模拟。
MM4:
MM力场的最新发展。
AMBER力场:
合并原子力场(基团合并计算,简化计算,精度下降),对有机小分子或无机分子缺乏有效的参数和原子类型定义,不适合用于小分子。
OPLS力场:
用于蛋白质和核酸大分子。
CHARMM力场:
适用于生物大分子,最新版本采用实验数据和从头算法数据,充分考虑溶剂和溶剂、溶质和溶质、溶剂和溶质之间的相互租用。
BIO+力场:
CHARMM力场的变种。
MMFF94力场:
运用凝聚态过程,适合大分子和小分子,目前最精确的力场之一。
Universal力场(UFF):
针对整个周期表的分子力学和动态模拟力场,静电电荷计算。
COMPASS力场:
可以计算孤立态和凝聚态分子,可用于有机、无机小分子,高分子,金属离子、金属氧化物与金属。
三、能量最小化
量子力学方法可以处理分子的键的断裂与生成并提供系统电子分布的信息,分子力学只能描述核间相互作用的能量变化。
能量最小化的计算按是否采用能量的导数分为两类:
非导数法(单纯型法)和导数法。
一阶能量导数的方向指向能量最小化的点,有最陡下降法(SD)、共轭梯度法(CONJ)、任意步长逼近法;二阶能量导数预测何处能量梯度方向发生变化,有牛顿-拉普森(Newton-Raphson)法。
SD:
梯度是进行搜索的方向,每次搜索之后旧的方向被新点处的梯度取代,对梯度依赖大,收敛慢,适合优化最初段,用于大分子。
CONJ:
不仅运用当前梯度,也采用先前的最小化历史来确定下一步,收敛比SD快,用于大分子。
Newton-Raphson法:
原则上可一步收敛,但存储导数的矩阵太大,不适用于大体系,用于小分子。
四、构象分析
构象搜寻采用适当的方法产生各种不同的构象并对这些狗系那个进行能量最小化,比较好泽泻构象并找出其中能量最低的构象。
根据产生构象的方法不同,可分为系统搜寻和非系统搜寻两类(模拟退火动力学,高温淬火动力学,遗传算法,距离几何法)。
(一)分子动力学
分子动力学是将系统运动方程组数值积分,得到N个经典例子系统的相轨道,进而研究该系统的平衡热力学性质、结构动力学性质、非平衡输运性质等。
分子动力学涉及Newton运动方程的积分,需要选取适当的时间步长,选取的时间步长和运动的频率有关。
分子动力学模拟经过三个阶段:
加热、模拟(平衡期和资料收集期)、冷却。
构建分子的时候是0K,缓慢加热(目的使每个步长里分子都接近平衡,避免高温区发生断裂或扭曲)使系统升温至模拟温度,平衡可以避免加热过程中引入人为因素,平衡所需时间(一般100ps)与体系所模拟的性质有关。
分子体系冷却可降低在较高温度时分子的张力,冷却过程也叫模拟退火。
分子动力学模拟溶剂的作用,可通过选择媒介的介电常数和周围边界条件模拟,水的介电常数。
采用周期边界条件模拟是为了消除刚性壁边界条件和自然边界条件的表面效应。
分子力学模拟一般都应先对分子进行几何优化,在几何优化基础上再对分子进行动力学模拟。
主要用于能量最小化,和分子力学的区别:
分子力学不能越过一定的能垒,只是局部优化;分子动力学模拟则是全局优化低能构象(越过能垒高度与模拟温度有关)。
1、分子动力学两种取样方式
(1)模拟退火动力学
首先体系升温,使分子体系有足够的能量,客服柔性分子中存在的各种旋转能垒和顺反异构能垒,搜寻全部构象空间,在构象空间中选出一些能量相对极小的构象,然后逐步降温,再进行分子动力学模拟。
(2)高温淬火动力学
对于相对比较深的局部能量最小点,需要提高温度600-1200K模拟,往往还需要结合模拟退火。
2、Langevin
Langevin模拟是随机动力学模拟,通过给各原子分配分子在溶剂中与其它分子相互碰撞并随着在溶剂中运动产生摩擦力的值进行模拟(可模拟所有与溶剂及相关的问题),计算速度比动力学快。
(二)随机构象搜索
MonteCarlo方法:
利用随机取样处理问题的方法称为MonteCarlo方法,它是一种通过的采取随机数和概率统计进行猜测来研究问题。
分子动力学不能越过的能垒,MonteCarlo构象搜寻可以是跳跃式的,其优点是取样的构象恰当,对低能构象取样几率大。
Metropolis取样法:
体系首先有一个起始构象,构象发生微小的随机变化产生新的构象,相应的能量发生变化。
如果能量差不大于零则接受构象变化,新构象成为下一步的初始构象。
如果能量大于零,选择一个随机数将其与原相比较,若能量小于则接受,新的构象成为下一次随机变化的起始点;否则拒绝变化,老的构象仍是下一次随机变化的起始点(过程中产生Boltzmann集合)。
(三)系统构象搜索
通过系统地搜寻分子的争购构象空间,寻找势能面上的极小点。
最基本的搜索方法是格点搜索,即在分子构象空间中以柔性键的旋转角度小间隔为变量逐渐搜索。
难以用于大分子体系,也不适用于环状结构。
(四)遗传算法
算法步骤:
随机产生初始群体,个体以二进制序列标记,计算适应值。
通过三个算子产生新一代,复制算子、杂交算子、变异算子。
对于构象搜索,适应值为能量的函数。
主要用于全局优化。
(五)距离几何法
第3章量子力学
量子化学可以用于:
(1)用于解释观察的实验结果,补充实验条件达不到因而不能观察的一些实验结果;
(2)预见分子的结构、过渡态和反应途径;(3)预见分子的电学、磁学和光学性质;(4)应用于从X射线到NMR等光谱中;(5)应用于大分子、晶体和溶液中分子间的相互作用。
一、量子化学理论简介
应用量子力学原理处理化学问题,形成分子轨道理论价键理论与配位场理论,自六十年代以来几乎所有的量子化学计算都是以分子轨道理论为基础的。
分子轨道理论在物理模型上有三个基本近似:
非相对论近似、Born-Oppenheimer近似(绝热近似)及单电子近似(轨道近似)。
非相对近似是电子在原子核附近运动而不被的原子核俘获,必须保持很高运动度。
近似认为电子质量等于电子静止的质量,即电子质量恒为1个单位。
Born-Oppenheimer近似是在计算分子总能量(电子总能量与核排斥能值和)时,把电子的运动和核的运动分开处理,忽略贡献很小的电子态之间偶合项,即非绝热项,故又称绝热近似。
上面两种近似主要是对Hamilton算符的处理。
单电子近似则认为电子被函数为n个电子所占据的轨道(单电子函数)的乘积。
单电子近似导致分子轨道理论的一个核心概念——自洽场的出现。
每一个电子的运动状态不仅受核力的作用,而且还要受其他n-1个电子所产生的势场的作用。
这样在描述与电子势场时,必须考虑受作用的那一个电子状态和它本身作为其他电子运动时对势场的贡献要一致起来,即自洽起来。
这样的势场,成为自洽场。
1、Hartree-Fock方法固有的