AMBER: tleap load pdb file

From: xiaowen <>
Date: Fri, 4 Feb 2005 15:18:09 -0600 (CST)

Thank David replying my question.

Based on David's suggestions, here is the new results and some explanations.
Briefly, I can say the problem is "leap cannot load pdb file".

1. double-check that the test.leap test suite works on your IBM computer;

cd $AMBERHOME/test/leap> ./Run.tleap
(cannot pass, some error with the input file or in sander)


Running leap tests on ./tleap -> ../../exe/teLeap

    building force field libraries:

        '91 ff.. "(all passed)"

    loadpdb and saveamberparm: 94 ff..

       - aminoan94
       comparing topology to earlier leap version..
diffing ff94/ with
possible FAILURE: check
       comparing sander energies
1525-001 The READ statement on the file all_aminoan94.crd cannot be completed
because the end of the file was reached. The program will stop.
Run.tleap: Program error in sander (could be leap input problem)


2. try to make a smaller version of nuc.pdb, to see if that will work.
I modify to:
A5 A A A A3
T5 T T T T3
$AMBERHOME/exe> nucgen -O -i -o nuc2.out -d
$AMBERHOME/dat/leap/parm/nucgen.dat -p nuc2.pdb

so far is fine, the nuc2.pdb can be open by internet Explorer.
But go to
test=loadpdb nuc2.pdb
the same result as nuc.pdb's

3 I try to use the exist unit to do next steps.
under tleap: ( try small unit ACE)
>saveamberparm ACE ACE.prmtop ACE.inpcrd (works well)
           (if under Xleap, Xwin still cannot show the graphy on unit editor
>test1=copy ACE
>test2=copy ACE
>solvatebox test1 TIP3PBOX 8.0 (WORKS WELL)
>solvateoct test2 TIP3PBOX8.0 (WOKRS WELL)

> Is there actually a chainid "T" in the input? How many atoms does it have?

Yes, the file: total number of atoms is 438 from the nuc.pdb file.

A5 A A A A A A A A A3
T5 T T T T T T T T T3


Your suggestions and ideas are really helpful to a new UNIX and amber user.
Best regards!

The AMBER Mail Reflector
To post, send mail to
To unsubscribe, send "unsubscribe amber" to
Received on Fri Feb 04 2005 - 21:53:00 PST
Custom Search