Leap does not add disulfide bonds automatically, you must
do this with the "bond" command.
jim
On Wed, 16 Apr 2003, Youyi Peng wrote:
> Hi,
> My protein contains a disulfide bond and a HIS is involved in the
> ligand binding. I renamed the corresponding CYS and HIS to CYX and HIE
> in the pdb file and then loadpdb into xleap and saved the prmtop and
> inpcrd files. S-S disulfide bond existed in inpcrd file. After energy
> minimization, S-S bond broke up and the distance is over 3 A in the
> min.restrt file. It seems to me that loadpdb treats CYS as a residue
> name but not the real CYS (including S-S bond). Am I Rright? what is the
> normal procedure to deal with CYX and HIE(D/P)?
> I have questions here:
> For CYS, even I changed the res. name, I still need to use BOND to make
> the bond in leap and loadpdbusingseq, Or loadpdb will treat CYX as the
> defined CYX by AMBER?
> For HIS, after I changed the name to HIE, when loaded into leap, it is
> treated as HIE or still the default HIS? Do I need to change the leaprc
> file?
> I am an AMBER beginner, If I ask some silly and simple questions,
> please don't feel mad at me.
> Thanks a lot.
> Youyi Peng
>
--
----------------------------------------------------------------------------
NOTE NEW MAILING ADDRESS & FAX NUMBER:
James W. Caldwell (voice) 415-476-8603
Department of Pharmaceutical Chemistry (fax) 415-514-3866
UCSF Box 2240 (email) caldwell.heimdal.ucsf.edu
600 16th Street
University of California
San Francisco, CA 94143-2240
----------------------------------------------------------------------------
Received on Wed Apr 16 2003 - 16:53:01 PDT
