AMBER: build a polymer with sequence command

From: Shuting Wei <gtg549i.mail.gatech.edu>
Date: Fri, 22 Sep 2006 16:41:43 -0400

Dear Amber,

I have to send this again cause i didn't find a solution. Attached
the input file of PAM monomer (prepin, and frcmod)

I tried to build a polymer chain in leap with sequence command, but i
got the following error:

> loadamberprep polyamine.prepin
Loading Prep file: ./polyamine.prepin
> loadamberparams polyamine.frcmod
Loading parameters: ./polyamine.frcmod
Reading force field mod type file (frcmod)
> set PAM tail PAM.1.12
> desc PAM
UNIT name: PAM
Head atom: .R<PAM 1>.A<C1 1>
Tail atom: .R<PAM 1>.A<N4 12>
Contents:
R<PAM 1>
> triPAM = sequence { PAM PAM PAM }
ERROR: Comparing atoms N4, H14, H16, H18 to atoms C2, N4, H18, H14
!FATAL ERROR----------------------------------------
!FATAL: In file [chirality.c], line 121
!FATAL: Message: Atom C2 is not in the first list
!
!ABORTING.

What's wrong with my input? Thanks a lot!

Shuting


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

Received on Sun Sep 24 2006 - 06:07:20 PDT
Custom Search