您所在位置： &
&nbsp&&nbsp
ELPs相变及其与木聚糖酶相互作用的分子动力学模拟.pdf82页
本文档一共被下载：
次 ,您可免费全文在线阅读后下载本文档
文档加载中...广告还剩秒
需要金币：50 &&
你可能关注的文档：
··········
··········
学校代码： 10385 分类号： 研究生学号： 密 级： 硕士学位论文 ELPs 相变及其与木聚糖酶相互作用的分子动力学模拟 The molecular dynamics simulation of ELPs phase transition and its interaction with xylanase 作者姓名： 黄可君 指导教师： 张光亚 教授 学
科：生物化学与分子生物学 研究方向： 生物信息学 所在学院： 化工学院 论文提交日期：二零一四年六月二日
万方数据 学位论文独创性声明 本人声明兹呈交的学位论文是本人在导师指导下完成的研究成 果。论文写作中不包含其他人已经发表或撰写过的研究内容，如参考 他人或集体的科研成果，均在论文中以明确的方式说明。本人依法享 有和承担由此论文所产生的权利和责任。 论文作者签名： 签名日期： 学位论文版权使用授权声明 本人同意授权华侨大学有权保留并向国家机关或机构送交学位论 文的复印件和电子版，允许学位论文被查阅和借阅。本人授权华侨大 学可以将本学位论文的全部内容或部分内容编入有关数据库进行检 索，可以采用影印、缩印或扫描等复制手段保存和汇编本学位论文。 论文作者签名： 指导教师签名： 签 名 日期： 签 名 日期： II
万方数据 摘要 摘 要 近年来，类弹性蛋白多肽 （Elastin-like polypeptides ，ELPs ）的可逆相变特 性逐渐开始受到学术界的关注，基于此，可逆相转换循环 （Inverse transition cycling,ITC ）的技术已经成为了一种高效低成本的分离ELPs
的方法。实验室设 计并表达了一种全新的ELPs 作为纯化标签，并通过可逆相转换循环完成了对包 括木聚糖酶等目标蛋白在内的分离纯化。这种ELPs 就是ELPs[KV F-20] 。但实 8 验中存在如下几个问题：1、所有的 ELPs
对不同盐（特别是阴离子）的响应存 在明显差异，但其机制不明
正在加载中，请稍后...您所在位置： &
&nbsp&&nbsp&nbsp&&nbsp
糖苷水解酶Fnce15A和GB1多肽分子动力学模拟.pdf59页
本文档一共被下载：
次 ,您可免费全文在线阅读后下载本文档
文档加载中...广告还剩秒
需要金币：200 &&
你可能关注的文档：
··········
··········
University
Technology
Adissertationformaster，S
glycosidehydrolase
Author’SName：
speciality：
Bioinformatics
Supervisor：Prof．Haiyan
Finishedtime：
中国科学技术大学学位论文原创性声明
本人声明所呈交的学位论文，是本人在导师指导下进行研究工作所取得的
成果。除已特别加以标注和致谢的地方外，论文中不包含任何他人已经发表或
撰写过的研究成果。与我一同工作的同志对本研究所做的贡献均已在论文中作
了明确的说明。
作者签名：
签字日期：―皇型瓣t／
中国科学技术大学学位论文授权使用声明
作为申请学位的条件之一，学位论文著作权拥有者授权中国科学技术大学
拥有学位论文的部分使用权，即：学校有权按有关规定向国家有关部门或机构
送交论文的复印件和电子版，允许论文被查阅和借阅，可以将学位论文编入《中
国学位论文全文数据库》等有关数据库进行检索，可以采用影印、缩印或扫描
等复制手段保存、汇编学位论文。本人提交的电子文档的内容和纸质论文的内
容相一致。
保密的学位论文
正在加载中，请稍后...您所在位置： &
&nbsp&&nbsp&nbsp&&nbsp
烟碱样乙酰胆碱受体门控机制的动力学研究.pdf144页
本文档一共被下载：
次 ,您可全文免费在线阅读后下载本文档。
文档加载中...广告还剩秒
需要金币：150 &&
烟碱样乙酰胆碱受体门控机制的动力学研究.pdf
你可能关注的文档：
··········
··········
大连理工大学 博士学位论文
烟碱样乙酰胆碱受体门控机制的动力学研究 姓名：刘信力 申请学位级别：博士 专业：工程力学 指导教师：王希诚;蒋华良 座机电话号码 大连理工大学博士学位论文 摘 要 自从20世纪70年代，分子动力学模拟首次在蛋白质分子应用以来，尤其在生物科
学领域变得越来越普遍。它提供了一种用计算机在原子水平上动态地研究分子内和分子 间反应的途径。直到今天，不仅是模拟体系的规模一直在快速增长，而且分子动力学方 molecular
法本身也一直在发展。例如，拉伸分子动力学 steered
分子动力学的一个分支，20世纪90年代以来在研究配体的结合／解离过程中起到了重要 的作用。在本文中共有三部分的工作。首先，用长时间 30纳秒 的常规分子动力学模
拟了探索烟碱样乙酰胆碱受体 nicotinicacetylcholine rotationmolecular
且我们发展了一种新的方法一一拉伸转动分子动力学 steered
作中揭示的烟碱样乙酰胆碱受体 hAChR 的细胞外结构区域转动的成因，我们使用自 己设计的三个拉伸模型，分别使用拉伸分子动力学在可以作为烟碱样乙酰胆碱受体 binding
应用在解离细胞色素P4503A4-甲吡酮的复合物上，比较结果显示优于传统的拉伸分子
动力学方法。 论文由五个部分组成，概述如下： ionchannel 的发 第一章概述了分子动力学理论和配体门控离子通道 1igand．gated
展现状和基本概念。此外，也简要地陈述了本文的主要工作和发现。 第二章作为预备知识，介绍了分子动力学的理论和应用。首先是回顾了分子动力学
的简史，并提到了牛顿方程、原子间势函数和有限差分算法。拉伸分子动力学和自由能
对计算也包括在这一章中。接着总结了使用和发展分子动力学的原因，并
正在加载中，请稍后... 上传我的文档
 下载
 收藏
