<p>I haven't worked on TFBSs but I really, really doubt you will be able to do this. In mutations, like in real estate, it's location, location, location. It isn't a question of how many mutations but what mutations and where. TFBSs are usually described using position weight matrices, precisely because that's how you can show the importance of different changes in different positions. </p>

<p>Here are a couple of random Google image search results for "transcription factor binding sites logo":</p>

<p><a href="https://i.sstatic.net/azBCo.png" rel="nofollow noreferrer"><img src="https://i.sstatic.net/azBCo.png" alt="CRP binding site logo"></a></p>

<p>(This was taken from <a href="http://compbio.pbworks.com/w/page/16252929/Transcription%20Regulation%20and%20Transcription%20Factor%20Motif%20Finding" rel="nofollow noreferrer">here</a> and you might want to read through it.)</p>

<p><a href="https://i.sstatic.net/1lDSd.jpg" rel="nofollow noreferrer"><img src="https://i.sstatic.net/1lDSd.jpg" alt="enter image description here"></a>
(source: Front. Genet., 23 February 2016 | <a href="https://doi.org/10.3389/fgene.2016.00024" rel="nofollow noreferrer">https://doi.org/10.3389/fgene.2016.00024</a>)</p>

<p>The main take home messages from these images are:</p>

<ol>
<li><p>Different positions have different importance and potential to disrupt binding. In the first image, I would expect that changing the T at position -3 to anything else would be pretty bad. However, changing an A to a T at -8 would likely be neutral. </p></li>
<li><p>Different TFs have very different binding sites and very different binding site lengths.  Look at panel B of the second image. That one appears to be a very conservative site and likely any change would be deleterious. Conversely, changing any of the nucleotides of the 17bp spacer shown in panel A would likely not be important as long as you don't change the length. </p></li>
</ol>

<p>Please note that I have not done more than glance at the sources I took these from, so my inferences may well be wrong, I only use them to illustrate the general point: you cannot really put a number to this. You will always need to look at each specific change and each specific position. In some cases, a single change can disrupt binding while in others multiple changes might have no effect at all. </p>

<hr>

<p>You might want to look into TFBS databases like JASPAR or TRANSFAC.</p>

CC BY-SA 4.0
2018-07-13T16:54:46.693Z
2018-07-13T16:54:46.693Z
terdon
<p>Every motif has a score. The score is calculated from the number of times each base was observed in that motif in the genome (I'm not sure how exactly these are calculated). Higher numbers indicate that this motif was found more frequently, which we take as a proxy for it providing better binding.</p>

<p>The <a href="http://www.ensembl.org/info/docs/tools/vep/index.html" rel="nofollow noreferrer">Ensembl VEP</a> can give you the score change for a mutated motif in the output, so this may be the easiest way to analyse this.</p>

CC BY-SA 4.0
2018-07-16T13:06:34.257Z
2018-07-16T13:06:34.257Z
Emily_Ensembl
<p>The motifbreakR package for R I believe does something like what you want to do;</p>

<p><a href="https://www.rdocumentation.org/packages/motifbreakR/versions/1.2.2/topics/motifbreakR" rel="nofollow noreferrer">https://www.rdocumentation.org/packages/motifbreakR/versions/1.2.2/topics/motifbreakR</a></p>

<p>If you're looking for a general rule as to what kind of mutations/where in the motif mutations will significantly affect motif binding, I suggest generating a lot of sequences with simulated SNP's (single nucleotide polymorphisms) and seeing how these score against the particular motif. Then, you will have a large(ish) database of sequences with one base pair difference from the motif and their corresponding score. This will allow you to judge what kinds of mutations affect the score the most. </p>

<p>Furthermore, you could somehow make a distribution of scores versus mutations and generate an empirical p-value to determine statistical significance (for instance, such and such mutation caused a score that was in the lowest 5% of scores for all mutations). You would know better than me how to do this, since your background is in CS and my PhD subject is genetics.</p>

<p>To make an assumption about the number of bases that need to be mutated, do this for one, two, three, etc. simulated SNP's.</p>

CC BY-SA 4.0
2020-05-10T21:08:09.610Z
2020-05-10T22:46:43.803Z
2020-05-10T22:46:43.803Z
SamChiappone
<p>How much does a single mutation/alteration of a nucleotide affect the presence of a transcription factor binding site (TFBS)? 
I am from computer science background(Obviously). 
I want to make a general assumption about the number of mutated bases in a TFBS to minimize the binding affinity?</p>

<fasta><sequence-analysis><phylogenetics>
Effect of mutation in DNA sequence on transcription factor binding sites
# Effect of mutation in DNA sequence on transcription factor binding sites

(Asked by: username_0 on 2018-07-13)

How much does a single mutation/alteration of a nucleotide affect the presence of a transcription factor binding site (TFBS)?
I am from computer science background(Obviously).
I want to make a general assumption about the number of mutated bases in a TFBS to minimize the binding affinity?

**Tags:** fasta, sequence-analysis, phylogenetics

---

## 3 answers

### Answer by username_1 on 2018-07-13. Score: 2

I haven't worked on TFBSs but I really, really doubt you will be able to do this. In mutations, like in real estate, it's location, location, location. It isn't a question of how many mutations but what mutations and where. TFBSs are usually described using position weight matrices, precisely because that's how you can show the importance of different changes in different positions.

Here are a couple of random Google image search results for "transcription factor binding sites logo":

[![CRP binding site logo](https://i.sstatic.net/azBCo.png)](https://i.sstatic.net/azBCo.png)

(This was taken from [here](http://compbio.pbworks.com/w/page/16252929/Transcription%20Regulation%20and%20Transcription%20Factor%20Motif%20Finding) and you might want to read through it.)

[![enter image description here](https://i.sstatic.net/1lDSd.jpg)](https://i.sstatic.net/1lDSd.jpg)
(source: Front. Genet., 23 February 2016 | <https://doi.org/10.3389/fgene.2016.00024>)

The main take home messages from these images are:

1. Different positions have different importance and potential to disrupt binding. In the first image, I would expect that changing the T at position -3 to anything else would be pretty bad. However, changing an A to a T at -8 would likely be neutral.
2. Different TFs have very different binding sites and very different binding site lengths. Look at panel B of the second image. That one appears to be a very conservative site and likely any change would be deleterious. Conversely, changing any of the nucleotides of the 17bp spacer shown in panel A would likely not be important as long as you don't change the length.

Please note that I have not done more than glance at the sources I took these from, so my inferences may well be wrong, I only use them to illustrate the general point: you cannot really put a number to this. You will always need to look at each specific change and each specific position. In some cases, a single change can disrupt binding while in others multiple changes might have no effect at all.

---

You might want to look into TFBS databases like JASPAR or TRANSFAC.


---

### Answer by username_2 on 2018-07-16. Score: 1

Every motif has a score. The score is calculated from the number of times each base was observed in that motif in the genome (I'm not sure how exactly these are calculated). Higher numbers indicate that this motif was found more frequently, which we take as a proxy for it providing better binding.

The [Ensembl VEP](http://www.ensembl.org/info/docs/tools/vep/index.html) can give you the score change for a mutated motif in the output, so this may be the easiest way to analyse this.


---

### Answer by username_3 on 2020-05-10. Score: 0

The motifbreakR package for R I believe does something like what you want to do;

<https://www.rdocumentation.org/packages/motifbreakR/versions/1.2.2/topics/motifbreakR>

If you're looking for a general rule as to what kind of mutations/where in the motif mutations will significantly affect motif binding, I suggest generating a lot of sequences with simulated SNP's (single nucleotide polymorphisms) and seeing how these score against the particular motif. Then, you will have a large(ish) database of sequences with one base pair difference from the motif and their corresponding score. This will allow you to judge what kinds of mutations affect the score the most.

Furthermore, you could somehow make a distribution of scores versus mutations and generate an empirical p-value to determine statistical significance (for instance, such and such mutation caused a score that was in the lowest 5% of scores for all mutations). You would know better than me how to do this, since your background is in CS and my PhD subject is genetics.

To make an assumption about the number of bases that need to be mutated, do this for one, two, three, etc. simulated SNP's.


---

