problem with GluP ACe CT3 patch
hello,
I am trying to patch a protonated glu as well as ACE and CT3. Most of the documentation on patches is for DISU. As for ACE or CT3 (which is what I tried first) I used:
class mymodel(automodel): def special_patches(self, aln): self.patch(residue_type='CT3', residues=(self.residues['473:A']))
a = mymodel(env, alnfile = 'all.pir', knowns = ('3ob6_A','3ob6_B'), sequence = 'Can1_E184H')
a.starting_model= 1 a.ending_model = 1
I got the following error:
patch___482E> Patching residue identifier not found: CT3
Same thing for ACE. I also tried with allhmodel but it fails as well.
Another question for GLUP patching in no hydrogen models how can one be sure that one got indeed GLUP?
Thanks in advance for your help,
Martine.
-------------------------------------------- Martine Prévost Structure et Fonction des Membranes Biologiques CP 206/2 Université Libre de Bruxelles Bld du Triomphe B-1050 Bruxelles Belgium Phone : 32-2-650 2049 Fax : 32-2-650 5382 E-mail: mprevost@ulb.ac.be --------------------------------------------
On 4/20/16 2:31 AM, Prevost Martine wrote: > I am trying to patch a protonated glu as well as ACE and CT3. > Most of the documentation on patches is for DISU. > As for ACE or CT3 (which is what I tried first) I used:
Your script looks fine to me (obviously I can't run it without the rest of your inputs though). And all three patches work fine for me with both automodel and allhmodel. Is it possible you are using a different topology file that doesn't have these patches in it, or you are doing something else unusual?
> Another question for GLUP patching in no hydrogen models how can one be sure that one got indeed GLUP?
Look in the log file for the "segment topology patched using RTF" message. But in no hydrogen models the GLUP patch doesn't do anything except change the charges, which aren't used anyway, so it doesn't really matter whether or not it gets applied.
Ben Webb, Modeller Caretaker
participants (2)
-
Modeller Caretaker
-
Prevost Martine