删除或更新信息,请邮件至freekaoyan#163.com(#换成@)

太赫兹波对钾离子通道蛋白二级结构影响的分子动力学模拟

本站小编 Free考研考试/2021-12-29

摘要:钾离子通道在神经细胞动作电位复极过程中起着重要作用. 钾离子通道蛋白种类繁多, 钾离子通道允许钾离子特异性穿过细胞膜, 从而维持神经细胞静息电位. 离子通道蛋白的二级结构决定其功能特性, 皮秒尺度内二级结构的波动会对离子通道蛋白的功能, 即离子通过速率有很大的影响. 本文使用分子动力学模拟方法, 模拟施加不同幅值的53.7 THz的太赫兹波对真实KcsA钾通道蛋白二级结构和钾离子通过速率的影响. 研究发现, 在53.7 THz的太赫兹波的作用下, KcsA钾通道蛋白中α螺旋数量减少, β折叠以及卷曲数量增加. 此外, 53.7 THz的太赫兹波能够加速钾离子通过KcsA钾通道. 本文从蛋白质的二级结构分析太赫兹波对钾离子通道蛋白的影响, 为太赫兹波和生物功能分子相互作用之间提供了新的观察角度.
关键词: 钾通道/
蛋白二级结构/
钾离子速率/
共振吸收

English Abstract


--> --> -->
钾离子通道广泛分布于神经、血管、血液等的细胞中, 在调节细胞, 尤其是神经细胞的膜电位和兴奋性中起重要作用[1-4]. 电压门控K+通道是典型的电压门控通道, 是最大的钾通道基因家族, 是神经元兴奋性的关键调节因子[5-7]. 人体中的电压门控钾离子通道包括40个不同的基因编码, 分为12个亚科, 从Kv1到Kv12. 哺乳动物电压门控钾通道是由围绕离子传导孔的4个相同α亚基组成的四聚体结构, 其中每个亚基有6个α螺旋跨膜结构域(S1—S6), 选择性过滤器由S5和S6之间的回路形成, 其中S1—S4序列对于通道电压感应和门控至关重要[5].
太赫兹波对研究生物系统的重要性在于其频率对应于脱氧核糖核酸(DNA)、核糖核酸(RNA)和蛋白质等生物功能分子的构象振荡范围, 以及其量子能量与氢键和范德瓦耳斯力等分子间相互作用的能量重合, 作用于生物分子上的太赫兹波和生物分子内部可能发生非线性共振, 导致其构象和功能发生剧烈变化, 从而对生物系统产生非热效应[8-13]. 因此, 太赫兹波有望具有广泛的潜在生物医学和安全相关应用. 已有研究表明, 太赫兹波和细胞内移动离子的相互作用主要是电作用[14]. 太赫兹波能够影响膜蛋白的功能从而对生理过程进行调制. 研究发现, 太赫兹波通过和离子通道蛋白某些基团形成共振, 影响离子通过速率, 进而影响生物体的记忆以及感知等功能[15-18]. 太赫兹波还可通过引起碱性磷酸酶的结构发生变化从而影响其活性[19]. 太赫兹波会引发各种细胞反应, 例如改变基因表达模式、改变基因组的结构和完整性以及影响细胞骨架的组装过程[20-22]. 近年来, 有研究报道了太赫兹波对细胞内遗传物质的影响. 太赫兹波能够通过与弱氢键相连的嘌呤振动之间形成共振, 破坏碱基对之间的氢键, 从而大大加速DNA的解螺旋速度[23]. 太赫兹波还可通过影响细胞内Zn2+离子的局部浓度变化来影响基因表达[24]. 上述研究表明, 具有特定频率的太赫兹波会对生物功能分子产生影响, 分布在细胞膜上的离子通道蛋白是太赫兹波的潜在分子靶标. 此外, 在模拟太赫兹波对生物功能分子的调制过程中, 大多使用仅包含主要基团的近似结构, 且对太赫兹波作用过程中离子通道蛋白二级结构的变化知之甚少. 本文采用分子动力学方法, 模拟太赫兹波对真实钾离子通道的作用, 并分析在此过程中钾离子通道蛋白二级结构的变化, 从而为太赫兹波和钾离子通道蛋白作用的动力学关系提供新的视角.
本文分析了太赫兹波对钾离子通道蛋白二级结构的影响, 计算不同幅值太赫兹波作用下钾离子通道蛋白氢键网络的变化, 以研究太赫兹波对蛋白质内部氢键的影响. 此外, 计算了不同幅值太赫兹波作用下钾离子通道蛋白各种二级结构的数目, 统计分析太赫兹波对各种二级结构的变化趋势的影响, 结果表明, 太赫兹波对蛋白质二级结构的影响主要体现在α螺旋、β折叠和卷曲这3种二级结构. 最后, 计算了不同幅值太赫兹波作用下, 钾离子通过钾离子通道的速率.
2
2.1.钾离子通道蛋白四聚体结构
-->从蛋白质结构数据库(RCSB PDB)中获得钾离子通道蛋白KcsA的单体结构, 并使用VMD将其转化成四聚体结构. 在计算前, 用CHARMMGUI网站补全原始结构中缺失的残基, 并用VMD检查补齐后的四聚体结构. 如图1(a)所示, KcsA通道结构包含1个胞内门, 将细胞质和充满水的空腔分开, 这个门是由4个TM2螺旋的羰基末端聚集在一起形成束交叉从而收缩通往孔道区的通道. 图1(b)是从正前方观察到的KcsA通道结构, 图中黄颜色为β折叠, 紫色为α螺旋, β折叠和α螺旋之间通过卷曲连接.
图 1 (a) KcsA四聚体结构(从上往下); (b) KcsA四聚体结构(从前往后); (c)钾离子通道模型
Figure1. (a) KcsA tetramer structure (from top to bottom); (b) KcsA tetramer structure (from front to back); (c) potassium channel model.

