24小时热门版块排行榜    

查看: 169  |  回复: 1
当前主题已经存档。

panhuairen

金虫 (小有名气)

[交流] 【求助】Amber 的Xleap 遇到的问题

为了把InsightII 产生的RNA的pdb文件进行分子模拟,就将该文件到xleap下进行跑动力学之前的准备,以便生成计算前的文件格式:
Here are commands:
________________________________________________________________________
>[root@localhost amber_tutorial]# $AMBERHOME/exe/tleap -s -f $AMBERHOME/dat/leap/cmd/leaprc.rna.ff99
-I: Adding /usr/local/amber9/dat/leap/prep to search path.
-I: Adding /usr/local/amber9/dat/leap/lib to search path.
-I: Adding /usr/local/amber9/dat/leap/parm to search path.
-I: Adding /usr/local/amber9/dat/leap/cmd to search path.
-s: Ignoring startup file: leaprc
-f: Source /usr/local/amber9/dat/leap/cmd/leaprc.rna.ff99.

Welcome to LEaP!
Sourcing: /usr/local/amber9/dat/leap/cmd/leaprc.rna.ff99
Log file: ./leap.log
Loading parameters: /usr/local/amber9/dat/leap/parm/parm99.dat
Reading title:
PARM99 for DNA,RNA,AA, organic molecules, TIP3P wat. Polariz.& LP incl.02/04/99
Loading library: /usr/local/amber9/dat/leap/lib/all_nucleic94.lib
Loading library: /usr/local/amber9/dat/leap/lib/all_amino94.lib
Loading library: /usr/local/amber9/dat/leap/lib/all_aminoct94.lib
Loading library: /usr/local/amber9/dat/leap/lib/all_aminont94.lib
Loading library: /usr/local/amber9/dat/leap/lib/ions94.lib
Loading library: /usr/local/amber9/dat/leap/lib/solvents.lib
> x = loadpdb "final_bf_dm.pdb"
Loading PDB file: ./final_bf_dm.pdb
Created a new atom named: P within residue: .R
Created a new atom named: O1P within residue: .R
Created a new atom named: O2P within residue: .R
Created a new atom named: H2* within residue: .R
Created a new atom named: HO2* within residue: .R
................................
.............
.......
Created a new atom named: HO2* within residue: .R
Created a new atom named: H2* within residue: .R
Created a new atom named: HO2* within residue: .R
Created a new atom named: H2* within residue: .R
Created a new atom named: HO2* within residue: .R
Created a new atom named: H2* within residue: .R
Created a new atom named: HO2* within residue: .R
Created a new atom named: H2* within residue: .R
Created a new atom named: HO2* within residue: .R
  total atoms in file: 1944
  Leap added 194 missing atoms according to residue templates:
       70 Heavy
       124 H / lone pairs
  The file contained 195 atoms not in residue templates
>
___________________________________________________________________
谁知道问题出在哪里了呢?

[ Last edited by zzgyb on 2008-7-20 at 18:59 ]
回复此楼
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

xxffliu

铜虫 (小有名气)

有些残基类型在参数文件中没有定义,可以试试leaprc.ff99或leaprc.ff03。amber主页上的第一个tutorial介绍了你出现的情况,要手动修改这些残基,可以看看。
2楼2008-07-12 22:44:15
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 panhuairen 的主题更新
普通表情 高级回复 (可上传附件)
信息提示
请填处理意见