1g_helixorient(1) GROMACS suite, VERSION 4.5 g_helixorient(1)
2
3
4
6 g_helixorient - calculates local pitch/bending/rotation/orientation
7 inside helices
8
9 VERSION 4.5
10
12 g_helixorient -s topol.tpr -f traj.xtc -n index.ndx -oaxis helix‐
13 axis.dat -ocenter center.dat -orise rise.xvg -oradius radius.xvg
14 -otwist twist.xvg -obending bending.xvg -otilt tilt.xvg -orot rota‐
15 tion.xvg -[no]h -[no]version -nice int -b time -e time -dt time -xvg
16 enum -[no]sidechain -[no]incremental
17
19 g_helixorient calculates the coordinates and direction of the average
20 axis inside an alpha helix, and the direction/vectors of both the alpha
21 carbon and (optionally) a sidechain atom relative to the axis.
22
23
24 As input, you need to specify an index group with alpha carbon atoms
25 corresponding to an alpha helix of continuous residues. Sidechain
26 directions require a second index group of the same size, containing
27 the heavy atom in each residue that should represent the sidechain.
28
29
30 Note that this program does not do any fitting of structures.
31
32
33 We need four Calpha coordinates to define the local direction of the
34 helix axis.
35
36
37 The tilt/rotation is calculated from Euler rotations, where we define
38 the helix axis as the local X axis, the residues/CA-vector as Y, and
39 the Z axis from their cross product. We use the Euler Y-Z-X rotation,
40 meaning we first tilt the helix axis (1) around and (2) orthogonal to
41 the residues vector, and finally apply the (3) rotation around it. For
42 debugging or other purposes, we also write out the actual Euler rota‐
43 tion angles as theta1-3.xvg
44
46 -s topol.tpr Input
47 Run input file: tpr tpb tpa
48
49 -f traj.xtc Input
50 Trajectory: xtc trr trj gro g96 pdb cpt
51
52 -n index.ndx Input, Opt.
53 Index file
54
55 -oaxis helixaxis.dat Output
56 Generic data file
57
58 -ocenter center.dat Output
59 Generic data file
60
61 -orise rise.xvg Output
62 xvgr/xmgr file
63
64 -oradius radius.xvg Output
65 xvgr/xmgr file
66
67 -otwist twist.xvg Output
68 xvgr/xmgr file
69
70 -obending bending.xvg Output
71 xvgr/xmgr file
72
73 -otilt tilt.xvg Output
74 xvgr/xmgr file
75
76 -orot rotation.xvg Output
77 xvgr/xmgr file
78
79
81 -[no]hno
82 Print help info and quit
83
84 -[no]versionno
85 Print version info and quit
86
87 -nice int 19
88 Set the nicelevel
89
90 -b time 0
91 First frame (ps) to read from trajectory
92
93 -e time 0
94 Last frame (ps) to read from trajectory
95
96 -dt time 0
97 Only use frame when t MOD dt = first time (ps)
98
99 -xvg enum xmgrace
100 xvg plot formatting: xmgrace, xmgr or none
101
102 -[no]sidechainno
103 Calculate sidechain directions relative to helix axis too.
104
105 -[no]incrementalno
106 Calculate incremental rather than total rotation/tilt.
107
108
110 gromacs(7)
111
112 More information about GROMACS is available at <http://www.gro‐
113 macs.org/>.
114
115
116
117 Thu 26 Aug 2010 g_helixorient(1)