小木虫 --- 500万硕博科研人员喜爱的学术科研平台
&&查看话题
用NAMD将pdb生成psf总是不成功，请问为什么呢？
我想将一个酶的PDB转换为psf，用了两种方法，一种是制作pgn指令运行，一种是使用Extensions→Modeling→Automatic PSF Builder菜单项。
*************************************
其中，制作pgn指令运行会出现以下错误：
Main console display active (Tcl8.5.6 / Tk8.5.6)
(VMD) 1 % source psf820.pgn
invalid command name &\rtf1\ansi\ansicpg936\deff0\nouicompat\deflang1033\deflangfe2052{\fonttbl{\f0\fswiss\fcharset134 SimS}{\f1\fswiss\fcharset0 SimS}{\f2\fnil\fcharset134 \'cb\'ce\'cc\'e5;}}
{\\red0\green0\blue0;}
{\*\generator Riched20 6.3.9600}\viewkind4\uc1
\pard\sb140\cf1\f0\fs20\lang9 package require psfgen \par
topology top all27_prot_lipid.inp \par
pdbalias residue HIS HSE \par
pdbalias atom ILE CD1 CD \par
segment U \{pdb \f1\lang1033 papainpdb\f0\lang9 .pdb\} \par
coordpdb \f1\lang1033 papainpdb\f0\lang9 .pdb U \par
guesscoord \par
writepdb \f1\lang1033 papainpdbpsf\f0\lang9 .pdb \par
writepsf \f1\lang1033 papainpdbpsf\f0\lang9 .psf \cf0\f2\fs22\par
&Main& (VMD) 2 %
---------------------------------
我的配置文件pgn是这样的：
package require psfgen
topology top all27_prot_lipid.inp
pdbalias residue HIS HSE
pdbalias atom ILE CD1 CD
segment U {pdb papainpdb.pdb}
coordpdb papainpdb.pdb U
guesscoord
writepdb papainpdbpsf.pdb
writepsf papainpdbpsf.psf
-------------------------------
***********
然后我又试了Automatic PSF Builder
具体的详情是：
Unable to load structure file papainpdb_autopsf_tmpfile.psf
Unable to load structure file papainpdb_autopsf_tmpfile.psf
& & while executing
&mol load psf ${basename}_tmpfile.psf pdb ${basename}_tmpfile.pdb&
& & (procedure &checkparams& line 4)
& & invoked from within
&checkparams&
& & (procedure &::autopsf::makepatches_gui& line 62)
& & invoked from within
&::autopsf::makepatches_gui&
& & invoked from within
&.autopsf.patches.finish invoke&
& & (&uplevel& body line 1)
& & invoked from within
&uplevel #0 &
& & (procedure &tk::ButtonUp& line 24)
& & invoked from within
&tk::ButtonUp .autopsf.patches.finish&
(command bound to event)
请问我的问题到底出在哪里呢？谢谢！！
请问是加在哪里？
研究生必备与500万研究生在线互动！
扫描下载送金币
浏览器进程
打开微信扫一扫
随时随地聊科研君，已阅读到文档的结尾了呢~~
介绍了如何利用NAMD软件在粗粒化水平上做分子模拟。
扫扫二维码，随身浏览文档
手机或平板扫扫即可继续访问
利用NAMD计算粗粒化模型
举报该文档为侵权文档。
举报该文档含有违规或不良信息。
反馈该文档无法正常浏览。
举报该文档为重复文档。
推荐理由：
将文档分享至：
分享完整地址
文档地址：
粘贴到BBS或博客
flash地址：
支持嵌入FLASH地址的网站使用
html代码：
&embed src='/DocinViewer-4.swf' width='100%' height='600' type=application/x-shockwave-flash ALLOWFULLSCREEN='true' ALLOWSCRIPTACCESS='always'&&/embed&
450px*300px480px*400px650px*490px
支持嵌入HTML代码的网站使用
您的内容已经提交成功
您所提交的内容需要审核后才能发布，请您等待！
3秒自动关闭窗口NAMD入门教程_图文-五星文库
免费文档下载
NAMD入门教程_图文
导读：solvateubq.psfubq.pdbCt5Coubq_wb等待运行结束，VMD就调用solvate将ubq.pdb和ubq.psf所储存的蛋白放入一个立方水体中。在图形窗口可以见到一个立方形的水体包围蛋白（如图）。参数Ct5通知程序如何确定立方体的各边长。方法是在每个坐标方向上选择坐标最大的那个原子，然后再延伸5A，即为该方向立方体面的边界。注意：生成
solvate ubq.psf ubq.pdb Ct 5 Co ubq_wb
等待运行结束，VMD就调用solvate将ubq.pdb和ubq.psf所储存的蛋白放入一个立方水体中。在图形窗口可以见到一个立方形的水体包围蛋白（如图）
参数 Ct 5 通知程序如何确定立方体的各边长。方法是在每个坐标方向上选择坐标最大的那个原子，然后再延伸5A，即为该方向立方体面的边界。注意：生成的立方水体并不一定是正立方体。各边长取决于坐标最大（距离原点最远）的原子的位置。还有一个参数 Co ubq_wb是为了通知程序生成的文件名。运行结束后我们得到的两个文件就是ubq_wb.psf 和 ubq_wb.pdb。
2、在VMD TkConsole中输入：
set everyone [atomselect top all]
measure minmax $everyone
这时返回的数值是整个体系中离原点最近的点和最远的点的坐标。我们需要的是整个立方体的中心，可以自己计算也可以用下面的命令：
measure center $everyone
这时返回的三个数值就是体系的中心。记下这三个数值，我们以后还会用到。
返回值如图：
在开始下一节之前，我们要将生成的文件拷贝到common公用目录下以方便访问。在Windows资源管理器中找到1-1-build目录，按Ctrl选择以下六个文件：ubq.pdb, ubq.psf, ubq_ws.pdb, ubq_ws.psf, ubq_wb.pdb, ubq_ws.psf，然后把它们拷贝到namd-tutorial/common目录下。 提示：
此处生成的立方水体事实上过小了。在实际应用时，应当保证水体足够大，以防止蛋白在拉伸运动时超出水环境。也要避免在使用周期边界条件时蛋白和四周各个单元的蛋白镜像相碰撞。周期边界条件的详细知识见下文。
此外，还应当注意的是在实际应用时应当在水环境中放入离子。特别是当蛋白质的净电荷不为0时，更应当设定离子数目以使得整个体系是中性的。在放入离子时应当将它们放在体系中静电势能的最低点，以节省计算时间。因为离子总会向势能最低点自发运动。
球状水体中泛素（Ubiquitin）的分子动力学模拟
在这一节中，我们将要对球状水体的泛素分子进行最简单的动力学模拟。
首先，我们要进行的分子动力学模拟的目的是什么？我们将泛素放入球状水体中，水体周围是真空，然后NAMD会根据我们设定好的温度值按照Boltzmann-Maxwell分子速
率分布给各原子赋予一定的初始速度，接下来就是要根据牛顿力学方程，求解个水分子以及蛋白质中各原子的运动轨迹。我们得到的结果，就是模拟泛素这一小肽在溶液状态下的运动状态。 知识链接：能量最小化和能量平衡（Minimization and Equilibration）
本次动力学模拟实际包括两个过程：能量最小化和能量平衡（Minimization and Equilibration）。能量最小化时，NAMD设定各原子的速度为0，然后不断改变各个原子的相对位置并计算体系总能量，搜索最低势能点，作为分子动力学模拟的初始状态。这一过程是不记录原子运动轨迹的。因为原子的位置改变只是因为NAMD需要搜索最低能量状态，而不是真实的相互作用引起的运动。
能量平衡是让蛋白质和水分子在设定好的环境温度（即原子的速度）下相互作用，达到能量平衡分配，整个体系达到稳定状态（熵达到极大值）。
为什么需要首先进行能量最小化？这是因为我们提供的体系有可能包含极度扭曲，拉伸或压缩变形的键和键角。它们是解析结构或同源模建时引入的错误结构，含有很高的能量。如果不首先进行最小化，直接进行能量平衡，蛋白质会和水分子相互作用，恢复伸展状态，释放掉这些错误结构中的高能量。这一过程是没有意义的――因为它是错误结构引起的反应，并不是蛋白质在溶液中的真实状态。从而就浪费了计算时间。不仅如此，能量释放引起的剧烈运动和相互作用最终可能使得蛋白质的行为不符合溶液中的真实行为。所以有必要在能量平衡之前，首先人为搜索能量最低点，作为分子动力学模拟的初始状态。
一般地，分子动力学模拟包括多个能量最小化和平衡过程。通常我们会首先将蛋白质固定而仅允许水分子运动，进行能量最小化和能量平衡；然后允许蛋白质和水分子同时运动，再次经历能量最小化和能量平衡这一循环。第一步的目的是使水分子达到能量最小，这通常是一个很快的过程。然后再放开蛋白质，使整个系统达到能量最小。这样可以减小计算量，并防止由于一开始蛋白结构很不稳定而结果产生假象。
在上一节中我们已经获得了所需要的ubq_ws.pdb 和ubq_ws.psf 两个文件。对照本章开始提到NAMD所需的四个文件知，还需要有配置文件就可以提交NAMD进行动力学模拟了。（力场参数文件在common文件夹中）。下面我们将首先得到配置文件。
2.3.1 配置文件
前面我们把使用NAMD比作写信，这里的“信”就是指的配置文件。配置文件记录了进行动力学模拟所需的全部参数和设置，NAMD只要得到这一文件就可以按照相应指令进行操作。对于本次动力学模拟，在namd-tutorial/1-2-sphere 目录下可以得到已经预先制作完毕的配置文件。下面我们将要仔细讲述文件的内容。
打开写字板，在菜单中选择 文件→打开，找到1-2-sphere 目录，文件类型选择全
部文件（图），然后打开文件ubq_ws_eq.conf。这个文件看起来好像很复杂，但是我们会
仔细分析讲解每一部分的含义。
注意：在配置文件中，每一行开头如果是“#”，则本行内容会被当作注释对待，NAMD会忽略其中的内容。因此为了便于区分，我们用#把文件分割成几大部分。如第一部分是：
############################################
## JOB DESCRIPTION
############################################
意思是这一部分是对所提交工作的描述。
1、大体浏览一下，可以发现整个文件被分成了以下几部分：
? 工作描述（job description）
? 可调参数 （adjustable parameters）
? 动力学模拟参数（simulation parameters）
? 附加参数（extra parameters）
? 执行脚本 （execution script）
它完整的记录了输入的蛋白质结构文件（pdb和psf文件）的位置，输出结果文件的文件名，以及动力学模拟时的环境温度，截止点，步长等各种参数。在进行动力学模拟时只需要提供给NAMD一个配置文件，NAMD就可以找到输入文件，调整好各种参数，按照要求进行动力学模拟之后输出结果。
2、然后，首先我们来看第一部分：Job Description。这一部分每一行开头都有#，因此只包括注释。它描述的是这一配置文件的目的：Minimization and Equilibration of Ubiquitin in a Water Sphere。就相当于一片文章的题目。
3、Adjustable Parameters 这一部分包括5项参数：
? structure: 给出调用的psf文件的位置
? coordinates：给出调用的坐标文件（即pdb文件）的位置
? set temperature 310: 定义一个变量temperature，并赋值310。以后如果要使
包含总结汇报、旅游景点、办公文档、文档下载、外语学习、资格考试、IT计算机以及NAMD入门教程_图文等内容。本文共25页
相关内容搜索小木虫 --- 500万硕博科研人员喜爱的学术科研平台
&&查看话题
既然NAMD是模拟多肽的结构，为什么有需要一个pdb文件作为input呢？
如题，我是今天才学namd的新手，想模拟一段10个氨基酸的序列的loop，但是在namd教程上看到namd第一步就是需要一个pdb文件，我既然不知道结构就当然不会有pdb文件了，那这是怎么回事呀，求高手解答:rol::rol::rol:
那我能不能在pymol里自己建一个pdb文件来进行模拟呢？求大神详细解释一下
当然可以&&可以画PDB文件的软件有很多
CHEMOFFICE也可以& &MS也可以
如果你要模拟的分子常见的话，也可以在蛋白质数据库里面下载
研究生必备与500万研究生在线互动！
扫描下载送金币
浏览器进程
打开微信扫一扫
随时随地聊科研}