计算化学公社

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

[Gaussian/gview] 使用Multiwfn生成复合物结合部位的电荷信息该怎么选择关键词?

[复制链接 Copy URL]

139

帖子

0

威望

449

eV
积分
588

Level 4 (黑子)

本帖最后由 YuniJ 于 2024-7-1 15:49 编辑

目的:将聚合物共价修饰到蛋白质上
操作:将聚合物反应端和蛋白质N端的连接区域单独抠出来(附件1:peg_asp.mol2),通过vmd进行了加氢操作,并且观测结构无重叠部分,搜索相关化学知识初步判断结构应该没有问题
问题:使用Multiwfn默认关键词计算RESP原子电荷,报错
Optimization stopped.
    -- Number of steps exceeded,  NStep= 138
    -- Flag reset to prevent archiving.

使用关键词:# B3LYP/def2SVP opt=loose scrf(solvent=water)

想请教一下各位老师,是否是因为关键字的原因?如果是的话我应该考虑哪些方面呢?

peg_asp.mol2

2.46 KB, 下载次数 Times of downloads: 2

连接区域mol2

RESP.sh

2.38 KB, 下载次数 Times of downloads: 0

6万

帖子

99

威望

6万

eV
积分
125156

管理员

公社社长

9#
发表于 Post on 2024-7-3 20:13:00 | 只看该作者 Only view this author
YuniJ 发表于 2024-7-3 19:30
请问老师我按照3.4形成PEG反应端-ASP-NME模型体系,拟合了peg和蛋白连接区域的静电,下一步我可以直接删 ...

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

139

帖子

0

威望

449

eV
积分
588

Level 4 (黑子)

8#
 楼主 Author| 发表于 Post on 2024-7-3 19:30:31 | 只看该作者 Only view this author
sobereva 发表于 2024-7-1 17:41
当前结构什么都不是,优化哪怕收敛了也毫无意义
怎么算单个残基的RESP电荷,在下文3.4节给了明确的例子
R ...

请问老师我按照3.4形成PEG反应端-ASP-NME模型体系,拟合了peg和蛋白连接区域的静电,下一步我可以直接删除NME,然后怼回原来的蛋白质结构中吗?

139

帖子

0

威望

449

eV
积分
588

Level 4 (黑子)

7#
 楼主 Author| 发表于 Post on 2024-7-2 08:27:13 | 只看该作者 Only view this author
student0618 发表于 2024-7-1 21:17
sob老師給的帖文有使用sobtop更詳盡的例子 faq5 只是大概的做法

The structure of modified residue,  ...

谢谢您!

886

帖子

4

威望

2063

eV
积分
3029

Level 5 (御坂)

A Student

6#
发表于 Post on 2024-7-1 21:17:54 | 只看该作者 Only view this author
本帖最后由 student0618 于 2024-7-1 21:58 编辑
YuniJ 发表于 2024-7-1 20:16
Thank you for your answer. The non-standard residue is polyethylene glycol (PEGylation). Are you s ...

sob老師給的帖文有使用sobtop更詳盡的例子 faq5 只是大概的做法

The structure of modified residue, how PEG covalently modified the protein N-terminal, its protonation etc 這些不確定就查一查 自己要做simulation的東西長什麼要清楚哦

PEG 的話也要check要用多少repeating units, 現在是否真的足夠?如果有多個repeating unit的話這可能要另外處理哦
敬仰一针见血的指责,厌倦别有用心的赞美。

139

帖子

0

威望

449

eV
积分
588

Level 4 (黑子)

5#
 楼主 Author| 发表于 Post on 2024-7-1 20:20:39 | 只看该作者 Only view this author
sobereva 发表于 2024-7-1 17:41
当前结构什么都不是,优化哪怕收敛了也毫无意义
怎么算单个残基的RESP电荷,在下文3.4节给了明确的例子
R ...

谢谢老师的回答,我想请教一下老师那我应该怎么构建这个结构呢?聚乙二醇的反应端是CH3O,蛋白质N端是ASP,因为不是化学背景所以可能在结构上有些搞不清楚,不好意思老师

139

帖子

0

威望

449

eV
积分
588

Level 4 (黑子)

4#
 楼主 Author| 发表于 Post on 2024-7-1 20:16:39 | 只看该作者 Only view this author
student0618 发表于 2024-7-1 16:37
Your structure has many problems.

Thank you for your answer. The non-standard residue is polyethylene glycol (PEGylation). Are you saying that I need to look up the protonation state of PEGylation and then proceed according to FAQ 5 in sobtop?

6万

帖子

99

威望

6万

eV
积分
125156

管理员

公社社长

3#
发表于 Post on 2024-7-1 17:41:12 | 只看该作者 Only view this author
当前结构什么都不是,优化哪怕收敛了也毫无意义
怎么算单个残基的RESP电荷,在下文3.4节给了明确的例子
RESP拟合静电势电荷的原理以及在Multiwfn中的计算
http://sobereva.com/441http://bbs.keinsci.com/thread-10880-1-1.html
北京科音自然科学研究中心http://www.keinsci.com)致力于计算化学的发展和传播,长期开办极高质量的各种计算化学类培训:初级量子化学培训班中级量子化学培训班高级量子化学培训班量子化学波函数分析与Multiwfn程序培训班分子动力学与GROMACS培训班CP2K第一性原理计算培训班,内容介绍以及往届资料购买请点击相应链接查看。这些培训是计算化学从零快速入门以及进一步全面系统性提升研究水平的高速路!培训各种常见问题见《北京科音办的培训班FAQ》
欢迎加入北京科音微信公众号获取北京科音培训的最新消息,并避免错过网上有价值的计算化学文章!
欢迎加入人气极高、专业性特别强的理论与计算化学综合交流群思想家公社QQ群(群号见此链接),合计达一万多人。北京科音培训班的学员在群中可申请VIP头衔,提问将得到群主Sobereva的最优先解答。
思想家公社的门口Blog:http://sobereva.com(发布大量原创计算化学相关博文)
Multiwfn主页:http://sobereva.com/multiwfn(十分强大、极为流行的量子化学波函数分析程序)
Google Scholar:https://scholar.google.com/citations?user=tiKE0qkAAAAJ
ResearchGate:https://www.researchgate.net/profile/Tian_Lu

886

帖子

4

威望

2063

eV
积分
3029

Level 5 (御坂)

A Student

2#
发表于 Post on 2024-7-1 16:37:50 | 只看该作者 Only view this author
本帖最后由 student0618 于 2024-7-1 16:58 编辑

Your structure has many problems.


If you know the name of this modified residue but you don't know what the problem is, try to google and check its structure and protonation.

If it is a modified residue within a protein, you also need to cap the two sides. Or if its modified terminal, you may still need to cap one end. Search for "非标准残基" in http://sobereva.com/soft/Sobtop/



(edited post to fix problems in figures)


敬仰一针见血的指责,厌倦别有用心的赞美。

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

GMT+8, 2026-2-25 02:08 , Processed in 0.255792 second(s), 24 queries , Gzip On.

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