Re: [AMBER] branched aminoacids

From: FyD <>
Date: Thu, 17 Dec 2009 10:49:16 +0100


> I'd like to derive a parameter file for an unnatural aminoacid having
> one amino head and two carboxy tails, but I have a couple of questions.
> First of all, is this possible with amber? How do I define the two
> carboxy tails in the off file?

If you look at the cystine (not cysteine) residue: CYX in the Amber
Force Field Topology DataBase I think you get an example of what you

xleap -f leaprc.ff99SB

desc CYX
> desc CYX
UNIT name: CYX
Head atom: .R<CYX 1>.A<N 1>
Tail atom: .R<CYX 1>.A<C 9>
R<CYX 1>
> desc CYX.1
RESIDUE sequence number: 1
RESIDUE PDB sequence number: 0
Type: protein
Connection atoms:
  Connect atom 0: A<N 1>
  Connect atom 1: A<C 9>
  Connect atom 2: A<SG 8>

  => connect0, connect1 & connect2 are defined

> Is it possible to use a standard version
> of REDIII for deriving RESP charges?

With R.E.D.-III.3, you will be able to derive the N-terminal,
C-terminal and Central fragments in three different jobs. With R.E.D.
Server/R.E.D.-IV all can be done in a single step.

See for the
global strategy for the three different fragments.

- If we only discuss about the central fragment:
   I think you need a fragment like:

     R with R = R'CO

to get this fragment, you could use the following 'whole' molecule:
      R with R = R'CO-NME

This means you simply use 3 intra-molecular charge constraints
(INTRA-MCC keyword) in the corresponding P2N file (for ACE & 2 NMEs)
and the job is done...

You follow exactly the same strategy for the N-terminal & C-terminal
fragments using Methylammonium & Acetate - but here you will have to
replace one intra-molecular charge constraint by one inter-molecular
charge constraint (INTER-MCC keyword). Once again this is

regards, Francois

AMBER mailing list
Received on Thu Dec 17 2009 - 02:00:02 PST
Custom Search