1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23
|
model_file='sdh.pdb'
e=io.LoadEntity(model_file)
def StripSidechain(edi, residue):
atoms = residue.GetAtomList()
for atom in atoms:
if atom.name not in ['CA', 'C', 'N', 'O']:
edi.DeleteAtom(atom)
def StripSidechains(entity):
edi=entity.EditXCS(mol.EditMode.BUFFERED_EDIT)
for c in e.GetChainList():
residues=c.GetResidueList()
for residue in residues:
StripSidechain(edi, residue)
StripSidechains(e)
g=gfx.Entity('bare backbone', e)
scene.Add(g)
scene.CenterOn(g)
|