[chimerax-users] Applying symmetry operations to generate an oligomeric model
Tom Goddard
goddard at sonic.net
Mon Nov 7 13:11:25 PST 2022
Hi Rob,
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.
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.
Tom
> On Nov 6, 2022, at 12:58 PM, Robert Harkness <r.harkness at utoronto.ca> wrote:
>
> Hi Tom,
>
> Thank you for your detailed response!
>
> 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.
>
> 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.
>
> 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?
>
> Best,
>
> Rob
>
>> On Nov 5, 2022, at 7:55 PM, Tom Goddard <goddard at sonic.net <mailto:goddard at sonic.net>> wrote:
>>
>> You don't often get email from goddard at sonic.net <mailto:goddard at sonic.net>. Learn why this is important <https://aka.ms/LearnAboutSenderIdentification>
>> Hi Rob,
>>
>> 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,
>>
>> measure center #1 <cxcmd:measure center #1>Center of mass grid index for emdb 15560 = (279.52, 279.92, 253.25)
>> Center of mass xyz coordinates for emdb 15560 = (289.20, 235.85, 185.32)
>>
>>
>> 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.
>>
>> 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".
>>
>> Sorry for the delayed answer. I was at a workshop all last week and did not have time for email.
>>
>> Tom
>>
>>> On Nov 3, 2022, at 4:35 AM, Robert Harkness via ChimeraX-users <chimerax-users at cgl.ucsf.edu <mailto:chimerax-users at cgl.ucsf.edu>> wrote:
>>>
>>> Hello,
>>>
>>> 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
>>>
>>> sym #0 (first asymmetric unit) group t,z3 axis <axis from measure rotation> center <axis points from measure rotation> coordinateSystem #1 (the density map)
>>>
>>> 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!
>>>
>>> Best,
>>>
>>> Rob
>>>
>>> _______________________________________________
>>> ChimeraX-users mailing list
>>> ChimeraX-users at cgl.ucsf.edu <mailto:ChimeraX-users at cgl.ucsf.edu>
>>> Manage subscription:
>>> https://www.rbvi.ucsf.edu/mailman/listinfo/chimerax-users <https://www.rbvi.ucsf.edu/mailman/listinfo/chimerax-users>
>>>
>>
>
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://www.rbvi.ucsf.edu/pipermail/chimerax-users/attachments/20221107/71a9bf7c/attachment.html>
More information about the ChimeraX-users
mailing list