AMBER Archive (2006)

Subject: Re: AMBER: Question about LEAP and twisting a dihedral angle

From: Ilyas Yildirim (yildirim_at_pas.rochester.edu)
Date: Sun Sep 03 2006 - 01:20:27 CDT


Dear Dave,

I have attached the .pdb file I am using. It is a tetramer, (ggcc.pdb).
Here is the xleap.in I am using:

------------ xleap.in ------------------
source leaprc.rna.ff99
ggcc = loadpdb ggcc.pdb
impose ggcc { 3 4 } {
  { C3' O3' P O5' 0 }
}
savepdb ggcc ggcc_new.pdb
quit
----------------------------------------

The tetramer has 2 strands and I am changing just one dihedral angle of
strand # 1 to zero. The new .pdb file (ggcc_new.pdb) is totally different.
Yes, it keeps this particular dihedral angle at zero degrees; but it
changes all of the atom positions of strand #1 while keeping the atom
positions of strand # 2 untouched. I know that what I am doing is not
right; I need to describe some more restraints to make leap understand
that strand #2 is part of the structure, and that it has to be helical,
but I am not sure how to do it.

Right now, I am writing a script which is going to rotate a structure in a
particular orientation, and plan on doing some complicated procedures
to just change a dihedral angle. Its pretty time consuming and
therefore I am wondering if I can solve this problem with leap, using the
'impose' command. Thanks in advance.

Best,

On Sat, 2 Sep 2006, David A. Case wrote:

> On Sat, Sep 02, 2006, Ilyas Yildirim wrote:
>
> > The impose command does not work; it somehow
> > changes the external coordinates of all of the atoms.
>
> It would be helpful to post details of the problem. I'm not sure what you
> mean by "external coordinates", but using the impose command is tricky, but
> people in Rob Woods' group have figured a lot out about this. Maybe someone
> could help if you gave details.
>
> I don't think "twist" works from the command line.
>
> ...dac
>
> -----------------------------------------------------------------------
> The AMBER Mail Reflector
> To post, send mail to amber_at_scripps.edu
> To unsubscribe, send "unsubscribe amber" to majordomo_at_scripps.edu
>
>

-- 
  Ilyas Yildirim
  ---------------------------------------------------------------
  - Department of Chemisty       -				-
  - University of Rochester      -				-
  - Hutchison Hall, # B10        -				-
  - Rochester, NY 14627-0216     - Ph.:(585) 275 67 66 (Office)	-
  - http://www.pas.rochester.edu/~yildirim/			-
  ---------------------------------------------------------------



----------------------------------------------------------------------- The AMBER Mail Reflector To post, send mail to amber_at_scripps.edu To unsubscribe, send "unsubscribe amber" to majordomo_at_scripps.edu