| 查看: 1776 | 回复: 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 |
回复此楼» 收录本帖的淘帖专辑推荐
 分子模拟相关的 |
» 猜你喜欢
一志愿华中农微生物，288分，三年实验经历
已经有7人回复
-
269求调剂
已经有3人回复
-
268分085602化学工程调剂
已经有16人回复
-
071000生物学调剂求助
已经有4人回复
-
296求调剂
已经有5人回复
-
求调剂
已经有3人回复
-
297求调剂
已经有25人回复
-
一志愿中国科学院上海有机所，有机化学356分找调剂
已经有10人回复
-
367求调剂
已经有9人回复
-
086003调剂求助
已经有17人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 用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人回复
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
5楼2011-04-18 10:19:14
6楼2011-04-18 10:19:25
★ ★
zh1987hs(金币+2): 谢谢 2011-04-18 18:33:22
zh1987hs(金币+2): 谢谢 2011-04-18 18:33:22
|
在产生psf文件时使用的脚本是： 产生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 不过在产生psf时会出现下面的问题： 这个问题是不是正常的，如果不正常，怎样避免？？  |
» 本帖已获得的红花（最新10朵）
宋梦华7楼2011-04-18 10:20:00
8楼2011-04-18 10:30:50
9楼2017-03-27 09:36:41
