Hi All,
Has anybody experienced problems with sleap and iwrap? (Is anybody using
sleap routinely?).
Attached is a pdb file and two scripts that 'should' do the same thing in
sleap and tleap. If you then run MD with the two resulting prmtop and inpcrd
files things appear to be fine without iwrap=0 but if you set iwrap=1 the
sleap generated prmtop run behaves VERY strangely.
You get the protein jumping out of the water box in a strange way. For this
case with just one molecule things appear to be ok, as in just a strange
imaging issue. However, if you set a system up with glycans attached then in
the sleap run only half the molecule gets imaged and then it crashes with
vlimit errors (due to bonds being massively stretched) as soon as you turn
on iwrap.
Hence it seems like there is something very wrong with the prmtop generated
by sleap. This is BAD since sleap is the only way to generate prmtops with
mixed 1-4 scaling parameters in right now.
Has anybody else seen such issues / has an idea of what might be wrong with
the sleap prmtop?
All the best
Ross
/\
\/
|\oss Walker
---------------------------------------------------------
| Assistant Research Professor |
| San Diego Supercomputer Center |
| Adjunct Assistant Professor |
| Dept. of Chemistry and Biochemistry |
| University of California San Diego |
| NVIDIA Fellow |
|
http://www.rosswalker.co.uk |
http://www.wmd-lab.org/ |
| Tel: +1 858 822 0854 | EMail:- ross.rosswalker.co.uk |
---------------------------------------------------------
Note: Electronic Mail is not secure, has no guarantee of delivery, may not
be read every day, and should not be used for urgent or sensitive issues.
_______________________________________________
AMBER-Developers mailing list
AMBER-Developers.ambermd.org
http://lists.ambermd.org/mailman/listinfo/amber-developers
Received on Fri Apr 08 2011 - 10:00:05 PDT