|
gmx cluster对结构进行团簇分析gmx cluster [-f [<.xtc/.trr/...>]] [-s [<.tpr/.tpb/...>]] [-n [<.ndx>]] [-dm [<.xpm>]] [-om [<.xpm>]] [-o [<.xpm>]] [-g [<.log>]] [-dist [<.xvg>]] [-ev [<.xvg>]] [-conv [<.xvg>]] [-sz [<.xvg>]] [-tr [<.xpm>]] [-ntr [<.xvg>]] [-clid [<.xvg>]] [-cl [<.xtc/.trr/...>]] [-nice ] [-b ] [-e ] [-dt ] [-tu ] [-[no]w] [-xvg ] [-[no]dista] [-nlevels ] [-cutoff ] [-[no]fit] [-max ] [-skip ] [-[no]av] [-wcl ] [-nst ] [-rmsmin ] [-method ] [-minstruct ] [-[no]binary] [-M ] [-P ] [-seed ] [-niter ] [-nrandom ] [-kT ] [-[no]pbc] gmx cluster可以使用几种不同的方法团簇化结构. 结构之间的距离可由轨迹来确定, 或使用-dm选项从.xpm矩阵文件读取. 结构间的距离可以由叠合后的RMS偏差或原子对距离的RMS偏差来定义. 确定团簇的方法有以下几种: single linkage(单连接): 当一个结构到团簇中任何一个原子的距离小于cutoff时, 就将此结构加入到团簇中. Jarvis Patrick: 当一个结构和团簇中的一个结构互为近邻结构, 并且至少有P个相同的邻近结构时, 将这个结构加入到团簇中. 一个结构的近邻结构是指距离它最近的M个结构或在cutoff以内的全部结构. Monte Carlo(蒙特卡洛): 利用蒙特卡洛方法重新排列RMSD矩阵, 以使帧的排列具有尽可能小的递增顺序. 这样做可以使从一个结构到另一个结构的动画尽量平滑, 并且结构之间的具有最大的可能(例如)RMSD, 但中间步骤应尽可能小. 这种方法可用于显示模拟的平均力势能系综或牵引模拟. 显然, 用户要仔细地准备轨迹文件(例如不能存在叠加的帧). 最终的结果可以通过.xpm矩阵文件进行直观的检查, 此文件从下到上都应该平滑地变化. diagonalization(对角化): 对角化RMSD矩阵. gromos: 利用Daura等介绍的算法(Angew. Chem. Int. Ed. 1999, 38, pp 236-240). 使用截断来数算近邻结构的个数, 把具有最多近邻的结构及其所有近邻作为一个团簇, 并从团簇池中将这个团簇移除. 然后对团簇池中剩下的结构重复以上算法. 当团簇化算法(single linkage, Jarvis Patrick and gromos) 将每个结构都精确地分配到了一个团簇, 并且提供了轨迹文件时, 在每一个团簇中, 相对于其他结构或平均结构或所有结构拥有最小平均距离的结构将被写入到轨迹文件中. 当输出所有结构时, 对每个团簇会使用单独编号的文件. 程序总会给出两个输出文件: -o: 输出矩阵左上半区域的RMSD值, 团簇图像的右下半区域. 当-minstruct = 1时, 若两个结构属于同一团簇, 相应的图像点为黑色; 当-minstruct > 1时, 对每一个团簇使用不同的颜色. -g: 输出所用选项的信息和所有团簇及其成员的详细列表. 此外, 程序也可以给出多个可选的输出文件: -dist: 输出RMSD的分布 -ev: 输出RMSD矩阵对角化的特征向量 -sz: 输出团簇的大小 -tr: 输出两个团簇之间的转变次数矩阵 -ntr: 输出从/到每个团簇的总转变次数 -clid: 输出团簇数随时间变化的函数 -cl: 输出每个团簇的平均(利用-av选项)或中心结构, 或对所选的一组团簇, 将团簇成员输出到带编号的文件(利用-wcl选项, 取决于-nst和-rmsmin选项). 一个团簇的中心是指团簇中与所有其他结构具有最小平均RMSD的结构.
|
|||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||