Init the library and needed functions.
In [6]:
import rna_tools.Seq as Seq
import rna_tools.BlastPDB
from rna_tools.BlastPDB import BlastPDB
reload(rna_tools.BlastPDB);
reload(Seq);
Create a RNASeqence object.
In [7]:
seq = Seq.RNASequence("GGGUCAGGCCGGCGAAAGUCGCCACAGUUUGGGGAAAGCUGUGCAGCCUGUAACCCCCCCACGAAAGUGGG")
print(seq)
In [8]:
print(seq.predict_ss())
In [9]:
print(seq.predict_ss(method='RNAsubopt'))
In [10]:
print(seq.predict_ss(method='ipknot'))
In [11]:
print(seq.predict_ss(method='centroid_fold'))
In [29]:
p = BlastPDB(seq.seq)
p.search()
print p.result
In [30]:
import rna_pdb_tools.RfamSearch as rf
#reload(rf)
#seq = Seq.Seq("GGGUCAGGCCGGCGAAAGUCGCCACAGUUUGGGGAAAGCUGUGCAGCCUGUAACCCCCCCACGAAAGUGGG")
rs = rf.RfamSearch()
print rs.cmscan(seq)
In [ ]:
from rna_pdb_tools.pdb_parser_lib import RNAStructure
fn = "rna_pdb_tools/data/260c8ff6-f24e-4eff-9760-1831407fc770_ALL_thrs5.30A_clust01-000001_AA.pdb"
s = RNAStructure(fn)
print s.get_report()
print s.get_info_chains()
print s.get_head()
#print s.view() # image paste here :-)
In [ ]:
%%bash
cd rna_pdb_tools
./rna-pdb-tools.py --no_hr --get_seq data/260c8ff6-f24e-4eff-9760-1831407fc770_ALL_thrs5.30A_clust01-000001_AA.pdb
In [ ]:
# model using SimRNA
#res = SimRNA(ss,seq.get_ss())
In [ ]:
# fake import, should be
res = "rna_pdb_tools/data/260c8ff6-f24e-4eff-9760-1831407fc770_ALL_thrs5.30A_clust01-000001_AA.pdb"
# view
view = nglview.show_structure_file(res)
view.add_representation(repr_type='cartoon')
view
In [ ]:
%%bash
cd ~/rna-bench/opt/xxxcx rna_pdb_tools
./rna-pdb-tools.py --no_hr --get_seq ~/rna-bench/examples/5k7c.pdb
In [ ]:
%%bash
cd rna_pdb_tools
./rna-pdb-tools.py --no_hr --get_seq input/5k7c.pdb
./rna-pdb-tools.py --no_hr --get_seq input/tetraloop.pdb
./rna-pdb-tools.py --get_seq input/1xjr.pdb
In [ ]: