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