为了方便使用可视化软件查看模拟的结果,经常需要将结果合并成pdb等文件,然后用shell远程下载文件后查看。
用trjconv 命令合并生产文件,-s 输入tpr能量文件,-f 输入xtc文件,最后-o表示输出文件。-pbc 后面的
是表示读取里面的数据,后面需要选择,命令如下:
trjconv -s dppc-md.tpr -f dppc-md.xtc -o dppc-md.pdb -pbc whole -conect
:-) trjconv (-: Option Filename Type Description ------------------------------------------------------------ -f dynamic.xtc Input Trajectory: xtc trr trj gro g96 pdb cpt -o dynamic.pdb Output Trajectory: xtc trr trj gro g96 pdb -s dynamic.tpr Input, Opt! Structure+mass(db): tpr tpb tpa gro g96 pdb -n index.ndx Input, Opt. Index file -fr frames.ndx Input, Opt. Index file -sub cluster.ndx Input, Opt. Index file -drop drop.xvg Input, Opt. xvgr/xmgr file Option Type Value Description ------------------------------------------------------ -[no]h bool no Print help info and quit -[no]version bool no Print version info and quit -nice int 19 Set the nicelevel -b time 0 First frame (ps) to read from trajectory -e time 0 Last frame (ps) to read from trajectory -tu enum ps Time unit: fs, ps, ns, us, ms or s -[no]w bool no View output .xvg, .xpm, .eps and .pdb files -xvg enum xmgrace xvg plot formatting: xmgrace, xmgr or none -skip int 1 Only write every nr-th frame -dt time 0 Only write frame when t MOD dt = first time (ps) -[no]round bool no Round measurements to nearest picosecond -dump time -1 Dump frame nearest specified time (ps) -t0 time 0 Starting time (ps) (default: don't change) -timestep time 0 Change time step between input frames (ps) -pbc enum whole PBC treatment (see help text for full description): none, mol, res, atom, nojump, cluster or whole -ur enum rect Unit-cell representation: rect, tric or compact -[no]center bool no Center atoms in box -boxcenter enum tric Center for -pbc and -center: tric, rect or zero -box vector 0 0 0 Size for new cubic box (default: read from input) -clustercenter vector 0 0 0 Optional starting point for pbc cluster option -trans vector 0 0 0 All coordinates will be translated by trans. This can advantageously be combined with -pbc mol -ur compact. -shift vector 0 0 0 All coordinates will be shifted by framenr*shift -fit enum none Fit molecule to ref structure in the structure file: none, rot+trans, rotxy+transxy, translation, transxy or progressive -ndec int 3 Precision for .xtc and .gro writing in number of decimal places -[no]vel bool yes Read and write velocities if possible -[no]force bool no Read and write forces if possible -trunc time -1 Truncate input trajectory file after this time (ps) -exec string Execute command for every output frame with the frame number as argument -[no]app bool no Append output -split time 0 Start writing new file when t MOD split = first time (ps) -[no]sep bool no Write each frame to a separate .gro, .g96 or .pdb file -nzero int 0 If the -sep flag is set, use these many digits for the file numbers and prepend zeros as needed -dropunder real 0 Drop all frames below this value -dropover real 0 Drop all frames above this value -[no]conect bool yes Add conect records when writing .pdb files. Useful for visualization of non-standard molecules, e.g. coarse grained ones
And the next is to choose which part you want to fill :
Reading file dynamic.tpr, VERSION 4.6.5 (single precision)
Select group for output
Group 0 ( System) has 906 elements
Group 1 ( Protein) has 163 elements
Group 2 ( Protein-H) has 163 elements
Group 3 ( C-alpha) has 0 elements
Group 4 ( Backbone) has 0 elements
Group 5 ( MainChain) has 0 elements
Group 6 ( MainChain+Cb) has 0 elements
Group 7 ( MainChain+H) has 0 elements
Group 8 ( SideChain) has 163 elements
Group 9 ( SideChain-H) has 163 elements
Group 10 ( Prot-Masses) has 163 elements
Group 11 ( non-Protein) has 743 elements
Group 12 ( Other) has 743 elements
Group 13 ( W) has 743 elements
Select a group:
Choose "System" mean everything of this design.
可以根据自己的需要进行选择生成。0代表系统里面的全部都转换。