Hello,

To define the disulfide bonds sometimes is upsetting, so I wrote a loop to set it up automatically:

Code:
set count = 0

label cysloop
incr count by 1
define residue sele ires @count end

if ?selresn .eq. CYS then
coor mind sele ires @count .and. type SG end sele type SG end
  if ?mind .lt. 2.5 then
    define CYS1 sele bynu ?minda1 end
    set CYS1 = ?selsegi ?selresi
    define CYS2 sele bynu ?minda2 end
    set CYS2 = ?selsegi ?selresi
    patch disu @CYS1 @CYS2 setup warn
  endif
endif

if count .lt. ?nres goto cysloop


I put the loop right after reading the coordinates/sequence of all chains.

If you all have any thoughts to improve this loop, it would be great.

Regards,

Antoniel.