This is a test version of Biostars. For the public version, visit https://www.biostars.org.
How to create protein-protein interaction network for specific genes (5000 genes) using biogrid?

Hi all,

I am doing chip-sea analysis. I have a set of 500 genes which correspond to the location of peaks. I was wondering how can I create a "protein-protein interaction network" for these genes and visualise with biogrid?

Kind regards

protein-protein-interaction cytoscape biogrid

If you're interested in direct interactions between your 5000 genes, then you could just filter out the biogrid interactions for which both partners are not in your 5000. If you're also interested in indirect interactions, that's a bit more complicated.

Thanks very much for your prompt response. Yes I am interested in direct interaction for this stage. So I will search about how to do the filtering in Biogrid. A link to appropriate tutorial would be highly appreciated. thanks again

If you download the mitab file corresponding to your organism (or the BIOGRID-ALL mitab if you have genes from multiple organisms), you'll see that the first two columns contain the Entrez gene IDs of the interaction partners, columns 3 and 4 contain the gene symbols. So you could use a short perl script like this (untested):

  open (FH,"<",$BIOGRID_file) or die ...
  while (my $line=<FH>) {
    chomp($line);
    my @tmpary = split(/\t/,$line);
    my @G1s = split(/\|/,$tmpary[2]); # gene symbol of partner 1
    my ($g1) = $G1s[0]=~/:(.+)/;
    my @G2s = split(/\|/,$tmpary[3]); # gene symbol of partner 2
    my ($g2) = $G2s[0]=~/:(.+)/;
    if ($wanted{$g1} && $wanted{$g2}) {
      # Assumes the symbols of the 5000 genes are keys to %wanted
      print "$g1\t$g2\n";
    }
  }
  close FH;

0 answers

No answers yet.

Log in to answer this question.