注册 登录  
 加关注
查看详情
   显示下一条  |  关闭
温馨提示!由于新浪微博认证机制调整,您的新浪微博帐号绑定已过期,请重新绑定!立即重新绑定新浪微博》  |  关闭

Wei Ding blog

 
 
 

日志

 
 

gromacs的基本使用  

2011-06-21 11:47:37|  分类: 默认分类 |  标签: |举报 |字号 订阅

  下载LOFTER 我的照片书  |
一:利用pdb文件产生gromacs能识别的文件。
pdb2gmx -f 2NW7.pdb  -o conf.gro -p topol.top -i posre.itp -n clean.ndx -q clean.pdb -water tip3p
二:决定进行模拟的空间大小(建议不使用该命令)
editconf -f conf.gro  -o editconf.gro -box 10nm
三:产生可以用pymol查看的文件
editconf -f editconf.gro -o pymol.pdb
四:进行能量最小化
先编写em.mdp
integrator= steep
nsteps= 200
nstlist= 10
rlist= 1.0
coulombtype= pme
rcoulomb = 1.0
vdw-type = cut-off
rvdw = 1.0
nstenergy = 10
然后
grompp -f em.mdp -po mdout.mdp -c conf.gro -p topol.top  -o topol.tpr
mdrun -s topol.tpr  -o traj.trr -x traj.xtc -c confout.gro -e ener.edr -g md.log
得到结果
Potential Energy  = -1.8466966e+05
Maximum force     =  8.7827275e+03 on atom 7711
Norm of force     =  1.4047519e+02
五:判断结构的合理性(RMSD)
 g_rms -s topol.tpr  -f traj.trr -o rmsd.xvg
xmgrace rmsd.xvg
六:比较结构扰动与结构因子的关系。
g_rmsf -f traj.trr -s topol.tpr -o rmsf.xvg
xmgrace rmsf.xvg
七:查看二级结构
do_dssp  -f traj.trr -s topol.tpr -o ss.xpm -sc scount.xvg
xpm2ps -f       root.xpm-o       plot.eps
八:Distance and hydrogen bonds
make_ndx -f conf.gro -o index.ndx
  评论这张
 
阅读(555)| 评论(0)
推荐 转载

历史上的今天

在LOFTER的更多文章

评论

<#--最新日志,群博日志--> <#--推荐日志--> <#--引用记录--> <#--博主推荐--> <#--随机阅读--> <#--首页推荐--> <#--历史上的今天--> <#--被推荐日志--> <#--上一篇,下一篇--> <#-- 热度 --> <#-- 网易新闻广告 --> <#--右边模块结构--> <#--评论模块结构--> <#--引用模块结构--> <#--博主发起的投票-->
 
 
 
 
 
 
 
 
 
 
 
 
 
 

页脚

网易公司版权所有 ©1997-2018