2
2.2.分子动力学模拟体系的构建
-->图1(c)所示, 将钾离子通道蛋白的四聚体结构正确嵌入二棕榈酰磷脂酰胆碱(DPPC)膜中, 并用SPC水填充整个盒子[25]. 为了模拟真实的生理环境, 在模拟体系中加入0.15 mmol的KCl, 对应于模拟体系中含有315个K+和331个Cl. 为了充分限制蛋白质和磷脂的位置, 在模拟之前, 先进行400 ps的预平衡. 之前的研究表明, 钾离子通道中负责钾离子特异性传输的基团是羰基, 羰基自发的振动频率为53.7 THz[16]. 因此, 在模拟盒子中应用53.7 THz的太赫兹波从上表面垂直入射, 模拟的时间步长为2 fs.
图1(a)图1(b)可以看到, KcsA过滤器有向内朝向孔中心的骨架羰基氧, 以协调通过的钾离子, 从而取代它们的溶剂化水. 已有研究表明, 钾离子通道的孔区域通过氢键网络固定在恰当的位置[26].
蛋白质内或蛋白质与溶剂分子之间的氢键网络的动态波动通过平均氢键对的相关函数[$ {C}_{hb}\left(\tau \right) $]来表征[27]:
$ {C}_{hb}\left(\tau \right) = \frac{\left\langle{h{b}_{i}\left(\tau \right) \cdot h{b}_{i}\left(0\right)}\right\rangle}{\big\langle{h{{b}_{i}(0)}^{2}}\big\rangle} , $
式中t时刻特定氢键i的存在函数是0或1 ($ hb (\tau) = $$ \left[\mathrm{0, 1}\right] $). 在分析中, 氢键是使用几何标准定义的, 即质心距离小于3.5 ?, O···H的距离小于2.6 ?, HO和O之间的角度小于30°. 蛋白质系统中其他弱相互作用(范德瓦耳斯力、静电等)被识别为近似截止距离内的接触.
二级结构, 如α螺旋和β折叠, 是以骨架氢键、骨架原子、Cα原子等连接[26]. 氢键和任何分子相互作用一样, 在短距离内很强并在无穷远处不断减小为零. 为了计算方便, 氢键被视为阶跃函数. 氢键的判定标准有能量准则、电子结构准则、几何准则等, 其中能量准则和电子结构准则的计算量大, 不易实现[27]. 在本文的计算中, 使用几何标准对氢键进行定义, 氢键是由氢原子-施体-受体所成角度的截断值与施体-受体之间距离的截断值共同决定的. OH和NH被作为氢键施体, O总是作为氢键受体, N默认为受体, 以便对KcsA四聚体中的氢键网络和各种二级结构进行区分.
钾离子通道蛋白质在脂质膜的复杂环境中起作用, 蛋白质原子不仅彼此相互作用并与水相互作用, 还和周围脂质分子的原子相互作用. 氢键可在蛋白质中氨基酸的残基内起作用, 是蛋白质三维结构稳定性和灵活性的驱动力之一, 对于膜蛋白的组装、结构和功能至关重要[28]. 天然膜环境中的膜蛋白不像在晶体结构中所观察到的是静态的, 相反, 它们是动态的. 理论和实验的研究结果都表明, 蛋白质二级结构中发生的皮秒时间尺度的波动对蛋白质的功能有很大影响[29-31]. 分子动力学模拟允许我们通过重建溶剂化膜蛋白的生理脂质膜环境扩展实验的观察范围, 从而研究原子细节的动力学[32]. 计算了不同幅值的太赫兹波辐射下钾离子通道蛋白氢键网络中氢键的数量, 结果如图2所示. 可以看出: 随着时间的增加, 氢键数量开始缓慢增加, 经过800 ps后, 氢键数量趋于稳定. 随着太赫兹波幅值从0增加到0.5 V/nm, 稳定后的氢键网络中的氢键数量在减少, 且随着幅值的增大, 减少趋势越来越明显, 下文也可以得出氢键数量的改变主要是由钾通道蛋白质中二级结构的数目变化引起的, 氢键的数量对蛋白质的稳定性具有重要作用, 氢键数量减少说明此时蛋白质的稳定性变弱, 钾离子通道蛋白的构象容易发生改变, 进而影响通道中钾离子通过速率.
图 2 太赫兹波强度不同时钾离子通道蛋白中氢键的变化
Figure2. Changes of hydrogen bonds in potassium channel proteins at different terahertz intensities.

