<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 Rob,<div class=""><br class=""></div><div class=""> Hopefully you are using ChimeraX, not Chimera -- we have only done maintenance on Chimera for the past 6 years. The ChimeraX commands all work in physical coordinates (Angstroms), not grid index units (pixels), and usually they use scene coordinates, not local model coordinates which makes a difference if you have moved individual models relative to one another in the scene.</div><div class=""><br class=""></div><div class=""> Your Dn symmetry problems are because the two-fold symmetry that flips one Cn ring to the other stacked Cn ring needs to be aligned on the x-axis. In your map apparently it is not aligned with the x-axis. In general it is helpful to have your map aligned in a standard way if you will be applying symmetry. For dihedral symmetry the z-axis is the cyclic symmetry axis, the origin should be on that cyclic symmetry axis mid-way between the two rings, and the x-axis should be a two-fold symmetry axis. If you don't use standard map orientations for symmetric maps it is going to be very painful to make symmetric copies of atomic models in those maps.</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><br class=""><blockquote type="cite" class=""><div class="">On Nov 6, 2022, at 12:58 PM, Robert Harkness <<a href="mailto:r.harkness@utoronto.ca" class="">r.harkness@utoronto.ca</a>> wrote:</div><br class="Apple-interchange-newline"><div class="">
<meta http-equiv="Content-Type" content="text/html; charset=utf-8" class="">
<div style="word-wrap: break-word; -webkit-nbsp-mode: space; line-break: after-white-space;" class="">
Hi Tom,
<div class=""><br class="">
</div>
<div class="">Thank you for your detailed response!</div>
<div class=""><br class="">
</div>
<div class="">I think I am using an older version of chimera. When I use the center <map> command, I just get one line of output that corresponds to the center of the map in grid coordinates. I have been converting these into angstrom coordinates using the
pixel size and applying those in sym.</div>
<div class=""><br class="">
</div>
<div class="">Sometimes when I open a model and a map, I have to rotate and translate the model into a reasonable position inside the map and then do a “fit in map” to get the initial placement for sym, so this means that the model’s coordinates will now be
different relative to the scene and the map. If I save the fitted model to a new PDB and then re-open it, will it now have the same coordinates of the scene (and map if I haven’t moved the map at all)? I think that way I can use scene coordinates for all maps
and models.</div>
<div class=""><br class="">
</div>
<div class="">I am finding for other cages I am working on that have Dn symmetry, in some cases they work beautifully if I just fit the model into the portion of the map that corresponds to one unit of the cage, then save that model as a PDB relative to the
map, then open it and use sym <model> group <Dn> center <map center in angstrom units> coordinateSystem <map>. In other cases doing this same procedure, the first half of the cage units are made perfectly as sym goes from the initial model around the Cn axis
placing copies in the density, but when it flips over the 2-fold axis to do the other half, the model units are out of alignment with the density. So I have been adding the axis option to get it closer but find it requires a little tweaking. I guess the axis
command doesn’t always give the perfect axis of rotation because the map is not perfect. Will this also be alleviated by just using everything in scene coordinates?</div>
<div class=""><br class="">
</div>
<div class="">Best,</div>
<div class=""><br class="">
</div>
<div class="">Rob<br class="">
<div class=""><br class="">
<blockquote type="cite" class="">
<div class="">On Nov 5, 2022, at 7:55 PM, Tom Goddard <<a href="mailto:goddard@sonic.net" class="">goddard@sonic.net</a>> wrote:</div>
<br class="Apple-interchange-newline">
<div class="">
<div style="word-wrap: break-word; -webkit-nbsp-mode: space; line-break: after-white-space;" class="">
<table border="0" cellspacing="0" cellpadding="0" width="100%" style="cellpadding:0;border:0;cellspacing:0;display:table;width:100%;table-layout:fixed;border-collapse:seperate;float:none;" align="left" class="">
<tbody style="display:block;" class="">
<tr class="">
<td valign="middle" width="1px" bgcolor="#A6A6A6" cellpadding="7px 2px 7px 2px" style="padding:7px 2px 7px 2px;background-color:#A6A6A6;valign:middle" class="">
</td>
<td valign="middle" width="100%" bgcolor="#EAEAEA" cellpadding="7px 5px 7px 15px" style="width:100%;background-color:#EAEAEA;padding:7px 5px 7px 15px;font-family:wf_segoe-ui_normal,Segoe UI,Segoe WP,Tahoma,Arial, sans-serif;font-size:12px;font-weight:normal;color:#212121;text-align:left;word-wrap:break-word;" class="">
<div class="">You don't often get email from <a href="mailto:goddard@sonic.net" class="">
goddard@sonic.net</a>. <a href="https://aka.ms/LearnAboutSenderIdentification" class="">
Learn why this is important</a></div>
</td>
<td valign="middle" align="left" width="75px" bgcolor="#EAEAEA" cellpadding="7px 5px 7px 5px" style="width:75px;background-color:#EAEAEA;padding:7px 5px 7px 5px;font-family:wf_segoe-ui_normal,Segoe UI,Segoe WP,Tahoma,Arial, sans-serif;font-size:12px;font-weight:normal;color:#212121;text-align:left;word-wrap:break-word;align:left" class="">
</td>
</tr>
</tbody>
</table>
<div class="">Hi Rob,
<div class=""><br class="">
</div>
<div class=""> It sounds like you did everything right. Using "measure center #1" where model #1 is the map should give the right center of the tetrahedron. When you use that command it reports the center in map grid index units and also in scene coordinates,
for example,</div>
<div class=""><br class="">
</div>
<div class="">
<div class="cxcmd" style="font-weight: bold; margin-top: 0.5em; background-color: rgb(221, 221, 221); font-family: "Times New Roman"; font-size: inherit; font-variant-ligatures: normal; orphans: 2; widows: 2; text-decoration-thickness: initial;">
<div class="cxcmd_as_cmd" style="display: inline;"><a href="cxcmd:measure center #1" class="">measure center #1</a></div>
</div>
<span style="font-family: "Times New Roman"; font-size: inherit; font-variant-ligatures: normal; orphans: 2; widows: 2; text-decoration-thickness: initial;" class="">Center of mass grid index for emdb 15560 = (279.52, 279.92, 253.25)</span><br style="font-family: "Times New Roman"; font-size: inherit; font-variant-ligatures: normal; orphans: 2; widows: 2; text-decoration-thickness: initial;" class="">
<span style="font-family: "Times New Roman"; font-size: inherit; font-variant-ligatures: normal; orphans: 2; widows: 2; text-decoration-thickness: initial;" class="">Center of mass xyz coordinates for emdb 15560 = (289.20, 235.85, 185.32)</span></div>
<div class=""><br class="">
</div>
<div class=""><br class="">
<div class="">In the sym command were you use that center to copy your atomic model with tetrahedral symmetry you should use the scene coordinates (the second line of output) not the grid index coordinates. The "sym" command expects scene coordinates for the
"center" option and you should not use the "sym" command coordinateSystem option. The coordinateSystem option is used only if the center coordinates are given in the coordinate system of a specific model rather than the whole scene. The models start out
with the same coordinate system as the whole scene, but if you move a model in the scene (e.g. when you fit an atomic model into the map) then that model has different local coordinates than the scene.</div>
<div class=""><br class="">
</div>
<div class=""> I see another potential problem in your sym command. The "axis" option is supposed to be a 3-fold symmetry axis of the tetrahedron. You are giving it the symmetry axis of the asymmetric unit trimer which may or may not be a symmetry axis of
the tetrahedron. To determine the symmetry axis of the tetrahedron open two copies of the map and rotate one copy by hand using the "Move Model" mouse mode (which moves only the selected model) to roughly align two different trimeric subunits then optimize
the alignment with "fitmap #2 in #1", then use measure rotation #2 to #1 and it will report a symmetry axis of the tetrahedron. Another way to go about this is if you density map already has a standard orientation, such as tetrahedron 2-fold symmetry axes
along x,y,z then you don't need to specify the axis in the sym command, just use "t,222" in the sym command. Or if it has a symmetry axis along z use "t,z3" with no axis specification. You can probably see if it has a standard symmetry using "measure symmetry
#1".</div>
<div class=""><br class="">
</div>
<div class=""> Sorry for the delayed answer. I was at a workshop all last week and did not have time for email.</div>
<div class=""><br class="">
</div>
<div class=""><span class="Apple-tab-span" style="white-space:pre"></span>Tom</div>
<div class=""><br class="">
<blockquote type="cite" class="">
<div class="">On Nov 3, 2022, at 4:35 AM, Robert Harkness via ChimeraX-users <<a href="mailto:chimerax-users@cgl.ucsf.edu" class="">chimerax-users@cgl.ucsf.edu</a>> wrote:</div>
<br class="Apple-interchange-newline">
<div class="">
<div class="">Hello,<br class="">
<br class="">
I am trying to generate a model of an oligomeric protein cage that consists of 4 trimers and has tetrahedral symmetry. I am using the asymmetric unit model consisting of one monomer (#0) and a density map (#1) as a starting point. My strategy has been to take
the asymmetric unit and fit it into a portion of the map corresponding to one of the monomers within one of the trimer faces of the tetrahedron. Then I generate a copy of that asymmetric unit (#2) and fit it into the map for an adjacent protomer within the
same trimer face. This results in the two asymmetric unit models being rotated 120 degrees relative to one another and part of the same trimer. Then I use the measure rotation command to measure the rotation of the second asymmetric unit relative to the first,
giving me a rotation matrix and axis parameters. I then do<br class="">
<br class="">
sym #0 (first asymmetric unit) group t,z3 axis <axis from measure rotation> center <axis points from measure rotation> coordinateSystem #1 (the density map)<br class="">
<br class="">
And this generates the right looking trimers, so the 3-fold axis appears to be correct, but the center is off such that only one of the four trimers matches the density map and is placed at the top of the tetrahedron, while the other 3 trimers are off into
the void of the scene but placed at the other 3 vertices. I am not sure how to appropriately define the center. My feeling is that it should be the center of the map, taken from the measure center #1 command, but implementing those numbers in sym doesn’t work
either. Any thoughts on how these center coordinates are defined/should be measured, or if there is a better approach for doing this would be awesome!<br class="">
<br class="">
Best,<br class="">
<br class="">
Rob<br class="">
<br class="">
_______________________________________________<br class="">
ChimeraX-users mailing list<br class="">
<a href="mailto:ChimeraX-users@cgl.ucsf.edu" class="">ChimeraX-users@cgl.ucsf.edu</a><br class="">
Manage subscription:<br class="">
<a href="https://www.rbvi.ucsf.edu/mailman/listinfo/chimerax-users" class="">https://www.rbvi.ucsf.edu/mailman/listinfo/chimerax-users</a><br class="">
<br class="">
</div>
</div>
</blockquote>
</div>
<br class="">
</div>
</div>
</div>
</div>
</blockquote>
</div>
<br class="">
</div>
</div>
</div></blockquote></div><br class=""></div></body></html>