Hello Users<br>
<br>
i am carrying out a dimeric protein simulation with G43a1 force field using Gromacs 3.3<br>
<br>
I&nbsp; have added some missing residues into my protein using a
seperate program. Now before adding water i want to do some simulation
in vaccum to minimized those added resideus and also the hydrogens
added via pdb2gmx command.<br>
<br>
For these i am freezing the protein and crystal waters leaving the
added residues so that they can freely move. for that i am using the
following input (.mdp)&nbsp; file:<br>
<br>
<br>
&nbsp;title&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; = minimizatio in vaccum<br>
integrator&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; = steep<br>
dt&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; = 0.002<br>
constraints&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; = none<br>
nsteps&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; = 500<br>
nstlist&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; = 10<br>
rcoulomb&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; = 1<br>
rvdw&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; = 1<br>
ns_type&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; = simple<br>
coulombtype&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; = PME<br>
vdwtype&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; = cut-off<br>
pbc&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; = no<br>
gen_vel&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; = no<br>
freezegrps&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; = protein&nbsp; solvent<br>
freezedim&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; = Y Y Y&nbsp;&nbsp;&nbsp;&nbsp; Y Y Y<br>
<br>
emtol&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; = 100.0<br>
emstep&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; = 0.01<br>
<br>
With that&nbsp; i have made an index file where i have mentioned two groups with names &#39;protein&#39; and &#39;solvent&#39;  <br>
<br>
but inspite of this i am getting the following problem :<br>
....................................................<br>
<br>
Fatal error:<br>
Group protein not found in indexfile.<br>
Maybe you have non-default goups in your mdp file, while not using the &#39;-n&#39; option of grompp.<br>
In that case use the &#39;-n&#39; option.<br>
................................................<br>
<br>
although i have mentioned the index file in the grompp :<br>
<br>
&quot;aquapo.ndx is the name of&nbsp; index file &quot;<br>
grompp -f 1min.mdp -c aquapo_grom.pdb -p aquapo_grom.top -n aquapo.ndx -o 1min.tpr <br>
<br>
<br>
index file used <br>
&nbsp;[protein]<br>
&nbsp; 1 2 3 4 5 6 7 8 9 11 12 13 14 15 16 17 18 19 20 21 22 23 24 26
27 28 29 30 31 32 33 34 35 37 38 39 40 41 42 43 44 45 47 48 49 50 51 52
53 55 56 57 58 59 60&nbsp; 61 64 65 66 67 68 69 70 72 73 74 75 76 77 78
79 80 81 83 85 86 87 88 89 91 92 93 94 95 96 97 98 99 100 101 102 103
104 105 106 107 108 110 111 112 113 114&nbsp;&nbsp; 115 116 117 118 119
120 121 122 123 125 126 127 128 129 130 131 132 133 134 136 137 138 139
140 141 142 143 144 146 147 148 149 150 151 152 154 155 156 157&nbsp;
158 159 160 163 164 165 166 167 168 169 171 172 173 174 175 176 177 178
179 180 182 184 185 186 187 188 190 191 192 193 194 195 196 197 198<br>
<br>
<br>
&nbsp; [solvent]<br>
&nbsp; 199 200 201 202 203 204 205 206 207 208 209 210 211 212 213 214
215 216 217 218 219 220 221 222 223 224 225 226 227 228 229 230 231 232
234 235 236 237 238&nbsp; 239 240 241 242&nbsp; 243 244 245 246 247 248
249 250 251 252 253 254 255 256 257 258 259 260 261 262 263 264 265 266
267 268 269 270 271 272 273 274 275<br>
<br>
<br>
<br>
please tell what is the reason for this error <br>
<br>
waiting for your reply <br>
With Regards<br>
Gurpreet<br>
&nbsp;<br>
<br>
<br>
<br>
<br>