10. Merging 3D Maps and Atomic Structures: Rigid Fitting
Once we have the predicted model of any structural element included in our map, to fit that model in the volume constitutes the next step in the modeling workflow. Two protocols have been included in Scipion with this purpose, phenix-dock in map (Appendix PHENIX Dock in map [Liebschner et al., 2019]) and chimerax-rigid fit (Appendix CHIMERAX Rigid fit). The first one allows automatic fitting of models in maps, while the second one only does it when model and map are quite close, thus requiring manual fitting in advance. Although there is no a general rule to fit map and model, because it will depend on the particular problem and on our previous knowledge, in this tutorial we are going to use phenix-dock in map application first, followed by the final Fit in map in chimerax-rigid fit. Observe these two new steps in the modeling Scipion workflow in Fig. 10.1.
10.1. Initial rigid fit with PHENIX dock in map
Open phenix-dock in map protocol (Fig. 10.2 (1)), and complete the form with the the extracted map asymmetric unit (2), the map resolution (3), the model of atomic structure previously saved in ChimeraX (4), and the number of copies of this atomic structure that we’d like to fit in the map, 1 in this case (5). As an additional exercise you can check the result of fitting two copies of this structure in the initial input map.
After executing the protocol phenix-dock in map (Fig. 10.2 (6)), you can check the docking results clicking in Analyze Results (7). ChimeraX graphics window will be opened (Fig. 10.3) showing the map and the atomic structure modeled in its initial location (pink) and fitted in the map (green) (Fig. 10.3 (1)).
A rough inspection of the placed_model in Fig. 10.3 (remark the location of the HEME group, for example, which should be moved slightly to the right side) shows that the fitting could be improved a little. The second protocol, chimerax-rigid fit, will help in this purpose.
10.2. Completing rigid fit with CHIMERAX rigid fit
NOTE before starting!!!:As we already advised previously, we are going to use a ChimeraX-derived protocol (chimerax-rigid fit, Appendix CHIMERA Rigid fit). Remark that this use of ChimeraX is completely different from the use of ChimeraX as a visualization tool. By using the ChimeraX graphics window, opening it from the Scipion button Analyze Results, we can observe protocol results but we CANNOT save anything in Scipion. However, using ChimeraX as a tool, as it is the case in Scipion ChimeraX-derived protocols, we can perform different tasks, taking advantage of the available ChimeraX tools and, finally, we CAN save the obtained results and the working session in Scipion.
Once opened the ChimeraX graphical window, we can complete the fitting of the model to the map, by command line or through the ChimeraX GUI.
- By ChimeraX command line, considering that map and model have ID numbers #2 and #3 (Fig. 10.5 (B)):
fitmap #3 inMap #2
By the ChimeraX GUI: Select in the upper main menu Tools -> Volume Data -> Fit in Map. A small window will be opened (Fig. 10.5 (A)). Select the appropriate model to fit in the map and press Fit (1) to allow the automatic rigid fitting.
scipionwrite #3 prefix Hgb_alpha_ exit