分子模拟与药物设计

1、分子模拟与药物设计分子模拟与药物设计 研究蛋白质与小分子药物之间的相互作用研究蛋白质与小分子药物之间的相互作用 基于配体的药物设计基于配体的药物设计QSARwithCoMFA可以建立分子性质(包括生物活性)与它们结构之间相关的统计学与图形化模型。这些模型可用于预测新化合物的性质或活性。定量构效关系（QSARs）将分子的化学性质或生物活性与它的结构相关以便设计出活性更好的新化合物凝血酶抑制剂的CoMSIA分子场等势图（左：立体场，绿色有利黄色不利；右：疏水场，蓝色有利红色不利） 基于配体的药物设计基于配体的药物设计叠合一系列分子从而共享出体现生物活性的公共模式，生长出相应的药效团假设采用先进的遗

2、传算法和多元化的打分函数，考虑分子能量、立体相似性和药效位点重合，同时还考虑构象的柔性，不确定的立体化学性质，可变的环结构 基于受体的药物设计基于受体的药物设计高通量虚拟筛选高通量虚拟筛选大量化合物数据库的筛选是一项昂贵费时的任务。计算高通量筛选可以增加筛选数据集中合适化合物的占有率并且能降低先导物发现的成本打分函数是一个蛋白-配体原子表面距离的非线性函数的线性组合。蛋白-配体相互作用包括立体作用、极性作用、熵和溶剂化作用。 基于受体的药物设计基于受体的药物设计同源模建是一个根据模板蛋白将一级序列转成3D结构的技术总称，包括了threading和homology两种人类Xa因子的晶体结构（PD

3、BID：2BOK，橙色）与同源模建结构（白色）的比较 分子动力学分子动力学主要内容：主要内容： Gromacs简介简介 Gromacs原理原理 Gromacs运算步骤运算步骤 Gromacs结果分析结果分析Lysozyme in WaterRMSD 分子模拟是在分子模型的基础上用计算机做实验，“计算机实验” 通过模拟微观粒子的运动来计算宏观性质温度压力黏度传递性质表面张力分子间的作用模型牛顿力学量子力学统计力学等 分子模拟具有理论和实验的双重性质 分子模拟不能完全取代实验理论实验模拟理论的正确性模拟参数的正确性模拟方法的选择理论的更新模拟分子在特定环境下一定时间内的构象能量变化趋势；分子动力学

4、模拟是根据分子力学性质建立的适用于生化体系，聚合物，金属或非金属材料的力场和牛顿的运动力学原理发展出来的计算方法；分子力场(forcefield)中的各项参数，包括键长，键角，电荷分布等都可以通过量子化学计算得到；与量化计算相比，分子力学在速度方面具有明显的优势；最早出现在上世纪70年代；分子动力学对势能函数的依赖性：所有从分子动力学计算出来得到的宏观性质最终都取决于势能模型分子动力学的核心：牛顿运动方程势能（位能）模型：( )( )U rf r=iiFtrm22i=1,2,3,NrUr例：甲烷，某些惰性气体质点处理Ur方阱模型Ur阶梯模型键的振动键角扭矩分子内部各原子（基团）之间的范德华力、

5、静电力一般要计算1-4（相隔超过两个键的原子或基团对）15432qqq分子之间的范德华力分子之间的静电力CCCHHHHHHHH 10根键长作用 18个键角作用 8个扭矩作用 27个范德华力作用 27个静电作用初始化能量优化平衡数据产出避免局部分子重叠，并不是动力学模拟 根据所有分子的当前坐标计算个分子的受力（位能函数） 根据受力更新分子的坐标 在此过程中收集用来计算宏观性质的有关信息读入模型参数，模拟控制参数 去除某些可能存在的原子重叠 去除某些严重扭曲的键长、键角、扭矩等 方法最速下降法牛顿拉夫森方法其他一般优化几千到几万步简单小型体系大型（复杂）体系和并行计算 气体的模拟 小分子体系，不需

6、要复杂的势能模型 几百到几千个分子，分子分布稀疏，大部分是短程作用 一般用一台微机就可以处理，计算时间几分钟几小时 简单的液体，不涉及太多的界面性质 小分子体系，势能模型不是很复杂 几百个分子，可能涉及到静电作用，可能需要长程校正 用微机也可以处理，计算时间一般几小时几天 必要性 体系越来越大 模拟时间越来越长 解决办法 制造更快的处理器 并行计算机 例子：50000原子的生物体系，1ns模拟 单个处理器：12天 16个并行处理器：1天或者 MessagePassingInterface 90年代初制定和完善的一套并行语法 支持Fortran,C,C+ 简单易学NAMDAMBERCHARMMT

7、INKERLAMMPSDL-POLYGROMACS 主要针对生物体系，也适当照顾一般化学体系 优点 算法好，计算效率高 界面友好 维护服务好 免费软件 缺点 兼容性不好/nSelect the Force Field: From /usr/local/gromacs/share/gromacs/top:n 1: AMBER03 force field (Duan et al., J. Comp. Chem. 24, 1999-2012, 2003)n 2: AMBER94 force field (Cornell et al., JACS 117, 5

8、179-5197, 1995)n 3: AMBER96 force field (Kollman et al., Acc. Chem. Res. 29, 461-469, 1996) n4: AMBER99 force field (Wang et al., J. Comp. Chem. 21, 1049-1074, 2000)n 5: AMBER99SB force field (Hornak et al., Proteins 65, 712-725, 2006)n 6: AMBER99SB-ILDN force field (Lindorff-Larsen et al., Proteins

9、 78, 1950-58, 2010) n 7: AMBERGS force field (Garcia & Sanbonmatsu, PNAS 99, 2782-2787, 2002)n 8: CHARMM27 all-atom force field (with CMAP) - version 2.0 n 9: GROMOS96 43a1 force fieldn10: GROMOS96 43a2 force field (improved alkane dihedrals) n11: GROMOS96 45a3 force field (Schuler JCC 2001 22 1

10、205) n12: GROMOS96 53a5 force field (JCC 2004 vol 25 pag 1656)n13: GROMOS96 53a6 force field (JCC 2004 vol 25 pag 1656) n14: OPLS-AA/L all-atom force field (2001 aminoacid dihedrals)n15: DEPRECATED Encad all-atom force field, using full solvent chargesn16: DEPRECATED Encad all-atom force field, usin

