/usr/share/pyshared/Bio/SCOP/Des.py is in python-biopython 1.58-1.
This file is owned by root:root, with mode 0o644.
The actual contents of the file can be viewed below.
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60 61 62 63 64 65 66 67 68 69 70 71 72 73 74 75 76 77 78 79 80 81 82 83 84 85 86 87 88 89 90 91 | # Copyright 2001 by Gavin E. Crooks. All rights reserved.
# This code is part of the Biopython distribution and governed by its
# license. Please see the LICENSE file that should have been included
# as part of this package.
""" Handle the SCOP DEScription file.
The file format is described in the scop
"release notes.":http://scop.berkeley.edu/release-notes-1.55.html
The latest DES file can be found
"elsewhere at SCOP.":http://scop.mrc-lmb.cam.ac.uk/scop/parse/
"Release 1.55":http://scop.berkeley.edu/parse/des.cla.scop.txt_1.55 (July 2001)
"""
class Record(object):
"""Holds information for one node in the SCOP hierarchy.
sunid -- SCOP unique identifiers
nodetype -- One of 'cl' (class), 'cf' (fold), 'sf' (superfamily),
'fa' (family), 'dm' (protein), 'sp' (species),
'px' (domain). Additional node types may be added.
sccs -- SCOP concise classification strings. e.g. b.1.2.1
name -- The SCOP ID (sid) for domains (e.g. d1anu1),
currently empty for other node types
description -- e.g. "All beta proteins","Fibronectin type III",
"""
def __init__(self, line=None):
self.sunid = ''
self.nodetype = ''
self.sccs = ''
self.name = ''
self.description =''
if line:
self._process(line)
def _process(self, line):
"""Parses DES records.
Records consist of 5 tab deliminated fields,
sunid, node type, sccs, node name, node description.
"""
#For example ::
#
#21953 px b.1.2.1 d1dan.1 1dan T:,U:91-106
#48724 cl b - All beta proteins
#48725 cf b.1 - Immunoglobulin-like beta-sandwich
#49265 sf b.1.2 - Fibronectin type III
#49266 fa b.1.2.1 - Fibronectin type III
line = line.rstrip() # no trailing whitespace
columns = line.split("\t") # separate the tab-delineated cols
if len(columns) != 5:
raise ValueError("I don't understand the format of %s" % line)
sunid, self.nodetype, self.sccs, self.name, self.description = columns
if self.name=='-': self.name =''
self.sunid = int(sunid)
def __str__(self):
s = []
s.append(self.sunid)
s.append(self.nodetype)
s.append(self.sccs)
if self.name:
s.append(self.name)
else:
s.append("-")
s.append(self.description)
return "\t".join(map(str,s)) + "\n"
def parse(handle):
"""Iterates over a DES file, returning a Des record for each line
in the file.
Arguments:
handle -- file-like object
"""
for line in handle:
if line.startswith('#'):
continue
yield Record(line)
|