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!