计算化学公社

 找回密码 Forget password
 注册 Register
Views: 11480|回复 Reply: 6
打印 Print 上一主题 Last thread 下一主题 Next thread

[GROMACS] 配體拓樸生成問題 (已解決)

[复制链接 Copy URL]

63

帖子

0

威望

452

eV
积分
515

Level 4 (黑子)

本帖最后由 HarrisLin 于 2017-5-10 17:26 编辑


參考6F Sob老師的教學

個人做法如下:
先啟動amber
將acpype.py和配體.mol2扔同一資料夾
終端機輸入python acpype.py -i lead4l.mol2
開始運算約五分鐘完成( i7 6700)
============
目前是要做蛋白-配體的Gromacs MD
參考這篇的教學 點我
配體透過網站運算後得到拓樸.itp .gro再送回蛋白拓樸中

目前使用了PROFRG和TPPMKTOP都碰上了點問題所以上來請教

1.PROFRG產生的拓樸需要手動修正
不過原文教學的說明比較攏統  
不太明白如何修正(個人是生物背景)
想問問看是否有比較詳細的教學資料

2.參考 點我 這篇使用TPPMKTOP產生結果
產生三個檔案 file.itp  file.rtp lack.itp
我看了下應該file.itp file.rtp(對應.gro)都能直接使用在Gromacs
那lack.itp這個要補充在哪裡還是直接忽略呢?

3.或是使用AmberTool是最佳方式?
因為我雖然安裝成功仍有些使用技巧尚未克服
如果要改用可能得上來再多請教一下

謝謝

336

帖子

2

威望

4531

eV
积分
4907

Level 6 (一方通行)

2#
发表于 Post on 2017-5-9 09:00:07 | 只看该作者 Only view this author
这个问题暗含力场选择、小分子电荷以及力场参数获得的问题。
建议使用antichamber
http://ambermd.org/antechamber/antechamber.html

5万

帖子

99

威望

5万

eV
积分
112496

管理员

公社社长

3#
发表于 Post on 2017-5-9 12:43:34 | 只看该作者 Only view this author
那个翻译来的教程用的做法不怎么样。

建议参看
几种生成有机分子GROMACS拓扑文件的工具
http://sobereva.com/266

PRODRG早已经过时,J. Chem. Inf. Model. 2010, 50, 2221–2235里早已经曝光PRODRG对charge group有明显问题。现在要产生GROMOS力场的都用ATB了
靠acpype调用ambertools产生GAFF力场的拓扑文件用起来极为方便,安装可参考
Amber14安装方法
http://sobereva.com/263

另外在我上面给的博文里还介绍了怎么产生CGenFF力场的拓扑文件,可以结合GROMOS里自带的CHARMM力场使用。

我不推荐那几个产生OPLS力场的工具。

北京科音自然科学研究中心http://www.keinsci.com)致力于计算化学的发展和传播,长期开办高质量的各种计算化学类培训:初级量子化学培训班中级量子化学培训班高级量子化学培训班量子化学波函数分析与Multiwfn程序培训班分子动力学与GROMACS培训班CP2K第一性原理计算培训班,内容介绍以及往届资料购买请点击相应链接查看。这些培训是计算化学从零快速入门以及进一步全面系统性提升研究水平的高速路!培训各种常见问题见《北京科音办的培训班FAQ》
欢迎加入“北京科音”微信公众号获取北京科音培训的最新消息、避免错过网上有价值的计算化学文章!
欢迎加入人气非常高、专业性特别强的综合性理论与计算化学交流QQ群“思想家公社QQ群”:1号:18616395,2号:466017436,3号:764390338,搜索群号能搜到哪个说明目前哪个能加,合计9000人。北京科音培训班的学员在群中可申请VIP头衔,提问将得到群主Sobereva的最优先解答。
思想家公社的门口Blog:http://sobereva.com(发布大量原创计算化学相关博文)
Multiwfn主页:http://sobereva.com/multiwfn(十分强大的量子化学波函数分析程序)
ResearchGate:https://www.researchgate.net/profile/Tian_Lu
Money and papers are rubbish, get a real life!

63

帖子

0

威望

452

eV
积分
515

Level 4 (黑子)

4#
 楼主 Author| 发表于 Post on 2017-5-9 14:59:10 | 只看该作者 Only view this author
本帖最后由 HarrisLin 于 2017-5-9 15:03 编辑
sobereva 发表于 2017-5-9 12:43
那个翻译来的教程用的做法不怎么样。

建议参看

感謝Sob老師回答

