<p>You could run <a href="http://meme-suite.org/doc/fimo.html" rel="nofollow noreferrer">FIMO</a> on your entire genome for TFs (transcription factors; DNA-binding proteins) of interest, which gives you binding sites: genomic intervals where those TFs bind.</p>

<p><a href="https://bioinformatics.stackexchange.com/a/2491/776">https://bioinformatics.stackexchange.com/a/2491/776</a></p>

<p>You can then map or intersect promoter regions with those TF binding sites with BEDOPS <code>bedmap</code> or <code>bedops</code>, respectively:</p>

<p><a href="https://bedops.readthedocs.io/en/latest/content/reference/statistics/bedmap.html" rel="nofollow noreferrer">https://bedops.readthedocs.io/en/latest/content/reference/statistics/bedmap.html</a></p>

<p><a href="https://bedops.readthedocs.io/en/latest/content/reference/set-operations/bedops.html" rel="nofollow noreferrer">https://bedops.readthedocs.io/en/latest/content/reference/set-operations/bedops.html</a></p>

<p>The use of <code>bedmap</code> is recommended to assign TFs to a promoter of interest. This tool will return both the promoter region and any TFs that overlap it by the specified overlap criteria.</p>

CC BY-SA 4.0
2019-05-31T08:50:00.413Z
Thanks a lot for your very helpful response. "FIMO scans a set of sequences for individual matches to each of the motifs you provide ". So how do I relate ouput from FIMO to TF binding sites and the genomic intervals? Is there any method to get the specific names of the genes whose promoter regions bind my TF? I have used the FIMO web tool http://meme-suite.org/tools/fimo.
Charles
CC BY-SA 4.0
2019-05-31T12:20:22.153Z
I don't understand. If you have promoters, those are derived from genes, no?
Alex Reynolds
CC BY-SA 4.0
2019-05-31T14:03:22.213Z
Sorry for the silly question, but which column of FIMO output corresponds to the TF binding sites? I now want to  intersect promoter regions with those TF binding sites using BEDOPS bedops.
Charles
CC BY-SA 4.0
2019-05-31T14:31:45.637Z
The output from the command-line tool should be in BED format, so the binding site of the TF would be the first three columns. I don't know what format comes out of the web tool, but I would start by reading the documentation.
Alex Reynolds
CC BY-SA 4.0
2019-05-31T14:59:04.243Z
I have used the command line tool by running: fimo [options] <motif file> <sequence file> as found here http://meme-suite.org/doc/fimo.html. But none of the output from the above command is a .bed file. The output files include fimo.html, fimo.tsv, fimo.gff, cisml.xml and fimo.xml,  files. I want to be able to compare results from the command line and  results from the online tool, that is why I run the two methods.
Charles
CC BY-SA 4.0
2019-05-31T15:15:34.590Z
Maybe read the SE post I linked to at the top of my answer. This has a `fimo` command that writes to a BED file. Hope this helps!
Alex Reynolds
CC BY-SA 4.0
2019-06-04T16:33:29.660Z
Thank you once more. But  I discover so many repeats after running FIMO. The repeats also occur when I discover the motifs on my transcription factor DNA sequence (obtained using Ensembl), using the MEME tool. I set number of motifs to detect to 3 and size of motifs to 10bp. I generated the promoter sequences of my genes fed into FIMO using genome browser after setting promoter sequence to 1000bp above TSS. Please is there a way I can solve the problem of the repeats at any of the levels? The genes I use actually come from a Single cell RNA seq experiment.
Charles
CC BY-SA 4.0
2019-06-04T21:10:15.667Z
Is it normal to have a motif look like this:CACACACACACACACACACACACACACACACACACACACACACA
Charles
CC BY-SA 4.0
2019-06-05T00:38:46.230Z
It seems like you might be doing things out of order: You would feed the whole genome (mappable regions) into FIMO against Jaspar or other published motif databases, limiting your search to motifs no longer than, say, 20nt. Then do your set operations against this result with your promoters, at a desired statistical threshold (1e-5 or less). If you're looking to predict putative motifs from your promoters directly, then you would use MEME, not FIMO. Try re-reading the answer I gave.
Alex Reynolds
CC BY-SA 4.0
2019-06-06T06:35:13.630Z
Thanks a lot. I figured out many things. After running intersecting the TF binding sites with promoter sequences 500 bp above TSS, I obtain results like this: chr1 2000801 2000900 . I have searched the internet on how to decode these genomic coordinates into gene names to no avail. Is there a tool that i can use to do this? Secondly is there a possibility to have graphical output from running bedops? Thanks
Charles
CC BY-SA 4.0
2019-05-27T04:53:42.147Z
2019-05-28T17:51:20.840Z
2019-05-28T17:51:20.840Z
Alex Reynolds
<p>I have a list of human genes that have been proven to be up-regulated in a disease condition of interest, using microarray analysis. I also have a protein which according to literature is said to be a DNA binding protein. I wish to be able to computationally predict if the protein mentioned binds to any of the promoter sequences of my differentially expressed genes and if yes , then at which position?</p>

