Python script for pdb databases

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

Popularne posty z tego bloga

Kubernetes

Helm

Ansible Tower / AWX