我原始檔案是PDB
參考這篇 點我 用Amber製作配體拓樸
目前已經照個步驟產生了 .inpcrd和.prmtop( 這篇 說需要這兩個檔案跑acpype
不過後面的操作我就有點困惑
Sob老師給的 acpype網站 我看完Readme仍不太明白如何下載 系統有有內建python 2.7
後來是抓這個版本 點我
以及安裝後是像gmx指令直接啟用還是像amber要先進資料夾source amber.sh啟用呢?

另外我這方案和老師描述的步驟似乎不太一樣?
推荐用acpype,用法很简单,要处理xxx.mol2就执行./acpype.py -i xxx.mol2,然后程序会自动调用Antechamber处理xxx.mol2,算完后acpype就会将之转化为带有_GMX后缀的.gro、.itp、.top,直接在Gromacs里用即可
如果老師的方案才是比較簡便的
想請問哪裡找的到類似的詳細教學呢(中英皆可)
因為對ubuntu系統還在摸索階段
比較需要詳細的教學

謝謝

545

帖子

0

威望

3123

eV
积分
3668

Level 5 (御坂)

5#
发表于 Post on 2017-5-9 19:27:53 | 只看该作者 Only view this author
HarrisLin 发表于 2017-5-9 14:59
感謝Sob老師回答

我原始檔案是PDB

话说你准备用什么力场?我记得PRODRG只支持GROMOS力场,Antechamber只支持AMBER力场,这两种方案没什么横向比较的价值啊…

不过acpype确实挺好用的。实在不行你就把acpype.py源码下载下来,然后python acpype.py -i xxx.mol2算了。

(其实py3也能用来跑acpype,源码里改几个小地方就行…不过我估计你用不上就是了)

评分 Rate

参与人数
Participants 1
eV +1 收起 理由
Reason
HarrisLin + 1 谢谢

查看全部评分 View all ratings

5万

帖子

99

威望

5万

eV
积分
112496

管理员

公社社长

6#
发表于 Post on 2017-5-9 19:44:27 | 只看该作者 Only view this author
HarrisLin 发表于 2017-5-9 14:59
感謝Sob老師回答

我原始檔案是PDB

你先把ambertools装好,然后直接用./acpype.py -i xxx.mol2让acpype脚本载入xxx.mol2文件,acpype就会自动调用机子里的antechamber先把xxx.mol2转化成amber格式拓扑文件,然后再转换成GROMACS等格式的。用起来非常方便。
附件里这个acpype肯定是没问题的
acpype.py (142.73 KB, 下载次数 Times of downloads: 15)

评分 Rate

参与人数
Participants 1
eV +1 收起 理由
Reason
HarrisLin + 1 謝謝 我明天試試看

查看全部评分 View all ratings

北京科音自然科学研究中心http://www.keinsci.com)致力于计算化学的发展和传播,长期开办高质量的各种计算化学类培训:初级量子化学培训班中级量子化学培训班高级量子化学培训班量子化学波函数分析与Multiwfn程序培训班分子动力学与GROMACS培训班CP2K第一性原理计算培训班,内容介绍以及往届资料购买请点击相应链接查看。这些培训是计算化学从零快速入门以及进一步全面系统性提升研究水平的高速路!培训各种常见问题见《北京科音办的培训班FAQ》
欢迎加入“北京科音”微信公众号获取北京科音培训的最新消息、避免错过网上有价值的计算化学文章!
欢迎加入人气非常高、专业性特别强的综合性理论与计算化学交流QQ群“思想家公社QQ群”:1号:18616395,2号:466017436,3号:764390338,搜索群号能搜到哪个说明目前哪个能加,合计9000人。北京科音培训班的学员在群中可申请VIP头衔,提问将得到群主Sobereva的最优先解答。
思想家公社的门口Blog:http://sobereva.com(发布大量原创计算化学相关博文)
Multiwfn主页:http://sobereva.com/multiwfn(十分强大的量子化学波函数分析程序)
ResearchGate:https://www.researchgate.net/profile/Tian_Lu
Money and papers are rubbish, get a real life!

63

帖子

0

威望

452

eV
积分
515

Level 4 (黑子)

7#
 楼主 Author| 发表于 Post on 2017-5-9 21:37:29 | 只看该作者 Only view this author
本帖最后由 HarrisLin 于 2017-5-9 21:38 编辑
k64_cc 发表于 2017-5-9 19:27
话说你准备用什么力场?我记得PRODRG只支持GROMOS力场,Antechamber只支持AMBER力场,这两种方案没什么横 ...

因為我是我們實驗室第一個做的(前人都用DS跑MD)
所以打算先弄懂能跑的基本流程
並沒有很堅持要用哪種力場
換力場等進階步驟之後再慢慢研究

本版积分规则 Credits rule

手机版 Mobile version|北京科音自然科学研究中心 Beijing Kein Research Center for Natural Sciences|京公网安备 11010502035419号|计算化学公社 — 北京科音旗下高水平计算化学交流论坛 ( 京ICP备14038949号-1 )|网站地图

GMT+8, 2024-11-27 05:23 , Processed in 0.179233 second(s), 25 queries , Gzip On.

快速回复 返回顶部 返回列表 Return to list