Hello everybody,

I am having a problem with my model. I am new to Modeller but I googled and couldn't find an answer. I am selecting the residues I want to form the disulphide bonds but there is a one cys that gets in a way and when I visualize the disulphides in Pymol, I can see three cysteines are bonded together. How can I prevent that?

Here are my patches:

class MyModel(automodel):     def special_patches(self, aln):         self.patch(residue_type='DISU', residues=(self.residues['76'], self. residues['241']))         self.patch(residue_type='DISU', residues=(self.residues['46'], self. residues['122']))         self.patch(residue_type='DISU', residues=(self.residues['188'], self.residues['206']))         self.patch(residue_type='DISU', residues=(self.residues['105'], self.residues['202']))         self.patch(residue_type='DISU', residues=(self.residues['265'], self.residues['228']))         self.patch(residue_type='DISU', residues=(self.residues['113'], self.residues['183']))         self.patch(residue_type='DISU', residues=(self.residues['127'], self.residues['171']))

Since Modeller only reads ATOM and HETEATOM records I dont really understand why is this happening.

Thank you for any help!