On Mon, Aug 04, 2008, Alex Spasic wrote:
>
> 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,
> com(:2.C1',:2.C2',:2.C3',:2.C4',:2.O4'),:3.P,com(:3.C1',:3.C2',:3.C3',:3.C4',:3.O4')]",
> r1 = -180, r2 = -150.0, r3 = -150.0, r4 = -120.0,
> rk2 = 30.0, rk3 = 30.0
> /
>
As I remember (and from browsing briefly through he code) only distance
restraints can use the center of mass idea. It looks like the new format for
restraint input doesn't enforce this restriction. Matthew Seetin should look
at this and confirm (or not), but my intial guess is that you can't do what
you are trying here.
...dac
-----------------------------------------------------------------------
The AMBER Mail Reflector
To post, send mail to amber.scripps.edu
To unsubscribe, send "unsubscribe amber" (in the *body* of the email)
to majordomo.scripps.edu
Received on Wed Aug 06 2008 - 06:07:22 PDT
