On Wed, Sep 09, 2009, Dean Cuebas wrote:
>
> One of the tutorials at the ambermd site shows the preparation of a ligand
> and it's solvation, but does not show how to properly setup a system when
> the coordinates of the ligand and protein are known.
It seems obvious after you know how....here is a typical tleap workflow:
source leaprc.ff99SB (for the protein part)
source leaprc.gaff (for the ligand part)
loadMol2 ligand.mol2 (where ligand.mol2 is created by antechamber) OR
loadAmberPrep ligand.prepi (where ligand.prepi is created by antechamber)
loadAmberParams frcmod (to get ligand-specific parameters, if any, created
by parmchk)
complex = loadPDB complex.pdb (where complex.pdb has consistent coordinates
of both the protein and the ligand, with the
ligand residue and atom names the same as in
the mol2 or prepi files loaded above. You
should have a TER card separating protein and
ligand.)
...solvate, etc.
saveAmberParm complex complex.top complex.crd
quit
Note that Amber does not create the complex.pdb file -- you get that from
a docking program, or from the PDB, or from wherever. LEaP will use whatever
coordinates are in that file, so check things visually before proceeding.
[Note: things are rather more complicated if there is a covalent bond between
the protein and the ligand.]
...hope this helps...dac
_______________________________________________
AMBER mailing list
AMBER.ambermd.org
http://lists.ambermd.org/mailman/listinfo/amber
Received on Mon Sep 14 2009 - 13:38:43 PDT
