# Detwiler: Questions Behe, Polyphen, and Ratchets

**URL:** <https://discourse.peacefulscience.org/t/detwiler-questions-behe-polyphen-and-ratchets/5088>\
**Category:** Conversation\
**Tags:** Science\
**Created:** [March 5, 2019, 3:31pm UTC](https://discourse.peacefulscience.org/t/detwiler-questions-behe-polyphen-and-ratchets/5088 "2019-03-05T15:31:27Z")\
**Posts on this page:** 20\
**Page:** 3

<div class="post-metadata">

**Author:** ![Rumraket](https://sea2.discourse-cdn.com/flex016/user_avatar/discourse.peacefulscience.org/rumraket/32/9328_2.png) [@Rumraket](https://discourse.peacefulscience.org/u/Rumraket)\
**Post date:** [March 6, 2019, 12:41am UTC](https://discourse.peacefulscience.org/t/detwiler-questions-behe-polyphen-and-ratchets/5088/41 "2019-03-06T00:41:53Z")

</div>

> [@swamidass](#):
>
> Well let’s test it then with some proteins, missense mutations, and stability predictors. Are you game?

Sure.

---

<div class="post-metadata">

**Author:** ![swamidass](https://sea2.discourse-cdn.com/flex016/user_avatar/discourse.peacefulscience.org/swamidass/32/3_2.png) [@swamidass](https://discourse.peacefulscience.org/u/swamidass)\
**Post date:** [March 6, 2019, 12:43am UTC](https://discourse.peacefulscience.org/t/detwiler-questions-behe-polyphen-and-ratchets/5088/42 "2019-03-06T00:43:11Z")

</div>

Gotta get back to work. I’ll move this into a new thread. See if you can pull a good set of proteins from swissprot or something. Perhaps we limit to to cytosolic proteins, because the predictors work best on them.

---

<div class="post-metadata">

**Author:** ![Rumraket](https://sea2.discourse-cdn.com/flex016/user_avatar/discourse.peacefulscience.org/rumraket/32/9328_2.png) [@Rumraket](https://discourse.peacefulscience.org/u/Rumraket)\
**Post date:** [March 6, 2019, 12:45am UTC](https://discourse.peacefulscience.org/t/detwiler-questions-behe-polyphen-and-ratchets/5088/43 "2019-03-06T00:45:46Z")

</div>

Sure, but it’s sleepy time for me now so we’ll have to wait till tomorrow.

---

<div class="post-metadata">

**Author:** ![Mercer](https://avatars.discourse-cdn.com/v4/letter/m/e274bd/32.png) [@Mercer](https://discourse.peacefulscience.org/u/Mercer)\
**Post date:** [March 6, 2019, 12:51am UTC](https://discourse.peacefulscience.org/t/detwiler-questions-behe-polyphen-and-ratchets/5088/44 "2019-03-06T00:51:35Z")

</div>

> [@Rumraket](#):
>
> But apparently mutations towards higher stability are less likely on average compared to mutations towards lower stability. It is true that stability increasing mutations are not necessarily function improving.

It’s also true that mutations that increase stability of a purified protein in solution can easily decrease the stability of the complex in which most that protein functions.

So absent encyclopedic knowledge, 50/50 is the best estimate.

---

<div class="post-metadata">

**Author:** ![Mercer](https://avatars.discourse-cdn.com/v4/letter/m/e274bd/32.png) [@Mercer](https://discourse.peacefulscience.org/u/Mercer)\
**Post date:** [March 6, 2019, 1:00am UTC](https://discourse.peacefulscience.org/t/detwiler-questions-behe-polyphen-and-ratchets/5088/45 "2019-03-06T01:00:17Z")

</div>

> [@John\_Detwiler](#):
>
> I will never not look at data.

[http://www.jbc.org/content/290/11/7003/F2.expansion.html](http://www.jbc.org/content/290/11/7003/F2.expansion.html)

Note that all of these cardiomyopathy mutants are shifted, in other words MORE stable than wild-type (normal). Two cause hypocontractility, the rest cause hypercontractility; opposite functional effects. None are good mutations to carry.

---

<div class="post-metadata">

**Author:** ![Rumraket](https://sea2.discourse-cdn.com/flex016/user_avatar/discourse.peacefulscience.org/rumraket/32/9328_2.png) [@Rumraket](https://discourse.peacefulscience.org/u/Rumraket)\
**Post date:** [March 6, 2019, 10:36am UTC](https://discourse.peacefulscience.org/t/detwiler-questions-behe-polyphen-and-ratchets/5088/46 "2019-03-06T10:36:51Z")

</div>

> [@Mercer](#):
>
> It’s also true that mutations that increase stability of a purified protein in solution can easily decrease the stability of the complex in which most that protein functions.
> 
> So absent encyclopedic knowledge, 50/50 is the best estimate.

I can see how a mutation with one effect in one environment, can have another effect in another environment, that makes perfect sense to me. But it’s not obvious to me why your latter statement should be assumed.

Attempts have been made at estimating the universal distribution of the stability effects of mutations on proteins. From figure 2 of the paper linked earlier:

 ![Fig2_Touriki_Tawfik_2009](https://us1.discourse-cdn.com/flex016/uploads/peacefulscience/original/2X/b/bdc9daae9b68f767b9efe9f73fc33bbead833482.png)

They’re clearly shifted towards more destabilizing, than towards more stabilizing. The data comes from reference 31, which is this paper: [https://www.sciencedirect.com/science/article/pii/S0022283607004342](https://www.sciencedirect.com/science/article/pii/S0022283607004342)

I’d just like to add that I don’t really have any dog in this fight. The conclusions aren’t mine, I’m just relaying what I read in these articles. If their predictions are contradicted by better data, I have no compunction accepting that and moving on. Don’t shoot the messenger.

---

<div class="post-metadata">

**Author:** ![John\_Detwiler](https://sea2.discourse-cdn.com/flex016/user_avatar/discourse.peacefulscience.org/john_detwiler/32/3683_2.png) [@John\_Detwiler](https://discourse.peacefulscience.org/u/John_Detwiler)\
**Post date:** [March 6, 2019, 3:33pm UTC](https://discourse.peacefulscience.org/t/detwiler-questions-behe-polyphen-and-ratchets/5088/47 "2019-03-06T15:33:30Z")

</div>

> [@Mercer](#):
>
> > Polyphen does run the results against pdb for structure and even a database for secondary structure. This obviously returns results that will highly correlate with stability.
> 
> That doesn’t seem obvious to me at all. It’s also perfectly clear from my own work that more stability does not correlate with improved function. That’s your most untenable assumption.

I don’t care about “improved function”. I’m perfectly fine with improved function or different function, or less function or no function. I’m even fine with mutations that increase stability. All I’m saying is that if you continue to lose stability on average, you eventually won’t be stable. That’s not really a controversial claim, and does not warrant extreme extrapolation in the opposite direction. Can’t I hypothesize this without also having to believe that natural selection will turn every protein into a neutron star?

If you are scoring for similarity toward things that by definition already fold, then the closer you are to that, the closer you are to something that folds and has been so far sustainable than you are to everything else. To characterize “everything else”, I think studies show that most mutations decrease stability. I read at least a few that seemed to agree that about 30% of mutations actually prevent folding on average. This would seem to indicate an extreme sensitivity toward sequence (most proteins are not 3 amino acids long) and gives me no reason to believe a mutation might not be expected to decrease stability in a protein it doesn’t break.

---

<div class="post-metadata">

**Author:** ![T\_aquaticus](https://sea2.discourse-cdn.com/flex016/user_avatar/discourse.peacefulscience.org/t_aquaticus/32/1249_2.png) [@T\_aquaticus](https://discourse.peacefulscience.org/u/T_aquaticus)\
**Post date:** [March 6, 2019, 4:00pm UTC](https://discourse.peacefulscience.org/t/detwiler-questions-behe-polyphen-and-ratchets/5088/48 "2019-03-06T16:00:39Z")

</div>

> [@Rumraket](#):
>
> I don’t think this relationship between protein function and protein stability is one of _necessity_ , it is however one of generality. And it is generally the case that there is some zone between **too rigid** , and **too unstable** , where most natural proteins sit and function “best”.

I would agree, in general. The two main properties are the stability of the tertiary structure and the chemical properties of amino acids in the active site.

What I am not sure of is if Polyphen can differentiate between stability and change in tertiary structure. I could see a scenario where Polyphen labels something as damaging that is in fact beneficial. There is also the strong possibility of compensatory mutations preventing damage from other mutations.

---

<div class="post-metadata">

**Author:** ![T\_aquaticus](https://sea2.discourse-cdn.com/flex016/user_avatar/discourse.peacefulscience.org/t_aquaticus/32/1249_2.png) [@T\_aquaticus](https://discourse.peacefulscience.org/u/T_aquaticus)\
**Post date:** [March 6, 2019, 4:08pm UTC](https://discourse.peacefulscience.org/t/detwiler-questions-behe-polyphen-and-ratchets/5088/49 "2019-03-06T16:08:18Z")

</div>

> [@John\_Detwiler](#):
>
> If you are scoring for similarity toward things that by definition already fold, then the closer you are to that, the closer you are to something that folds and has been so far sustainable than you are to everything else. To characterize “everything else”, I think studies show that most mutations decrease stability.

In the most extreme general sense, this is probably true. Of course, the reality is complicated. Different parts of the protein are going to have different sensitivities to changes, and it also depends heavily on what the new amino acid is. Swapping out a leucine for an isoleucine may not have any affect, but swapping out a leucine for a proline could have a drastic effect.

---

<div class="post-metadata">

**Author:** ![Mercer](https://avatars.discourse-cdn.com/v4/letter/m/e274bd/32.png) [@Mercer](https://discourse.peacefulscience.org/u/Mercer)\
**Post date:** [March 6, 2019, 4:16pm UTC](https://discourse.peacefulscience.org/t/detwiler-questions-behe-polyphen-and-ratchets/5088/50 "2019-03-06T16:16:00Z")

</div>

> [@John\_Detwiler](#):
>
> I don’t care about “improved function”. I’m perfectly fine with improved function or different function, or less function or no function. I’m even fine with mutations that increase stability.

I don’t think you looked at the data.

> [@John\_Detwiler](#):
>
> All I’m saying is that if you continue to lose stability on average, you eventually won’t be stable.

There are two problems with that:

1. We’re not talking about **continuing** to lose stability.
2. Instability plays an essential role in protein function, particularly in regulation, so not being stable is not necessarily nonfunctional. The disease-causing conformation of prions (PrP) is incredibly stable. That’s not good.

> [@John\_Detwiler](#):
>
> If you are scoring for similarity toward things that by definition already fold,

I wouldn’t. Even random proteins fold. It’s what they do.

> [@John\_Detwiler](#):
>
> then the closer you are to that, the closer you are to something that folds and has been so far sustainable than you are to everything else.

Again, “something that folds” is a _non sequitur_. You didn’t look at the data, did you?

> [@John\_Detwiler](#):
>
> To characterize “everything else”, I think studies show that most mutations decrease stability.

Then how did I just show you 8 killer mutations that increased stability?

> [@John\_Detwiler](#):
>
> I read at least a few that seemed to agree that about 30% of mutations actually prevent folding on average.

I think that you’re not understanding the basis of folding.

> [@John\_Detwiler](#):
>
> This would seem to indicate an extreme sensitivity toward sequence (most proteins are not 3 amino acids long) and gives me no reason to believe a mutation might not be expected to decrease stability in a protein it doesn’t break.

That’s because you are discounting the importance of instability in regulation. Proteins don’t just sit around in solution to be scored by a computer program. To understand these principles, I would suggest that you do a deep dive into the sarcomere, because it’s one of our best-understood systems and the system is fine-tuned to sit on a razor’s edge (instability), to be turned on or off by minute changes in intracellular calcium. Even the most structurally stable one, tropomyosin, is “detuned” to help it function as part of that razor’s edge.

---

<div class="post-metadata">

**Author:** ![John\_Detwiler](https://sea2.discourse-cdn.com/flex016/user_avatar/discourse.peacefulscience.org/john_detwiler/32/3683_2.png) [@John\_Detwiler](https://discourse.peacefulscience.org/u/John_Detwiler)\
**Post date:** [March 6, 2019, 4:19pm UTC](https://discourse.peacefulscience.org/t/detwiler-questions-behe-polyphen-and-ratchets/5088/51 "2019-03-06T16:19:49Z")

</div>

> [@Rumraket](#):
>
> However, since it actually IS possible for a mutation to increase stability, or somehow enhance function, but the software doesn’t know how to determine whether the mutation will have this effect, it has been coded to just assume it will be damaging.

> [@Rumraket](#):
>
> Because positive, or purifying selection will actually _some times_ get rid of damaging mutations, and therefore only allow neutral or even function-improving ones. For that reason reliance on this software is a mistake, and that is why we actually have to look at what genes these mutations happen in, and what these genes actually do.

> [@Rumraket](#):
>
> In contrast though is the APOB gene, involved in fatty acid transport in blood.

> [@Rumraket](#):
>
> But this time, because of what we know about what the gene actually does, transporting fatty acid out of blood and into important tissues, it is more likely that the mutations actuallly improved the function of the gene, because polar bears eat so much fat.

> [@Rumraket](#):
>
> We simply can’t take it at face value when we are actually dealing with mutations evolving under strong selection in a new environment, as the polar bear’s have been.

I agree with all of this. I’m really interested in what you guys discover in the simulations you are going to run (by the way, it’s really encouraging to see people who want to experiment instead of just pontificate). I feel like Joshua’s hypothesis accurately represents the evolutionary perspective, that if sustainability has been achieved, then there should be roughly a 50/50 split, but I feel like you are more likely correct and the literature seems to support the idea that even conserved mutations tend to reduce stability in contrast to ancestral sequences. The question then becomes, does purifying selection have the strength to offset these reductions. I have backup hypotheses which are harder to falsify and so are weaker. One of those is that many proteins seem to be designed with a certain amount of binding affinity that tunes them to behave in beneficial ways, like fibrin or fibrinogen or thrombin or whatever the one is that falls apart seconds after being activated so blood only clots locally.

---

<div class="post-metadata">

**Author:** ![John\_Detwiler](https://sea2.discourse-cdn.com/flex016/user_avatar/discourse.peacefulscience.org/john_detwiler/32/3683_2.png) [@John\_Detwiler](https://discourse.peacefulscience.org/u/John_Detwiler)\
**Post date:** [March 6, 2019, 4:21pm UTC](https://discourse.peacefulscience.org/t/detwiler-questions-behe-polyphen-and-ratchets/5088/52 "2019-03-06T16:21:02Z")

</div>

crap I just realized that’s not what you said. You do believe conserved sequences offset that reduction.

---

<div class="post-metadata">

**Author:** ![Mercer](https://avatars.discourse-cdn.com/v4/letter/m/e274bd/32.png) [@Mercer](https://discourse.peacefulscience.org/u/Mercer)\
**Post date:** [March 6, 2019, 4:23pm UTC](https://discourse.peacefulscience.org/t/detwiler-questions-behe-polyphen-and-ratchets/5088/53 "2019-03-06T16:23:47Z")

</div>

> [@Rumraket](#):
>
> But it’s not obvious to me why your latter statement should be assumed.

It’s because proteins interact with themselves, other proteins, substrates, etc. All of these things are missed when scoring stability of a protein, either _in vitro_ or _in silico_, when it’s not interacting with anything.

> [@Rumraket](#):
>
> Attempts have been made at estimating the universal distribution of the stability effects of mutations on proteins. From figure 2 of the paper linked earlier:

But FoldX is only being run on the protein doing absolutely nothing in aqueous solution. What if wild-type protein becomes less stable when bound to substrate, but the mutant, while less stable in solution, becomes more stable but less functional when bound to substrate? What if the converse occurs? Can you think of any biochemical reason why all four combos wouldn’t occur?

---

<div class="post-metadata">

**Author:** ![T\_aquaticus](https://sea2.discourse-cdn.com/flex016/user_avatar/discourse.peacefulscience.org/t_aquaticus/32/1249_2.png) [@T\_aquaticus](https://discourse.peacefulscience.org/u/T_aquaticus)\
**Post date:** [March 6, 2019, 4:29pm UTC](https://discourse.peacefulscience.org/t/detwiler-questions-behe-polyphen-and-ratchets/5088/54 "2019-03-06T16:29:01Z")

</div>

> [@John\_Detwiler](#):
>
> I agree with all of this. I’m really interested in what you guys discover in the simulations you are going to run (by the way, it’s really encouraging to see people who want to experiment instead of just pontificate).

Protein structure predictions are notoriously difficult, so I’m not sure how far we can really get. I have done very little work in this area, but have worked with others who have. Good protein modeling often requires time on supercomputers, so the best we may be able to do is simple comparisons for alpha helices and beta sheets.

---

<div class="post-metadata">

**Author:** ![John\_Detwiler](https://sea2.discourse-cdn.com/flex016/user_avatar/discourse.peacefulscience.org/john_detwiler/32/3683_2.png) [@John\_Detwiler](https://discourse.peacefulscience.org/u/John_Detwiler)\
**Post date:** [March 6, 2019, 4:30pm UTC](https://discourse.peacefulscience.org/t/detwiler-questions-behe-polyphen-and-ratchets/5088/55 "2019-03-06T16:30:51Z")

</div>

One thing I noticed in the original paper and also the wikipedia writeup on ApoB, " Mice overexpressing mApoB have increased levels of LDL “bad cholesterol” and decreased levels of [HDL] ‘good cholesterol’. Mice containing only one functional copy of the mApoB gene show the opposite effect, being resistant to [hypercholesterolemia]." This would seem to indicate that more is going on here than a simple A takes care of B type analysis. It would almost seem to argue in the opposite direction from a gene dosing perspective. In other words, it may be that blunting the effects of the protein product actually improve the ability of polar bears to deal with the required increase in dietary fat. Of course that would seem to vindicate Behe’s rule, but doesn’t exactly address my concerns or the way they overlap with his in terms of stability.

So I’m having trouble keeping up with all the quotes and redundancy here. If I’m replying at the wrong levels and duplicating responses, I’m terribly sorry. This format is killing me.

---

<div class="post-metadata">

**Author:** ![Mercer](https://avatars.discourse-cdn.com/v4/letter/m/e274bd/32.png) [@Mercer](https://discourse.peacefulscience.org/u/Mercer)\
**Post date:** [March 6, 2019, 4:35pm UTC](https://discourse.peacefulscience.org/t/detwiler-questions-behe-polyphen-and-ratchets/5088/56 "2019-03-06T16:35:36Z")

</div>

> [@John\_Detwiler](#):
>
> This would seem to indicate that more is going on here than a simple A takes care of B type analysis.

Well, yes, starting with the important fact that the adjectives “good” and “bad” only have that meaning in the context of one’s susceptibility to atherosclerosis, not for all the rest of biology involving cholesterol.

---

<div class="post-metadata">

**Author:** ![swamidass](https://sea2.discourse-cdn.com/flex016/user_avatar/discourse.peacefulscience.org/swamidass/32/3_2.png) [@swamidass](https://discourse.peacefulscience.org/u/swamidass)\
**Post date:** [March 6, 2019, 4:42pm UTC](https://discourse.peacefulscience.org/t/detwiler-questions-behe-polyphen-and-ratchets/5088/57 "2019-03-06T16:42:09Z")

</div>

> [@John\_Detwiler](#):
>
> This would seem to indicate that more is going on here than a simple A takes care of B type analysis.

That is right, which is why the whole notion of degradation of genes is poorly defined. What you are saying with stability is different. It has better definition, but still suffers, as @Mercer has explained.

As for @Rumraket clarifies he means there is about. 65%-70% chance that a muation will reduce stability in a folded protein. I don’t object to this. That is essentially 50/50, and it is self limiting too. Unstable proteins are more likely to get a mutation that increases stability. I had thought he meant that mutations that increased stability were exceedingly rare (that isn’t true).

If you keep negative selection in mind, we expect that proteins will remain as stable as they need to in order to work. I’m not seeing how this is a problem for evolution.

---

<div class="post-metadata">

**Author:** ![Mercer](https://avatars.discourse-cdn.com/v4/letter/m/e274bd/32.png) [@Mercer](https://discourse.peacefulscience.org/u/Mercer)\
**Post date:** [March 6, 2019, 4:43pm UTC](https://discourse.peacefulscience.org/t/detwiler-questions-behe-polyphen-and-ratchets/5088/58 "2019-03-06T16:43:51Z")

</div>

> [@swamidass](#):
>
> If you keep negative selection in mind, we expect that proteins will remain as stable as they need to in order to work.

And they often need to be less stable.

---

<div class="post-metadata">

**Author:** ![Rumraket](https://sea2.discourse-cdn.com/flex016/user_avatar/discourse.peacefulscience.org/rumraket/32/9328_2.png) [@Rumraket](https://discourse.peacefulscience.org/u/Rumraket)\
**Post date:** [March 6, 2019, 4:45pm UTC](https://discourse.peacefulscience.org/t/detwiler-questions-behe-polyphen-and-ratchets/5088/59 "2019-03-06T16:45:07Z")

</div>

> [@Mercer](#):
>
> It’s because proteins interact with themselves, other proteins, substrates, etc. All of these things are missed when scoring stability of a protein, either _in vitro_ or _in silico_ , when it’s not interacting with anything.

Right, but why should that skewer the distribution of stability effects of mutations towards more destabilizing?

You appear to be suggesting that if we put a protein back in it’s natural environment in a cell, this will have the effect of reversing the stability effects of mutations more toward stabilizing, for destabilizing mutations, than it will have the effect of reversing stabilizing mutations towards more destabilizing. I see zero reason why this should be the case.

I would agree that putting any given protein in another context will have the effect of reversing SOME of the stability effects of certain mutations, but I don’t see why the cellular environment would push or bias this potential reversal effect of mutations more towards being more likely stabilizing compared to the distribution outside of this context.

> [@Mercer](#):
>
> But FoldX is only being run on the protein doing absolutely nothing in aqueous solution. What if wild-type protein becomes less stable when bound to substrate, but the mutant, while less stable in solution, becomes more stable but less functional when bound to substrate? What if the converse occurs? Can you think of any biochemical reason why all four combos wouldn’t occur?

No I agree all four combos can occur, but they can occur to both stabilizing and destablizing mutation.

Suppose mutations are ~30% more likely to be destabilizing than stabilizing in vitro. Now we move the protein to it’s “natural” environment where it evolved, inside a cell, and we predict that some of the mutations will have reversed effects. Some mutations that are destabilizing in vitro, are stabilizing in vivo, but why would it not also be the case that some mutations, in at least equal proportion, are stabilizing in vitro, but become destabilizing in vivo?

I don’t see why the intracellular environment should produce a counterbiasing effect on on the stability effects of mutationsmutations.

---

<div class="post-metadata">

**Author:** ![swamidass](https://sea2.discourse-cdn.com/flex016/user_avatar/discourse.peacefulscience.org/swamidass/32/3_2.png) [@swamidass](https://discourse.peacefulscience.org/u/swamidass)\
**Post date:** [March 6, 2019, 4:46pm UTC](https://discourse.peacefulscience.org/t/detwiler-questions-behe-polyphen-and-ratchets/5088/60 "2019-03-06T16:46:50Z")

</div>

> [@Rumraket](#):
>
> Suppose mutations are ~30% more likely to be destabilizing than stabilizing in vitro.

Then it lends no support to @John_Detwiler.

[Previous page](https://discourse.peacefulscience.org/t/detwiler-questions-behe-polyphen-and-ratchets/5088.md?page=2)

[Next page](https://discourse.peacefulscience.org/t/detwiler-questions-behe-polyphen-and-ratchets/5088.md?page=4)