<p>Are there any Bioinformatic tools available for this task?</p>

<computation><homer>
How do I computationally predict the binding position of a DNA binding protein on the promoters of differentially expressed genes?
# How do I computationally predict the binding position of a DNA binding protein on the promoters of differentially expressed genes?

(Asked by: username_0 on 2019-05-26)

I have a list of human genes that have been proven to be up-regulated in a disease condition of interest, using microarray analysis. I also have a protein which according to literature is said to be a DNA binding protein. I wish to be able to computationally predict if the protein mentioned binds to any of the promoter sequences of my differentially expressed genes and if yes , then at which position?

Are there any Bioinformatic tools available for this task?

**Tags:** computation, homer

---

## 1 answers

### Answer by username_1 on 2019-05-27. Score: 4

You could run [FIMO](http://meme-suite.org/doc/fimo.html) on your entire genome for TFs (transcription factors; DNA-binding proteins) of interest, which gives you binding sites: genomic intervals where those TFs bind.

<https://bioinformatics.stackexchange.com/a/2491/776>

You can then map or intersect promoter regions with those TF binding sites with BEDOPS `bedmap` or `bedops`, respectively:

<https://bedops.readthedocs.io/en/latest/content/reference/statistics/bedmap.html>

<https://bedops.readthedocs.io/en/latest/content/reference/set-operations/bedops.html>

The use of `bedmap` is recommended to assign TFs to a promoter of interest. This tool will return both the promoter region and any TFs that overlap it by the specified overlap criteria.

> **Comments:**
> **username_0**: Thanks a lot for your very helpful response. "FIMO scans a set of sequences for individual matches to each of the motifs you provide ". So how do I relate ouput from FIMO to TF binding sites and the genomic intervals? Is there any method to get the specific names of the genes whose promoter regions bind my TF? I have used the FIMO web tool http://meme-suite.org/tools/fimo.  
> **username_1**: I don't understand. If you have promoters, those are derived from genes, no?  
> **username_0**: Sorry for the silly question, but which column of FIMO output corresponds to the TF binding sites? I now want to  intersect promoter regions with those TF binding sites using BEDOPS bedops.  
> **username_1**: The output from the command-line tool should be in BED format, so the binding site of the TF would be the first three columns. I don't know what format comes out of the web tool, but I would start by reading the documentation.  
> **username_0**: I have used the command line tool by running: fimo [options] <motif file> <sequence file> as found here http://meme-suite.org/doc/fimo.html. But none of the output from the above command is a .bed file. The output files include fimo.html, fimo.tsv, fimo.gff, cisml.xml and fimo.xml,  files. I want to be able to compare results from the command line and  results from the online tool, that is why I run the two methods.  
> **username_1**: Maybe read the SE post I linked to at the top of my answer. This has a `fimo` command that writes to a BED file. Hope this helps!  
> **username_0**: Thank you once more. But  I discover so many repeats after running FIMO. The repeats also occur when I discover the motifs on my transcription factor DNA sequence (obtained using Ensembl), using the MEME tool. I set number of motifs to detect to 3 and size of motifs to 10bp. I generated the promoter sequences of my genes fed into FIMO using genome browser after setting promoter sequence to 1000bp above TSS. Please is there a way I can solve the problem of the repeats at any of the levels? The genes I use actually come from a Single cell RNA seq experiment.  
> **username_0**: Is it normal to have a motif look like this:CACACACACACACACACACACACACACACACACACACACACACA  
> **username_1**: It seems like you might be doing things out of order: You would feed the whole genome (mappable regions) into FIMO against Jaspar or other published motif databases, limiting your search to motifs no longer than, say, 20nt. Then do your set operations against this result with your promoters, at a desired statistical threshold (1e-5 or less). If you're looking to predict putative motifs from your promoters directly, then you would use MEME, not FIMO. Try re-reading the answer I gave.  
> **username_0**: Thanks a lot. I figured out many things. After running intersecting the TF binding sites with promoter sequences 500 bp above TSS, I obtain results like this: chr1 2000801 2000900 . I have searched the internet on how to decode these genomic coordinates into gene names to no avail. Is there a tool that i can use to do this? Secondly is there a possibility to have graphical output from running bedops? Thanks  

---

