Re: [AMBER] Analysis of minimization stage

From: Valentina Romano <valentina.romano.unibas.ch>
Date: Mon, 12 May 2014 08:47:42 +0000

Dear Amber users

I am still having problems with the equilibration stage.

My system is a complex of a protein kinase binds to the 6amino-purine.
The whole system contains 247 residues.

First of all, I minimized the solvent (water +ions) keeping the protein's heavy atoms frozen:

 &cntrl
  imin=1,
  maxcyc=1000,
  ntb=1,
  ntr=1,
  cut=10,
  restraint_wt = 10.0,
  restraintmask=':1-247.CA,C,O,N,H'
 /

Afterwards, I wanted to run a short MD keeping frozen the protein's heavy atoms to relax solvent molecules before running a minimization of the whole system:

 &cntrl
  imin=0,
  irest=0,
  ntx=1,
  ig=-1,
  ntb=1,
  ntr=1,
  cut=10,
  ntc=2,
  ntf=2,
  tempi=300.0,
  temp0=300.0,
  ntt=3,
  gamma_ln=5.0,
  nstlim=15000, dt=0.001,
  ntpr=100, ntwx=100, ntwr=1000
  restraint_wt = 10.0,
  restraintmask=':1-247.CA,C,O,N'

Sander immediately stops with the following error message:

vlimit exceeded for step 2; vmax = 144.9409

     Coordinate resetting (SHAKE) cannot be accomplished,
     deviation is too large
     NITER, NIT, LL, I and J are : 0 2 1927 3834 3844

     Note: This is usually a symptom of some deeper
     problem with the energetics of the system.

If I decrease the dt from 0.001 to 0.0005, sander runs the MD but between each step i got the following messages:

e.g.:
vlimit exceeded for step 4001; vmax = 22.8241
vlimit exceeded for step 4003; vmax = 25.4000
vlimit exceeded for step 4004; vmax = 27.7180
vlimit exceeded for step 4006; vmax = 29.4842
vlimit exceeded for step 4007; vmax = 28.0214
vlimit exceeded for step 4008; vmax = 56.6682
vlimit exceeded for step 4009; vmax = 45.4394
vlimit exceeded for step 4010; vmax = 31.1217
vlimit exceeded for step 4011; vmax = 24.5171
vlimit exceeded for step 4012; vmax = 34.7059

 Moreover the TEMP(K) is not around 300.0 at each step but it changes a lot.

All this things suggested me that there is a problem in the structure I am using to run the MD.
After the mininimization stage I checked the TOTAL energy and it was relatively flat at the end of the mninimization. Thus I have no idea how to find out the problem with the minimized structure i am using.

Any suggestions?
I would like the use cheackoverlap in cptraj but is not really clear how it works.

Best,
Valentina


~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
Valentina Romano | PhD Student | Biozentrum, University of Basel & SIB Swiss Institute of Bioinformatics
Klingelbergstrasse 61 | CH-4056 Basel |

Phone: +41 61 267 15 80






_______________________________________________
AMBER mailing list
AMBER.ambermd.org
http://lists.ambermd.org/mailman/listinfo/amber
Received on Mon May 12 2014 - 02:00:03 PDT
Custom Search