AMBER: coordinate change of pdb structure using xleap

From: Henrik Öberg <>
Date: Wed, 27 Aug 2008 08:26:42 +0200

Dear all,
I am using AMBER 9 simulating an enzyme docking during a 1 ns period. My problem is not the actual simulation but the initial settings leading up to it. When I load my ligand pdb file into leap, after using antechamber and parmchk, through the following order:
source leaprc.gaff
loadamberprep ligand.prepin
loadamberparams ligand.frcmod
loadamberparams frcmod.acn
DOCK=loadpdb dock.pdb
check DOCK (here I recieve some warnings but overall the unit is ok)
addions DOCK Na+ 0
solvateoct DOCK2 ACNBOX 8.0
saveamberparm dock.prmtop dock.inpcrd

It is now I observe the problem. When I visualize the dock.inpcrd file the ligand has changed coordinates, meaning the starting structure before minimizing the system is not what I initially obtained from the docking program. Has anyone encountered similar problems, or has anyone any suggestion?
Looking forward to replies!
Best regards
Henrik Oberg
KTH, Sweden
The AMBER Mail Reflector
To post, send mail to
To unsubscribe, send "unsubscribe amber" (in the *body* of the email)
Received on Sun Aug 31 2008 - 06:07:03 PDT
Custom Search