Re: AMBER: peoblemwith RED II

From: FyD <>
Date: Fri, 27 Aug 2004 12:02:37 +0200


> Thanks for the reply. I was able to rectify the problem
> without modification to the RED source code. I have calculated the RESP
> charges and then created the frcmod and prepin files for the new residue.
> But now I face another problem.


> 1. I made the modifications to the protein structure by attaching an
> organic molecule to an aminoacid residue in INSIGHT II. In the saved pdb
> file the attached organic molecule is present in the end of the file as
> the last residue. Now I load my prepin and frcmod files followed by this
> molecule in leap and make proper bonds between the organic molecule and
> the aminoacid. Now when I check this structure it shows a few warnings of
> close contacts and an error:

If I understand you, you have two part in your residue... The aminoacids and the
org. fragment.
This means you have two residues and you have a branched molecule...

> Created a new atom named: OXT within residue: .R<ALA 582>

Did you check the atom names of this ALA 582? It must match the ones in the ALA
OFF library of xleap.

> How to overcome this problem?

Check carefully the atom name AND residue name in your PDB file and in your
xleap OFF libraries...
They have to be the same...

> 2. Next I tried to make modifications at 2 aminoacids. Now when I tried to
> load this structurte in leap, leap aborts with this error:
> ERROR: Comparing atoms C15, H161, H162, H163 to atoms C15, O11, H161, H162
> !FATAL ERROR----------------------------------------
> !FATAL: In file [chirality.c], line 120 !FATAL: Message: Atom O11 is not
> in the first list

First see

I wonder if your problem does not come from the fact that your have two residues
in this new aminoacid and consequently that you have a brached molecule.
I had this problem once and was able to solve it checking carefully the 'head'
and 'tail' of each residue...
See the xleap documentation for this.
Moreover, I wonder if you should not had a TER character between the two parts
of your new aminoacids in your PDB file as your molecule is branched...

Regards, Francois
The AMBER Mail Reflector
To post, send mail to
To unsubscribe, send "unsubscribe amber" to
Received on Thu Sep 02 2004 - 01:15:11 PDT
Custom Search