11、g scaled-down vacuum chargesn17: DEPRECATED Gromacs force field (see manual) n18: DEPRECATED Gromacs force field with hydrogens for NMR pdb2gmx-f1AKI.pdb-o1AKI_processed.gro-waterspce-f蛋白质名称-o坐标water水分子类型Step One: Prepare the Topology#include“oplsaa.ff/forcefield.itp“立场文件;Namenrexcl包含的分子名称Protein_A3

12、atoms每个原子的坐标、电荷、质量、类型信息;nrtyperesnrresidueatomcgnrchargemasstypeBchargeBmassB;residue1LYSrtpLYSHq+2.01opls_2871LYSN1-0.314.0067;qtot-0.32opls_2901LYSH110.331.008;qtot0.033opls_2901LYSH210.331.008;qtot0.364opls_2901LYSH310.331.008;qtot0.695opls_293B1LYSCA10.2512.011;qtot0.94Step Two: Examine the Topo

13、logyStep Three: Defining the Unit Cell & Adding Solvent1. Define the box dimensions using editconf.2. Fill the box with water using genbox.editconf -f 1AKI_processed.gro -o 1AKI_newbox.gro -c -d 1.0 -bt cubic-c 居中 d 盒子边缘离质心的距离 bt 盒子类型genbox -cp 1AKI_newbox.gro -cs spc216.gro -o 1AKI_solv.gro -p

14、topol.top-cs 水分子类型Step Four: Adding Ionsgrompp -f ions.mdp -c 1AKI_solv.gro -p topol.top -o ions.tpr-f 输入mdp文件 o生成tpr文件mdp 即molecule dynamics parameter文件，起控制动力学模拟过程的作用genion -s ions.tpr -o 1AKI_solv_ions.gro -p topol.top -pname NA -nname CL -nn 8-pname 阳离子类型 nname 阴离子类型 nn 离子数目Step Five: Energy Mini

15、mizationgrompp -f minim.mdp -c 1AKI_solv_ions.gro -p topol.top -o em.tprmdrun -v -deffnm em-v 输出每步的运算信息 deffnm 定义输出文件的名称em.log: EM过程的记录文件em.edr: 二进制能量文件em.trr: 二进制轨迹文件em.gro: 能量最小化后的结构Step Five: Energy Minimization能量分析：g_energy -f em.edr -o potential.xvgStep Six: Equilibration目的：1.使蛋白质分子周围的溶剂和离子达到平衡

16、状态,防止整个体系坍塌；2.先等温使得溶剂分子达到既定的模拟温度；3.再等压使得整个体系的密度达到均衡一致；grompp -f nvt.mdp -c em.gro -p topol.top -o nvt.tpr（等温恒容过程）mdrun -deffnm nvtg_energy -f nvt.edrStep Six: Equilibrationgrompp -f npt.mdp -c nvt.gro -t nvt.cpt -p topol.top -o npt.tpr（等温恒压过程）mdrun -deffnm nptg_energy -f npt.edr -o pressure.xvg（分析压力

17、）g_energy -f npt.edr -o density.xvg（分析密度）Step Eight: Production MDgrompp -f md.mdp -c npt.gro -t npt.cpt -p topol.top -o md_0_1.tpr mdrun -deffnm md_0_1mdrun -deffnm md_0_1mpirun -np X mdrun_mpi -deffnm md_0_1（并行）Step Nine: Analysistrjconv -s md_0_1.tpr -f md_0_1.xtc -o md_0_1_noPBC.xtc -pbc mol -ur

18、 compact-pbc周期性边界条件 -xtc文件压缩后的坐标文件g_rms -s md_0_1.tpr -f md_0_1_noPBC.xtc -o rmsd.xvg -tu ns-tu 时间单位nsg_rms -s em.tpr -f md_0_1_noPBC.xtc -o rmsd_xtal.xvg -tu ns与结晶结构相比整个体系的RMSD值g_gyrate -s md_0_1.tpr -f md_0_1_noPBC.xtc -o gyrate.xvg回转半径计算，代表整个蛋白质体系的紧致程度，评价体系是否稳定的指标。Step Nine: AnalysisStep Nine: An

19、alysisngmx f md.trr (or md.xtc) s md.tpr（GROMACS自带图形分析查看模块）Step Nine: Analysismake_ndx f xxx.pdb o xxx.ndx（体系分组）g_confrms f1 1OMB.pdb f2 md.gro o fit.pdb（叠合结构并输出RMSD值）g_rmsf -f xxx.xtc s xxx.tpr b e o xxx.xvg ox xxx.pdb（均方根波动）-b开始时间 e结束时间 ox 输出开始时间到结束时间的平均PDB结构g_hbond f xxx.xtc s xxx.tpr num xxx.xvg输出氢键的数目 r0.35nm 角度小于30trjconv f xxx.xtc s xxx.tpr o xxx.pdb dump xxx在特定的时间点输出该时间点的pdb结构文件 Ste