On Fri, Apr 03, 2009, chenquan wrote:
>
> I am doing a simulation of protein. Ester bond forms between the
> carboxyl of one residue and hydroxy of threonine in the protein. when
> I run the command "saveamberparm AA cmp.prmtop cmp.inpcrd", there is
> always a error："Could not find angle parameter: CT - OH - C". The
> force fields I specfied is leaprc.ff03 and leaprc.glycam04EP.
There are no esters in standard amino acids, and hence no relevant
parameters in the standard parameter files. You will need to develop
your own force field parameters here. Check out tutorial B4 for one way
to begin.
>
> In additon, after I load the protein and run "savepdb", the ester bond
> is breaked which is treated as the N-terminal and C-terminal. In fact,
> the bond is connected before I load it.
It's hard to help, since you give so few details of what you actually
did. My guess is that you have not yet created a unit with the ester
bond in it; this will have a special residue name (not one of the
standard names). If you use that, the bonds should remain intact.
...dac
_______________________________________________
AMBER mailing list
AMBER.ambermd.org
http://lists.ambermd.org/mailman/listinfo/amber
Received on Sun Apr 05 2009 - 01:08:27 PDT
