𧬠superdarks: Code files : 3-post_query_informatics
- π 7-clusterupdate_branch/
- π library/
- π 1-listing_parse.py β Download
- π 10-goFo_one_node.py β Download
- π 11-reduce_for_cytoscape_strict.py β Download
- π 12-clustalo_pairwise_traceback_map_con_to_structure.py β Download
- π 2-uniprot_parse.py β Download
- π 3-listing_uniprot_merge.py β Download
- π 4-uniprot_interpro_merge.py β Download
- π 5-gpcr_filtering.py β Download
- π 6-sequence_chunkify_and_collect.py β Download
- π 8-chunkify_and_blastp.py β Download
- π 9-blastp_chunks_to_network_fast.py β Download
π Scripts in this folder
1-listing_parse.pyβ Parsesrank / coverage / TM-score / UniProt-IDtuples out of therank--coverage-tm--UNIPROT-β¦-full.pdbfilenames produced upstream.2-uniprot_parse.pyβ Streaming parser for Swiss-Prot / TrEMBL.datflat-file entries, one protein record at a time (memory-safe).3-listing_uniprot_merge.pyβ Streams the full Swiss-Prot and TrEMBL.dat.gzarchives viapigz/gzipand merges them with the hit listing into a unified TSV (gene, organism, superkingdom, and protein-function columns).4-uniprot_interpro_merge.pyβ Joins the merged UniProt TSV with InterProβsprotein2ipr.dat.gz, emitting one combined TSV per hit (UPPERCASE-normalized for downstream matching).5-gpcr_filtering.pyβ Pandas-based 7TM + GPCR keyword filter over the InterPro-enriched TSV (with a specialGCR-\d+family rule), writing per-match Excel workbooks up to the 1,048,576-row sheet limit.6-sequence_chunkify_and_collect.pyβ Rewrites FASTA headers into the canonicalsuperkingdom|uniprot_id|gene|organism|rankform and chunks sequences into downstream-ready files.7-clusterupdate_branch/β Optional pre-BLASTP reduction: splits the FASTA by superkingdom and runs MMseqs2linclust+clusterupdateto produce priority representatives (Archaea β Bacteria β Eukaryota).8-chunkify_and_blastp.pyβ Splits one FASTA into N whole-record chunks and submits one all-vs-all BLASTP LSF job per chunk (no CLI β edit the CONFIG block).9-blastp_chunks_to_network_fast.pyβ Unified network builder (FAST in-process, or TURBO / LSF cluster) that turns per-chunk BLASTP TSVs into gzippednodes+edgestables ready for Cytoscape.10-goFo_one_node.pyβ Isom-Lab βgo-forwardβ BFS walk over the thresholded BLAST similarity graph (by % identity plus e-value or bitscore), producing a no-rings, strongest-edges-first subgraph.11-reduce_for_cytoscape_strict.pyβ Reduces a very large network to a Cytoscape-friendly subset by keeping the largest components, extracting a Maximum Spanning Forest backbone, and pruning by node-importance score.12-clustalo_pairwise_traceback_map_con_to_structure.pyβ Clustal-Omega pairwise traceback that maps sequence conservation onto a reference PDB, enabling 3D visualization of evolutionary signal.
β οΈ Disclaimer
This is research-grade code provided for academic and exploratory purposes only. It is not intended for professional or clinical deployment.