I have a little more information about my disulfide bond problem. Dr.
Case mentioned that the problem with LEaP crashing after the "bond"
command to form disulfides might mean that these bonds are already
present. I did a little investigating. When I load my pdb file into LEaP
and view it in the edit window, all of the disulfide bonds are properly
connected. If I just avoid using the "bond" command, saveAmberParm
writes the *.top and *.crd files normally. If, while still in LEaP, I
save the pdb file (savePdb), when I look at the pdb in Chimera the
disulfide bonds are broken. So, somewhere in the process of writing that
pdb file from LEaP the information that is needed to connect the
disulfide bonds is lost. The pdb file that is saved from LEaP is
stripped of SSBOND and CONNECT lines, so maybe it's no surprise that
these S - S connectivities are lost.
My question is this- should I just assume, since the disulfide bonds are
present in the edit screen in LEaP, that these bonds are really there
and that I don't need to worry about using the "bond" procedure for SG -
SG?
Another problem is that when I convert the *.rst files output by
minimization or MD to pdb files, once again the disulfides show up as
not connected. The two sulfur atoms, however, seem to stay in close
relationship to each other.
Thanks,
John
-----------------------------------------------------------------------
The AMBER Mail Reflector
To post, send mail to amber.scripps.edu
To unsubscribe, send "unsubscribe amber" to majordomo.scripps.edu
Received on Sun Mar 04 2007 - 06:07:37 PST
