Re: AMBER: MD problems/shake-many thanks

From: <>
Date: Wed, 21 Dec 2005 00:04:24 +0100 (CET)

Dear Professor Cheatham!

Thank you for helping me.
The calculations are good now, and everything is right.
Many many thanks once again!

Pawel Gruszczynski
University of Gdansk
Chemistry Department

> Yes, this is a ligand (ATP) problem. The atoms HO-OH in the first
> phosphoreus group. The parameters to ATP are from parm94.dat.
> I have tried to do dynamics on isolated ligand having shake on,
> and that is a problem. I got the SHAKE error again on the same
> atoms:
> "Coordinate resetting (SHAKE) cannot be accomplished,
> deviation is too large
> NITER, NIT, LL, I and J are : 0 0 1 1 2

As a followup to the list regarding what was happening and what to
for... If you have an isolated ligand and SHAKE failures are
evident in
MD simulation, something is likely funny with the parameters that
atoms to move within vdw spheres leading to large forces such that
integrator or SHAKE fails. Working with Pawel, I obtained the
prmtop/inpcrd and did some experimenting.

(1) If you get SHAKE failures on your drug/ligand, try first
minimization of the isolated compound in-vacuo. If atoms are
vdw spheres (leading to overlap), this will be readily apparent
looking at the structure. In this case, the HO atom (with no vdw
parameters) overlapped with a nearby oxygen. This is occurring
since the
HO atom has a relatively large (~ +0.4) charge and strongly
favorably with the nearby oxygens (~ -0.6 charge). As there is no
repulsion, this attraction overcomes the bond and angle forces
it away.

(2) Running MD simulation on the isolated system in-vacuo shows
temperature cannot be maintained (it quickly rises) which is
evidence of the instability.

I got around this problem (assuming the charges were reasonable) by
adding a vdw potential to this hydrogen. I changed the OH
(for this compound) to have a small vdw like HS (0.6 0.0157) and
minimization and MD were stable (although the hydrogen in question
still very strongly attracted to the nearby oxygen, it does not
the vdw sphere).

The AMBER Mail Reflector
To post, send mail to
To unsubscribe, send "unsubscribe amber" to
Received on Tue Dec 20 2005 - 23:53:00 PST
Custom Search