Dear George,
> I am experiencing the following problem when trying to save the parameter
> file of my protein.
>
> <saveamberparm> generates the following output:
>
> Could not find bond parameter for: SH - SH
>
> Building angle parameters.
> Could not find angle parameter: HS - SH - SH
> Could not find angle parameter: SH - SH - CT
> etc.
>
> I'm dealing with a small protein with 3 disulfide bonds.
>
> I'm using leaprc.ff99SB as a forcefield
Cysteines involved in disulphide bonds have to be renamed CYX in the pdb.
A CYS residue has types CT-SH-HS
while CYX has CT-S
So not only do you have additional hydrogens here that shouldn't be there
the types are incorrect for disulphide bonds. Edit your pdb and reload if
and all should be good.
All the best
Ross
/\
\/
|\oss Walker
| Assistant Research Professor |
| San Diego Supercomputer Center |
| Tel: +1 858 822 0854 | EMail:- ross.rosswalker.co.uk |
|
http://www.rosswalker.co.uk | PGP Key available on request |
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 mailing list
AMBER.ambermd.org
http://lists.ambermd.org/mailman/listinfo/amber
Received on Sun Jan 25 2009 - 01:11:59 PST