二级结构允许简单描述蛋白质的复杂三维结构. 分子动力学软件Gromacs采用DSSP程序将通过X射线晶体学或核磁共振确定的原子坐标转换为供体和受体之间的主链氢键, 然后将重复键合模式转换为螺旋、片和匝等类别, 进一步可细分为α螺旋、卷曲、β折叠、β桥、弯曲、转角、回折等二级结构. 本文计算了在不同幅值太赫兹波作用下钾离子通道蛋白二级结构(总数 = α螺旋 + β折叠 + β桥 + 转角)的变化趋势. 如图3所示, 在经过近700 ps的稳定后, 钾离子通道蛋白的二级结构数目也趋于稳定, 从图3可以看出, 当外加太赫兹波时, 钾离子通道蛋白二级结构数目比不加太赫兹电磁波时多, 外加太赫兹电磁波会影响蛋白质二级结构的变化从而影响其功能. 并且在外加太赫兹波幅值为0.2 V/nm时, 钾离子通道蛋白二级结构的数目增加最多, 此时我们认为是钾离子通道蛋白和太赫兹波形成频率和幅值匹配才造成钾离子通道蛋白结构变化显著.
图 3 太赫兹波强度不同时钾离子通道蛋白中二级结构数目的变化
Figure3. Changes in the number of secondary structures in potassium channel proteins under different terahertz intensities.

