4 受控分子动力学模拟(Steered Molecular Dynamics)
所谓受控分子动力学模拟（Steered Molecular Dynamics，SMD），就是指在进行分子动力学模拟时，人为地给分子中的某个或某几个原子施加一个假想的外力，或者人为地固定某个或某几个原子的位置。从而控制整个分子的行为。我们前面进行的动力学模拟都是通过各种参数设置，尽量逼近蛋白质分子在溶液体系中的真实状态，以研究分子的各种性质和行为。但受控分子动力学模拟却要用一个假想的外力干扰控制生物大分子的行为，这是为什么呢？（需要查一下更多的SMD的应用实例，详细了解SMD究竟可以用于研究什么问题，才能做出回答）
在下面的例子中，我们将首先固定泛素分子中一个原子的位置，然后用一个假想的外力牵拉另一个原子。球形的泛素分子会因此而被逐渐拉开，最终成为伸展状态的肽链。
进行SMD时，我们需要用已经完成能量最小化和能量平衡，达到稳定状态的蛋白结构。如果用含有扭曲、拉伸、变形构象的原始蛋白结构，我们将无法分清蛋白各部分的运动是由于蛋白内部的形变张力引起的，还是由于SMD实验附加的外力造成的。因此，进行SMD之前必须预先进行一次平衡态分子动力学模拟，获得稳定结构。
这里，我们使用2.4节球状水体分子动力学模拟输出的恢复文件（.restart）进行SMD。恢复文件输出时，泛素已在水体中完成了能量最小化和能量平衡，达到了稳定状态。
4.1 除去水分子
为了节省计算时间，我们在进行本次SMD之前将除去体系中所有的水分子。但读者必须要注意：在真正进行SMD实验的时候决不可以将水分子除去！
1、打开VMD，选择File→New Molecule菜单项，载入common目录下的文件ubq_ws.psf。不要关闭窗口，此时窗口“Load file for:”一项应该显示“0:ubq_ws.psf”。再次单击按钮Browser，找到1-2-sphere目录下的文件ubq_ws_eq.restart.coor，载入该文件。关闭Molecule File Browser窗口。
现在我们载入了2.4节球状水体动力学模拟时输出的恢复文件(.restart.coor)。 2、选择Extension → tk菜单项，打开tk，首先用cd命令改变当前目录至common目录下，然后输入：
set selprotein [atomselect top protein] $selprotein writepdb ubq_ww_eq.pdb
这样我们便在common文件夹下新建了一个pdb文件ubq_ww_eq.pdb，储存已经达到平衡态的蛋白质分子，但没有水分子。
3、删除当前分子，但不要关闭VMD。 4.2 恒速牵拉
在这一节我们将使用恒速牵拉进行SMD动力学模拟。进行恒速牵拉时，我们需要首先设定一个假想的原子作为施力原子，这个实际并不存在的原子就叫做假原子（dummy atom）。假原子通过一根假想的弹簧与系统中真实存在的某个原子（称为SMD原子）相连，
牵拉时，假原子以恒定的速度运动，因此称这种牵拉方式为恒速牵拉。实际上被牵拉原子（SMD原子）的速度不是恒定的。SMD原子所受拉力由弹簧的形变和弹性系数决定，满足胡克定律：。
???F??k?x
又有：
?????x?vt?(r?r0)
????x是假想其中F是SMD原子所受的力，
?弹簧的形变量。它等于假原子的位移vt减 ?去SMD原子的位移，其中r0是SMD原子
?图 恒速牵拉的图示。图中SMD原子是的初始位置，r是SMD原子的当前位置。
蓝色的，假原子（dummy atom）是红色的，二者之间有一根假想的弹簧相连。假原子以恒定的速度运动，弹簧发生形变后作用于SMD原子，于是SMD原子开始运动。 因此，SMD原子的受力可以计算出来，然后根据牛顿定律计算它的运动情况。
4.2.1 设定SMD原子和固定原子（Fixed atom）
NAMD使用pdb文件中的B因子一
栏区别哪些原子是被固定的，哪些原子未被固定：如果某个原子在该栏中的对应值是1（非0值即可），那么该原子被固定，在整个动力学模拟过程中，它的空间坐标不会变化；反之如果为0，那么该原子不受影响。 知识链接：PDB文件的格式 使用写字板打开common目录下我们刚刚制作的ubq_ww_eq.pdb，可以看到文件的内容分成几栏（column）： 在上图中，给我们提供有用信息的共有10栏数据，已在图中标注出： （1）表示所指示为原子 （2）该原子序列号 （3）IUPAC标准格式的原子名称 （4）残基名称 （5）残基序列号 （6）原子的X坐标 （7）原子的Y坐标 （8）原子的Z坐标 （9）位置（Occupancy） （10）温度因子或B因子（beta factor） 其中后两栏在动力学模拟的时候被NAMD用于标记两类特殊的原子：位置（Occupancy）栏用于标记SMD原子，B因子一栏用于标记固定原子（fixed atom）。
4.2.1 设定固定原子（fixed atom）和SMD原子
下面我们将使用VMD创建我们所需的pdb文件。我们在本例中将使用pdb文件中B因子一栏（即温度因子）指定哪个原子将被固定，使用位置（occupancy）一栏指定哪个原子作为SMD原子被牵拉。
1、在VMD中选择File→New Molecule 菜单项，单击Browse按钮找到common目录下的文件ubq.psf，载入它。不要关闭Molecule File Browser窗口，注意Load file for 一项应当显示1:ubq.psf。再次单击Browse按钮，载入common目录下的ubq_ww_eq.pdb。关闭Molecule File Browse窗口，在VMD图形窗口中应当可以看到没有水分子的泛素分子（图）。
图 载入除去水分子后的泛素分子 2、在VMD tk中输入以下命令，固定第一个氨基酸残基的α碳： set allatoms [atomselect top all] $allatoms set beta 0
set fixedatom [atomselect top “resid 1 and name CA”] $fixedatom set beta 1
下面是对命令的解释： ?
set allatoms [atomselect top all] 新建了一个变量allatoms，它代表体系中的所有原子。 ? ?
$allatoms set beta 0 所有原子在pdb文件中的B栏（即beta）全部设为0。 set fixedatom [atomselect top “resid 1 and name CA”] 新建一个变量fixedatom，代表第一个氨基酸残基的α碳。 ?
$fixedatom set beta 1 fixedatom所代表的原子——第一个氨基酸残基的α碳的B栏被设置为1。这样NAMD会在分子动力学模拟时保持该原子固定。 同样地，我们还需要设定哪一个原子将被牵拉。前面提到过，被牵拉的原子就叫做SMD原子。对于这一种原子，NAMD使用pdb文件中的“位置”（occupancy）一栏进行标记。同样地，0代表不受影响的原子，标记为1的原子将被牵拉。设置方法仍然是在tk中输入命令进行。
$allatoms set occupancy 0
set smdatom [atomselect top “resid 76 and name CA”] $smdatom set occupancy 1 以下是对上述命令的解释： ?
$allatoms set occupancy 0 当前pdb文件中所有原子的“occupancy”一栏设置为0。 ?
set smdatom [atomselect top “resid 76 and name CA”] 新建了一个变量smdatom，代表第76个氨基酸（即最后一个氨基酸）的α碳。 ?
$smdatom set occupancy 1 变量smdatom所代表的原子的occupancy 一栏设定为1，这样NAMD会在分子动力学模拟时将假想的力作用于该原子上。 完成上述设定之后，我们需要让VMD输出我们设置好的文件。如果当前目录不是common，用cd命令改变当前目录到common中，然后输入：$allatoms writepdb ubq_ww_eq.ref 回车后会在当前目录生成文件ubq_ww_eq.ref，这就是我们进行SMD所需的文件，它记录了泛素的结构，以及哪一个原子作为固定原子（fixed atom），哪一个作为牵拉原子。
不要关闭VMD，将所有窗口最小化。使用写字板打开common目录下的文件ubq_ww_eq.ref，注意打开时“文件类型”依然要选择“所有文件（*.*）”否则将看不到该文件。在文件开头第六行即可找到固定原子——第一个氨基酸甲硫氨酸的α碳。这个原子