AMBER Archive (2009)

Subject: Re: [AMBER] broken Ligand

From: David A. Case (case_at_biomaps.rutgers.edu)
Date: Sun May 31 2009 - 07:50:57 CDT


On Sat, May 30, 2009, Shubhra Gupta wrote:
>
> Sorry for sparse information. I have attached the files. a.pdb file I got
> after saving topology (a.prmtop) and parameters (a.inpcrd). I Sourced
> leaprc.GLYCAM_06 and leaprc.gaff. For solvation I used TIP3PBOX and 7.0
> buffer. min.in was input file to execute sander for minimization. Resulted
> a_min.pdb is after minimization, where the ligand is broken.

I don't see what is "broken" in the a_min.pdb file, but that is OK:
since you know what the problem seems to be, you can debug for yourself.
Your files show a periodic system with explicit water, but the min.in file is
for a non-periodic system using GB. Try doing just a gas-phase minimization
on the 4YA residue all by itself. If you have problems, make is smaller until
things work, and then try to build back up to the full ligand molecule.

...dac

_______________________________________________
AMBER mailing list
AMBER_at_ambermd.org
http://lists.ambermd.org/mailman/listinfo/amber