极化电荷的发展、应用以及生物大分子体系分子动力学模拟的研究

【摘要】分子动力学模拟可以用来研究生物大分子体系的各种性质,如：预测受体和配体之间的亲和力,模拟蛋白质的折叠过程,研究蛋白质与细胞膜之间的相互作用、离子通道以及酶的催化反应等等,因此逐渐受到科学家的重视。分子动力学模拟理论将每个原子看成一个球,原子与原子之间通过弹簧连接,模拟体系在高维的势能函数面上运动。目前,存在多种不同的势能函数,如：AMBER,CHARMM,OPLS等,所有的势能函数都可划分成两大类：通过共价键相互作用的势能和非共价键作用的势能,共价键作用能包含键的伸缩能,角的弯曲能,二面角的扭转势能,非共价键作用能包含静电作用能和范德华作用能,势能函数中不同的参数形成不同的分子力场。静电作用是一种长程作用力,对动力学模拟的结果具有重大的影响,静电作用力的正确与否直接关系到氢键、盐桥等作用力的合理性,而大量研究证明氢键和盐桥对蛋白质折叠、配体和受体相互作用等都具有非常重要的意义。目前所有的力场都是采用库伦定律来计算静电作用能,分子理论假设每个原子上的电荷都集中在原子的中心,通过拟合分子周围的电势得到每个原子上的电荷。传统分子力场中,每个氨基酸上的原子电荷为一个定值,不受周围环境的影响。而实际上,生物大分子的一切运动都离不开水溶液,水分子为极性分子与生物大分子之间存在明显的极化效应,因此我们在分子动力学模拟的时候必须考虑生物大分子与周围环境之间的极化作用。经过长期的努力,目前已发展了多种极化电荷模型,主要有：Drude振子模型,涨落电荷模型和诱导多极模型。但是这些方法仍然存在一些无法避免的缺陷,如：Drude振子模型引入了人为的参数,诱导多极模型的计算量非常大,目前根本无法与分子动力学模拟结合,而且力场的可迁移性也很难实现。我们课题组完全放弃了力场的可迁移性,运用分子碎片共轭帽法,将量化计算的方法引入到生物大分子的电势计算过程中,在计算电势的同时考虑了周围环境与生物大分子之间的极化效应,最后对电势进行拟合得到生物大分子的极化电荷。借助RESP方法拟合的蛋白质内部原子的极化电荷随构象的变化波动幅度较大,非常不利于分子动力学模拟,特别是在QM/MM模拟等过程中,可能存在大构象的变化,因此大幅度的电荷涨落,对模拟非常的不利。基于2011年张等人选择参考电荷将偶极矩进行分割的方法,本文在RESP方法的基础上发展了一套在数值上更稳定的拟合极化电荷的方法：dRESP方法。dRESP方法中,每个原子上的电荷被划分成两部分：参考电荷和需要拟合的电荷。此方法的前提是参考电荷能够近似描述原子的电荷分布,所以参考电荷一般为MBER、CHARM、OPLS等力场中常用的电荷。量化计算的电势与参考电荷在相同格点上产生的电势之间的差量就是我们需要拟合的电荷在这个格点上产生的电势,拟合这个电势就能得到一个差值电荷。在拟合的过程中,我们在每个原子上所加的限制都是不同,与其参考电荷平方的倒数成正比,最后参考电荷与拟合得到的差值电荷之和就是整个原子的电荷。通过这种方式拟合电荷的优势是：拟合的电荷比原子电荷小了几个数量级,因此造成的数值误差比直接拟合原子的电荷要小很多。计算结果证明参考电荷的选取不影响拟合电荷的波动性,这种分步拟合电荷的方法的确很大程度上避免了拟合过程中出现的数值问题。Streptavidin与biotin之间的结合自由能高达-18.1kcal/mol,是目前发现的自然界中亲和力最强的化合物,对合理药性物设计具有很高的研究价值。大量研究证明氢键网络对biotin与streptavidin之间的结合非常重要,然后Baugh等人发现streptavidin末端的130号氨基酸由PHE突变成LEU之后,结合自由能减少了4.2kcal/mol,但是结构却没有明显的变化,突变前后A链上Ca原子的RMSD值仅为0.377A。我们同时运用AMBER电荷和dRESP方法拟合的极化电荷对streptavidin突变前后与biotin的复合物进行分子动力学模拟,并借助MM/PBSA方法预测原生态和130氨基酸突变体系的结合自由能。AMBER电荷的计算结果发现biotin与streptavidin之间的亲和力为正值,两者无法自由结合,同时计算的突变引起的能量差与实验值相比也过低。dRESP极化电荷加强了biotin与streptavidin之间的结合能,计算的绝对结合能和原生态与突变态之间的相对能量都与实验值一致。我们的计算证明了Baugh等人的猜测：静电极化效应是引起130号氨基酸突变能量变化的主要原因。CYP是一大类单氧酶,涉及到70-80%药物的前期代谢反应,因此对药物-药物相互作用具有重大的影响。CYP2J2是CYP酶家族的一类子族,主要表达在小肠,心血管等肝外器官中,大量研究证明CYP2J2能够催化大部分的药物,是药物前期代谢必不可少的一种酶,也可能会涉及到药物-药物相互作用,但是CYP2J2药物代谢的机理还不明了。分子对接是计算机辅助药物设计的主要手段,结合实验、分子对接、分子动力学模拟等技术,我们从69种已知的药物中筛选到两种CYP2J2有效的抑制剂：telmisartan、flunarizine、telmisartan和flunarizine对CYP2J2的半抑制剂浓度分别为0.42声μM和0.94μM,而对CYP家族其它成员的活性几乎没有影响,同时我们还证明这两种抑制剂不是CYP2J2的底物,是研究CYP2J2药物代谢机理最理想的工具。分子对接和动力学模拟进一步证实了实验的结论,并在原子结构的级别上展示telmisartan是非竞争性抑制剂,flunarizine是竞争性的抑制剂。
【作者】曾娟；
【导师】张增辉；梅晔；
【作者基本信息】华东师范大学，原子与分子物理，2014，博士
【关键词】MD模拟；分子力场；静电作用；极化效应；分子碎片共轭帽法；PPC；RESP；dRESP；AMBER；Streptavidin；Biotin；溶剂模型；MM/PBSA；结合自由能；CYP2J2；分子对接；药物-药物作用；

