24小时热门版块排行榜    

北京石油化工学院2026年研究生招生接收调剂公告
查看: 682  |  回复: 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的回帖

我爱问路

新虫 (初入文坛)

引用回帖:
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的回帖
查看全部 4 个回答

qianxiong

铜虫 (小有名气)

【答案】应助回帖

★ ★
感谢参与,应助指数 +1
我爱问路(月只蓝代发): 金币+2, 感谢指导! 2014-07-08 16:32:31
干吗要手动添加COO?
程序自动会添加。
EM不充分,跑动力学约束失败。
戒微博
2楼2014-07-06 15:44:32
已阅   回复此楼   关注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的回帖
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 285求调剂 +5 AZMK 2026-04-02 6/300 2026-04-03 08:41 by yulian1987
[考研] 一志愿陕西师范大学生物学317分 +5 1563日。 2026-04-02 5/250 2026-04-03 06:58 by ilovexiaobin
[考研] 一志愿北京科技大学材料工程085601,求调剂 +11 cdyw 2026-04-02 11/550 2026-04-02 23:06 by 马儿快快地跑
[考研] 309分085801求调剂 +10 学员Gtwj7W 2026-03-31 10/500 2026-04-02 22:42 by yunlongyang
[考研] 296求调剂 +4 sdhu 2026-04-02 4/200 2026-04-02 21:29 by baoball
[考研] 279求调剂 +5 傅文秋 2026-04-02 5/250 2026-04-02 18:10 by 笔落锦州
[考研] 环境工程297分求调剂一志愿杭高院 +15 GENJIOW 2026-03-31 16/800 2026-04-02 17:56 by cyh—315
[考研] 354求调剂 +4 lxb598 2026-03-31 5/250 2026-04-02 09:55 by Jaylen.
[考研] 08工科275分求调剂 +13 AaAa7420 2026-03-31 13/650 2026-04-02 06:19 by 286640313
[考研] 085410 一志愿211 22408分数359求调剂 +3 123456789qw 2026-03-31 4/200 2026-04-02 00:06 by 义文wang
[考研] 291求调剂 +20 Y-cap 2026-03-29 25/1250 2026-04-01 23:49 by 欣喜777
[考研] 303分 0807学硕求调剂 +3 TYC3632 2026-04-01 3/150 2026-04-01 19:24 by lwk2004
[考研] 285求调剂 +5 FZAC123 2026-03-30 5/250 2026-04-01 15:50 by 韩雨涵
[考研] 085410人工智能 初试316分 求调剂 +3 残星拂曙 2026-03-31 3/150 2026-04-01 11:09 by 小熊raider
[考研] 材料与化工272求调剂 +25 阿斯蒂芬2004 2026-03-28 25/1250 2026-03-31 16:27 by hypershenger
[考研] 求收留 +8 1943443204 2026-03-28 8/400 2026-03-31 15:00 by -迷了路啊路
[考研] 105500药学求调剂,一志愿山东大学药学,348分 +3 gr哈哈哈 2026-03-28 3/150 2026-03-30 18:56 by 源_2020
[考研] 293求调剂 +3 末未mm 2026-03-30 5/250 2026-03-30 17:23 by 王保杰33
[考研] 一志愿南京航空航天大学材料学硕求调剂 +3 @taotao 2026-03-28 3/150 2026-03-28 10:26 by JourneyLucky
[考研] 272求调剂 +7 脚滑的守法公民 2026-03-27 7/350 2026-03-27 17:23 by laoshidan
信息提示
请填处理意见