生物分子模拟论坛

 找回密码
 我想注册

QQ登录

只需一步,快速开始

扫一扫,访问微社区

查看: 7941|回复: 3

[其他] charmm分子动力学模拟尝试

[复制链接]
发表于 2014-11-30 15:41:12 | 显示全部楼层 |阅读模式

马上注册,结交更多好友,下载更多分子模拟资源。

您需要 登录 才可以下载或查看,没有帐号?我想注册

x
本帖最后由 lrf1980 于 2014-11-30 16:20 编辑

当把charmm安装好后,就来试试跑一下动力学模拟。基本的步骤参照了下面网站的步骤。
http://www.ch.embnet.org/MD_tutorial/
对于charmm来说由于程序本身用fortran写的,同时对pdb的要求比较多。所以对于初学者,有的时候甚至会一直困在如何得到一个能被charmm识别的pdb。前段时间,以为朋友教了我一个很好的办法,我觉得简单,而且不需要很多专业知识,现在介绍给大家。该方法就是用vmd软件来生成charmm能识别的pdb和psf文件。
现在安装好vmd后,直接用vmd打开我们建好的蛋白质模型。然后选extensions菜单,-> modeling -> automatic psf builder
这个时候程序弹出另外一个对话框。load topology file,这个基本上就是加载charmm力场。这个在charmm的安装包里面可以找到。
load完成后,就可以进行step2, 点击guess 菜单,进入step3, 点击create chains就能生成psf,进入step4, apply patches,pdb和psf文件就生成了。需要注意一下,如果你发现你找不到新生成的pdb和psf文件,在vmd的终端下看看目前的操作目录是什么,然后到那个目录下去找。
如果这些都顺利完成,恭喜你,你的pdb文件基本上charmm可以识别,你可以可以来开始运行charmm脚本了。

现在以上次我帮woniuzhixing做的不全残基的蛋白质为例来运行一下charmm的build, minimize,heat, equip几个脚本。

直接拷贝Home MD网站上的build脚本,然后将相关的信息修改。build脚本内容如下:


* Generate PSF and CHARMM coordinate set for
* protein 1qs4b4c
*

bomb -1

! Read in Topology and  Parameter files

open unit 1 card read name "./partop/top_all32_prot_lipid-all36_carb.rtf"
read RTF card unit 1
close unit 1


open unit 1 card read name "./partop/par_all32_prot_lipid-all36_carb.prm"
read PARA card unit 1
close unit 1

! Read sequence from the PDB coordinate file
open unit 1 card read name 1qs4b4c.pdb
read sequ pdb unit 1



! now generate the PSF and also the IC table (SETU keyword)

generate prot setu

! read in pdb coordinate file

open unit 1 card read name 1qs4b4c.pdb
read coor pdb unit 1
close unit 1


! build in hydrogens if using a crystal structure
hbuild  sele all end

! build in missing coordinates using values in
! the parameter set
ic fill
ic para
ic build

! print all coordinates:  undefined coordinates are indicated as
! 9999.00

print coor



! write out the protein structure file (psf) and
! the coordinate file in pdb format.  

open write formatted unit 27 name 1qs4b4c.psf
write psf card unit 27
* PSF for 1qs4
*


open unit 1 card write name 1qs4b4c-allh.pdb
write coor pdb unit 1
* 1qs4  Structure with hydrogens
*
close unit 1

stop


相比home MD网站,我在脚本开始设置了bomb -1. 另外我也把立场文件放在上层目录里的一个partop文件夹里。所以可以看到我在读立场时,用了“ ”指示立场文件的具体位置。
build.inp 完成后,运行 charmm <build.inp> build.out 就可以生成添加h原子的pdb和psf文件了。

用相同的方法可以运行minimize,heat, equip脚本得到相关的文件。至于每个脚本里面的具体内容,就需要各位自己去读charmm的相关文件了。

我把运行的脚本及结果贴在附件里,供大家参考。
有问题站内联系。


1QS4_Charmm.zip

870.47 KB, 下载次数: 44

charmm 脚本

评分

参与人数 1金币 +20 收起 理由
川大-灰太狼 + 20 很给力!

查看全部评分

发表于 2014-12-1 21:35:15 | 显示全部楼层
抛个charmm的license出来啊。
发表于 2015-10-25 19:56:58 | 显示全部楼层
对脚本的事情一直感到很神奇
发表于 2015-10-27 13:37:50 | 显示全部楼层
霸气                        
您需要登录后才可以回帖 登录 | 我想注册

本版积分规则

QQ|分迪科技|小黑屋|手机版|Archiver|生物分子模拟论坛 ( 蜀ICP备14009200号-3 )

GMT+8, 2024-3-29 09:52 , Processed in 0.071192 second(s), 23 queries .

Powered by Discuz! X3.4

Copyright © 2001-2020, Tencent Cloud.

快速回复 返回顶部 返回列表