计算化学公社
标题: 使用MolAICal计算纳米管和蛋白质孔道半径的教程 [打印本页]
作者Author: MolAICal 时间: 2020-8-12 00:48
标题: 使用MolAICal计算纳米管和蛋白质孔道半径的教程
本帖最后由 MolAICal 于 2020-9-2 23:15 编辑
使用MolAICal计算纳米管和蛋白质孔道半径的教程
更多教程（含英文教程）请见如下：
1.简介
本教程介绍使用
MolAICal （https://doi.org/10.1093/bib/bbaa161）计算纳米管和蛋白质半径的方法。共分为三个部分：纳米管半径计算，蛋白质孔道半径计算和肽通道半径的计算。最后一个教程是肽通道半径的测量，若果你熟悉
VMD和
NAMD，可以使用由
CHARMM力场产生的
PDB和
PSF文件来测量肽通道的半径，当然也可以只使用肽段的
PDB文件进行肽通道半径的测量。
2.工具
2.1. 所需软件下载地址
2.2. 操作示例文件
所有用到的操作教程文件均可在下面的网站下载：
3.操作流程
3.1. 纳米管半径计算
1)在VMD软件中构建纳米管：Extensions-->Modeling-->Nanotube Builder (如图1)。
(, 下载次数 Times of downloads: 110)
图1.构建纳米管
2）使用类似图1的方法打开VMD Tk Console: Extensions-->Tk Console。在Tk Console中使用cd命令切换到含有文件“nanotube.pdb”和“parameter.dat”的目录，例如：
#>cd d:/005-radiiCal/nanotube
3）使用如下命令在Tk Console中保存生成的纳米管文件，并命名为“nanotube.pdb”
#>set all [atomselect top all]
#>$all writepdb nanotube.pdb
4)选择已构建纳米管内的任意点。你可以在纳米管孔道内表面选择2个不同的原子，然后将这两个原子的连线中心作为“cpoint”。在本教程中，选择的点坐标为-0.2015 0.4185 30.147。打开“005-radiiCal\nanotube”文件夹中的“parameter.dat”，将以上所选点坐标添加到“cpoint”。然后按下文所示修改 “vector”参数：
-----------------------------------------------------------------
cpoint -0.2015 0.4185 30.147
vector 0.00 0.00 1.00
-----------------------------------------------------------------
“0.00 0.00 1.00” 表示沿着Z轴方向的半径测量。“0.00 1.00 0.00” 表示沿着Y轴方向的半径测量。“1.00 0.00 0.00” 表示沿着X轴方向的半径测量。 通道可大致沿任意轴方向放置，即和vector的方向大致一致。
5)在Windows DOS或Linux console中运行如下命令计算半径：
#> molaical.exe -channel radii -cpp parameter.dat
命令运行会生成 “channel_radii.dat”,“dot.vmd_plot” 和“surf.vmd_plot”文件。“dot.vmd_plot”和“surf.vmd_plot”可通过VMD软件展示通道表层。类似图1的方式打开VMD Tk Console： Extensions-->Tk Console。然后运行以下命令:
#> source dot.vmd_plot
本教程省略了纳米管卡通图的做法，你可以根据自己的偏好自行设置。你将看到如图2所示的通道点曲面：
(, 下载次数 Times of downloads: 116)
图2.纳米管通道表面
文件 “channel_radii.dat” 包含了反应坐标和半径值。文件“channel_radii.dat”中的第一列是反应坐标，第二列是半径值。可以使用 OriginLab, Microsoft Excel等工具将其绘制成图。绘制结果如图3所示:
(, 下载次数 Times of downloads: 118)
图3.半径对反应坐标作图结果
3.2. 蛋白孔道半径计算
转至教程所在目录下：
#> cd 005-radiiCal/KcsA
在蛋白通道中选择任意点，然后在参数文件“parameter.dat”中，将参数“cpoint”设置成为所选点的坐标。按照下文所示设置参数“cpoint”和“vector”：
-----------------------------------------------------------------
cpoint 0.001 0.006 1.927
vector 0.00 0.00 1.00
-----------------------------------------------------------------
“0.00 0.00 1.00” 表示沿着Z轴方向的半径测量。“0.00 1.00 0.00” 表示沿着Y轴方向的半径测量。“1.00 0.00 0.00” 表示沿着X轴方向的半径测量。通道可大致沿任意轴方向放置，即和vector的方向大致一致。
1）在Windows DOS或Linux console中运行如下命令计算半径：
#> molaical.exe -channel radii -cpp parameter.dat
2）本运算也会生成“channel_radii.dat”, “dot.vmd_plot”和 “surf.vmd_plot”文件。类似图1的方式打开VMD Tk Console：Extensions-->Tk Console。然后运行以下命令:
#> mol load pdb KcsA.pdb
#> source dot.vmd_plot
本教程省略了蛋白卡通图的做法，你可以根据自己的偏好自行设置。你将看到图4所示点曲面：
(, 下载次数 Times of downloads: 123)
图4.蛋白通道的点曲面
文件 “channel_radii.dat” 包含了反应坐标和半径值。文件“channel_radii.dat”中的第一列是反应坐标，第二列是半径值。可以使用 OriginLab, Microsoft Excel等工具将其绘制成图。半径绘制如图5所示:
(, 下载次数 Times of downloads: 123)
图5.半径对反应坐标作图结果
注意事项：文件“parameter.dat”中的参数“conpar”是一个控制参数，参数“conpar”的默认值是0.15，其数值的增加可以提升随机测量的概率。在这种情况下，有可能出现一些奇怪的测量路线。如果你的蛋白孔道比较规则且大致沿着X, Y, Z轴的某一个方向，你可以减少参数“conpar”的数值，比如设置成0.04，但是参数“conpar”不能设置成0，这样你就能得到较为规整的测量通道。
3.3. 半径计算的高级教程
本部分示例利用由 CHARMM 力场产生的PDB和PSF文件计算多肽孔道半径。转至教程所在目录：
#> cd 005-radiiCal/GramicidinA
选择多肽孔道中的任意点。将参数“cpoint”设置为所选任意点的坐标。按照下文所示设置参数“pdbpath”, “psfpath”,“cpoint” 和 “vector” :
-----------------------------------------------------------------
pdbpath 1JNO.pdb
psfpath 1JNO.psf
cpoint 0.1625 -0.629 -1.838
vector 0.00 0.00 1.00
-----------------------------------------------------------------
“0.00 0.00 1.00”表示沿着Z轴方向的半径测量。“0.00 1.00 0.00”表示沿着Y轴方向的半径测量。“1.00 0.00 0.00”表示沿着X轴方向的半径测量。通道应大致沿任意轴方向放置，即和vector的方向大致一致。
1）在Windows DOS或Linux console中运行如下命令计算半径：
#> molaical.exe -channel radii -cpp parameter.dat -fc charmm
2）本次运算也会生成 “channel_radii.dat”,“dot.vmd_plot” 和 “surf.vmd_plot”文件。类似图1的方式打开VMD Tk Console：Extensions-->Tk Console。运行以下命令:
#> mol load pdb 1JNO.pdb
#> source surf.vmd_plot
本教程省略了多肽卡通图的做法，你可以根据自己的偏好自行设置。你将看到图6所示多肽通道 (如图6所示):
(, 下载次数 Times of downloads: 129)
图6.多肽通道
文件 “channel_radii.dat” 包含了反应坐标和半径值。文件“channel_radii.dat”中的第一列是反应坐标，第二列是半径值。可以使用 OriginLab, Microsoft Excel等工具将其绘制成图。半径绘制如图7所示：
(, 下载次数 Times of downloads: 107)
图7.半径对反应坐标作图结果
作者Author: Eva_Winter 时间: 2021-8-16 22:00
请问下，对于非对称的蛋白，孔道每个位置的横截面形状都不一样，也可能不是圆，那么这里的半径值是将孔道近似为圆得到的值吗？还是取的孔半径最小值？
作者Author: amber菜鸟 时间: 2021-8-28 11:16
非常Nice的软件，请问这款软件在windos环境下如何安装呢？搜了很多发现教程非常详细，但是安装教程确一致没有找到。
作者Author: MolAICal 时间: 2021-10-17 02:16
windows下直接解压就可以用，不需要安装
安装的方式看一下这个： https://molaical.gitee.io/install.html
作者Author: MolAICal 时间: 2021-11-10 15:37
近似值，采取蒙特卡洛退火模拟的算法获取最优值
作者Author: zuozuo 时间: 2024-2-2 22:14
3.2. 蛋白孔道半径计算
转至教程所在目录下：
#> cd 005-radiiCal/KcsA
在蛋白通道中选择任意点，然后在参数文件“parameter.dat”中，将参数“cpoint”设置成为所选点的坐标。按照下文所示设置参数“cpoint”和“vector”：
-----------------------------------------------------------------
cpoint 0.001 0.006 1.927
您好，请问在计算蛋白通道半径的时候，“转至教程所在目录下：#> cd 005-radiiCal/KcsA”和参数文件“parameter.dat”在哪里可以找到呀，我下载出来的压缩包里没有这些文件
作者Author: MolAICal 时间: 2024-2-7 16:24
你看一下这个链接：https://github.com/MolAICal/tutorials/tree/master/005-radiiCal/KcsA 或者 https://gitee.com/molaical/tutorials/tree/master/005-radiiCal/KcsA 你可能下错压缩包，或者没有下载完全，parameter.dat就在里面啊，你在检查一下。实在不行，你加一下MolAICal讨论的QQ群 1151656349，我传给你
作者Author: ZZZephyr 时间: 2025-8-7 13:21
你好，请问我在DOS输入计算半径指令时显示“此应用无法在你的电脑上运行”是什么问题呀？版本不兼容吗
| 欢迎光临 计算化学公社 (http://ccc.keinsci.com/) |
Powered by Discuz! X3.3 |