分子动力学学习笔记3

4--肾上腺素受体的同源模建及与的对接研究
在pdb库中搜索结构，在insightii中模建，affinity 对接。

monte carlo用方法预优化搜索10 个能量最小化的合理的起始构象；以得到的10个最低能量构象为起始点， 用多极单元非键力方法cell muitipole 进行模拟退火动力学优化，起始温度500k 经50步控温冷却到300k。
利用discoveer软件cvff力场 先分步对氢原子进行最陡下降法优化及共轭梯度法优化。再对重原子进行束缚限制优化，经优化后的主链二面角值处于正常区域范围内的高于91%。考虑到对接结果与同源模建条件的一致性，计
算中仍然采用cvff力场， 每个计算过程得到十个最低能量构象，代表十种可能的受体配体结合方式。我们进行了十次对接研究。

5--


6--蛋白质结构和功能的紫外共振喇曼光谱研究
研究蛋白质分子结构的方法分为两类:一类是测定蛋白质晶体的分子结构,如X 射线衍射分析法和小角中子衍射法;另一类是研究溶液状态下的蛋白质分子,如核磁共振法、圆二色谱法、荧光光谱法和激光喇曼光谱法等。因为生命体中的蛋白质是以溶液状态存在,所以第二类方法具有更重要的意义。
理论研究预测,α2螺旋结构的变性发生在亚纳秒量级,而二级结构模体发生在纳秒量级,近来荧光光谱和红外光谱常用来在纳秒量级范围内研究蛋白质的折叠和去折叠过程。用泵浦2探测技术可以研究蛋白质与肽的折叠和去折叠现象

7--中国对虾热休克蛋白70 (HSP70) 的三维结构优化
运用GROMOS96 的43A1 力场〔6〕,采用充分含水模型. 在线预测得到两个结构域,即ATP 酶结构域(N 端) 和肽结合结构域(C 端) ,将其分别放入立方体周期性盒子里,确定溶质到盒边的最小间距( 见表1) , 加入水分子( SPC 模型) , 利用SHAKE算法,进行400 步(每步0. 002ps ,下同) 的能量优化. 为模拟溶液中性环境,加入适当的离子(见表1) ,再进行同样的能量优化. 取随机种子数为210185 , 玻耳兹曼常数为8. 31441 ×10 - 3KJMol - 1K- 1 ,压强为1 个大气压,采用逐步升温技术,从50K升至300K(见表2) ,进行总计150ps 的约束动力学模拟,在室温(300K) 和常压(1bar) ,进行4ns 的自由动力学模拟,系统达到平衡.
自由动力学模拟残基Cα原子的均方根偏差
(Root - Mean - Square Deviation ,RMSD) 如图1 所示. 图1 (a) 为N 末端结构域,3～5ns 的RMSD 的平均值为0. 537nm , 最大相对偏差为15. 2 % ,RMSD 的均方根涨落Rdc 为5. 7 %;图1 (b) 为C 末端结构域, 最后500ps 的RMSD 的平均值为0. 354nm ,最大相对偏差为8 % ,RMSD 的均方根涨
落Rdc 为2. 5 %. 由于受计算能力的限制及热运动造成原子位置的涨落等因素的影响,要使系统达到完全平衡需花费较长的时间,根据准平衡态概念〔7〕,可以认为已达到准平衡态. 由最后500ps的轨迹,求得其平均坐标,得到如图2 (a) 和(c) 所示中国对虾热休克蛋白70 的优化结构.
8--鱼类Ⅱ型抗冻蛋白冰水复合体系的理论模拟研究
　鱼类Ⅱ型抗冻蛋白冰水复合体系的构建　以鱼类Ⅱ型抗冻蛋白实验测定的结构为模型,在已有的理论模拟结果基础上,选定将鱼类Ⅱ型抗冻蛋白与冰晶的{1121}晶面作为最佳束缚冰面[6 ] ,如图1 所示,将在真空条件下所获得的这种最佳作用模型,放入体积为5. 5 nm ×5. 0 nm ×5. 6 nm 的周期性水溶剂盒子中,重新构建出抗冻蛋白冰水复合体系,如图2 所示. 溶剂水分子选用HyperChem 数据库提供的TIP3P 模型[7 ] . 复合体系中共包含1 874 个蛋白原子,1 380 个冰晶原子以及作为介质的1 695 个水分子
9--的中间态构象是α核被破坏而β核仍保持天然构象
离散分子动力学(DMD)方法近来得到了广泛的应用. DMD 方法采用类似于MD 方法的非格点模型. 但不同于MD 方法中使用的连续势(例如Lennard2Jones 势) ,DMD 方法使用方势阱来表示基元间的相互作用. 通过变MD 方法中的求解经典运动方程为DMD 方法中的求解瞬间碰撞方程,从而提高运算效率

标签：

版权申明：本站文章部分自网络，如有侵权，请联系：west999com@outlook.com
特别注意：本站所有转载文章言论不代表本站观点，本站所提供的摄影照片，插画，设计作品，如需使用，请与原作者联系，版权归原作者所有