AMBER Archive (2005)

Subject: RE: AMBER: xleap with modified residue

From: Ross Walker (ross_at_rosswalker.co.uk)
Date: Thu Oct 20 2005 - 13:58:53 CDT


Dear Whitney,
 
This sounds like it could be a bug with leap. Can you try repeating all the
steps but this time call the unit MGG instead of 2MG. The number at the
beginning of the unit name may be giving problems. Also replace the unit
name mol with test. So
 
test = loadpdb test.pdb
 
Please make sure you have applied all the bugfixes from
http://amber.scripps.edu/bugfixes80.html first since several of these deal
with leap crashes.
 
If you still have problems please email me back with the library and pdb
files and I'll see if I can reproduce the error.
 
All the best
Ross
 

/\
\/
|\oss Walker

| Department of Molecular Biology TPC15 |
| The Scripps Research Institute |
| Tel: +1 858 784 8889 | EMail:- ross_at_rosswalker.co.uk |
| http://www.rosswalker.co.uk <http://www.rosswalker.co.uk/> | PGP Key
available on request |

Note: Electronic Mail is not secure, has no guarantee of delivery, may not
be read every day, and should not be used for urgent or sensitive issues.

 

  _____

From: owner-amber_at_scripps.edu [mailto:owner-amber_at_scripps.edu] On Behalf Of
Whitney Allen
Sent: Thursday, October 20, 2005 11:13
To: amber_at_scripps.edu
Subject: AMBER: xleap with modified residue

Hi Amber users

I am trying to run MD on a TRNA structure (6tna.pdb). This molecule has
several modified bases and so I have used Antechamber to compute a molecule
library for each modified base. However, after loading the library file
into xleap (which is successful) I try and load the TRNA structure and xleap
completely shuts down.

Below is the procedure that I follow. Note that I have tried to follow
tutorials 4 & 5 as applicable. Here, for simplicity, I have cut the TRNA
PDB file so that it contains only one modified residue (2MG or residue 10 in
the 6tna.pdb file) and 15 total residues (residues 1-15 in the 6tna.pdb
file) and have called it test.pdb. If anyone can help please reply and I can
send you the 2MG.pdb, 2MG.mol2, 2mg.lib, and test.pdb files.

 

(Antechamber- note that I change the residue to a molecule)

 

% $AMBERHOME/exe/antechamber -i 2MG.pdb -fi pdb -o 2MG.mol2 -fo mol2 -c bcc
-s 2 -mk "CARTESIAN AM1 STANDARD DIRECT OPT=LBFGS XTEST=0.0001"

 

% $AMBERHOME/exe/parmchk -i 2MG.mol2 -f mol2 -o 2MG.frcmod

 

% xleap -f leaprc.rna.ff99

source leaprc.gaff

loadmol2 2MG.mol2

check 2MG

edit 2MG

loadamberparms 2MG.frcmod

check 2MG

saveoff 2MG.lib

desc 2MG

 

 

% xleap -f leaprc.rna.ff99

source leaprc.gaff

loadoff 2mg.lib

check 2MG

mol=loadpdb test.pdb

 

(At this point xleap closes and I get the following error)

 

!FATAL ERROR----------------------------------------

!FATAL: In file [chirality.c], line 120

!FATAL: Message: Atom P is not in the first list

!

!ABORTING.

 

Please help me understand what is going on here and let me know how I can
get this to work. It looks like all of the files are of normal format so I
am a bit confused.

Thanks in advance

Whitney Hastings

  _____

Yahoo!
<http://pa.yahoo.com/*http://us.rd.yahoo.com/evt=36035/*http://music.yahoo.c
om/unlimited/> Music Unlimited - Access over 1 million songs. Try it free.

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