AMBER: Problem with setting the NMR restraints

From: Alex Spasic <>
Date: Mon, 04 Aug 2008 10:15:25 -0400

Hello everyone,

I am trying to run a simulation of a small single chain DNA fragment
using Amber 10 while holding fixed the torsion angle formed by the:

(phosphate)i-(center of mass of sugar)i-(phosphate)i+1-(center of mass
of sugar)i+1

The restraint is defined by the following lines:

# psps restraint
 &rst restraint = "torsion[:2.P,
          r1 = -180, r2 = -150.0, r3 = -150.0, r4 = -120.0,
          rk2 = 30.0, rk3 = 30.0

However this restraint doesn't produce the angles centered around -150.
If I use just one of the of the atoms from the sugar (instead of their
center of mass) the results are fine. My question is whether I have made
a mistake in setting up these restrains or if this is a bug in the code?
Has anyone else had similar problems?

Thanks in advance,
The AMBER Mail Reflector
To post, send mail to
To unsubscribe, send "unsubscribe amber" (in the *body* of the email)
Received on Wed Aug 06 2008 - 06:07:21 PDT
Custom Search