24小时热门版块排行榜    

查看: 11790  |  回复: 210

bay__gulf

金虫 (著名写手)

刘苏州

★ ★
zeoliters(金币+2,VIP+0):多谢! 12-25 16:42
那一个帮助页上面有个Paratool's user guide 的链接
http://www.ks.uiuc.edu/Research/ ... ool/usersguide.html
81楼2009-09-18 21:19:13
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

bay__gulf

金虫 (著名写手)

刘苏州

★ ★
zeoliters(金币+2,VIP+0):多谢! 12-25 16:42
我整理了一份paratool 的中文指南,在http://www.mdbbs.org/thread-16875-1-1.html 的第21-23楼
82楼2009-09-20 10:24:50
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

arthurii

木虫 (正式写手)


小木虫(金币+0.5):给个红包,谢谢回帖交流
引用回帖:
Originally posted by bay__gulf at 2009-9-20 10:24:
我整理了一份paratool 的中文指南,在http://www.mdbbs.org/thread-16875-1-1.html 的第21-23楼

万分感谢!!
83楼2009-09-22 21:57:59
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

involute

新虫 (初入文坛)

请教PSF文件中关于IMPHI的写法


小木虫(金币+0.5):给个红包,谢谢回帖交流
您好!
     我最近用topology_all27_na.inp生成DNA的PSF文件的时候发现,对于如下的IMPROPER:

    A—B—C
         ▏
        D
     在psf文件的IMPHI行,是写成 B A C D的格式。但我查了一下NAMD关于bionano的一个tutorial,发现它给的sample DNA,对以上结构,采取的写法是 A C B D。
     另外,我查看了一下,这个tutorial里面,PDB和PSF中,atom name,atom type均采取Amber的命名格式,因为它最后计算采取的parameter file是cornell.prm,里面的atom type是AMBER格式的。

     所以,我现在碰到一个问题,直接采取top_all27_na.inp生成的IMPHI,在用cornell.prm参数运行NAMD时,它报错,因为不认识B A C D,所以找不到B A C D的参数。

     我想请教一下,是不是因为top_all27_na.inp采取的是CHARMM格式,而cornell.prm采取的是Amber格式,所以它们对IMPROPER的四个原子排列顺序不一样?

     如果是那样,我个人打算采取的解决方案是,直接修改top_all27_na.inp里面每个residue对IMPR的格式定义,使其跟cornell.prm一致,但我不知道这样会不会有别的意想不到的问题。

     ps:我的专业背景是微电子类,对生物知道的很少。请多多赐教!

    谢谢!
84楼2009-09-24 11:15:41
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

bay__gulf

金虫 (著名写手)

刘苏州

★ ★ ★ ★ ★ ★ ★ ★
gwdavid(金币+5,VIP+0):谢谢! 9-29 16:07
zeoliters(金币+3,VIP+0):多谢回帖交流! 12-25 16:43
amber 力场用于处理核酸,charmm力场用于处理蛋白质,这个算是个"潜规则"吧,但不是绝对要遵循的。

amber的形式同charmm 还有些不同,bionano例中应该是在charmm 的形式上借用了amber的参数。那个impr 的确有问题,需要手工修改。作者可能修改过了,在大多数的psf 文件开头都有 REMARKS topology cornell.rtf ,可惜这个神秘的cornell.rtf 没有在例中出现。

一套力场参数是经过反复拟合得到的,不同力场间的参数并不是完全通用。而charmm力场是namd 的原生力场,支持最好,且charmm 中同样有全套的核酸的力场参数,可以用这个。

amber 用在namd 中要经过比较复杂的设置。
http://www.mdbbs.org/thread-17292-1-1.html

[ Last edited by bay__gulf on 2009-9-24 at 12:39 ]
85楼2009-09-24 12:38:29
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

involute

新虫 (初入文坛)

回复:85


小木虫(金币+0.5):给个红包,谢谢回帖交流
非常感谢您的及时回答!

我仔细看过了,在topology_all27_na里面,对于A,T,C,G,只定义了两个Improper;而在bionano-tutorial里面,那个sample DNA的psf中,定义了5个,即凡是跟周围三个原子成键的原子,它都以此为中心,定义了一个improper。

这样看来,CHARMM和Amber对参数的定义,还有拟合方式是很不一样的。

所以,我决定还是修改topology_all27这个文件,直到最后得到的psf文件跟tutorial里面给出的sample DNA一致。

感谢你的指教!
86楼2009-09-24 14:00:31
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

bay__gulf

金虫 (著名写手)

刘苏州

