Skip to content
Snippets Groups Projects
Commit 249e94b8 authored by masikora's avatar masikora
Browse files

fix bad atom selection for tutorial

parent 66bba25d
Branches main
No related tags found
No related merge requests found
...@@ -132,13 +132,13 @@ class glycoshield: ...@@ -132,13 +132,13 @@ class glycoshield:
if self.initialsugarframes[isugar] == 0: if self.initialsugarframes[isugar] == 0:
self.initialsugarframes[isugar] = self.usugar.trajectory.n_frames self.initialsugarframes[isugar] = self.usugar.trajectory.n_frames
# select (part of) the protein that is connected to the sugar # select (part of) the protein that is connected to the sugar
tripep = self.usugar.select_atoms('protein and resid {} and name N CA C O'.format(" ".join([str(i) for i in resids_on_sugar]))) tripep = self.usugar.select_atoms('protein and resid {} and name in N CA C O OT1'.format(" ".join([str(i) for i in resids_on_sugar])))
# Test whether sequon is in the 2:N-1 range: # Test whether sequon is in the 2:N-1 range:
self._test_sequon(resids_on_protein, protchain) self._test_sequon(resids_on_protein, protchain)
# select sequon on protein # select sequon on protein
sequon = self.uprot.select_atoms('segid {} and resid {} and name N CA C O'.format(protchain, " ".join([str(i) for i in resids_on_protein]))) sequon = self.uprot.select_atoms('segid {} and resid {} and name in N CA C O OT1'.format(protchain, " ".join([str(i) for i in resids_on_protein])))
# the selection has to have same length for protein and peptide with the sugar! # the selection has to have same length for protein and peptide with the sugar!
# Now select part of the peptide that should be transplanted onto protein. By default we do not transplant the part which gets selected. # Now select part of the peptide that should be transplanted onto protein. By default we do not transplant the part which gets selected.
......
0% Loading or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment