User talk:ElNando888/Blog/SHAPE?

From Eterna Wiki

< User talk:ElNando888‎ | Blog

Revision as of 18:53, 22 July 2013 by Omei (talk | contribs)

Latest comment: 22 July 2013 by ElNando888

== 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]

 

----

Sorry.  As phrased, that was a really dumb question. You stated very clearly that you were comparing two different sequences, but somehow I managed to forget that as I was writing my comment.

What I was actually wondering was how you had figured out that the configuration in the two 3D models you show represented the configuration the sequence took in your design.  I too have spent a lot of time looking at the 3D structure of tetraloop hairpins in PDB, and my general impression is that a 6-base sequence will take on quite different configurations (for reasons unknown), just as the SHAPE results for the same sequence in the same position in the same lab can vary widely.  So while I have hypothesized that certain 3D configurations correspond to certain SHAPE score patterns, I really have no way of confirming or denying this.  I was hoping you had come up with something to help do that. 

Have you gathered many instances of the GGUAAC and CGUAAG loops in lab to see how consistent the SHAPE patterns are?  If they are reasonably consistent, and the PDB structures are consistent for this sequence across distinct molecules (there are a lot of duplicatations of the same molecule in PDB ), then your explanation bears a lot of weight.