[gmx-users] problem with energy minimization

[Geraudis Mustelier]

Hi

I would like to refine a model obtained with modeling by homology technique
from X-ray structures. So, I am treating to do an energy minimization. The
parameter values in .mdp file are:
;
cpp                 =  /lib/cpp
define              = -DPOSRES
constraints       =  none
morse              =  no
integrator         =  steep
nsteps              =  2000
;
; Energy minimizing stuff
emtol               =  100
emstep             =  0.00001
;
nstcomm           =  1
ns_type            =  grid
rlist                   =  1
rcoulomb          =  1.0
rvdw                =  1.0
Tcoupl              =  no
Pcoupl              =  no
gen_vel             =  no
 
 But the last results showed in log file are:

Reading file r3md01.tpr, VERSION 3.0.5 (single precision)

Back Off! I just backed up ener.edr to #ener.edr#
STEEPEST DESCENTS
   Tolerance         =          100

Stepsize too small Converged to machine precision,
but not to the requested precision (100)
You might need to increase your constraint accuracy

writing lowest energy coordinates.

Back Off! I just backed up r3md01.trr to #r3md01.trr#

Back Off! I just backed up r3md02.pdb to #r3md02.pdb#
Maximum force:  2.14823e+03

STEEPEST DESCENTS did not converge in 82 steps
  Potential Energy  = -2.11852e+05

I have explored with different values for these parameters, but always the
final structure is identical to initial. So, someone can help me to give
correct value to these parameters, or introduce others, in order to obtain
another conformation?
 
 ====================================================================
 MSc. Geraudis Mustelier-Portuondo

 Center of Molecular Immunology
 Calle 216 esq. A 15, Atabey, Playa
 A. Postal 16040, Havana 11600, CUBA

 email: gera at ict.cim.sld.cu
 http://bioinfo.cim.sld.cu

 Tel: 53-7-214335
 Fax: 53-7-335049
 ==================================================================