★ ★
gwdavid(金币+2,VIP+0):谢谢! 9-29 16:07
你可以改改看,或是问问bionano的作者那个cornell.rtf 文件长得啥样子,
其实用charmm的核酸力场也不错。

还有,阁下是不是姓何

[ Last edited by bay__gulf on 2009-9-24 at 15:33 ]
87楼2009-09-24 14:51:24
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

arthurii

木虫 (正式写手)


小木虫(金币+0.5):给个红包,谢谢回帖交流
请教,NAMD跑完动力学,怎么算平均结构?还有用VMD算的SASA面积要求录入的那个参数是什么来的?看user guide不是很明白,怎么可以分别算极性和非极性的部分?
谢谢
88楼2009-09-26 09:12:32
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

bay__gulf

金虫 (著名写手)

刘苏州

★ ★ ★
gwdavid(金币+3,VIP+0):谢谢! 9-29 16:07
measure avpos 和measure sasa可以分别完成你的要求
参考ug107和109页,很清楚的。
我在这里再怎么表述也不会比ug更清楚。
atomselect 中有一个keyword:polar,不知道能否判断极性。
89楼2009-09-26 10:16:58
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

arthurii

木虫 (正式写手)


小木虫(金币+0.5):给个红包,谢谢回帖交流
引用回帖:
Originally posted by bay__gulf at 2009-9-26 10:16:
measure avpos 和measure sasa可以分别完成你的要求
参考ug107和109页,很清楚的。
我在这里再怎么表述也不会比ug更清楚。
atomselect 中有一个keyword:polar,不知道能否判断极性。

Thank you very much. 版主动作真快,辛苦了
90楼2009-09-26 14:55:47
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 bay__gulf 的主题更新
普通表情 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 材料学硕318求调剂 +10 February_Feb 2026-03-01 12/600 2026-03-01 23:34 by L135790
[考研] 265分求调剂不调专业和学校有行学上就 +6 礼堂丁真258 2026-02-28 8/400 2026-03-01 22:50 by jian_
[考研] 272求调剂 +6 田智友 2026-02-28 6/300 2026-03-01 21:40 by 公瑾逍遥
[考研] 0856求调剂285 +10 吕仔龙 2026-02-28 10/500 2026-03-01 21:37 by 公瑾逍遥
[考研] 299求调剂 +3 Y墨明棋妙Y 2026-02-28 5/250 2026-03-01 21:01 by tangxiaotian
[考研] 291分工科求调剂 +9 science饿饿 2026-03-01 10/500 2026-03-01 18:55 by 18137688336
[考研] 328求调剂 +3 aaadim 2026-03-01 5/250 2026-03-01 17:29 by njzyff
[考研] 281求调剂 +4 2026计算机_诚心 2026-03-01 7/350 2026-03-01 17:20 by 2026计算机_诚心
[考研] 321求调剂一志愿东北林业大学材料与化工英二数二 +4 虫虫虫虫虫7 2026-03-01 7/350 2026-03-01 16:52 by caszguilin
[考研] 化工专硕342,一志愿大连理工大学,求调剂 +3 kyf化工 2026-02-28 4/200 2026-03-01 16:49 by yywzz
[考研] 304求调剂 +6 曼殊2266 2026-02-28 7/350 2026-03-01 15:14 by wjLi2017
[考研] 303求调剂 +4 今夏不夏 2026-03-01 4/200 2026-03-01 14:46 by 嘟嘟小浣熊
[考研] 课题组接收材料类调剂研究生 +3 gaoxiaoniuma 2026-02-28 4/200 2026-03-01 14:30 by jjj三跨
[考研] 295复试调剂 +3 简木ChuFront 2026-03-01 3/150 2026-03-01 14:27 by zzxw520th
[考研] 284求调剂 +6 天下熯 2026-02-28 6/300 2026-03-01 14:19 by Ducount.Y
[考研] 317一志愿华南理工电气工程求调剂 +6 Soliloquy_Q 2026-02-28 11/550 2026-03-01 11:14 by 歌liekkas
[论文投稿] 求助coordination chemistry reviews 的写作模板 10+3 ljplijiapeng 2026-02-27 4/200 2026-03-01 09:07 by babero
[论文投稿] Optics letters投稿被拒求助 30+3 luckyry 2026-02-26 4/200 2026-03-01 09:06 by babero
[考研] 材料调剂 +4 爱擦汗的可乐冰 2026-02-28 4/200 2026-03-01 00:38 by 猫猫球alter
[硕博家园] 【博士招生】太原理工大学2026化工博士 +4 N1ce_try 2026-02-24 8/400 2026-02-26 08:40 by N1ce_try
信息提示
请填处理意见