2015-04-17 1 views
0
from Bio import PDB 
pdb1 = PDB.PDBList() 
pdb1.retrieve_pdb_file("1CRK") #downloads the .pdb file from the internet 
parser = PDB.PDBParser(PERMISSIVE=1) 
structure = parser.get_structure("1CRK",r'C:\Users\Sam\AppData\Local\Enthought\Canopy\User\Lib\site-packages\Bio\cr\pdb1crk.ent') 
print(structure) 
print(dir(structure)) 
for model in structure: 
    for chain in model: 
     for residue in chain: 
      for atom in residue:         
       N = atom.get_name() 
       I = atom.get_id() 
       Y = atom.get_coord()     
       V = atom.get_vector() 
       O = atom.get_occupancy() 
       B = atom.get_bfactor()   

Je suis nouveau à python et essayé un aller au module biopython, Ce que je veux faire est pour fondamentalement, imprimer (Y, B) seulement pour les valeurs qui ont B> 10 et B < 50. Mais je n'ai aucune idée de la façon d'aborder cela. Aucune suggestion? :(this is where i got all the code from, this might helpBiopython extrait PDB coordonnées (x, y, z) et les données de valeur b si b valeur> 10 && <50

Répondre

0
if 10 < B < 50: 
    print(Y, B) 

coincé dans cette dernière boucle

+0

Je ne savais pas comment à l'époque, tout compris it out: D –