bio-octopus 0.1.0 → 0.1.1

Sign up to get free protection for your applications and to get access to all the features.
@@ -2,6 +2,15 @@
2
2
 
3
3
  A bioruby plugin for running and parsing octopus, available from http://octopus.cbr.su.se/ Currently, the plugin is not Ruby 1.9 (only 1.8) compatible, because of the dependency on the rio gem.
4
4
 
5
+ To run it as a script, a fasta file of sequences to predict, and a database formatted in legacy blast database format needs to be specified:
6
+
7
+ $ /usr/local/rvm/gems/ruby-1.8.7-p334/gems/bio-octopus-0.1.0/bin/biooctopus my.fasta /blastdb/UniProt2011_04/uniprot_sprot.fasta
8
+
9
+ Output are topology predictions, including signal peptide and transmembrane domain predictions. If that is not the best format for you, look at the script file and modify it - it is just a front-end. Also required are:
10
+ * a working install of SPOCTOPUS
11
+ * NCBI legacy blast
12
+ Binaries for these must be accessible in the PATH.
13
+
5
14
  == Contributing to bio-octopus
6
15
 
7
16
  * Check out the latest master to make sure the feature hasn't been implemented or the bug hasn't been fixed yet
data/VERSION CHANGED
@@ -1 +1 @@
1
- 0.1.0
1
+ 0.1.1
@@ -1,7 +1,7 @@
1
1
  #!/usr/bin/env ruby
2
2
  #
3
3
  # BioRuby bio-octopus Plugin
4
- # Version 0.0.1
4
+ # Version 0.1.1
5
5
  # Author:: Ben J Woodcroft
6
6
  # Copyright:: 2011
7
7
  # License:: The Ruby License
@@ -21,12 +21,12 @@ if $0 == __FILE__
21
21
  ].flatten
22
22
  end
23
23
  o.parse!
24
-
24
+
25
25
  unless ARGV.length == 2
26
- puts o.help
27
- exit
26
+ puts o.help
27
+ exit
28
28
  end
29
-
29
+
30
30
  fasta_path = ARGV[0]
31
31
  blastdb_path = ARGV[1]
32
32
 
@@ -36,6 +36,21 @@ if $0 == __FILE__
36
36
  result = runner.calculate(seq.seq, blastdb_path)
37
37
  name = seq.definition
38
38
 
39
+ # Output singal peptide if it has one
40
+ if result.signal?
41
+ puts [
42
+ name,
43
+ 'Signal peptide',
44
+ result.signal_peptide.start,
45
+ result.signal_peptide.stop
46
+ ].join("\t")
47
+ else
48
+ puts [
49
+ name,
50
+ 'No signal peptide found'
51
+ ].join("\t")
52
+ end
53
+
39
54
  if result.has_domain?
40
55
  # At least one TMD found. Output each on a separate line
41
56
  result.transmembrane_domains.each do |tmd|
@@ -50,7 +65,7 @@ if $0 == __FILE__
50
65
  else
51
66
  puts [
52
67
  name,
53
- 'No Transmembrane Domain Found'
68
+ 'No transmembrane domain found'
54
69
  ].join("\t")
55
70
  end
56
71
  end
@@ -1,7 +1,7 @@
1
1
  #!/usr/bin/env ruby
2
2
  #
3
3
  # BioRuby bio-octopus Plugin
4
- # Version 0.0.1
4
+ # Version 0.1.1
5
5
  # Author:: Ben J Woodcroft
6
6
  # Copyright:: 2011
7
7
  # License:: The Ruby License
@@ -21,12 +21,12 @@ if $0 == __FILE__
21
21
  ].flatten
22
22
  end
23
23
  o.parse!
24
-
24
+
25
25
  unless ARGV.length == 2
26
- puts o.help
27
- exit
26
+ puts o.help
27
+ exit
28
28
  end
29
-
29
+
30
30
  fasta_path = ARGV[0]
31
31
  blastdb_path = ARGV[1]
32
32
 
@@ -36,6 +36,21 @@ if $0 == __FILE__
36
36
  result = runner.calculate(seq.seq, blastdb_path)
37
37
  name = seq.definition
38
38
 
39
+ # Output singal peptide if it has one
40
+ if result.signal?
41
+ puts [
42
+ name,
43
+ 'Signal peptide',
44
+ result.signal_peptide.start,
45
+ result.signal_peptide.stop
46
+ ].join("\t")
47
+ else
48
+ puts [
49
+ name,
50
+ 'No signal peptide found'
51
+ ].join("\t")
52
+ end
53
+
39
54
  if result.has_domain?
40
55
  # At least one TMD found. Output each on a separate line
41
56
  result.transmembrane_domains.each do |tmd|
@@ -50,7 +65,7 @@ if $0 == __FILE__
50
65
  else
51
66
  puts [
52
67
  name,
53
- 'No Transmembrane Domain Found'
68
+ 'No transmembrane domain found'
54
69
  ].join("\t")
55
70
  end
56
71
  end
@@ -5,11 +5,11 @@
5
5
 
6
6
  Gem::Specification.new do |s|
7
7
  s.name = %q{bio-octopus}
8
- s.version = "0.1.0"
8
+ s.version = "0.1.1"
9
9
 
10
10
  s.required_rubygems_version = Gem::Requirement.new(">= 0") if s.respond_to? :required_rubygems_version=
11
11
  s.authors = ["Ben J Woodcroft"]
12
- s.date = %q{2011-04-17}
12
+ s.date = %q{2011-05-02}
13
13
  s.description = %q{Running and parsing of the protein transmembrane domain predictor octopus}
14
14
  s.email = %q{gmail.com after donttrustben}
15
15
  s.executables = ["biooctopus.rb", "biooctopus"]
metadata CHANGED
@@ -1,13 +1,13 @@
1
1
  --- !ruby/object:Gem::Specification
2
2
  name: bio-octopus
3
3
  version: !ruby/object:Gem::Version
4
- hash: 27
4
+ hash: 25
5
5
  prerelease:
6
6
  segments:
7
7
  - 0
8
8
  - 1
9
- - 0
10
- version: 0.1.0
9
+ - 1
10
+ version: 0.1.1
11
11
  platform: ruby
12
12
  authors:
13
13
  - Ben J Woodcroft
@@ -15,7 +15,7 @@ autorequire:
15
15
  bindir: bin
16
16
  cert_chain: []
17
17
 
18
- date: 2011-04-17 00:00:00 +10:00
18
+ date: 2011-05-02 00:00:00 +10:00
19
19
  default_executable:
20
20
  dependencies:
21
21
  - !ruby/object:Gem::Dependency