User talk:ElNando888/Blog/SHAPE?: Difference between revisions

From Eterna Wiki

< User talk:ElNando888‎ | Blog

Latest comment: 22 July 2013 by ElNando888

(Request for clarification)
(replying)
Line 17: Line 17:
<p>I'm not claiming the latter is the right explanation; I'm just trying to better understand your chain of reasoning.</p>
<p>I'm not claiming the latter is the right explanation; I'm just trying to better understand your chain of reasoning.</p>
<p>[[User:Omei|Omei]] ([[User talk:Omei|talk]]) 16:51, 22 July 2013 (UTC)</p>
<p>[[User:Omei|Omei]] ([[User talk:Omei|talk]]) 16:51, 22 July 2013 (UTC)</p>
<p>----</p>
<p>No, you could not switch them, because the sequences are clearly different. GGUAAC vs CGUAAG (closing pair swapped)</p>
<p>FRABASE returned in each case a few hits, all of them presenting a similar GNRA-like pattern (which is natural). I didn't verify if the ones I selected were more representative than the others, but I'm fairly sure that all GGUAAC looked very much alike, that all CGUAAG also looked alike, and that the main difference between the two sets was indeed to be found in the position of the first G in the loop relative to the bases forming the closing pair, in other words, their stacking relationships.</p>
<p>And using solved structures seems to me like the best possible option here. Would you rather trust a software simulation than whatever we can manage to measure with XRD or NMR?</p>
<p><em>-- [[User:ElNando888|ElNando888]] ([[User talk:ElNando888|talk]]) 17:45, 22 July 2013 (UTC)</em></p>
<p><em><br /></em></p>
<p><em><br /></em></p>

Revision as of 17:45, 22 July 2013

== General ==

Great topic! I agree with the basic premise that the SHAPE results are affected by numerous things besides the presence or absence of the three base pairings that the game's energy model acknowledges.  Working to figure out those additional factors is perhaps the most cool aspect of EteRNA.

Omei (talk) 20:34, 21 July 2013 (UTC)Reply[reply]

== Toolset ==

What tools did you use to create these images?  Something other than RNA Composer to predict the 3D structure?  And the renderings don't look like options I have seen in Chimera, though there could be plugins I don't know about.  I'm especially interested in what you used to predict stacking interactions.

Omei (talk) 20:34, 21 July 2013 (UTC)Reply[reply]

----

The 3D renderings are no simulations, they are all segments of solved structures from PDB. First, I used FRABASE to identify the PDB entries containing the sequences I was interested in. Then I fetched the structures in Chimera and worked them with it.

Atoms and "normal" bonds are easy to color any way you want, but you need to define your own color to get the the translucid white I applied on the backbone (select/structure/backbone/full). For the hydrogen bonds, I currently like to make them look like springs, but it's an option relatively hard to access: tools/general controls/pseudobond panel/hydrogen bonds/attributes/component pseudobond attribute/bond style

And you're right, there are no built-in command to detect and represent stacking interactions. First, and as I was saying, it took me a very long while (months) before I could find the parameters that validate a stacking interaction. I found them on the website of another 3D rendering tool, PyMol, at http://www.schrodinger.com/kb/1556

Then it's an arduous work in Chimera, but at least it's possible. There are tools in Structure Analysis which allow to define centroids, planes, axis, show or hide these objects, and also to calculate distances and plane angles. Once a pi stacking interaction is validated, I represent the interaction by defining an axis based on the atoms of both aromatic rings.

The other difficulty I had related to this topic was to ascertain whether the imidazole (pentagonal ring in purines) is aromatic or not. I don't recall exactly where I read that they actually are, but considering that reference and the various examples I saw in solved structures, I believe that they are.

-- ElNando888 (talk) 05:15, 22 July 2013 (UTC)Reply[reply]

----

The 3D renderings are no simulations, they are all segments of solved structures from PDB. First, I used FRABASE to identify the PDB entries containing the sequences I was interested in. Then I fetched the structures in Chimera and worked them with it.

Does that mean the association of the specific 3D structures with the 2 hairpins is simply your choice?  If so, could I equally well switch them, and attribute the lower SHAPE reactivity at 26 to the fact that the G has two hydrogen bonds with the opposing A?

I'm not claiming the latter is the right explanation; I'm just trying to better understand your chain of reasoning.

Omei (talk) 16:51, 22 July 2013 (UTC)Reply[reply]

----

No, you could not switch them, because the sequences are clearly different. GGUAAC vs CGUAAG (closing pair swapped)

FRABASE returned in each case a few hits, all of them presenting a similar GNRA-like pattern (which is natural). I didn't verify if the ones I selected were more representative than the others, but I'm fairly sure that all GGUAAC looked very much alike, that all CGUAAG also looked alike, and that the main difference between the two sets was indeed to be found in the position of the first G in the loop relative to the bases forming the closing pair, in other words, their stacking relationships.

And using solved structures seems to me like the best possible option here. Would you rather trust a software simulation than whatever we can manage to measure with XRD or NMR?

-- ElNando888 (talk) 17:45, 22 July 2013 (UTC)Reply[reply]