Hi Prathvi,
Assuming that 'm' is your protein model, then the code would be:

from ChemGroup import findGroup
disulfides = findGroup("disulfide", [m])
num_disulfides = len(disulfides)

--Eric

Eric Pettersen
UCSF Computer Graphics Lab


On Nov 7, 2022, at 12:32 AM, Prathvi Singh via Chimera-users <chimera-users@cgl.ucsf.edu> wrote:

Hi Elaine,

I have found out a way to select the disulphide bonded CYS residues in UCSF chimera using the command: select disulphide.

However, is there a way to count the number of disulphide bridges in a protein and save in a variable?

--
Prathvi Singh,
Research Fellow,
Department of Biological Sciences & Bioengineering,
Indian Institute of Technology, Kanpur-208016
_______________________________________________
Chimera-users mailing list: Chimera-users@cgl.ucsf.edu
Manage subscription: https://www.rbvi.ucsf.edu/mailman/listinfo/chimera-users