<html><head><meta http-equiv="Content-Type" content="text/html charset=utf-8"></head><body style="word-wrap: break-word; -webkit-nbsp-mode: space; -webkit-line-break: after-white-space;" class=""><div dir="auto" style="word-wrap: break-word; -webkit-nbsp-mode: space; -webkit-line-break: after-white-space;" class="">Hi Zongli,<div class=""><br class=""></div><div class=""> Assign icosahedral symmetry to the map using the BIOMT matrices in the fit PDB coordinates then use vop cover to extend the octant map using that symmetry:</div><div class=""><br class=""></div><div class=""><span class="Apple-tab-span" style="white-space:pre"> </span>open emdbID:5199</div><div class=""><span class="Apple-tab-span" style="white-space:pre"> </span>open 4v7q</div><div class=""><span class="Apple-tab-span" style="white-space:pre"> </span>volume #0 symmetry #1</div><div class=""><span class="Apple-tab-span" style="white-space:pre"> </span>vop cover #0 ibox -499,-499,-499,499,499,499 cellsize 2000,2000,2000</div><div class=""><br class=""></div><div class="">The octant map is size 500,500,500 and I assume the 0,0,0 grid point is at the center of symmetry so the vop cover command says to cover grid points -499 to 499. The vop cover command is intended to be used with crystallography so I also tell it the unit cell is bigger than the full map (2000 instead of 1000) so it does not use periodicity of the map. The vop command took about 45 minutes to execute — very slow. </div><div class=""><br class=""></div><div class=""> The result of this symmetric covering is not ideal because the edges of the octant map could introduce artifacts in the full map. Also it is interpolating the one octant at all symmetry equivalent points and averaging, and this trilinear interpolation is going to introduce a little bit of blurring. Unfortunately because of the symmetry you can’t simply put 8 copies of the octant map together to form a full virus. The x,y,z axes are 2-fold symmetry axes so the given octant can be mapped to 4 of the 8 octants with no grid interpolations, but the other four require other symmetry operators, not just 90 degree rotations. So you can’t reproduce the full virus map with complete fidelity. The only way to avoid the interpolation is to get a full map from the Harrison lab. It is unfortunate that only this inadequate octant map is what is in the EMDB. In the attached image I show the EMDB 5199 octant, and a neighbor octant which obviously is not just a rotated copy since the original octant has 9 spikes while the neighbor octant has only 6.</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="10A59312-FF66-41BF-9C88-89D288122AD1" src="cid:FE798923-04E2-46DF-B9F1-CDE4D99131FF@cgl.ucsf.edu" class=""></div><div class=""><br class=""><div><blockquote type="cite" class=""><div class="">On Mar 27, 2017, at 9:49 AM, 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 Zongli,<br class="">I believe you can use “vop cover” to expand the map if you can supply the needed symmetry information and box size. Please see this previous post:<br class=""><br class=""><<a href="http://plato.cgl.ucsf.edu/pipermail/chimera-users/2014-March/009654.html" class="">http://plato.cgl.ucsf.edu/pipermail/chimera-users/2014-March/009654.html</a>><br class=""><br class="">I hope this helps,<br class="">Elaine<br class="">----------<br class="">Elaine C. Meng, Ph.D. <br class="">UCSF Computer Graphics Lab (Chimera team) and Babbitt Lab<br class="">Department of Pharmaceutical Chemistry<br class="">University of California, San Francisco<br class=""><br class=""><blockquote type="cite" class="">On Mar 26, 2017, at 8:01 PM, Li, Zongli <<a href="mailto:Zongli_Li@hms.harvard.edu" class="">Zongli_Li@hms.harvard.edu</a>> wrote:<br class=""><br class="">Hello,<br class="">I am wonder if there is a way in Chimera to create an entire icosahedral virus map from an octant of the virus density map?<br class="">Thanks,<br class="">Zongli<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="http://plato.cgl.ucsf.edu/mailman/listinfo/chimera-users" class="">http://plato.cgl.ucsf.edu/mailman/listinfo/chimera-users</a><br class=""><br class=""></div></div></blockquote></div><br class=""></div></div></body></html>