Skip to content
New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

Fail to read P00533 #47

Open
bkellman opened this issue May 11, 2018 · 0 comments
Open

Fail to read P00533 #47

bkellman opened this issue May 11, 2018 · 0 comments

Comments

@bkellman
Copy link

My code:


PROJECT = 'GTspec_query'
LIST_OF_GENES = ['P01106','P01229','P01374']
PDB_FILE_TYPE = 'mmtf'

# Create the GEM-PRO project
my_gempro = GEMPRO(gem_name=PROJECT, root_dir=ROOT_DIR, genes_list=LIST_OF_GENES, pdb_file_type=PDB_FILE_TYPE)

# UniProt mapping
my_gempro.uniprot_mapping_and_metadata(model_gene_source='ACC+ID') 

fails and it appears to be due to a bad model in P00533. Error below

<ipython-input-2-d2f5408c2c3e> in <module>()
----> 1 execfile('02_structure_annote.py')

/media/ben/9c17f1c9-a45e-49ec-b547-8fbd2f25ccc6/GTspecificity/02_structure_annote.py in <module>()
    43
    44 # UniProt mapping
---> 45 my_gempro.uniprot_mapping_and_metadata(model_gene_source='ACC+ID')  ## [Nathan] can leave this as ACC+ID which means your list o
    46 print('Missing UniProt mapping: ', my_gempro.missing_uniprot_mapping)
    47 my_gempro.df_uniprot_metadata.head()

/home/ben/anaconda2/lib/python2.7/site-packages/ssbio/pipeline/gempro.pyc in uniprot_mapping_and_metadata(self, model_gene_source, custo
   546                         uniprot_prop = g.protein.load_uniprot(uniprot_id=mapped_uniprot, download=True, outdir=outdir,
   547                                                               set_as_representative=set_as_representative,
--> 548                                                               force_rerun=force_rerun)
   549                     except HTTPError as e:
   550                         log.error('{}, {}: unable to complete web request'.format(g.id, mapped_uniprot))

/home/ben/anaconda2/lib/python2.7/site-packages/ssbio/core/protein.pyc in load_uniprot(self, uniprot_id, uniprot_seq_file, uniprot_xml_f
   378             if download:
   379                 uniprot_prop.download_metadata_file(outdir=outdir, force_rerun=force_rerun)
--> 380                 uniprot_prop.download_seq_file(outdir=outdir, force_rerun=force_rerun)
   381
   382             # Also check if UniProt sequence matches a potentially set representative sequence

/home/ben/anaconda2/lib/python2.7/site-packages/ssbio/databases/uniprot.pyc in download_seq_file(self, outdir, force_rerun)
   183                                                    force_rerun=force_rerun)
   184
--> 185         self.sequence_path = uniprot_fasta_file
   186
   187     def download_metadata_file(self, outdir, force_rerun=False):

/home/ben/anaconda2/lib/python2.7/site-packages/ssbio/protein/sequence/seqprop.pyc in sequence_path(self, fasta_path)
Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment
Labels
None yet
Projects
None yet
Development

No branches or pull requests

1 participant