24小时热门版块排行榜    

北京石油化工学院2026年研究生招生接收调剂公告
查看: 683  |  回复: 3
当前只显示满足指定条件的回帖,点击这里查看本话题的所有回帖

我爱问路

新虫 (初入文坛)

[求助] 蛋白质CTER中两个氧原子异常接近 已有1人参与

手动在蛋白质的C端加了一个CTER,但是在跑分子动力学的时候CTER中的两个氧原子会靠的非常近(0.271埃),进程会因为这一段上的某些原子无法束缚而报错终止。错误信息如下:
ERROR: Constraint failure in RATTLE algorithm for atom 24416!
ERROR: Constraint failure; simulation has become unstable.

我所用的CTER如下:
PRES CTER -1.00 ! standard C-terminus
GROUP ! use in generate statement
ATOM C CC 0.34 ! OT2(-)
ATOM OT1 OC -0.67 ! /
ATOM OT2 OC -0.67 ! -C
DELETE ATOM O ! \\
BOND C OT2 ! OT1
DOUBLE C OT1
IMPR C CA OT2 OT1
ACCEPTOR OT1 C
ACCEPTOR OT2 C
IC N CA C OT2 0.0000 0.0000 180.0000 0.0000 0.0000
IC OT2 CA *C OT1 0.0000 0.0000 180.0000 0.0000 0.0000
其中OC 在VDW中的参数如下:
OC 0.000000 -0.120000 1.700000 ! ALLOW POL ION
                ! JG 8/27/89

不知道大家有没有遇到类似的问题,如何解决?
回复此楼

» 猜你喜欢

» 本主题相关商家推荐: (我也要在这里推广)

已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

qianxiong

铜虫 (小有名气)

【答案】应助回帖

★ ★ ★
我爱问路(月只蓝代发): 金币+2, 鼓励交流! 2014-07-08 16:32:37
我爱问路: 金币+1, ★★★★★最佳答案, 感谢您的解答,确实是角度问题,现已解决,多谢!!! 2014-07-27 01:08:44
charmm 力场  默认给出COO与PERS CTER的参数是一样的,

这和电量有吗关系,它们是angle关系,又不计算它们之间的coul 作用
戒微博
4楼2014-07-07 00:16:05
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 4 个回答

qianxiong

铜虫 (小有名气)

【答案】应助回帖

★ ★
感谢参与,应助指数 +1
我爱问路(月只蓝代发): 金币+2, 感谢指导! 2014-07-08 16:32:31
干吗要手动添加COO?
程序自动会添加。
EM不充分,跑动力学约束失败。
戒微博
2楼2014-07-06 15:44:32
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

我爱问路

新虫 (初入文坛)

引用回帖:
2楼: Originally posted by qianxiong at 2014-07-06 15:44:32
干吗要手动添加COO?
程序自动会添加。
EM不充分,跑动力学约束失败。

你好,感谢你的回复。
我的系统分为两部分,一部分是蛋白,另一部分是一种生物大分子。对于蛋白是需要加COO的,但对于生物大分子则不需要。
所以,我把原来程序中的默认加CTER和NTER的语句注释掉了,手动给每个蛋白的两端加CTER和NTER。但是效果和程序自动添加是一样的。如果程序默认添加的话生物大分子也会被加上CTER和NTER。

请问常用的默认CTER和我上面给出的是一样的吗?有没有替代的方法?
现在是在优化的过程中就会出现两个氧原子靠的太近的情况,我在尝试改CTER中原子的电荷量,不知道是否可行。
3楼2014-07-06 19:11:43
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 专硕 351 086100 也是考的材科基 本科也是材料 +7 202451007219 2026-04-02 7/350 2026-04-03 08:06 by 岸上的一条鱼
[考研] 一志愿华中农业071010,总分320求调剂 +6 困困困困坤坤 2026-04-02 6/300 2026-04-02 21:28 by dongzh2009
[考研] 一志愿武汉理工0856,初试334 +3 26考研材料 2026-04-02 3/150 2026-04-02 21:22 by dongzh2009
[考研] 调剂 +3 好好读书。 2026-04-01 6/300 2026-04-02 15:49 by liumengping
[考研] 一志愿北京科技大学085601材料工程英一数二初试总分335求调剂 +8 双马尾痞老板2 2026-04-02 9/450 2026-04-02 14:45 by 5896
[考研] 学硕化学工程与技术,一志愿中国海洋大学320+求调剂 +8 披星河 2026-04-02 8/400 2026-04-02 14:12 by oooqiao
[考研] 一志愿北京科技大学材料学硕328分求调剂 +6 1段时间 2026-03-31 7/350 2026-04-02 13:57 by 3041
[考研] 能源动力 调剂 +3 不破不立0 2026-04-02 3/150 2026-04-02 12:46 by ffffjjjj
[考研] 311求调剂 +9 勇敢的小吴 2026-04-02 9/450 2026-04-02 11:37 by Sammy2
[考研] 266分,一志愿电气工程,本科材料,求材料专业调剂 +10 哇呼哼呼哼 2026-04-01 11/550 2026-04-02 11:31 by lnilvy
[考研] 材料考研调剂 +9 Gs大王 2026-04-02 9/450 2026-04-02 10:24 by olim
[考研] 304求调剂 +12 素年祭语 2026-03-31 15/750 2026-04-01 22:41 by peike
[考研] 273求调剂 +19 李芷新1 2026-03-31 19/950 2026-04-01 21:49 by chyhaha
[考研] 求调剂0703 +5 周嘉尧 2026-03-31 8/400 2026-04-01 20:32 by ltltkkk
[考研] 085600 一志愿9 总分351 求调剂学校 +7 czhcz 2026-03-31 9/450 2026-04-01 19:24 by 无际的草原
[考研] 291求调剂 +3 迷蒙木木 2026-04-01 4/200 2026-04-01 11:07 by 逆水乘风
[考研] 一志愿南昌大学324求调剂 +6 hanamiko 2026-03-29 6/300 2026-03-31 16:35 by hypershenger
[基金申请] 面上5B能上会吗? +8 redcom 2026-03-29 8/400 2026-03-31 15:53 by niuailing
[考研] 266求调剂 +3 哇呼哼呼哼 2026-03-29 3/150 2026-03-31 10:06 by cal0306
[考研] 262求调剂 +7 ZZ..000 2026-03-30 8/400 2026-03-31 10:05 by cal0306
信息提示
请填处理意见