<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 Yeping,<div class=""><br class=""></div><div class=""> Chimera does not need high map resolution for fitting -- in fact high resolution (e.g. 4 Angstroms) makes it harder to find the correct orientation of the molecule because there can be many locally optimal fits and Chimera only finds locally optimal fits and has limited search capabilities. So for a high resolution map it is often helpful to first Gaussian smooth it (Tools / Volume Data / Volume Filter or equivalent "vop gaussian" command) to get a lower resolution blob, fit the atomic model in that, then switch to fitting into the high resolution map once the approximately correct orientation is found at lower resolution.</div><div class=""><br class=""></div><div class=""> The main difficult fitting problem you are likely to encounter is that Chimera is only fitting the molecule as a rigid body. So if the shape of your atomic model does not match the conformation seen in the map then you won't find any good fit. </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 Aug 20, 2020, at 6:49 PM, sunyeping <<a href="mailto:sunyeping@aliyun.com" class="">sunyeping@aliyun.com</a>> wrote:</div><br class="Apple-interchange-newline"><div class=""><div class="__aliyun_email_body_block"><div style="clear:both;" class="">Hi Greg,</div><div style="clear:both;" class=""><br class=""></div><div style="clear:both;" class="">Thank you for your help.</div><div style="clear:both;" class="">The atomic model of the Ab was built with Rosetta. Maybe one reason it cannot be automatically fitted into the map is that it does not has a correct H-L orientation.</div><div style="clear:both;" class="">But is high local resolution of the map a prerequisite for successful automatical fitting with chimera? </div><div style="clear:both;" class=""><br class=""></div><div style="clear:both;" class="">Best,</div><div style="clear:both;" class="">Yeping</div><blockquote style="margin-right:.0px;margin-top:.0px;margin-bottom:.0px;" class=""><div style="clear:both;" class=""><span style="font-family: Tahoma, Arial, STHeiti, SimSun; font-size: 14px;" class="">------------------------------------------------------------------</span></div><div style="clear:both;" class=""><span style="font-family: Tahoma, Arial, STHeiti, SimSun; font-size: 14px;" class="">From:Tom Goddard <<a href="mailto:goddard@sonic.net" class="">goddard@sonic.net</a>></span></div><div style="clear:both;" class=""><span style="font-family: Tahoma, Arial, STHeiti, SimSun; font-size: 14px;" class="">Sent At:2020 Aug. 21 (Fri.) 00:50</span></div><div style="clear:both;" class=""><span style="font-family: Tahoma, Arial, STHeiti, SimSun; font-size: 14px;" class="">To:孙业平 <<a href="mailto:sunyeping@aliyun.com" class="">sunyeping@aliyun.com</a>></span></div><div style="clear:both;" class=""><span style="font-family: Tahoma, Arial, STHeiti, SimSun; font-size: 14px;" class="">Cc:chimera-users <<a href="mailto:chimera-users@cgl.ucsf.edu" class="">chimera-users@cgl.ucsf.edu</a>></span></div><div style="clear:both;" class=""><span style="font-family: Tahoma, Arial, STHeiti, SimSun; font-size: 14px;" class="">Subject:Re: [Chimera-users] Cannot dock molecule to targeted region?</span></div><div style="clear:both;" class=""><span style="font-family: Tahoma, Arial, STHeiti, SimSun; font-size: 14px;" class=""><br class=""></span></div>Hi Yeping,<div class=""><br class=""></div><div class=""> You have a difficult problem and Chimera is not going to be able to automatically produce the right fit, and neither am I. The preprint on the map you are fitting EMDB 22158 is here</div><div class=""><br class=""></div><div class=""> <a class="" href="https://www.biorxiv.org/content/10.1101/2020.06.17.153486v2.full" target="_blank">https://www.biorxiv.org/content/10.1101/2020.06.17.153486v2.full</a></div><div class=""><br class=""></div><div class="">They do not provide a PDB for the 22158 spike + antibody, although they do provide PDB 6xey for the related EMDB 22156 which is at higher resolution. Why is that? Probably because they are not confident of their antibody fit to 22158 or could not even make a trustworthy model of that antibody given the low resolution. What is the resolution? If you look at figure 13e where they show the local resolution, the antibody region has map resolution in the range 9.5 to 11.5 A. That is definitely not good enough to build an antibody model, nor even good enough to see the beta sheets of the antibody. By comparison look at the higher resolution EMDB 22156 and PDB 6xey they determined shown in figure 12 -- map resolution of around 4-6 A near the antibody, and they show their antibody model and density in figure 12e and you can see it is pretty reasonable.</div><div class=""><br class=""></div><div class=""> In summary, I don't know if your antibody atomic model is in the right conformation and you are trying to fit it into low ~10 Angstrom density. Chimera is not able to exhaustively sample configurations and tell you the best and even if it did it could very well be wrong with these uncertainties. Chimera is primarily used to interactively do fitting allowing you to understand whether a reliable fit is even possible.</div><div class=""><br class=""></div><div class=""> Good luck!</div><div class=""><br class=""></div><div class=""> Tom</div><div class=""><br class=""></div><div class=""><br class=""></div><div class="">Fig 13 from paper.</div><div class=""><span id="cid:__aliyun159797454317359771"><F13.large.jpg></span><br class=""><div class="">Fig 12 from paper</div><div class=""><br class=""></div><div class=""><span id="cid:__aliyun159797454317359772"><F12.large.jpg></span><br class=""><div class="">On Aug 19, 2020, at 8:59 PM, sunyeping <<a class="" href="mailto:sunyeping@aliyun.com" target="_blank">sunyeping@aliyun.com</a>> wrote:</div><br class="Apple-interchange-newline"><div class=""><div class="" style="clear:both;">Hi Tom,</div><div class="" style="clear:both;"><br class=""></div><div class="" style="clear:both;">Thank you for your suggestions and instructions. I have tried your method. Yes, the <span class=" __aliyun_node_has_bgcolor" style="font-family:Tahoma,Arial,STHeiti,SimSun;font-size:14.0px;font-style:normal;font-variant-ligatures:normal;font-variant-caps:normal;font-weight:400;text-align:start;text-indent:.0px;text-transform:none;widows:2;background-color:#ffffff;float:none;display:inline;"> "Optimize correlation" option in the "Fit in Map" tool can avoid large shit of the molecule. However, the problem with this method is that the final fitted result is determined by the initial position of the antibody molecule. I have no idea how the antibody binds the antigen, and my expectation is to find the correct binding positon of the antibody to the antigen by using the fitting tools of chimera. The only thing I know is that the CDRs face of the antibody should be placed towards the antigen. I have put the coordinations for the antigen (Ag, S.pdb) and the antibody (Ab, 4-8_H3-0175.pdb) at the following link:</span></div><div class="" style="clear:both;"><span class=" __aliyun_node_has_bgcolor" style="font-family:Tahoma,Arial,STHeiti,SimSun;font-size:14.0px;font-style:normal;font-variant-ligatures:normal;font-variant-caps:normal;font-weight:400;text-align:start;text-indent:.0px;text-transform:none;widows:2;background-color:#ffffff;float:none;display:inline;"> <a class="" href="https://drive.google.com/drive/folders/16ELIDsUxiNwRGp9ely2SRn-zjAdHCudH?usp=sharing" target="_blank">https://drive.google.com/drive/folders/16ELIDsUxiNwRGp9ely2SRn-zjAdHCudH?usp=sharing</a></span></div><div class="" style="clear:both;"><br class=""></div><div class="" style="clear:both;">Firstly, I fitted the Ag to the map for the ag-ab complex (emd_22158.map), then I put the antibody to target area of the map, and then used Fit in Map tool to fit it and get the complex coordinate fit-in-map_1.pdb. After that, I had another try. I put the Ab at the same initial position in the target area but rotated it by 180 degree around its long axis, then fit it, and get <span class=" __aliyun_node_has_bgcolor" style="font-family:Tahoma,Arial,STHeiti,SimSun;font-size:14.0px;font-style:normal;font-variant-ligatures:normal;font-variant-caps:normal;font-weight:400;text-align:start;text-indent:.0px;text-transform:none;widows:2;background-color:#ffffff;float:none;display:inline;">complex coordinate </span><span class=" __aliyun_node_has_bgcolor" style="font-family:Tahoma,Arial,STHeiti,SimSun;font-size:14.0px;font-style:normal;font-variant-ligatures:normal;font-variant-caps:normal;font-weight:400;text-align:start;text-indent:.0px;text-transform:none;widows:2;background-color:#ffffff;float:none;display:inline;">fit-in-map_1.pdb. Both of the two </span> <span class=" __aliyun_node_has_bgcolor" style="font-family:Tahoma,Arial,STHeiti,SimSun;font-size:14.0px;font-style:normal;font-variant-ligatures:normal;font-variant-caps:normal;font-weight:400;text-align:start;text-indent:.0px;text-transform:none;widows:2;background-color:#ffffff;float:none;display:inline;">complex coordinate can also be found at the above link. In these two complex structures, the heavy chain (H) and light chain (L) of the Ab are totally opposite (a rotation of 180 degree). </span></div><div class="" style="clear:both;"><span class=" __aliyun_node_has_bgcolor" style="font-family:Tahoma,Arial,STHeiti,SimSun;font-size:14.0px;font-style:normal;font-variant-ligatures:normal;font-variant-caps:normal;font-weight:400;text-align:start;text-indent:.0px;text-transform:none;widows:2;background-color:#ffffff;float:none;display:inline;"> Maybe you could try to do the fitting and find the most reasonable binding mode of the Ab? One clue is that the CDR3 of the H chain (residue 93-107) should usually be involved in contact with the antigen.</span></div><div class="" style="clear:both;"><span class=" __aliyun_node_has_bgcolor" style="font-family:Tahoma,Arial,STHeiti,SimSun;font-size:14.0px;font-style:normal;font-variant-ligatures:normal;font-variant-caps:normal;font-weight:400;text-align:start;text-indent:.0px;text-transform:none;widows:2;background-color:#ffffff;float:none;display:inline;"><br class=""></span></div><div class="" style="clear:both;"><span class=" __aliyun_node_has_bgcolor" style="font-family:Tahoma,Arial,STHeiti,SimSun;font-size:14.0px;font-style:normal;font-variant-ligatures:normal;font-variant-caps:normal;font-weight:400;text-align:start;text-indent:.0px;text-transform:none;widows:2;background-color:#ffffff;float:none;display:inline;">Best regards,</span></div><div class="" style="clear:both;"><span class=" __aliyun_node_has_bgcolor" style="font-family:Tahoma,Arial,STHeiti,SimSun;font-size:14.0px;font-style:normal;font-variant-ligatures:normal;font-variant-caps:normal;font-weight:400;text-align:start;text-indent:.0px;text-transform:none;widows:2;background-color:#ffffff;float:none;display:inline;">Yeping</span></div><div class="" style="clear:both;"><span class="" style="font-family:Tahoma,Arial,STHeiti,SimSun;font-size:14.0px;">------------------------------------------------------------------</span></div><div class="" style="clear:both;"><span class="" style="font-family:Tahoma,Arial,STHeiti,SimSun;font-size:14.0px;">From:Tom Goddard <<a class="" href="mailto:goddard@sonic.net" target="_blank">goddard@sonic.net</a>></span></div><div class="" style="clear:both;"><span class="" style="font-family:Tahoma,Arial,STHeiti,SimSun;font-size:14.0px;">Sent At:2020 Aug. 20 (Thu.) 01:15</span></div><div class="" style="clear:both;"><span class="" style="font-family:Tahoma,Arial,STHeiti,SimSun;font-size:14.0px;">To:孙业平 <<a class="" href="mailto:sunyeping@aliyun.com" target="_blank">sunyeping@aliyun.com</a>></span></div><div class="" style="clear:both;"><span class="" style="font-family:Tahoma,Arial,STHeiti,SimSun;font-size:14.0px;">Cc:chimera-users <<a class="" href="mailto:chimera-users@cgl.ucsf.edu" target="_blank">chimera-users@cgl.ucsf.edu</a>>; Greg Pintilie <<a class="" href="mailto:gregdp@gmail.com" target="_blank">gregdp@gmail.com</a>></span></div><div class="" style="clear:both;"><span class="" style="font-family:Tahoma,Arial,STHeiti,SimSun;font-size:14.0px;">Subject:Re: [Chimera-users] Cannot dock molecule to targeted region?</span></div><div class="" style="clear:both;"><span class="" style="font-family:Tahoma,Arial,STHeiti,SimSun;font-size:14.0px;"><br class=""></span></div>Hi Yeping,<div class=""><br class=""></div><div class=""> As Greg says Chimera is trying to move your antibody into higher density. Besides the masking he suggests, another thing you can do is just use the Fit in Map tool (not Fit to Segments), place your anti-body by hand in approximately the correct position using the "Active" checkbuttons in model panel and the mouse, and under options in the Fit in Map tool choose "Optimize correlation" instead of "Optimize overlap". The default optimize overlap means to move the molecule into the highest density, but optimizing correlation instead moves it to get the best agreement between the shape of the density. That option will be grayed out until you check "Use map simulated from atoms" since correlation needs a simulated map, and you should enter the approximate resolution of your map. I don't think Fit to Segments offers that option to optimize correlation. When you use Fit in Map it just does a local optimization so your initial placement and orientation of the antibody has to be approximately right.</div><div class=""><br class=""></div><div class=""> Tom</div><div class=""><br class=""></div><div class=""><span id="cid:__aliyun159789594465151061" class=""><fit_in_map.png></span></div><div class=""><div class=""><br class=""><div class="">On Aug 19, 2020, at 5:13 AM, Greg Pintilie <<a class="" href="mailto:gregdp@gmail.com" target="_blank">gregdp@gmail.com</a>> wrote:</div><br class="Apple-interchange-newline"><div class=""><div class="" style="overflow-wrap:break-word;"><div class=""><br class=""></div>Hi Yeping,<div class=""><br class=""></div><div class="">It seems that the part of the map you are fitting to has lower values, and the Fit-To-Segments procedure is moving the model (your antibody) towards higher values. To prevent this, there is an option in the Fit To Segments dialog ‘Mask map with region(s) to prevent large drifts’. If you enable this option you should not see the model move outside the region. You may have to press the Options button at the bottom of the dialog first to reveal this option if you don’t see it. Also make sure to try ‘Rotational search’ if the default/faster ‘Align principal axes’ doesn’t seem to find the right fit.</div><div class=""><br class=""></div><div class="">There is also the ‘fitmap’ command in Chimera which does a more global search, perhaps the Chimera team may have some advice on it in this scenario, but I think you have the right idea with segmenting first and then doing the fit - it’s a bit more work but you can target the search better with it.</div><div class=""><a class="" href="https://www.cgl.ucsf.edu/chimera/current/docs/UsersGuide/midas/fitmap.html" target="_blank">https://www.cgl.ucsf.edu/chimera/current/docs/UsersGuide/midas/fitmap.html</a></div><div class=""><br class=""></div><div class="">Greg</div><div class=""><br class=""><div class=""><br class=""><div class="">On Aug 19, 2020, at 1:53 AM, 孙业平 <<a class="" href="mailto:sunyeping@aliyun.com" target="_blank">sunyeping@aliyun.com</a>> wrote:</div><br class="Apple-interchange-newline"><div class=""><div class="" style="clear:both;"><span class=" __aliyun_node_has_color" style="font-family:Tahoma,Arial,STHeiti,SimSun;font-size:14.0px;">Dear all,</span></div><div class="" style="clear:both;"><span class=" __aliyun_node_has_color" style="font-family:Tahoma,Arial,STHeiti,SimSun;font-size:14.0px;"><br class=""></span></div><div class="" style="clear:both;"><span class=" __aliyun_node_has_color" style="font-family:Tahoma,Arial,STHeiti,SimSun;font-size:14.0px;"> I am trying to fit an antibody molecule to the cryo-EM density map for an antigen-antibody complex. The related files can be found at the following link</span></div><div class="" style="clear:both;"><span class=" __aliyun_node_has_color" style="font-family:Tahoma,Arial,STHeiti,SimSun;font-size:14.0px;"> <a class="" href="https://drive.google.com/drive/folders/16ELIDsUxiNwRGp9ely2SRn-zjAdHCudH?usp=sharing" target="_blank">https://drive.google.com/drive/folders/16ELIDsUxiNwRGp9ely2SRn-zjAdHCudH?usp=sharing</a>.</span></div><div class="" style="clear:both;"><span class=" __aliyun_node_has_color" style="font-family:Tahoma,Arial,STHeiti,SimSun;font-size:14.0px;">I first segmented the map ("emd_22158.map") and get the segment map ("emd_22158.seg") in chimera. Then I selected the regions corresponding to the antibody and use the "Fit to segments" tool of chimera to fit the antibody molecule. The "image1.png" shows the inital antibody postion and the expected docking postion in the segmented complex density (indicated by the red arrow). However, the antibody molecule cannot be fitted into the selected target region after this operation, but moved to the position shown in "image2.png (only the target regions are shown). </span></div><div class="" style="clear:both;"><span class=" __aliyun_node_has_color" style="font-family:Tahoma,Arial,STHeiti,SimSun;font-size:14.0px;"> I have tried to ungoup and regroup the target region to make its overall shape more simial to the antibody, and place the antibody into the region and ajust its orienation, but whatever I tried, the results are the same: the antibody will move to the position shown in "image2.png" after performing the <span class=" __aliyun_node_has_bgcolor" style="font-family:Tahoma,Arial,STHeiti,SimSun;font-size:14.0px;font-style:normal;font-variant-ligatures:normal;font-variant-caps:normal;font-weight:400;text-align:start;text-indent:.0px;text-transform:none;widows:2;background-color:#ffffff;float:none;display:inline;">"Fit to segments" action</span>.</span></div><div class="" style="clear:both;"><span class=" __aliyun_node_has_color" style="font-family:Tahoma,Arial,STHeiti,SimSun;font-size:14.0px;"> I wonder whether it is impossible to dock the antibody to the expected position with chimera or there is anything wrong with my operation. Could you help me with this issue? I will appreciate any help from you.</span></div><div class="" style="clear:both;"><span class=" __aliyun_node_has_color" style="font-family:Tahoma,Arial,STHeiti,SimSun;font-size:14.0px;"><br class=""></span></div><div class="" style="clear:both;"><span class=" __aliyun_node_has_color" style="font-family:Tahoma,Arial,STHeiti,SimSun;font-size:14.0px;">Best regards,</span></div><div class="" style="clear:both;"><span class=" __aliyun_node_has_color" style="font-family:Tahoma,Arial,STHeiti,SimSun;font-size:14.0px;">Yeping Sun</span></div>_______________________________________________<br class="">Chimera-users mailing list: <a class="" href="mailto:Chimera-users@cgl.ucsf.edu" target="_blank">Chimera-users@cgl.ucsf.edu</a><br class="">Manage subscription: <a class="" href="https://plato.cgl.ucsf.edu/mailman/listinfo/chimera-users" target="_blank">https://plato.cgl.ucsf.edu/mailman/listinfo/chimera-users</a><br class=""></div></div><br class=""></div></div>_______________________________________________<br class="">Chimera-users mailing list: <a class="" href="mailto:Chimera-users@cgl.ucsf.edu" target="_blank">Chimera-users@cgl.ucsf.edu</a><br class="">Manage subscription: <a class="" href="https://plato.cgl.ucsf.edu/mailman/listinfo/chimera-users" target="_blank">https://plato.cgl.ucsf.edu/mailman/listinfo/chimera-users</a><br class=""></div></div><br class=""></div><div class=""><br class=""></div>_______________________________________________<br class="">Chimera-users mailing list: <a class="" href="mailto:Chimera-users@cgl.ucsf.edu" target="_blank">Chimera-users@cgl.ucsf.edu</a><br class="">Manage subscription: <a class="" href="https://plato.cgl.ucsf.edu/mailman/listinfo/chimera-users" target="_blank">https://plato.cgl.ucsf.edu/mailman/listinfo/chimera-users</a><br class=""></div></div><br class=""></div></blockquote><div class=""><br class=""></div></div>_______________________________________________<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=""></div></blockquote></div><br class=""></div></body></html>