NAMD入门教程三

1、4受控分子动力学模拟 (Steered Molecular Dynamics)所谓受控分子动力学模拟( Steered M olecular Dynamics , SMD ),就是指在进行 分子动力学模拟时，人为地给分子中的某个或某几个原子施加一个假想的外力，或者人为地固定某个或某几个原子的位置。从而控制整个分子的行为。我们前面进行的动力学模拟都是通过各种参数设置， 尽量逼近蛋白质分子在溶液体系中的真实状态，以研究分子的各种性质和行为。但受控分子动力学模拟却要用一个假想的外力干扰控制生物大分子的行为，这是为什么呢？(需要查一下更多的 SMD的应用实例，详细了解SMD究竟可以用于研究什么问题，

2、才能做出回答)在下面的例子中，我们将首先固定泛素分子中一个原子的位置，然后用一个假想的外 力牵拉另一个原子。球形的泛素分子会因此而被逐渐拉开，最终成为伸展状态的肽链。进彳t SMD时，我们需要用已经完成能量最小化和能量平衡，达到稳定状态的蛋白结构。如果用含有扭曲、拉伸、变形构象的原始蛋白结构，我们将无法分清蛋白各部分的运动是由 于蛋白内部的形变张力引起的，还是由于SMD实验附加的外力造成的。因此，进行 SMD之前必须预先进行一次平衡态分子动力学模拟，获得稳定结构。这里，我们使用2.4节球状水体分子动力学模拟输出的恢复文件(.restart )进行SMD。恢复文件输出时，泛素已在水体中完成了能量

3、最小化和能量平衡，达到了稳定状态。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

4、 窗口。现在我们载入了 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动力学模拟。进行恒速牵拉时，我

5、们需要首先dummySMD原子）相连,设定一个假想的原子作为施力原子，这个实际并不存在的原子就叫做假原子（atom ）。假原子通过一根假想的弹簧与系统中真实存在的某个原子（称为牵拉时，假原子以恒定的速度运动，因此 称这种牵拉方式为恒速牵拉。实际上被牵 拉原子（SMD原子）的速度不是恒定的。SMD原子所受拉力由弹簧的形变和弹性系数决定，满足胡克定律:LT F喻4又有:VLr vtr (rrTo)ur其中F是SMD原子所受的力,X是假想弹簧的形变量。它等于假原子的位移 Vt减r去SMD原子的位移，其中r0是SMD原子图恒速牵拉的图示。图中SMD原子是蓝色的，假原子（dummyatom ）是红色的，

