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.scripps.edu
> To unsubscribe, send "unsubscribe amber" to majordomo.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.scripps.edu
To unsubscribe, send "unsubscribe amber" to majordomo.scripps.edu
Received on Wed Sep 06 2006 - 06:07:02 PDT
