Re: [AMBER] disulfide bonds braking during MD simulation

From: David A Case <>
Date: Tue, 24 Oct 2017 20:39:09 -0400

On Tue, Oct 24, 2017, Michael Shokhen wrote:
> I have observed braking of disulfide bonds between CYX residues

This implies that the bonds were never properly created.
> 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
                 ^^^^ ???

What is "prot"? It doesn't seem to have been defined, and it would not
be surprising if the expected bond were breaking.

As Carlos pointed out, this problem is almost certainly noted in the tleap
output, and it's good to learn what to look for. Another good idea (TM)
is to use the "checkvalidity" command in parmed: it searches for missing
disulfide bonds.

But to developers: this points up two genuine problems with tleap:

  1. The output is so voluminous that people understandably have a hard
       time reading it.
  2. tleap really should quit when errors like this are encounterd. I
       understand that this would break some existing scripts, but the
       benefits are likely to outweigh the disadvantages.

I suppose I should be the one to attack this, but it would be really great
if someone would step up to the plate (this being World Series time....)
It's really pretty trivial coding to find these and other error messages
and make them fatal. Removing un-needed messages is more of an art, but
is still super-straightforward.


AMBER mailing list
Received on Tue Oct 24 2017 - 18:00:02 PDT
Custom Search