蛋白质的二级结构又可进一步分为α螺旋、β折叠、卷曲和转角等结构. α螺旋是蛋白质肽链主链骨架借助氢键卷曲形成的一种周期性螺旋结构, α螺旋跨膜结构域的相互作用和动力学在膜蛋白的折叠、组装和功能中起着重要作用[33-35]. 氢键是α螺旋稳定的主要次级键, 当氢键破坏时, α螺旋构象即被破坏. 借助VMD, 我们发现α螺旋主要集中于S5和S6的孔道处, 即包围了钾离子通道最窄的地方. 不同幅值的太赫兹波作用下α螺旋的变化趋势如图4所示, 可以看出, 在体系平衡后, 外加太赫兹波时, 钾离子通道蛋白中α螺旋的数目要比未加太赫兹波时少, 并且外加太赫兹波的幅值增大时, 对应的α螺旋的数目也在缓慢减小. 因为太赫兹波的能量极低, 不能直接破坏共价键, 只能通过共振形式破坏氢键. 图4中的结果表明, 太赫兹波会破坏钾离子通道蛋白中α螺旋中的氢键, 使其数目减少.
图 4 太赫兹波强度不同时钾离子通道蛋白中α螺旋数目的变化
Figure4. Changes in the number of α-helices in potassium channel protein under different terahertz intensities.

β折叠是由两条以上肽链或同一条肽链之间的不同部分平行或反平行排列形成股, 股与股之间通过氢键固定, 氢键主要在股间. β折叠主要存在于S1—S4序列, 影响通道电压感应和门控. 不同幅值的太赫兹波作用下β折叠的变化趋势如图5所示, 可以看出, 外加太赫兹波时的β折叠数量比未加太赫兹波时多. 说明太赫兹波可能会促进不同肽链之间或同一条肽链之间不同部分形成新的β折叠, 也有可能是太赫兹波导致部分α螺旋展开形成β折叠. 从图4图5可知, 太赫兹波造成α螺旋减少, β折叠增多. β折叠往往含有较多的非极性残基, 并埋在蛋白质内部形成疏水核心, 而α螺旋通常是两性的, 亲水面位于表面, 疏水一侧朝向蛋白质内部. α螺旋的减少和β折叠的增多会导致钾离子通道中疏水核的暴露和亲水面的减少, 从而影响钾离子通过速率.
图 5 太赫兹波强度不同时钾离子通道蛋白中β折叠数目的变化
Figure5. Changes in the number of β-sheet in potassium channel protein under different terahertz intensities.

卷曲通过主链间的氢键或主链与侧链间的氢键稳定其构象, 其是柔性的, 可以使肽链改变走向, 且有利于连接结构相对刚性的α螺旋和β折叠. 从图6可以看出, 外加太赫兹波时钾离子通道蛋白中的卷曲数目减少, 并且在外加太赫兹波幅值为0.2 V/nm时, 卷曲的数目变化最大, 此时蛋白质的构象变化可能最大.
图 6 太赫兹波强度不同时钾离子通道蛋白中卷曲数目的变化
Figure6. Changes in the number of coil in potassium channel protein under different terahertz intensities.

钾离子通过钾通道的速率对于调节神经细胞的静息电位和动作电位至关重要, 并且钾离子通道对一些疾病起着调节作用, 是治疗的新靶点. 本文分别计算了无太赫兹电场、太赫兹电场强度为0.1, 0.2, 0.3, 0.4, 0.5 V/nm情况下钾离子通过电压门控钾离子通道的速率, 从图7可知, 最开始随着外加太赫兹波幅值的增大, 钾离子通过钾离子通道的速率也随之增大, 但是当幅值增加到0.3 V/nm及以上时, 钾离子通过的速率不再增加, 我们推测已经达到了钾离子通道的饱和, 这是由于孔道蛋白分子不是电中性的, 其对通过的钾离子具有一定的电吸引, 限定了它的最大运输速度. 并且从图4图6可知, 当外加太赫兹波的幅值为0.2 V/nm时, 对应的二级结构变化最大, 但是此时的钾离子通过速率却没有达到最大. 静电分析表明, 当钾离子沿着狭窄的孔道穿过膜时它必须穿过一个能量闸, 太赫兹波除了和钾离子通道蛋白中的羰基形成共振外, 还对模拟盒子中的钾离子有库仑力的作用, 两者共同作用使钾离子快速通过钾离子通道.
图 7 太赫兹波强度不同时钾离子通过钾离子通道的数目
Figure7. Number of potassium ions passing through potassium channels under different intensities of terahertz.