该文档贡献者很忙，什么也没留下。
 下载此文档
正在努力加载中...
分子动力学研究2-乙酸苯并噻吩在HIV-1蛋白酶与抑制剂结合中的作用
下载积分：216
内容提示：分子动力学研究2-乙酸苯并噻吩在HIV-1蛋白酶与抑制剂结合中的作用
文档格式：PDF|
浏览次数：3|
上传日期： 16:17:52|
文档星级：
全文阅读已结束，如果下载本文需要使用
 216 积分
下载此文档
该用户还上传了这些文档
分子动力学研究2-乙酸苯并噻吩在HIV-1蛋白酶与抑制剂
官方公共微信AMBER分子动力学简例(一)
09:55&3384人阅读&(0)&&
&&&&&&&&&&&&&&&&&&&
AMBER分子动力学简例（一）
以下是使用AMBER包的简单教程，希望对开始学习分子动力学的同学有用处。申明一下，以下教程原版来自网上，是最最基本的教程，同时也非常实用，有非常好的借鉴意义。
AMBER分子动力学程序包是加州圣弗兰西斯科大学（University of California San
Francisco,UCSF)的Peter
A Kollman和其同事编的，程序很全，现在已经发展到版本9.0。AMBER功能涵盖种类非常多的生物分子，包括蛋白、核算以及药物小分子。软件详细情况请浏览amber.scripps.edu.
以下是AMBER软件包中四个主要的大程序：
--------------------------------------------------------------------------------
Leap：用于准备分子系统坐标和参数文件，有两个程序：
xleap：X－windows版本的leap，带GUI图形界面。
tleap：文本界面的Leap。
Antechamber：用于生成少见小分子力学参数文件的。有的时候一些小分子Leap程序不认识，
&&&&&&&&&&&&&&&&&&&&&&&&需要加载其力学参数，这些力学参数文件就要antechamber生成。
Sander：MD数据产生程序，即MD模拟程序，被称做AMBER的大脑程序。
Ptraj：MD模拟轨迹分析程序。
--------------------------------------------------------------------------------
本教程研究的题目是脑下垂体荷尔蒙之一的oxytocin，需要X光衍射晶体结构文件1NPO.PDB。该文件包含了该荷尔蒙和其运载蛋白的复合物，可以从蛋白数据库下载。
PDB文件是不包含氢原子的，Leap程序会自动的加上PDB文件缺少的东西。当第一次使用PDB文件的时候，要十分留意文件包含的信息，所以PDB文件缺少的残疾、侧链或者添加的变异都在这个地方记录。可以用文本阅读程序阅读PDB文件头部的信息，即以REMARKS开头的信息文本行。PDB一个重要的信息是SSBOND记录，该记录说明结构中二硫键的位置，这样的信息在使用Leap程序建立分子结构的时候需要。在本教程中，我们将比较oxytocin在真空和溶液中分子动力学的差异，如果没有二硫键，将影响整个结果。整个过程在Linux系统下完成，如对Linux不熟悉，请翻阅有关Linux书籍。
建立项目目录，并进入该目录：
下载1NPO.PDB文件到该目录下，使用文本阅读器阅读该PDB文件。从在文件的开头部分可以得到该文件是一个二聚物的PDB文件，删除文件中A、C、和D链对应的文本行，剩下的就是oxytocin的晶体机构了，保存为oxyt.pdb文件。
使用Swiss PDB viewer分析以下oxyt.pdb文件，可以得知该pdb文件缺少了一个侧链。如果使用Deep View，它会自动给文件加上侧链。重新保存pdb文件为oxyt.pdb文件。
xleap和tleap程序功能是一样的，都是准备分子结构的坐标文件和拓扑文件。xleap启动一个X界面，比较慢。tleap纯文本，要快一些，我们选择tleap：
leaprc.ff03
其中leaprc.ff03是AMBER的2003力场文件。力场是一个很重要的文件，定义了分子、原子和残疾等的信息，请查阅有关资料。
oxy=loadpdb
该命令读入oxyt.pdb文件到oxy变量中。
oxy.1.SG& oxy.6.SG
连接oxy变量的第一和第六个残疾的SG原子，即是连接二硫键。可以使用check命令检查oxy是否完好，没有特殊情况的话，最后应该是“OK”，表示分子系统状态是好的，可以保存。
接下来就是保存分子系统了：
saveamberparm
oxy_vac.top
oxy_vac.crd
其中oxy_vac.top和oxy_vac.crd分别是分子系统真空状态下的拓扑和坐标文件。
接下来要给分子系统添加水环境，
solvateoct
该命令让Leap程序使用TIP3PBOX水模型,水环境距离分子为9纳米。也可以用solvatebox命令，但是solvatebox添加的水环境是一个立方体，不是八面体。一般要求水表面到蛋白距离要达到8.5纳米，避免MD过程中蛋白冲出水环境，但是水环境越大计算时间将越长。如果MD过程内含PME，那么水箱长度要大于2倍截矩（cutoff），截矩在MD配置文件中定义。加溶剂之后，要计算系统是否为电中性，使用命令：
charge oxy
如果现实不是电中性，要使用addions添加反性电荷，如Cl－或者Na＋等。
最后保存溶剂环境的系统拓扑结构和坐标文件，并退出Leap程序：
saveamberparm
也可以把命令集中到一个，让tleap读取。如将上面的命令集中到以下文件中
"oxy.leaprc"
--------------------------------------------------------------------------------
leaprc.ff03
oxy=loadpdb
bond& oxy.1.SG
saveamberparm oxy oxy_vac.top
oxy_vac.crd
solvateoct
saveamberparm oxy
oxy.top& oxy.crd
--------------------------------------------------------------------------------
将以上两横线之间的命令保存为"oxy.leaprc"，然后使用以下命令一步搞定：
oxy.leaprc
到此，分子系统准备已经完成，目录中会多了一个leap.log文件，这是Leap程序的记录文件，它包含了leap程序所做的一切，包括自动的和手动命令。
现在已经有了可以进行分子动力学的基本文件了，再编写sander的控制文件，就可以进行模拟了。这些将在后文中介绍。&&&&
AMBER分子动力学简例（二）
分子动力学（1）
真空模式分子动力学模拟将使用NVT系宗分两步进行，即系统能量最优化和分子动力学过程。
1、系统能量最优化。我们将使用淬火能量最优化解除系统内原子之间的不正常相互作用，这些原子之间的高能量相互作用如果不消除，可能影响后续的分子动力学过程。因为动力学过程是能量梯度变化的，太高的能量壁垒可能让MD局限在某一个能量局部最小化位置中。
Sander程序是分子动力学模拟程序，它的主要功能是能量最优化，动力学模拟和NMR优化计算。我们必须给Sander一个运行的配置文件，使其按照我们的要求进行计算。能量最优化的配置文件如下：
"min_vac.in"
--------------------------------------------------------------------------------
oxytocin:& initial minimization prior to MD
maxcyc =& 500,
ncyc = 250,
--------------------------------------------------------------------------------
Sander程序的基本命令参数如下：
sander& &O &i in &o out &p prmtop &c inpcrd &r
restrt [-ref refc &x mdcrd &v mdvel &e& mden &inf
所以我们的命令可以如下：
sander &O &i& min_vac.in &o min_vac.out &p
oxy_vac.top &c oxy_vac.crd &r oxy_vacmin.rst&
可以使用more命令或者tail命令查看min_vac.out的输出内容，那是能量最优化的记录。
2、分子动力学模拟。
Sander程序的配置文件为：
--------------------------------------------------------------------------------
oxytocin& MD in-vacuo, 12 angstrom cut off, 250
imin = 0, ntb =& 0,
igb = 0, ntpr = 100, ntwx = 500,
ntt = 3, gamma_ln = 1.0,
tempi =& 300.0, temp0 = 300.0,
nstlim = 125000, dt = 0.002,
cut =& 12.0
--------------------------------------------------------------------------------
sander& &O &i md_vac.in &o md_vac.out &p
oxy_vac.top &c oxy_vacmin.rst &r oxy_vacmd.rst& &x
oxy_vacmd.mdcrd &ref oxy_vacmin.rst && &
MD的计算过程一般比较久，真空相对与溶剂中要快以下，250ps的模拟大概在一个主频为2。0GHz的Linux单机上运行5分钟。
以上配置文件中用到很多参数，这些参数这是sander程序参数的一小部分，以下将相应解释。如果要了解sander的其他参数，请阅读AMBER用户指南。
&ctrl和"/"：sander的参数一般要求出现在这两个标识符号之间，参数以及这两个标识符被称做控制模块。
cutoff：以纳米为单位的截矩。即超出截矩范围的非键连接相互作用将不计。
ntr：原子位置能量抑制位，1表示抑制，0表示不抑制。
imin：能量最优化标志位。1表示sander将进行能量最优化，0表示让sander进行分子动力学模拟。
macyc：能量最优化次数。
ncyc：便是经过多少次能量优化以后，能量优化从淬火过程变为梯度变化过程。
ntmin：能量优化方法标志位。0表示前10个能量最优化为淬火过程，然后进行梯度能量优化；1表示ncyc次淬火过程，然后进行梯度能量优化，为默认值；2表示只进行淬火过程。
dx0：表示启动模拟步长。
dxm：最大优化步数。
drms：梯度能量优化标准，默认值为1.0E-4& kcal/mol.A。
更多参数将在后文中解释。&&&
AMBER分子动力学简例（三）
分子动力学（2）&&
水环境中的分子动力学模拟
溶剂环境中的分子动力学模拟分为以下四步进行：
1、溶剂环境能量最优化。这一步保持溶质（蛋白）不变，去除溶剂中能量不正常的范德华相&&&
2、整系统能量最优化。去除整个系统中能量不正常的相互作用。
3、有限制的分子动力学。保持蛋白质不动，溶解溶剂的不同层，同时逐渐将系统温度从0K提升到300K。
4、整系统分子动力学模拟。在一个大气压，300K的环境下整个系统分子动力学模拟。可以得到成果的分子动力学模拟。
#############################
1、溶剂环境能量最优化。
该步骤的配置文件min1.in如下：
--------------------------------------------------------------------------------
oxytocin:& initial minimisation solvent +
maxcyc =& 1000,
ncyc = 500,
Hold the& protein fixed
--------------------------------------------------------------------------------
该过程保持肽链不动，其中500.0单位是kcal/mol，表示作用在肽链上使其不动的力。“RES&
1 9”表示肽链残基数目，因为我们学习使用的oxytocin有9个残基。
模拟命令如下：
sander &O &i min1.in& &o min1.out &p oxy.top &c
oxy.crd &r oxy_min1.rst &ref oxy.crd& &
2、整系统能量最优化。
配置文件min2.in如下：
--------------------------------------------------------------------------------
oxytocin:& initial minimisation whole system
maxcyc =& 2500,
ncyc = 1000,
--------------------------------------------------------------------------------
命令如下：
sander& &O &i min2.in &o min2.out &p oxy.top &c
oxy_min1.rst &r oxy_min2.rst& &
3、有限制的分子动力学。
第一步分子动力学保持蛋白分子位置不变，但是不是完全固定每个原子，同时缓解蛋白分子周围的水分子，是溶剂环境能量优化。在这个步骤中，我们将主要目的是对特定的原子使用作用力使其能量优化。我们要优化溶剂环境，至少需要10ps，我们将使用20ps用来优化我们上两步制作的分子系统的周期性边界的溶剂环境。
命令配置文件md1.in如下：
--------------------------------------------------------------------------------
oxytocin:& 20ps MD with res on protein
irest = 0,
tempi =& 0.0,
temp0 = 300.0,
gamma_ln = 1.0,
nstlim = 10000, dt =& 0.002,
ntpr = 100, ntwx = 500, ntwr = 1000
Keep protein fixed with& weak restraints
--------------------------------------------------------------------------------
上述参数解释如下：
ntb& = 1：表示分子动力学过程保持体积固定。
imin = 0：表示模拟过程为分子动力学，不是能量最优化。
nstlim =& #：＃表示计算的步数。
dt = 0.002：表示步长，单位为ps，0.002表示2fs。
temp0 =& 300：表示最后系统到达并保持的温度，单位为K。
tempi = 100：系统开始时的温度。&
gamma_ln =&
1：表示当ntt＝3时的碰撞频率，单位为ps-1（请参考AMBER手册）
ntt = 3：温度转变控制，3表示使用兰格氏动力学。
tautp =& 0.1：热浴时间常数，缺省为1.0。小的时间常数可以得到较好的耦联。
20.0：保持分子动力学稳定性速度极限。20.0为缺省值，当动力学模拟中原子速度大于极限值时，程序将其速度降低到极限值以下。
comp = 44.6：& 溶剂可压缩单位。
ntc = 2：Shake算法使用标志位。1表示不实用使用，2表示氢键将被计算，3表示所有键都将被计算在内。
tol =& #.#####：坐标位置重新设置的几何位置相对容忍度。
我们将使用一个较小的作用力，10kcal/mol。在分子动力学中，当ntr＝1时，作用力只需要5-10kcal/mol（我们需要引用一个坐标文件做分子动力学过程的比较，我们需要使用"-ref"参数）。太大的作用力同时使用Shake算法和2fs步长将使整个系统变得不稳定，因为大的作用力使系统中的原子产生大频率的振动，模拟过程并步需要。
运行命令如下：
sander& &O &i md1.in &o md1.out &p oxy.top &c
oxy_min2.rst &r oxy_md1.rst &x& oxy_md1.mdcrd &ref
oxy_min2.rst &inf& &
进行MD的运行时间一般较长，可以使用程序的并行版本提交集群计算。在主频位2.0GHz的P4单机上，大概需要一个钟头。可以随时查看md1.in文件的程序输出。
4、整系统分子动力学模拟。
这一步中，我们将进行整个系统的分子动力学模拟，而不对某些特定原子位置进行限制。因为知识一个小例子，我们将只进行250ps的MD计算。配置文件md2.in如下：
--------------------------------------------------------------------------------
oxytocin:& 250ps MD
imin = 0, irest = 1, ntx = 7,
ntb = 2, pres0 = 1.0,& ntp = 1,
taup = 2.0,
cut = 10, ntr = 0,
ntc = 2, ntf = 2,
tempi =& 300.0, temp0 = 300.0,
ntt = 3, gamma_ln = 1.0,
nstlim = 125000, dt =& 0.002,
ntpr = 100, ntwx = 500, ntwr =& 1000
--------------------------------------------------------------------------------
上面一些参数解释如下:
ntb=2：表示分子动力学过程的压力常数。
ntp＝1：表示系统动力学过程各向同性。
taup& = 2.0：压力缓解时间，单位为ps。
pres＝1：引用1个单位的压强。
使用以下命令进行MD：
sander &O &i md2.in& &o md2.out &p oxy.top &c
oxy_md1.rst & r oxy_md2.rst &x oxy_md2.mdcrd &ref&
oxy_md1.rst && &
模拟时间较长，在P4单机2.0GHz的上需要7.5个小时。
到此，模拟全部完成，接下来要对得到的数据进行分析。主要数据文件.out文件，包含系统能量、温度，压力等等；.mdcrd文件，是分子动力学轨迹文件，可以求系统蛋白的RMSD，回转半径等等。数据分析根据不同的研究目的不同而不同，我们将在后文中进行一些简单的分析。&&&
使用pdbviewer软件 load structure的时候，假设侧链缺失会 有提示错误信息。
直接使用xleap产生该结构的top&
crd后，再使用ambpdb重新产生structure即可加上侧链
已投稿到：
以上网友发言只代表其个人观点，不代表新浪网的观点或立场。