SHAKE problem

From: Rosalia Pascual <>
Date: Wed 07 Aug 2002 17:13:12 +0200


I have a problem with SHAKE ( coordinate resetting cannot be ...)
My system is a protein-ligand solvated with a cap of waters of 20
(total charge of the
system is 6).

First, I have minimized it in the following way ( AMBER 7):

Minimization with Cartesian restraints
    imin=1, maxcyc=5000,
    scee=1.2, igb=0, ntb=0,
    ntr=0, ibelly=0,
    ntmin=1, ncyc=1000, nsnb=25, cut=10,

and then I have run a dynamics for 100 ps:

    ntx=1, ntpr=100, ntwx=100, ntwv=100, ntwe=100,

    ntf=2, ntc=2,
    ntb=0, igb=0, cut=10, scnb=2, scee=1.2,
    ibelly=0, ntr=1,

    nstlim=50000, ntt=1, tautp=0.2,
   type='TEMP0', istep1=0, istep2=10000,
                 value1=0, value2=300.0,
   type='TEMP0', istep1=10000, istep2=50000,
                 value1=300.0, value2=300.0,
fixed residues (away from the active center)
RES 1 2
RES 4 6
RES 65 65
RES 68 86
RES 116 139
RES 172 172
RES 174 187

After 100 ps, the system is equilibrated ( ETOT constant, RMS
Then, I try to run a production ( same as above) for 300 ps.
However I have a SHAKE problem.

I have seen different solutions in similar e-mails, but it doesn`t help

1) VLIMIT is 20

2) The cutoff is smaller ( 10 ) than the diameter of the water cap (20 )

happens if I use cutoff 9 or 8 ).

3) The problem happens 104 ps after the dynamics ( so, maybe is not a
from the initial structure - I have checked it and it seems OK).

4) If I turn off the SHAKE OPTION ( although I have a cap of
WATBOX216 ), temperature raises a lot.

Could anybody help me?
Thanks in advance.
Received on Wed Aug 07 2002 - 08:13:12 PDT
Custom Search