使用Multiwfn图形化研究弱相互作用

内容发布更新时间 : 2025/7/4 8:37:15星期一 下面是文章的全部内容请认真阅读。

[图15]

与B3LYP/6-31G*波函数下的散点图相对比,在Promolecule近似下,主要区别是对应位阻效应的spike明显靠右,这说明芳香环中间的电子密度在形成分子后显著降低以减小互斥效应。并且Promolecule的散点图整体分布比较靠下。所以在观看Promolecule近似下的等值面就必须用与之前不同的等值面数值和电子密度屏蔽范围,如果仍只保留

sign(λ_2)ρ在[-0.05,0.05]区域,则最右边的spike将会被截掉很大部分;如果观看的等值面仍为0.5,则y=0.5的直线在散点图中将与一大片范围相交。所以建议保留sign(λ_2)ρ在[-0.1,0.1]区域的点,并观看RDG=0.25的等值面,这样只有这四个spike对应的区域被保留了下来,并且y=0.25的直线只与它们相交。这样得到的填色等值面如图15下方所示,色彩刻度范围设为了[-0.04,0.03],可见和B3LYP/6-31G*下的图没有太大差别,只是芳香环中间的梭形区域变成椭圆。

虽然靠实际电子密度分析时需要保留的sign(λ_2)ρ的范围有比较通用的值[-0.05,0.05],但是在Promolecule近似下情况比较复杂,spike的位置变化范围较大,应保留的范围不好确定,所以最好每次都考察一下散点图确定保留范围。在Multiwfn中默认的保留范围是[-0.1,0.1],可以通过修改settings.ini里的RDGprodens_maxrho参数来改变,当sign(λ_2)ρ小于-RDGprodens_maxrho或者大于RDGprodens_maxrho的点的RDG值将被自动设为100。若此参数设为0.0,则不自动作屏蔽。

7. 用于生物大分子体系

由于在Promolecule密度下就能得到可靠结果,计算函数时速度很快,也不再需要量子化学软件计算波函数,使得这种分析方法又快又方便,能够容易地用于蛋白质、核酸等生物大分子体系。

此例将绘制DNA双螺旋结构的某一部分的RDG填色等值面图。由于这种分析方法所用的分子结构应当尽量接近平衡构型,获得这样的结构最简单的方法就是从RCSB数据库(http://www.rcsb.org/pdb/home/home.do)里面找相应的晶体结构,然后加氢。也可以自行建立DNA结构,这样序列可以自定义,可以用比如Ambertools里的NAB工具实现,之后再做能量最小化。附件里的DNA.pdb是已建好的含10个腺嘌呤-胸腺嘧啶碱基对儿的DNA片段,已通过Amber在ff99SB力场、GB溶剂模型下优化,细节见http://ambermd.org/tutorials/basic/tutorial1/。为了视觉上比较清楚,将要研究的区域是分子结构边缘区域,如下图灰色方框所示。其中黄色圆球所示的是第84、565号原子,将它们的中点作为网格中心并往外延展一定区域就能使网格涵盖感兴趣的范围。如果不清楚延展多少比较合适,可以先设定一个预期的延展范围,然后用粗糙格点计算一下,看看延展范围是否合适,再做高质量格点的计算。

[图16]

操作步骤如下

DNA.pdb //此文件已在压缩包里。假设读入某个pdb文件时提示发现未知元素,说明原子名有问题,应进行修改。 100 //功能100 1 //子功能1

16,14 //Promolecule近似下的sign(λ_2)ρ函数与RDG函数 -10 //设定延展范围 9 //延展9 bohr

7 //通过两原子中心定义网格中心

84,565 //两原子序号为84和565。注意这个序号是在Multiwfn里的序号,如果在pdb里原子序号是连贯的,比如此例,则在Multiwfn里和pdb文件里原子序号是一致的。如果pdb文件里序号不连贯,建议直接输入网格中心的坐标(网格选项中的第6项)。

120,120,120 //由于网格范围较大,所以用较多格点。

>>展开全文<<
12@gma联系客服:779662525#qq.com(#替换为@) 苏ICP备20003344号-4 ceshi