6、二者之间有一根假想的弹簧相连。假原的初始位置，r是SMD原子的当前位置。因此，SMD原子的受力可以计算出 来，然后根据牛顿定律计算它的运动情况。子以恒定的速度运动，弹簧发生形变后作用4.2.1设定SMD原子和固定原子于SMD原子，于是SMD原子开始运动。(Fixed atom )NAMD使用pdb文件中的 B因子1 (非栏区别哪些原子是被固定的，哪些原子未被固定：如果某个原子在该栏中的对应值是o值即可），那么该原子被固定，在整个动力学模拟过程中，它的空间坐标不会变化；反之 如果为0,那么该原子不受影响。知识链接：PDB文件的格式使用写字板打开common目录下我们刚刚制作的ubq_ww_eq.

7、pdb ,可以看到文件的内容分成几栏（column ）：M阍固.MSB0曾nbq >w eq -写字板CRYST10.0000.0000,000 90.00 00.0090,00P 11ATOffl1NJ127.70124,3583.冏1.000.0DuATOm2HTlMETJ127.29123.4V93. 041kOO0.00uATOffi3HT2METu128.24224,2254.033hOO0.00uATOM4HT3METIT128, 34324, 4972. 348kOO0.00uATOM5CAMETIT126, 59925. 2fi03. 2521.000.00uAiorn6

8、HAMETu126, 24925. seo2. 2851.000.00uATOM7CBTiffiTJ125.3552% 6453. 99。1,000,0 口uATOffl8HEIMETJ124.53825.3S3工364LOO0.00uATOMgHB2METJ125.09723.7113,4451.000.00uATOfll£0OGMETIT125. S8624. 3935.423kOO0.00uATOJI11HGLHETIT126. 5S823, 9£25. 6141. 000.00uATOM12HGZTT：J125 68225,4155,359LOO0.00uATOM1

9、33DMETJ124.36923,466. 1561.000.00uATOffi14C£METUi24.T2T23, S23LSI 21.000.00uATOHl15HEkMETITi25, 53222. 9958, 2391.000.00uATOM16METIT 川一筌里14K3%,(3)(4)(6)(7)(8)(9)(1t>)文件褊辑直看插入(1)辂式帮助此在上图中，给我们提供有用信息的共有10栏数据，已在图中标注出：(1 )表示所指示为原子(2)该原子序列号(3) IUPAC标准格式的原子名称(4)残基名称(5)残基序列号(6)原子的X坐标(7)原子的Y坐标(8)原子的Z

10、坐标(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 菜单项，单击 Brows

11、e 按钮找到common 目 录下的文件ubq.psf ,载入它。不要关闭 Molecule File Browser 窗口，注意Load file for 一项应当显示1:ubq.psf 。再次单击Browse 按钮，载入common 目录下的 ubq_ww_eq.pdb 。关闭Molecule File Browse 窗口，在VMD图形窗口中应当可以看到没 有水分子的泛素分子(图)。图载入除去水分子后的泛素分子-I匚| X|H Vllf LELS 抑 eiiGL Ui splay2、在VMD tk中输入以下命令，固定第一个氨基酸残基的a碳:set allatoms atomselect t

12、op allresid 1 and name CA”$allatoms set beta 0set fixedatom atomselect top$fixedatom set beta 1卜面是对命令的解释:set allatoms atomselect top all新建了一个变量 allatoms , 它代表体系中的所有原子。$allatoms set beta 0 所有原子在pdb文件中的B栏（即beta ）全部设为0。set fixedatom atomselect top伺sid 1 and name CA” 新建一个变量fixedatom ，代表第一个氨基酸残基的a碳。$fixe

13、datom set beta 1 fixedatom 所代表的原子第一个氨基酸残基的a碳的B栏被设置为1。这样NAMD会在分子动力学模拟时保持该原子固定。同样地，我们还需要设定哪一个原子将被牵拉。前面提到过，被牵拉的原子就叫做 SMD原子。对于这一种原子，NAMD使用pdb文件中的"位置"（occupancy ） 一栏进行标记。同样地，0代表不受影响的原子，标记为1的原子将被牵拉。设置方法仍然是在tk中输入命令进行。$allatoms set occupancy 0set smdatom atomselect top“ resid 76 and name CA ”$smda

14、tom set occupancy 1以下是对上述命令的解释：$allatoms set occupancy 0 当前 pdb 文件中所有原子的 “ occupancy ” 栏设置为0。set smdatom atomselect top“resid 76 and name CA ” | 新建了一个变量smdatom ，代表第76个氨基酸（即最后一个氨基酸）的a碳。$smdatom set occupancy 1 变量 smdatom 所代表的原子的 occupancy 一栏设定为1 ,这样NAMD会在分子动力学模拟时将假想的力作用于该原子上。完成上述设定之后，我们需要让VMD输出我们设置好的

15、文件。如果当前目录不是common ,用 cd命令改变当前目录至Ucommon 中，然后输入：$allatoms writepdbubq_ww_eq.ref回车后会在当前目录生成文件ubq_ww_eq.ref ,这就是我们进行 SMD所需的文件，它记录了泛素的结构，以及哪一个原子作为固定原子（ fixed atom ）,哪一个 作为牵拉原子。不要关闭VMD ,将所有窗口最小化。使用写字板打开common目录下的文 件ubq_ww_eq.ref ,注意打开时“文件类型”依然要选择“所有文件（*.* ）”否则将看不到该文件。在文件开头第六行即可找到固定原子一一第一个氨基酸甲硫氨酸的“碳。这个原子的

16、B因子一栏为1.00 ,而其它原子都为0.00 (图)。拖动滚动条向下，一直到倒数第四行，可以看到我们定义的 SMD原子一一最后一个氨 基酸甘氨酸的a碳。这个原子的Occupancy 一栏是1.00 ,而其它原子都是0.00 (图)。需要说明一点：在上面两张图中出现的原子坐标数值(第 7、8、9栏数值)可能和读 者看到的不一致。这是因为动力学模拟中，能量平衡( equilibration ) 一步中各个原子的 初始速度是随机选择的。我们使用的pdb文件是各个原子能量平衡结束后的末速度，因为初速度是随机的，末速度也不一定相同。4.2.2设定拉力的方向现在，关闭写字板。我们已经设定了固定原子( f

17、ixed atom )和SMD原子。下面我 们需要做的是设定拉力的方向。我们定义拉力的方向为沿固定原子( fixed atom )和SMD原子的连线方向（即由固定原子到SMD原子的矢量的方向）。在VMD tk 中输入：set smdpos lindex $smdatom get x y z 0set fixedpos lindex $fixedatom get x y z 0vecnorm vecsub $smdpos $fixedpos输入以上命令后，tk中给出三个数值，这三个值就是固定原子（fixed atom ）到SMD原子的矢量坐标，分别为 x, y, z坐标。记下这三个值（图），我们

18、后面还会用到。然后关 闭VMD即可。图假想的拉力矢量的坐标现在我们已经得到了文件： ubq_ww_eq.ref ,这个文件储存了原子的位置，以及哪个原子是固定原子，哪个原子是SMD原子。下一步我们将制作 NAMD配置文件。注意在下面的步骤中一定要避免输入错误，否则很可能导致动力学模拟无法正常进行。我们将对一个样本配置文件进行修改，制作我们所需要的配置文件。在实际工作中， 读者也可以修改NAMD教程中提供的样本配置文件，得到自己所需要的文件。1、使用Windows资源管理器获得我们所需的样本配置文件sample.conf 。文件在common目录下。将这一文件粘贴到 3-1-pullcv 目录下

19、，然后改名为ubq_ww_pcv.conf 以 和其他文件相区分。2、使用写字板打开 ubq_ww_pcv.conf 。现在我们看到的是一个配置文件的样本 （图）。浏览一下可以发现很多参数是用户需要根据实际情况进行设定的。事实上，这个样 本配置文件就是NAMD动力学模拟配置文件的基本蓝本，用户可以在此基础上量身定做符合 自己需要的配置文件。图NAMD配置文件样本3、在 Job Description一栏力口入：# N- C- Termini Constant Velocity Pulling当然这只是一个注释，不加也不会影响我们的动力学模拟。4、在Adjustable Parameters一部

20、分进行如下改变，以符合我们本次动力学模拟的实际条件：structure mypsf.psf-structure ./common/ubq.psfcoordinates mypdb.pdb-coordinates ./common/ubq_ww_eq.pdboutputName myoutputfoutputName ubq_ww_pcv以上各项的含义参见|2.3.1 配置文干|5、在Input 一部分进行如下改变，以符合我们本次动力学模拟的实际条件:parameters par_all27_prot_lipid.inpf parameters ./common/par_all27_prot_l

21、ipid.inp我们没有使用周期性边界条件，因此不需要设置Periodic Boundary Conditions部分。此夕卜 Force-Field Parameters,Integrator Parameters和 PME三项也不需改变。6、Constant Temperature Control运动：一项应该关闭，因为温度控制可能干扰原子的langevin on f langevin off7、Constant Pressure Control8、Fixed Atoms Constraint不需要进行改变，默认是不启用。需要启用。请对前两行作出如下改变:if 0 - if 1 fixed

22、AtomsFile myfixedatoms.pdbf fixedAtomsFile ./common/ubq wweq.ref(注：和C语言中判断表达式真假类似，if1表示条件成立，if0表示条件不成立)第三行不需要改变，因为我们是使用B值这一栏设定固定原子的。如果在实际工作中使用了其他栏，就需要在这里作出相应改变以通知NAMD那一栏用于标记固定原子。9、IMD Settings一项不需要改变，默认是不启用10、在Extra Parameters一部分输入以下内容：SMDonSMDFile./common/ubq_ww_eq.refSMDk7SMDVel0.005上面的内容的作用是：通知NA

23、MD我们将使用SMD原子。SMD原子和假原子之间的弹簧弹性常数为7 kcal/mol/A 2 (1 kcal/mol = 69.479 pN /A )。SMDVel 一项设定的就是假原子的速度。这里设定为0.005A /timestep 。因为我们设定 Itimestep = 2fs ,因此假原子的速度为2.5A /ps。在实际动力学模拟过程中，假原子将以该恒定速度运动，并通过一个假想的弹簧同 SMD原子相连，牵拉SMD原子运动。11、在Extra Parameters部分我们还需要制定拉力的方向。请输入SMDDirnx ny nz注意nx, ny, n z三个参数就是我们在 4.2.2 一节

24、求出的矢量的坐标。12、接下来继续输入：SMDOutputFreq 10这个值指定每隔多少步输出一次SMD过程产生的数据13、最后，在Execution Script部分将动力学模拟时间改为20000步，相当于40Ps :run 50000 一 run 20000然后保存文件，关闭写字板。我们已经完成了配置文件的设定。提示请再次注意，我们的设置仅仅是为了使得动力学模拟能够迅速完成，因此采取了许多 不恰当的体系设定。比如，具有科学意义的动力学模拟应当至少计算几个ns （几百万步），牵拉速度应该尽可能慢一些，在0.1A/ps左右。并且正如我们提到的，绝不可以忽略水分子。4.2.4 SMD 模拟下面

25、我们可以进行动力学模拟了。为了保证模拟正常进行，请确认 3-1-pullcv 目录下 是否有文件 ubq_ww_pcv.conf , 并确认 common 目录下有以下文件：ubq.psfubq_ww_eq.pdbubq_ww_eq.refpar_all27_prot_lipid.inp下面，打开terminal ,使用cd命令改变当前目录到 namd-tutorial/namd 目录下（即 NAMD所在目录），然后输入：namd2 ./3-1-pullcv/ubq_ww_pcv.conf > ./3-1-pullcv/ubq_ww_pcv.log实测：CPU: Pentium M740

26、 1.73GHz,内存：512M 可以在 10min 之内完成。结果的获得与分析将在4.4节讲解。4.3恒力牵拉下面我们将进行另一种类型的 SMD模拟：恒力牵拉。我们仍需要固定一个原子并指定 另一个原子被牵拉，但是拉力是恒定的，直接作用于SMD原子上。因此不需要指定假原子以及假想弹簧的弹性系数。4.3.1设定固定原子(fixed atom )和SMD原子我们在本例中仍将使用 pdb文件中B因子一栏(即温度因子)指定哪个原子将被固定， 使用位置(occupancy ) 一栏指定哪个原子作为SMD原子被牵拉，并设定拉力的大小。1、在VMD中选择File 一 New Molecule菜单项，单击 B

27、rowse 按钮找到common 目录下的文件ubq.psf ,载入它。不要关闭 Molecule File Browser 窗口，注意Load file for 一项应 当显示0:ubq.psf 。再 次单击Browse 按钮，载入common 目 录下的 ubq_ww_eq.pdb 。关闭Molecule File Browse 窗口，在VMD图形窗口中应当可以看到没 有水分子的泛素分子2、在VMD tk中输入以下命令，固定第一个氨基酸残基的“碳：set allatoms atomselect top all$allatoms set beta 0set fixedatom atomsel

28、ect topresid 1 and name CA ”$fixedatom set beta 13、输入以下命令，设定 SMD原子和拉力的大小：$allatoms set occupancy 0set smdatom atomselect topresid 76 and namd CA”$smdatom set occupancy 11.54这样第76个残基的a碳的 Occupancy 一栏被设定为11.54 ,其他原子都为 0。这 样不仅设定了 SMD原子，而且设定了拉力大小为 11.54 kcal/mol/A °。相当于800pN (8.0 x 10 -10N)4、下面设定拉力

29、的方向。$smdatom set x nx$smdatom set y ny$smdatom set z n注意nx, ny , nz三个参数就是我们在4.2.2 一节求出的矢量的坐标 (0.3671,0.3783,0.8498)。设定好之后图形窗口显示一个极度变形的分子。不过没有关系，这是因为VMD将我们设定的拉力坐标按照原子的坐标显示出来了。5、下面在tk中首先使用cd命令改变当前目录到 common目录下，然后输入$allatoms writepdb ubq_ww_eq2.ref保存我们的设定。然后打开Windows资源管理器，找到common目录下，检查一下是否已经有ubq_ww_eq

30、2.ref 文件？如果有了，就可以关闭VMD。4.3.2 配置文件和上一节我们进行的过程一样，我们还需要修改样本配置文件sample.conf ,制作我 们需要的配置文件。首先打开资源管理器，将 common目录下的sample.conf 拷贝到3-2-pullcf目录下，并改名为ubq_ww_pcf.conf ,然后用写字板打开该文件开始编辑：3、在 Job Description一栏加入：# N- C- Termini Constant Force Pulling当然这只是一个注释，不加也不会影响我们的动力学模拟。4、在Adjustable Parameters一部分进行如下改变，以符合我

31、们本次动力学模拟的实际条件：structure mypsf.psfcoordinates mypdb.pdboutputName myoutput-以上各项的含义参见2J1配置文件5、在 Inputstructure ./common/ubq.psfcoordinates ./common/ubq_ww_eq.pdboutputName ubq_ww_pcf部分进行如下改变，以符合我们本次动力学模拟的实际条件:parameters par_all27_prot_lipid.inpf parameters ./common/par_all27_prot_lipid.inp同上一节，我们不需要设置

32、 Periodic Boundary Conditions一部分。此夕卜Force-FieldParameters,Integrator Parameters 和 PME 三项也不需改变。6、同上一节，Constant Temperature Control一项应该关闭，因为温度控制可能干扰原子的运动：langevin on f langevin off不需要进行改变，默认是不启用。需要启用。请对前两行作出如下改变：- if 1 f fixedAtomsFile ./common/ubq ww默认是不启用7、Constant Pressure Control8、Fixed Atoms Cons

33、traintif 0 fixedAtomsFile myfixedatoms.pdb eq2.ref9、IMD Settings一项不需要改变部分输入以下内容：10、在 Extra Parametersconstantforceyesconsforcefile./common/ubq_ww_eq2.ref上面的内容的作用是：通知NAMD我们将对SMD原子施加一个恒定的拉力。拉力的大O拉力的方向不需要在此指定,小就是Occupancy 一栏我们设定的值 11.54 kcal/mol/A因为我们刚才已经记录到了ubq_ww_eq2.ref文件中。11、最后，在Execution Script部分将

34、动力学模拟时间改为20000步，相当于40Ps :run 50000 一 run 20000然后保存文件，关闭写字板。我们已经完成了配置文件的设定。4.3.3 动力学模拟下面我们可以进行动力学模拟了。为了保证模拟正常进行，请确认 3-2-pullcf目录下 是否有文件 ubq_ww_pcf.conf , 并确认 common 目录下有以下文件：ubq.psfubq_ww_eq.pdbubq_ww_eq2.refpar_all27_prot_lipid.inp下面，打开terminal ,使用cd命令改变当前目录到 namd-tutorial/namd 目录下（即 NAMD所在目录），然后输入：

35、namd2 ./3-2-pullcf/ubq_ww_pcf.conf> ./3-2-pullcf/ubq_ww_pcf.log实测：CPU: Pentium M740 1.73GHz,内存：512M 可以在 10min 之内完成。4.4结果分析4.4.1 恒速拉伸SMD结果分析我们将在本节中分析恒速拉伸SMD的结果。如果读者没有成功完成动力学模拟，可以使用3-1-pullcv/example-output中的结果文件。4.4.1.1 日志文件首先我们看一下输出的日志文件。打开写字板，找到日志文件 3-1-pullcWubq_ww_pcv.log打开它。向下拖动滚动条，是否发现比我们第一次

36、动力学模拟的日志文件多了一些内容？读者可以找到1-2-sphere/ubq_ws_eq.log对照一下增加的内容（图）。置ifbq_ir_puT -写字板Jn文件旧 褊辑旧 查看（V）插入co 格式® 帮助oo口|¥旧|管I® 闯 相嘲电H 曳jETITLE:TSB0MDAJTGLEDIKED0PEK吟 EKEEND切 SYSTEM 工鼬FILESWIt 39. 2904 38. 5471 32.8393.3831 13.7S94 3L 0079, 3933561： 316C38 32. S5 H 2601; W. 1205 22. 7573SfflD3C 39.

37、 2642 38. 6804 32. 813； ?7-8025 28. 635S 64. 3927M 33. 27S1 38. 7534 32. S04： 32.3131 33.2363 74.S512SID51 39.253 38.7697 32.7475 5).3944 5 L 0049 116.717ENERGY; 0233.5451707. 7219411. 1340图输出文件中新增的内容可以看到，在每次输出ENERGY之后，紧接着增加以 SMD开头的几行文字 （图）， 这是NAMD输出的关于SMD原子（即我们定义的第 76氨基酸的a碳）的数据。图中（1）记录的是步数，（2）记录的是输

38、出时 SMD原子的坐标，（3）记录的是输出时 SMD原子所受拉力的大小（x,y,z三个方向上的分力）4.4.1.2 轨迹动画下面我们将看一看动力学模拟生成的原子运动轨迹，直观地感受一下泛素肽链是如何 被拉力拉开的。1、打开 VMD ,选择 File 一 New Molecule菜单项，单击按钮 Browse 找 Uncommon目录下的ubq.psf ,载入该文件； 不要关闭Molecule File Browser 窗口，注意Load file for 一项应当显示 0:ubq.psf 。再次单击 Browse按钮，载入 3-1-pullcv目录下的轨迹文件ubq_ww_pcv.dcd 。关

39、闭Molecule File Browse窗口，在VMD图形窗口中应当可以看到 泛素分子。2、为了更形象的观察蛋白质分子的变化，我们需要改变一下显示模式。选才i Graphic s一 Representation 菜单项，打开显示模式设置，将Drawing Method一项由Lines改为Cartoon （图）以便直观地观察到二级结构的变化。然后调节主窗口右下角的速度滑块至中央，单击最右边的播放按钮，就可以看到整个过程的动画了:-Jni2slVBD lain速度滑块 革上此调节到中部鬻某 始播放loop 为注意Cartoon 模式下二级结构的表示：“螺旋为圆柱，3片层为片层箭头,线绳。但是因为

40、 VMD不会每一帧都将二级结构重新渲染，因此我们看到，在动画的整个过 程中肽链都是线绳状（图）。这是因为我们从最后一帧开始渲染。最后一帧时，整个肽链都 处于伸展状态，相当于全是loop；因此VMD在渲染其它帧的时候也认为整条肽链都是loop,从而以线绳表示3、单击播放控制滑块，将它向右推到第一帧，这时的蛋白还是以loop形式表示的（图）。我们需要让VMD重新渲染一下以正确地看到二级结构。方法是打开tk ,输入命令：mol ssrecalc top图：伸展状态下 VMD认为整个肽链都是loop将播放控制滑块推到开头第一帧图将播放控制滑块推到第一帧ran i .产三的"L叫。.t可以看到

41、VMD重新对二级结构进行了渲染(图)。按照此方法可以逐帧地观察肽链被拉开的过程中二级结构的变化 ，研究一下在拉伸的过程中哪一种二级结构先消失，哪一种后消失。下面显示的是调整其他参数后渲染出的三帧。图中左下灰色圆球为固定原子，右上 角的灰色圆球是SMD原子。知识链接：蛋白质的解折叠途径( unfolding pathway )做到这一步，可能有的读者会问，我们用一个假想的拉力将蛋白质拉开解折叠有什么实际意义呢？事实上，蛋白质的解折叠途径是生命体中很重要的一个过程。我们知道，除少数在线粒体和叶绿体中所合成的蛋白质之外， 游离核糖体和内置网膜上的核糖体所合成的蛋白质，-般都含有分拣信号(sortin

42、g signal ),决定它们的最终去向和定位。分拣型号可以引导蛋白质跨膜转运，将蛋白质从细胞质溶质装移到细胞内的不同部位。在跨膜的过程中，蛋白质必须解折叠，以便于跨膜转运。解折叠的过程是蛋白质信号序列或前导序列一端固定在膜上， 整个肽链signal ）序列。这一信号序列决定了蛋白质的趋向和最终定位（470页）4.4.1.3 牵引力-时间作图下面我们将分析牵引拉力随时间的变化情况，并使用 Excel进行作图。1、在tk中使用cd命令改变当前目录到3-1-pullcv ,然后输入：mkdir analysis这样我们在3-1-pullcv目录下新建了一个文件夹analysis ,用于储存我们分析

43、时生成的数据文件。2、下面我们将使用一个脚本文件从日志文件中提取NAMD输出文件中关于 SMD受力大小的信息。脚本文件就在3-1-pullcv目录下，读者可以输入dir命令查看当前目录下的文件，应该有一个文件为ft.tcl。下面在tk中输入：source ft.tcl这样就调用了脚本文件ft.tcl。调用之后，脚本会要求用户输入三个值：n_x , n_y ,n_z。分别输入我们在4.2.2一节计算得到的拉力矢量坐标（0.3671,0.3783,0.8498）（图）。该脚本最后会在3-1-pullcv/analysis目录下新建文件ft.dat储存提取出的数据。?il« Cons.ol

44、4 Edi t I nt er p Prefs History Helploading history file ,. 43 events addedMain consols display active (Tel8.4.1 / Tk8.4.1)(VND) 49 % cd . ./nt,u/3-1-pullCV>PIain< (3rpu工up) 5 0 % source ft. tcl Enter a value for n x: .3671Elite rvaluefor输入我们定义拄力0.3733Bute 匚valuefor向的矢量的坐标0,6498 n<(3-i-pullcv) 51 %图ft.tcl脚本的调用3、 VMD ,ExcelExcel ,选择菜单项 File 打开，文件类型一栏要选择“所有文件”，找到3-1-pullcv/analysis目录下的文件ft.dat,打开它。在弹出的“文本导入向导”中，先单击“下一项”，然后一定要选中“空格” 一项，再