import os
import sys
import urllib2
import string
#url = 'http://www.rcsb.org/pdb/rest/search'
chain = sys.argv[2]
id = sys.argv[1]
d1 = id+"."+chain
url1 = 'http://www.pdb.org/pdb/rest/describeMol?structureId='
url = url1+d1
req = urllib2.Request(url)
f = urllib2.urlopen(req)
result = f.read()
#print "Found number of PDB entries:", result.count('\n')
wynik = []
wynik = result
b= string.split(wynik,'\n')
for i in range (len(b)):
if '<polymerDescription description' in b[i]:
wynik = b[i]
a=[]
c= string.join(wynik,'')
for i in range (len(c)):
#print c[i]
if '"' is c[i]:
a.append(i)
end=[]
for i in range (len(c)):
if i<a[1] and i>a[0]:
end.append(c[i])
print string.join(end,'')
Komentarze
Prześlij komentarz