4
受控分子动力学模?/p>
(Steered Molecular Dynamics)
所谓受控分子动力学模拟?/p>
S
teered
M
olecular
D
ynamics
?/p>
SMD
?/p>
,就是指在进?/p>
分子动力学模拟时?/p>
人为地给分子中的某个或某几个原子施加一个假想的外力?/p>
或者人为地
固定某个或某几个原子的位置?/p>
从而控制整个分子的行为?/p>
我们前面进行的动力学模拟都是
通过各种参数设置?/p>
尽量逼近蛋白质分子在溶液体系中的真实状态,
以研究分子的各种性质
和行为?/p>
但受控分子动力学模拟却要用一个假想的外力干扰控制生物大分子的行为?/p>
这是?/p>
什么呢?/p>
(需要查一下更多的
SMD
的应用实例,
详细了解
SMD
究竟可以用于研究什么问题,
才能做出回答?/p>
在下面的例子中,我们将首先固定泛素分子中一个原子的位置,然后用一个假想的?/p>
力牵拉另一个原子。球形的泛素分子会因此而被逐渐拉开,最终成为伸展状态的肽链?/p>
进行
SMD
时,
我们需要用已经完成能量最小化和能量平衡,
达到稳定状态的蛋白结构?/p>
如果用含有扭曲?/p>
拉伸?/p>
变形构象的原始蛋白结构,
我们将无法分清蛋白各部分的运动是?/p>
于蛋白内部的形变张力引起的,还是由于
SMD
实验附加的外力造成的。因此,进行
SMD
?/p>
前必须预先进行一次平衡态分子动力学模拟,获得稳定结构?/p>
这里?/p>
我们使用
2.4
节球状水体分子动力学模拟输出的恢复文?/p>
?/p>
.restart
?/p>
进行
SMD
?/p>
恢复文件输出时,泛素已在水体中完成了能量最小化和能量平衡,达到了稳定状态?/p>
4.1
除去水分?/p>
为了节省计算时间?/p>
我们在进行本?/p>
SMD
之前将除去体系中所有的水分子?/p>
但读者必
须要注意?/p>
在真正进?/p>
SMD
实验的时候决不可以将水分子除去!
1
、打开
VMD
,选择
File
?/p>
New
Molecule
菜单项,载入
common
目录下的文件
ubq_ws.psf
?/p>
不要关闭窗口?/p>
此时窗口
?/p>
Load file for:
?/p>
一项应该显?/p>
?/p>
0:ubq_ws.psf
?/p>
?/p>
再次单击按钮
Browser
,找?/p>
1-2-sphere
目录下的文件
ubq_ws_eq.restart.coor
,载?/p>
该文件。关?/p>
Molecule File Browser
窗口?/p>
现在我们载入?/p>
2.4
节球状水体动力学模拟时输出的恢复文件
(.restart.coor)
?/p>
2
、选择
Extension
?/p>
tk
菜单项,打开
tk
,首先用
cd
命令改变当前目录?/p>
common
目录下,然后输入?/p>
set selprotein [atomselect top protein]
$selprotein writepdb ubq_ww_eq.pdb
这样我们便在
common
文件夹下新建了一?/p>
pdb
文件
ubq_ww_eq.pdb
?/p>
储存已经?/p>
到平衡态的蛋白质分子,但没有水分子?/p>