Just to inform of the way I found to deal with this problem. Basically, my
the ligand in my protein/ligand setup does not have atom types.
1. Extract ligand and protein coordinates separately in pdb (can use mol2).
--> ligand.pdb, protein.pdb
2. In sleap
>LIG = loadpdb ligand.pdb
>fixbond LIG
>addhydr LIG
>savemol2 LIG ligand.mol2
3. Using antechamber:
>antechamber -i ligand.pdb -fi pdb -o ligand+types.mol2 -fo mol2 -c bcc
-s 2
4. In tleap
>LIG = loadmol2 ligand+types.mol2
>PROT = loadpdb protein.pdb
>COMP = combine { LIG PROT }
fixbond and addhydr seem to work fine in the adding of the hydrogens. Since
antechamber can export as a mol2, we avoid the internal coordinates of the
prepi, and original coordinates are conserved.
Thanks to all who helped me out here!
-Sean
Sean Johnston
johnstos.stolaf.edu
sean.johnst.gmail.com
On Mon, Oct 20, 2008 at 5:39 PM, Bill Ross <ross.cgl.ucsf.edu> wrote:
> Once you have a prepin, the secret is to load it and then load your
> pdb with protein and ligand coordinates. These pdb coordinates will
> be applied to all residues, so the relative locations of protein and
> ligand will be preserved. Ultimately the whole system will be shifted
> when you solvate it, so do not get too attached to the coords in the
> pdb.
>
> Bill
> -----------------------------------------------------------------------
> The AMBER Mail Reflector
> To post, send mail to amber.scripps.edu
> To unsubscribe, send "unsubscribe amber" (in the *body* of the email)
> to majordomo.scripps.edu
>
-----------------------------------------------------------------------
The AMBER Mail Reflector
To post, send mail to amber.scripps.edu
To unsubscribe, send "unsubscribe amber" (in the *body* of the email)
to majordomo.scripps.edu
Received on Sun Nov 02 2008 - 05:08:34 PST
