| 查看: 504 | 回复: 0 | ||
straysme铁虫 (初入文坛)
|
[求助]
求助：怎样用ambertools 添加保守水分子。。
|
|
在amber邮件求助中找到：Optionally, you can improve the positions of the hydrogen atoms with pol_h and gwh. Is so, you must put crystal water molecules in a separate file called watpdb, before running leap, and then run leap again with the final result file from gwh. Thus: 1. Move crystal waters to a separate file wat.pdb 2. Run tleap as above, but do not set pdbwritecharges, do not solvate and do not check: x=loadpdb pdb_file bond x.22.SG x.65.SG savepdb x pdbout saveamberparm x prmtop prmcrd quit 3. Run pol_h (it gives the position of the polar hydrogens, but not on water molecules): pol_h ppdbout This usually takes ~1 minute. This is the protonated file for the protein. Add a proper title to this file (a row starting with REMARK) 4. Convert the prmtop file to old format: new2oldparm prmtop.old 5. Run gwh (it calculates improved positions of the water hydrogens) gwh -w wat.pdb -elstat -p prmtop.old wpdb It usually takes about 5 minutes, and it may crash. Add this file to the protonated pdb file: cat wpdb >> ppdbout remove then END record before the water molecules 但是找不到pol_h and gwh 这两个程序。求各路大神帮助。。。谢谢@smutao@oxox6085@oxox6085@oxox6085 |
回复此楼» 猜你喜欢
2026国自然函评费到账
已经有24人回复
-
售SCI一区T0P文章，我:8.O.55.1.O.5.4,科目全,可+急
已经有5人回复
-
2026年8月25日国自然放榜前突然收到列入评审专家邮件，有关系吗？
已经有25人回复
-
售SCI一区文章，我:8.O.551.O.5.4,科目全,可伽急
已经有6人回复
-
售SCI一区T0P文章，我:8.O.55.1.O.5.4,科目全,可+急
已经有5人回复
-
售SCI一区T0P文章，我:8O.55.1.O.5.4,科目齐全,可+急
已经有4人回复
-
今天务委会开完了，明天出结果吗
已经有22人回复
-
在坚冰还盖着北海的时候，我看到了怒放的梅花。
已经有9人回复
-
明天应该可查了！？
已经有6人回复
-
有没有大神帮我看看基金代码
已经有10人回复
» 本主题相关商家推荐: (我也要在这里推广)
