甘油水溶液氢键特性的分子动力学模拟3
目 录
摘要................................................................................................................................ 1
关键词............................................................................................................................ 1 Abstract ...................................................................................................................... 1
引言................................................................................................................................ 1
1.一些研究者通过经典的结晶理论建立了胞内冰.................................................. 2
2.模拟结果.................................................................................................................. 2
2.1径向分布函数.......................................................................................................... 2
2.2甘油分子的构型分布.............................................................................................. 3
2.3氢键分析.................................................................................................................. 3
2.3.1氢键的选择........................................................................................................... 3
2.3.2氢键结构分析....................................................................................................... 4
2.3.3氢键动力分析....................................................................................................... 4
3.结论.......................................................................................................................... 5
参考文献........................................................................................................................ 5
甘油水溶液氢键特性的分子动力学模拟
摘要:为了研究低温保护剂溶液的结构和物理化学特性,以甘油为保护剂,采用分子动力学方法,对不同浓度的甘油和水的二元体系进行了模拟。得到了不同浓度的甘油水溶液在2 ns 内的分子动力学运动轨迹,通过对后1 ns 内运动轨迹的分析,得到了各个原子对的径向分布函数和甘油分子的构型分布。根据氢键的图形定义,分析了氢键的结构和动力学特性。计算了不同浓度下体系中平均每个原子(O和H) 和分子(甘油和水) 参与氢键个数的百分比分布及其平均值。同时还计算了所有氢键、水分子之间的氢键以及甘油与水分子之间的氢键的生存周期。 关键词:分子动力学模拟; 径向分布函数; 低温保护剂; 水溶液
Abstract :To study the structure and physicochemical characteristics of cry protective agent (CPA) solutions, Glycerol has been chosen as a CPA and the molecular dynamics method was used to simulate glycerol and water binary Systems with different concentrations. Molecular dynamics trajectories of aqueous glycerol solutions within 2 ns were obtained. After a detailed analysis of trajectories within the last 1 ns, the intermolecular radial distribution functions for C-C,C-O,C-H,O-H,O-O and H-H pairs and the backbone conformation distributions of glycerol molecules were calculated. Based on geometrical criteria, structural and dynamics characteristics of the hydrogen bonding network were analyzed. Distribution percentages and average values of the number of hydrogen bonds per atom (O and H atoms) and per molecule(glycerol and water molecules)were calculated. The lifetimes of total hydrogen bonds, hydrogen bonds between water molecules and hydrogen bonds between glycerol and water molecules were also studied.
Key Words :Molecular dynamics simulation; Radial distribution function; Cry protective agent; Aqueous Solution
引言
迄今为止,人类虽然已经成功保存多种细胞和组织,但是对细胞在冷冻和复温过程中的损伤机理还不是十分清楚[1]。为了保存更大体积的生物材料(如人的肾脏、肝脏和心脏) ,有不少研究者仍在从事相关的基础研究工作。胞内冰晶的
形成将导致微细结构的破坏,使细胞死亡。因此最理想的保存方案是玻璃化降温保存,但此过程要求很高的降温速率。目前的降温手段以及样品内温度均匀分布要求还不允许有这么高的降温速率。若能从理论上预测胞内冰晶的形成和生长规律,在此基础上寻找抑制冰晶生成和生长的方法,则可以实现正常降温速率下的玻璃化保存,最大限度地提高细胞冷藏后的存活率。
1、一些研究者通过经典的结晶理论建立了胞内冰
晶生成和生长模型,如Karlsson 的扩散控制模型[2]以及Zhao 的统一模型[3]。我们曾经确定了冷冻过程中胞内溶液的均相成核温度下降值与平衡凝固点下降值的关系[4],并根据此关系,建立了基于软冲突的扩散控制生长模型[5]。虽然我们的模型预测结果与实验结果基本吻合,但是在某些情况下,尤其是在较低的冷冻保护剂(cry protective agent, CPA) 浓度下,模型预测结果与实验结果仍有不小的偏差。同时低温环境下CPA 溶液的物理化学参数,如扩散系数、黏度、成核温度、成核速率等都缺乏精确的数据[5]。氢键的存在使得CPA 溶液的结构和物理化学行为更为复杂,因此,对CPA 溶液内部氢键特性的研究就显得更为迫切和必要。微尺度数值模拟方法的不断成熟和完善,为展开此项研究提供了一种强有力的手段。分子动力学模拟(molecular dynamics simulation, MDS)方法已经被广泛用来研究物质和体系的微观结构及其氢键特性。甘油分子小且柔性强,易于被玻璃化,被认为是最有前途的低温保护溶剂,因此引起了广泛的关注[6-11]。本文以甘油为CPA ,利用分子动力学方法研究不同浓度下甘油水溶液的氢键特性。系中加入不同数量的甘油分子,得到不同浓度的甘油水溶液。
2、模拟结果
2.1径向分布函数
对l ns的分子动力学轨迹进行了分析,得到了各个原子对(包括甘油分子间,水分子间以及甘油与水分子间) 的径向分布函数,如图1所示。甘油浓度对径向分布函数存在一定的影响,但对波峰波谷处的距离基本没有影响。原子对的径向分布函数对了解溶液内部的微观结构有重要的意义。O-H 和O-O 原子对的径向分布函数分别在0.19和0.29 nm 处存在波峰表明着氢键的存在。同时,O-H 和O-O 原子对的径向分布函数第一个波峰的距离恰好相差0.1 nm,表明氢键基本上是线性的。O-O 和O-H 原子对的径向分布函数的第一个波谷的距离dO-O 和dO-H
通常被用来定义氢键[6,21-23]。随甘油浓度不同,dO-H 在0.224-0.243 nm之间变化,dO-O 在0.324-0.343 nm之间变化。
2.2甘油分子的构型分布
甘油分子含有两个CCCO 二面角,根据其大小二面角可以分为α(-60°,60°) ,β(60°,180°)和γ(180°,300°)三种,相应地甘油分子的构型分为6种[6,24]:αα,αβ,αγ,ββ,βγ,γγ。相对于氢键,对甘油分子构型的研究比较少,且主要是采用数值的方法,各种分析结果的偏差比较大[6,24,25]。根据本文的模拟结果,计算了各个浓度下各种构型的分布,结果汇于表2。同时, 表中还列出了各个浓度下统计结果的算术平均值及其它研究者的结果。根据我们的研究结果,αα构型最稳定, 甘油分子构型分布随浓度有微小的变化,但考虑到甘油分子个数等误差因素,可以视为基本没有变化,这也符合其它的研究结果[6]。αα构型广泛存在于甘油晶体态,这也表明了其稳定性[24]。NMR 研究[25]表明水溶液中αα,αβ,αγ,βγ的分布比率要大于ββ,γγ的分布比率, 这与我们的研究结果相符。
2.3氢键分析
2.3.1氢键的选择
在经典分子动力学模拟中,氢键的定义比较模糊[24],通常采用两种方式:能量准则和图形准则。根据能量准则,分子对间的相互作用能强于一定值,便视为该分子对之间存在氢键作用[26];分子间是否存在氢键取决于两个分子间的相对位置[9,24,25,27,28]。氢键X —H…A对应的图形约束包含三个距离和三个角。三个距离分别为X —H 、X —A 和H —A ,其中X —H 为共价键长,取决于采用的势能模型, 因此在氢键的图形准则中,一般采用后两个。三个角分别为X —H…A、H —X…A和X…A…H ,其中前两个角常用来定义氢键。两个分子间的相对位置满足一个(或多个) 距离约束, 该两个分子间就存在氢键作用。角约束通常作为附属准则来限制氢键,因为氢键大体呈线性[9,21,23,24,27,29,30]。本文采用距离H —A 和角H —X…A双重约束来定义氢键。本文研究的甘油水溶液系统可能存在两种形式氢键:O—H…O和C —H…O.C—H…O是相对于O —H…O更弱的一种相互作用
[6],主要存在于蛋白溶液中[31]。
实验发现甘油分子CH 组红外光谱频率不随温度变化,因此有的研究者认为CH 组没有参与氢键作用[32],但随后其它作者认为CH 组红外光谱频率不随温度
变化并不能表明甘油中不存在C —H…O作用,C —H…O作用对角变形的不敏感性也能解释这一现象[6]。C —H…O是否存在于甘油水溶液体系中需要进一步的实验和理论研究,本文主要研究O —H…O形式的氢键。根据O-H 原子对的径向分布函数第一个波谷的距离dO-H ,H…O的特征长度Rc0O-H 选为0.24nm ,分子间氢键O —H…O的特征角度选为30°, 而分子内氢键的特征角度不进行约束
[21,27]。
2.3.2氢键结构分析
大部分的O 原子(大约90%)参与一个或者个氢键作用。在较低的甘油浓度下,50%的O 原子参与两个氢键作用,随着甘油浓度的增加,该比例有所下降, 当甘油浓度为3.83 mol·L-1时, 下降到不足40%。而参与单个氢键作用的O 原子的比例随着甘油浓度的增加逐渐增加,从最初的38.4%增加到46.6%。O 原子参与两个氢键作用的比例的下降与O 原子参与单个氢键作用比例的上升导致平均每个O 原子参与氢键的数目随甘油浓度的增加不断减少(从1.51减少到1.34) 。75%以上的H 原子参与了氢键作用,且这种比例基本不随甘油浓度的变化而变化,甘油浓度对平均每个H 原子参与的氢键数目的影响很微弱。
纯水和纯甘油体系的模拟结果与文献中的结果非常一致[23,27]。由于我们这里采用的H…O的特征长度Rc0O-H 比较小, 因此得到的平均每个分子参与的氢键数目相对较小。甘油分子由于含有3个O —H 键, 每个分子参与的氢键数目大约是每个水分子参与氢键数目的2倍。在较低的浓度下, 参与4、5、6、7、8个氢键的甘油分子比例均在10%以上, 随着浓度的增加, 甘油分子参与高数目氢键的比例有所降低, 当浓度增加到2.99 mol·L-1时,参与8个氢键的甘油分子的比例降到10%以下。除了纯甘油体系,参与6个氢键的甘油分子的比例最高, 大约在25%。大约有70%的水分子参与3或者4个氢键作用。随着甘油浓度的增加,参与4个氢键的水分子所占的比例逐渐下降,而参与3个氢键的水分子的比例逐渐增加。平均每个水分子参与的氢键数目随甘油浓度增加缓慢下降, 但是考虑到计算误差,这种下降趋势并不是很明显。随着甘油浓度的增加, 平均每个分子参与的氢键数目不断增加。
2.3.3氢键动力分析
由于氢键的快速破裂和重建, 溶液的氢键构型不断发生变化[27,33]。氢键的生
存周期τ对了解氢键的动力有着至关重要的作用, 人们采用不同的方法来研究t
[34-37]。我们采用自相关函数CHB(t)来研究氢键的生存周期[9,27,33,38],如果O 原子i 和H 原子j 在时刻0和t 均存在氢键作用,且破裂时间从未超过t*,则n ij (t)=1,否则,n ij (t)=0。其中,t*为特征时间,t*=0和t*=∞分别对应两种极限情况,即连续周期和截断周期[27]。这里我们只研究连续周期。但是由于模拟的时间步长永远不可能为0,因此不可能真正达到t*=0,本文的时间步长为2 fs,分子轨迹每两个步长记录一次,对应于t*=4 fs。但是数值模拟表明,氢键的长时间松弛行为并不是严格的指数关系,而是与扩散系数有关[39,40]。我们这里主要是研究氢键的生存周期随甘油浓度变化的趋势,并不需要精确的τ值,因此仍然采用指数关系来确定τ。水中氢键的生存周期为0.68 ps,与其它文献中的值非常接近[27,33,38]。甘油中氢键的生存周期为2.30 ps,与文献中的计算值3.0 ps相近[27]。虽然在浓度为2.99和3.34 mol·L-1时,t 、t w-w 和t g-w 均小于浓度为2.60 mol·L-1时的值, 但是氢键的生存周期t 、水分子之间相互氢键作用(w-w)的生存周期t w-w , 和甘油与水分子之间相互氢键作用(g-w)的生存周期t g-w 随甘油浓度的增加均呈现增加的趋势。在大部分甘油浓度下(box1除外), 水分子间相互氢键作用(w-w)的周期t w-w 大于甘油与水分子之间相互氢键作用(g-w)的周期t g-w 。
3、结论
利用分子动力学方法,模拟了不同浓度下的甘油和水的二元体系。分析了C-C ,C-O ,C-H ,O-O ,O-H 和H-H 原子对的分子间径向分布函数,并根据O-H 原子对的径向分布函数确定了氢键的图形准则。计算了参与i 个氢键作用的原子(O和H) 和分子(甘油和水) 的百分比分布,并根据其百分比分布计算了平均每个原子和分子参与的氢键个数。研究表明,平均每个O 原子参与的氢键个数随甘油浓度的增大而减小,但平均每个H 原子参与的氢键个数随甘油浓度变化并没有明显的变化。平均每个甘油和水分子参与的氢键个数随甘油浓度的增加呈现减小的趋势,但是考虑到误差, 这种趋势并不明显。如果不考虑分子的类别,则平均每分子参与的氢键个数随甘油浓度的增加而增加。同时,本文还报道了各个浓度下甘油分子的构型分布和氢键的生存周期。
参考文献
[1] Mazur, P.; Pinn,I.L.;Kleinhans,F.W.Cryobiology,2007,54(2):223
[2] Karlsson,J.O.M.;Cravalho,E.G.;Toner,M.J.Appl.Phys.,1994,75(9):4442
[3] Zhao,G.;Luo,D.;Gao,D.AIChE J.,2006,52:2596
[4]Chen,C.;Li,W.Z.Acta Phys.-Chim.Sin.,2008,24(1):74[陈聪, 李维仲. 物理化学学报,2008,24(1):74]
[5] Chen,C.;Li,W.Z.Cryoletters,2008,29(5):371
[6] Dashnau,J.L.;Nucci,N.V.;Sharp,K.A.;Vanderkooi,J.M.J.Phys.Chem.B,2006,110(27):13670
[7] Kuttel,M.;Brady,J.W.;Naidoo,K.J.J.Comput.Chem.,2002,23:1236
[8] Chelli,R.;Procacci,P.;Cardini,G.;Valle,R.G .D.;Califano,S.Phys.Chem.Chem.Phys.,1999,1:871
[9] Root,L.J.;Berne,B.J.J.Chem.Phys.,1997,107(11):4350
[10] Reiling,S.;Schlenkrich,M.;Brickmann,J.J.Comput.Chem.,
1996,17(4):450
[11] Blieck,J.;Affouard,F.;Bordat,P.;Lerbret,A.;Descamps,M.Chem.Phys.,2005,317:253
[12]Phillips,J.C.;Braun,R.;Wang,W.;Gumbar,J.;Tajkhorshid,E.;Villa,E.;Chipot,C.;Skeel,R.D.;Kale,L.;chulten,K.J.Comput.Chem.,2005,26:1781
[13] Berendsen,H.J.C.;Grigera,J.R.;Straatsma,T.P.J.Phys.Chem.,1987,91:6269
[14] Ryckaert,J.P.Mol.Phys.,1985,55(3):549
[15] Darden,T.;York,D.;Pedersen,L.J.Chem.Phys.,1993,98:10089
[16] Procacci,P.;Marchi,M.J.Chem.Phys.,1996,104:3003
[17] Martyna,G.J.;Tobias,D.J.;Klein,M.L.J.Chem.Phys.,1994,101:4177
[18] Feller,S.E.;Zhang,Y.;Pastor,R.W.;Brooks,B.R.J.Chem.Phys.,1995,103:4613
[19] Brunger,A.T.X-PLOR.The Howard Higher Medical Institute andDepartment of Molecular Biophysics and Biochemistry,YaleUniversity.p.A System for X.1992
[20] Hayward,J.A.;Reimers,J.R.J.Chem.Phys.,1997,106(4):1518
[21] Guardia,E.;Marti,J.;Garcia-Tarres,L.;Laria,D.J.Mol.Liq.,2005,117:63
[22] Jeffrey,G.A.Crystallography Reviews,2003,9(2-3):135
[23] Padro,J.A.;Saiz,L.;Guardia,E.J.Mol.Struct.,1997,416(1-3):243
[24] Chelli,R.;Procacci,P.;Cardini,G.;Califano,S.Phys.Chem.Chem.Phys.,1999,1:879
[25] Callam,C.S.;Singer,S.J.;Lowary,T.L.;Hadad,C.M.J.Am.Chem.Soc.,2001,123:11743
[26] Matsumoto,M.J.Chem.Phys.,2007,126:054503
[27] Guardia,E.;Marti,J.;Padro,J.A.;Saiz,L.;Komolkin,A.V.J.Mol.Liq.,2002,96-97:3
[28] Sarkar,S.;Joarder,R.N.Phys.Lett.A,1996,222(3):195
[29] Loof,H.D.;Nilsson,L.;Rigler,R.J.Am.Chem.Soc.,1992,114:4028
[30] Nag,A.;Chakraborty,D.;Chandra,A.J.Chem.Sci.,2008,120(1):71
[31] Jiang,L.;Lai,L.J.Biol.Chem.,2002,277(40):37732
[32] Zelent,B.;Nucci,N.V.;V anderkooi,J.M.J.Phys.Chem.A,2004,108(50):11141
[33] Elola,M.D.;Ladanyi,B.M.J.Chem.Phys.,2006,125:184506
[34] Rapaport, D.C.Mol.Phys., 1983, 50(5):1151
[35] Matsumoto,M.;Gubbins,K.E.J.Chem.Phys.,1990,93:1981
[36] Luzar,A.;Chandler,D.J.Chem.Phys.,1993,98:8160
[37] Ferrario,M.;Haughney,M.;McDonald,I.R.;Klein,M.L.J.Chem.Phys.,1990,93:5156
[38] Lee,H.S.;Tuckerman,M.E.J.Chem.Phys.,2007,126:164501
[39] Luzar,A.;Chandler,D.Phys.Rev.Lett.,1996,76(6):928
[40] Luzar,A.;Chandler,D.Nature,1996,379:55
相关文章
- 20**年执业兽医考试[动物生物化学]全真模拟题
- 食品化学重点
- 分子印迹技术及其应用
- 生化名词解释 (自动保存的)
- 纺织用酶制剂复配技术研究与进展
- 高三生物月考试题
- 氢键对吡啶与水络合物的影响
- 生物化学知识点总结
- 考研生物化学复习笔记.doc
- 生物化学试题库重点........
1.构成生物分子的元素有27种,其中约16种是所有生命所共有,下面说法正确的是 A.生物体中含量最丰富的元素是H.0.C.N四种元素,约占细胞重量的99% B生物体中含量最丰富的元素是H.O.C.S四种元素,约占细胞重量的99% C.生物体 ...
第一章 绪论 1.食品科学(food Science)--一门将基础学科和工程学的理论用于研究食品基本的物理. 化学和生物化学性质以及食品加工原理的学问. 它是一门涉及到食品的特性及其变化. 保藏和改性原理的科学 . 2.食品化学(Food ...
[摘要]分子印迹技术是在近年来发展起来的一项新兴技术.本文介绍了其基本原理.印迹聚合物的制备及分子印迹技术在膜制备和有机合成等方面的应用. [关键词]分子印迹技术:分子印迹聚合物:应用 [中图分类号]O658.9[文献标识码]A[文章编号] ...
第一章 糖类 α-及β-异头物(αor βanomer ):是指葡萄糖分子形成环状半缩醛结构以后,C1原子也变成了不对称碳原子,半缩醛羟基可产生两种不同的排列方位,因此形成了α-及β-两种异头物,α-型的羟基位于决定构型的羟基的同侧,β-型 ...
摘要:本文主要介绍了纺织用酶制剂复配技术的研究与进展.探讨了各类稳定剂.防腐剂.表面活性剂等的添加对酶制剂的影响,以及单酶复配与多酶复配技术等. 关键词:纺织用酶制剂:复配:定剂:防腐剂 Abstract:This paper introd ...
高三生物月考试题 考查范围:细胞的物质基础.结构基础及物质的输入和输出 注意事项: 1.本试卷分第I卷(选择题)和第卷(非选择题)两部分.答题前,考生务必将自己 的姓名.准考证号填写在答题卡上. 2.回答第I卷时,选出每小题答案后,用2B铅 ...
本科生毕业论文(设计) 中文题目 氢键对吡啶与水络合物的影响 英文题目 Hydrogen bonding with water of 学生姓名 李国徽 班级 物理学3班 学号11070328 学 院 吉林大学物理学院 专 业 指导教师 姜永 ...
生物化学复习题 第一章 绪论 1. 名词解释 生物化学: 生物化学指利用化学的原理和方法,从分子水平研究生物体的化学组成,及其在体内的代谢转变规律,从而阐明生命现象本质的一门科学.其研究内容包括①生物体的化学组成,生物分子的结构.性质及功能 ...
第一篇 生物大分子的结构与功能 第一章 氨基酸和蛋白质 一.组成蛋白质的20种氨基酸的分类 1.非极性氨基酸 包括:甘氨酸.丙氨酸.缬氨酸.亮氨酸.异亮氨酸.苯丙氨酸.脯氨酸 2.极性氨基酸:极性中性氨基酸:色氨酸.酪氨酸.丝氨酸.半胱氨酸 ...
生物化学试题库 蛋白质化学 答 案: 一.填空 1. 20 非极性 极性 疏水 亲水 赖 精 天 冬 2.色 苯丙 酪 3. -OH -SH -COOH N NH 4.氨基 紫红 亮黄 5.肽 氢键 二硫 键 疏水作用(键) 范德华力 二硫 ...