VMD加载xtb优化或动力学轨迹时读入/保留pdb信息的方法

tjuptz

        xtb可以拿来做蛋白等生物体系，这种体系往往用pdb格式的坐标文件，但用xtb做优化或动力学后的轨迹文件都是xyz格式的，丢失了残基等信息。虽然优化之后还顺带输出了一个优化好的原格式结构文件，但就这个文件可能还有bug，比如把ligand每一个原子都价格TER给分开-_-！所以为了一些特殊目的，比如用VMD做一些分析啥的，方法如下：
1、就像gromacs的轨迹载入一样，先载入一个pdb，再删除第一帧(xtb轨迹应该是不包含输入构象，但用了--cma的话轨迹可能jump)，然后把xtb轨迹文件加载到这个mol中

1. mol new original_input.pdb
   
2. mol addfile xtb.trj type xyz waitfor all

复制代码

2、基于上述可以把整个轨迹都转化成pdb格式，可能对其他可视化程序有用。通过$sel可以自定义输出组、还可以选择输出的帧、频率等。里面还定义了个分帧输出的功能，可自行更改起始帧、终止帧、每隔多少帧输出。

1. # convert xtb opt/md trjactory to pdb file
   
2. 
   
3. # load template pdb and trajectory into one molid
   
4. # if --cma used in xtb, delete frame 0 to avoid jump
   
5. mol new pro_IA.pdb
   
6. animate delete  beg 0 end 0 skip 0 0
   
7. mol addfile xtb.trj type xyz first 0 last -1 step 1 waitfor all 0
   
8. 
   
9. set sel [atomselect 0 all]
   
10. set n [molinfo 0 get numframes]
    
11. animate write pdb trj.pdb beg 0 end [incr n -1] skip 1 waitfor all sel $sel 0
    
12. 
    
13. puts "All Done!"
    
14. 
    
15. # function to separate each frame into individual pdbfile
    
16. proc seppdb {} {
    
17. set sel [atomselect 0 all]
    
18. set n [molinfo 0 get numframes]
    
19. for { set i 0 } { $i < $n } { incr i 50 } {
    
20.     $sel frame $i
    
21.     $sel update
    
22.     $sel writepdb trj_$i.pdb
    
23. }
    
24. }

复制代码

注：VMD1.9.3输出的pdb应该是没有MODEL字段；xtb6.3.3之后已经可以自己输出pdb格式的轨迹了，且含有MODEL字段。另外ligand被TER分割是因为是用了本来给蛋白用的ATOM标引而非HETATM。