论文部分内容阅读
本文采用分子动力学(MD)方法对三种管径分别为6、8和10 A的跨膜环肽纳米管8xcyclo(WL)n=3,4,5-PNT/POPE,采用平衡及非平衡分子动力学方法模拟研究水分子在此类跨膜纳米通道中的结构特点和传输特性(包括水链结构、偶极取向、密度分布和渗透与扩散特性等),主要研究成果如下:
1.采用模拟退火和MD模拟的构象搜索方法获得了三种不同管径的环肽纳米管8xcyclo(WL)n=3,4,s-PNT的最低能量构象。考察了模拟过程中环肽纳米管骨架结构的稳定性。完全水化的POPE脂膜经过预平衡后,获得了每个脂分子的面积和有效厚度基本保持不变的POPE脂膜平衡构型。将环肽纳米管的平衡构型插入平衡的水--脂膜中,预平衡后获得了相应的平衡构型。
2.对跨膜环八肽纳米管8xcyclo(WL)4-PNT/POPE体系进行50 ns的MD模拟,详细考察了纳米管内水链的结构及其形成机制。水分子在该纳米管内存在波状模式的水链,以“1-2-1-2”方式排列,即每个α-plane区域有一个水分子,每个mid-plane区域有两个水分子,波长恰好为两相邻环肽间的距离4.8 A。
(1)管内水分子的柱形分布函数曲线和PMF(Potential of Mean Force)计算说明,跨膜环八肽纳米管中形成“1-2-1-2”水链模式的主要原因是该纳米管的空间限制作用。管内水分子的径向分布函数(O-O和H-H)合理解释了环八肽纳米管中“1-2-1-2”水链结构。管内mid-plane区域的每个水分子与骨架间形成强氢键的几率不到0.1,表明环八肽纳米管与水分子间的氢键作用不是影响管内“1-2-1-2”水链结构形成的主要因素。环八肽纳米管中水分子形成的氢键网络模式为:在α-plane区域的一个水分子与相邻mid-plane区域的两个水分子形成氢键,mid-plane区域的每个水分子与相邻α-plane区域的水分子和纳米管骨架C=O形成氢键。
(2)纳米管内水分子的偶极取向和水分子沿管轴的势能曲线表明,水分子的偶极翻转发生在第三和第四个环之间,来自于环肽纳米管、脂膜和本体水分子的合力对管内靠近管口的水分子(|z|>10 A)有很强的定位效应;在管中间部分的水分子(|z|<10 A),外界的影响不够强,不足以克服热运动。
3.采用平衡MD和非平衡MD模拟,分析了水分子在三种跨膜环肽纳米管8×cyclo(WL)n=3,4,5-PNT/POPE中的传输性质,得到三种纳米管中水分子传输的特点和特性参数diffusion permeability(Pd)osmotic Prmeability(Pf).
(1)跨膜环六肽纳米管中,水分子呈单链排列,大多数水分子占据着mid-Plane区域的位置;跨膜环八肽纳米管中,水分子为“1-2-1-2”模式排列;纳米管中,水分子的排列近似为本体水中的结构。
(2)水分子在三种跨膜环肽纳米管8xcyclo(WL)n=3,4,5-PNT/POPE中运动是一维扩散运动,沿纳米管管轴方向运动,满足菲克运动。
(3)采用平衡MD模拟首次获得了水分子在三种跨膜环肽纳米管8xcyclo(WL)n=3,4,5-PNT/POPE中的扩散特性参数diffusion permeability(Pd)分别为0.488、21.17和94.08 A3/ns。采用平衡和非平衡MD模拟两种方法首次获得了水分子在上述三种纳米通道中的渗透特性参数osmotic permeability(Pf)分别为3.97、66.15和496.86 A3/ns(平衡MD),2.30、71.85和473.31 A3/ns(非平衡MD)。本文模拟获得的Pf值与AQP1的pf实测值有一定的可比性。比较了三种环肽纳米管的Pf/Pd值与+1(为管内水分子的平均数目)的关系,发现只有环六肽纳米管Pf/Pd值近似等于+1,说明水分子在管中协同运动,没有交换和穿越,符合单链传输的特点。随着纳米管管径变大,水链不再以单链形式存在,(+1)/(Pf/Pd)值随着管径变大。
综上,本文选定具有合适疏水性侧链的环肽纳米管8xcyclo(WL)n=3,4,5-PNT嵌入到完全水化的POPE脂膜中组成的系统作为生物水通道的简化模型,获得的水链结构及其形成机制的分析和水分子在三种跨膜通道中的扩散和渗透特性的分析,将有助于在原子水平揭示纳米尺度水通道的本质特征。