<br>Hi!<br><br>I&#39;m am trying to visualize the energy minimization trajectory path of my protein+ligand in VMD. I imported the .trr file into the ending energy minimized pdb file, and I was able to see the frames of the trr (53 all together with nstxout=1). But when I ran the trajectory I saw barely any change in the position on my protein. It just shifted slightly before the last frame. The rest of the time, the protein was stationary. My protein is 243 residues, so not very large, and it includes the ligand (small as well) and the protein; I have constrained nothing. I&#39;m using Gromacs 4.0.4.<br>
<br>When I did the energy minimization, the em.log output file did say that the protein converged to machine precision but not to my specified precision (Fmax&lt;100), and that my stepsize might be too small or the energy change not significant. Is this message a problem? Is it because maybe my energy change isn&#39;t significant, that my protein is not fluctuating much? The em.log output file showed that it took 79 steps, steepest descent, to minimize my protein..the beginning potential was -1.5E5 kj/mol and the ending potential was -4.9E5....I graphed the energies vs step and it did show a significant, stable convergence of potential.<br>
 <br>Anyway, my em.mdp file is below. If I can get some suggestions on how to more easily see my protein fluctuating to its minimized state, in vmd, that would be great! Should I maybe use the .xtc trajectory file instead of the .trr? Or should I add the lines nstvout or nstfout?<br>
<br>title                           =  BR6_trp<br>cpp                             =  /usr/bin/cpp<br>define                          =  -DFLEXIBLE<br>constraints                     =  none<br>integrator                      =  steep<br>
dt                              =  0.002 ; ps !<br>nsteps                          =  500<br>nstxout                         =  1<br>nstlist                         =  10<br>ns_type                         =  grid<br>rlist                           =  0.9<br>
coulombtype                     =  PME<br>rcoulomb                        =  0.9<br>vdwtype                         =  cut-off<br>rvdw                            =  0.9<br>fourierspacing                  =  0.12<br>fourier_nx                    =  0<br>
fourier_ny                    =  0<br>fourier_nz                    =  0<br>pme_order                     =  4<br>ewald_rtol                    =  1e-5<br>optimize_fft                  =  yes<br><br><br>Thanks in advance,<br>
Halie Shah<br>University of Houston<br>Computational Biology Lab<br>