本文采用分子动力学的方法模拟了太赫兹波对真实钾离子通道蛋白二级结构的影响, 以及太赫兹波对钾离子通过钾离子通道的影响. 基于DSSP程序, 统计了在不同幅值太赫兹波的作用下, 钾离子通道蛋白主要二级结构皮秒时间尺度随时间的波动. 结果表明, 太赫兹波的作用下, α螺旋的数目减少, β折叠、卷曲等二级结构的数目有小幅增加. 特别是卷曲数目的改变表明了钾离子通道蛋白构象在太赫兹作用下发生了改变. 太赫兹波还能够促进钾离子通过钾离子通道. 开始阶段随着外加太赫兹波幅值的增加, 钾离子通过速率也跟着增加, 当增加到0.3 V/nm及以上时, 钾离子的通过速率不再变化, 这主要是由于已经达到了钾离子通道的饱和. 本文的研究为理解太赫兹波和离子通道蛋白集体振动的关系提供了一个新的视角.
相关话题/结构 计算 生物 网络 电压

  • 领限时大额优惠券,享本站正版考研考试资料!
    大额优惠券
    优惠券领取后72小时内有效,10万种最新考研考试考证类电子打印资料任你选。涵盖全国500余所院校考研专业课、200多种职业资格考试、1100多种经典教材,产品类型包含电子书、题库、全套资料以及视频,无论您是考研复习、考证刷题,还是考前冲刺等,不同类型的产品可满足您学习上的不同需求。 ...
    本站小编 Free壹佰分学习网 2022-09-19
  • 水的氢键网络动力学与其太赫兹频谱的关系
    摘要:水是万物生命之源,认识水的太赫兹吸收谱是太赫兹技术在生物医学上应用的前提,太赫兹频率的选择对高效、低能耗地实现太赫兹的生物效应至关重要.水的复杂氢键网络使得其具有较宽的太赫兹吸收峰,因此有必要研究水的氢键网络动力学与其太赫兹吸收谱之间的关系,然而这方面的研究仍然非常缺乏.采用分子动力学模拟方法 ...
    本站小编 Free考研考试 2021-12-29
  • COP-PDMS微流控芯片的制备及在太赫兹对肠道上皮细胞生物效应中的应用
    摘要:用于细胞培养的常规培养皿无法直接用于激光共聚焦显微镜的观察,在很大程度上限制了细胞水平上的太赫兹生物效应研究的开展.本文采用了对太赫兹吸收率低的材料—环烯烃聚合物(COP),与聚二甲基硅氧烷(PDMS)进行共键合,借用软蚀刻、光刻、等离子清洗和高温高压孵育技术将两种材料整合成新型微流控芯片,并 ...
    本站小编 Free考研考试 2021-12-29
  • 海洋上空折射率结构常数廓线估算
    摘要:湍流问题从提出到现在已困扰人们300多年.虽然提出了一些可行的方案,但在可预见的未来这一问题仍将困扰人们.湍流主要由浮力热气泡和风切变产生,在地球表面和大气之间传递物质和能量.2019年6月开展了第二次海洋季风实验(seamonsoonexperiment-II:SMEX-II),实验过程中通 ...
    本站小编 Free考研考试 2021-12-29
  • 二维<i>X</i>O<sub>2</sub> (<i>X</i> = Ni, Pd, Pt)弹性、电子结构和热导率
    摘要:基于第一性原理计算方法,研究了二维XO2(X=Ni,Pd,Pt)的稳定性、弹性、电子结构和热导率.计算结果显示,二维XO2同时具备较好的机械和动力学稳定性.此外,二维NiO2,PdO2和PtO2的杨氏模量分别为124.69N/m,103.31N/m和116.51N/m,泊松比分别为0.25,0 ...
    本站小编 Free考研考试 2021-12-29
  • 磁电电压可调电感器的理论设计与可调性优化
    摘要:磁电电压可调电感器(ME-VTI)是基于磁电效应实现电场对电感的调控.与其他可调电感器相比,具有能耗低、体积小、可调性大且连续等特点.而以往对ME-VTI的研究主要针对结构和磁致伸缩材料,致使电感器制备工艺复杂或可调性提高不大.本文通过构建理论模型,着眼于压电材料场致应变对电感可调性的影响.采 ...
    本站小编 Free考研考试 2021-12-29
  • ENSO气温关联网络结构特征差异及成因分析
    摘要:基于全球地表气温资料,分别构建了ElNi?o和LaNi?a事件对应的全球气温关联网络,并分析网络结构特征的差异,探究可能的成因.结果表明,与LaNi?a事件气温网络(简称LaNi?a网络)相比,ElNi?o事件气温网络(简称ElNi?o网络)中格点温度序列间的相关性减弱,气温网络的连通度显著降 ...
    本站小编 Free考研考试 2021-12-29
  • 电子束对ZnO和TiO<sub>2</sub>辐照损伤的模拟计算
    摘要:电子辐照在材料中产生的缺陷主要是相互独立的空位-间隙原子对,由于不同靶原子的离位阈能不同,通过改变电子束的能量可以调控在材料中产生的缺陷类型,同时,电子的注量又可以决定电子辐照产生的缺陷的浓度.ZnO和TiO2的磁光电特性受Zn空位、Ti空位、O空位、Zn间隙原子、Ti间隙原子等缺陷的影响,因 ...
    本站小编 Free考研考试 2021-12-29
  • 多孔结构体材料热整流效应
    摘要:基于傅里叶导热定律,两种具有不同热导率温度依赖特性的材料组合而成的两段式组合材料可以实现热整流效应.本文提出在体材料上均匀布置多孔结构,通过多孔结构孔隙率调整材料的热导率参数,进而强化热整流效应.基于有限元方法和有效介质理论,计算并分析了温差和孔隙率等参数对体材料热整流系数的影响.计算结果表明 ...
    本站小编 Free考研考试 2021-12-29
  • 腔体结构参数对毛细管放电型脉冲等离子体推力器放电特性的影响
    摘要:毛细管放电型脉冲等离子体推力器在微纳卫星在轨机动应用中具备良好的发展前景.本文系统研究了单次放电能量5J条件下,不同毛细管腔体尺寸参数对推力器能量沉积特性、烧蚀特性、输出推力参数和等离子体羽流参数的影响规律.实验结果表明,增大毛细管腔体内径会显著降低放电电流密度,减小沉积能量和等效功率;增大腔 ...
    本站小编 Free考研考试 2021-12-29
  • 基于太赫兹石墨烯等离激元的多参数相位可调谐结构及其应用
    摘要:太赫兹波的振幅和相位进行主动调控由于在太赫兹功能器件方面的广泛应用,受到了广泛关注.目前采用的金属-介质-金属超表面结构结合石墨烯等二维材料可实现太赫兹振幅/相位的动态调控,但存在调控自由度少(电压或光强)以及超表面结构加工工艺复杂及价格昂贵等缺点.本文提出了一种棱镜耦合石墨烯等离激元结构的相 ...
    本站小编 Free考研考试 2021-12-29