1g_rms(1) GROMACS suite, VERSION 4.5 g_rms(1)
2
3
4
6 g_rms - calculates rmsd's with a reference structure and rmsd matrices
7
8 VERSION 4.5
9
11 g_rms -s topol.tpr -f traj.xtc -f2 traj.xtc -n index.ndx -o rmsd.xvg
12 -mir rmsdmir.xvg -a avgrp.xvg -dist rmsd-dist.xvg -m rmsd.xpm -bin
13 rmsd.dat -bm bond.xpm -[no]h -[no]version -nice int -b time -e time -dt
14 time -tu enum -[no]w -xvg enum -what enum -[no]pbc -fit enum -prev int
15 -[no]split -skip int -skip2 int -max real -min real -bmax real -bmin
16 real -[no]mw -nlevels int -ng int
17
19 g_rms compares two structures by computing the root mean square devia‐
20 tion (RMSD), the size-independent 'rho' similarity parameter (rho) or
21 the scaled rho (rhosc), see Maiorov & Crippen, PROTEINS 22, 273
22 (1995). This is selected by -what.
23
24 Each structure from a trajectory ( -f) is compared to a reference
25 structure. The reference structure is taken from the structure file (
26 -s).
27
28
29 With option -mir also a comparison with the mirror image of the refer‐
30 ence structure is calculated. This is useful as a reference for 'sig‐
31 nificant' values, see Maiorov & Crippen, PROTEINS 22, 273 (1995).
32
33
34 Option -prev produces the comparison with a previous frame the speci‐
35 fied number of frames ago.
36
37
38 Option -m produces a matrix in .xpm format of comparison values of
39 each structure in the trajectory with respect to each other structure.
40 This file can be visualized with for instance xv and can be converted
41 to postscript with xpm2ps.
42
43
44 Option -fit controls the least-squares fitting of the structures on
45 top of each other: complete fit (rotation and translation), translation
46 only, or no fitting at all.
47
48
49 Option -mw controls whether mass weighting is done or not. If you
50 select the option (default) and supply a valid tpr file masses will be
51 taken from there, otherwise the masses will be deduced from the atom‐
52 mass.dat file in the GROMACS library directory. This is fine for pro‐
53 teins but not necessarily for other molecules. A default mass of 12.011
54 amu (Carbon) is assigned to unknown atoms. You can check whether this
55 happend by turning on the -debug flag and inspecting the log file.
56
57
58 With -f2, the 'other structures' are taken from a second trajectory,
59 this generates a comparison matrix of one trajectory versus the other.
60
61
62 Option -bin does a binary dump of the comparison matrix.
63
64
65 Option -bm produces a matrix of average bond angle deviations analo‐
66 gously to the -m option. Only bonds between atoms in the comparison
67 group are considered.
68
70 -s topol.tpr Input
71 Structure+mass(db): tpr tpb tpa gro g96 pdb
72
73 -f traj.xtc Input
74 Trajectory: xtc trr trj gro g96 pdb cpt
75
76 -f2 traj.xtc Input, Opt.
77 Trajectory: xtc trr trj gro g96 pdb cpt
78
79 -n index.ndx Input, Opt.
80 Index file
81
82 -o rmsd.xvg Output
83 xvgr/xmgr file
84
85 -mir rmsdmir.xvg Output, Opt.
86 xvgr/xmgr file
87
88 -a avgrp.xvg Output, Opt.
89 xvgr/xmgr file
90
91 -dist rmsd-dist.xvg Output, Opt.
92 xvgr/xmgr file
93
94 -m rmsd.xpm Output, Opt.
95 X PixMap compatible matrix file
96
97 -bin rmsd.dat Output, Opt.
98 Generic data file
99
100 -bm bond.xpm Output, Opt.
101 X PixMap compatible matrix file
102
103
105 -[no]hno
106 Print help info and quit
107
108 -[no]versionno
109 Print version info and quit
110
111 -nice int 19
112 Set the nicelevel
113
114 -b time 0
115 First frame (ps) to read from trajectory
116
117 -e time 0
118 Last frame (ps) to read from trajectory
119
120 -dt time 0
121 Only use frame when t MOD dt = first time (ps)
122
123 -tu enum ps
124 Time unit: fs, ps, ns, us, ms or s
125
126 -[no]wno
127 View output xvg, xpm, eps and pdb files
128
129 -xvg enum xmgrace
130 xvg plot formatting: xmgrace, xmgr or none
131
132 -what enum rmsd
133 Structural difference measure: rmsd, rho or rhosc
134
135 -[no]pbcyes
136 PBC check
137
138 -fit enum rot+trans
139 Fit to reference structure: rot+trans, translation or none
140
141 -prev int 0
142 Compare with previous frame
143
144 -[no]splitno
145 Split graph where time is zero
146
147 -skip int 1
148 Only write every nr-th frame to matrix
149
150 -skip2 int 1
151 Only write every nr-th frame to matrix
152
153 -max real -1
154 Maximum level in comparison matrix
155
156 -min real -1
157 Minimum level in comparison matrix
158
159 -bmax real -1
160 Maximum level in bond angle matrix
161
162 -bmin real -1
163 Minimum level in bond angle matrix
164
165 -[no]mwyes
166 Use mass weighting for superposition
167
168 -nlevels int 80
169 Number of levels in the matrices
170
171 -ng int 1
172 Number of groups to compute RMS between
173
174
176 gromacs(7)
177
178 More information about GROMACS is available at <http://www.gro‐
179 macs.org/>.
180
181
182
183 Thu 26 Aug 2010 g_rms(1)