AMBER Archive (2003)

Subject: Re: CYS and HIS

From: caldwell_at_heimdal.compchem.ucsf.edu
Date: Wed Apr 16 2003 - 09:55:55 CDT


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_at_heimdal.ucsf.edu 600 16th Street University of California San Francisco, CA 94143-2240 ----------------------------------------------------------------------------