24小时热门版块排行榜    

查看: 2483  |  回复: 18
当前只显示满足指定条件的回帖,点击这里查看本话题的所有回帖

fanc232

新虫 (小有名气)

[求助] 使用gromacs模拟失真——关键氢键无法形成 已有3人参与

如附件图片所示,图中Tyr侧链上羟基H atom本应当和配体上的羰基O atom形成氢键;这是配体-受体的关键作用力之一。
但实际模拟的时候,Tyr侧链的这个H atom很快转了方向,指向Phe的侧链。
我怀疑是不是Phe侧链上的大π电子云把H atom吸引过去的,加上配体附近那个Zn离子排斥,就把这个H atom弄走了。
但总之,模拟中如果连这个关键作用力都会丢失,我就肯定没办法获得准确的动力学数据。
在此请各位大神指点迷津!
大家有什么想法请尽管说,我现在已经不知道该往哪个方向查资料了。求助求助!

使用gromacs模拟失真——关键氢键无法形成
ask.png
回复此楼
我很喜欢计算机辅助药物设计
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

awaken2013

铜虫 (小有名气)

引用回帖:
6楼: Originally posted by fanc232 at 2015-05-12 08:16:59
谢谢您帮忙!这个起始结构是结晶体。模拟之前保留了结晶水,但起始结构中,这个位置没有结晶水。模拟的过程是,起始结构——能量最小化——限制性nvt——限制性npt——放开production。

看了您的头像我忍不住想 ...

哈哈 金毛可不好养哦
不用NVT,直接NPT,AN模拟退火到从0到你需要的的温度就可以了

如果确定模拟过程中水没有参与到TYR的氢键漂移的话,你可以尝试尝试ff99SB*-ildn力场(可以在gromacs上下载到),或者直接DB下载编译个AMBER14 用amber去跑,ff12是最新也是目前最好amber力场。

QMMM也是另外一个选择。但是不涉及到化学键的生成和断裂,能量计算等,就没必要
9楼2015-05-12 18:16:27
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 19 个回答

jerkwin

专家顾问 (正式写手)

【答案】应助回帖

★ ★
感谢参与,应助指数 +1
fanc232(abinitio代发): 金币+2, 鼓励交流 2015-05-12 08:05:01
这和所用的力场有关, 你的立场中没有氢键的项, 不出现氢键也很正常
2楼2015-05-11 20:37:47
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

fanc232

新虫 (小有名气)

引用回帖:
2楼: Originally posted by jerkwin at 2015-05-11 20:37:47
这和所用的力场有关, 你的立场中没有氢键的项, 不出现氢键也很正常

首先谢谢您回复!我使用的AMBER03力场,自己也没什么把握。我查查它的氢键项。
关键作用力中包含几个氢键,貌似就这个没法形成。我再仔细查查。
我很喜欢计算机辅助药物设计
3楼2015-05-11 21:11:50
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

awaken2013

铜虫 (小有名气)

★ ★
fanc232(abinitio代发): 金币+2, 鼓励交流 2015-05-12 08:05:11
ff03模拟的精确度还可以。
光从你这张图来看不能说明什么。如果你的复合物结构中存在有结晶水,确定是不是与水形成了氢键,我认为这个可能性更大。另外你的蛋白是否是复合物晶体,若不是的话,最好先加个限制势,进行一小段时间的模拟,再放开production
4楼2015-05-12 00:11:36
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 欢迎发来filecode的Mz6后的代码验证其规律 +32 医学老男孩 2026-08-13 74/3700 2026-08-16 21:05 by 医学老男孩
[论文投稿] 岩土工程学报什么时候才能终审完呐 25+3 yeager111 2026-08-10 4/200 2026-08-16 20:42 by tfang
[基金申请] 时间戳又变了8-15 +13 archvillain 2026-08-15 25/1250 2026-08-16 20:13 by zhaosm1982
[基金申请] 2027广东省杰青 +3 奶牛小黑 2026-08-15 6/300 2026-08-16 20:07 by 奶牛小黑
[基金申请] 我的国基提前知道中了,可是同事的操作让我实在接受不了,怎么会有这样的人 +11 家与远方 2026-08-10 16/800 2026-08-16 10:28 by ray43
[基金申请] 有时候,自然基金真的不能太认真 (我的申报经验) +10 majunge000 2026-08-11 12/600 2026-08-16 08:18 by xli1984
[精细化工] 招聘 金属平磨液,抛光液研发工程师 +3 小天0311 2026-08-14 3/150 2026-08-16 07:31 by H9PLUS
[基金申请] 小木虫上这么多卖论文的,真有人买论文么?感觉没必要啊 +12 Tide man 2026-08-10 13/650 2026-08-15 16:34 by 氺木
[基金申请] 关于Filecode分析方法 +10 majunge000 2026-08-10 13/650 2026-08-15 12:26 by hanpeng972
[基金申请] 各位道友,我要去昆明玩几天,回来见。 +7 Tide man 2026-08-14 8/400 2026-08-15 01:11 by arzu_hma
[基金申请] 是这周出结果还是下周出结果? +4 yuleib84 2026-08-11 4/200 2026-08-14 23:05 by lfy8008
[基金申请] filecode +15 documentary 2026-08-10 17/850 2026-08-14 10:08 by kissu88
[基金申请] 重要来源:本周末出结果 +10 瞬息宇宙 2026-08-12 10/500 2026-08-13 15:46 by likettle
[基金申请] Filecode 又变了,巨变 +3 WH3796 2026-08-12 4/200 2026-08-13 14:13 by 小木虫6752397
[基金申请] 结合人工智能,周易传统文化,filecode打分制来了,3分以上希望很大。 +3 Tide man 2026-08-12 4/200 2026-08-13 08:35 by ZJTJZ
[基金申请] 帮忙看看fileCode +7 wwncly 2026-08-10 13/650 2026-08-11 19:36 by 冰心玉壶晴
[基金申请] 为什么网上很多人说本周 12号出结果 +6 瞬息宇宙 2026-08-10 7/350 2026-08-11 19:25 by Tide man
[基金申请] 什么时候出结果,有咨询渠道??? +3 Tide man 2026-08-11 3/150 2026-08-11 17:54 by kudofaye
[基金申请] 确定了,国自然21号放榜 +6 布布和一二 2026-08-10 7/350 2026-08-10 19:15 by 2000zf36392
[基金申请] 国自然结果 +4 Vierhys 2026-08-10 8/400 2026-08-10 15:06 by Vierhys
信息提示
请填处理意见