[AMBER] solvate with crystal water

From: Hashem Taha <htahaphd.gmail.com>
Date: Mon, 11 Apr 2011 16:09:58 -0400

Hi,

I am trying to solvate a protein that has a probe attached to a cysteine
residue and includes one crystal water molecule. The leap input file I am
using is included below. Basically, when I try to solvate the system with
the crystal water included in the initial pdb, I get the message
"ContainerCenterAt(): Skipping atom w/ unknown position: O:O" and the system
does not solvate properly (the number of water residues that leap adds
is 169969). If I try removing the crystal water, solvateOct works fine, and
solvateOct also works without problems on the protein without the probe on
the cysteine residue. Please advise!

source /usr/local/amber/amber10/dat/leap/cmd/leaprc.gaff
source /usr/local/amber/amber10/dat/leap/cmd/leaprc.ff03
loadamberparams /nfs/01/osu6030/libs/CY3.frcmod
loadOff /nfs/01/osu6030/libs/CY3.lib
loadamberparams /nfs/01/osu6030/libs/linkers/LK5.frcmod
loadOff /nfs/01/osu6030/libs/linkers/LK5.lib
loadamberparams /nfs/01/osu6030/libs/TAMs/4S2.frcmod
loadOff /nfs/01/osu6030/libs/TAMs/4S2.lib
loadamberparams /nfs/01/osu6030/libs/linkers/linktoCYS.frcmod

p = loadpdb test.pdb
check p
set p head p.63.SG

b = sequence {4S2 LK5 p}
check b

addions b Cl- 0
solvateOct b TIP3PBOX 8
savepdb b test.pdb
saveamberparm b test.prmtop test.inpcrd

quit
_______________________________________________
AMBER mailing list
AMBER.ambermd.org
http://lists.ambermd.org/mailman/listinfo/amber
Received on Mon Apr 11 2011 - 13:30:04 PDT
Custom Search