基于构象空间离散化的复杂体系自由能计算方法开发

来源 :吉林大学 | 被引量 : 2次 | 上传用户:lconan
下载到本地 , 更方便阅读
声明 : 本文档内容版权归属内容提供方 , 如果您对本文有版权争议 , 可与客服联系进行内容授权或下架
论文部分内容阅读
研究生物大分子结构和功能的关系在生物和化学领域是非常重要的课题。在对蛋白质分子的研究过程中,无论是对蛋白质折叠过程分析、蛋白质与小分子的对接,还是蛋白质结构预测等过程都离不开对蛋白质构象转化的分析。明确构象转化过程中的自由能的变化是研究蛋白质分子构象转化的核心问题。由于生物大分子结构有着非常高的复杂度,因此对生物大分子高阶相关性的研究方法并不多,而且在复杂体系计算自由能时,现有方法效率和准确度不能同时满足。因此,从一个高维整体角度来分析蛋白质的构象空间以及开发新的自由能算法都具有非常重要的价值。本篇论文首先通过分层高维聚类方法来研究了溶菌酶蛋白质的构象转化,然后提出了基于构象空间的离散化的一种新型的复杂体系计算自由能方法,并且从理论推导和分子模拟数据中证明了该方法的可行性,最后我们还对比了分子模拟轨迹中能量与自由能在蛋白质不同二级结构上的相关性。自由能地貌图被广泛用于研究自然折叠状态下的蛋白质结构与功能的关系。蛋白质等生物大分子的复杂度非常高,其构象空间十分复杂,经过分子动力学模拟等采样方法,一般我们可以得到较为广泛的样本构象空间,然而,人们普遍从一维或者低维的角度来勾勒蛋白质的自由能地貌,以规避高阶相关性带来的困难。这也使得我们得到的蛋白质构象变化的结果局限在低维分析中,缺少在高阶相关方面上的分析,因此我们期望在一个高维的角度去描述自由能地貌,从而了解生物大分子的构象变化。在本文中,我们通过毫秒级别的分子动力学模拟方法研究了溶菌酶蛋白质在自然折叠状态下的构象空间,通过选择主链二面角为自由度,依据转换频率,区分自由度的快与慢,在与时间相关的不同层次从高维的角度研究了溶菌酶构象变化。经过与实验数据对比,晶体结构的构型在采样轨迹的分组中为主导构型,而且有着较高的能垒,相对难于向其他构型转化,并且蛋白质的结构和功能与非作用位点二面角的构象转化息息相关。自由能计算在生物化学等领域是一个亟待解决的非常重要的问题。在自由能计算的先河中,主要分为两类计算方法,一种是依据实验数据采取的打分函数;一种是通过物理模型与公式推导开发的方法。由于生物大分子的复杂度,第一种方法,需要大量的实验数据,且难于对未知体系进行预测,因此利用严格的公式推导更适合应用到更多的体系中。主要的理论计算自由能的方法分为两种模式,第一种为需要对两种状态之间路径采样的热力学积分、自由能微扰等方法;第二种为只需对两个状态无需对过程采样的MM/P(G)BSA、LIE等方法。第一类方法理论推导严密,能够准确的计算不同状态之间的自由能差,局限性是对于复杂的生物大分子体系,其操作过程较为复杂而且需要大量的采样,会消耗大量的计算机资源。第二种方法,依据各状态分子结构特点,通过对模型的近似求解两个状态之间的自由能差,能够较快的计算出结果,相对节省计算机资源,但是,由于这几种方法并没有解决体系中熵焓补偿的作用,对熵的计算非常不准确,导致计算得到的自由能误差很大。因此,开发新的方法还是非常有必要的。本篇论文中我们提出了一种基于显性恒定统计权重分布(ECISWD)假设,将构象空间离散化的在复杂体系中计算自由能的新方法。当对构象空间离散化足够精细的时候,不同宏观状态之间的自由能变可以通过对不同宏观状态下的conformer计数得到。另外,我们发现依据此时足够精细的conformer标准作为微观状态时,构象熵的变化和自由能的变化仅相差一个负的温度因子,不同宏观状态焓变和宏观状态内conformer内部的熵的变化近似相等。因此这种方法,很好地解决了熵焓补偿带来的误差,可以较为准确的预测自由能差。经过理论推导,这个方法在简单体系下是不成立的,对于复杂体系才适用,并且我们在磷脂分子和蛋白质分子分子动力学模拟轨迹中加以验证,证明了这种方法是一种非常有效的计算不同状态之间的自由能差的方法。最后,我们分析了分子动力学模拟轨迹中宏观状态之间不同能量项与自由能的关系,提出了能量跨度在计算自由能差中的重要地位,并且对比了多种蛋白质不同二级结构上的不同能量项与自由能的相关性差异,我们发现,在稳定的二级结构中,自由能与最小能量、最大能量以及能量跨度的相关性更高,其中能量跨度是相关性最高的参数。
其他文献
金属硫族化合物具有诸多优异的性能,在众多领域具有良好的应用前景,而成为人们研究的热点。利用其优越的热电性质,制备出热电器件己得到广泛应用;该族化合物具有良好的拓扑绝缘性,成为实现量子反常霍尔效应以及寻找马约拉纳费米子的候选材料,而得到广泛关注;此外,在压力作用下,金属硫族化合物呈现出丰富的压致结构(晶体结构与电子结构)变化信息,并带来全新的物理性质。如高压下金属硫族化合物发生的电子拓扑相变导致热电
考虑到理论上研究蛋白质生物化学过程中的动力学和热力学对理解其生物功能和本质及其重要,而实验上所能给出的测量通常是宏观的热力学或者动力学量,缺乏在分子层次上理解生命附加活动的过程。本文主要研究兴趣是:以分子动力学模拟为基础手段,建立起一套可以捕捉到某些由于实验精度限制而不能观测到的蛋白质生物学过程的模拟方法,在分子层次上去理解这些生物学过程。从而实现理论和实验桥梁的作用,深化对生物蛋白分子功能的理解
现代教育发展中,人文素养教育是重要内容,人文素养主要是人具有的人文精神素养。人们在阅读过程中,一方面能获取信息知识,同时也是读者和作者心灵交流的过程,因此,阅读对提升人们的人文素养有积极的意义。现当代文学作品内容丰富,与时代发展贴近,文字通俗易懂,且其中蕴含着丰富的人文精神和价值理念。基于阅读视角下,从多个维度对现当代文学作品中所蕴含的人文素养进行分析,对培养和提升学生的人文精神来说意义深远。
期刊
纳米硫化镉的同质异相转变引起电子结构的改变,影响相应的激子复合动力学和自旋极化机制。激子复合动力学可以探测激子寿命和对应的弛豫或者复合路径,进而理解其相应的电子结构。在室温下制备具有高自旋极化率的材料,为发展自旋半导体器件奠定了基础。尺寸的变化会影响激子的束缚状态,能带结构和材料的磁畴结构,影响着激子动力学和自旋极化率。掺杂可以改变半导体中载流子的种类并引入杂质能级进而改变激子复合的路径和自旋极化
随着晶体X射线衍射和核磁共振(NMR)等技术的发展,越来越多的蛋白质结构得到解析。截止到2017年,蛋白质(PDB)数据库中已存有约120000个蛋白质晶体结构。但蛋白质的静态结构通常不能为我们研究蛋白质功能提供足够的信息。要将蛋白质的结构与功能联系起来就需要我们研究蛋白质的运动和结构变化。计算机技术和分子动力学模拟方法的发展,为我们探索蛋白质分子的微观运动提供了一条可行途径。随着计算机速度的提高
高压是一种清洁的、强大的实验方法,是可调节的物理参量。高压能够非常有效地缩短原子间距离、增加相邻电子轨道重叠,进而改变原子(分子)间的相互作用和电子结构,形成常规条件下难以形成的具有新结构与新性质的高压新相。氢元素(H)是相对原子质量最小的元素,因此在化学元素周期表中排于第一位,由于其简单的电子结构也被称作是最简单的元素,同时也是宇宙中分布最广的元素。氢元素独特的电子排布及成键特征,使其既可归类为
目的 探讨加巴喷丁联合背根神经节脉冲射频治疗带状疱疹后神经痛(PHN)的临床疗效及安全性。方法 收集2019年1月至2021年1月河北北方学院附属第一医院疼痛科收治的68例PHN患者的临床资料,按照治疗方法不同分为对照组和研究组,每组34例。对照组采用常规治疗+背根神经节脉冲射频治疗,研究组采用常规治疗+背根神经节脉冲射频+加巴喷丁胶囊(第1天每次0.3 g,每日1次;第2~3天每次0.3 g,每
稀土离子具有独特的电子层结构和丰富的能级,使其被广泛地用于荧光材料的研究。随着上世纪末纳米科学技术研究的兴起,稀土掺杂的纳米材料逐渐引起了科研工作者们的兴趣。近年来,高温热分解法被广泛地用于制备纳米材料。相比于其它制备方法,这种方法的优势是制备周期短并且容易获得更小尺寸的纳米粒子。另外,基质对材料的发光性能有很大的影响。为了减少稀土发光材料中非辐射跃迁造成的能量损失以及最大化辐射能量跃迁,基质材料
主持人语核心素养时代呼唤大单元教学和跨学科融合教学。《普通高中课程方案(2017年版2020年修订)》在前言中提出,要“重视以学科大概念为核心,使课程内容结构化,以主题为引领,使课程内容情境化,促进学科核心素养的落实”。人们常说“政史不分家”,思想政治学科强调与其他学科的相互配合。面对当前社会变革和实践创新中的新挑战、新问题,一方面,教师要用历史唯物主义的眼光引领学生通过观察、辨析、反思和实
期刊