|
|||||||||||||||||||||||||||||||||
AMBER Archive (2003)Subject: AMBER: Restrain of dihedral angle!
From: Marco Aurelio Correia Preto (mcpreto_at_fc.up.pt)
Dear Amber users:
We've been testing some parameters for non-coded amino acids!
I know Amber does not support any option to do this automatically, but I was hopping we could restrict the dihedral to a fixed value, perform an optimization on the hole unit and then (after frreiing the restriction) perform a single point to have an MM energy read!
It would be necessary to do this to each point of the profile!
We've been trying to use restraintorsion unit a1 a2 a3 a4 999999999 danglevalue 1 in xleap to keep the dihedral on the disered value but we are not achiving this objective.
Unless we are wrongly difining the atoms (we use unit.residuenumber.atomname) we don't understand the obtained values! We would expect a constrain en value of zero for the first point of the optimization (because the starting geometry has the desired value), but we don't.
We belive that the restrain is not being correctly aplyed by us.
Sincerely yours
Marco Preto
-----------------------------------------------------------------------
| |||||||||||||||||||||||||||||||||
|