Kegg soap enzymes by pw.pl

 
#!/usr/bin/perl
 
use SOAP::Lite;
use strict;
 
# this script is largely inspired from the sample scripts published by KEGG
# it will return all yeast 'sce' pathways
 
my $wsdl = 'http://soap.genome.jp/KEGG.wsdl';
my $server = SOAP::Lite -> service($wsdl);
my $pathcount=0;
 
foreach my $path (@{$server -> list_pathways("sce")}) {
 $pathcount++;
 my $keggID = $path -> {entry_id};
 $keggID =*tilde* s/path://;
 print "Pathway: $pathcount\t$keggID\t$path->{definition}\n";
 my $genecount=0;
 foreach my $enzyme ( @{$server -> get_enzymes_by_pathway($path -> {entry_id})}) {
  $genecount++;
  my @gene_array=@{$server -> get_genes_by_enzyme($enzyme, 'sce')}; 
  my $genes = join (";", @gene_array);
  $genes =*tilde* s/sce//g;
  print "Pathway-member: $genecount\t$enzyme\t$genes\n";
 }
}