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.ambermd.org
http://lists.ambermd.org/mailman/listinfo/amber
Received on Mon Jun 01 2009 - 01:08:31 PDT
