Did you read the leap.log file?
On Oct 24, 2017 11:00 AM, "Michael Shokhen" <michael.shokhen.biu.ac.il>
wrote:
> Dear Amber experts,
>
>
> I have observed braking of disulfide bonds between CYX residues
>
> in protein during MD simulation by Amber16.
>
> I used ff14SB force field.
>
> During protein preparing I have defined in tleap input script file all CYX
>
> pairs connected by disulfide bonds (see below).
>
> I would appreciate your advice how to resolve the problem.
>
>
> Thank you,
>
> Michael
>
>
> Script for tleap:
>
>
>
> source leaprc.protein.ff14SB
> source leaprc.water.tip3p
> mc = loadPDB 1fle_nh.pdb
> bond mc.31.SG mc.47.SG
> bond mc.128.SG mc.195.SG
> bond mc.159.SG prot.175.SG
> bond mc.185.SG mc.215.SG
> bond mc.248.SG mc.277.SG
> bond mc.255.SG prot.281.SG
> bond mc.264.SG mc.276.SG
> bond mc.270.SG prot.285.SG
> addions2 mc Cl- 0
> solvateoct mc TIP3PBOX 15
> saveamberparm mc mc.prmtop mc.inpcrd
> quit
>
> _______________________________________________
> AMBER mailing list
> AMBER.ambermd.org
> http://lists.ambermd.org/mailman/listinfo/amber
>
_______________________________________________
AMBER mailing list
AMBER.ambermd.org
http://lists.ambermd.org/mailman/listinfo/amber
Received on Tue Oct 24 2017 - 08:30:03 PDT
