Idealize does not set ideality to zero

Started by Susume

bkoep Staff Lv 1

The "Idealize Secondary Structure" tool sets the torsion angles of the residue to that of a common helix or beta sheet.

The "Idealize Peptide Bonds" tool sets the bond lengths, bond angles, and bond torsions of the peptide bond between two residues. This is generally useful for correctly closing cut-points, which disrupts the peptide bond between residues.

The "Ideality" subscore of a residue reflects all bond lengths, angles, and torsions of a residue, and will is affected by each of these tools.

Susume Lv 1

Seeing similar behavior on puzzle 1208 in new devprev. Many residues have slightly negative ideality (-1 to -2.5), and the idealize tool does not change them at all. There are no prolines in my solution. Are these small negative scores considered tolerable, and is idealize supposed to ignore them? The overall ideality score for the 75-residue protein, after idealizing the whole protein, is -75. This is pretty bad for a small protein with no prolines that has already been idealized (from experience I would expect this number to be in the -30 to -50 range).

Shared solution as "1208 idealize does not change ideality".

bkoep Staff Lv 1

Thanks for sharing, Susume!

The Ideality score component has been up-weighted in recent design puzzles like 1208, so I'm not surprised this penalty is more severe than what you're used to. We've noticed some peculiar geometry in Foldit designs and are trying to figure out what the problem is. There could be an issue with the Idealize tool, and your shared solution should be helpful in tracking down the problem.

Susume Lv 1

Seeing this again (or still) in 1285, the No Rebuild puzzle. Shared solution "1285 Idealize does not set ideality to zero". This protein has already been idealized once and then wiggled on low wiggle power, and now has a total ideality of -88, which is pretty poor for a protein this size (even among recent design puzzles). If you select all and hit Idealize a second time, the total ideality goes to -44, not zero. I still think the Idealize tool should set ideality to zero.

There are no prolines. The protein has never had cutpoints, and never been rebuilt, tweaked, or local wiggled. I have used rama map drag, Ideal SS, shake, global wiggle (always on auto wiggle power, with and without bands, with and without portions frozen), remix (by hand and in script) and mutate (by hand and in script). Somehow that subset of tools was enough to introduce a total ideality score of -100 (after several scripts, before I idealized the first time), and now Idealize is not able to get it back to zero.

Susume Lv 1

I reset the puzzle (1285 No Rebuild), and the total ideality on the extended chain is -204. This is approximately -2.7 per segment, but all the segments have slightly different ideality scores. I was assuming the extended chain we were given had zero ideality, and that all the segments were the same, but not so.

I selected the whole protein, set all to sheet, and hit Ideal SS. Now total ideality is -205. The ideality scores are now uniform, but not zero.

I thought the reason for the No Rebuild puzzle was that we kept coming up with poor ideality. Starting us off with poor ideality seems counterproductive.

inkycatz Lv 1

We're thinking of a number of solutions here to clear up some of the ambiguity, and possibly updating the tooltip. Nothing firm yet though.

Susume Lv 1

Puzzle 1290 (No Rebuilding) starts off with terrible ideality too (-247). Idealize on the extended chain only brings it up to -245. Setting the whole thing to either sheet or helix and using Ideal SS (which I'm pretty sure used to fix ideality as well as creating the chosen structure) also only gets it to -245. In end game on one track I was able to get it up to -49 (not a stable pose though). If this subscore measures bond angles and bond lengths that are not supposed to change, why give us a starting position with lots of room for improvement in those angles and lengths?

bkoep Staff Lv 1

Sorry again for the delayed reply. I've chased down a few experts, and have done a little digging myself, but still haven't quite gotten to the bottom of it…

It looks like the issue stems from improper hydrogen torsions. (In addition to bond lengths and angles, the ideality scoring also reports on bad dihedral angles.)

There are a few cases in which hydrogens along the protein backbone can move slightly out of place. It's not a huge problem, but there are tools in Rosetta to address it, and we should be able to incorporate those tools into Foldit sometime soon.

However, I'm also seeing that particular hydrogens on ILE and VAL sidechains are systematically misplaced. As above, the difference is slight enough that it shouldn't affect any model's validity from a scientific perspective. Nevertheless, there seem to be some small differences between the "ideal" torsion values and the coordinates used by Rosetta to create an ILE residue from scratch. I'm still unsure why there is any difference to begin with (and someone might have had a good reason for making it so), but this issue is very unlikely to be resolved anytime soon.

I'd like to reiterate that this is no cause for concern among players. The Idealize tool was never meant to completely level the ideality score. An ideality score in the range of 1-3 points per residue is entirely reasonable and lies well within the range of geometries observed in high-resolution crystal structures.

Susume Lv 1

Thank you for the detailed response!

It does seem related to certain sidechains - I had not noticed that before. The extended chain at the start of 1290 has ideality -247 (all isoleucine). When I change all to serine, ideality jumps to -9. All leucine: ideality -58. All valine: ideality -232.