Re: AMBER: not equivalent oxygens in COO- / how to force order of atoms in IMPROPER?

From: Francesco Pietra <chiendarret.yahoo.com>
Date: Wed, 19 Mar 2008 10:26:24 -0700 (PDT)

Hi Edyta:
Please, put me in the list of those very much interested in your Python script.
Sincerely

francesco pietra

--- Edyta Malolepsza <em427.cam.ac.uk> wrote:

> Dear David,
>
> > > The reason can be in order of atoms in IMPROPER angles for NH2 and COO-
> groups.
> >
> > This will only really be fixed by a change in the way the impropers are
> > defined, I think. Someone will have to take this on as a project; it will
> > be easier once we implement the CHARMM force field equations (but don't
> hold
> > your breath.)
>
> It can be done also by reordering atom positions in topology file
> without any code or energy formula change. I have written program in
> python which reads original topology file and AMBER libraries for
> amino acids and nucleic to create new topology file with correct order
> of atoms for IMPROPER's. With new topology file any permutations of
> equivalent atoms give the same results - potential is symmetric.
>
> I can share this program, please let me know if somebody is interested.
>
> With best regards,
>
> Edyta Malolepsza
> -----------------------------------------------------------------------
> The AMBER Mail Reflector
> To post, send mail to amber.scripps.edu
> To unsubscribe, send "unsubscribe amber" to majordomo.scripps.edu
>



      ____________________________________________________________________________________
Looking for last minute shopping deals?
Find them fast with Yahoo! Search. http://tools.search.yahoo.com/newsearch/category.php?category=shopping
-----------------------------------------------------------------------
The AMBER Mail Reflector
To post, send mail to amber.scripps.edu
To unsubscribe, send "unsubscribe amber" to majordomo.scripps.edu
Received on Fri Apr 18 2008 - 21:12:19 PDT
Custom Search