【参考文献】
[1]孟琳.北京城市轨道交通商业物业价值评估[D].北京交通大学,交通运输规划与管理,2013,硕士.
[2]薛美东,赵波,张雪松,高小淇,江全元.基于分布式控制的独立型光储水柴微网调度策略[J].电力系统自动化,2014,04:1-7.
[3]宋亚男.基层民警职业道德建设的思考[D].吉林大学,公共管理,2014,硕士.
[4]张胡佳.“枢纽型”社会组织运作机制研究[D].华东理工大学,公共管理（专业学位）,2013,硕士.
[5]王磊.天麻流浸膏提取工艺的优化及质量控制的研究[D].黑龙江中医药大学,药物分析学,2013,硕士.
[6]刘静.卡通元素融入现代中国画创作的探索与研究[D].首都师范大学,美术,2014,硕士.
[7]鲁俊毅.企业核心竞争力及核心技术创新研究——以爱生技术集团为例[D].西北工业大学,企业管理,2004,硕士.
[8]南楠.中国传统造型艺术中鹿形象研究[D].西安美术学院,美术学,2014,硕士.
[9]李琨.蓝宝石回音壁模温度计的探索研究[D].河北大学,仪器仪表工程,2014,硕士.
[10]黄初生.白细胞滤过在体外循环中应用的Meta分析[D].广西医科大学,胸心外科（专业学位）,2013,硕士.
[11]李莎莎.网络消费对职业女性生活方式的影响研究[D].郑州大学,社会学,2013,硕士.
[12]马长信.赣北中志留统风暴岩的初步研究[J].江西地质科技.1995(02)
[13]苟文蓉.操控论观照下葛浩文英译《天堂蒜薹之歌》研究[D].西北师范大学,英语语言文学,2014,硕士.
[14]田静华.MPEG-4标准的研究及其运动估计部分的实现[D].浙江工业大学,通信与信息系统,2004,硕士.
[15]程胜,段绪果,吴敬.重组蔗糖异构酶的制备及应用条件优化[J].食品与发酵工业.
[16]杨楠.薪酬满意度对组织承诺影响的实证研究[D].辽宁大学,企业管理,2012,硕士.
[17]侯永慧.唐代诗刻文献研究[D].广西师范大学,中国古典文献学,2014,硕士.
[18]张建伟.高超声速再入体气动热数值模拟研究[D].山东大学,2012.
[19]许云光.水力发电机组振动监测与信号分析系统的研究[D].哈尔滨工业大学,控制科学与工程,2013,硕士.
[20]何家洪,张胜涛,邢少华,宋仲容.L-半胱氨酸/普鲁士蓝复合修饰玻碳电极的制备及电化学性能[J].分析化学,2013,10:1561-1565.
[21]张云程.带式输送机中间带条驱动用于下运的探讨[J].煤炭科学技术.1989(03)
[22]贺进,王东瑞,党智敏.乳液混合法制备低渗流阈值、高导电聚苯乙烯/石墨烯复合颗粒[A].中国化学会高分子学科委员会.2013年全国高分子学术论文报告会论文摘要集——主题J：高分子复合体系[C].中国化学会高分子学科委员会:,2013:1.
[23]刘改.思想政治教育的社会管理功能研究[D].中国矿业大学,思想政治教育,2014,硕士.
[24]徐焕祥,李道飞,王雷,叶锦,樊之鹏,俞小莉.回收余热能的气动/内燃复合循环效率[J].浙江大学学报(工学版),2014,04:649-652+703.
[25]路甬祥.关于材料科技与产业的回眸与思考[J].科技导报,2012,30:15-19.
[26]邢志广.基于新国标的汽车用LED前照灯设计[D].天津理工大学,光学工程,2013,硕士.
[27]谢佳.基于WebService的B/S架构的在线考试系统的设计与实现[D].西安科技大学,计算机技术,2013,硕士.
[28]栾慧鑫.从跨文化视角看《孙子兵法》在德国的传播[D].黑龙江大学,德语语言文学,2013,硕士.
[29]胡粤.区域旅游经济增长的差距与收敛性研究[D].湖南师范大学,旅游管理,2014,硕士.
[30]刘春梅.基于预制语块的中职英语博客写作教学研究[D].河北大学,现代教育技术,2013,硕士.
[31]肖兆权.腐败成因的经济学分析和国家廉政体系构建的研究[D].河海大学,产业经济学,2004,硕士.
[32]卢利平.美科学家首次用碳纳米管制造出石墨烯带[J].功能材料信息,2009,03:51.
[33]张立志.农村高中地理与经济社会协调发展关系研究[D].河北师范大学,学科教学,2014,硕士.
[34]王巍.大学生创业模式研究[D].吉林大学,工商管理,2004,硕士.
[35]张淑华.新闻传播专业学位研究生教育的困境与对策研究[J].郑州大学学报(哲学社会科学版),2014,03:187-192.
[36]张新.基于CONVEF的扩散模型和图像去噪方法的研究[D].天津理工大学,计算机技术,2013,硕士.
[37]申权.异形坯连铸结晶器流场与锥度的研究[D].燕山大学,重型装备设计理论及其数字化技术,2013,硕士.
[38]白夏夏.延安市产业结构现状及优化模型[D].西安科技大学,应用数学,2013,硕士.
[39]韩周琨.半个世纪以来李清照形象中西对比研究[D].西南民族大学,外国语言学及应用语言学,2014,硕士.
[40]王珊珊.单通道角跟踪接收机基带信号处理技术[D].西安电子科技大学,电子与通信工程,2012,硕士.
[41]岳琦.董事会特征与财务舞弊相关性研究[D].内蒙古财经大学,会计学,2012,硕士.
[42]邓帅奇.全空间弹性波场数值模拟与逆时偏移成像方法研究[D].中国矿业大学,2012.
[43]杜波.EDXRF分析中吸收—增强效应的蒙特卡罗模拟[D].成都理工大学,辐射防护及环境保护,2013,硕士.
[44]丁琳,吴长永.数据挖掘在远程教育个性化服务中的应用[J].电化教育研究,2002,09:43-46.
[45]孙琦.高校马克思主义大众化教育主体合力机制研究[D].哈尔滨理工大学,思想政治教育,2012,硕士.
[46]郭彦杰.我国优秀运动员退役安置政策变迁的研究[D].北京体育大学,体育人文社会学,2013,硕士.
[47]韩亚斌.基于积件思想的初中物理声、光现象教学研究[D].内蒙古师范大学,学科教学（专业学位）,2012,硕士.
[48]李明.基于决策树的数据检测模式获取[J].计算机时代,2005,07:37-38.
[49]杨柳.镁合金方盒形件温热拉深成形试验研究[D].燕山大学,材料加工工程,2014,硕士.
[50]黄平,路中侃.川东石炭系气水分布及预测[J].石油地球物理勘探,1996,S2:95-100+151.