This is a test version of Biostars. For the public version, visit https://www.biostars.org.
transcription binding site prediction from fasta sequence

Hi all I have 2 kb long sequence which is promoter region of a gene. From that sequence I would like to pinpoint transcription factor binding sites of smad3. How can I do that and is it possible with JASPAR?

Thanks

transcription site jaspar
• 1,720 views • link
updated by ATpoint ATpoint • written by newbio • newbio

Hi, you can use python to find all the smad3 binding sites in your query sequence, for example if smad3 binding site is 5'-GTCTAGAC-3' you can use the following dirty python code to get all your sites

#usr/bin/env python3
from re import finditer
query="your 2 kb DNA sequence"
smad3_binding_site="GTCTAGAC"
for matches in finditer(smad3_binding_site,query):
    print(matches.span(), matches.group())

hope it helps.

• 0 views • link
written by Prosad • Prosad

1 answer

JASPAR allows exporting motifs in MEME format which you can use along with FIMO from the MEME suite to look for specific motif occurrences in your sequence(s). As opposed to exact string matching it takes into consideration the probability of each base occurring per-position since chromatin binding proteins can be promiscuous with their binding sites.

• 0 views • link
written by rpolicastro rpolicastro

Log in to answer this question.

More posts like this