[Chimera-users] Visualising disulfide bonds
meng at cgl.ucsf.edu
Sat Nov 12 08:42:46 PST 2011
If the bonds were literally present in the trajectory (defined in the prmtop), they would be shown if you were showing Cys atoms and bonds, and you could select them using "Select... Chemistry... functional group... disulfide" or command: sel disulfide
However, if they are not explicitly defined as bonds and you simply want to identify when cysteine sidechains come close to each other, you would need to do distance measurements instead. It wouldn't work to find cysteines within some zone of a cysteine since of course, each would be within the zone of itself, so you would need to either (a) use residue numbers and measure the SG-SG distance for each possible pair, or (b) use a python script that does the zone "except self." Neither of these would check for reasonable angles for bonding, however.
For (b) there is a similar script "selLys.py" available on our wiki
You would need to change it from using lysine NZ and 8 angstroms to cysteine SG and whatever cutoff you feel appropriate.
The command or script could be used on your trajectory on a per-frame basis (see Per-Frame menu of MD Movie).
I hope this helps,
Elaine C. Meng, Ph.D.
UCSF Computer Graphics Lab (Chimera team) and Babbitt Lab
Department of Pharmaceutical Chemistry
University of California, San Francisco
On Nov 12, 2011, at 3:33 AM, George Tzotzos wrote:
> Hi everybody,
> I'm trying to find a way to visualise disulfide bonds in AMBER trajectory files to no avail.
> Any help how to do so will be very much appreciated.
More information about the Chimera-users