| 查看: 2206 | 回复: 8 | ||||
| 当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖 | ||||
[交流]
【求助】产生psf文件时出现的Missing atoms for conformation definition 已有2人参与
|
||||
在产生psf文件时，结束后。shell出现了这样的信息： 这是什么原因？？是原始pdb文件的原因吗？？ 这是我产生psf的脚本： 进行psf生成前的文件准备(prepfiles.tcl)： mol load pdb 1J4N-ALL.pdb foreach S {A B C D} { set seg [atomselect top "segname $S and chain A"] $seg writepdb seg$S.pdb $seg delete } foreach S {A B C D} { set wat [atomselect top "segname $S and resname HOH and within 100 of chain A"] $wat writepdb crystwat$S.pdb $wat delete } exit 产生psf文件（build.psf）: package require psfgen topology ../top_all27_prot_lipid.rtf pdbalias residues HIS HSE pdbalias atom ILE CD1 CD pdbalias atom HOH O OH2 pdbalias resid HOH TIP3 foreach S {A B C D} { segment $S { pdb seg$S.pdb } coordpdb seg$S.pdb $S regenerate angles dihedrals segment WC$S { auto none pdb crystwat$S.pdb } coordpdb crystwat$S.pdb WC$S } guesscoord writepdb aqp0.pdb writepdb aqp0.psf exit |
回复此楼» 收录本帖的淘帖专辑推荐
 分子模拟相关的 |
» 猜你喜欢
哪位老哥知道今年的国自然具体哪一天放榜？
已经有8人回复
-
咱们一起用铁证分析2026国家社科基金中标与否
已经有25人回复
-
我的国基提前知道中了，可是同事的操作让我实在接受不了，怎么会有这样的人
已经有16人回复
-
欢迎发来filecode的Mz6后的代码验证其规律
已经有70人回复
-
有时候，自然基金真的不能太认真 （我的申报经验）
已经有12人回复
-
招聘 金属平磨液，抛光液研发工程师
已经有3人回复
-
时间戳又变了8-15
已经有23人回复
-
filecode=后面第一个是大写字母
已经有7人回复
-
filecode
已经有12人回复
-
小木虫上这么多卖论文的，真有人买论文么？感觉没必要啊
已经有13人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 用VASP计算Mg的PDOS的时候出现d电子是什么原因？
已经有6人回复
- 求立方结构（cubic）的La1-xSrxCoO的prc或cif文件
已经有4人回复
- FORTRAN新手 求助主程序循环问题
已经有10人回复
- 利用概率密度函数生成随机变量过程中出现超几何分布函数问题
已经有4人回复
- 植物体产生的气体怎么收集?
已经有10人回复
- 【the song I want to sing for u】My love by 小猫（重在参与）O(∩_∩)O~
已经有94人回复
- 请教Atoms too close错误
已经有8人回复
- 请问在做药物分子对接时药物分子结构文件的准备
已经有9人回复
- 【求助】fortran计算时出现NaN
已经有11人回复
6楼2011-04-18 10:19:25
2楼2011-04-16 16:32:00
★ ★ ★ ★ ★
御剑江湖(金币+5): 谢谢 2011-04-18 10:12:08
御剑江湖(金币+5): 谢谢 2011-04-18 10:12:08
|
找到原因了，也就是我在运行prepfiles.tcl这一步后，产生的segA.pdb segB.pdb segC.pdb segD.pdb，在这些文件的通道里还是含有水分子（这一步按照教程的理解,应该是蛋白和水分子要分别在不同的pdb文件里），而产生的crystwatA.pdb crystwatB.pdb crystwatC.pdb crystwatD.pdb是正常的，也就是只含有通道里的水分子 。 这是到prefiles.tcl这一步的操作流程，大家看看我哪一步需要修改： 1. 通过vmd先找到通道中的水分子，获取它们的编号 2. 然后生成这些水分子的pdb文件（例子）： mol load pdb 1J4N.pdb set model [atomselect top "water and (resid 301 or resid 302 or resid 303 or resid 304)"] $model writepdb wat.pdb exit 3. 生成蛋白的pdb文件： mol load pdb 1J4N.pdb set model [atomselect top protein] $model writepdb 1j4n.pdb exit 4. 合并wat.pdb 与1j4n.pdb 这一步不用写脚本，直接把wat.pdb的内容（即去头去尾）拷贝到1j4n.pdb的END前即可。 5. over。 构建四聚体(buildtetra.tcl): mol load pdb 1j4n.pdb set all [atomselect top all] $all set segname A $all writepdb 1J4N-A.pdb $all delete set sel [atomselect top all] $sel set segname B $sel move {{-1.0 0.0 0.0 93.331} {0.0 -1.0 0.0 93.331} {0.0 0.0 1.0 0.0} {0.0 0.0 0.0 1.0}} $sel writepdb 1J4N-B.pdb $sel delete mol delete top mol load pdb 1j4n.pdb set sel [atomselect top all] $sel set segname C $sel move {{0.0 -1.0 0.0 93.331} {1.0 0.0 0.0 0.0} {0.0 0.0 1.0 0.0} {0.0 0.0 0.0 1.0}} $sel writepdb 1J4N-C.pdb $sel delete mol delete top mol load pdb 1j4n.pdb set sel [atomselect top all] $sel set segname D $sel move {{0.0 1.0 0.0 0.0} {-1.0 0.0 0.0 93.331} {0.0 0.0 1.0 0.0} {0.0 0.0 0.0 1.0}} $sel writepdb 1J4N-D.pdb $sel delete cat 1J4N-A.pdb 1J4N-B.pdb 1J4N-C.pdb 1J4N-D.pdb > 1J4N-ALL.pdb 然后打开1J4N-ALL.pdb,搜索END,把中间的三个END及相应的杂项去掉，保存为1J4N-ALL.pdb。此时就生成了四聚体pdb文件。 进行psf生成前的文件准备(prepfiles.tcl)： mol load pdb 1J4N-ALL.pdb foreach S {A B C D} { set seg [atomselect top "segname $S and chain A"] $seg writepdb seg$S.pdb $seg delete } foreach S {A B C D} { set wat [atomselect top "segname $S and resname HOH and within 100 of chain A"] $wat writepdb crystwat$S.pdb $wat delete } exit |
3楼2011-04-18 09:58:21
4楼2011-04-18 10:01:51
