<html><head><meta http-equiv="Content-Type" content="text/html; charset=utf-8"></head><body style="word-wrap: break-word; -webkit-nbsp-mode: space; line-break: after-white-space;" class="">Hi Sadeem,<div class=""><br class=""></div><div class=""> As Elaine points out the trick is to get the helix axis in a form that the sym command can use. Unfortunately the sym command does not know about axes made with "define axis". But here is a pretty easy way to get the axis. I open 2 copies of the RNA structure, align the second copy RNA to the first copy shifted by one nucleotide, then use the axis of rotation defined by that alignment:</div><div class=""><br class=""></div><div class=""><span class="Apple-tab-span" style="white-space:pre"> </span>open 6g19</div><div class=""><span class="Apple-tab-span" style="white-space:pre"> </span>open 6g19</div><div class=""><span class="Apple-tab-span" style="white-space:pre"> </span>match #0:2-14.X,1-13.Y@P #1:1-13.X,2-14.Y@P </div><div class=""><span class="Apple-tab-span" style="white-space:pre"> </span>measure rotation #0 #1 </div><div class=""><span class="Apple-tab-span" style="white-space:pre"> </span>sym #0 group H,46,74,3 axis #2 center #2 </div><div class=""><br class=""></div><div class="">Below is what the result looks like.</div><div class=""><br class=""></div><div class=""><span class="Apple-tab-span" style="white-space:pre"> </span>Tom</div><div class=""><br class=""></div><div class=""><img apple-inline="yes" id="C1D83AC6-F71D-43C5-9389-1CC26D76F3BA" width="628" height="620" src="cid:AEA61E44-BDE8-42F9-B14E-76501765D94A" class=""><br class=""><div><br class=""><blockquote type="cite" class=""><div class="">On Sep 10, 2020, at 1:31 PM, Elaine Meng <<a href="mailto:meng@cgl.ucsf.edu" class="">meng@cgl.ucsf.edu</a>> wrote:</div><br class="Apple-interchange-newline"><div class=""><div class="">Hi Sadeem,<br class="">If the axis defined from the RNA is what you really need, here are some thoughts on a few steps of how to convert it to something that could be used by "sym."<br class=""><br class="">If there are two atoms that fall very close to the axis, you can just specify those as the axis value in the "sym" command. You would probably want to make the axis with a small radius in the define command and show the nucleic acid as atoms so that you can evaluate this, e.g.<br class=""><br class="">define axis radius 0.2 #0:.X,.Y<br class="">~nuc<br class="">~ribbon nucleic acid<br class=""><br class="">Then specify those two atoms as the "sym" axis value, e.g. select them with ctrl-click, Shift-ctrl click, and then specify the in the command as "sel", e.g.<br class=""><br class="">sym #0 group H,46,74,3 axis sel<br class=""><br class="">I tried that on some other structure with a RNA double helix and the result didn't look very good. It looks somewhat better if I also specify a center value, e.g. maybe in your case something like<br class=""><br class="">sym #0 group H,46,74,3 axis sel center #0:.X,.Y<br class=""><br class="">However, I don't know if that gives the appropriate center point for your system. See the "sym" pages for possible ways to specify a center.<br class=""><br class="">A longer way is to add new atoms at the axis-object endpoints:<br class="">Use the "Save" button on the Axes/Planes/Centroids table to save the axis information to a text file. Then calculate the endpoint coordinates of the axis according to the formulae given right next to the bottom of this page:<br class=""><<a href="http://www.rbvi.ucsf.edu/chimera/docs/ContributedSoftware/structuremeas/structuremeas.html#axes" class="">http://www.rbvi.ucsf.edu/chimera/docs/ContributedSoftware/structuremeas/structuremeas.html#axes</a>><br class=""><br class="">Then use the Build Structure tool, Start Structure to add one atom at the coordinates of one end of the axis and then again to add another atom at the coordinates of the other end of the axis.<br class=""><br class=""><<a href="http://www.rbvi.ucsf.edu/chimera/docs/ContributedSoftware/editing/editing.html" class="">http://www.rbvi.ucsf.edu/chimera/docs/ContributedSoftware/editing/editing.html</a>><br class=""><<a href="http://www.rbvi.ucsf.edu/chimera/docs/ContributedSoftware/editing/editing.html#start" class="">http://www.rbvi.ucsf.edu/chimera/docs/ContributedSoftware/editing/editing.html#start</a>><br class=""><br class="">Then in the "sym" command you can specify the two atoms you just added as the axis value. Maybe the axis center would be an appropriate point for the "sym" center value.<br class=""><br class="">I hope this helps,<br class="">Elaine<br class=""><br class=""><blockquote type="cite" class="">On Sep 10, 2020, at 12:23 PM, Elaine Meng <<a href="mailto:meng@cgl.ucsf.edu" class="">meng@cgl.ucsf.edu</a>> wrote:<br class=""><br class="">Hi Sadeem,<br class="">Axis objects from the "define" command cannot be used as axis specifications for the sym command.<br class=""><br class="">See the "sym" command documentation for how to specify axes:<br class=""><<a href="http://www.rbvi.ucsf.edu/chimera/docs/UsersGuide/midas/sym.html" class="">http://www.rbvi.ucsf.edu/chimera/docs/UsersGuide/midas/sym.html</a>><br class=""><br class="">"Specify axis of symmetry (default z), where the axis value can be:<br class=""><span class="Apple-tab-span" style="white-space:pre"> </span>• x - X-axis<br class=""><span class="Apple-tab-span" style="white-space:pre"> </span>• y - Y-axis<br class=""><span class="Apple-tab-span" style="white-space:pre"> </span>• z - Z-axis<br class=""><span class="Apple-tab-span" style="white-space:pre"> </span>• x,y,z (three values separated by commas only) - an arbitrary vector in the reference coordinate system <br class=""><span class="Apple-tab-span" style="white-space:pre"> </span>• an atom-spec of exactly two atoms (not necessarily bonded or in the same model) or one bond. A bond can only be specified by selecting it and using the word selected, sel, or picked; any atoms also selected at the time will be ignored."<br class=""><br class="">I hope this helps,<br class="">Elaine<br class="">-----<br class="">Elaine C. Meng, Ph.D. <br class="">UCSF Chimera(X) team<br class="">Department of Pharmaceutical Chemistry<br class="">University of California, San Francisco<br class=""><br class=""><br class=""><blockquote type="cite" class="">On Sep 10, 2020, at 12:18 PM, sadeem ahmad <<a href="mailto:sadeemahmad@gmail.com" class="">sadeemahmad@gmail.com</a>> wrote:<br class=""><br class="">Hi all,<br class="">I am trying to generate a helical trimer from the monomer pdb around a central double-stranded RNA. The rise and twist are 46 A and 74.0 deg respectively. I opened my monomer pdb in chimera (Screenshot-0) and used the following command to generate an axis around the central RNA double helix:<br class=""><br class="">define axis #0:.X,.Y<br class=""><br class="">where X and Y are the chain IDs for the 2 RNA strands<br class=""><br class="">This command helped me generate an axis named a1 along the RNA duplex (Screenshot-1). Now I want to generate a helical trimer along the RNA duplex with rise=46A and twist=74deg. For this I used the command:<br class=""><br class="">sym #0 group H,46,74,3 axis a1<br class=""><br class="">However, it is not placing the symmetric copies around axis a1 but around a distant axis (Screenshot-2). I am not sure how to generate the helical trimer along the a1 axis. Any help in this regard would be appreciated. Thanks.<br class=""><br class="">-Sadeem<br class=""></blockquote><br class=""><br class="">_______________________________________________<br class="">Chimera-users mailing list: <a href="mailto:Chimera-users@cgl.ucsf.edu" class="">Chimera-users@cgl.ucsf.edu</a><br class="">Manage subscription: <a href="https://plato.cgl.ucsf.edu/mailman/listinfo/chimera-users" class="">https://plato.cgl.ucsf.edu/mailman/listinfo/chimera-users</a><br class=""><br class=""></blockquote><br class=""><br class="">_______________________________________________<br class="">Chimera-users mailing list: <a href="mailto:Chimera-users@cgl.ucsf.edu" class="">Chimera-users@cgl.ucsf.edu</a><br class="">Manage subscription: <a href="https://plato.cgl.ucsf.edu/mailman/listinfo/chimera-users" class="">https://plato.cgl.ucsf.edu/mailman/listinfo/chimera-users</a><br class=""><br class=""></div></div></blockquote></div><br class=""></div></body></html>