On Wed, Aug 03, 2011, Maura Catherine Mooney wrote:
>
> The prmtop was constructed in xleap from a stripped pdb file. I created
> a lib file using RESP charges generated with amber then neutralized and
> solvated. The four crystal waters, I think, are important for metal
> coordination and I wanted to eventually treat this region with qm (along
> with a ligand). Therefore I thought it best to model these 4 waters
> with the 'normal' geometry, as opposed to the three-bond model.
This is wrong: no matter what you eventually intend to do with QM, in the MM
force field, if you are using TIPx waters, you need to make sure that the
waters are rigid.
> Also I attempted the rdparm command on my
> prmtop and I thought the output was strange...see below:
>
> ==========================================
> RDPARM MENU: angles :WAT
> Mask [:WAT] represents 28557 atoms
>
> Angle Kthet degrees atom names (numbers)
>
> ===========================================
This looks correct: LEaP removes all the H-O-H angle force field terms in
WAT residues, unless you issue the command "set default flexibleWater on".
...good luck....dac
_______________________________________________
AMBER mailing list
AMBER.ambermd.org
http://lists.ambermd.org/mailman/listinfo/amber
Received on Thu Aug 04 2011 - 05:30:03 PDT